Papers for Monday, Sep 18 2023

We seek to identify the source regions of the slow solar wind (SSW) through combining models with in situ observations. We leverage an opportune conjunction between Solar Orbiter and Parker Solar Probe (PSP) during PSP Encounter 11 to include compositional diagnostics from the Solar Orbiter heavy ion sensor (HIS) as these variations provide crucial insights into the origin and nature of the solar wind. We use Potential Field Source Surface (PFSS) and Magnetohydrodynamic (MHD) models to connect the observed plasma at PSP and Solar Orbiter to its origin footpoint in the photosphere, and compare these results with the in situ measurements. A very clear signature of a heliospheric current sheet (HCS) crossing as evidenced by enhancements in low FIP elements, ion charge state ratios, proton density, low-Alfvenicity, and polarity estimates validates the combination of modeling, data, and mapping. Fast wind from a small equatorial coronal hole (CH) with low ion charge state ratios, low FIP bias, high-Alfvenicity, and low footpoint brightness mostly fits together, yet includes anomalously low alpha particle abundance. We identify slow wind from many different sources, with broad variation in composition and Alfvenicity. We distinguish between classical non-Alfvenic and high-Alfvenicity SSW and make the association of low-Alfvenicity, decreased alpha-to-proton abundance, high charge state ratios, and FIP bias with streamer wind while intermediate alpha abundance, high-Alfvenicity, and dips in ion charge state ratios correspond to CH boundaries. Through this comprehensive analysis, we highlight the power of multi-instrument conjunction studies in assessing the sources of the solar wind.

We present results from a pair of high resolution, long timescale ($\sim10^5 GM/c^3$), global, three dimensional magnetohydrodynamical accretion disk simulations with differing initial magnetic plasma $\beta$ in order to study the effects of initial toroidal field strength on production of large-scale poloidal field. We initialize our disks in approximate equilibrium with purely toroidal magnetic fields of strength $\beta_0=5$ and $\beta_0=200$. We also perform a limited resolution study. We find that simulations of differing field strength diverge early in their evolution and remain distinct over the time studied, indicating that initial magnetic conditions leave a persistent imprint in our simulations. Neither simulation enters the Magnetically Arrested Disk (MAD) regime. Both simulations are able to produce poloidal fields from initially-toroidal fields, with the $\beta_0=5$ simulation evolving clear signs of a large-scale poloidal field. We make a cautionary note that computational artifacts in the form of large-scale vortices may be introduced in the combination of initially-weak field and disk-internal mesh refinement boundaries, as evidenced by the production of an $m=1$ mode overdensity in the weak field simulation. Our results demonstrate that the initial toroidal field strength plays a vital role in simulated disk evolution for the models studied.

Multiphase galactic outflows, generated by supernova feedback, are likely to be more metal-rich than the interstellar media from which they are driven due to incomplete mixing between supernova ejecta and the ambient ISM. This enrichment is important for shaping galactic metallicities and metallicity gradients, but measuring it quantitatively from simulations requires resolution high enough to resolve mass, momentum and energy exchanges between the different phases of the outflows. In this context, we present simulations of outflows, driven by SN feedback, conducted using \textsc{Quokka}, a new GPU-optimised AMR radiation-hydrodynamics code. This code allows us to reach combinations of resolution, simulation volume, and simulation duration larger than those that have previously been possible, and to resolve all gas phases from cold neutral medium, $T \sim 100$ K, to hot ionised gas, $T \gtrsim 10^7$ K. In this, a first of a series of papers exploring generation and evolution of multiphase outflows from a wide range of galactic environments and star formation rates, we quantify the extent of selective metal loading in Solar neighbourhood-like environments. We explain the selective metal loading we find as a result of the transport of metals within and between phases, a phenomenon we can study owing to the parsec-scale resolution that our simulations achieve. We also quantify the sensitivity of metal loading studies to numerical resolution, and present convergence criteria for future studies.

We exploit moderately resolved [OIII], [CII] and dust continuum ALMA observations to derive the gas density ($n$), the gas-phase metallicity ($Z$) and the deviation from the Kennicutt-Schmidt (KS) relation ($\kappa_s$) on ~sub-kpc scales in the interstellar medium (ISM) of five bright Lyman Break Galaxies at the Epoch of Reionization ($z\approx 7$). To do so, we use GLAM, a state-of-art, physically motivated Bayesian model that links the [CII] and [OIII] surface brightness ($\Sigma_{\rm [CII]}$, $\Sigma_{\rm [OIII]}$) and the SFR surface density ($\Sigma_{\rm SFR}$) to $n$, $\kappa_s$, and $Z$. All five sources are characterized by a central starbursting region, where the $\Sigma_{\rm gas}$ vs $\Sigma_{\rm SFR}$ align ~10x above the KS relation ($\kappa_s\approx10$). This translates into gas depletion times in the range $t_{\rm dep}\approx 80-250$ Myr. The inner starbursting centers are characterized by higher gas density ($\log (n/{\rm cm^{-3}}) \approx 2.5-3.0$) and higher metallicity ($\log (Z/Z_{\odot}) \approx -0.5$) than the galaxy outskirts. We derive marginally negative radial metallicity gradients ($\nabla \log Z \approx -0.03 \pm 0.07$dex/kpc), and a dust temperature ($T_d\approx$32-38 K) that anticorrelates with the gas depletion time.

The metal-poor stars in the Galactic halo are thought to show the imprints of the first (PopIII) stars, and thus provide a glance at the first episodes of star formation. In this work, we aim at understanding whether all very metal-poor stars formed in environments polluted by PopIII supernovae (SNe) and at what level. With a general parametric model for early metal enrichment, we study the chemical abundances (from C to Zn) of an environment imprinted by a single PopIII SN. We investigate how these abundances depend on the initial mass and internal mixing of PopIII stars, as well as on their SN explosion energy. We then study how subsequent generations of normal (PopII) SNe affect the PopIII chemical signatures. By comparing the observed chemical abundances with our model predictions, we show that stars with [C/Fe]$>+2.5$ form in environments polluted purely by low-energy PopIII SNe ($E_{\rm SN}<2\times 10^{51}$erg). At lower [C/Fe], stars can be imprinted either by PopIII only, or also by normal PopII SNe. The probability of being enriched by PopII SNe increases as [C/Fe] decreases. When PopII stars contribute more to the pollution, they wash out the diverse chemical peculiarities left by the different PopIII SNe, and the chemical dispersion between their descendants decreases. We conclude that C-normal stars ($\rm [C/Fe] \leq +0.7$) have likely been enriched by PopII SNe at a $\geq 50\%$ level and we identify in the abundance scatter a key diagnostic to pinpoint the signature of PopIII SNe.

A number of recent studies have found evidence for a tension between observations of large-scale structure (LSS) and the predictions of the standard model of cosmology with the cosmological parameters fit to the cosmic microwave background (CMB). The origin of this '$S_8$ tension' remains unclear, but possibilities include new physics beyond the standard model, unaccounted for systematic errors in the observational measurements and/or uncertainties in the role that baryons play. Here we carefully examine the latter possibility using the new FLAMINGO suite of large-volume cosmological hydrodynamical simulations. We project the simulations onto observable harmonic space and compare with observational measurements of the power and cross-power spectra of cosmic shear, CMB lensing, and the thermal Sunyaev-Zel'dovich (tSZ) effect. We explore the dependence of the predictions on box size and resolution, cosmological parameters including the neutrino mass, and the efficiency and nature of baryonic 'feedback'. Despite the wide range of astrophysical behaviours simulated, we find that baryonic effects are not sufficiently large to remove the $S_8$ tension. Consistent with recent studies, we find the CMB lensing power spectrum is in excellent agreement with the standard model, whilst the cosmic shear power spectrum, tSZ effect power spectrum, and the cross-spectra between shear, CMB lensing, and the tSZ effect are all in varying degrees of tension with the CMB-specified standard model. These results suggest that some mechanism is required to slow the growth of fluctuations at late times and/or on non-linear scales, but that it is unlikely that baryon physics is driving this modification.

Recent works have suggested that energy balance spectral energy distribution (SED) fitting codes may be of limited use for studying high-redshift galaxies for which the observed ultraviolet and far-infrared emission are offset (spatially `decoupled'). It has been proposed that such offsets could lead energy balance codes to miscalculate the overall energetics, preventing them from recovering such galaxies' true properties. In this work, we test how well the SED fitting code Magphys can recover the stellar mass, star formation rate (SFR), specific SFR, dust mass and luminosity by fitting 6,706 synthetic SEDs generated from four zoom-in simulations of dusty, high-redshift galaxies from the FIRE project via dust continuum radiative transfer. Comparing our panchromatic results (using wavelengths 0.4-500$\mu$m, and spanning $1<z<8$) with fits based on either the starlight ($\lambda_\mathrm{eff} \le 2.2\,\mu$m) or dust ($\ge 100\,\mu$m) alone, we highlight the power of considering the full range of multi-wavelength data alongside an energy balance criterion. Overall, we obtain acceptable fits for 83 per cent of the synthetic SEDs, though the success rate falls rapidly beyond $z \approx 4$, in part due to the sparser sampling of the priors at earlier times since SFHs must be physically plausible (i.e. shorter than the age of the Universe). We use the ground truth from the simulations to show that when the quality of fit is acceptable, the fidelity of Magphys estimates is independent of the degree of UV\FIR offset, with performance very similar to that previously reported for local galaxies.

Supermassive black holes in active galactic nuclei (AGN) are known to launch relativistic jets, which are observed across the entire electromagnetic spectrum and are thought to be efficient particle accelerators. Their primary radiation mechanism for radio emission is polarized synchrotron emission produced by a population of non-thermal electrons. In this Letter, we present a global general relativistic magnetohydrodynamical (GRMHD) simulation of a magnetically arrested disk (MAD). After the simulation reaches the MAD state, we show that waves are continuously launched and propagate along the shear layer, where they shear magnetic field lines and alter the shear layer geometry. We then perform polarized radiation transfer calculations of our GRMHD simulation and find signatures of the waves in both total intensity and linear polarization, effectively lowering the fully resolved polarization fraction. The tell-tale polarization signatures of the waves could be observable by future Very Long Baseline Interferometric observations, e.g., by the next-generation Event Horizon Telescope.

We simulate a black-hole accretion disk system with full-transport general relativistic neutrino radiation magnetohydrodynamics (GR$\nu$RMHD) for 1.2 seconds. This system is likely to form after the merger of two compact objects and is thought to be a robust site of $r$-process nucleosynthesis. We consider the case of a black-hole accretion disk arising from the merger of two neutron stars. Our simulation time coincides with the nucleosynthesis timescale of the $r$ process ($\sim$ 1 second). Because these simulations are time consuming, it is common practice to run for `short' duration of approximately 0.1 to 0.3 seconds. We analyze the nucleosynthetic outflow from this system and compare the results between stopping at 0.12 and 1.2 seconds respectively. We find that the addition of mass ejected in the longer simulation as well as more favorable thermodynamic conditions from emergent viscous ejecta greatly impacts the nucleosynthetic outcome. We quantify the error in nucleosynthetic outcomes between short and long cuts.

For neutron stars, there exist universal relations insensitive to the equation of states, the so called I-Love-Q relations, which show the connections among the moment of inertia, tidal Love number and quadrupole moment. In this paper, we show that these relations also applies to dark stars, bosonic or fermionic. The relations can be extended to a higher range of the variables, as those curves all approximate the ones generated by a polytropic equation of state, when taking the low density (pressure) limit.

We construct complex seismic models of two high-amplitude delta Sct stars, AE UMa and RV Ari, each pulsating in two radial modes: fundamental and first overtone. The models reproduce, besides the frequencies of two radial modes, also the amplitude of bolometric flux variations (the parameter f) for the dominant mode. Applying the Monte Carlo-based Bayesian analysis, we derive strong constraints, on the parameters of the model as well as on the free parameters of the theory. A vast majority of seismic models of the two stars are just at the beginning of hydrogen-shell burning and a small fraction is at the very end of an overall contraction. The stars have a similar age of about 1.6 Gyr for the hydrogen-shell burning phase. Both stars have unusual low overshooting from the convective core; about 0.02 and 0.004 of the pressure scale height for AE UMa and RV Ari, respectively. This result presumably indicates that overshooting should vary with time and scale with a decreasing convective core. The efficiency of convection in the envelope of both stars is rather low and described by the mixing length parameter alphaMLT of about 0.3-0.6. The third frequency of RV Ari, confirmed by us in the TESS photometry, can only be associated with mixed nonradial modes l=1, g4-g8 or l=2, g10-g12. We include the dipole mode into our Bayesian modelling and demonstrate its huge asteroseismic potential.

The absence of a breakthrough in directly observing dark matter (DM) through prominent large-scale detectors motivates the development of novel tabletop experiments probing more exotic regions of the parameter space. If DM contains ultralight bosonic particles, they would behave as a classical wave and could manifest through an oscillating force on baryonic matter that is coherent over $\sim 10^6$ periods. Our Helium ultraLIght dark matter Optomechanical Sensor (HeLIOS) uses the high-$Q$ acoustic modes of superfluid helium-4 to resonantly amplify this signal. A superconducting re-entrant microwave cavity enables sensitive optomechanical readout ultimately limited by thermal motion at millikelvin temperatures. Pressurizing the helium allows for the unique possibility of tuning the mechanical frequency to effectively broaden the DM detection bandwidth. We demonstrate the working principle of our prototype HeLIOS detector and show that future generations of HeLIOS could explore unconstrained parameter space for both scalar and vector ultralight DM after just an hour of integration time.

Gamma-Ray Bursts are the most luminous events in the Universe, and are excellent laboratories to study extreme physical phenomena in the cosmos. Despite a long trajectory of progress in understanding these highly energetic events, there are still many observed features that are yet to be fully explained. Observations of the jet opening angle of long gamma-ray bursts (LGRBs) suggest that LGRB jets are narrower for those GRBs at higher redshift. This phenomenon has been explained in the context of collimation by the stellar envelope, with denser (lower metallicity) stars at higher redshifts able to collimate the jet more effectively. However, until now, the dependence of jet opening angle on the properties of the central engine has not been explored. We investigate the effect of black hole spin on the jet collimation angle for a magnetically launched jet, using the General Relativistic Radiation Magnetohydrodynamical (GRRMHD) code nu-bhlight. We present 3D results for a range of spin values. The simulations show that higher spinning black holes tend to create narrower jets. If indeed LGRB progenitors in the early universe are able to produce black hole central engines with higher spin, this could account for at least some of the observed jet opening angle-redshift correlation.

A new generation of radio telescopes with excellent sensitivity, instantaneous {\it uv} coverage, and large fields of view, are providing unprecedented opportunities for performing commensal transient searches. Here we present such a commensal search in deep observations of short gamma-ray burst fields carried out with the MeerKAT radio telescope in South Africa at 1.3 GHz. These four hour observations of eight different fields span survey lengths of weeks to months. We also carry out transient searches in time slices of the full observations, at timescales of 15 minutes, and 8 seconds. We find 122 variable sources on the long timescales, of which 52 are likely active galactic nuclei, but there are likely also some radio flaring stars. While the variability is intrinsic in at least two cases, most of it is consistent with interstellar scintillation. In this study, we also place constraints on transient rates based on state-of-the-art transient simulations codes. We place an upper limit of $2\times10^{-4}$ transients per day per square degree for transients with peak flux of 5 mJy, and an upper limit of $2.5\times10^{-2}$ transients per day per square degree for transients with a fluence of 10 Jy ms, the minimum detectable fluence of our survey.

We intend to study a modified version of the planar Circular Restricted Three-Body Problem (CRTBP) by incorporating several perturbing parameters. We consider the bigger primary as an oblate spheroid and emitting radiation while the small primary has an elongated body. We also consider the perturbation from a disk-like structure encompassing this three-body system. First, we develop a mathematical model of this modified CRTBP. We have found there exist five equilibrium points in this modified CRTBP model, where three of them are collinear and the other two are non-collinear. Second, we apply our modified CRTBP model to the Sun-Haumea system by considering several values of each perturbing parameter. Through our numerical investigation, we have discovered that the incorporation of perturbing parameters has resulted in a shift in the equilibrium point positions of the Sun-Haumea system compared to their positions in the classical CRTBP. The stability of equilibrium points is investigated. We have shown that the collinear equilibrium points are unstable and the stability of non-collinear equilibrium points depends on the mass parameter $\mu$ of the system. Unlike the classical case, non-collinear equilibrium points have both a maximum and minimum limit of $\mu$ for achieving stability. We remark that the stability range of $\mu$ in non-collinear equilibrium points depends on the perturbing parameters. In context of the Sun-Haumea system, we have found that the non-collinear equilibrium points are stable.

We present an extension of the Coronal Reconstruction Onto B-Aligned Regions (CROBAR) method to Linear Force Free Field (LFFF) extrapolations, and apply it to the reconstruction of a set of AIA, MDI, and STEREO EUVI data. The results demonstrate that CROBAR can not only reconstruct coronal emission structures, but also that it can help constrain the coronal field extrapolations via the LFFF's helicity $\alpha$ parameter. They also provide a real-world example of how CROBAR can easily incorporate information from multiple perspectives to improve its reconstructions, and we also use the additional perspectives to help validate the reconstructions. We furthermore touch on the use of real-world emission passbands rather than idealized power-law type ones using DEMs. We conclude with a comparison of CROBAR generated emission to observed emission and those produced with idealized DEM based power-laws. These results further illustrate the promise of CROBAR for real-world applications, and we make available a preliminary release of the software available for download.

As the maximum RF input and output frequencies of the integrated data converters in RFSoC increase, it becomes practical to digitize and synthesize RF signals in the majority of C band directly without analogue up and down mixing circuits. The elimination of the mixer circuits can significantly simplify the architecture of the receivers or readouts for radio astronomy telescopes. For the systems with large bandwidth or high channel counts, direct sampling can dramatically reduce the size and cost of overall system. This paper with focus on summarising part of the preliminary characterization results for direct sampling with RFSoC data converters in higher order Nyquist zones.

Optical quasi-periodic oscillations (QPOs) are the most preferred signs of sub-pc binary black hole (BBH) systems in AGN. In this manuscript, robust optical QPOs are reported in quasar SDSS J1609+1756 at $z=0.347$. In order to detect reliable optical QPOs, four different methods are applied to analyze the 4.45 years-long ZTF g/r/i-band light curves of SDSS J1609+1756, direct fitting results by sine function, Generalized Lomb-Scargle periodogram, Auto-Cross Correlation Function and Weighted Wavelet Z-transform method. The Four different methods can lead to well determined reliable optical QPOs with periodicities $\sim340$ days with confidence levels higher than 5$\sigma$, to guarantee the robustness of the optical QPOs in SDSS J1609+1756. Meanwhile, based on simulated light curves through CAR process to trace intrinsic AGN activities, confidence level higher than $3\sigma$ can be confirmed that the optical QPOs are not mis-detected in intrinsic AGN activities, re-confirming the robust optical QPOs and strongly indicating a central sub-pc BBH system in SDSS J1609+1756. Furthermore, based on apparent red-shifted shoulders in broad Balmer emission lines in SDSS J1609+1756, space separation of the expected central BBH system can be estimated to be smaller than $107\pm60$ light-days, accepted upper limit of total BH mass $\sim(1.03\pm0.22)\times10^8{\rm M_\odot}$. Therefore, to detect and report BBH system expected optical QPOs with periodicities around 1 year is efficiently practicable through ZTF light curves, and combining with peculiar broad line emission features, further clues should be given on space separations of BBH systems in broad line AGN in the near future.

Early dark energy is a promising potential resolution of the Hubble tension. Unfortunately, many models suffer from the need to fine-tune their initial conditions to ensure that the epoch of early dark energy coincides with matter-radiation equality. We propose a class of attractive early dark energy models where this coincidence arises naturally as a saddle point of a dynamical system that attracts a large volume of phase-space trajectories regardless of the initial conditions. The system approaches a global dark energy attractor at late-times. Our framework therefore unifies early and late dark energy using a single scalar degree of freedom. We analyze a fiducial attractive early dark energy model and find that it is disfavored by cosmological data due to the presence of a long-lived saddle point in the matter era where the scalar plays the role of an additional component of (non-clustering) dark matter. Our investigations provide lessons for future model-building efforts aimed at constructing viable attractive early dark energy models.

Jetted active galactic nuclei (AGNs) are the principal extragalactic $\gamma$-ray sources. Fermi-detected high-redshift ($z>3$) blazars are jetted AGNs thought to be powered by massive, rapidly spinning supermassive black holes (SMBHs) in the early universe ($<2$ Gyr). They provide a laboratory to study early black hole (BH) growth and super-Eddington accretion -- possibly responsible for the more rapid formation of jetted BHs. However, previous virial BH masses of $z>3$ blazars were based on C IV in the observed optical, but C IV is known to be biased by strong outflows. We present new Gemini/GNIRS near-IR spectroscopy for a sample of nine $z>3$ Fermi $\gamma$-ray blazars with available multi-wavelength observations that maximally sample the spectral energy distributions (SEDs). We estimate virial BH masses based on the better calibrated broad H$\beta$ and/or Mg II . We compare the new virial BH masses against independent mass estimates from SED modeling. Our work represents the first step in campaigning for more robust virial BH masses and Eddington ratios for high-redshift Fermi blazars. Our new results confirm that high-redshift Fermi blazars indeed host overly massive SMBHs as suggested by previous work, which may pose a theoretical challenge for models of the rapid early growth of jetted SMBHs.

We study a compact nucleus embedded in an early-type dwarf galaxy, MATLAS-167, which is in the process of disruption by the tidal force of the neighboring giant S0 galaxy, NGC 936, in a group environment. Using the imaging data of the MATLAS survey, we analyze the stellar tidal tail of MATLAS-167 and its central compact nucleus, designated as NGC 936_UCD. We find that NGC 936_UCD has a luminosity of M$_{g}$ = $-$11.43$\pm$0.01 mag and a size of 66.5$\pm$17 pc, sharing the global properties of Ultra Compact Dwarf galaxies (UCDs) but significantly larger and brighter compared to the typical UCD populations observed in the Virgo cluster. By integrating the total luminosity of both the tidal stream and MATLAS-167, we estimate that the disrupted dwarf progenitor possesses a luminosity of M$_{g}$ = $-$15.92$\pm$0.06 mag, a typical bright dE luminosity. With the help of the optical spectrum observed by the SDSS survey, we derive the simple stellar population properties of NGC 936_UCD: a light-weighted age of 5.6$\pm$0.7 Gyr and metallicity of [Z/H] = $-$0.83$\pm$0.3 dex. Our findings suggest that tidal threshing is a possible formation mechanism of bright UCD populations in close proximity to giant galaxies.

In this paper, we present an analysis of the pulsating behavior of Kepler target KIC 9845907. Using the data from Kepler, we detected 85 significant frequencies, including the first overtone $f_{1}$ = 17.597 day$^{-1}$ as the dominant frequency, the non-radial independent frequency $f_{3}$ = 31.428 day$^{-1}$ ($\ell$=1), as well as two modulation terms $f_{m1}$ = 0.065 day$^{-1}$ and $f_{m2}$ = 1.693 day$^{-1}$. We found fourteen pairs of triplet structures with $f_{m1}$ or $f_{m2}$, four pairs of which can further form quintuplet structures. We note these are the most intriguing features discovered in this study and they were recognized for the first time in $\delta$ Scuti stars. We discussed several possible explanations, i.e., beating, the Blazhko effect, combination mode hypothesis, nonlinear mode coupling, large separation, and stellar rotational splitting for these equidistant structures. Our asteroseismic models indicate this modulation with $f_{m1}$ might be related to the rotational splitting. The study of more $\delta$ Scuti stars with triplet and/or quintuplet structures using high-precision space photometry would be helpful to further explore its origin.

The Internal Linear Combination (ILC) method is commonly employed to extract the cosmic microwave background (CMB) signal from multi-frequency observation maps. However, the performance of the ILC method tends to degrade when the signal-to-noise ratio (SNR) is relatively low, especially when measuring the primordial $B$-modes for detecting the primordial gravitational waves. To address this issue, we propose an enhanced version of needlet ILC (NILC) method on $B$ map called constrained NILC which is more suitable for situations with low SNR by adding additional prior foreground information. We illustrate our methods using mock data generated from the combination of WMAP, Planck and a ground-based experiment in the northern hemisphere, and the chosen noise level for the ground-based experiment are very conservative which can be easily achieved in the very near future. The results show that the level of foreground residual can be well controlled. In comparison to the standard NILC method, which introduces a bias to the tensor-to-scalar ratio ($r$) of approximately $0.05$, the constrained NILC method exhibits a significantly reduced bias of only $5\times10^{-3}$ towards $r$ which is much smaller than the statistical error.

Proton auroras are widely observed on the day side of Mars, identified as a significant intensity enhancement in the hydrogen Ly alpha (121.6 nm) emission between 120 and 150~km altitudes. Solar wind protons penetrating as energetic neutral atoms into the Martian thermosphere are thought to be responsible for these auroras. Understanding proton auroras is therefore important for characterizing the solar wind interaction with the atmosphere of Mars. Recent observations of spatially localized "patchy" proton auroras suggest a possible direct deposition of protons into the atmosphere of Mars during unstable solar wind conditions. Here, we develop a purely data-driven model of proton auroras using Mars Atmosphere and Volatile EvolutioN (MAVEN) in situ observations and limb scans of Ly alpha emissions between 2014 and 2022. We train an artificial neural network that reproduces individual Ly alpha intensities with a Pearson correlation of 0.95 along with a faithful reconstruction of the observed Ly alpha emission altitude profiles. By performing a SHapley Additive exPlanations (SHAP) analysis, we find that Solar Zenith Angle, seasonal CO2 atmosphere variability, solar wind temperature, and density are the most important features for the modelled proton auroras. We also demonstrate that our model can serve as an inexpensive tool for simulating and characterizing Ly alpha response under a variety of seasonal and upstream solar wind conditions.

GPU computing is popular due to the calculation potential of a single card. The N-body integrator GENGA is built to for this, but it suffers a performance penalty on consumer-grade GPUs due to their truncated double precision (FP64) performance. We aim to speed up GENGA on consumer-grade cards by harvesting their high single-precision performance (FP32). We modified GENGA to be able to compute the long-distance forces between bodies in FP32 precision and tested this with 5 experiments. We ran simulations with similar initial conditions of 6600 planetesimals in both FP32 and FP64 precision. We also ran simulations that i) began with a mixture of planetesimals and planetary embryos, ii) planetesimal-driven giant planet migration, and iii) terrestrial planet formation with a gas disc. Second, we ran the same simulation beginning with 40 000 planetesimals using both FP32 and FP64 precision forces on a variety of consumer-grade and Tesla GPUs to measure the performance boost of FP32 computing. There are no statistical differences when running in FP32 or FP64 precision that can be attributed to the force prescription rather than stochastic effects. The uncertainties in energy are almost identical when using both precisions. However, the uncertainty in the angular momentum using FP32 rather than FP64 precision long-range forces is about two orders of magnitude greater, but still very low. Running the simulations in single precision on consumer-grade cards decreases running time by a factor of three and becomes within a factor of three of a Tesla A100 GPU. Additional tuning speeds up the simulation by a factor of two across all types of cards. The option to compute the long-range forces in single precision in GENGA when using consumer-grade GPUs dramatically improves performance at a little penalty to accuracy. There is an additional environmental benefit because it reduces energy usage.

This index contains the proceedings submitted to the 38th International Cosmic Ray Conference (ICRC 2023) in the name of the CTA consortium.

We present the analysis of spatially-resolved spectroscopic observations of the born-again planetary nebula (PN) A 78 that are used to investigate the chemistry and physical properties of its three main morphological components, namely the inner knots, its eye-like structure, and the low surface-brightness outer nebula. The H-poor chemical abundances of the inner knots confirm the born-again nature of A 78, with a N/O abundances ratio consistent with the predictions of very late thermal pulses (VLTP). On the other hand, the high Ne/O is not expected in VLTP events, which prompts a possible different evolutionary path may be involving a binary system. The low N/O ratio and He/H abundances of the outer shell are indicative of a low-mass progenitor, whereas the chemical abundances of the eye-like structure, which results from the interaction between the H-poor born-again material and the outer nebula, evidence their mixture. Unlike previous works, the extinction is found to be inhomogeneous, being much higher towards the H-poor inner knots, where the presence of large amounts of C-rich dust has been previously reported. Dust-rich material seems to diffuse into outer nebular regions, resulting in zones of enhanced extinction.

In this paper we develop a framework for studying unstratified, magnetised eccentric discs and compute uniformly precessing eccentric modes in a cylindrical annulus which provide convenient initial conditions for numerical simulations. The presence of a magnetic field in an eccentric disc can be described by an effective gas with a modified equation of state. At magnetic field strengths relevant to the magneto-rotational instability the magnetic field has negligible influence on the evolution of the eccentric disc, however the eccentric disc can significantly enhance the magnetic field strength over that in the a circular disc. We verify the suitability of these eccentric disc solutions by carrying out 2D simulations in RAMSES. Our simulated modes (in 2D) follow a similar evolution to the purely hydrodynamical modes, matching theoretical expectations, provided they are adequately resolved. Such solutions will provide equilibrium states for studies of the eccentric magneto-rotational instability and magnetised parametric instability in unstratified discs and are useful for exploring the response of disc turbulence on top of a fluid flow varying on the orbital timescale.

The large-scale distribution of baryons is sensitive to gravitational collapse, mergers, and galactic feedback. Known as the Cosmic Web, its large-scale structure (LSS) can be classified as halos, filaments, and voids. Fast Radio Bursts (FRBs) are extragalactic sources that undergo dispersion along their propagation paths. They provide insight into ionised matter along their sightlines via their dispersion measures (DMs), and have been investigated as probes of the LSS baryon fraction, the diffuse baryon distribution, and of cosmological parameters. We use the cosmological simulation IllustrisTNG to study FRB DMs accumulated while traversing different types of LSS. We combine methods for deriving electron density, classifying LSS, and tracing FRB sightlines. We identify halos, filaments, voids, and collapsed structures along random sightlines and calculate their DM contributions. We analyse the redshift-evolving cosmological DM components of the Cosmic Web. We find that the filamentary contribution dominates, increasing from ~71% to ~80% on average for FRBs originating at z=0.1 vs z=5, while the halo contribution falls, and the void contribution remains consistent to within ~1%. The majority of DM variance originates from halos and filaments, potentially making void-only sightlines more precise probes of cosmological parameters. We find that, on average, an FRB originating at z=1 will intersect ~1.8 foreground collapsed structures, increasing to ~12.4 structures for a z=5 FRB. The impact parameters between our sightlines and TNG structures of any mass appear consistent with those reported for likely galaxy-intersecting FRBs. However, we measure lower average accumulated DMs from these structures than the $\sim90\;{\rm pc\;cm^{-3}}$ DM excesses reported for these literature FRBs, indicating some DM may arise beyond the structures themselves.

The study of phosphorus chemistry in the interstellar medium has become a topic of growing interest in astrobiology, because it is plausible that a wide range of P-bearing molecules were introduced in the early Earth by the impact of asteroids and comets on its surface, enriching prebiotic chemistry. Thanks to extensive searches in recent years, it has become clear that P mainly appears in the form of PO and PN in molecular clouds and star-forming regions. Interestingly, PO is systematically more abundant than PN by factors typically of $\sim1.4-3$, independently of the physical properties of the observed source. In order to unveil the formation routes of PO and PN, in this work we introduce a mathematical model for the time evolution of the chemistry of P in an interstellar molecular cloud and analyze its associated chemical network as a complex dynamical system. By making reasonable assumptions, we reduce the network to obtain explicit mathematical expressions that describe the abundance evolution of P-bearing species and study the dependences of the abundance of PO and PN on the system's kinetic parameters with much faster computation times than available numerical methods. As a result, our model reveals that the formation of PO and PN is governed by just a few critical reactions, and fully explains the relationship between PO and PN abundances throughout the evolution of molecular clouds. Finally, the application of Bayesian methods constrains the real values of the most influential reaction rate coefficients making use of available observational data.

PSR B1259-63 is a gamma-ray binary system with a radio pulsar orbiting an O9.5Ve star, LS 2883, with a period of ~3.4 yr. Close to the periastron the system is detected at all wavelengths, from radio to the TeV energies. The emission in this time period is believed to originate from the interaction of LS 2883 and pulsar's outflows. The observations of 4 periastra passages taken in 2010-2021 show strong correlation of the radio and X-ray lightcurves with two peaks just before and after the periastron. The observations of the latest 2021 periastron passage reveal the presence of the 3rd X-ray peak and subsequent disappearance of radio/X-ray flux correlation. In this paper we present the results of our optical, radio and X-ray observational campaigns on PSR B1259-63 performed in 2021 accompanied with the analysis of the publicly available GeV FERMI/LAT data. We compare the properties of different periastron passages, discuss the obtained results and show that they can be explained in terms of the 2-zone emission cone model proposed by us previously.

We investigate the propagation of spherically symmetric shocks in relativistic homologously expanding media with density distributions following a power-law profile in their Lorentz factor. That is, $\rho_{ej} \propto t^{-3}\gamma_{e}(R,t)^{-\alpha}$, where $\rho_{ej}$ is the medium proper density, $\gamma_{e}$ is its Lorentz factor, $\alpha>0$ is constant and $t$, $R$ are the time and radius from the center. We find that the shocks behavior can be characterized by their proper velocity, $U'=\Gamma_s'\beta_s'$, where $\Gamma_s'$ is the shock Lorentz factor as measured in the immediate upstream frame and $\beta_s'$ is the corresponding 3-velocity. While generally, we do not expect the shock evolution to be self-similar, for every $\alpha>0$ we find a critical value $U'_c$ for which a self-similar solution with constant $U'$ exists. We then use numerical simulations to investigate the behavior of general shocks. We find that shocks with $U'>U'_c$ have a monotonously growing $U'$, while those with $U'<U'_c$ have a decreasing $U'$ and will eventually die out. Finally, we present an analytic approximation, based on our numerical results, for the evolution of general shocks in the regime where $U'$ is ultra-relativistic.

We present chemical composition and fundamental parameters (the effective temperature, surface gravity and radius) for four sharp-lined A-type stars ${\gamma}$ Gem (HD 41705), o Peg (HD 214994), ${\theta}$ Vir (HD 114330) and ${\nu}$ Cap (HD 193432). Our analysis is based on a self-consistent model fitting of high-resolution spectra and spectrophotometric observations over a wide wavelength range. We refined the fundamental parameters of the stars with the SME package and verified their accuracy by comparing with the spectral energy distribution and hydrogen line profiles. We found Teff/log g = 9190+/-130 K/3.56+/-0.08, 9600+/-50 K/3.81+/-0.04, 9600+/-140 K/3.61+/-0.12, and 10200+/-220 K/3.88+/-0.08 for ${\gamma}$ Gem, o Peg, ${\theta}$ Vir and ${\nu}$ Cap, respectively. Our detailed abundance analysis employs a hybrid technique for spectrum synthesis based on classical model atmospheres calculated in local thermodynamic equilibrium (LTE) assumption together with the non-LTE (NLTE) line formation for 18 of 26 investigated species. Comparison of the abundance patterns observed in A stars of different types (normal A, Am, Ap) with similar fundamental parameters reveals significant abundance diversity that cannot be explained by the current mechanisms of abundance peculiarity formation in stellar atmospheres. We found a rise of the heavy element (Zn, Sr, Y, Zr, Ba) abundance excess up to +1 dex with Teff increasing from 7200 to 10000 K, with a further decrease down to solar value at Teff = 13000 K, indicating that stars with solar element abundances can be found among late B-type stars.

We present the results of new photometric and spectroscopic observations of WN5+O6 binary V444 Cyg and a detailed analysis of extant spectroscopy and photometry. Using elements of the spectroscopic orbit and assuming e = 0, i = 78{\deg} we determined the masses and orbit sizes of the components of V444 Cyg M_O6 = 26.4 M_Sun, M_WN5 = 10.7 M_Sun, a_O6 = 10.6 R_Sun, a_WN5 = 26.1 R_Sun. Based on new and archival light curves and by applying the Hertzsprung method, we improved the photometrical estimate of secular increase rate of the orbital period in V444 Cyg dP/dt_ph = 0.119+-0.003 s/yr. From the comparison of the new and archival radial velocity curves of V444 Cyg we independently derived the secular orbital period change rate dP/dt_sp = 0.147+-0.032 s/yr, in agreement with the photometrical dP/dt_ph. The obtained secular increase rate of the binary orbital period dP/dt and the mean radii of the components enabled us to estimate the stellar wind mass-loss rate from WR star dM/dt_WN5 = -(6.0+-0.4)*10^-6 M_Sun/yr.

Planetary mass loss is governed by several physical mechanisms, including photoionisation that may impact the evolution of the atmosphere. Stellar radiation energy deposited as heat depends strongly on the energy of the primary electrons following photoionisation and on the local fractional ionisation. All these factors affect the model-estimated atmospheric mass loss rates and other characteristics of the outflow in ways that have not been clearly elucidated. The shape of the XUV stellar spectra influences strongly the photoionisation and heating deposition on the atmosphere. We elaborate on the local and planet-wise effects, to clearly demonstrate the significance of such interactions. Using the PLUTO code, we performed 1D hydrodynamics simulations from Neptune to Jupiter size planets and stars from M dwarfs to Sun-like. Our results indicate a significant decrease of the planetary mass loss rate for all planetary systems when secondary ionisation is taken into account. The mass loss rate is found to decrease by 43$\%$ for the more massive exoplanet to 54$\%$ for the less massive exoplanet orbiting solar-like stars, and up to 52$\%$ for a Jupiter-like planet orbiting a M type star. Our results also indicate much faster ionisation of the atmosphere due to photoelectrons. We built a self-consistent model including secondary ionisation by photoelectron to evaluate its impact on mass loss rates. We find that photoelectrons affect the mass loss rates by factors that are potentially important for planetary evolution theories. We also find that enhanced ionisation occurs at altitudes that are often probed with specific atomic lines in transmission spectroscopy. Future modelling of these processes should include the role of photoelectrons. Finally, we make available a simple yet accurate parameterisation for atomic hydrogen atmospheres.

As the millimeter wavelength range remains a largely unexplored spectral region for galaxies, the IMEGIN large program aims to map the millimeter continuum emission of 22 nearby galaxies at 1.15 and 2 mm. Using the high-resolution maps produced by the NIKA2 camera, we explore the existence of very cold dust and take possible contamination by free-free and synchrotron emission into account. We study the IR-to-radio emission coming from different regions along the galactic plane and at large vertical distances. New observations of NGC 891, using the NIKA2 camera on the IRAM 30m telescope, along with a suite of observations at other wavelengths were used to perform a multiwavelength study of the spectral energy distribution in the interstellar medium in this galaxy. This analysis was performed globally and locally, using the advanced hierarchical Bayesian fitting code, HerBIE, coupled with the THEMIS dust model. Our dust modeling is able to reproduce the near-IR to millimeter emission of NGC 891, with the exception of an excess at a level of 25% obtained by the NIKA2 observations in the outermost parts of the disk. The radio continuum and thermal dust emission are distributed differently in the disk and galaxy halo. Different dusty environments are also revealed by a multiwavelength investigation of the emission features. Our detailed decomposition at millimeter and centimeter wavelengths shows that emission at 1 mm is purely originated by dust. Radio components become progressively important with increasing wavelengths. Finally, we find that emission arising from small dust grains accounts for ~ 9.5% of the total dust mass, reaching up to 20% at large galactic latitudes. Shock waves in the outflows that shatter the dust grains might explain this higher fraction of small grains in the halo.

We apply state-of-the-art, likelihood-free statistical inference (machine-learning-based) techniques for reconstructing the spectral shape of a gravitational wave background (GWB). We focus on the reconstruction of an arbitrarily shaped signal by the LISA detector, but the method can be easily extended to either template-dependent signals, or to other detectors, as long as a characterisation of the instrumental noise is available. As proof of the technique, we quantify the ability of LISA to reconstruct signals of arbitrary spectral shape (blind reconstruction), considering a diversity of frequency profiles, and including astrophysical backgrounds in some cases. As a teaser of how the method can reconstruct signals characterised by a parameter-dependent template (template reconstruction), we present a dedicated study for power-law signals. While our technique has several advantages with respect to traditional MCMC methods, we validate it with the latter for concrete cases. This work opens the door for both fast and accurate Bayesian parameter estimation of GWBs, with essentially no computational overhead during the inference step. Our set of tools will be integrated into the package GWBackFinder, which will be made publicly available in due time.

We re-examine the interactions of ultra-high energy cosmic rays (UHECRs) with photons from the cosmic microwave background (CMB) under a changed, locally non-linear temperature redshift relation $T(z)$. This changed temperature redshift relation is motivated by the postulate of subjecting thermalised and isotropic photon gases such as the CMB to an SU(2) rather than a U(1) gauge group. This modification of $\Lambda$CDM is called SU(2)$_{\rm CMB}$, and some cosmological parameters obtained by SU(2)$_{\rm CMB}$ seem to be in better agreement with local measurements of the same quantities. In this work, we apply the reduced CMB photon density under SU(2)$_{\rm CMB}$ to the propagation of UHECRs. This leads to a higher UHECR flux just below the ankle in the cosmic ray spectrum and slightly more cosmogenic neutrinos under otherwise equal conditions for emission and propagation. Most prominently, the proton flux is significantly increased below the ankle ($5\times10^{18}$ eV) for hard injection spectra and without considering the effects of magnetic fields. The reduction in CMB photon density also favours a decreased cosmic ray source evolution than the best fit using $\Lambda$CDM. In consequence, it seems that SU(2)$_{\rm CMB}$ favours sources that evolve as the star formation rate (SFR), like star burst galaxies (SBG) and gamma ray bursts (GRB), over active galactic nuclei (AGNs) as origins of UHECRs. We conclude that the question about the nature of primary sources of UHECRs is tightly-knit with the actual temperature redshift relation of the CMB.

Molecular clouds in the Galactic center (GC) reprocess radiation from past outbursts of nearby high-energy sources, generating a bright Fe K-alpha fluorescence at 6.4 keV. The closest clouds to the GC are only $\simeq 1.5$ pc from Sgr A*, forming a torus-like structure known as the circumnuclear disk (CND). The study of fluorescence emission can lead to a characterization of the illuminating source(s), the reflecting clouds, and the global geometry of such a system lying in the GC. The primary purpose of our study is to analyze possible fluorescence signals arising in the CND. This signal would allow us to constrain the CND's physical properties and the source-reflector system's geometry. By exploiting the last $\simeq 20$ years of XMM-Newton observations of the GC, we studied the variability of the Fe K-alpha line in the region around Sgr A*. We identified regions with a flux excess and computed the spectrum therein. We then derived the hydrogen column density of the CND after relating the intensity of the 6.4 keV line to the total energy emitted by known transient sources in the region. Starting from data collected in 2019, we find significant line excesses in a region compatible with the eastern portion of the CND. The echo radiation can be linked to the 2013 outburst of the magnetar SGR J1745-2900. We derive a mean effective hydrogen column density of the CND in the eastern region of $\simeq 10^{23}$ cm$^{-2}$. The scenario depicted is physically plausible, given the luminosity, the position of the illuminating source, and the expected density of the CND.

Correlations between intrinsic properties of gamma-ray burst (GRB) light curves provide clues to the nature of the central engine, the jet, and a possible means to standardise GRBs for cosmological use. Here we report on the discovery of a correlation between the intrinsic early time luminosity, $L_{G,\rm 10s}$, measured at rest frame 10s, and the average decay rate measured from rest frame 10s onward, $\alpha_{G,\rm avg>10s}$, in a sample of 13 Fermi Large Array Telescope (LAT) long GRB light curves. We note that our selection criteria, in particular the requirement for a redshift to construct luminosity light curves, naturally limits our sample to energetic GRBs. A Spearman's rank correlation gives a coefficient of -0.74, corresponding to a confidence level of 99.6%, indicating that brighter afterglows decay faster than less luminous ones. Assuming a linear relation with $\log(L_{G,\rm 10s})$, we find $\alpha_{G,\rm avg>10s} = -0.31_{-0.09}^{+0.12}\log(L_{G,\rm 10s}) + 14.43_{-5.97}^{+4.55}$. The slope of -0.31 is consistent at $1\sigma$ with previously identified correlations in the optical/UV and X-ray light curves. We speculate that differences in the rate at which energy is released by the central engine or differences in observer viewing angle may be responsible for the correlation.

Stars in short-period binaries typically have spins that are aligned and synchronized with the orbit of their companion. In triple systems, however, the combination of spin and orbital precession can cause the star's rotation to evolve to a highly misaligned and sub-synchronous equilibrium known as a Cassini state. We identify a population of recently discovered stars that exhibit these characteristics and which are already known to have tertiary companions. These third bodies have a suitable orbital period to allow the inner binary to evolve into the sub-synchronous Cassini state, which we confirm with orbital evolution models. We also compute the expected stellar obliquity and spin period, showing that the observed rotation rates are often slower than expected from equilibrium tidal models. However, we show that tidal dissipation via inertial waves can alter the expected spin-orbit misalignment angle and rotation rate, potentially creating the very slow rotation rates in some systems. Finally, we show how additional discoveries of such systems can be used to constrain the tidal physics and orbital evolution histories of stellar systems.

Recent asteroseismic measurements have revealed a small population of stars in close binaries, containing primaries with extremely slow rotation rates. Such stars defy the standard expectation of tidal synchronization in such systems, but they can potentially be explained if they are trapped in a spin-orbit equilibrium known as Cassini state 2 (CS2). This state is maintained by orbital precession due to an outer tertiary star, and it typically results in a very sub-synchronous rotation rate and high degree of spin-orbit misalignment. We examine how CS2 is affected by magnetic braking and different types of tidal dissipation. Magnetic braking results in a slower equilibrium rotation rate, while tidal dissipation via gravity waves can result in a slightly higher rotation rate than predicted by equilibrium tidal theory, and dissipation via inertial waves can result in much slower rotation rates. For seven binary systems with slowly rotating primaries, we predict the location of the outer tertiary predicted by the CS2 theory. In five of these systems, a tertiary companion has already been detected, although it closer than expected in three of these, potentially indicating tidal dissipation via inertial waves. We also identify a few new candidate systems among a population of eclipsing binaries with rotation measurements via spot modulation.

We study the stellar population properties of 182 spectroscopically-confirmed (MUSE/VLT) Lyman-$\alpha$ emitters (LAEs) and 450 photometrically-selected Lyman-Break galaxies (LBGs) at z = 2.8 - 6.7 in the Hubble eXtreme Deep Field (XDF). Leveraging the combined power of HST and JWST NIRCam and MIRI observations, we analyse their rest-frame UV-through-near-IR spectral energy distributions (SEDs) with MIRI playing a crucial role in robustly assessing the LAE's stellar mass and ages. Our LAEs are low-mass objects (log$_{10}$(M$_\star$[M$_\odot$]) ~ 7.5), with little or no dust extinction (E(B - V) ~ 0.1) and a blue UV continuum slope ($\beta$ ~ -2.2). While 75% of our LAEs are young (< 100 Myr), the remaining 25% have significantly older stellar populations (> 100 Myr). These old LAEs are statistically more massive, less extinct and have lower specific star formation rate (sSFR) compared to young LAEs. Besides, they populate the M$_\star$ - SFR plane along the main-sequence (MS) of star-forming galaxies, while young LAEs populate the starburst region. The comparison between the LAEs properties to those of a stellar-mass matched sample of LBGs shows no statistical difference between these objects, except for the LBGs redder UV continuum slope and marginally larger E(B - V) values. Interestingly, 48% of the LBGs have ages < 10 Myr and are classified as starbursts, but lack detectable Ly$\alpha$ emission. This is likely due to HI resonant scattering and/or selective dust extinction. Overall, we find that JWST observations are crucial in determining the properties of LAEs and shedding light on the properties and similarities between LAEs and LBGs.

The detection of a stochastic Gravitational Wave (GW) background sourced by a cosmological phase transition would allow us to see the early Universe from a completely new perspective, illuminating aspects of Beyond the Standard Model (BSM) physics and inflationary cosmology. In this study, we investigate whether the evolution of the scalar potential of a minimal SM extension after inflation can lead to a strong first-order phase transition. In particular, we focus on a BSM spectator scalar field that is non-minimally coupled to gravity and has a dynamical double-well potential. As inflation ends, the potential barrier diminishes due to the evolution of the curvature scalar. Therefore, a phase transition can proceed through the nucleation of true-vacuum bubbles that collide as they fill the Universe and produce GWs. We consider high and low scales of inflation, while also taking into account a kination period between inflation and the onset of radiation domination. With this prescription, we showcase a proof-of-concept study of a new triggering mechanism for BSM phase transitions in the early Universe, whose GW signatures could potentially be probed with future detectors.

The stability of hybrid stars with first-order phase transitions as determined by calculating fundamental radial oscillation modes is known to differ from the predictions of the widely-used Bardeen--Thorne--Meltzer criterion. We consider the effects of out-of-chemical-equilibrium physics on the radial modes and hence stability of these objects. For a barotropic equation of state, this is done by allowing the adiabatic sound speed to differ from the equilibrium sound speed. We show that doing so extends the stable branches of stellar models, allowing stars with rapid phase transitions to support stable higher-order stellar multiplets similarly to stars with multiple slow phase transitions. We also derive a new junction condition to impose on the oscillation modes at the phase transition. Termed the reactive condition, it is physically motivated, consistent with the generalized junction conditions between two phases, and has the common rapid and slow conditions as limiting cases. Unlike the two common cases, it can only be applied to nonbarotropic stars. We apply this junction condition to hybrid stellar models generated using a two-phase equation of state consisting of nuclear matter with unpaired quark matter at high densities joined by a first-order phase transition and show that like in the slow limiting case, stars that are classically unstable are stabilized by a finite chemical reaction speed.

Third-order lensing statistics contain a wealth of cosmological information that is not captured by second-order statistics. However, the computational effort for estimating such statistics on forthcoming stage IV surveys is prohibitively expensive. We derive and validate an efficient estimation procedure for the three-point correlation function (3PCF) of polar fields such as weak lensing shear. We then use our approach to measure the shear 3PCF and the third-order aperture mass statistics on the KiDS-1000 survey. We construct an efficient estimator for third-order shear statistics which builds on the multipole decomposition of the 3PCF. We then validate our estimator on mock ellipticity catalogs obtained from $N$-body simulations. Finally, we apply our estimator to the KiDS-1000 data and present a measurement of the third-order aperture statistics in a tomographic setup. Our estimator provides a speedup of a factor of $\sim$ 100-1000 compared to the state-of-the-art estimation procedures. It is also able to provide accurate measurements for squeezed and folded triangle configurations without additional computational effort. We report a significant detection of the tomographic third-order aperture mass statistics in the KiDS-1000 data $(\mathrm{S/N}=6.69)$. Our estimator will make it computationally feasible to measure third-order shear statistics in forthcoming stage IV surveys. Furthermore, it can be used to construct empirical covariance matrices for such statistics.

This paper performs the first cosmological parameter analysis of the KiDS-1000 data with second- and third-order shear statistics. This work builds on a series of papers that describe the roadmap to third-order shear statistics. We derive and test a combined model of the second-order shear statistic, namely the COSEBIs and the third-order aperture mass statistics $\langle M_\mathrm{ap}^3\rangle$ in a tomographic set-up. We validate our pipeline with $N$-body simulations that mock the fourth Kilo Degree survey data release. To model the second- and third-order statistics, we use the latest version of \textsc{HMcode2020} for the power spectrum and \textsc{BiHalofit} for the bispectrum. Furthermore, we use an analytic description to model intrinsic alignments and hydro-dynamical simulations to model the effect of baryonic feedback processes. Lastly, we decreased the dimension of the data vector significantly by considering for the $\langle M_\mathrm{ap}^3\rangle$ part of the data vector only equal smoothing radii, making a data analysis of the fourth Kilo Degree survey data release using a combined analysis of COSEBIs third-order shear statistic possible. We first validate the accuracy of our modelling by analysing a noise-free mock data vector assuming the KiDS-1000 error budget, finding a shift in the maximum-a-posterior of the matter density parameter $\Delta \Omega_m< 0.02\, \sigma_{\Omega_m}$ and of the structure growth parameter $\Delta S_8 < 0.05\, \sigma_{S_8}$. Lastly, we performed the first KiDS-1000 cosmological analysis using a combined analysis of second- and third-order shear statistics, where we constrained $\Omega_m=0.248^{+0.062}_{-0.055}$ and $S_8=\sigma_8\sqrt{\Omega_m/0.3}=0.772\pm0.022$. The geometric average on the errors of $\Omega_\mathrm{m}$ and $S_8$ of the combined statistics increased compared to the second-order statistic by 2.2.

We summarize common notations and concepts in the field of Intrinsic Alignments (IA). IA refers to physical correlations involving galaxy shapes, galaxy spins, and the underlying cosmic web. Its characterization is an important aspect of modern cosmology, particularly in weak lensing analyses. This resource is both a reference for those already familiar with IA and designed to introduce someone to the field by drawing from various studies and presenting a collection of IA formalisms, estimators, modeling approaches, alternative notations, and useful references.

Gravitational Waves (GWs) have been detected in the $\sim 100$ Hz and nHz bands, but most of the gravitational spectrum remains unobserved. A variety of detector concepts have been proposed to expand the range of observable frequencies. In this work, we study the capability of GW detectors in the "mid-band", the $\sim 30$ mHz -- 10 Hz range between LISA and LIGO, to measure the signals from and constrain the properties of $\sim 1-100 M_\odot$ compact binaries. We focus on atom-interferometer-based detectors. We describe a Fisher matrix code we use to evaluate their capabilities, and present numerical results for two benchmarks: terrestrial km-scale detectors, and satellite-borne detectors in medium Earth orbit. One unique capability of mid-band GW detectors is pinpointing the location of GW sources on the sky. We demonstrate that a satellite-borne detector could achieve sub-degree sky localization for any detectable source with chirp mass $\mathcal{M}_c \lesssim 50 M_\odot$. We also compare different detector configurations, including different locations of terrestrial detectors and various choices of the orbit of a satellite-borne detector. As we show, a network of only two terrestrial single-baseline detectors or one single-baseline satellite-borne detector would each provide close-to-uniform sky-coverage, with signal-to-noise ratios varying by less than a factor of two across the entire sky. We hope that this work contributes to the efforts of the GW community to assess the merits of different detector proposals.

The next generation of space- and ground-based facilities promise to reveal an entirely new picture of the gravitational wave sky: thousands of galactic and extragalactic binary signals, as well as stochastic gravitational wave backgrounds (SGWBs) of unresolved astrophysical and possibly cosmological signals. These will need to be disentangled to achieve the scientific goals of experiments such as LISA, Einstein Telescope, or Cosmic Explorer. We focus on one particular aspect of this challenge: reconstructing an SGWB from (mock) LISA data. We demonstrate that simulation-based inference (SBI) - specifically truncated marginal neural ratio estimation (TMNRE) - is a promising avenue to overcome some of the technical difficulties and compromises necessary when applying more traditional methods such as Monte Carlo Markov Chains (MCMC). To highlight this, we show that we can reproduce results from traditional methods both for a template-based and agnostic search for an SGWB. Moreover, as a demonstration of the rich potential of SBI, we consider the injection of a population of low signal-to-noise ratio supermassive black hole transient signals into the data. TMNRE can implicitly marginalize over this complicated parameter space, enabling us to directly and accurately reconstruct the stochastic (and instrumental noise) contributions. We publicly release our TMNRE implementation in the form of the code saqqara.

Symmetries and conservation laws associated with the ideal Einstein-Euler system, for stationary and axisymmetric stars, can be utilised to define a set of \emph{flow constants}. These quantities are conserved along flow lines in the sense that their gradients are orthogonal to the four-velocity. They are also conserved along surfaces of constant magnetic flux, making them powerful tools to identify general features of neutron star equilibria. One important corollary of their existence is that mixed poloidal-toroidal fields are inconsistent with the absence of meridional flows except in some singular sense, a surprising but powerful result first proven by Bekenstein and Oron. In this work, we revisit the flow constant formalism to rederive this result together with several new ones concerning both non-linear and perturbative magnetic equilibria. Our investigation is supplemented by some numerical solutions for multipolar magnetic fields on top of a Tolman-VII background, where strict power-counting of the flow constants is used to ensure a self-consistent treatment.

We point out the theoretical possibility of detecting purely gravitationally interacting dark matter using the very high sensitivity of gravitationally mediated quantum phase shift.

A curved space-time is known to act as a source for particle production in environments where gravity plays a significant role. We explore this effect in a minimal setup of cosmic inflation on the production of a scalar dark matter candidate during and after the inflationary stage of the universe. We consider the production of dark matter via direct coupling to the inflaton field or from pure gravitational interactions. Cosmological constraints from structure formation and dark matter isocurvature perturbations are discussed. A new analytical expression for the isocurvature power spectrum is provided.

In this letter, we demonstrate how to use the generalized $\delta N$ formalism, which enables us to compute all the large scale fluctuations, including the gravitational waves, solely by solving the evolution of the background homogeneous Universe. Using the Noether charge density, we derive an analytic formula which describes the mapping between the fluctuations at the horizon crossing and the sourced gravitational waves at the end of inflation. This formula can apply also to an inflation model with an anisotropic background.