Papers for Wednesday, Oct 18 2023

The search for rocky planet atmospheres with JWST has focused on planets transiting M dwarfs. Such planets have favorable planet-to-star size ratios, enhancing the amplitude of atmospheric features. Since the expected signal strength of atmospheric features is similar to the single-transit performance of JWST, multiple observations are required to confirm any detection. Here, we present two transit observations of the rocky planet GJ 1132 b with JWST NIRSpec G395H, covering 2.8-5.2 $\mu$m. Previous HST WFC3 observations of GJ 1132 b were inconclusive, with evidence reported for either an atmosphere or a featureless spectrum based on analyses of the same dataset. Our JWST data exhibit substantial differences between the two visits. One transit is consistent with either a H$_2$O-dominated atmosphere containing ~1% CH$_4$ and trace N$_2$O ($\chi^{2}_{\nu}$ = 1.13) or stellar contamination from unocculted starspots ($\chi^{2}_{\nu}$ = 1.36). However, the second transit is consistent with a featureless spectrum. Neither visit is consistent with a previous report of HCN. Atmospheric variability is unlikely to explain the scale of the observed differences between the visits. Similarly, our out-of-transit stellar spectra show no evidence of changing stellar inhomogeneity between the two visits - observed 8 days apart, only 6.5% of the stellar rotation rate. We further find no evidence of differing instrumental systematic effects between visits. The most plausible explanation is an unlucky random noise draw leading to two significantly discrepant transmission spectra. Our results highlight the importance of multi-visit repeatability with JWST prior to claiming atmospheric detections for these small, enigmatic planets.

We use a suite of hydrodynamics simulations of the interstellar medium (ISM) within a galactic disk, which include radiative transfer, a non-equilibrium model of molecular hydrogen, and a realistic model for star formation and feedback, to study the structure of the ISM and H$_2$ abundance as a function of local ISM properties. We show that the star formation rate and structure of the ISM are sensitive to the metallicity of the gas with a progressively smoother density distribution with decreasing metallicity. In addition to the well-known trend of the HI-H$_2$ transition shifting to higher densities with decreasing metallicity, the maximum achieved molecular fraction in the interstellar medium drops drastically at $Z \lesssim 0.2 \, Z_\odot$ as the formation time of H$_2$ becomes much longer than a typical lifetime of dense regions of the ISM. We present accurate fitting formulae for both volumetric and projected $f_\mathrm{H_2}$ measured on different scales as a function of gas metallicity, UV radiation field, and gas density. We show that when the formulae are applied to the patches in the simulated galaxy the overall molecular gas mass is reproduced to better than a factor of $\lesssim 1.5$ across the entire range of metallicities and scales. We also show that the presented fit is considerably more accurate than any of the previous $f_\mathrm{H_2}$ models and fitting formulae in the low-metallicity regime. The fit can thus be used for modeling molecular gas in low-resolution simulations and semi-analytic models of galaxy formation in the dwarf and high-redshift regimes.

Observations and theory suggest that core-collapse supernovae can span a range of explosion energies, and when sub-energetic, the shockwave initiating the explosion can decelerate to speeds comparable to the escape speed of the progenitor. In these cases, gravity will complicate the explosion hydrodynamics and conceivably cause the shock to stall at large radii within the progenitor star. To understand these unique properties of weak explosions, we develop a perturbative approach for modeling the propagation of an initially strong shock into a time-steady, infalling medium in the gravitational field of a compact object. This method writes the shock position and the post-shock velocity, density, and pressure as series solutions in the (time-dependent) ratio of the freefall speed to the shock speed, and predicts that the shock stalls within the progenitor if the explosion energy is below a critical value. We show that our model agrees very well with hydrodynamic simulations, and accurately predicts (e.g.) the time-dependent shock position and velocity and the radius at which the shock stalls. Our results have implications for black hole formation and the newly detected class of fast X-ray transients (FXTs). In particular, we propose that a ``phantom shock breakout'' -- where the outer edge of the star falls through a stalled shock -- can yield a burst of X-rays without a subsequent optical/UV signature, similar to FXTs. This model predicts the rise time of the X-ray burst, $t_{\rm d}$, and the mean photon energy, $kT$, are anti-correlated, approximately as $T \propto t_{\rm d}^{-5/8}$.

Strong gravitational lensing can be used to find otherwise invisible dark matter subhaloes. In such an analysis, the lens galaxy mass model is a significant source of systematic uncertainty. In this paper we analyse the effect of angular complexity in the lens model. We use multipole perturbations which introduce low-order deviations from pure ellipticity in the isodensity contours, keeping the radial density profile fixed. We find that, in HST-like data, multipole perturbations consistent with those seen in galaxy isophotes are very effective at causing false positive substructure detections. We show that the effectiveness of this degeneracy depends on the deviation from a pure ellipse and the lensing configuration. We find that, when multipoles of one per cent are allowed in the lens model, the area in the observation where a subhalo could be detected drops by a factor of three. Sensitivity away from the lensed images is mostly lost. However, the mass limit of detectable objects on or close to the lensed images does not change. We do not expect the addition of multipole perturbations to lens models to have a significant effect on the ability of strong lensing to constrain the underlying dark matter model. However, given the high rate of false positive detections, angular complexity beyond the elliptical power-law should be included for such studies to be reliable. We discuss implications for previous detections and future work.

The first all-sky maps of the diffuse emission of high ionization lines observed in X-rays by SRG/eROSITA, provide an excellent probe for the study of the warm-hot phase (T~10^6 K) of the circumgalactic medium (CGM) of the Milky Way (MW). In this work we analyse the O VIII line detected in the first eROSITA All-Sky Survey data (eRASS1). We fit a sky map made in a narrow energy bin around this line, with physical emission models embedded in a 3D geometry to constrain the density distribution of the warm-hot gas around our Galaxy, with a focus on mid and high (absolute) Galactic latitudes. By masking out the eROSITA bubbles and other bright extended foreground sources, we find that an oblate geometry of the warm-hot gas (T~0.15-0.17 keV), flattened around the Galactic disk with scale height z_h~1-3 kpc, best describes the eRASS1 O VIII map, with most of the observed emission resulting to be produced within a few kpc from the Sun. The additional presence of a large scale warm-hot spherical halo, while providing a minor contribute to the X-ray emission, accounts for the high O VII absorption column densities detected with XMM-Newton, as well as most of the baryon budget of the CGM of the MW. The eROSITA data carry the largest amount of information and detail of O VIII CGM intensities to date, allowing for a significant reduction of the statistical uncertainties of the inferred physical parameters.

We present DE-VAE, a variational autoencoder (VAE) architecture to search for a compressed representation of dynamical dark energy (DE) models in observational studies of the cosmic large-scale structure. DE-VAE is trained on matter power spectra boosts generated at wavenumbers $k\in(0.01-2.5) \ h/\rm{Mpc}$ and at four redshift values $z\in(0.1,0.48,0.78,1.5)$ for the most typical dynamical DE parametrization with two extra parameters describing an evolving DE equation of state. The boosts are compressed to a lower-dimensional representation, which is concatenated with standard cold dark matter (CDM) parameters and then mapped back to reconstructed boosts; both the compression and the reconstruction components are parametrized as neural networks. Remarkably, we find that a single latent parameter is sufficient to predict 95% (99%) of DE power spectra generated over a broad range of cosmological parameters within $1\sigma$ ($2\sigma$) of a Gaussian error which includes cosmic variance, shot noise and systematic effects for a Stage IV-like survey. This single parameter shows a high mutual information with the two DE parameters, and these three variables can be linked together with an explicit equation through symbolic regression. Considering a model with two latent variables only marginally improves the accuracy of the predictions, and adding a third latent variable has no significant impact on the model's performance. We discuss how the DE-VAE architecture can be extended from a proof of concept to a general framework to be employed in the search for a common lower-dimensional parametrization of a wide range of beyond-$\Lambda$CDM models and for different cosmological datasets. Such a framework could then both inform the development of cosmological surveys by targeting optimal probes, and provide theoretical insight into the common phenomenological aspects of beyond-$\Lambda$CDM models.

Despite powerful X-ray emission, some AGNs are known to either lack optical emission lines (so-called 'optically dull' AGNs) or have lines that fall on the star-forming branch of the BPT diagram ('misclassified' AGNs). Aperture effects have been proposed to explain such atypical spectra, especially when based on SDSS (3'') fibers. We use observations from VLT-MUSE with Adaptive Optics to explore the spatially resolved optical emission line properties of 4 optically dull and 1 misclassified X-ray AGN candidates. VLT-MUSE IFU spectra allow us to investigate the extent to which the aperture size affects the emission line measurements. The optically dull AGNs become detectable in deeper VLT-MUSE spectroscopic apertures having the same size (3'') as SDSS fibers, suggesting no AGN is truly lineless. However, in no case does the line become more detectable as the aperture decreases, as would be expected if dilution by strong continuum was responsible for making the lines appear weak. We also show that the misclassified X-ray AGN retains the same position on the BPT diagram in smaller apertures (down to 0.''5), demonstrating that its misclassification is not the result of the dilution by HII regions. Thus, we conclude that continuum swamping or star formation dilution, i.e., aperture effects, are not responsible for atypical lines. Rather, the AGN lines are intrinsically weak.

Large-scale cosmological galaxy formation simulations typically prevent gas in the interstellar medium (ISM) from cooling below $\approx 10^4$ K. This has been motivated by the inability to resolve the Jeans mass in molecular gas (>>$10^5\,\mathrm{M}_{\odot}$) which would result in undesired artificial clumping. We show that the classical Jeans criteria derived for Newtonian gravity are not applicable in the simulated ISM if the spacing of resolution elements representing the dense ISM is below the gravitational force softening length and gravity is therefore softened and not Newtonian. We re-derive the Jeans criteria for softened gravity in Lagrangian codes and use them to analyse gravitational instabilities at and below the hydrodynamical resolution limit for simulations with adaptive and constant gravitational softening lengths. In addition, we define criteria for which a numerical runaway collapse of dense gas clumps can occur caused by over-smoothing of the hydrodynamical properties relative to the gravitational force resolution. This effect is illustrated using simulations of isolated disk galaxies with the smoothed particle hydrodynamics code Swift. We also demonstrate how to avoid the formation of artificial clumps in gas and stars by adjusting the gravitational and hydrodynamical force resolutions.

The effect of angular momentum on galaxy formation and evolution has been studied for several decades. Our recent two papers using IllustrisTNG-100 simulation have revealed the acquisition path of the angular momentum from large-scale environment (satellites within hundreds of kpc) through the circum-galactic medium (CGM) to the stellar discs, putting forward the co-rotation scenario across the three distance scales. In real observations, although the rotation signature for the CGM and environmental three-dimensional (3d) angular momentum are difficult to obtain, line-of-sight kinematics of group member galaxies and stellar disc kinematics of central galaxies are available utilizing existing group catalogue data and integral field unit (IFU) data. In this paper, we use (1) the group catalogue of SDSS DR7 and MaNGA IFU stellar kinematic maps and (2) the group catalogue of GAMA DR4 data and SAMI IFU stellar kinematic maps, to test if the prediction above can be seen in real data. We found the co-rotation pattern between stellar discs and satellites can be concluded with 99.7 percent confidence level ($\sim 3\sigma$) when combining the two datasets. And the random tests show that the signal can be scarcely drawn from random distribution.

Recent studies of high-redshift galaxies such as GN-z11 at z=10.6 with JWST show significant amounts of nitrogen (N) in their spectra, which is unexpected. As this phenomenology appears to extend to gravitionally-lensed galaxies at Cosmic noon such as the Sunburst Arc at z=2.37, as well as Globular Clusters, we suggest the common ingredient to be very massive stars (VMS) with zero-age main sequence (ZAMS) masses in the range 100-1000 Msun. The He II in the Sunburst Arc might also be the result of the disproportionally large contribution of VMS to the total stellar contribution. We analyse the pros and cons of the previous suggestions, including classical Wolf-Rayet (cWR) stars and supermassive stars (SMS), and we conclude that only our VMS alternative ticks all the relevant boxes. We discuss the VMS mass loss history via their peculiar vertical evolution in the HR diagram resulting from a self-regulatory effect of these wind-dominated VMS, and we estimate that the large amounts of N present in star-forming galaxies may indeed result from VMS. We conclude that VMS should be included in population synthesis and chemical evolution models, and moreover that it is critical that this is done self-consistently, as a small error in their mass-loss rates would have dramatic consequences for their stellar evolution, as well as their ionising & chemical feedback.

Observing a supernova explosion shortly after it occurs can reveal important information about the physics of stellar explosions and the nature of the progenitor stars of supernovae (SNe). When a star with a well-defined edge explodes in vacuum, the first photons to escape from its surface appear as a brief shock-breakout flare. The duration of this flare can extend to at most a few hours even for nonspherical breakouts from supergiant stars, after which the explosion ejecta should expand and cool. Alternatively, for stars exploding within a distribution of sufficiently dense optically thick circumstellar material, the first photons escape from the material beyond the stellar edge, and the duration of the initial flare can extend to several days, during which the escaping emission indicates photospheric heating. The difficulty in detecting SN explosions promptly after the event has so far limited data regarding supergiant stellar explosions mostly to serendipitous observations that, owing to the lack of ultraviolet (UV) data, were unable to determine whether the early emission is heating or cooling, and hence the nature of the early explosion event. Here, we report observations of SN 2023ixf in the nearby galaxy M101, covering the early days of the event. Using UV spectroscopy from the Hubble Space Telescope (HST) as well as a comprehensive set of additional multiwavelength observations, we trace the photometric and spectroscopic evolution of the event and are able to temporally resolve the emergence and evolution of the SN emission.

Advanced primordial chemistry networks have been developed to model the collapse of metal-free baryonic gas within the gravitational well of dark matter (DM) halos and its subsequent collapse into Population III stars. At the low densities of 10^-26-10^-21 g cm-3 (10-3-10^2 cm-3) the collapse is dependent on H2 production, which is a function of the compressional heating provided by the DM potential. Once the gas decouples from the DM, the temperature-density relationship follows a well established path dictated by various chemical reactions until the formation of the protostar at 10^-4 g cm-3 (10^19 cm-3). Here we explore the feasibility of replacing the chemical network (CN) with a barotropic equation of state (EoS) just before the formation of the first protostar, to reduce the computational load of simulating the further fragmentation, evolution and characteristics of the very high density gas. We find a significant reduction in fragmentation when using the EoS. The EoS method produces a protostellar mass distribution that peaks at higher masses when compared to CN runs. The change in fragmentation behaviour is due to a lack of cold gas falling in through the disc around the first protostar when using an EoS. Despite this, the total mass accreted across all sinks was invariant to the switch to an EoS, hence the star formation rate (Msun yr^-1) is accurately predicted using an EoS. The EoS routine is approximately 4000 times faster than the CN, however this numerical gain is offset by the lack of accuracy in modelling secondary protostar formation and hence its use must be considered carefully.

We highlight the role of weak lensing measurements from current and upcoming stage-IV imaging surveys in the search for cosmic inflation, specifically in measuring the scalar spectral index $n_s$. To do so, we combine the Dark Energy Survey 3 years of observation weak lensing and clustering data with BICEP/Keck, Planck and Sloan Digital Sky Survey data in $r\Lambda$CDM where $r$ is the tensor-to-scalar ratio. While there is no significant improvement in constraining power, we obtain a 1$\sigma$ shift on $n_s$. Additionally, we forecast a weak lensing and clustering data vector from the 10-year Legacy Survey of Space and Time by the Vera C. Rubin Observatory and show its combination with current data would improve their $n_s$ constraints by 25$\%$ in $r\Lambda$CDM.

The assembly of galaxies over cosmic time is tightly connected to the assembly of their host dark matter halos. We investigate the stellar mass growth history and the chemical enrichment history of central galaxies in SDSS-MaNGA. We find that the derived stellar metallicity of passive central galaxies is always higher than that of the star-forming ones. This stellar metallicity enhancement becomes progressively larger towards low-mass galaxies (at a given epoch) and earlier epochs (at a given stellar mass), which suggests strangulation as the primary mechanism for star formation quenching in central galaxies not only in the local universe, but also very likely at higher redshifts up to $z\sim3$. We show that at the same present-day stellar mass, passive central galaxies assembled half of their final stellar mass $\sim 2$ Gyr earlier than star-forming central galaxies, which agrees well with semi-analytic model. Exploring semi-analytic model, we find that this is because passive central galaxies reside in, on average, more massive halos with a higher halo mass increase rate across cosmic time. As a consequence, passive central galaxies are assembled faster and also quenched earlier than their star-forming counterparts. While at the same present-day halo mass, different halo assembly history also produces very different final stellar mass of the central galaxy within, and halos assembled earlier host more massive centrals with a higher quenched fraction, in particular around the "golden halo mass" at $10^{12}\mathrm{M_\odot}$. Our results call attention back to the dark matter halo as a key driver of galaxy evolution.

In the flourishing field of gravitational-wave astronomy, accurately inferring binary black hole merger formation channels is paramount. The Bayesian hierarchical model selection analysis offers a promising methodology (see, e.g., One Channel to Rule Them All, Zevin et al. 2021). However, recently, Cheng et al. (2023) highlighted a critical caveat: observed channels absent in known models can bias branching fraction estimates. In this research note, we introduce a test to detect missing channels in such analyses. Our findings show a commendable success rate in identifying these elusive channels. Yet, in scenarios where missing channels closely overlap with recognized ones, discerning the difference remains challenging.

(Abridged) Electron-molecule interaction is a fundamental process in radiation-driven chemistry in space, from the interstellar medium to comets. Therefore, knowledge of interaction cross-sections is key. While there has been a plethora of studies of total ionization cross-sections, data is often spread over many sources, or not public or readily available. We introduce the Astrochemistry Low-energy Electron Cross-Section (ALeCS) database, a public database for electron interaction cross-sections and ionization rates for molecules of astrochemical interest. In this work, we present the first data release comprising total ionization cross-sections and ionization rates for over 200 neutral molecules. We include optimized geometries and molecular orbital energies at various levels of theory, and for a subset of the molecules, the ionization potentials. We compute total ionization cross-sections using the binary-encounter Bethe model and screening-corrected additivity rule, and ionization rates and reaction network coefficients for molecular cloud environments for $>$200 neutral molecules ranging from diatomics to complex organics. We demonstrate that our binary-encounter Bethe cross-sections agree well with experimental data. We show that the ionization rates scale roughly linearly with the number of constituent atoms in the molecule. We introduce and describe the public ALeCS database. For the initial release, we include total ionization cross-sections for $>$200 neutral molecules and several cations and anions calculated with different levels of quantum chemistry theory, the chemical reaction rates for the ionization, and network files in the formats of the two most popular astrochemical networks, the KIDA and UMIST. The database will be continuously updated for more molecules and interactions.

We investigate sharp structures visible in solar magnetic field tracers. It is shown that the sunspot magnetic boundaries do not coincide with the photometric ones. Moreover, there is no clear boundary of the magnetic field in the vicinity of sunspots. Thus, the widely accepted concept of magnetic tubes with sharp edges is not always correct and should be used with caution. It is also shown that even in the moments of complete absence of visible spots on the Sun, there are magnetic fields over 800 Gauss. The nature of these strong magnetic fields remains unclear; they may originate at relatively small depths under the photosphere.

Most Active Galactic Nuclei (AGN) in the local Universe are obscured. In these obscured AGN an excess is usually observed in the soft X-rays below ~2 keV above the absorbed X-ray continuum. This spectral component is associated with the scattering of X-ray photons off free electrons in the Narrow Line Region (NLR), and/or to photoionised lines. Recent studies have found that in highly obscured AGN this component has lower flux relative to the primary X-ray continuum than in less obscured AGN. This is measured by the scattering fraction, or fscatt, which is the ratio of the scattered flux to the continuum. Here, we use the ray-tracing platform RefleX to perform simulations of scattered X-ray radiation to test two possible explanations for this phenomenon: (1) sources with lower fscatt are viewed at higher inclinations or (2) low fscatt sources are characterized by larger covering factors. We consider a conical NLR of free electrons, while allowing the column density and opening angle (and hence covering factor) to vary. We also consider electron densities inferred from observations, and from simulations carried out with the spectral synthesis code Cloudy. Our simulations show fscatt is expected to be related to both the inclination angle and covering factor of the torus; however, the observed negative correlation between fscatt and NH can only be explained by a positive relation between the column density and the covering factor of the obscuring material. Additional contributions to fscatt can come from unresolved photoionised lines and ionised outflowing gas.

Pulsar timing is a powerful tool that, by accounting for every rotation of a pulsar, precisely measures the spin frequency, spin frequency derivatives, astrometric position, binary parameters when applicable, properties of the ISM, and potentially general relativistic effects. Typically, this process demands fairly stringent scheduling requirements for monitoring observations as well as deep domain knowledge to "phase connect" the timing data. We present an algorithm that automates the pulsar timing process for binary pulsars, whose timing solutions have an additional level of complexity, although the algorithm works for isolated pulsars as well. Using the statistical F-test and the quadratic dependence of the reduced $\chi ^2$ near a minimum, the global rotation count of a pulsar can be determined efficiently and systematically. We have used our algorithm to establish timing solutions for two newly discovered binary pulsars, PSRs J1748$-$2446aq and J1748$-$2446at, in the globular cluster Terzan 5, using $\sim$70 Green Bank Telescope observations from the last 13 years.

We present a study of six dusty and gaseous pillars (containing the HH 1004 and HH 1010 objects) and globules (that contain the HH 666, HH 900, HH 1006, and HH 1066 objects) localized in the Carina nebula using sensitive and high angular resolution ($\sim$0.3$''$) Atacama Large Millimeter/Sub-millimeter Array (ALMA) observations. This is a more extensive study that the one presented in \citet{Cortes}. As in this former study, we also analyzed the 1.3 mm continuum emission and C$^{18}$O(2$-$1), N$_2$D$^+$(3$-$2) and $^{12}$CO(2$-$1) spectral lines. These new observations revealed the molecular outflows emanating from the pillars, the dusty envelopes$+$disks that are exciting them, and the extended HH objects far from their respective pillars. We reveal that the masses of the disks$+$envelopes are in a range of 0.02 to 0.38 M$_\odot$, and those for the molecular outflows are of the order of 10$^{-3}$ M$_\odot$, which suggests that their exciting sources might be low- or intermediate-mass protostars as already revealed in recent studies at infrared and submillimeter bands. In the regions associated with the objects HH 900 and HH 1004, we report multiple millimeter continuum sources, from where several molecular outflows emanate.

Although the majority of hot magnetic stars have extremely stable, $\sim$kG strength surface magnetic fields with simple topologies, a subset undergo small-scale explosions due to centrifugal breakout (CBO). The resulting small-scale flares are typically below the sensitivity of current magentospheric diagnostics and do not generate detectable transient signatures. However, a recently reported radio flare from the hot magnetic star CU Vir suggests that some of the most energetic events do reach detectable levels. Motivated by this, we searched for transient radio sources in the first two epochs of the VLITE Commensal Sky Survey (VCSS) at the position of 761 hot magnetic stars. We report three detections. A false association analysis shows a less than 1% probability that the sources are imaging artifacts. We then examine the stellar parameters of the three stars to understand if they are likely to produce flares. We conclude that while at this stage we cannot make a definitive association of the detections with the stars, the current data are consistent with the hypothesis that the flares originate in the stellar magnetospheres.

Large-amplitude accretion outbursts in young stars are expected to play a central role in proto-stellar assembly. Yet most outburst identifications and detailed studies have resulted from searches in optical time domain surveys, which are not sensitive to events located in heavily dust-obscured regions of the Galactic plane. Here, we present the discovery of WTP$\,$10aaauow, a large-amplitude mid-infrared (MIR) outburst identified in a systematic search of NEOWISE data using image subtraction techniques designed to recover transients in dense Galactic stellar fields. The new outburst is located towards the RCW$\,$49 star-forming region, and estimated to be at a distance of $\approx 4\,$kpc. Concurrent with the MIR brightening, we show that the source underwent a $\gtrsim5\,$mag outburst in the optical and near-infrared (NIR) bands around 2014-2015 reaching a peak luminosity of $\approx260\,$L$_\odot$, followed by a slow decline over the next 7 years. Analysis of the pre- and post-outburst spectral energy distributions reveal a pre-outburst stellar photosphere at a temperature of $3600-4400\,$K, surrounded by a likely two-component dust structure similar to that of a flat spectrum or Class I type YSO. We present optical and NIR follow-up spectroscopy of the source, that show a GK-type spectrum in the optical bands exhibiting complex line profiles in strong absorption features, and evidence for a wind reaching a terminal velocity of $\approx 400\,$km$\,$s$^{-1}$. The NIR bands are characterized by a cooler M-type spectrum exhibiting a forest of atomic and molecular features. All together, the spectra demonstrate that WTP\,10aaauow is a new FU Ori outburst. Ongoing systematic infrared searches will continue to reveal the extent of this population in the Galactic disk.

The cosmic infrared background (CIB) is the accumulated infrared (IR) radiation mainly from interstellar dust heated up by early stars. In this work, we measure the cross-correlation between galaxies from the unWISE catalog and the CIB maps from the Planck satellite to simultaneously constrain the cosmic star formation rate (SFR), dust spectral energy distribution (SED), and the halo occupation distribution (HOD). The unWISE galaxy catalog is divided into three tomographic bins centered at $z\sim 0.6, 1.1, 1.5$, and the CIB maps are at 353, 545, and 857 $\mathrm{GHz}$. We measure the cross-correlations between these galaxy samples and CIB maps and get a 194$\sigma$ signal within an angular scale $100<\ell<2000$, from which we constrain two CIB halo models from previous literature and one new model. The SFR, SED, and HOD model parameters are constrained consistently among the three models. Specifically, the dust temperature at $z=0$ is constrained $T_0=21.44^{+1.03}_{-1.06} \,\mathrm{K}$, which is slightly lower than $T_0=24.4\pm1.9 \,\mathrm{K}$ measured by the Planck collaboration. The halo mass that gives the most efficient star formation is around $10^{11.8}M_{\odot}$. From the model parameters, combined with the SFR density at $z=0$ synthesized from multi-wavelength observations, we break the degeneracy between SED and SFR and recover the cosmic star formation history that is consistent with multi-wavelength surveys. We also constrain the graybody SED model in agreement with previous measurements from infrared flux stacking. From the HOD constraints, we derive an increasing trend of galaxy linear bias along redshifts that agrees with the results from cross- and auto-correlation with unWISE galaxies. This study indicates the power of using CIB-galaxy cross-correlation to study star formation, dust, and abundance of galaxies across cosmic time.

Time-Division Multiplexing is the readout architecture of choice for many ground and space experiments, as it is a very mature technology with proven outstanding low-frequency noise stability, which represents a central challenge in multiplexing. Once fully populated, each of the two BICEP Array high frequency receivers, observing at 150GHz and 220/270GHz, will have 7776 TES detectors tiled on the focal plane. The constraints set by these two receivers required a redesign of the warm readout electronics. The new version of the standard Multi Channel Electronics, developed and built at the University of British Columbia, is presented here for the first time. BICEP Array operates Time Division Multiplexing readout technology to the limits of its capabilities in terms of multiplexing rate, noise and crosstalk, and applies them in rigorously demanding scientific application requiring extreme noise performance and systematic error control. Future experiments like CMB-S4 plan to use TES bolometers with Time Division/SQUID-based readout for an even larger number of detectors.

Relations between neutron star properties that do not depend on the nuclear equation of state offer insights on neutron star physics and have practical applications in data analysis. Such relations are obtained by fitting to a range of phenomenological or nuclear physics equation of state models, each of which may have varying degrees of accuracy. In this study we revisit commonly-used relations and re-assess them with a very flexible set of phenomenological nonparametric equation of state models that are based on Gaussian Processes. Our models correspond to two sets: equations of state which mimic hadronic models, and equations of state with rapidly changing behavior that resemble phase transitions. We quantify the accuracy of relations under both sets and discuss their applicability with respect to expected upcoming statistical uncertainties of astrophysical observations. We further propose a goodness-of-fit metric which provides an estimate for the systematic error introduced by using the relation to model a certain equation-of-state set. Overall, the nonparametric distribution is more poorly fit with existing relations, with the I--Love--Q relations retaining the highest degree of universality. Fits degrade for relations involving the tidal deformability, such as the Binary-Love and compactness-Love relations, and when introducing phase transition phenomenology. For most relations, systematic errors are comparable to current statistical uncertainties under the nonparametric equation of state distributions.

High-fidelity simulated astronomical images are an important tool in developing and measuring the performance of image-processing algorithms, particularly for high precision measurements of cosmic shear -- correlated distortions of images of distant galaxies due to weak gravitational lensing caused by the large-scale mass distribution in the Universe. For unbiased measurements of cosmic shear, all other sources of correlated image distortions must be modeled or removed. One such source is the correlated blurring of images due to optical turbulence in the atmosphere, which dominates the point-spread function (PSF) for ground-based instruments. In this work, we leverage data from weather forecasting models to produce wind speeds and directions, and turbulence parameters, that are realistically correlated with altitude. To study the resulting correlations in the size and shape of the PSF, we generate simulated images of the PSF across a ~10 square-degree field of view -- the size of the camera focal plane for the Vera C. Rubin Observatory in Chile -- using weather data and historical seeing for a geographic location near the Observatory. We make quantitative predictions for two-point correlation functions (2PCF) that are used in analyses of cosmic shear. We observe a strong anisotropy in the two-dimensional 2PCF, which is expected based on observations in real images, and study the dependence of the orientation of the anisotropy on dominant wind directions near the ground and at higher altitudes. The code repository for producing the correlated weather parameters for input to simulations (psf-weather-station) is public at https://github.com/LSSTDESC/psf-weather-station.

GMagAO-X is a visible to NIR extreme adaptive optics (ExAO) system that will be used at first light for the Giant Magellan Telescope (GMT). GMagAO-X is designed to deliver diffraction-limited performance at visible and NIR wavelengths (6 to 10 mas) and contrasts on the order of $10^{-7}$. The primary science case of GMagAO-X will be the characterization of mature, and potentially habitable, exoplanets in reflected light. GMagAO-X employs a woofer-tweeter system and includes segment phasing control. The tweeter is a 21,000 actuator segmented deformable mirror (DM), composed of seven individual 3,000 actuator DMs. This new ExAO framework of seven DMs working in parallel to produce a 21,000 actuator DM significantly surpasses any current or near future actuator count for a monolithic DM architecture. Bootstrapping, phasing, and high order sensing are enabled by a multi-stage wavefront sensing system. GMT's unprecedented 25.4 m aperture composed of seven segments brings a new challenge of co-phasing massive mirrors to 1/100th of a wavelength. The primary mirror segments of the GMT are separated by large >30 cm gaps so there will be fluctuations in optical path length (piston) across the pupil due to vibration of the segments, atmospheric conditions, etc. We have developed the High Contrast Adaptive-optics Testbed (HCAT) to test new wavefront sensing and control approaches for GMT and GMagAO-X, such as the holographic dispersed fringe sensor (HDFS), and the new ExAO parallel DM concept for correcting aberrations across a segmented pupil. The CoDR for GMagAO-X was held in September 2021 and a preliminary design review is planned for early 2024. In this paper we will discuss the science cases and requirements for the overall architecture of GMagAO-X, as well as the current efforts to prototype the novel hardware components and new wavefront sensing and control concepts for GMagAO-X on HCAT.

In the last two decades many people have been searching for the optimal wavefront sensor as it can boost the performance of high-contrast imagining by orders of magnitude on the ELTs. According classical information theory, the optimal sensitivity of a wavefront sensor is 1/2 radian rms per photon. We show that classical limit is also the quantum metrology limit for starlight, which means that 1/2 radian rms per photon is really the limit. This proceeding introduces the Phase Induced Amplitude Apodized Zernike Wavefront sensor. The PIAA-ZWFS modifies a standard ZWFS with a set of aspheric lenses to increase its sensitivity. The optimized system reaches the fundamental limit for all spatial frequencies >1.7 cycles/pupil and is very close to the limit for the spatial frequencies <1.7 cycles/pupil. The PIAA-ZWFS can be seamlessly integrated with the PIAA-CMC coronagraphy. This makes the PIAA-ZWFS an ideal candidate as wavefront sensor for high-contrast imaging.

Uncorrected wavefront errors create speckle noise in high-contrast observations at small inner-working angles. These speckles can be sensed and controlled by using coronagraph integrated wavefront sensors. Here, we will present how the Phase Induced Amplitude Apodized Complex Mask Corongraph (PIAACMC) can be integrated with both a Self-Coherent Camera (SCC) for focal plane wavefront sensing and an extremely sensitivity high-order pupil plane Zernike wavefront sensor (ZWFS). Non-common path aberrations can be completely erased by integrating both sensors into the PIAACMC, which is of extremely high importance in high-contrast imaging.

In data from the New Horizons encounter with Pluto in 2015, attention was called to a crater named Kiladze and its surroundings because of the water ice spectral properties, which contrast with the primarily methane ice regional surface composition. The water ice carries the spectral signature of an ammoniated compound, similar to that seen at two other sites on Pluto where cryovolcanism has been identified. The faulted structure of Kiladze, including shaping by numerous collapse pits and the distorted shape of the crater, are compatible with the surroundings in Hayabusa Terra, east of Spunik Planitia. They are further compatible with an interpretation as a resurgent caldera formed during a past era of active cryovolcanic period that appears to be significantly more recent than the overall age of the planet's surface, possibly in the last several million years. In view of the size of the caldera and the large scale of the surrounding distribution of water ice, we propose that Kiladze is a "supervolcano" in which one or more explosive events has scattered more than ~1000 km$^{3}$ of icy cryomagma erupted from the interior onto the surface.

RelativisticDynamics.jl is an open-source Julia package for relativistic spin-orbital dynamics in the gravitational strong-field of a Kerr spacetime. Existing codes for modelling the dynamics of spinning objects like pulsars in the strong-field regime are generally lacking, since such systems occupy an intermediate regime that is generally overlooked. At the "low" end of this regime there are post-Newtonian descriptions which neglect the influence of the pulsar spin on the underlying spacetime metric ("spin-curvature" coupling). At the "high" end there are the full numerical relativity solutions which are primarily applicable toe two black holes with a mass ratio $\mathcal{O}(1)$, and are computationally intractable for pulsar systems observed over a large number of orbital cycles. RelativisticDynamics.jl aims to bridge this gap by providing a modern, fast code for accurate numerical evolution of spinning relativistic systems via the Mathisson-Papetrou-Dixon formalism. Julia is a modern language that solves the "two language problem", enabling fast dynamic typing and JIT compilation on conjunction with petaflop performance, comparable with numerical languages that are better known in the astrophysics community such as C or Fortran. RelativisticDynamics.jl is written to be fully type flexible, being able to support arbitrary number formats, and fully differentiable via automatic differentiation.

While the most exciting explanation of the observed dust asymmetries in protoplanetary disks is the presence of protoplanets, other mechanisms can also form the dust features. This paper presents dual-wavelength Atacama Large Millimeter/submillimeter Array (ALMA) observations of a large asymmetric dusty ring around the M-type star CIDA 9A. We detect a dust asymmetry in both 1.3 mm and 3.1 mm data. To characterize the asymmetric structure, a parametric model is used to fit the observed visibilities. We report a tentative azimuthal shift of the dust emission peaks between the observations at the two wavelengths. This shift is consistent with a dust trap caused by a vortex, which may be formed by an embedded protoplanet or other hydrodynamical instabilities, such as a dead zone. Deep high-spatial observations of dust and molecular gas are needed to constrain the mechanisms that formed the observed millimeter cavity and dust asymmetry in the protoplanetary disk around CIDA 9A.

Observing the first generation of stars, Population III (Pop III), is still a challenge even with the James Webb Space Telescope (JWST) due to their faintness. Instead, searching for fossil records of Pop III stars in nearby dwarf galaxies provides an alternative method for studying their physical properties. It is intriguing that a star recently discovered in the Sculptor dwarf galaxy, named AS0039, is considered to show the unique signature of a Pop~III star. The detailed abundance patterns of AS0039 are well-matched with those predicted by nucleosynthesis models for Pop~III exploding as an energetic hypernova (HN), confirming its potential to provide insight into the properties of the first stars. This study aims to explore the environmental conditions required for the formation of such a unique star using cosmological hydrodynamic zoom-in simulations on dwarf galaxies with a mass of M_vir~10^8 solar mass at z=0 while varying the fraction of Pop~III stars that undergo HNe. Our simulations identify rapid gas inflow (~0.08 solar mass/yr) as a possible factor in facilitating the formation of stars similar to AS0039. Alternatively, the delayed formation of subsequent Pop~II stars in the gas-enriched environment may lead to low-metallicity stars like AS0039. Additionally, using the A-SLOTH code, we investigate the probability of finding remnants of Pop II stars with HN signatures in nearby dwarf satellite galaxies. We suggest that the most likely dwarf galaxies to contain HN signatures are massive satellites with a probability of 40% in the range of M_peak~10^{10}-10^{11} solar mass and M_star~10^7-10^8 solar mass, considering observational limitations.

The dichotomy referred to as a partition or separation of a whole into two parts and specifically, the dichotomy is very important feature of Mars between the Southern and Northern regions of Mars, and another thing that makes Mars very special that is the occurrence of Dust Devils on Mars. So, we studied and survey the dust devils occurrence on Mars in different Martian Years on the whole Mars. We create a 2D map of Martian Surface and plot the coordinates where the dust devils are captured during their activity and use those locations where they leave the tracks behind them after passing from those locations and those tracks commonly referred to as a Dust Devils Tracks. So, we plot them in two different categories Direct Observations and Indirect Observations of Dust Devils and in the map, we have plotted the locations (coordinates) of DDs shows a variation in locations of occurrences with the Dichotomy the serpent like variation we observed and we find most of the dust devil are occurred on the Dichotomy and the nearby regions of it which follows the serpent like trajectory of dichotomy of Mars and another observation shows that these locations lie on the remanent magnetic fields zones of mars which referred to as crustal magnetic fields of Mars this previously unknown relationship between crustal magnetic fields, dichotomy of mars and occurrence of dust devils is being examined here.

We investigate the explicit role of negative local non-Gaussianity, $f_{\rm NL}$, in suppressing the abundance of primordial black holes (PBHs) in the single-field model of Galileon inflation. PBH formation requires significantly enhancing the scalar power spectrum, which greatly affects their abundance. The associated frequencies in the nHz regime are also sensitive to the generation of scalar-induced gravitational waves (SIGWs) which may explain the current data from the pulsar timing arrays (PTAs). Our analysis using the threshold statistics on the compaction function demonstrates that Galileon theory not only avoids PBH overproduction using the curvature perturbation enhancements that give $f_{\rm NL} \sim {\cal O}(-6)$, but also generates SIGWs that conform well with the PTA data.

We present the first results of the eXtreme UV Environments (XUE) James Webb Space Telescope (JWST) program, that focuses on the characterization of planet forming disks in massive star forming regions. These regions are likely representative of the environment in which most planetary systems formed. Understanding the impact of environment on planet formation is critical in order to gain insights into the diversity of the observed exoplanet populations. XUE targets 15 disks in three areas of NGC 6357, which hosts numerous massive OB stars, among which some of the most massive stars in our Galaxy. Thanks to JWST we can, for the first time, study the effect of external irradiation on the inner ($< 10$ au), terrestrial-planet forming regions of proto-planetary disks. In this study, we report on the detection of abundant water, CO, CO$_2$, HCN and C$_2$H$_2$ in the inner few au of XUE 1, a highly irradiated disk in NGC 6357. In addition, small, partially crystalline silicate dust is present at the disk surface. The derived column densities, the oxygen-dominated gas-phase chemistry, and the presence of silicate dust are surprisingly similar to those found in inner disks located in nearby, relatively isolated low-mass star-forming regions. Our findings imply that the inner regions of highly irradiated disks can retain similar physical and chemical conditions as disks in low-mass star-forming regions, thus broadening the range of environments with similar conditions for inner disk rocky planet formation to the most extreme star-forming regions in our Galaxy.

We present the results of the first X-ray polarimetric observation of the low-mass X-ray binary 4U 1957+115, performed with the Imaging X-ray Polarimetry Explorer IXPE in May 2023. The binary system resides in a high/soft spectral state since discovery and is thought to host a black hole. The $\sim$571 ks IXPE observation reveals a linear polarization degree of $1.9\% \pm 0.4$\% and a polarization angle of $-41.\!\!^\circ8 \pm 5.\!\!^\circ7$ in the 2-8 keV energy range. Spectral modelling is consistent with the dominant contribution from the standard accretion disc, while polarimetric data suggest a significant role of returning radiation, i.e. photons that are bent by strong gravity effects and forced to return to the disc surface, where they can be reflected before eventually reaching the observer. In this setting, we find that models with a black hole spin lower than 0.96 and inclination lower than 50$\hbox{$^\circ$}$ are disfavoured.

In warm (equilibrium temperature <1000 K) gas giant exoplanet atmospheres, the observation of trace species in abundances deviating from thermochemical equilibrium predictions could be used as an indicator of disequilibrium chemical processes, such as photochemistry. To predict which compounds could be used as such tracers, it is therefore essential to study how photochemical processes affect their abundances. For this purpose, we investigated experimentally the efficiency of the photochemical formation of hydrocarbons in gas mixtures representative of warm gas giant atmospheres as a function of the gas temperature at millibar pressures. We find that, compared to thermal reactions alone, photochemistry efficiently promotes, under the studied conditions, the formation of hydrocarbons, with the detection of acetylene, ethane, and propane, as well as carbon monoxide. Therefore, our results confirm the importance of photochemistry in exoplanet atmospheres as a disequilibrium process. Ethane is the major hydrocarbon formed in our experiments, in apparent contradiction with the prediction by thermophotochemical models that acetylene should be the main hydrocarbon product. We also observe an evolution of the hydrocarbon production efficiency as a function of the temperature, a behavior not reproduced by a 0D thermophotochemical model. Additional studies are necessary to definitively understand the origin of the differences between the experimental and modeling results and to infer the importance of our results for understanding hydrocarbon formation in warm gas giant exoplanet atmospheres. Finally, our work demonstrates the importance of experimental studies together with modeling studies to accurately interpret, understand, and predict observations of exoplanet atmospheres.

The Didymos binary asteroid was the target of the Double Asteroid Redirection Test (DART) mission, which intentionally impacted Dimorphos, the smaller member of the binary system. We used the Near-Infrared Spectrograph (NIRSpec) and Mid-Infrared Instrument (MIRI) instruments on JWST to measure the 0.6-5 $\mu$m and 5-20 $\mu$m spectra of Didymos approximately two months after the DART impact. These observations confirm that Didymos belongs to the S asteroid class and is most consistent with LL chondrite composition as was previously determined from its 0.6-2.5-$\mu$m reflectance spectrum. Measurements at wavelengths $>$ 2.5 $\mu$m show Didymos to have thermal properties typical for an S-complex asteroid of its size and to be lacking absorptions deeper than $\sim$2\% due to OH or H2O. Didymos' mid-infrared emissivity spectrum is within the range of what has been observed on S-complex asteroids observed with Spitzer Space Telescope and is most consistent with emission from small ($<$ 25 $\mu$m) surface particles. We conclude that the observed reflectance and physical properties make the Didymos system a good proxy for the type of ordinary chondrite asteroids that cross near-Earth space, and a good representative of likely future impactors.

Although anomalous surface abundances are often observed in A-type main-sequence stars (known as chemically peculiar stars; e.g., metallic line stars or Am stars), our understanding about the behavior of aluminium is still insufficient. Actually, even whether Al is overabundant or underabundant in Am stars is not clarified. This is presumably because most of the previous studies employed the Al I 3944/3961 lines with the assumption of LTE, despite that a considerable non-LTE effect is expected in this resonance doublet. With an aim to shed light on this issue, extensive statistical-equilibrium calculations on Al I/Al II were carried out for a wide range of atmospheric parameters, based on which the non-LTE Al abundances were determined by applying the spectrum-fitting technique to the Al I 3944/3961 lines for 63 A-type dwarfs (7000 < Teff < 10000 K) of comparatively lower rotational velocities (vsini < 100 km/s). The following results were obtained. (1) The non-LTE corrections (Delta) are positive (reflecting the importance of overionization) and significantly large (0.3 < Delta < 1.0 dex depending on Teff; generally Delta_3944 < Delta_3961). (2) By applying these corrections (and indispensable inclusion of Balmer line wings as background opacity), consistent non-LTE abundances for both lines could be obtained, and the serious zero-point discrepancy (considerably negative [Al/H] for normal metallicity stars of [Fe/H]~ 0) found in old studies has been settled. (3) Al abundances of A-type stars are almost in proportion to [Fe/H] (tending to be overabundant in Am stars) with an approximate relation of [Al/H]~1.2[Fe/H]. which is qualitatively consistent with the prediction of the diffusion theory (suggesting an Al excess in the photosphere of Am stars).

Radial velocity surveys suggest that the Solar System may be unusual and that Jupiter-like planets have a frequency <20% around solar-type stars. However, they may be much more common in one of the closest associations in the solar neighbourhood. Young moving stellar groups are the best targets for direct imaging of exoplanets and four massive Jupiter-like planets have been already discovered in the nearby young beta Pic Moving Group (BPMG) via high-contrast imaging, and four others were suggested via high precision astrometry by the European Space Agency's Gaia satellite. Here we analyze 30 stars in BPMG and show that 20 of them might potentially host a Jupiter-like planet as their orbits would be stable. Considering incompleteness in observations, our results suggest that Jupiter-like planets may be more common than previously found. The next Gaia data release will likely confirm our prediction.

We report the results of the stellar occultation by (UII) Umbriel on September 21st, 2020. The shadow crossed the USA and Canada, and 19 positive chords were obtained. A limb parameter accounted for putative topographic features in the limb fittings. Ellipse fittings were not robust - only upper limits were derived for the true size/shape of a putative Umbriel ellipsoid. The adopted spherical solution gives radius = 582.4 +/- 0.8 km, smaller/close to 584.7 +/- 2.8 km from Voyager II. The apparent ellipse fit results in a true semi-major axis of 584.9 +/- 3.8 km, semi-minor axes of 582.3 +/- 0.6 km and true oblateness of 0.004 +/- 0.008 for a putative ellipsoid. The geometric albedo was pV = 0.26 +/- 0.01. The density was rho = 1.54 +/- 0.04 g cm-3. The surface gravity was 0.251 +/- 0.006 m s-2 and the escape velocity 0.541 +/- 0.006 km s-1 . Upper limits of 13 and 72 nbar (at 1 sigma and 3 sigma levels, respectively) were obtained for the surface pressure of a putative isothermal CO2 atmosphere at T = 70 K. A milliarcsecond precision position was derived: RA = 02h 30m 28.84556s +/- 0.1 mas, DE = 14o 19' 36.5836" +/- 0.2 mas. A large limb parameter of 4.2 km was obtained, in striking agreement with opposite southern hemisphere measurements by Voyager II in 1986. Occultation and Voyager results indicate that the same strong topography variation in the surface of Umbriel is present on both hemispheres.

Cosmic shear statistics, such as the two-point correlation function (2PCF), can be evaluated with the PDF-SYM method instead of the traditional weighted-sum approach. It makes use of the full PDF information of the shear estimators, and does not require weightings on the shear estimators, which can in principle introduce additional systematic biases. This work presents our constraints on $S_8$ and $\Omega_m$ from the shear-shear correlations using the PDF-SYM method. The data we use is from the z-band images of the Dark Energy Camera Legacy Survey (DECaLS), which covers about 10000 deg$^2$ with more than 100 million galaxies. The shear catalog is produced by the FQ method, and well tested on the real data itself with the field-distortion effect. Our main approach is called quasi-2D as we do use the photo-$z$ information of each individual galaxy, but without dividing the galaxies into redshift bins. We mainly use galaxy pairs within the redshift interval between 0.2 and 1.3, and the angular range from $4.7$ to $180$ arcmin. Our analysis yields $S_8=0.762 \pm 0.026$ and $\Omega_m=0.234 \pm 0.075$, with the baryon effects and the intrinsic alignments included. The results are robust against redshift uncertainties. We check the consistency of our results by deriving the cosmological constaints from auto-correlations of $\gamma_1$ and $\gamma_2$ separately, and find that they are consistent with each other, but the constraints from the $\gamma_1$ component is much weaker than that from $\gamma_2$. It implies a much worse data quality of $\gamma_1$, which is likely due to additional shear uncertainties caused by CCD electronics (according to the survey strategy of DECaLS). We also perform a pure 2D analysis, which gives $S_8=0.81^{+0.03}_{-0.04}$ and $\Omega_m=0.25^{+0.06}_{-0.05}$. Our findings demonstrate the potential of the PDF-SYM method for precision cosmology.

We aim to reveal the sizes of the continuum and broad emission line (BEL) emitting regions in the gravitationally lensed quasar SDSS J1004+4112 by analyzing the unique signatures of microlensing in this system. Through a comprehensive analysis of 20 spectroscopic observations acquired between 2003 and 2018, we studied the striking deformations of various BEL profiles and determined the sizes of their respective emitting regions. Our approach involves a detailed analysis of the magnitude differences in the BEL wings and their adjacent continua, and the implementation of a statistical model to quantify the distribution and impact of microlensing magnifications. To ensure a reliable baseline for no microlensing, we used the emission line cores as a reference. We then applied a Bayesian estimate to derive the size lower limits of the Ly$\alpha$, Si IV, C IV, C III], and Mg II emitting regions, as well as the sizes of the underlying continuum-emitting sources. We analyzed the outstanding microlensing-induced distortions in the line profiles of various BELs in the quasar image A, characterized by a prominent magnification of the blue part and a strong demagnification of the red part. From the statistics of microlensing magnifications and using Bayesian methods, we estimate the lower limit to the overall size of the regions emitting the BELs to be a few lt-days across, which is significantly smaller than in typically lensed quasars. The asymmetric deformations in the BELs indicate that the broad-line region is generally not spherically symmetric, and is likely confined to a plane and following the motions of the accretion disk. Additionally, the inferred continuum-emitting region sizes are larger than predictions based on standard thin-disk theory by a factor of $\sim$3.6 on average. The size-wavelength relation is consistent with that of a geometrically thin and optically thick accretion disk.

Evaporation of chromospheric plasma by particle beams has been a standard feature of models of solar flares for many decades, supported both by observations of strong hard X-ray bremsstrahlung signals, and detailed 1D hydrodynamic radiative transfer models with near-relativistic electron beams included. However in multi-dimensional models, evaporation, if included, has only been driven by heat conduction and by the impact and reflection of fast plasma outflows on the lower atmosphere. Here we present the first multi-dimensional flare simulation featuring evaporation driven by energetic electrons. We use a recent magnetohydrodynamic model that includes beam physics, but decrease the initial anomalous resistivity to create a gentler precursor phase, and improve on the dynamic resistivity treatment that determines where beams are injected. Beam-driven evaporation is achieved. The relevant factors are thermal conduction and electron beams, with the beam electrons more than doubling the kinetic energy flux, and adding 50% to the upward mass from the chromosphere. These findings finally pave the way for integrating detailed 1D flare modelling within a self-consistent 2D and 3D context. The beam fluxes from these self-consistent models can be used to directly compare multi-dimensional results with those from the externally injected beam fluxes of 1D models, as well as understand further evaporation-driven phenomena relating to beams of particles.

The recent detection of a stochastic signal in the NANOGrav 15-year data set has aroused great interest in uncovering its origin. However, the evidence for the Hellings-Downs correlations, a key signature of the gravitational-wave background (GWB) predicted by general relativity, remains inconclusive. In this paper, we search for an isotropic non-tensorial GWB, allowed by general metric theories of gravity, in the NANOGrav 15-year data set. Our analysis reveals a Bayes factor of approximately 2.5, comparing the quadrupolar (tensor transverse, TT) correlations to the scalar transverse (ST) correlations, suggesting that the ST correlations provide a comparable explanation for the observed stochastic signal in the NANOGrav data. We obtain the median and the $90\%$ equal-tail amplitudes as $\mathcal{A}_\mathrm{ST} = 7.8^{+5.1}_{-3.5} \times 10^{-15}$ at the frequency of 1/year. Furthermore, we find that the vector longitudinal (VL) and scalar longitudinal (SL) correlations are weakly and strongly disfavoured by data, respectively, yielding upper limits on the amplitudes: $\mathcal{A}_\mathrm{VL}^{95\%} \lesssim 1.7 \times 10^{-15}$ and $\mathcal{A}_\mathrm{SL}^{95\%} \lesssim 7.4 \times 10^{-17}$. Lastly, we fit the NANOGrav data with a general transverse (GT) correlations parameterized by a free parameter $\alpha$. Our analysis yields $\alpha=1.74^{+1.18}_{-1.41}$, thus excluding both the TT ($\alpha=3$) and ST ($\alpha=0$) models at the $90\%$ confidence level.

We investigate the efficacy of using the cosmic web nodes identified by the DisPerSE topological filament finder to systematically identify galaxy groups in the infall regions around massive clusters. The large random motions and infall velocities of galaxies in the regions around clusters complicate the detection and characterisation of substructures through normal group-finding algorithms. Yet understanding the co-location of galaxies within filaments and/or groups is a key part of understanding the role of environment on galaxy evolution, particularly in light of next-generation wide-field spectroscopic surveys. Here we use simulated massive clusters from TheThreeHundred collaboration and compare the derived group catalogues, (haloes with $\sigma_{v} > 300 h^{-1}$ km/s) with the critical points from DisPerSE, ran on haloes with more than 100 particles. We find that in 3D, 56\% of DisPerSE nodes are correctly identified as groups (purity) while 68\% of groups are identified as nodes (completeness). The fraction of matches increases with group mass and with distance from the host cluster centre. This rises to a completeness of 100\% for the most massive galaxy groups ($M>10^{14}$ M$_{\odot}$) in 3D, or 63\% when considering the projected 2D galaxy distribution. When a perfect match occurs between a cosmic web node and a galaxy group, the DisPerSE node density ($\delta$) serves as an estimate of the group's mass, albeit with significant scatter. We conclude that the use of a cosmic filament finder shows promise as a useful and straightforward observational tool for disentangling substructure within the infall regions of massive clusters.

We present 0$.\!\!^{\prime\prime}22$-resolution Atacama Large Millimeter/submillimeter Array (ALMA) observations of CO(2$-$1) emission from the circumnuclear gas disk in the red nugget relic galaxy PGC 11179. The disk shows regular rotation, with projected velocities near the center of 400 km s$^{-1}$. We assume the CO emission originates from a dynamically cold, thin disk and fit gas-dynamical models directly to the ALMA data. In addition, we explore systematic uncertainties by testing the impacts of various model assumptions on our results. The supermassive black hole (BH) mass ($M_\mathrm{BH}$) is measured to be $M_\mathrm{BH} = (1.91\pm0.04$ [$1\sigma$ statistical] $^{+0.11}_{-0.51}$ [systematic])$\times 10^9$ $M_\odot$, and the $H$-band stellar mass-to-light ratio $M/L_H=1.620\pm0.004$ [$1\sigma$ statistical] $^{+0.211}_{-0.107}$ [systematic] $M_\odot/L_\odot$. This $M_\mathrm{BH}$ is consistent with the BH mass$-$stellar velocity dispersion relation but over-massive compared to the BH mass$-$bulge luminosity relation by a factor of 3.7. PGC 11179 is part of a sample of local compact early-type galaxies that are plausible relics of $z\sim2$ red nuggets, and its behavior relative to the scaling relations echoes that of three relic galaxy BHs previously measured with stellar dynamics. These over-massive BHs could suggest BHs gain most of their mass before their host galaxies do. However, our results could also be explained by greater intrinsic scatter at the high-mass end of the scaling relations, or by systematic differences in gas- and stellar-dynamical methods. Additional $M_\mathrm{BH}$ measurements in the sample, including independent cross-checks between molecular gas- and stellar-dynamical methods, will advance our understanding of the co-evolution of BHs and their host galaxies.

Analysis of f-mode frequencies has provided a measure of the radius of the Sun which is lower, by a few hundredths per cent, than the photospheric radius determined by direct optical measurement. Part of this difference can be understood by recognizing that it is primarily the variation of density well beneath the photosphere of the star that determines the structure of these essentially adiabatic oscillation modes, not some aspect of radiative intensity. In this paper we attempt to shed further light on the matter, by considering a differently defined, and dynamically more robust, seismic radius, namely one determined from p-mode frequencies. This radius is calibrated by the distance from the centre of the Sun to the position in the subphotospheric layers where the first derivative of the density scale height changes essentially discontinuously. We find that that radius is more- or-less consistent with what is suggested by the f modes. In addition, the interpretation of the radius inferred from p modes leads us to understand more deeply the role of the total mass constraint in the structure inversions. This enables us to reinterpret the sound-speed inversion, suggesting that the positions of the photosphere and the adiabatically stratified layers in the convective envelope differ nonhomologously from those of the standard solar model.

Within a simulated Milky Way-like galaxy, we identify and analyse analogues of the Gaia-Enceladus (GE), Kraken and Sequoia mergers that each matches remarkably well observational results, including in velocity and chemical abundance space, and their distributions in the $j_{z}$-Energy plane. The Kraken analogue is the earliest merger and has the highest total mass ratio. Consistent with previous studies, it is chemically indistinguishable from old in-situ stars at the time of its accretion. The GE and Sequoia analogue events accrete at similar times in our simulation, both along filaments but from opposite sides of the main galaxy. The mean stellar ages of the GE and Sequoia analogues are both similar and, from our simulation results, we see that they can be separate entities and still naturally reproduce the observed properties of their stellar remnants at the present day, including the significant retrograde velocities of the Sequoia analogue remnant stars and the difference in the tracks of the two galaxies through chemical abundance space. Our results provide supporting information about the properties of these three merger events, and show for the first time that they can all be reproduced with a fully cosmological simulation, providing a possible self consistent evolutionary pathway for the Milky Way's formation.

Context. Large spectroscopic surveys devoted to the study of the Milky Way, including Gaia, use automated pipelines to massively determine the atmospheric parameters of millions of stars. The Gaia FGK Benchmark Stars are reference stars with Teff and log g derived through fundamental relations, independently of spectroscopy, to be used as anchors for the parameter scale. The first and second versions of the sample have been extensively used for that purpose, and more generally to help constrain stellar models. Aims. We provide the third version of the Gaia FGK Benchmark Stars, an extended set intended to improve the calibration of spectroscopic surveys, and their interconnection. Methods. We have compiled about 200 candidates which have precise measurements of angular diameters and parallaxes. We determined their bolometric fluxes by fitting their spectral energy distribution. Masses were determined using two sets of stellar evolution models. In a companion paper we describe the determination of metallicities and detailed abundances. Results. We provide a new set of 192 Gaia FGK Benchmark Stars with their fundamental Teff and logg, and with uncertainties lower than 2% for most stars. Compared to the previous versions, the homogeneity and accuracy of the fundamental parameters are significantly improved thanks to the high quality of the Gaia data reflecting on distances and bolometric fluxes.

In this work we analyse 3 average-luminosity hard X-ray selected AGN: ESO 506-G27, IGR J19039+3344 and NGC 7465. They have simultaneous Swift/XRT and NuSTAR data never published before and have been poorly studied at X-ray energies. These sources make for interesting targets both from a methodological and scientific point of view. Scientifically, they are of interest since they are possibly heavily absorbed objects, belong to a peculiar class and are variable both in flux and in spectral shape. Methodologically, because it is an interesting exercise to understand how existing spectral models can be applied to faint sources and how the use of NuSTAR data alone and then simultaneous and/or average data impacts on the spectral parameters determination. In this work we demonstrate that simultaneous data are not sufficient if their statistical quality is poor. Moreover, we show that also the use of time-averaged data when dealing with faint AGN does not always provide confident results as for brighter AGN. Regardless of the poor data quality employed in our analysis, we are able to provide insights into the spectral characteristics of each source. We analyse in detail for the first time the iron line complex of ESO 506-G27, finding not only the presence of the iron K$\alpha$ and K$\beta$ lines, but also of the iron K edge around 7 keV in the NuSTAR data. We also highlight changes in the absorption properties of IGR J19039+3344 and confirm NGC 7465 to be an unabsorbed type 1 LINER.

We want to investigate the influence of the powerful outflow driven by the hot core Sgr B2(N1) on the gas molecular inventory of the surrounding medium. We used the data taken as part of the 3 mm imaging spectral-line survey ReMoCA (Re-exploring Molecular Complexity with ALMA). Integrated intensity maps of SO and SiO emission reveal a bipolar structure with blue-shifted emission dominantly extending to the SE from the centre of the hot core and red-shifted emission to the NW. This is also prominently observed in emission of other S-bearing molecules and species that only contain N as a heavy element, including COMs, but also CH3OH, CH3CHO, HNCO, and NH2CHO. For a selection of COMs and simpler species, spectra were modelled under the assumption of LTE and population diagrams were derived at two positions, one in each outflow lobe. From this analysis, we obtained rotational temperatures, which are in a range of ~100-200K, and column densities. Abundances were subsequently compared to predictions of astrochemical models and to observations of L1157-B1, a position located in the well-studied outflow of the low-mass protostar L1157, and the source G+0.693-0.027, located in the Sgr B2 molecular cloud complex. Given the short distance of the analysed outflow positions to the centre of Sgr B2(N1), we propose a scenario in which a phase of hot-core chemistry (i.e. thermal desorption of ice species and high-temperature gas-phase chemistry) preceded a shock wave. The subsequent compression and further heating of the material resulted in the accelerated destruction of (mainly O-bearing) molecules. Gas-phase formation of cyanides seems to be able to compete with their destruction in the post-shock gas. Abundances of HCnN (n=3,5) are enhanced in the outflow component pointing to (additional) gas-phase formation. To confirm such a scenario, appropriate chemical shock models need to be run.

A magnetised flux rope model, "Flux Rope in 3D" (FRi3D) is used in the framework of European Heliospheric Forecasting Information Asset (EUHFORIA) for studying the evolution and propagation of coronal mass ejections (CME). In this paper, we rectify the mistake in the mentioned magnetic field profile of the FRi3D model used in Maharana et al., 2022, and we clarify the actual profile used in that work. In addition, we provide the recent updates introduced to the FRi3D implementation in EUHFORIA like optimising the "ai.fri3d" python package to reduce computational time and exploring different CME leg disconnection methods to make the numerical implementation more stable.

Void galaxies are essential to understand the physical processes that drive galaxy evolution as they are less affected by external factors than galaxies in denser environments, i.e. filaments, walls, and clusters. The stellar metallicity of a galaxy traces the accumulated fossil record of star formation through its entire life. Comparing the stellar metallicity of galaxies in various environments, including voids, filaments, walls, and clusters, can provide valuable insights into how the large-scale environment impacts galaxy chemical evolution. We present the first comparison of the total stellar mass vs. central stellar metallicity relation between galaxies in voids, filaments, walls, and clusters with different star formation history (SFH) types, morphologies, and colours, for stellar masses between 10^8.0 to 10^11.5 solar masses and redshift 0.01 < z < 0.05. We aim to better understand how the large-scale structure affects galaxy evolution by studying the stellar mass-metallicity relation of thousands of galaxies, which allows us to make a statistically sound comparison between galaxies in voids, filaments, walls, and clusters. We apply non-parametric full spectral fitting techniques (pPXF and STECKMAP) to 10807 spectra from the SDSS-DR7 (987 in voids, 6463 in filaments and walls, and 3357 in clusters) and derive their central mass-weighted average stellar metallicity. We find that galaxies in voids have on average slightly lower stellar metallicities than galaxies in filaments and walls (by 0.1 dex), and much lower than galaxies in clusters (by 0.4 dex). These differences are more significant for low-mass (10^9.25) than for high-mass galaxies, for long-timescale SFH (LT-SFH, extended along time) galaxies than for short-timescale SFHs (ST-SFH, concentrated at early times) galaxies, for spiral than for elliptical galaxies, and for blue than for red galaxies.

In our study, we have adopted the framework of Horava-Lifshitz gravity to model the Universe's dark matter and dark energy components. Specifically, we have considered two recent parametrizations for dark energy models: the CBDRM-type and CADMM-type parameterizations. In our analysis, we have explicitly expressed the Hubble parameter, denoted as $H(z)$, for these two distinct dark energy models. By doing so, we have aimed to investigate and quantify the accelerated cosmic expansion rate characterizing the late-time Universe. Our study uses a wide range of datasets. This dataset consists of recent measurements of baryon acoustic oscillations (BAO) collected over a period of twenty years with the Cosmic Chronometers (CC) dataset, Type Ia supernovae (SNIa) dataset, the Hubble diagram of gamma-ray bursts (GRBs), quasars (Q), and the latest measurement of the Hubble constant (R22). Consequently, we present a crucial aspect of our study by plotting the $r_{d}$ vs. $H_{0}$ plane. In the context of the $\Lambda$CDM model, after incorporating all the datasets, including the R22 prior, we obtain the following results: $H_{0}$ = $71.674089 \pm 0.734089$ $km s^{-1} Mpc^{-1}$ and $r_d = 143.050380 Mpc \pm 3.702038$. For the CBDRM model, we find $H_{0}$ = $72.355058 \pm 1.004604$ $km s^{-1} Mpc^{-1}$ and $r_d = 144.835069 Mpc \pm 2.378848$. In the case of the CADMM model, our analysis yields $H_{0}$ = $72.347804 \pm 0.923328$ $km s^{-1} Mpc^{-1}$ and $r_d = 144.466836 Mpc \pm 4.288758$. We have conducted cosmographic analyses for both of the proposed parameterizations in comparison to the $\Lambda$CDM paradigm. Additionally, we have applied Diagnostic tests to investigate the evolution of both models. Finally, the Information Criteria test demonstrates that the $\Lambda$CDM model emerges as the preferred choice among the models we have considered.

We present a study of the ratio of visible mass to total mass in spiral galaxies to better understand the relative amount of dark matter present in galaxies of different masses and evolutionary stages. Using the velocities of the H-alpha emission line measured in spectroscopic observations from the SDSS MaNGA DR17, we evaluate the rotational velocity of 5522 disk galaxies at their 90% elliptical Petrosian radii, R90. We compare this to the velocity expected from the total visible mass, which we compute from the stellar, HI, molecular hydrogen, heavy metal, and dust masses. Molecular hydrogen mass measurements are available for only a small subset of galaxies observed in SDSS MaNGA DR17, so we derive a parameterization of the molecular hydrogen mass as a function of absolute magnitude in the r-band using galaxies observed as part of SDSS DR7. With these parameterizations, we calculate the fraction of visible mass within R90 that corresponds to the observed velocity. Based on statistically analyzing the likelihood of this fraction, we conclude that the null hypothesis (no dark matter) cannot be excluded at a confidence level better than 95% within the visible extent of the disk galaxies. We also find that, by including all of these mass components, star-forming disk galaxies contain statistically the same ratio of visible-to-total mass, independent of magnitude.

In this paper, we investigate the scalar-induced gravitational waves in single-field non-attractor inflation for the Pulsar Timing Arrays data. Our model comprises three phases of inflation: the first and third phases are slow-roll inflation, while the second phase is a period of non-attractor inflation. We analyze the model's predictions for various values of the sound speed $c_s$ and examine the sharp transitions to the final attractor phase. Furthermore, we study the model's predictions for NANOGrav observations and future gravitational wave observations. We also calculate the non-Gaussianity parameter $f_{NL}$ for the non-attractor setup with a general sound speed and the sharpness parameter.

Aqueous chemistry within carbonaceous planetesimals is promising for synthesizing prebiotic organic matter essential to all life. Meteorites derived from these planetesimals delivered these life building blocks to the early Earth, potentially facilitating the origins of life. Here, we studied the formation of vitamin B$_3$ as it is an important precursor of the coenzyme NAD(P)(H), which is essential for the metabolism of all life as we know it. We propose a new reaction mechanism based on known experiments in the literature that explains the synthesis of vitamin B$_3$. It combines the sugar precursors glyceraldehyde or dihydroxyacetone with the amino acids aspartic acid or asparagine in aqueous solution without oxygen or other oxidizing agents. We performed thermochemical equilibrium calculations to test the thermodynamic favorability. The predicted vitamin B$_3$ abundances resulting from this new pathway were compared with measured values in asteroids and meteorites. We conclude that competition for reactants and decomposition by hydrolysis are necessary to explain the prebiotic content of meteorites. In sum, our model fits well into the complex network of chemical pathways active in this environment.

The recurring transient outbursts in low-mass X-ray binaries (LMXBs) provide ideal laboratories to study the accretion process. Unlike their supermassive relatives, LMXBs are far too small and distant to be imaged directly. Fortunately, phase-resolved spectroscopy can provide an alternative diagnostic to study their highly complex, time-dependent accretion discs. The primary spectral signature of LMXBs are strong, disc-formed emission lines detected at optical wavelengths. The shape, profile, and appearance/disappearance of these lines change throughout a binary orbit, and thus, can be used to trace how matter in these discs behaves and evolves over time. By combining a \textit{Swift} multi-wavelength monitoring campaign, phase-resolved spectroscopy from the Gran Telescopio Canarias (GTC) and Liverpool Telescope, and modern astrotomography techniques, we find a clear empirical connection between the line emitting regions and physical properties of the X-rays heating the disc in the black hole LMXB MAXI J1820+070 during its 2018 outburst. In this paper, we show how these empirical correlations can be used as an effective observational tool for understanding the geometry and structure of a LMXB accretion disc and present further evidence for an irradiation-driven warped accretion disc present in this system.

This proceeding is an introduction to cosmological applications of the Lorentz gauge theory. It provides the ingredients for a unique, though yet tentative $\Lambda$CDM theory of cosmology. The emergence of spacetime is described by the spontaneous symmetry breaking called here the khronogenesis. Space is then associated with the field strength of the antiself-dual gauge potential, and gravity is associated with the self-dual field strength. In the cosmological setting, khronogenesis seems to predict inflation. It is shown that the Lorentz gauge theory allows the consistent description of spin currents which could have important roles in cosmological phenomenology.

General Relativity (GR) yields a number of gravitational wave memory effects that correspond to symmetries of spacetime infinitely far away from gravitational fields. These symmetries and memory effects hint at the fundamental mathematical connection between gravity and quantum fields in the low-energy "infrared" regime. In this study, we propose to shift the paradigm of memory from merely a target in gravitational wave searches to a unique way of measuring symmetries of nature. Thus, we extend previous memory detection efforts to the proof-of-principle parameter estimation and model selection. Through simulating binary black hole (BBH) mergers in scenarios where spacetime exhibits certain sets of Bondi-Metzner-Sachs symmetries, we point out that both the current and the future gravitational wave observatories are excellent probes of spacetime symmetries. In particular, the design sensitivity of the proposed Einstein Telescope (ET) allows to constrain the strain amplitude of the leading-order displacement memory, associated with superrotational symmetries, to a 2% level in one year. Whereas the weaker spin memory amplitude can be constrained to a 22% level, providing a pathway to probe superrotational symmetries. Finally, there is almost no doubt among GR experts that displacement memory exists in nature, although its effect on parameter estimation was largely neglected in the predictions of the science output of future experiments such as LISA and ET. We find that it may lead to an overestimation of the measurement uncertainties for inferred parameters of the loudest BBH mergers by the order of 10%.

Neutrinos -- amongst the lightest known particles -- can mediate a force driving dark matter self-interaction and the small scale structure of the universe. We explore such a possibility in the simplest neutrino portal dark sector model where neutrino has a Yukawa coupling with a scalar $\phi$ and fermion $\chi$ that are degenerate in mass and together comprise 100\% of dark matter in the universe. We derive the non-relativistic potential generated by single-neutrino exchange which is of the monopole-dipole form and explore $\chi\phi\to\chi\phi$ scattering based on phase shift analysis. Our result shows that Born approximation continues to be valid in the low energy regime and the scattering cross section scales as $1/v^2$ over a wide range of dark matter velocities. Such a velocity-dependent self-interacting cross section can be large enough to explain the shallow density of dwarf galaxy cores and consistent with the upper limit from colliding galaxy clusters. The $1/v^2$ behavior persists down to rather low velocities $v\sim m_\nu/m$ where $m$ is the dark matter mass, leaving the opportunity for further astrophysical probes. Through the neutrino portal, the self-interacting dark matter parameter space can be tested by searches for $Z$-boson and light mesons decaying into the dark sector, as well as low-mass dark matter direct detection experiments.

We consider the optical appearance under a thin accretion disk of a regular black hole with a central de Sitter core implementing $\mathcal{O}(l^2/r^2)$ far-corrections to the Schwarzschild black hole. We use the choice $l=0.25M$, which satisfies recently found constraints from the motion of the S2 star around Sgr A$^*$ in this model, and which leads to thermodynamically stable black holes. As the emission model, we suitably adapt ten samples of the Standard Unbound emission profile for a monochromatic intensity in the disk's frame, which have been previously employed in the literature within the context of reproducing General Relativistic Magneto-Hydrodynamic simulations of the accretion flow. We find the usual central brightness depression surrounded by the bright ring cast by the disk's direct emission as well as two non-negligible photon ring contributions. As compared to the usual Schwarzschild solution, the relative luminosities of the latter are significantly boosted, while the size of the former is strongly decreased. We discuss the entanglement of the background geometry and the choice of emission model in generating these black hole images, as well as the capability of these modifications of Schwarzschild solution to pass present and future tests based on their optical appearance when illuminated by an accretion disk.

We present a calculation of the continuum part of the gamma-ray spectra resulting from Dark Matter annihilation in the framework of the MSSM taking into account Sommerfeld effects. Concentrating on pure wino and pure higgsino scenarios we compare our calculation to existing work and explore the numerical impact of the features not captured by previous approximative descriptions. We find that, in particular for large neutralino masses, when the Sommerfeld enhancement is very large, chargino-antichargino annihilation processes, which have not been considered before, lead to sizable differences with respect to existing calculations. In scenarios with neutralinos in the intermediate-mass range, we find that the role of the charginos is crucial in the endpoint regime. Our calculation provides the currently most accurate prediction for the continuum gamma-ray spectra.

The focus of dark matter searches to date has been on Weakly Interacting Massive Particles (WIMPs) in the GeV/$c^2$-TeV/$c^2$ mass range. The direct, indirect and collider searches in this mass range have been extensive but ultimately unsuccessful, providing a strong motivation for widening the search outside this range. Here we describe a new concept for a dark matter experiment, employing superfluid $^3$He as a detector for dark matter that is close to the mass of the proton, of order 1 GeV/$c^2$. The QUEST-DMC detector concept is based on quasiparticle detection in a bolometer cell by a nanomechanical resonator. In this paper we develop the energy measurement methodology and detector response model, simulate candidate dark matter signals and expected background interactions, and calculate the sensitivity of such a detector. We project that such a detector can reach sub-eV nuclear recoil energy threshold, opening up new windows on the parameter space of both spin-dependent and spin-independent interactions of light dark matter candidates.

The large-scale magnetic structure of interplanetary coronal mass ejections (ICMEs) has been shown to affect the galactic cosmic ray (GCR) flux measured in situ by spacecraft, causing temporary decreases known as Forbush decreases (Fds). In some ICMEs, the magnetic ejecta exhibits a magnetic flux rope (FR) structure; the strong magnetic field strength and closed field line geometry of such ICME FRs has been proposed to act as a shield to GCR transport. In this study, we identify four ICMEs near Earth that drove Fds with similar mean magnetic field strengths (20 - 25 nT); two ICMEs with more typical mean speeds (~400 km/s), and two fast (~750 km/s) ICMEs. Within each speed pairing, we identify an ICME that exhibited an open magnetic field line topology and compare its effect on the GCR flux to that which exhibited a mostly closed topology. We investigate the different mechanisms that contribute to the resulting ICME-related Fds and their recovery, and determine which properties, if any, play a more important role than others in driving Fds. We find that much of the GCR response to the ICME events in this study is independent of the open or closed magnetic field line topology of the flux rope, and that features such as the fluctuations in speed, magnetic field structure, and expansion within the FR may play more of a role in determining the smaller-scale structure of the Fd profile.

How thermal particles are accelerated to suprathermal energies is an unsolved issue, crucial for many astrophysical systems. We report novel observations of irregular, dispersive enhancements of the suprathermal particle population upstream of a high-Mach number interplanetary shock. We interpret the observed behavior as irregular "injections" of suprathermal particles resulting from shock front irregularities. Our findings, directly compared to self-consistent simulation results, provide important insights for the study of remote astrophysical systems where shock structuring is often neglected.