Papers for Friday, Apr 05 2024

This STScI Working Group (WG) was charged with soliciting community feedback and evaluating the strategic planning for exoplanet science with JWST and HST given the high quality of exoplanet observations, the significantly lengthened mission lifetime for JWST, and the pronounced expansion of the field over the last decade. We were charged with identifying key science themes, providing recommendations on issues associated with optimal timing and scale of resources, as well as providing a recommended DDT concept achievable with 500 hours of JWST time. We recommend a DDT concept to survey the atmospheres of rocky-M dwarf exoplanets. It is critical to quickly survey a wide sample of such targets to ascertain if they indeed host significant atmospheres, i.e., define the cosmic shoreline, and to identify high priority targets for future follow-up. It is important for this effort to occur early in the mission lifetime. In the context of strategic planning of exoplanet observations, it is useful to estimate the expected exoplanet observational commitment over JWST's lifetime. Given the current usage associated with exoplanets, extended over 20 cycles, it is anticipated that JWST will dedicate $\approx$30,000 hours to exoplanet observations. We recommend efforts to support GO-driven programs that will contribute to this unprecedented data product of JWST. Of the $\approx$30,000 hours of anticipated JWST full-mission time dedicated to exoplanets, we expect that 1/3 of it could, and perhaps inevitably would, form a comprehensive, high S/N, panchromatic, 10$^4$ hour atmospheric survey of planets. Such an observational sample would be a legacy archive that would address a broad range of science questions across various populations of planets. It would also bridge the direct imaging and transit communities and involve a multitude of techniques to detect and characterize exoplanets.

In the cold dark matter paradigm, our Galaxy is predicted to contain >10000 dark matter subhaloes in the $10^5-10^8M_\odot$ range which should be completely devoid of stars. Stellar streams are sensitive to the presence of these subhaloes, which can create small-scale features in streams if they pass closely enough. Modelling these encounters can therefore, potentially recover the subhalo's properties. In this work, we demonstrate this for streams generated in numerical simulations, modelled on eccentric orbits in a realistic Milky Way potential, which includes the Large Magellanic Cloud and the subhalo itself. We focus on a mock model of the ATLAS-Aliqa Uma stream and inject a $10^7 M_\odot$ subhalo, creating a similar discontinuous morphology to current observations. We then explore how well subhalo properties are recovered using mock stream observations, consisting of no observational errors, as well as assuming realistic observational setups. These setups include present day style observations, and what will be possible with 4MOST and Gaia DR5 in the future. We show that we can recover all parameters describing the impact even with uncertainties matching existing data, including subhalo positions, velocities, mass and scale radius. Modelling the subhalo on an orbit instead of assuming an impulse approximation, we greatly reduce the degeneracy between subhalo mass and velocity seen in previous works. However, we find a slight bias in the subhalo mass (~0.1 dex). This demonstrates that we should be able to reliably extract the properties of subhaloes with stellar streams in the near future.

Protostellar discs are mostly modelled as circular structures of gas and dust orbiting a protostar. However, a number of physical mechanisms, e.g. the presence of a (sub)stellar companion or initial axial asymmetry, can cause the gas and dust orbital motion to become eccentric. Theoretical studies have revealed that, when present, disc eccentricity is expected to occur with predictable profiles that can be long-lasting and potentially observable in protostellar systems. We construct an analytical model predicting the typical features of the kinematics and morphology of eccentric protostellar discs, with the final goal of characterising the observational appearance of eccentricity in discs. We validate the model using a numerical simulation of a circumbinary disc (where the binary makes the disc eccentric). We finally post-process the simulation with Monte Carlo Radiative Transfer to study how eccentric features would appear through the "eyes" of ALMA. Besides the motion of the material on eccentric Keplerian orbits in the disc orbital plane, the most characteristic eccentric feature emerging from the analytical model is strong vertical motion with a typical anti-symmetric pattern (with respect to the disc line of pericentres). A circumbinary disc with a $\approx 40$ au eccentric cavity ($e_{\rm cav}=0.2$), carved by an $a_{\rm bin}=15$ au binary, placed at a distance $d=130$ pc, is expected to host in its upper emission surface vertical oscillations up to $v_{z}\sim 400\, {\rm ms}^{-1}$ close to the cavity edge, i.e. well within ALMA spectral and spatial resolution capabilities. A residual spiral pattern in the vertical velocity $\Delta v_{z}\sim 150\, {\rm ms}^{-1}$ of the simulation cannot be captured by the theoretical model, we speculate it to be possibly linked to the presence of a companion in the system.

In this work, we compare the SFRs and SFHs of AGN and non-AGN galaxies. We explore these aspects across different density fields and over three orders of magnitude in L$_X$. For that purpose, we employ X-ray AGN detected in the XMM-XXL field and construct a galaxy control sample, using sources from the VIPERS catalogue. Our final samples consist of 149 X-ray AGN with $\rm 42<log,[L_{X,2-10keV}(ergs^{-1})]<45$ and 3\,488 non-AGN systems. The sources span a redshift range of $\rm 0.5<z<1.0$ and $10.5<\rm log\,[M_*(M_\odot)]<11.5$. For these systems, there are available measurements for their local densities and their spectral lines from the VIPERS catalogue. To compare the SFR of these two populations, we calculate the SFR$_{norm}$ parameter. The latter is defined as the ratio of the SFR of AGN to the SFR of non-AGN galaxies with similar M$_*$ and redshift. Our findings reveal that low and moderate L$_X$ AGN that live in low density fields have a nearly flat SFR$_{norm}-$L$_X$ relation. In contrast, AGN of similar L$_X$ that live in high density environments present an increase of SFR$_{norm}$ with L$_X$. Notably, our results suggest that the most luminous of the AGN exhibit increased SFR relative to non-AGN galaxies, and this trend appears to be independent of the density of the environment. Furthermore, for AGN with similar L$_X$, those in high-density regions tend to have higher SFR$_{norm}$ values compared to their counterparts in low-density areas. Comparison of the D$_n$4000 spectral index, which serves as a proxy for the age of the stellar population, reveals that low-to-moderate L$_X$ AGN live in galaxies with comparable stellar populations with non-AGN systems, regardless of the density field they live in. However, the most luminous X-ray sources tend to live in galaxies that have younger stellar populations than non-AGN galaxies, regardless of the galaxy's environment.

Theoretical models predict that z~6 quasars are hosted in the most massive halos of the underlying dark matter distribution and thus would be immersed in protoclusters of galaxies. However, observations report inconclusive results. We investigate the 1.1 pMpc2 environment of the z = 7.54 luminous quasar ULAS J1342+0928. We search for Lyman-break galaxy candidates (LBG) using deep imaging from the Hubble Space Telescope (HST) in the ACS/F814W, WFC3/F105W/F125W bands, and Spitzer/IRAC at 3.6 $\mu$m and 4.5 $\mu$m. We report a zphot = $7.69^{+0.33}_{-0.23}$ LBG with magF125W = 26.41 at 223 projected-pkpc from the quasar. We find no HST counterpart to one [CII]-emitter previously found with ALMA at 27 projected-pkpc and $z[CII]=7.5341\pm0.0009$ (Venemans et al. 2020). We estimate the completeness of our LBG candidates using results from CANDELS/GOODS deep blank field searches sharing a similar filter setup. We find that >50% of the z~7.5 LBGs with magF125W >25.5 are missed due to the absence of a filter redward of the Lyman-break in F105W, hindering the UV color accuracy of the candidates. We conduct a QSO-LBG clustering analysis revealing a low LBG excess of $0.46^{+1.52}_{-0.08}$ in this quasar field, consistent with an average or low-density field. Consequently, this result does not present strong evidence of an LBG overdensity around ULAS J1342+0928. Furthermore, we identify two LBG candidates with a zphot matching a confirmed z=6.84 absorber along the line-of-sight to the quasar. All these galaxy candidates are excellent targets for follow-up observations with JWST and/or ALMA to confirm their redshift and physical properties.

Galaxies are complex objects, yet the number of independent parameters to describe them remains unknown. We present here a non-parametric method to estimate the intrinsic dimensionality of large datasets. We apply it to wide-band photometric data drawn from the COSMOS2020 catalogue and a comparable mock catalogue from the Horizon-AGN simulation. Our galaxy catalogues are limited in signal-to-noise ratio in all optical and NIR bands. Our results reveal that most of the variance in the wide-band photometry of this galaxy sample can be described with at most $4.3\pm0.5$ independent parameters for star-forming galaxies and $2.9\pm0.2$ for passive ones, both in the observed and simulated catalogues. We identify one of these parameters to be noise-driven, and recover that stellar mass and redshift are two key independent parameters driving the magnitudes. Our findings support the idea that wide-band photometry does not provide more than one additional independent parameter for star-forming galaxies. Although our sample is not mass-limited and may miss some passive galaxies due to our cut in SNR, our work suggests that dimensionality reduction techniques may be effectively used to explore and analyse wide-band photometric data, provided the used latent space is at least four-dimensional.

Multi-planet systems exhibit remarkable architectural diversity. However, short-period giant planets are typically isolated. Compact systems like TOI-5398, with an outer close-orbit giant and an inner small-size planet, are rare among systems containing short-period giants. TOI-5398's unusual architecture coupled with its young age (650 $\pm$ 150 Myr) make it a promising system for measuring the original obliquity between the orbital axis of the giant and the stellar spin axis in order to gain insight into its formation and orbital migration. We collected in-transit (plus suitable off-transit) observations of TOI-5398 b with HARPS-N at TNG on March 25, 2023, obtaining high-precision radial velocity time series that allowed us to measure the Rossiter-McLaughlin (RM) effect. By modelling the RM effect, we obtained a sky-projected obliquity of $\lambda = 3.0^{+6.8}_{-4.2}$ deg for TOI-5398 b, consistent with the planet being aligned. With knowledge of the stellar rotation period, we estimated the true 3D obliquity, finding $\psi = (13.2\pm8.2)$ deg. Based on theoretical considerations, the orientation we measure is unaffected by tidal effects, offering a direct diagnostic for understanding the formation path of this planetary system. The orbital characteristics of TOI-5398, with its compact architecture, eccentricity consistent with circular orbits, and hints of orbital alignment, appear more compatible with the disc-driven migration scenario. TOI-5398, with its relative youth (compared with similar compact systems) and exceptional suitability for transmission spectroscopy studies, presents an outstanding opportunity to establish a benchmark for exploring the disc-driven migration model.

We present the results from a complex study of an eclipsing O-type binary (Aa+Ab) with the orbital period $P_{A}=3.2254367$ days, that forms part of a higher-order multiple system in a configuration (A+B)+C. We derived masses of the Aa+Ab binary $M_{1}= 19.02 \pm 0.12 \,M_\odot$, $M_{2}= 17.50 \pm 0.13 \,M_\odot$, radii $R_{1}= 7.70 \pm 0.05 \,R_\odot$, $R_{2}= 6.64 \pm 0.06 \,R_\odot$, and temperatures $T_1 = 34250 \pm 500 $ K, $T_2 = 33750 \pm 500 $ K. From the analysis of radial velocities, we found a spectroscopic orbit of A in the outer A+B system with $P_{A+B}=195.8$ days ($P_{A+B}/P_{A}\approx 61$). In the O-C analysis, we confirmed this orbit and found another component orbiting the A+B system with $P_{AB+C}=2550$ days ($P_{AB+C}\,/P_{A+B}\approx 13$). From the total mass of the inner binary and its outer orbit, we estimated the mass of the third object, $M_B \gtrsim 10.7 M_\odot$. From the light-travel time effect fit to the O-C data, we obtained the limit for the mass of the fourth component, $M_C \gtrsim 7.3 M_\odot$. These extra components contribute to about 20% to 30% (increasing with wavelength) of the total system light. From the comparison of model spectra with the multiband photometry, we derived a distance modulus of 18.59 $\pm$ 0.06 mag, a reddening of 0.16 $\pm$ 0.02 mag, and an $R_V$ of $3.2$. This work is part of our ongoing project, which aims to calibrate the surface brightness-color relation for early-type stars.

Long-period transiting exoplanets bridge the gap between the bulk of transit- and Doppler-based exoplanet discoveries, providing key insights into the formation and evolution of planetary systems. The wider separation between these planets and their host stars results in the exoplanets typically experiencing less radiation from their host stars; hence, they should maintain more of their original atmospheres, which can be probed during transit via transmission spectroscopy. Although the known population of long-period transiting exoplanets is relatively sparse, surveys performed by the Transiting Exoplanet Survey Satellite (TESS) and the Next Generation Transit Survey (NGTS) are now discovering new exoplanets to fill in this crucial region of the exoplanetary parameter space. This study presents the detection and characterisation of NGTS-30 b/TOI-4862 b, a new long-period transiting exoplanet detected by following up on a single-transit candidate found in the TESS mission. Through monitoring using a combination of photometric instruments (TESS, NGTS, and EulerCam) and spectroscopic instruments (CORALIE, FEROS, HARPS, and PFS), NGTS-30 b/TOI-4862 b was found to be a long-period (P = 98.29838 day) Jupiter-sized (0.928 RJ; 0.960 MJ) planet transiting a 1.1 Gyr old G-type star. With a moderate eccentricity of 0.294, its equilibrium temperature could be expected to vary from 274 K to 500 K over the course of its orbit. Through interior modelling, NGTS-30 b/TOI-4862 b was found to have a heavy element mass fraction of 0.23 and a heavy element enrichment (Zp/Z_star) of 20, making it metal-enriched compared to its host star. NGTS-30 b/TOI-4862 b is one of the youngest well-characterised long-period exoplanets found to date and will therefore be important in the quest to understanding the formation and evolution of exoplanets across the full range of orbital separations and ages.

Here we present speckle observations of 16 low-separation ($s<30$ AU) high probability candidate binaries from the catalog by Medan et al., where secondaries typically lack astrometric solutions in Gaia. From these speckle observations, we find a second component is always detected within the field of view. To determine if the detection is consistent with a physical companion or a chance alignment with a background source, we utilize a statistic from Tokovinin & Kiyaeva that compares the apparent motion of the systems to the expected orbital motion ($\mu^\prime$). Using simulated binary orbits, we construct likelihood distributions of $\mu^\prime$ assuming various total errors on the measurements. With the hypothesis that the system is a true binary, we show that large measurement errors can result in $\mu^\prime$ values higher than expected for bound systems. Using simulated chance alignments, we also create similar likelihoods to test this alternative hypothesis. By combining likelihoods of both true binaries and chance alignments, we find that 15 of the 16 candidates are physical systems regardless of the level of measurement error. Our findings also accommodate all 16 as physical systems if the average, relative measurement error on the binary separations and position angles is $\sim4.3\%$, which is consistent with our knowledge of the Gaia and Gemini speckle pipelines. Importantly, beyond assessing the likelihood of a true binary vs. chance alignment, this quantitative assessment of the true average measurement error will allow more robust error estimates of mass determinations from short separation binaries with Gaia and/or Gemini speckle data.

Optical turbulence modelling and simulation are crucial for developing astronomical ground-based instruments, laser communication, laser metrology, or any application where light propagates through a turbulent medium. In the context of spectrum-based optical turbulence Monte-Carlo simulations, we present an alternative approach to the methods based on the Fast Fourier Transform (FFT) using a quasi-random frequency sampling heuristic. This approach provides complete control over the spectral information expressed in the simulated measurable, without the drawbacks encountered with FFT-based methods such as high-frequency aliasing, low-frequency under-sampling, and static sampling statistics. The method's heuristics, implementation, and an application example from the study of differential piston fluctuations are discussed.

We present the DESI 2024 galaxy and quasar baryon acoustic oscillations (BAO) measurements using over 5.7 million unique galaxy and quasar redshifts in the range 0.1<z<2.1. Divided by tracer type, we utilize 300,017 galaxies from the magnitude-limited Bright Galaxy Survey with 0.1<z<0.4, 2,138,600 Luminous Red Galaxies with 0.4<z<1.1, 2,432,022 Emission Line Galaxies with 0.8<z<1.6, and 856,652 quasars with 0.8<z<2.1, over a ~7,500 square degree footprint. The analysis was blinded at the catalog-level to avoid confirmation bias. All fiducial choices of the BAO fitting and reconstruction methodology, as well as the size of the systematic errors, were determined on the basis of the tests with mock catalogs and the blinded data catalogs. We present several improvements to the BAO analysis pipeline, including enhancing the BAO fitting and reconstruction methods in a more physically-motivated direction, and also present results using combinations of tracers. We present a re-analysis of SDSS BOSS and eBOSS results applying the improved DESI methodology and find scatter consistent with the level of the quoted SDSS theoretical systematic uncertainties. With the total effective survey volume of ~ 18 Gpc$^3$, the combined precision of the BAO measurements across the six different redshift bins is ~0.52%, marking a 1.2-fold improvement over the previous state-of-the-art results using only first-year data. We detect the BAO in all of these six redshift bins. The highest significance of BAO detection is $9.1\sigma$ at the effective redshift of 0.93, with a constraint of 0.86% placed on the BAO scale. We find our measurements are systematically larger than the prediction of Planck-2018 LCDM model at z<0.8. We translate the results into transverse comoving distance and radial Hubble distance measurements, which are used to constrain cosmological models in our companion paper [abridged].

We present cosmological results from the measurement of baryon acoustic oscillations (BAO) in galaxy, quasar and Lyman-$\alpha$ forest tracers from the first year of observations from the Dark Energy Spectroscopic Instrument (DESI), to be released in the DESI Data Release 1. DESI BAO provide robust measurements of the transverse comoving distance and Hubble rate, or their combination, relative to the sound horizon, in seven redshift bins from over 6 million extragalactic objects in the redshift range $0.1<z<4.2$. DESI BAO data alone are consistent with the standard flat $\Lambda$CDM cosmological model with a matter density $\Omega_\mathrm{m}=0.295\pm 0.015$. Paired with a BBN prior and the robustly measured acoustic angular scale from the CMB, DESI requires $H_0=(68.52\pm0.62)$ km/s/Mpc. In conjunction with CMB anisotropies from Planck and CMB lensing data from Planck and ACT, we find $\Omega_\mathrm{m}=0.307\pm 0.005$ and $H_0=(67.97\pm0.38)$ km/s/Mpc. Extending the baseline model with a constant dark energy equation of state parameter $w$, DESI BAO alone require $w=-0.99^{+0.15}_{-0.13}$. In models with a time-varying dark energy equation of state parametrized by $w_0$ and $w_a$, combinations of DESI with CMB or with SN~Ia individually prefer $w_0>-1$ and $w_a<0$. This preference is 2.6$\sigma$ for the DESI+CMB combination, and persists or grows when SN~Ia are added in, giving results discrepant with the $\Lambda$CDM model at the $2.5\sigma$, $3.5\sigma$ or $3.9\sigma$ levels for the addition of Pantheon+, Union3, or DES-SN5YR datasets respectively. For the flat $\Lambda$CDM model with the sum of neutrino mass $\sum m_\nu$ free, combining the DESI and CMB data yields an upper limit $\sum m_\nu < 0.072$ $(0.113)$ eV at 95% confidence for a $\sum m_\nu>0$ $(\sum m_\nu>0.059)$ eV prior. These neutrino-mass constraints are substantially relaxed in models beyond $\Lambda$CDM. [Abridged.]

The first year of data from the Dark Energy Spectroscopic Instrument (DESI) contains the largest set of Lyman-$\alpha$ (Ly$\alpha$) forest spectra ever observed. This data, collected in the DESI Data Release 1 (DR1) sample, has been used to measure the Baryon Acoustic Oscillation (BAO) feature at redshift $z=2.33$. In this work, we use a set of 150 synthetic realizations of DESI DR1 to validate the DESI 2024 Ly$\alpha$ forest BAO measurement. The synthetic data sets are based on Gaussian random fields using the log-normal approximation. We produce realistic synthetic DESI spectra that include all major contaminants affecting the Ly$\alpha$ forest. The synthetic data sets span a redshift range $1.8<z<3.8$, and are analysed using the same framework and pipeline used for the DESI 2024 Ly$\alpha$ forest BAO measurement. To measure BAO, we use both the Ly$\alpha$ auto-correlation and its cross-correlation with quasar positions. We use the mean of correlation functions from the set of DESI DR1 realizations to show that our model is able to recover unbiased measurements of the BAO position. We also fit each mock individually and study the population of BAO fits in order to validate BAO uncertainties and test our method for estimating the covariance matrix of the Ly$\alpha$ forest correlation functions. Finally, we discuss the implications of our results and identify the needs for the next generation of Ly$\alpha$ forest synthetic data sets, with the top priority being to simulate the effect of BAO broadening due to non-linear evolution.

We present an optimized way of producing the fast semi-analytical covariance matrices for the Legendre moments of the two-point correlation function, taking into account survey geometry and mimicking the non-Gaussian effects. We validate the approach on simulated (mock) catalogs for different galaxy types, representative of the Dark Energy Spectroscopic Instrument (DESI) Data Release 1, used in 2024 analyses. We find only a few percent differences between the mock sample covariance matrix and our results, which can be expected given the approximate nature of the mocks, although we do identify discrepancies between the shot-noise properties of the DESI fiber assignment algorithm and the faster approximation used in the mocks. Importantly, we find a close agreement (<~ 5% relative differences) in the projected errorbars for distance scale parameters for the baryon acoustic oscillation measurements. This confirms our method as an attractive alternative to simulation-based covariance matrices, especially for non-standard models or galaxy sample selections, in particular, relevant to the broad current and future analyses of DESI data.

JWST has created a new era of terrestrial exoplanet atmospheric characterization, and with it the possibility to detect potential biosignature gases like CH$_{4}$. Our interpretation of exoplanet atmospheric spectra, and the veracity of these interpretations, will be limited by our understanding of atmospheric processes and the accuracy of input modeling data. Molecular cross-sections are essential inputs to these models. The photochemistry of temperate planets depends on photolysis reactions whose rates are governed by the dissociation cross-sections of key molecules. H$_{2}$O is one such molecule; the photolysis of H$_{2}$O produces OH, a highly reactive and efficient sink for atmospheric trace gases. We investigate the photochemical effects of improved H$_{2}$O cross-sections on anoxic terrestrial planets as a function of host star spectral type (FGKM) and CH$_{4}$ surface flux. Our results show that updated H$_{2}$O cross-sections, extended to wavelengths $>$200 nm, substantially impact the predicted abundances of trace gases destroyed by OH. The differences for anoxic terrestrial planets orbiting Sun-like host stars are greatest, showing changes of up to three orders of magnitude in surface CO levels, and over an order of magnitude in surface CH$_{4}$ levels. These differences lead to observable changes in simulated planetary spectra, especially important in the context of future direct-imaging missions. In contrast, the atmospheres of planets orbiting M-dwarf stars are substantially less affected. Our results demonstrate a pressing need for refined dissociation cross-section data for H$_{2}$O, where uncertainties remain, and other key molecules, especially at mid-UV wavelengths $>$200 nm.

The soft excess in active galactic nuclei (AGNs) may arise through a combination of relativistic reflection and the effects of a warm corona at the surface of the accretion disc. Detailed examination of the soft excess can therefore constrain models of the transport and dissipation of accretion energy. Here, we analyze 34 XMM-Newton observations from 14 Type I AGNs with the reXcor spectral model which self-consistently combines emission from a warm corona with relativistic reflection assuming a lamppost corona. The model divides accretion energy between the disc, the warm corona, and the lamppost. The XMM-Newton observations span a factor of 188 in Eddington ratio ($\lambda_{\mathrm{obs}}$) and 350 in black hole mass, and we find that a warm corona is a significant contributor to the soft excess for 13 of the 14 AGNs with a mean warm corona heating fraction of $0.51$. The reXcor fits reveal that the fraction of accretion energy dissipated in the lamppost is anti-correlated with $\lambda_{\mathrm{obs}}$. In contrast, the relationship between $\lambda_{\mathrm{obs}}$ and both the optical depth and heating fraction of the warm corona appears to transition from an anti-correlation to a correlation at $\lambda_{\mathrm{obs,t}} \approx 0.15$. Therefore, at least one other physical process in addition to the accretion rate is needed to explain the evolution of the warm corona. Overall, we find that a warm corona appears to be a crucial depository of accretion energy in AGNs across a broad range of $\lambda_{\mathrm{obs}}$ and black hole mass.

Some electromagnetically observed ultra-compact binaries will be strong gravitational wave sources for space-based detectors like the Laser Interferometer Space Antenna (LISA). These sources have historically been referred to as "verification binaries" under the assumption that they will be exploited to assess mission performance. This paper quantitatively interrogates that scenario by considering targeted analyses of known galactic sources in the context of a full simulation of the galactic gravitational wave foreground. We find that the analysis of the best currently known LISA binaries, even making maximal use of the available information about the sources, is susceptible to ambiguity or biases when not simultaneously fitting to the rest of the galactic population. While galactic binaries discovered electromagnetically in advance of, or during, the LISA survey are highly valuable multimessenger systems, the need for a global treatment of the galactic gravitational wave foreground calls into question whether or not they are the best sources for data characterization.

An outstanding gap in our knowledge of the solar wind is the relationship between switchbacks and solar wind turbulence. Switchbacks are large fluctuations, even reversals, of the background magnetic field embedded in the solar wind flow. It has been proposed that switchbacks may form as a product of turbulence and decay via coupling with the turbulent cascade. In this work, we examine how properties of solar wind magnetic field turbulence vary in the presence or absence of switchbacks. Specifically, we use in-situ particle and fields measurements from Parker Solar Probe to measure magnetic field turbulent wave power, separately in the inertial and kinetic ranges, as a function of switchback magnetic deflection angle. We demonstrate that the angle between the background magnetic field and the solar wind velocity in the spacecraft frame ($\theta_{vB}$) strongly determines whether Parker Solar Probe samples wave power parallel or perpendicular to the background magnetic field. Further, we show that $\theta_{vB}$ is strongly modulated by the switchback magnetic deflection angle. In this analysis, we demonstrate that switchback deflection angle does not correspond to any significant increase in wave power in either the inertial range or at kinetic scales. This result implies that switchbacks do not strongly couple to the turbulent cascade in the inertial or kinetic ranges via turbulent wave-particle interactions. Therefore, we do not expect switchbacks to contribute significantly to solar wind heating through this type of energy conversion pathway, although contributions via other mechanisms, such as magnetic reconnection may still be significant.

Mid-infrared emission features probe the properties of ionized gas, and hot or warm molecular gas. The Orion Bar is a frequently studied photodissociation region (PDR) containing large amounts of gas under these conditions, and was observed with the MIRI IFU aboard JWST as part of the "PDRs4All" program. The resulting IR spectroscopic images of high angular resolution (0.2") reveal a rich observational inventory of mid-IR emission lines, and spatially resolve the substructure of the PDR, with a mosaic cutting perpendicularly across the ionization front and three dissociation fronts. We extracted five spectra that represent the ionized, atomic, and molecular gas layers, and measured the most prominent gas emission lines. An initial analysis summarizes the physical conditions of the gas and the potential of these data. We identified around 100 lines, report an additional 18 lines that remain unidentified, and measured the line intensities and central wavelengths. The H I recombination lines originating from the ionized gas layer bordering the PDR, have intensity ratios that are well matched by emissivity coefficients from H recombination theory, but deviate up to 10% due contamination by He I lines. We report the observed emission lines of various ionization stages of Ne, P, S, Cl, Ar, Fe, and Ni, and show how certain line ratios vary between the five regions. We observe the pure-rotational H$_2$ lines in the vibrational ground state from 0-0 S(1) to 0-0 S(8), and in the first vibrationally excited state from 1-1 S(5) to 1-1 S(9). We derive H$_2$ excitation diagrams, and approximate the excitation with one thermal (~700 K) component representative of an average gas temperature, and one non-thermal component (~2700 K) probing the effect of UV pumping. We compare these results to an existing model for the Orion Bar PDR and highlight the differences with the observations.

H-poor superluminous supernovae (SLSNe-I) are characterized by O II lines around 4,000 - 4,500 A in pre-/near-maximum spectra, so-called W-shaped O II lines. As these lines are from relatively high excitation levels, they have been considered a sign of non-thermal processes, which may give a hint of power sources of SLSNe-I. However, the conditions for these lines to appear have not been understood well. In this work, we systematically calculate synthetic spectra to reproduce observed spectra of eight SLSNe-I, parameterizing departure coefficients from the nebular approximation in the SN ejecta (expressed as b_neb). We find that most of the observed spectra can be reproduced well with b_neb ~< 10, which means that no strong departure is necessary for the formation of the W-shaped O II lines. We also show that the appearance of the W-shaped O II lines is sensitive to temperature; only spectra with temperatures T ~ 14,000 - 16,000 K can produce the W-shaped O II lines without large departures. Based on this, we constrain the non-thermal ionization rate near the photosphere. Our results suggest that spectral features of SLSNe-I can give independent constraints on the power source through the non-thermal ionization rates.

We present strong constraints on the spacetime variation of the fine-structure constant $\alpha$ using the Dark Energy Spectroscopic Instrument (DESI). In this pilot work, we utilize $\sim110,000$ galaxies with strong and narrow O III $\lambda\lambda$4959,5007 emission lines to measure the relative variation $\Delta\alpha/\alpha$ in space and time. The O III doublet is arguably the best choice for this purpose owing to its wide wavelength separation between the two lines and its strong emission in many galaxies. Our galaxy sample spans a redshift range of $0<z<0.95$, covering half of all cosmic time. We divide the sample into subsamples in 10 redshift bins ($\Delta z=0.1$), and calculate $\Delta\alpha/\alpha$ for the individual subsamples. The uncertainties of the measured $\Delta\alpha/\alpha$ are roughly between $2\times10^{-6}$ and $2\times10^{-5}$. We find an apparent $\alpha$ variation with redshift at a level of $\Delta\alpha/\alpha=(2\sim3)\times10^{-5}$. This is highly likely to be caused by systematics associated with wavelength calibration, since such small systematics can be caused by a wavelength distortion of $0.002-0.003$ \AA, which is beyond the accuracy that the current DESI data can achieve. We refine the wavelength calibration using sky lines for a small fraction of the galaxies, but it does not change our main results. We further probe the spatial variation of $\alpha$ in small redshift ranges, and do not find obvious, large-scale structures in the spatial distribution of $\Delta\alpha/\alpha$. As DESI is ongoing, we will include more galaxies, and by improving the wavelength calibration, we expect to obtain a better constraint that is comparable to the strongest current constraint.

Photometric wide-area observations in the next decade will be capable of detecting a large number of galaxy-scale strong gravitational lenses, increasing the gravitational lens sample size by orders of magnitude. To aid in forecasting and analysis of these surveys, we construct a flexible model based on observed distributions for the lens and source properties and test it on the results of past lens searches, including SL2S, SuGOHI and searches on the COSMOS HST and DES fields. We use this model to estimate the expected yields of some current and planned surveys, including Euclid Wide, Vera Rubin LSST, and Roman High Latitude Wide Area. The model proposed includes a set of free parameters to constrain on the identifiability of a lens in an image, allowing construction of prior probability distributions for different lens detection methods. The code used in this work is made publicly available.

The detection of GW170817/AT2017gfo inaugurated an era of multimessenger astrophysics, in which gravitational wave and multiwavelength photon observations complement one another to provide unique insight on astrophysical systems. A broad theoretical consensus exists in which the photon phenomenology of neutron star mergers largely rests upon the evolution of the small amount of matter left on bound orbits around the black hole or massive neutron star remaining after the merger. Because this accretion disk is far from inflow equilibrium, its subsequent evolution depends very strongly on its initial state, yet very little is known about how this state is determined. Using both snapshot and tracer particle data from a numerical relativity/MHD simulation of an equal-mass neutron star merger that collapses to a black hole, we show how gravitational forces arising in a non-axisymmetric, dynamical spacetime supplement hydrodynamical effects in shaping the initial structure of the bound debris disk. The work done by hydrodynamical forces is ${\sim}10$ times greater than that due to time-dependent gravity. Although gravitational torques prior to remnant relaxation are an order of magnitude larger than hydrodynamical torques, their intrinsic sign symmetry leads to strong cancellation; as a result, hydrodynamical and gravitational torques have comparable effect. We also show that the debris disk's initial specific angular momentum distribution is sharply peaked at roughly the specific angular momentum of the merged neutron star's outer layers, a few $r_g c$, and identify the regulating mechanism.

Astrophysically motivated population models for binary black hole observables are often insufficient to capture the imprints of multiple formation channels. This is mainly due to the strongly parametrized nature of such investigations. Using a non-parametric model for the joint population-level distributions of binary black hole component masses and effective inspiral spins, we find hints of multiple subpopulations in the third gravitational wave transient catalog. The higher (more positive) spin subpopulation is found to have a mass-spectrum without any feature at the $30-40M_{\odot}$, which is consistent with the predictions of isolated stellar binary evolution, simulations for which place the pile up due to pulsational pair-instability supernovae near $50M_{\odot}$ or higher. The other sub-population with effective spins closer to zero shows a feature at $30-40M_{\odot}$ and is consistent with binary black holes formed dynamically in globular clusters, which are expected to peak around $30M_{\odot}$. We also compute merger rates for these two subpopulations and find that they are consistent with the theoretical predictions of the corresponding formation channels. We validate our results by checking their robustness against variations of several model configurations and by analyzing large simulated catalogs with the same model.

Gamma-ray bursts (GRBs) are believed to launch relativistic jets, which generate prompt emission by their internal processes and drive external shocks into surrounding medium, accounting for the long-lasting afterglow emission. However, how the jet powers the external shock is an open question. The unprecedented observations of the keV-MeV emission with GECAM and the TeV emission with LHAASO of so far the brightest burst, GRB 221009A, offer a great opportunity to study the prompt-to-afterglow transition and the early dynamical evolution of the external shock. In this letter, we find that the cumulative light curve of keV-MeV emission could well fit the rising stage of the TeV light curve of GRB 221009A, with a time delay of $4.45^{+0.26}_{-0.26}$\,s for TeV emission. Moreover, both the rapid increase in the initial stage and the excess from about \T+260\,s to 270\,s in the TeV light curve could be interpreted by inverse Compton (IC) scatterings of the inner-coming photons by the energetic electrons in external shock. Our results not only reveal a close relation between the keV-MeV and TeV emission, but also indicate a continuous, rather than impulsive, energy injection to the external shock. Assuming an energy injection rate proportional to the keV-MeV flux, we build a continuous energy injection model which well fits the TeV light curve of GRB 221009A, and provides an estimate of the Lorentz factor of the jet.

Flux ropes erupting from the vicinity of the black hole are thought to be a potential model for the flares observed in Sgr\,A$^*$. In this study, we examine the radiative properties of flux ropes that emerged from the vicinity of the black hole. We have performed three-dimensional two-temperature General Relativistic Magnetohydrodynamic (GRMHD) simulations of magnetized accretion flows with alternating multiple magnetic loops, and General Relativistic Radiation Transfer (GRRT) calculations. In GRMHD simulations, two different sizes of initial magnetic loops are implemented. In the small loop case, magnetic dissipation leads to a weaker excitement of magneto-rotational instability inside the torus which generates a lower accretion rate compared to the large loop case. However, it makes more generation of flux ropes due to frequent reconnection by magnetic loops with different polarities. By calculating the thermal synchrotron emission, we found that the variability of light curves and emitting region are tightly related. At $230\,\rm GHz$ and higher frequency, the emission from the flux ropes is relatively stronger compared with the background, which is responsible for the filamentary structure in the images. At lower frequencies, e.g. $43\,\rm GHz$, emission comes from more extended regions, which have a less filamentary structure in the image. Our study shows self-consistent electron temperature models are essential for the calculation of thermal synchrotron radiation and the morphology of the GRRT images. Flux ropes contribute considerable emission at frequencies $\gtrsim 230\,\rm GHz$.

For the last decade, several probes have pointed to a cosmological tension between the amplitude of density fluctuations extrapolated from the cosmic microwave background within the standard cosmological model and the one encapsulated by the $S_8$ parameter from large scale structure. The origin of this $S_8$ tension has not yet been elucidated and may hint at systematics in the data, unaccounted effects from baryonic physics, or new physics beyond the standard model of cosmology. Baryonic physics may in principle provide a nonlinear solution to the tension by suppressing the matter power spectrum more strongly on nonlinear scales than is traditionally assumed. Such a solution would not worsen the Hubble tension, contrary to many other proposed solutions to the $S_8$ tension. However, no realistic baryonic feedback in hydrodynamical simulations provides the needed suppression as a function of redshift. Here, we point out that a scale-dependence of local-type primordial non-Gaussianities (PNG), with significant PNG at scales of a few Mpc, can provide the needed suppression, since such PNG can suppress the power spectrum at slightly larger scales than baryons do. We demonstrate this by devising collisionless numerical simulations of structure formation in boxes of 0.5 Gpc/h with scale-dependent local-type PNG. Our simple models show that, as a proof of principle, scale-dependent PNG, with a Gaussian random field for primodial density fluctuations on large scales and $f_{\rm NL} \simeq -300$ at $\lesssim 10$ Mpc scales, together with state-of-the-art baryonification of the matter power spectrum, can in principle solve the $S_8$ tension. The $S_8$ tension would then be a smoking-gun of non-trivial inflationary physics.

Thanks to the precise astrometric measurements of proper motions by the Gaia mission, a new tidal stellar stream has been discovered in the northern hemisphere. The distribution of star count shows that the stream is approximately $80$ degrees long and $1.70$ degrees wide. Observations of $21$ member stars, including 14 RR Lyrae stars, indicate that the stream has an eccentric and retrograde orbit with $e=0.58$. The low metallicity, high total energy, and large angular momentum suggest that it is associated with the merging event Sequoia. This discovery suggests the possibility of finding more substructures with high eccentricity orbits, even in the inner halo.

MACS J0600.1-2008 (MACS0600) is an X-ray luminous, massive galaxy cluster at $z_{\mathrm{d}}=0.43$, studied previously as part of the REionization LensIng Cluster Survey (RELICS) and ALMA Lensing Cluster Survey (ALCS) projects which revealed a complex, bimodal mass distribution and an intriguing high-redshift object behind it. Here, we report on the results of an extended strong-lensing (SL) analysis of this system. Using new JWST and ground-based Gemini-N and Keck data, we obtain 13 new spectroscopic redshifts of multiply imaged galaxies and identify 12 new photometric multiple-image systems and candidates, including two multiply imaged $z\sim7$ objects. Taking advantage of the larger areal coverage, our analysis reveals a new bimodal, massive SL structure adjacent to the cluster which we measure spectroscopically to lie at the same redshift and whose existence was implied by previous SL-modeling analyses. While based in part on photometric systems identified in ground-based imaging requiring further verification, our extended SL model suggests that the cluster may have the second-largest critical area and effective Einstein radius observed to date, $A_{\mathrm{crit}}\simeq2.16\,\mathrm{arcmin}^2$ and $\theta_{\mathrm{E}}=49.7''\pm5.0''$ for a source at $z_{\mathrm{s}}=2$, enclosing a total mass of $M(<\theta_{\mathrm{E}})=(4.7\pm0.7)\times10^{14}\,\mathrm{M}_{\odot}$. Yet another, probably related massive cluster structure, discovered in X-rays $5'$ (1.7 Mpc) further north, suggests that MACS0600 is in fact part of an even larger filamentary structure. This discovery adds to several recent detections of massive structures around SL galaxy clusters and establishes MACS0600 as a prime target for future high-redshift surveys with JWST.

For many transitions, atomic data, such as the oscillator strength (log(gf)) and the central wavelength of the line, are poorly constrained or even unknown. We present and test a new inversion method that infers atomic line parameters and the height stratification of the atmospheric parameters from spatially resolved spectropolarimetric observations of the Sun. This method is implemented in the new inversion code $\texttt{globin}$. The new method imposes a spatial coupling in inversion parameters common to all pixels, such as the atomic parameters of the observed spectral lines, and infers atmospheric parameters for each spatial pixel individually. The uniqueness of this method lies in its ability to retrieve reliable atomic parameters even for heavily blended spectral lines. We tested the method by applying it to a set of 18 spectral lines between 4015 \r{A} and 4017 \r{A}, synthesized from a 3D magnetohydrodynamic simulation containing a sunspot and the quiet Sun region around it. The results were then compared with a previously used inversion method where atomic parameters were determined for every pixel independently (pixel-by-pixel method). The new method was able to retrieve the log(gf) values of all lines to an accuracy of 0.004 dex, while the pixel-by-pixel method retrieved the same parameter to an accuracy of only 0.025 dex. The largest differences between the two methods are evident for the heavily blended lines, with the former method performing better than the latter. In addition, the new method is also able to infer reliable atmospheric parameters in all the inverted pixels by successfully disentangling the degeneracies between the atomic and atmospheric parameters. The new method is well suited for the reliable determination of both atomic and atmospheric parameters and works well on all spectral lines, including those that are weak and/or severely blended.

Absorption line spectroscopy using background quasars can provide strong constraints on galactic outflows. In this paper, we investigate possible scaling relations between outflow properties, namely outflow velocity \Vout, the mass ejection rate $\dot M_{\rm out}$, and the mass loading factor $\eta$ and the host galaxy properties, such as star formation rate (SFR), SFR surface density, redshift, and stellar mass using galactic outflows probed by background quasars from MEGAFLOW and other surveys. We find that $V_{\rm out}$ ($\eta$) is (anti-)correlated with SFR and SFR surface density. We extend the formalism of momentum-driven outflows of Heckman et al. to show that it applies not only to down the barrel studies but also to winds probed by background quasars, suggesting a possible universal wind formalism. Under this formalism, we find a clear distinction between ``strong'' and ``weak'' outflows where ``strong'' outflows seem to have tighter correlations with galaxy properties (SFR or galaxy stellar mass) than ``weak'' outflows.

HAT-P-67 b is one of the lowest-density gas giants known to date, making it an excellent target for atmospheric characterization through the transmission spectroscopy technique. In the framework of the GAPS large programme, we collected four transit events, with the aim of studying the exoplanet atmosphere and deriving the orbital projected obliquity. We exploited the high-precision GIARPS observing mode of the TNG, along with additional archival TESS photometry, to explore the activity level of the host star. We performed transmission spectroscopy, both in the VIS and in the nIR wavelength range, and analysed the RML effect both fitting the RVs and the Doppler shadow. Based on the TESS photometry, we redetermined the transit parameters of HAT-P-67 b. By modelling the RML effect, we derived a sky-projected obliquity of ($2.2\pm0.4$){\deg} indicating an aligned planetary orbit. The chromospheric activity index $\log\,R^{\prime}_{\rm HK}$, the CCF profile, and the variability in the transmission spectrum of the H$\alpha$ line suggest that the host star shows signatures of stellar activity and/or pulsations. We found no evidence of atomic or molecular species in the VIS transmission spectra, with the exception of pseudo-signals corresponding to Cr I, Fe I, H$\alpha$, Na I, and Ti I. In the nIR range, we found an absorption signal of the He I triplet of 5.56$^{+0.29}_{-0.30}$%(19.0$\sigma$), corresponding to an effective planetary radius of $\sim$3$R_p$ (where $R_p\sim$2$R_J$) which extends beyond the planet's Roche Lobe radius. Owing to the stellar variability, together with the high uncertainty of the model, we could not confirm the planetary origin of the signals found in the optical transmission spectrum. On the other hand, we confirmed previous detections of the infrared He I triplet, providing a 19.0$\sigma$ detection. Our finding indicates that the planet's atmosphere is evaporating.

In the recent second part of the third observation run by the LIGO-Virgo-KAGRA collaboration, a candidate with sub-solar mass components was reported, which we labelled as SSM200308. This study investigates the premise that primordial black holes (PBHs), arising from Gaussian perturbation collapses, could explain SSM200308. Through Bayesian analysis, we obtain the primordial curvature power spectrum that leads to the merger rate of PBHs aligning with observational data as long as they constitute $3.54^{+7.66}_{-2.82}\times 10^{-3}$ of the dark matter. However, while the gravitational wave (GW) background from binary PBH mergers is within current observational limits, the scalar-induced GWs associated with PBH formation exceed the constraints imposed by pulsar timing arrays, challenging the Gaussian perturbation collapse PBH model as the source of SSM200308.

The planet's mass loss is important for the planet's formation and evolution. The radius valley (RV) is believed to be triggered by evaporation-induced mass loss. As an alternative mechanism for the RV, the mass loss of post-impact planets is thoroughly investigated in this work. The impact energy is converted to the planet's internal energy, enhancing its core energy and accelerating mass loss and orbital migration. As the host star changes from K-type to F-type, the planet's mass loss and orbital migration increase. When the initial gas-to-core mass ratio (GCR) is small, the migration efficiency for planets around K-type stars will increase, which helps to suppress mass loss and retain the planet's mass and radius within a specific range. On the contrary, planets around more massive F-type stars experience more substantial mass loss, potentially leading to complete mass loss, and migrate to orbits with longer periods. Our calculation shows that planets around different spectral types of host stars give rise to an RV ranging from 1.3-2.0 $R_{\oplus}$, consistent with the observed range of 1.3-2.6 $R_{\oplus}$. Despite the presence of uncertain parameters, the planetesimal impact can promote the RV establishment for planets around host stars of different spectral types.

We report the first detection in the interstellar medium of $N$-cyanomethanimine (H$_2$CNCN), the stable dimer of HCN of highest energy, and the most complex organic molecule identified in space containing the prebiotically relevant NCN backbone. We have identified a plethora of $a$-type rotational transitions with 3 $\leq J_\text{up} \leq$ 11 and $K_\text{a} \leq$ 2 that belong to this species towards the Galactic Center G+0.693-0.027 molecular cloud, the only interstellar source showing the three cyanomethanimine isomers (including the $Z$- and $E$- isomers of $C$-cyanomethanimine, HNCHCN). We have derived a total column density for H$_2$CNCN of (2.9$\, \pm \,$0.1)$\times$10$^{12}$ cm$^{-2}$, which translates into a total molecular abundance with respect to H$_2$ of (2.1$\, \pm \,$0.3)$\times$10$^{-11}$. We have also revisited the previous detection of $E$- and $Z$-HNCHCN, and found a total $C/N$-cyanomethanimine abundance ratio of 31.8$\, \pm \,$1.8 and a $Z/E$-HNCHCN ratio of 4.5$\, \pm \,$0.2. While the latter can be explained on the basis of thermodynamic equilibrium, chemical kinetics are more likely responsible for the observed $C/N$-cyanomethanimine abundance ratio, where the gas-phase reaction between methanimine (CH$_2$NH) and the cyanogen radical (CN) arises as the primary formation route.

In this article, we introduce an innovative parametric representation of the Hubble parameter, providing a model-independent means to explore the dynamics of an accelerating cosmos. The model's parameters are rigorously constrained through a Markov Chain Monte Carlo (MCMC) approach, leveraging a comprehensive dataset consisting of 31 data points from cosmic chronometers (CC), 1701 updated observations of Pantheon supernovae type Ia (SNeIa), and 6 data points from baryonic acoustic oscillations (BAO). Our analysis delves into the behavior of various cosmological parameters within the model, including the transition from a decelerating phase to an accelerating one, as well as the density parameters and the equation of state (EoS) parameter. The outcomes of our investigation reveal that the equation of state parameter aligns with characteristics reminiscent of the phantom model, supporting the prevailing understanding of our universe's current state of acceleration. This research contributes valuable insights into the ongoing cosmic expansion and underscores the utility of our novel parametric approach.

The abundance of galaxy clusters with mass and redshift is a well-known cosmological probe. The cluster mass is a key parameter for studies that aim to constrain cosmological parameters using galaxy clusters, making it critical to understand and properly account for the errors in its estimates. Subsequently, it becomes important to correctly calibrate scaling relations between observables like the integrated Compton parameter and the mass of the cluster. The NIKA2 Sunyaev-Zeldovich Large program (LPSZ) enables one to map the intracluster medium profiles in the mm-wavelength band with great details (resolution of $11 \ \mathrm{\&}\ 17^{\prime \prime}$ at $1.2 \ \mathrm{\&}\ 2 $ mm, respectively) and hence, to estimate the cluster hydrostatic mass more precisely than previous SZ observations. However, there are certain systematic effects which can only be accounted for with the use of simulations. For this purpose, we employ THE THREE HUNDRED simulations which have been modelled with a range of physics modules to simulate galaxy clusters. The so-called twin samples are constructed by picking synthetic clusters of galaxies with properties close to the observational targets of the LPSZ. In particular, we use the Compton parameter maps and projected total mass maps of these twin samples along 29 different lines of sight. We investigate the scatter that projection induces on the total masses. Eventually, we consider the statistical values along different lines of sight to construct a kind of 3D scaling law between the integrated Compton parameter, total mass, and overdensity of the galaxy clusters to determine the overdensity that is least impacted by the projection effect.

Context. In recent years, the R-Process Alliance (RPA) has conducted a successful search for stars enhanced in elements produced by the rapid neutron-capture (r-)process. In particular, the RPA has uncovered a number of stars strongly enriched in light r-process elements, such as Sr, Y and Zr, the so-called limited-r stars, in order to investigate the astrophysical production site(s) of these elements. Aims. With this paper, we aim to investigate the possible formation sites for light neutron-capture elements, by deriving detailed abundances for neutron-capture elements from high-resolution, high signal-to-noise spectra of three limited-r stars. Methods. We conducted a 1D local thermodynamic equilibrium spectroscopic abundance analysis of three stars, as well as a kinematic analysis. Further, we calculated the lanthanide mass fraction (XLa) of our stars and of limited-r stars from the literature. Results. We found that the neutron-capture element abundance pattern of limited-r stars behaves differently depending on their [Ba/Eu]] ratios, and suggest that this should be taken into account in future investigations of their abundances. Furthermore, we found that the XLa of limited-r stars is lower than that of the kilonova AT2017gfo. The latter seems to be in the transition zone between limited-r XLa and that of r-I, r-II stars. Finally, we found that, unlike r-I and r-II stars, the current sample of limited-r stars are largely born in the Galaxy rather than being accreted.

We present a comprehensive characterization of 26 CME-driven compressive waves known as Coronal Bright Fronts (CBFs) observed in the low solar corona between 2010 and 2017. These CBFs have been found to be associated with SEP events near Earth, indicating their importance in understanding space weather phenomena. The aim of this study is to analyze and describe the early dynamics of CBFs using a physics-based heliospheric SEP forecasting system known as the SPREAdFAST framework. This framework utilizes a chain of data-driven analytic and numerical models to predict SEP fluxes at multiple locations in the inner heliosphere by considering their acceleration at CMEs near the Sun and subsequent interplanetary transport. To estimate the time-dependent plasma and compression parameters of the CBFs, we utilized sequences of base-difference images obtained from the AIA instrument on board the SDO satellite, and measurements of the height-time profiles of the CMEs obtained from the LASCO instrument on board the SOHO satellite. We employed kinematic measurements and plasma model results to derive these parameters. The SPREAdFAST framework facilitated the analysis and correlation of these observations with SEP events near Earth. Our analysis yielded statistical relations and distributions for both the shocks and plasma parameters associated with the 26 CBFs investigated. By combining the observations from the AIA and LASCO instruments, as well as the data products from the SPREAdFAST framework, we obtained a comprehensive understanding of the early dynamics of CBFs, including their temporal evolution, plasma properties, and compressional characteristics. These findings contribute to the growing body of knowledge in the field and have implications for space weather forecasting and the study of SEP events.

Herbig Ae/Be (HAeBe) stars are intermediate-mass pre-main sequence stars, characterized by infrared excess and emission lines. They are observed to emit X-rays, whose origin is a matter of discussion and not settled yet. X-ray emission is not expected in HAeBe stars, as they lack the sub-surface convective zone. In this study, we retrieved observations from the Chandra archive for 62 HAeBe stars, among which 44 sources (detection fraction $\sim$71%) were detected in X-rays, with 7 being new detections. We use this sample as a test bed to conduct a comparative analysis of the X-ray properties of HAeBe stars and their low-mass counterparts, T Tauri Stars (TTSs). Further, we compare the X-ray properties of HAeBe stars and TTSs with optical and IR properties to constrain the X-ray emission mechanism in HAeBe stars. We found no correlation between X-ray emission and disk properties of HAeBe stars, confirming that X-rays are not related to accretion shocks. About 56% of HAeBe stars without any known sub-arcsec companions have lower plasma temperatures (kT $\leq$ 2 keV). We observe flaring/variability in HAeBe stars with confirmed low-mass companions. These stars show plasma temperatures > 2 keV, similar to TTSs. Guided by this information we discuss the role of a T Tauri companion for X-ray emission seen in our sample of HAeBe stars. From the results obtained in this paper, we suggest that X-ray emission from HAeBe stars may not be related to accretion shocks or hidden TTS, but rather can be due to magnetically driven coronal emission.

We present new H- and K-band spectroscopy for the bulge of M31, taken with the LUCI spectrograph at the Large Binocular Telescope (LBT). We studied radial trends of CO absorption features (namely, CO1.58, CO1.60, CO1.64, CO1.66, CO1.68, CO2.30, CO2.32, CO2.35) in the bulge of M31, out to a galactocentric distance of 100'' (380pc). We find that most COs do not exhibit a strong radial gradient, despite the strong metallicity gradient inferred from the optical spectral range, except for CO1.64, showing a steep increase in the center. We compared the observed line strengths to predictions of different state-of-the-art stellar population models, including an updated version of EMILES models, which also uses the extended IRTF spectral library. The observed COs are close to models' predictions, but in some models they turn out to be underestimated. We find that the lack of radial gradients is due to the combination of increasing CO strength with metallicity and C abundance, and decreasing CO strength with IMF slope and O abundance. We speculate that the steep gradient of CO1.64 might be due to Na overabundance. Remarkably, we were able to fit, at the same time, optical indices and all the NIR COs except for CO1.68, leaving abundance ratios (i.e., [C/Fe], [O/Fe], and [Mg/Fe]) as free-fitting parameters, imposing age and metallicity constraints from the optical, with no significant contribution from intermediate-age populations. For the majority of the bulge, we find [Mg/Fe]~0.15dex, [O/Fe] larger than [Mg/Fe] (by ~0.1dex), and C abundance consistent with that of Mg. In the central (few arcsec) region, we still find an enhancement of O and Mg, but significantly lower [C/Fe]. We find that the COs' line strengths of the bulge are significantly lower than those of massive galaxies, possibly because of a difference in carbon abundance, as well as, to some extent, total metallicity.

Millisecond pulsars (MSPs) are believed to be very old neutron stars (NSs) whose age may exceed significantly $10^8$ yrs. Although cooling scenarios of isolated NSs predict for that age a surface temperature $T_s\sim 10^4$ K, observations of the nearest MSP J0437-4715 indicate $T_s$ well above that value. Besides the heating of the polar cap surface by backflowing charged particles, Joule heating in the crust can contribute to the overall heat budget of MSPs. Since the dipolar field component, derived from $P$ and $\dot{P}$ measurements, is much too weak for remarkable heating, smaller-scale structures should be analysed whether they can supply the demanded heat. For this purpose we study the small scale field structure of radio pulsars. Magnetic field components, significantly stronger than the dipolar one, may exist especially at the surface of MSPs. We assign upper limits to the strength of single field components up to a multipolarity of $l=10$ and the corresponding deviations from axial symmetry $m \le l$. Arguments are provided that the decay of the small-scale components with $l=3$ or $l=4$ of the crustal magnetic field may cause the relatively high surface temperature of isolated MSPs.

To enhance our understanding of the impact of galaxy mergers on the kinematics of early-type galaxies (ETGs), we examine differences in specific stellar angular momentum within the half-light radius ($\lambda_{R_e}$) among ETGs with different types of tidal features and those without such features. This is accomplished by categorizing tidal features, which serve as direct evidence of recent mergers, into shells, streams, and tails, through deep images from the DESI Legacy Survey, and by using MaNGA data for the analysis of the kinematics of 1244 ETGs at $z<0.055$. We find that ETGs with tidal features typically have reduced $\lambda_{R_e}$ values that are lower by 0.12 dex than ETGs without tidal features. ETGs with shells contribute most to the reduction in $\lambda_{R_e}$. Consequently, nearly half of ETGs with shells are classified as slow rotators, a fraction that is more than twice as high as that of ETGs with tails or streams, and over three times higher than that of ETGs without tidal features. These trends generally remain valid even when ETGs are divided into several mass bins. Our findings support the idea that radial mergers, which are more effective at reducing $\lambda_{R_e}$ than circular mergers, are more closely associated with the formation of shells rather than streams or tails. The detection of shells in slightly more massive ETGs compared to streams and tails may be attributed to the fact that massive satellite galaxies are more likely to be accreted through radial orbits, due to the nature of dynamical friction.

Hyperspectral images are data cubes with two spatial dimensions and a third spectral dimension, providing a spectrum for each pixel, and thus allow the mapping of extended sources' physical properties. In this article, we present the Semi-blind Unmixing with Sparsity for Hyperspectral Images (SUSHI), an algorithm for non-stationary unmixing of hyperspectral images with spatial regularization of spectral parameters. The method allows for the disentangling of physical components without the assumption of a unique spectrum for each component. Thus, unlike most source separation methods used in astrophysics, all physical components obtained by SUSHI vary in spectral shape and in amplitude across the data cube. Non-stationary source separation is an ill-posed inverse problem that needs to be constrained. We achieve this by training a spectral model and applying a spatial regularization constraint on its parameters. For the spectral model, we used an Interpolatory Auto-Encoder, a generative model that can be trained with limited samples. For spatial regularization, we applied a sparsity constraint on the wavelet transform of the model parameter maps. We applied SUSHI to a toy model meant to resemble supernova remnants in X-ray astrophysics, though the method may be used on any extended source with any hyperspectral instrument. We compared this result to the one obtained by a classic 1D fit on each individual pixel. We find that SUSHI obtains more accurate results, particularly when it comes to reconstructing physical parameters. We applied SUSHI to real X-ray data from the supernova remnant Cassiopeia A and to the Crab Nebula. The results obtained are realistic and in accordance with past findings but have a much better spatial resolution. Thanks to spatial regularization, SUSHI can obtain reliable physical parameters at fine scales that are out of reach for pixel-by-pixel methods.

More than a century ago, Albert Einstein presented his general theory of gravitation (GR) to the Prussian Academy of Sciences. One of the predictions of the theory is that not only particles and objects with mass, but also the quanta of light, photons, are tied to the curvature of space-time, and thus to gravity. There must be a critical mass density, above which photons cannot escape. These are black holes (henceforth BH). It took fifty years after the theory was announced before possible candidate objects were identified by observational astronomy. And another fifty years have passed, until we finally have in hand detailed and credible experimental evidence that BHs of 10 to 10^10 times the mass of the Sun exist in the Universe. Three very different experimental techniques, but all based on Michelson interferometry or Fourier-inversion spatial interferometry have enabled the critical experimental breakthroughs. It has now become possible to investigate the space-time structure in the vicinity of the event horizons of BHs. We briefly summarize these interferometric techniques, and discuss the spectacular recent improvements achieved with all three techniques. Finally, we sketch where the path of exploration and inquiry may go on in the next decades.

The flux of ultra-high energy cosmic rays reaching Earth above the ankle energy (5 EeV) can be described as a mixture of nuclei injected by extragalactic sources with very hard spectra and a low rigidity cutoff. Extragalactic magnetic fields existing between the Earth and the closest sources can affect the observed CR spectrum by reducing the flux of low-rigidity particles reaching Earth. We perform a combined fit of the spectrum and distributions of depth of shower maximum measured with the Pierre Auger Observatory including the effect of this magnetic horizon in the propagation of UHECRs in the intergalactic space. We find that, within a specific range of the various experimental and phenomenological systematics, the magnetic horizon effect can be relevant for turbulent magnetic field strengths in the local neighbourhood of order $B_{\rm rms}\simeq (50-100)\,{\rm nG}\,(20\rm{Mpc}/{d_{\rm s})( 100\,\rm{kpc}/L_{\rm coh}})^{1/2}$, with $d_{\rm s}$ the typical intersource separation and $L_{\rm coh}$ the magnetic field coherence length. When this is the case, the inferred slope of the source spectrum becomes softer and can be closer to the expectations of diffusive shock acceleration, i.e., $\propto E^{-2}$. An additional cosmic-ray population with higher source density and softer spectra, presumably also extragalactic and dominating the cosmic-ray flux at EeV energies, is also required to reproduce the overall spectrum and composition results for all energies down to 0.6~EeV.

Analysis of microwave sky signals, such as the cosmic microwave background, often requires component separation with multi-frequency methods, where different signals are isolated by their frequency behaviors. Many so-called "blind" methods, such as the internal linear combination (ILC), make minimal assumptions about the spatial distribution of the signal or contaminants, and only assume knowledge of the frequency dependence of the signal. The ILC is a minimum-variance linear combination of the measured frequency maps. In the case of Gaussian, statistically isotropic fields, this is the optimal linear combination, as the variance is the only statistic of interest. However, in many cases the signal we wish to isolate, or the foregrounds we wish to remove, are non-Gaussian and/or statistically anisotropic (in particular for Galactic foregrounds). In such cases, it is possible that machine learning (ML) techniques can be used to exploit the non-Gaussian features of the foregrounds and thereby improve component separation. However, many ML techniques require the use of complex, difficult-to-interpret operations on the data. We propose a hybrid method whereby we train an ML model using only combinations of the data that $\textit{do not contain the signal}$, and combine the resulting ML-predicted foreground estimate with the ILC solution to reduce the error from the ILC. We demonstrate our methods on simulations of extragalactic temperature and Galactic polarization foregrounds, and show that our ML model can exploit non-Gaussian features, such as point sources and spatially-varying spectral indices, to produce lower-variance maps than ILC - eg, reducing the variance of the B-mode residual by factors of up to 5 - while preserving the signal of interest in an unbiased manner. Moreover, we often find improved performance when applying our model to foreground models on which it was not trained.

Lyman break galaxies (LBGs) are promising probes for clustering measurements at high redshift, $z>2$, a region only covered so far by Lyman-$\alpha$ forest measurements. In this paper, we investigate the feasibility of selecting LBGs by exploiting the existence of a strong deficit of flux shortward of the Lyman limit, due to various absorption processes along the line of sight. The target selection relies on deep imaging data from the HSC and CLAUDS surveys in the $g,r,z$ and $u$ bands, respectively, with median depths reaching 27 AB in all bands. The selections were validated by several dedicated spectroscopic observation campaigns with DESI. Visual inspection of spectra has enabled us to develop an automated spectroscopic typing and redshift estimation algorithm specific to LBGs. Based on these data and tools, we assess the efficiency and purity of target selections optimised for different purposes. Selections providing a wide redshift coverage retain $57\%$ of the observed targets after spectroscopic confirmation with DESI, and provide an efficiency for LBGs of $83\pm3\%$, for a purity of the selected LBG sample of $90\pm2\%$. This would deliver a confirmed LBG density of $\sim 620$ deg$^{-2}$ in the range $2.3<z<3.5$ for a $r$-band limiting magnitude $r<24.2$. Selections optimised for high redshift efficiency retain $73\%$ of the observed targets after spectroscopic confirmation, with $89\pm4\%$ efficiency for $97\pm2\%$ purity. This would provide a confirmed LBG density of $\sim 470$ deg$^{-2}$ in the range $2.8<z<3.5$ for a $r$-band limiting magnitude $r<24.5$. A preliminary study of the LBG sample 3d-clustering properties is also presented and used to estimate the LBG linear bias. A value of $b_{LBG} = 3.3 \pm 0.2 (stat.)$ is obtained for a mean redshift of 2.9 and a limiting magnitude in $r$ of 24.2, in agreement with results reported in the literature.

We present a sample of 341 "little red dots" (LRDs) spanning the redshift range $z\sim2-11$ using data from the CEERS, PRIMER, JADES, UNCOVER and NGDEEP surveys. These sources are likely heavily-reddened AGN that trace a previously-hidden phase of dust-obscured black hole growth in the early Universe. Unlike past use of color indices to identify LRDs, we employ continuum slope fitting using shifting bandpasses to sample the same rest-frame emission blueward and redward of the Balmer break. This approach allows us to identify LRDs over a wider redshift range and is less susceptible to contamination from galaxies with strong breaks that otherwise lack a rising red continuum. The redshift distribution of our sample increases at $z<8$ and then undergoes a rapid decline at $z\sim4.5$, which may tie the emergence, and obscuration, of these sources to the inside-out growth that galaxies experience during this epoch. We find that LRDs are 2-3 dex more numerous than bright quasars at $z\sim5-7$, but their number density is only 0.6-1 dex higher than X-ray and UV selected AGN at these redshifts. Within our sample, we have identified the first X-ray detected LRDs at $z=3.1$ and $z=4.66$. An X-ray spectral analysis confirms that these AGN are moderately obscured with $\log\,(N_{\rm H}/{\rm cm}^{2}$) of $23.3^{+0.4}_{-1.3}$ and $22.72^{+0.13}_{-0.16}$. Our analysis reveals that reddened AGN emission dominates their rest-optical light, while the rest-UV originates from their host galaxies. We also present NIRSpec follow-up spectroscopy of 17 LRDs that show broad emission lines consistent with AGN activity. The confirmed AGN fraction of our sample is $71\%$ for sources with F444W$<26.5$. In addition, we find three LRDs with narrow blue-shifted Balmer absorption features in their spectra, suggesting an outflow of high-density, low ionization gas from near the central engine of these faint, red AGN.

Weak magnetic fields must have existed in the early Universe, as they were sourced by the cross product of electron density and temperature gradients through the Biermann-battery mechanism. In this paper we calculate the magnetic fields generated at cosmic dawn by a variety of small-scale primordial perturbations, carefully computing the evolution of electron density and temperature fluctuations, and consistently accounting for relative velocities between baryons and dark matter. We first compute the magnetic field resulting from standard, nearly scale-invariant primordial adiabatic perturbations, making significant improvements to previous calculations. This "standard" primordial field has a root mean square (rms) of $\sim10^{-15}$ nG at $20\lesssim z \lesssim 100$, with fluctuations on $\sim$ kpc comoving scales, and could serve as the seed of present-day magnetic fields observed in galaxies and galaxy clusters. In addition, we consider early-Universe magnetic fields as a possible probe of non-standard initial conditions of the Universe on small scales $k \sim 1-10^3$ Mpc$^{-1}$. To this end, we compute the maximally-allowed magnetic fields within current upper limits on small-scale adiabatic and isocurvature perturbations. Under the current Cosmic Microwave Background spectral-distortion constraints magnetic fields could be produced with a rms of $\sim 5\times 10^{-11}$ nG at $z = 20$. Uncorrelated small-scale isocurvature perturbations within current Big-Bang Nucleosynthesis bounds could potentially enhance the magnetic field to $\sim 10^{-14}-10^{-10}$ nG at $z = 20$, depending on the specific isocurvature mode considered. While these very weak fields remain well below current observational capabilities, our work points out that magnetic fields could potentially provide an interesting window into the poorly constrained small-scale initial conditions of the Universe.

Understanding the non-linear dynamics of coupled scalar fields often necessitates simulations on a three-dimensional mesh. These simulations can be computationally expensive if a large scale separation is involved. A common solution is adaptive mesh refinement which, however, greatly increases a simulation's complexity. In this work, we present sledgehamr, an AMReX-based code package to make the simulation of coupled scalar fields on an adaptive mesh more accessible. Compatible with both GPU and CPU clusters, sledgehamr offers a flexible and customizable framework. While the code had been primarily developed to evolve axion string networks, this framework enables various other applications, such as the study of gravitational waves sourced by the dynamics of scalar fields.

We examine the descent via membrane nucleation through a landscape of vacua where the cosmological constant is given by a combination of four-form fluxes. It has been shown that this descent can be slowed exponentially for very low curvature vacua close to Minkowski space in a wide class of models satisfying certain parametric conditions, providing a possible solution to the cosmological constant problem. We explore in detail whether or not those parametric conditions are compatible with the membrane weak gravity conjecture. Whilst it is true that there is often tension, we show that this is not always the case and present an explicit model where Minkowski space is absolutely stable and the weak gravity conjecture is satisfied. This corresponds to an extension of the Bousso-Polchinski model into a generalised DBI action for four-forms. We also clarify how the landscape should be populated in a consistent model.

We present the first systematic investigation of supervised scaling laws outside of an ImageNet-like context - on images of galaxies. We use 840k galaxy images and over 100M annotations by Galaxy Zoo volunteers, comparable in scale to Imagenet-1K. We find that adding annotated galaxy images provides a power law improvement in performance across all architectures and all tasks, while adding trainable parameters is effective only for some (typically more subjectively challenging) tasks. We then compare the downstream performance of finetuned models pretrained on either ImageNet-12k alone vs. additionally pretrained on our galaxy images. We achieve an average relative error rate reduction of 31% across 5 downstream tasks of scientific interest. Our finetuned models are more label-efficient and, unlike their ImageNet-12k-pretrained equivalents, often achieve linear transfer performance equal to that of end-to-end finetuning. We find relatively modest additional downstream benefits from scaling model size, implying that scaling alone is not sufficient to address our domain gap, and suggest that practitioners with qualitatively different images might benefit more from in-domain adaption followed by targeted downstream labelling.

A possible extension of the Standard Model able to explain the recent measurement of the anomalous magnetic moment of the muon consists in adding a gauged $U(1)_{L_{\mu}-L_{\tau}}$ symmetry. If the dark matter particle is charged under this symmetry, the kinetic mixing between the new gauge boson and the photon induces dark matter-electron interactions. We derive direct detection constraints on light dark matter charged under a $U(1)_{L_{\mu}-L_{\tau}}$ symmetry with electron recoil experiments, and explore prospects with XLZD and OSCURA to close in the parameter space able to explain simultaneously the recent measurement on the anomalous magnetic moment of the muon and the observed relic density of dark matter. We further discuss the spin-dependent scattering contribution arising in this model, which was ignored previously in the literature.

We investigate the neutrino elastic differential cross-section (NDCS) and corresponding mean free path for neutral current scattering in the dense matter of a neutron star. A wide range of observed neutron star (NS) masses is considered, including the presence of $\Lambda$, $\Xi^{-}$, and $\Xi^{0}$ hyperons in the heaviest stars. Their presence significantly decreases the total neutrino mean free path in the heavier stars.

We consider head-on collisions of two particles near the event horizon. Particle 1 is outgoing, particle 2 is ingoing. We elucidate, in which case the energy $E_{c.m.}$ in the center of mass frame can grow unbounded. If the proper time between the horizon and an arbitrary point outside it for particle 1 is finite, we deal with a white hole. If it is infinite, we deal with a black hole. Particles can be either free or experience the action of a finite force. Our results are complementary to those for the standard BSW effect when particles move in the same direction. The results rely on classification of particles developed in our previous work H.V. Ovcharenko, O.B. Zaslavskii, Phys. Rev. D 108, 064029 (2023).

The equation of state of asymmetric nuclear matter as well as the neutron and proton effective masses and their partial-wave and spin-isospin decomposition are analyzed within the Brueckner--Hartree--Fock approach. Theoretical uncertainties for all these quantities are estimated by using several phase-shift-equivalent nucleon-nucleon forces together with two types of three-nucleon forces, phenomenological and microscopic. It is shown that the choice of the three-nucleon force plays an important role above saturation density, leading to different density dependencies of the energy per particle. These results are compared to the standard form of the Skyrme energy-density functional and we find that it is not possible to reproduce the BHF predictions in the $(S,T)$ channels in symmetric and neutron matter above saturation density, already at the level of the two-body interaction, and even more including the three-body interaction.

The Boundary-to-Bound (B2B) correspondence, which connects orbital and radiative observables between bound and unbound orbits, has recently been introduced and demonstrated in the perturbative regime. We produce a large number of numerical relativity simulations of bound and unbound encounters between two nonspinning equal mass black holes in order to test this correspondence in the non-perturbative regime. We focus on testing the radiated energy and angular momentum, as well as orbital parameters such as the period and periastron advance. We find that, across a wide range of eccentricities, the B2B relationships do not hold in the non-perturbative regime, thereby placing a clear limit on the applicability of these relationships. We also approximate the separatrix between bound and unbound relativistic encounters as a function of their initial energies and angular momenta.