Papers for Thursday, Jun 22 2023

Using the Hubble Ultraviolet Globular Cluster Survey (HUGS) and additional HST archival data, we have carried out a search for optical counterparts to the low-luminosity Chandra X-ray sources in the globular cluster M4 (NGC 6121). We have also searched for optical or X-ray counterparts to radio sources detected by the VLA. We find 24 new confident optical counterparts to Chandra sources for a total of 40, including the 16 previously identified. Of the 24 new identifications, 18 are stellar coronal X-ray sources (active binaries, ABs), the majority located along the binary sequence in a V-I colour-magnitude diagram and generally showing an H-alpha excess. In addition to confirming the previously detected cataclysmic variable (CV, CX4), we identify one confident new CV (CX76), and two candidates (CX81 and CX101). One MSP is known in M4 (CX12), and another strong candidate has been suggested (CX1); we identify some possible MSP candidates among optical and radio sources, such as VLA20, which appears to have a white dwarf counterpart. One X-ray source with a sub-subgiant optical counterpart and a flat radio spectrum (CX8, VLA31) is particularly mysterious. The radial distribution of X-ray sources suggests a relaxed population of average mass ~ 1.2 - 1.5 Msun. Comparing the numbers of ABs, MSPs, and CVs in M4 with other clusters indicates that AB numbers are proportional to cluster mass (primordial population), MSPs to stellar encounter rate (dynamically formed population), while CVs seem to be produced both primordially and dynamically.

Time-delay interferometry (TDI) is a crucial step in the on-ground data processing pipeline of the Laser Interferometer Space Antenna (LISA), as it reduces otherwise overwhelming laser noise and allows for the detection of gravitational waves. This being said, several laser noise couplings have been identified that limit the performance of TDI. First, on-board processing, which is used to decimate the sampling rate from tens of MHz down to a few Hz, requires careful design of the anti-aliasing filters to mitigate folding of laser noise power into the observation band. Furthermore, the flatness of those filters is important to limit the effect of the flexing-filtering coupling. Secondly, the post-processing delays applied in TDI are subject to ranging and interpolation errors. All of these effects are partially described in the literature. In this paper, we present them in a unified framework and give a more complete description of aliased laser noise and the coupling of interpolation errors. Furthermore, for the first time, we discuss the impact of laser locking on laser noise residuals in the final TDI output. To verify the validity of the analytic PSD models we derive, we run numerical simulations using LISA Instrument and calculate second-generation TDI variables with PyTDI. We consider a setup with six independent lasers and with locked lasers (locking configuration N1-12). We find that laser locking indeed affects the laser noise residual in the TDI combination as it introduces correlations among the six lasers inducing slight modulations of the PSD compared to the case of six independent lasers. This implies further studies on laser noise residuals should consider the various locking configurations to produce accurate results.

We study the motion of neutron superfluid vortices in a spinning-down neutron star, assuming axisymmetry of the flow and ignoring motion of vortices about the rotation axis. We find that the vortex array, if initially rectilinear, is soon substantially deformed as the star spins down; vortices are swept outward by the Magnus force, accumulating in regions of the inner crust where they pin, accompanied by significant bending of the vortex array. As the star spins down to below a spin rate of ~20 Hz (twice the spin rate of the Vela pulsar), the Magnus and pinning forces gradually compress the vortex array into dense sheets that follow spherical shells. In some cases, the vortex array bends on itself and reconnects, forming one or more tori of vortex rings that contain superfluid ``rivers" with significant angular momentum. Vortex sheets are likely to form near the base of the inner crust, in the regime of nuclear pasta.

Research efforts that require observations of high solar activity, such as multiwavelength studies of large solar flares and CMEs, must contend with the 11-year solar cycle to a degree unparalleled by other segments of heliophysics. While the "fallow" years around each solar minimum can be a great time frame to build the next major solar observatory, the corresponding funding opportunity and any preceding technology developments would need to be strategically timed. Even then, it can be challenging for scientists on soft money to continue ongoing research efforts instead of switching to other, more consistent topics. The maximum of solar cycle 25 is particularly concerning due to the lack of a US-led major mission targeting high solar activity, which could result in significant attrition of expertise in the field. We recommend the development of a strategic program of missions and analysis that ensures optimal science return for each solar maximum while sustaining the research community between maxima.

It is essential that there be coordinated and co-optimized observations in X-rays, gamma-rays, and EUV during the peak of solar cycle 26 (~2036) to significantly advance our understanding of impulsive energy release in the corona. The open questions include: What are the physical origins of space-weather events? How are particles accelerated at the Sun? How is impulsively released energy transported throughout the solar atmosphere? How is the solar corona heated? Many of the processes involved in triggering, driving, and sustaining solar eruptive events -- including magnetic reconnection, particle acceleration, plasma heating, and energy transport in magnetized plasmas -- also play important roles in phenomena throughout the Universe. This set of observations can be achieved through a single flagship mission or, with foreplanning, through a combination of major missions (e.g., the previously proposed FIERCE mission concept).

Earth is deficient in carbon and nitrogen by up to ${\sim}4$ orders of magnitude compared with the Sun. Destruction of (carbon- and nitrogen-rich) refractory organics in the high-temperature planet forming regions could explain this deficiency. Assuming a refractory cometary composition for these grains, their destruction enhances nitrogen-containing oxygen-poor molecules in the hot gas ($\gtrsim 300$K) after the initial formation and sublimation of these molecules from oxygen-rich ices in the warm gas (${\sim}150$K). Using observations of $37$ high-mass protostars with ALMA, we find that oxygen-containing molecules (CH$_3$OH and HNCO) systematically show no enhancement in their hot component. In contrast, nitrogen-containing, oxygen-poor molecules (CH$_3$CN and C$_2$H$_3$CN) systematically show an enhancement of a factor ${\sim} 5$ in their hot component, pointing to additional production of these molecules in the hot gas. Assuming only thermal excitation conditions, we interpret these results as a signature of destruction of refractory organics, consistent with the cometary composition. This destruction implies a higher C/O and N/O in the hot gas than the warm gas, while, the exact values of these ratios depend on the fraction of grains that are effectively destroyed. This fraction can be found by future chemical models that constrain C/O and N/O from the abundances of minor carbon, nitrogen and oxygen carriers presented here.

Abell 3827 is a unique galaxy cluster with a dry merger in its core causing a highly-resolved multiple-image configuration of a blue spiral galaxy at $z_\mathrm{s}=1.24$. The surface brightness profiles of four merging galaxies around $z_\mathrm{d}=0.099$ complicate a clear identification of the number of images and finding corresponding small-scale features across them. The entailed controversies about offsets between luminous and dark matter have never been settled and dark-matter characteristics in tension with bounds from complementary probes and simulations seemed necessary to explain this multiple-image configuration. We resolve these issues with a systematic study of possible feature matchings across all images and their impact on the reconstructed mass density distribution. From the local lens properties directly constrained by these feature matchings without imposing any global lens model, we conclude that none of them are consistent with expected local characteristics from standard single-lens-plane lensing, nor can they be motivated by the light distribution in the cluster. Inspecting complementary spectroscopic data, we show that all these results originate from an insufficient constraining power of the data and seem to hint at a thick lens and not at exotic forms of dark matter or modified gravity. If the thick-lens hypothesis can be corroborated with follow-up multi-plane lens modelling, A3827 suffers from a full three-dimensional degeneracy in the distribution of dark matter because combinations of shearings and scalings in a single lens plane can also be represented by an effective shearing and a rotation caused by multiple lens planes.

Dark matter structures within strong gravitational lens galaxies and along their line of sight leave a gravitational imprint on the multiple images of lensed sources. Strong gravitational lensing provides, therefore, a key test of different dark matter models in a way that is independent of the baryonic content of matter structures on subgalactic scales. In this chapter, we describe how galaxy-scale strong gravitational lensing observations are sensitive to the physical nature of dark matter. We provide a historical perspective of the field, and review its current status. We discuss the challenges and advances in terms of data, treatment of systematic errors and theoretical predictions, that will enable one to deliver a stringent and robust test of different dark matter models in the near future. With the advent of the next generation of sky surveys, the number of known strong gravitational lens systems is expected to increase by several orders of magnitude. Coupled with high-resolution follow-up observations, these data will provide a key opportunity to constrain the properties of dark matter with strong gravitational lensing.

What would a parity-violating universe look like? We present a numerical and theoretical study of mirror asymmetries in the late universe, using a new suite of $N$-body simulations: QUIJOTE-Odd. These feature parity-violating initial conditions, injected via a simple ansatz for the imaginary primordial trispectrum and evolved into the non-linear regime. We find that the realization-averaged power spectrum, bispectrum, halo mass function, and matter PDF are not affected by our modifications to the initial conditions, deep into the non-linear regime, which we argue arises from rotational and translational invariance. In contrast, the parity-odd trispectrum of matter (measured using a new estimator), shows distinct signatures proportional to the parity-violating parameter, $p_{\rm NL}$, which sets the amplitude of the primordial trispectrum. We additionally find intriguing signatures in the angular momentum of halos, with the primordial trispectrum inducing a non-zero correlation between angular momentum and smoothed velocity field, proportional to $p_{\rm NL}$. Our simulation suite has been made public to facilitate future analyses.

The Nancy Grace Roman Space Telescope is capable of delivering an unprecedented all-sky, high-spatial resolution, multi-epoch infrared map to the astronomical community. This opportunity arises in the midst of numerous ground- and space-based surveys that will provide extensive spectroscopy and imaging together covering the entire sky (such as Rubin/LSST, Euclid, UNIONS, SPHEREx, DESI, SDSS-V, GALAH, 4MOST, WEAVE, MOONS, PFS, UVEX, NEO Surveyor, etc.). Roman can uniquely provide uniform high-spatial-resolution (~0.1 arcsec) imaging over the entire sky, vastly expanding the science reach and precision of all of these near-term and future surveys. This imaging will not only enhance other surveys, but also facilitate completely new science. By imaging the full sky over two epochs, Roman can measure the proper motions for stars across the entire Milky Way, probing 100 times fainter than Gaia out to the very edge of the Galaxy. Here, we propose NANCY: a completely public, all-sky survey that will create a high-value legacy dataset benefiting innumerable ongoing and forthcoming studies of the universe. NANCY is a pure expression of Roman's potential: it images the entire sky, at high spatial resolution, in a broad infrared bandpass that collects as many photons as possible. The majority of all ongoing astronomical surveys would benefit from incorporating observations of NANCY into their analyses, whether these surveys focus on nearby stars, the Milky Way, near-field cosmology, or the broader universe.

Very massive stars (VMS) up to 200-300 $M_\odot$ have been found in the Local Universe. If they would lose little mass they produce intermediate-mass black holes or pair-instability supernovae (PISNe). Until now, VMS modellers have extrapolated mass-loss vs. metallicity ($Z$) exponents from optically-thin winds, resulting in a range of PISN thresholds that might be unrealistically high in $Z$, as VMS develop optically-thick winds. We utilize the transition mass-loss rate of Vink and Gr\"afener (2012) that accurately predicts mass-loss rates of Of/WNh ("slash") stars that characterize the morphological transition from absorption-dominated O-type spectra to emission-dominated WNh spectra. We develop a wind efficiency framework, where optically thin winds transition to enhanced winds, enabling us to study VMS evolution at high redshift where individual stars cannot be resolved. We present a MESA grid covering $Z_\odot/2$ to $Z_\odot/100$. VMS above the transition evolve towards lower luminosity, skipping the cool supergiant phase but directly forming pure He stars at the end of hydrogen burning. Below the transition, VMS evolve as cooler luminous blue variables (LBVs) or yellow hypergiants (YHGs), naturally approaching the Eddington limit. Strong winds in this YHG/LBV regime -- combined with a degeneracy in luminosity -- result in a mass-loss runaway where a decrease in mass increases wind mass loss. Our models indicate an order-of-magnitude lower threshold than usually assumed, at $Z_\odot/20$ due to our mass-loss runaway. While future work on LBV mass loss could affect the PISN threshold, our framework will be critical for establishing definitive answers on the PISN threshold and galactic chemical evolution modelling.

As the nearest accessible massive early-type galaxy, NGC 5128 presents an exceptional opportunity to measure dark matter halo parameters for a representative elliptical galaxy. Here we take advantage of rich new observational datasets of large-radius tracers to perform dynamical modeling of NGC 5128, using a discrete axisymmetric anisotropic Jeans approach with a total tracer population of nearly 1800 planetary nebulae, globular clusters, and dwarf satellite galaxies extending to a projected distance of $\sim250$ kpc from the galaxy center. We find that a standard NFW halo provides an excellent fit to nearly all the data, excepting a subset of the planetary nebulae that appear to be out of virial equilibrium. The best-fit dark matter halo has a virial mass of ${\rm M}_{vir}=4.4^{+2.4}_{-1.4}\times10^{12} {\rm M}_{\odot}$, and NGC 5128 appears to sit below the mean stellar mass--halo mass and globular cluster mass--halo mass relations, which both predict a halo virial mass closer to ${\rm M}_{vir} \sim 10^{13} {\rm M}_{\odot}$. The inferred NFW virial concentration is $c_{vir}=5.6^{+2.4}_{-1.6}$, nominally lower than $c_{vir} \sim 9$ predicted from published $c_{vir}$--${\rm M}_{vir}$ relations, but within the $\sim 30\%$ scatter found in simulations. The best-fit dark matter halo constitutes only $\sim10\%$ of the total mass at 1 effective radius but $\sim50\%$ at 5 effective radii. The derived halo parameters are relatively insensitive to reasonable variations in the tracer population considered, tracer anisotropies, and system inclination. Our analysis highlights the value of comprehensive dynamical modeling of nearby galaxies, and the importance of using multiple tracers to allow cross-checks for model robustness.

We present nebular spectroscopy of SN 2020hvf, a Type Ia supernova (SN Ia) with an early bump in its light curve. SN 2020hvf shares many spectroscopic and photometric similarities to the carbon-rich high-luminosity "03fg-like" SNe Ia. At $>$240 days after peak brightness, we detect unambiguous emission from [Ca II] $\lambda\lambda$7291, 7324 which is never observed in normal-SNe Ia and only seen in peculiar subclasses. SN 2020hvf displays "saw-tooth" emission profiles near 7300 A that cannot be explained with single symmetric velocity components of [Fe II], [Ni II], and [Ca II], indicating an asymmetric explosion. The broad [Ca II] emission is best modeled by two velocity components offset by 1,220 km s$^{-1}$, which could be caused by ejecta associated with each star in the progenitor system, separated by their orbital velocity. For the first time in a SN Ia, we identify narrow (${\rm FWHM} = 180\pm40$ km s$^{-1}$) [Ca II] emission, which we associate with a wind from a surviving, puffed-up companion star. Few published spectra have sufficient resolution and signal-to-noise ratio necessary to detect similar narrow [Ca II] emission, however, we have detected similar line profiles in other 03fg-like SNe Ia. The extremely narrow velocity width of [Ca II] has only otherwise been observed in SNe Iax at late times. Since this event likely had a double-degenerate "super-Chandrasekhar" mass progenitor system, we suggest that a single white dwarf (WD) was fully disrupted and a wind from a surviving companion WD is producing the observed narrow emission. It is unclear if this unique progenitor and explosion scenario can explain the diversity of 03fg-like SNe Ia, potentially indicating that multiple progenitor channels contribute to this subclass.

Star-forming galaxies like the Milky Way are surrounded by a hot gaseous halo at the virial temperature - the so-called galactic corona - that plays a fundamental role in their evolution. The interaction between the disc and the corona has been shown to have a direct impact on accretion of coronal gas onto the disc with major implications for galaxy evolution. In this work, we study the gas circulation between the disc and the corona of star-forming galaxies like the Milky Way. We use high-resolution hydrodynamical N-body simulations of a Milky Way-like galaxy with the inclusion of an observationally-motivated galactic corona. In doing so, we use SMUGGLE, an explicit interstellar medium (ISM) and stellar feedback model coupled with the moving-mesh code Arepo. We find that the reservoir of gas in the galactic corona is sustaining star formation: the gas accreted from the corona is the primary fuel for the formation of new stars, helping in maintaining a nearly constant level of cold gas mass in the galactic disc. Stellar feedback generates a gas circulation between the disc and the corona (the so-called galactic fountain) by ejecting different gas phases that are eventually re-accreted onto the disc. The accretion of coronal gas is promoted by its mixing with the galactic fountains at the disc-corona interface, causing the formation of intermediate temperature gas that enhance the cooling of the hot corona. We find that this process acts as a positive feedback mechanism, increasing the accretion rate of coronal gas onto the galaxy.

Radio-mode feedback is a key ingredient in galaxy formation and evolution models, required to reproduce the observed properties of massive galaxies in the local Universe. We study the cosmic evolution of radio-AGN feedback out to $z\sim2.5$ using a sample of 9485 radio-excess AGN. We combine the evolving radio luminosity functions with a radio luminosity scaling relationship to estimate AGN jet kinetic powers and derive the cosmic evolution of the kinetic luminosity density, $\Omega_{\rm{kin}}$ (i.e. the volume-averaged heating output). Compared to all radio-AGN, low-excitation radio galaxies (LERGs) dominate the feedback activity out to $z\sim2.5$, with both these populations showing a constant heating output of $\Omega_{\rm{kin}} \approx 4-5 \times 10^{32}\,\rm{W\,Mpc^{-3}}$ across $0.5 < z < 2.5$. We compare our observations to predictions from semi-analytical and hydrodynamical simulations, which broadly match the observed evolution in $\Omega_{\rm{kin}}$, although their absolute normalisation varies. Comparison to the Semi-Analytic Galaxy Evolution (SAGE) model suggests that radio-AGN may provide sufficient heating to offset radiative cooling losses, providing evidence for a self-regulated AGN feedback cycle. We integrate the kinetic luminosity density across cosmic time to obtain the kinetic energy density output from AGN jets throughout cosmic history to be $\sim 10^{50}\,\rm{J\,Mpc^{-3}}$. Compared to AGN winds, the kinetic energy density from AGN jets dominates the energy budget at $z \lesssim 2$; this suggests that AGN jets play an important role in AGN feedback across most of cosmic history.

Massive stars are progenitors of supernovae, neutron stars and black holes. During the hydrogen-core burning phase their convective cores are the prime drivers of their evolution, but inferences of core masses are subject to unconstrained boundary mixing processes. Moreover, uncalibrated transport mechanisms can lead to strong envelope mixing and differential radial rotation. Ascertaining the efficiency of the transport mechanisms is challenging because of a lack of observational constraints. Here we deduce the convective core mass and robustly demonstrate non-rigid radial rotation in a supernova progenitor, the $12.0^{+1.5}_{-1.5}$ solar-mass hydrogen-burning star HD 192575, using asteroseismology, TESS photometry, high-resolution spectroscopy, and Gaia astrometry. We infer a convective core mass ($M_{\rm cc} = 2.9^{+0.5}_{-0.8}$ solar masses), and find the core to be rotating between 1.4 and 6.3 times faster than the stellar envelope depending on the location of the rotational shear layer. Our results deliver a robust inferred core mass of a massive star using asteroseismology from space-based photometry. HD 192575 is a unique anchor point for studying interior rotation and mixing processes, and thus also angular momentum transport mechanisms inside massive stars.

We investigate the incidence and properties of ionized gas outflows in a sample of 52 galaxies with stellar mass between $10^7$ M$_{\odot}$ and $10^9$ M$_{\odot}$ observed with ultra-deep JWST/NIRSpec MSA spectroscopy as part of the JWST Advanced Deep Extragalactic Survey (JADES). The high-spectral resolution (R2700) NIRSpec observations allowed us to identify for the first time the signature of outflows in the rest-frame optical nebular lines in low-mass galaxies at $z>3$. The incidence fraction of ionized outflows, traced by broad components, is about 25-40$\%$ depending on the intensity of the emission lines. The low incidence fraction might be due to both the sensitivity limit and the fact that outflows are not isotropic but have a limited opening angle which results in a detection only when this is directed toward our line of sight. Evidence for outflows increases slightly with stellar mass and star-formation rate. The median velocity and mass loading factor (i.e., the ratio between mass outflow rate and star formation rate) of the outflowing ionized gas are 500 km s$^{-1}$ and $\eta=2.1^{+2.5}_{-1.6}$, respectively. These are two and 100 times higher, respectively than the typical values observed in local dwarf galaxies. These outflows are able to escape the gravitational potential of the galaxy and enrich the circum-galactic medium and, potentially, the inter-galactic medium. Our results indicate that outflows can significantly impact the star formation activity in low-mass galaxies within the first 2 Gyr of the Universe.

The origin of astrophysical neutrinos is one of the most debated topics today. Perhaps the most robust evidence of neutrino counterpart comes from supermassive black holes in active galactic nuclei associated with strongly collimated outflows, or jets, that can accelerate particles to relativistic energies and produce neutrinos through hadronic interactions. Similar outflows can also be found from X-ray binaries, or `microquasars', that consist of a neutron star or a stellar-mass black hole accreting matter from a non-degenerate companion star. In some cases, these systems can accelerate particles up to GeV energies implying an efficient acceleration mechanism in their jets. Neutrino production in microquasar jets can be expected with suitable conditions and a hadronic particle population. Microquasar Cyg X-3 is a unique, short orbital period X-ray binary hosting a Wolf-Rayet companion star with a strong stellar wind. The interaction of the dense stellar wind with a relativistic jet leads to particle collisions followed by high-energy gamma-ray and potentially neutrino emission. Here, using the 10-year neutrino candidate sample of the IceCube neutrino observatory, we find that the events with the highest spatial association with Cyg X-3 occur during short-lived high-energy gamma-ray flaring periods indicating the possible astrophysical nature of these events.

White dwarfs stars are known to be polluted by their active planetary systems, but little attention has been paid to the accretion of wind from low-mass companions. The capture of stellar or substellar wind by white dwarfs is one of few methods available to astronomers which can assess mass-loss rates from unevolved stars and brown dwarfs, and the only known method to extract their chemical compositions. In this work, four white dwarfs with closely-orbiting, L-type brown dwarf companions are studied to place limits on the accretion of a substellar wind, with one case of a detection, and at an extremely non-solar abundance $m_{\rm Na}/m_{\rm Ca}>900$. The mass-loss rates and upper limits are tied to accretion in the white dwarfs, based on limiting cases for how the wind is captured, and compared with known cases of wind pollution from close M dwarf companions, which manifest in solar proportions between all elements detected. For wind captured in a Bondi-Hoyle flow, mass-loss limits $\dot M\lesssim 5\times10^{-17}$ M$_\odot$ yr$^{-1}$ are established for three L dwarfs, while for M dwarfs polluting their hosts, winds in the range $10^{-13} - 10^{-16}$ M$_\odot$ yr$^{-1}$ are found. The latter compares well with the $\dot M\sim 10^{-13} - 10^{-15}$ M$_\odot$ yr$^{-1}$ estimates obtained for nearby, isolated M dwarfs using Ly$\alpha$ to probe their astropsheres. These results demonstrate that white dwarfs are highly-sensitive stellar and substellar wind detectors, where further work on the actual captured wind flow is needed.

We present the first JWST spectral energy distribution of a Y dwarf. This spectral energy distribution of the Y0 dwarf WISE J035934.06$-$540154.6 consists of low-resolution ($\lambda$/$\Delta\lambda$ $\sim$ 100) spectroscopy from 1$-$12 $\mu$m and three photometric points at 15, 18, and 21 $\mu$m. The spectrum exhibits numerous fundamental, overtone, and combination rotational-vibrational bands of H$_2$O, CH$_4$, CO, CO$_2$, and NH$_3$, including the previously unidentified $\nu_3$ band of NH$_3$ at 3 $\mu$m. Using a Rayleigh-Jeans tail to account for the flux emerging at wavelengths greater than 21 $\mu$m, we measure a bolometric luminosity of $1.523\pm0.090\times10^{20}$ W. We determine a semi-empirical effective temperature estimate of $467^{+16}_{-18}$ K using the bolometric luminosity and evolutionary models to estimate a radius. Finally, we compare the spectrum and photometry to a grid of atmospheric models and find reasonably good agreement with a model having $T_{\mathrm{eff}}$=450 K, log $g$=3.25 [cm s$^{-2}$], [M/H]=$-0.3$. However, the low surface gravity implies an extremely low mass of 1 $M_{\rm{Jup}}$ and a very young age of 20 Myr, the latter of which is inconsistent with simulations of volume-limited samples of cool brown dwarfs.

Recent observations demonstrate that misalignments and other out-of-plane structures are common in protoplanetary discs. Many of these have been linked to a central host binary with an orbit that is inclined with respect to the disc. We present simulations of misaligned circumbinary discs with a range of parameters to gain a better understanding of the link between those parameters and the disc morphology in the wave-like regime of warp propagation that is appropriate to protoplanetary discs. The simulations confirm that disc tearing is possible in protoplanetary discs as long as the mass ratio, $\mu$, and disc-binary inclination angle, $i$, are not too small. For the simulations presented here this corresponds to $\mu > 0.1$ and $i \gtrsim 40^\circ$. For highly eccentric binaries, tearing can occur for discs with smaller misalignment. Existing theoretical predictions provide an estimate of the radial extent of the disc in which we can expect breaking to occur. However, there does not seem to be a simple relationship between the disc properties and the radius within the circumbinary disc at which the breaks appear, and furthermore the radius at which the disc breaks can change as a function of time in each case. We discuss the implications of our results for interpreting observations and suggest some considerations for modelling misaligned discs in the future.

We study the formation of spinning primordial black holes during an early matter-dominated era. Using non-linear 3+1D general relativistic simulations, we compute the efficiency of mass and angular momentum transfer in the process -- which we find to be $\mathcal{O}(10\%)$ and $\mathcal{O}(5\%)$, respectively. We show that subsequent evolution is important due to the seed PBH accreting non-rotating matter from the background, which decreases the dimensionless spin. Unless the matter era is short, we argue that the final dimensionless spin will be negligible.

Measurements of the star formation activity on cloud scales are fundamental to uncovering the physics of the molecular cloud, star formation, and stellar feedback cycle in galaxies. Infrared (IR) emission from small dust grains and polycyclic aromatic hydrocarbons (PAHs) are widely used to trace the obscured component of star formation. However, the relation between these emission features and dust attenuation is complicated by the combined effects of dust heating from old stellar populations and an uncertain dust geometry with respect to heating sources. We use images obtained with NIRCam and MIRI as part of the PHANGS--JWST survey to calibrate dust emission at 21$\rm \mu m$, and the emission in the PAH-tracing bands at 3.3, 7.7, 10, and 11.3$\rm \mu m$ as tracers of obscured star formation. We analyse $\sim$ 16000 optically selected HII regions across 16 nearby star-forming galaxies, and benchmark their IR emission against dust attenuation measured from the Balmer decrement. We model the extinction-corrected H$\alpha$ flux as the sum of the observed H$\alpha$ emission and a term proportional to the IR emission, with $a_{IR}$ as the proportionality coefficient. A constant $a_{IR}$ leads to extinction-corrected H$\alpha$ estimates which agree with those obtained with the Balmer decrement with a scatter of $\sim$ 0.1 dex for all bands considered. Among these bands, 21$\rm \mu m$ emission is demonstrated to be the best tracer of dust attenuation. The PAH-tracing bands underestimate the correction for bright HII regions, since in these environments the ratio of PAH-tracing bands to 21$\rm \mu m$ decreases, signalling destruction of the PAH molecules. For fainter HII regions all bands suffer from an increasing contamination from the diffuse infrared background.

JWST observations of the strongly lensed galaxy The Sparkler have revealed a population of gravitationally bound globular cluster (GC) candidates. Different analyses have resulted in broadly similar ages but significantly different metallicities, questioning the assembly history that has led to the formation of such a population. In this letter, we re-analyse the two sets of photometry available in the literature with the code MCMAME especially tailored to fit physical properties of GCs. We find the ages and metallicities from both datasets are consistent within 1 $\sigma$ uncertainties. A significant group of GCs is consistent with being old and metal poor ([Fe/H] $\sim -1.7$). For this group, the ages do not converge, hence, we conclude that they are definitively older than 1 Gyr and can be as old as the age of the Universe. The remaining GCs have younger ages and a metallicity spread. The ages and metallicities distribution of GCs in the Sparkler are consistent with those observed in Local Group's galaxies at similar lookback times. Comparing with predictions from E-MOSAICS simulations we confirm that the Sparkler GC population traces the self-enrichment history of a galaxy which might become a few times $10^9$ M$_{\odot}$ massive system at redshift $z = 0$

The recently discovered transiting super-Earth LP 890-9 c is potentially one of the best rocky exoplanets for atmospheric characterization. Orbiting an ultracool M-dwarf at the inner edge of the habitable zone, LP 890-9 c offers a new opportunity to study the climate of rocky planets at the inner edge of the habitable zone. We investigate the molecular detectability with simulated JWST transmission spectra for five potential atmospheres of LP 890-9 c. We find that a small three-transit JWST program can infer evidence of H2O (at 3$\sigma$ confidence) for a full runaway greenhouse scenario. Alternatively, CO2-dominated atmospheres resembling Venus without high-altitude terminator clouds can be identified with eight transits. However, these predictions could be complicated by the impact of clouds and/or unocculted starspots. Nevertheless, JWST observations of LP 890-9 c could provide critical insights and potentially distinguish between models of rocky planets at the inner edge of the habitable zone.

The study of light lensed by cosmic matter has yielded much information about astrophysical questions. Observations are explained using geometrical optics following a ray-based description of light. After deflection the lensed light interferes, but observing this diffractive aspect of gravitational lensing has not been possible due to coherency challenges caused by the finite size of the sources or lack of near-perfect alignment. In this article, we report on the observation of these wave effects of gravitational lensing by recreating the lensing conditions in the laboratory via electro-optic deflection of coherent laser light. The lensed light produces a beam containing regularities, caustics, and chromatic modulations of intensity that depend on the symmetry and structure of the lensing object. We were also able to observe previous and new geometric-optical lensing situations that can be compared to astrophysical observations. This platform could be a useful tool for testing numerical/analytical simulations, and for performing analog simulations of lensing situations when they are difficult to obtain otherwise. We found that laboratory lensed beams constitute a new class of beams, with long-range, low expansion, and self-healing properties, opening new possibilities for non-astrophysical applications.

Mercury, the innermost planet, formed under highly reduced conditions, based mainly on surface Fe, S, and Si abundances determined from MESSENGER mission data. The minor element Cr may serve as an independent oxybarometer, but only very limited Cr data have been previously reported for Mercury. We report Cr/Si abundances across Mercury's surface based on MESSENGER X-Ray Spectrometer data throughout the spacecraft's orbital mission. The heterogeneous Cr/Si ratio ranges from 0.0015 in the Caloris Basin to 0.0054 within the high-magnesium region, with an average southern hemisphere value of 0.0008 (corresponding to about 200 ppm Cr). Absolute Cr/Si values have systematic uncertainty of at least 30%, but relative variations are more robust. By combining experimental Cr partitioning data along with planetary differentiation modeling, we find that if Mercury formed with bulk chondritic Cr/Al, Cr must be present in the planet's core and differentiation must have occurred at log fO2 in the range of IW-6.5 to IW-2.5 in the absence of sulfides in its interior, and a range of IW-5.5 to IW-2 with an FeS layer at the core-mantle boundary. Models with large fractions of Mg-Ca-rich sulfides in Mercury's interior are more compatible with moderately reducing conditions (IW-5.5 to IW-4) owing to the instability of Mg-Ca-rich sulfides at elevated fO2. These results indicate that if Mercury differentiated at a log fO2 lower than IW-5.5, the presence of sulfides whether in the form of a FeS layer at the top of the core or Mg-Ca-rich sulfides within the mantle would be unlikely.

Mercury's orbit can destabilize, resulting in a collision with either Venus or the Sun. Chaotic evolution can cause $g_1$ to decrease to the approximately constant value of $g_5$ and create a resonance. Previous work has approximated the variation in $g_1$ as stochastic diffusion, which leads to a model that can reproduce the Mercury instability statistics of secular and $N$-body models on timescales longer than 10~Gyr. Here we show that the diffusive model underpredicts the Mercury instability probability by a factor of 3-10,000 on timescales less than 5~Gyr, the remaining lifespan of the Solar System. This is because $g_1$ exhibits larger variations on short timescales than the diffusive model would suggest. To better model the variations on short timescales, we build a new subdiffusive model for $g_1$ including a quadratic spring potential above a certain value of $g_1$, which we refer to as a soft upper boundary. Subdiffusion is similar to diffusion, but exhibits larger displacements on short timescales and smaller displacements on long timescales. We choose model parameters based on the short-time behavior of the $g_1$ trajectories in the $N$-body simulations, leading to a tuned model that can reproduce Mercury instability statistics from 1-40~Gyr. This work motivates several questions in planetary dynamics: Why does subdiffusion better approximate the variation in $g_1$ than standard diffusion? Why is a soft upper boundary condition on $g_1$ an appropriate approximation? Why is there an upper bound on $g_1$, but not a lower bound that would prevent it from reaching $g_5$?

To determine the nature of the PeVatron's emission (hadronic or leptonic), it is essential to characterize the physical parameters of the environment from where it originates. We unambiguously confirm the association of molecular gas with the PeVatron candidate LHAASO J2108+5157 using unprecedented high angular-resolution (17$^{\prime \prime}$) $^{12,13}$CO($J$=1$\rightarrow$0) observations carried out with the Nobeyama 45m radio telescope. We characterize a molecular cloud in the vicinity of the PeVatron candidate LHAASO J2108+5157 by determining its physical parameters from our $^{12,13}$CO($J$=1$\rightarrow$0) line observations. We use an updated estimation of the distance to the cloud, which allows us to obtain a more reliable result. The molecular emission is compared with excess gamma-ray images obtained with Fermi--LAT at energies above 2 GeV to search for spatial correlations and test a possible hadronic ($\pi^0$ decay) origin for the gamma-ray emission. We find that the morphology of the spatial distribution of the CO emission is strikingly similar to that of the Fermi--LAT excess gamma-ray. By combining our observations with archival 21cm HI line data, the nucleons (HI + H$_2$) number density of the target molecular cloud is found to be 133.0 $\pm$ 45.0 cm$^{-3}$, for the measured angular size of 0.55 $\pm$ 0.02$^\circ$ at a distance of 1.6 $\pm$ 0.1 kpc. The resulting total mass of the cloud is M(HI +H$_2$) = 7.5$\pm$2.9$\times$10$^3$ M$_{\odot}$. Under a hadronic scenario, we obtain a total energy of protons of W$_p$ = 4.3$\pm$1.5 $\times$ 10$^{46}$ erg with a cutoff of 700$\pm$300 TeV, which reproduces the sub-PeV gamma-ray emission. We identified a molecular cloud in the vicinity of LHAASO J2107+5157 as the main target where cosmic rays from an unknown PeVatron produce the observed gamma-ray emission via $\pi^0$ decay.

We present optical and infrared imaging and spectroscopy of the R Coronae Borealis-type (R Cor Bor) star IRAS 00450+7401. Optical spectra further confirm its classification as a cool R Cor Bor system, having a hydrogen-deficient carbon star spectral sub-class of HdC5 or later. Mid-infrared spectroscopy reveals the typical ~8 um ``hump'' seen in other R Cor Bor stars and no other features. A modern-epoch spectral energy distribution shows bright emission from hot dust having Tdust>600 K. Historical infrared data reveal generally cooler dust color temperatures combined with long-term fading trends, but provide no discernible correlation between flux level and temperature. Investigating the most mid-infrared variable R Cor Bor stars found in IRAS, AKARI, and WISE data reveals similar fading trends, bursts that can show a factor of up to 10 change in flux density between epochs, and blackbody-fit dust color temperatures that span 400-1300 K. While some R Cor Bor stars such as IRAS 00450+7401 appear to undergo fade/burst cycles in the mid-infrared, significant gaps in temporal coverage prevent conclusively identifying any preferred timescale for their mid-infrared variability and circumstellar dust temperature changes.

Utilising archival Chandra X-ray Observatory data and Hubble Space Telescope globular cluster catalogues, we probe the time-domain properties of the low mass X-ray binary population in the elliptical galaxy NGC 4261. Of the 98 unique X-ray sources identified in this study, 62 sources are within the optical field of view and, of those, 33% are aligned with an optical cluster counterpart. We find twenty X-ray sources coincident with globular clusters; two are previously discovered ultra-luminous X-ray sources (ULXs) and eighteen are low mass X-ray binaries (GCLMXBs) with $L_X < 10^{39}$ erg s$^{-1}$. ULXs are a heterogeneous class of extremely bright X-ray binaries ($L_X > 10^{39}$ erg s$^{-1}$) and ULXs located in globular clusters (GCULXs) and may be indicators of black holes. Identifying these unusually X-ray bright sources and measuring their optical properties can provide valuable constraints on the progenitors of gravitational wave sources. We compare observations of these sources to the twenty previously-studied GCULXs from five other early-type galaxies, and find that GCULXs in NGC 4261 are of similar colour and luminosity and do not significantly deviate from the rest of the sample in terms of distance from the galaxy centre or X-ray luminosity. Both the GCULX and low mass X-ray binary (GCLMXB) populations of NGC 4261 show long term variability; the former may have implications for fast radio bursts originating in globular clusters and the latter will likely introduce additional scatter into the low mass end of GCLMXB X-ray luminosity functions.

The evolved massive star populations of the Local Group galaxies are generally thought to be well-understood. However, recent work suggested that the Wolf-Rayet (WR) content of M31 may have been underestimated. We therefore began a pilot project to search for new WRs in M31 and re-examine the completeness of our previous WR survey finished almost a decade prior. Our improved imaging data and spectroscopic follow-up confirmed 19 new WRs across three small fields in M31. These newly discovered WRs are generally fainter than the previously known sample due to slightly increased reddening as opposed to intrinsic faintness. From these findings, we estimate that there are another ~60 WRs left to be discovered in M31; however, the overall ratio of WN-type (nitrogen-rich) to WC-type (carbon-rich) WRs remains unchanged with our latest additions to the M31 WR census. We are in the process of extending this pilot WR survey to include the rest of M31, and a more complete population will be detailed in our future work.

We present in this Letter the first global comparison between traditional line-tied steady state magnetohydrodynamic models and a new, fully time-dependent thermodynamic magnetohydrodynamic simulation of the global corona. The maps are scaled to the approximate field distributions and magnitudes around solar minimum using the Lockheed Evolving Surface-Flux Assimilation Model to incorporate flux emergence and surface flows over a full solar rotation, and include differential rotation and meridional flows. Each time step evolves the previous state of the plasma with a new magnetic field input boundary condition. We find that this method is a significant improvement over steady-state models, as it closely mimics the constant photospheric driving on the Sun. The magnetic energy levels are higher in the time-dependent model, and coronal holes evolve more along the following edge than they do in steady-state models. Coronal changes, as illustrated with forward-modeled emission maps, evolve on longer timescales with time-dependent driving. We discuss implications for active and quiet Sun scenarios, solar wind formation, and widely-used steady state assumptions like potential field source surface calculations.

We report on a remarkable change in the spectral energy distribution (SED) of Mrk 841, providing new insights on how the soft X-ray excess emission in active galactic nuclei (AGNs) is produced. By Swift monitoring of a sample of Seyfert-1 galaxies, we found an X-ray spectral hardening event in Mrk 841. We thereby triggered our XMM-Newton, NuSTAR, and HST observations in 2022 to study this event. Our previous investigations of such events in other AGNs had shown that they are caused by obscuring winds. However, the event in Mrk 841 has different spectral characteristics and origin. We find it is the soft X-ray excess component that has become dimmer. This is, importantly, accompanied by a similar decline in the optical/UV continuum, suggesting a connection to the soft X-ray excess. In contrast, there is relatively little change in the X-ray power-law and the reflection components. Our SED modeling suggests that the soft X-ray excess is the high-energy extension of the optical/UV disk emission, produced by warm Comptonization. We find the temperature of the disk dropped in 2022, explaining the observed SED dimming. We then examined the Swift data, taken over 15 years, to further decipher the UV and X-ray variabilities of Mrk 841. A significant relation between the variabilities of the X-ray spectral hardness and that of the UV continuum is found, again suggesting that the soft excess and the disk emission are interlinked. This is readily explicable if the soft excess is produced by warm Comptonization.

The analysis of core-collapse supernova (CCSN) environments can provide important information on the life cycle of massive stars and constrain the progenitor properties of these powerful explosions. The MUSE instrument at the VLT enables detailed local environment constraints of the progenitors of large samples of CCSNe. Using a homogeneous SN sample from the ASAS-SN survey has enabled us to perform a minimally biased statistical analysis of CCSN environments. We analyze 111 galaxies observed by MUSE that hosted 112 CCSNe detected or discovered by the ASAS-SN survey between 2014 and 2018. The majority of the galaxies were observed by the the AMUSING survey. Here we analyze the immediate environment around the SN locations and compare the properties between the different CCSN types and their light curves. We used stellar population synthesis and spectral fitting techniques to derive physical parameters for all HII regions detected within each galaxy, including the star formation rate (SFR), H$\alpha$ equivalent width (EW), oxygen abundance, and extinction. We found that stripped-envelope (SE) SNe occur in environments with a higher median SFR, H$\alpha$ EW, and oxygen abundances than SNe II and SNe IIn/Ibn. The distributions of SNe II and IIn are very similar, indicating that these events explode in similar environments. For the SESNe, SNe Ic have higher median SFRs, H$\alpha$ EWs, and oxygen abundances than SNe Ib. SNe IIb have environments with similar SFRs and H$\alpha$ EWs to SNe Ib, and similar oxygen abundances to SNe Ic. We also show that the postmaximum decline rate, $s$, of SNe II correlates with the H$\alpha$ EW, and that the luminosity and the $\Delta m_{15}$ parameter of SESNe correlate with the oxygen abundance, H$\alpha$ EW, and SFR at their environments. This suggests a connection between the explosion mechanisms of these events to their environment properties.

Core-collapse supernovae (CCSNe) are widely accepted to be caused by the explosive death of massive stars with initial masses $\gtrsim 8$M$_\odot$. There is, however, a comparatively poor understanding of how properties of the progenitors -- mass, metallicity, multiplicity, rotation etc. -- manifest in the resultant CCSN population. Here we present a minimally biased sample of nearby CCSNe from the ASAS-SN survey whose host galaxies were observed with integral-field spectroscopy using MUSE at the VLT. This dataset allows us to analyze the explosion sites of CCSNe within the context of global star formation properties across the host galaxies. We show that the CCSN explosion site oxygen abundance distribution is offset to lower values than the overall HII region abundance distribution within the host galaxies. We further show that within the subsample of low-metallicity host galaxies, the CCSNe unbiasedly trace the star-formation with respect to oxygen abundance, while for the sub-sample of higher-metallicity host galaxies, they preferentially occur in lower-abundance star-forming regions. We estimate the occurrence of CCSNe as a function of oxygen abundance per unit star formation, and show that there is a strong decrease as abundance increases. Such a strong and quantified metallicity dependence on CCSN production has not been shown before. Finally, we discuss possible explanations for our result and show that each of these has strong implications for our understanding of CCSNe and massive star evolution.

The identification of molecules in exoplanetary atmospheres is only possible thanks to the availability of high-resolution molecular spectroscopic data. However, due to its intensive and time-consuming generation process, at present, only on order 100 molecules have high-resolution spectroscopic data available, limiting new molecular detections. Using routine quantum chemistry calculations (i.e., scaled harmonic frequency calculations using the B97-1/def2-TZVPD model chemistry with median errors of 10cm-1), here we present a complementary high-throughput approach to rapidly generate approximate vibrational spectral data for 2743 molecules made from the biologically most important elements C, H, N, O, P and S. Though these data are not accurate enough to enable definitive molecular detections and does not seek to replace the need for high-resolution data, it has powerful applications in identifying potential molecular candidates responsible for unknown spectral features. We explore this application for the 4.1 micron (2439cm-1) feature in the atmospheric spectrum of WASP-39b, listing potential alternative molecular species responsible for this spectral line, together with SO2. Further applications of this big data compilation also include identifying molecules with strong absorption features that are likely detectable at quite low abundances, and training set for machine learning predictions of vibrational frequencies. Characterising exoplanetary atmospheres through molecular spectroscopy is essential to understand the planet's physico-chemical processes and likelihood of hosting life. Our rapidly generated quantum chemistry big data set will play a crucial role in supporting this understanding by giving directions into possible initial identifications of the more unusual molecules to emerge.

We performed $42$ simulations of the first star formation with initial supersonic gas flows relative to the dark matter at the cosmic recombination era. Increasing the initial streaming velocities led to delayed halo formation and increased halo mass, enhancing the mass of the gravitationally shrinking gas cloud. For more massive gas clouds, the rate of temperature drop during contraction, in other words, the structure asymmetry, becomes more significant. When the maximum and minimum gas temperature ratios before and after contraction exceed about ten, the asymmetric structure of the gas cloud prevails, inducing fragmentation into multiple dense gas clouds. We continued our simulations until $10^5$ years after the first dense core formation to examine the final fate of the massive star-forming gas cloud. Among the $42$ models studied, we find the simultaneous formation of up to four dense gas clouds, with a total mass of about $2254\,M_\odot$. While the gas mass in the host halo increases with increasing the initial streaming velocity, the mass of the dense cores does not change significantly. The star formation efficiency decreases by more than one order of magnitude from $\epsilon_{\rm III} \sim 10^{-2}$ to $10^{-4}$ when the initial streaming velocity, normalised by the root mean square value, increases from 0 to 3.

The first series of observations by the James Webb Space Telescope (JWST) discovered unexpectedly abundant luminous galaxies at high redshift, posing possibly a serious challenge to popular galaxy formation models. We study early structure formation in a cosmological model with a blue, tilted power spectrum (BTPS) given by $P(k) \propto k^{m_{\rm s}}$ with $m_{\rm s} > 1$ at small length scales. We run a set of cosmological $N$-body simulations and derive the abundance of dark matter halos and of galaxies under simplified assumptions on star formation efficiency. The enhanced small-scale power allows rapid formation of nonlinear structure at $z>7$, and galaxies with stellar mass exceeding $10^{10}\,M_\odot$ can be formed by $z=9$. Because of frequent mergers, the structure of galaxies and galaxy groups appears overall clumpy. The BTPS model reproduces the observed stellar mass density at $z=7-9$, and thus eases the claimed tension between galaxy formation theory and recent JWST observations. Large-scale structure of the present-day Universe is largely unaffected by the modification of the small-scale power spectrum. Finally, we discuss the formation of the first stars and early super-massive black holes in the BTPS model.

We collected a total of 4,058 X-ray binary stars, out of which 339 stars had three atmospheric parameters for optical companions from Gaia and LAMOST, while 264 stars had masses and radii of optical companions determined using stellar evolution models. We conducted a thorough discussion on the reliability of each parameter. The statistical analysis revealed a noticeable bimodal distribution in the mass, radius, and age of the optical components. Our findings led to the proposal of a new quantitative classification criterion for X-ray binary stars. In this classification, one type is categorized as high-mass, high-temperature, and young, while the other type is classified as low-mass, low-temperature, and old, corresponding to High-Mass X-Ray Binaries (HMXBs) and Low-Mass X-Ray Binaries (LMXBs), respectively. The dividing lines were established at 11,500 K, 1.7 $M_{\odot}$, and 0.14 Gyr. The classification results of the three parameters showed consistency with one another in 90\% of the cases and were in agreement with previous classifications in 80\% of the cases. We found that the parameters of the two types had well-defined boundaries and distinct patterns. Based on our findings, we suggest that temperature is the best parameter for classification. Therefore, we propose that X-ray binary stars should be classified into high-temperature X-ray binaries (HTXBs) and low-temperature X-ray binaries (LTXBs). We believe that this classification is more convenient in practice and aligns well with physics.

Ground-based solar image restoration is a computationally expensive procedure that involves nonlinear optimization techniques. The presence of atmospheric turbulence produces perturbations in individual images that make it necessary to apply blind deconvolution techniques. These techniques rely on the observation of many short exposure frames that are used to simultaneously infer the instantaneous state of the atmosphere and the unperturbed object. We have recently explored the use of machine learning to accelerate this process, with promising results. We build upon this previous work to propose several interesting improvements that lead to better models. As well, we propose a new method to accelerate the restoration based on algorithm unrolling. In this method, the image restoration problem is solved with a gradient descent method that is unrolled and accelerated aided by a few small neural networks. The role of the neural networks is to correct the estimation of the solution at each iterative step. The model is trained to perform the optimization in a small fixed number of steps with a curated dataset. Our findings demonstrate that both methods significantly reduce the restoration time compared to the standard optimization procedure. Furthermore, we showcase that these models can be trained in an unsupervised manner using observed images from three different instruments. Remarkably, they also exhibit robust generalization capabilities when applied to new datasets. To foster further research and collaboration, we openly provide the trained models, along with the corresponding training and evaluation code, as well as the training dataset, to the scientific community.

A recent ultraviolet luminosity function (UVLF) analysis in the Hubble Frontier Fields, behind foreground lensing clusters, has helped solidify estimates of the faint-end of the $z \sim 5-9$ UVLF at up to five magnitudes fainter than in the field. These measurements provide valuable information regarding the role of low luminosity galaxies in reionizing the universe and can help in calibrating expectations for JWST observations. We fit a semi-empirical model to the lensed and previous UVLF data from Hubble. This fit constrains the average star formation efficiency (SFE) during reionization, with the lensed UVLF measurements probing halo mass scales as small as $M \sim 2 \times 10^9 {\rm M}_\odot$. The implied trend of SFE with halo mass is broadly consistent with an extrapolation from previous inferences at $M \gtrsim 10^{10} {\rm M}_\odot$, although the joint data prefer a shallower SFE. This preference, however, is partly subject to systematic uncertainties in the lensed measurements. Near $z \sim 6$ we find that the SFE peaks at $\sim 20 \%$ between $\sim 10^{11}-10^{12} {\rm M}_\odot$. Our best fit model is consistent with Planck 2018 determinations of the electron scattering optical depth, and most current reionization history measurements, provided the escape fraction of ionizing photons is $f_{\rm esc} \sim 10-20\%$. The joint UVLF accounts for nearly $80\%$ of the ionizing photon budget at $z \sim 8$. Finally, we show that recent JWST UVLF estimates at $z \gtrsim 11$ require strong departures from the redshift evolution suggested by the Hubble data.

We present an extended grid of multi-epoch 1D nonlocal thermodynamic equilibrium radiative transfer calculations for nebular-phase Type Ibc supernovae (SNe) from He-star explosions. Compared to Dessart+21, we study the spectral evolution from 100 to about 450d and augment the model set with progenitors that were evolved without wind mass loss. Models with the same final, preSN mass have similar yields and produce essentially the same emergent spectra. Hence, the uncertain progenitor mass loss history compromises the inference of the initial, main sequence mass. This shortcoming does not affect Type IIb SNe. However, our 1D models with a different preSN mass tend to yield widely different spectra, as seen through variations in the strong emission lines due to [NII]6548-6583, [OI]6300-6364, [CaII]7291-7323, [NiII]7378, and the forest of FeII lines below 5500A. At the lower mass end, the ejecta are He rich and at 100d cool through HeI, NII, CaII, and FeII lines, with NII and FeII dominating at 450d. These models, associated with He giants, conflict with observed SNe Ib, which typically lack strong NII emission. Instead they may lead to SNe Ibn or, because of additional stripping by a companion star, ultra-stripped SNe Ic. In contrast, for higher preSN masses, the ejecta are progressively He poor and cool at 100d through OI, CaII, and FeII lines, with OI and CaII dominating at 450d. Nonuniform, aspherical, large-scale mixing rather than composition differences likely determines the SN type at intermediate preSN masses. Variations in clumping, mixing, as well as departures from spherical symmetry would increase the spectral diversity but also introduce additional degeneracies. More robust predictions from spectral modeling require a careful attention to the initial conditions informed by physically-consistent 3D explosion models [abridged].

The homochirality of biological molecules on the Earth is a long-standing mystery regarding the origin of life. Circularly polarized ultraviolet (UV) light could induce the enantiomeric excess of biological molecules in the interstellar medium, leading to the homochirality on the earth. By performing 3D radiation transfer simulations with multiple scattering processes in interstellar dusty slabs, we study the generation of circular polarization (CP) of ultraviolet light at Lyman $\alpha$ ($\lambda = 0.1216~{\rm \mu m}$) as well as in the near-infrared (NIR, $\lambda = 2.14~{\rm \mu m}$) wavelengths. Our simulations show that the distributions of CP exhibit a symmetric quadrupole pattern, regardless of wavelength and viewing angle. The CP degree of scattered light from a dusty slab composed of aligned grains is $\sim 15$ percent for Ly$\alpha$ and $\sim 3$ percent at NIR wavelengths in the case of oblate grains with an MRN size distribution. We find that the CP degree of Ly$\alpha$ is well correlated with that in the NIR regardless of viewing angles, whilst being a factor of $\sim 5$ higher. Thus, high CP of Ly$\alpha$ is expected in sites where NIR CP is detected. We suggest that such circularly polarized Ly$\alpha$ may initiate the enantiomeric excess of biological molecules in space.

Four pulsars, PSRs J1838+1523, J1901+0510, J1909+0007 and J1929+1844, are found to exhibit bright and weak emission states from sensitive FAST observations. New FAST observations have measured their polarization properties for the two states, and revealed that the polarization profiles, linear polarization percentage, and polarization position angle curves, as well as circular polarization percentage are partially or entirely different in the two emission states. Remarkably, PSR J1838+1523 has very different slopes for the polarization position angle curves. PSR J1901+0510 has a wider profile and a higher linear polarization in the weak state than those in the bright state. PSR J1909+0007 has very distinct polarization angle curves for the two modes. While in the case of PSR J1929+1844, the central profile component evolves with frequency in the bright state, and the senses of circular polarization are opposite in the two modes. The different polarization properties of the two emission states provide valuable insights into the physical processes and emission conditions in pulsar magnetosphere.

Determining the behaviour of convection and clouds is one of the biggest challenges in our understanding of exoplanetary climates. Given the lack of in situ observations, one of the most preferable approaches is to use cloud-resolving or cloud-permitting models (CPM). Here we present CPM simulations in a quasi-global domain with high spatial resolution (4$\times$4 km grid) and explicit convection to study the cloud regime of 1 to 1 tidally locked rocky planets orbiting around low-mass stars. We show that the substellar region is covered by deep convective clouds and cloud albedo increases with increasing stellar flux. The CPM produces relatively less cloud liquid water concentration, smaller cloud coverage, lower cloud albedo, and deeper H2O spectral features than previous general circulation model (GCM) simulations employing empirical convection and cloud parameterizations. Furthermore, cloud streets--long bands of low-level clouds oriented nearly parallel to the direction of the mean boundary-layer winds--appear in the CPM and substantially affect energy balance and surface precipitation at a local level.

We performed photometric observations of the S255, S257, S140, NGC7358 and the Orion~Bar photo-dissociation regions (PDRs) at 2 micron using narrow-band filters centered on the Br-gamma, H2 and [FeII] lines, as well as the narrow-band Kcont and the broad-band H filters for continuum subtraction. The observations were done with the 2.5-m telescope of the SAI Caucasian Mountain Observatory and the near-infrared camera and spectrograph ASTRONIRCAM. We find several high-density arc-like structures in the Br-gamma and [FeII] images of the ionized gas in NGC7538 and extended shells and arcs visible through the H2 emission. The H ionization front and H2 dissociation front are merged in NGC7538. In S255 and S257 we detected only Br-gamma emission from the HII regions and bright H2 emission from the PDRs. The projected distance between the H ionization and H2 dissociation fronts are approx. 0.3-0.4 pc, which cannot be explained using models of a uniform medium. Most probably, the ionized and neutral gas in these PDRs is clumpy. The H-to-H2 transitions in the NGC7538, S255, S257 and S140 PDRs are gradual with no sharp borders. This conclusion also confirms the suggestion of a clumpy medium.

Exoplanet detection by direct imaging is a difficult task: the faint signals from the objects of interest are buried under a spatially structured nuisance component induced by the host star. The exoplanet signals can only be identified when combining several observations with dedicated detection algorithms. In contrast to most of existing methods, we propose to learn a model of the spatial, temporal and spectral characteristics of the nuisance, directly from the observations. In a pre-processing step, a statistical model of their correlations is built locally, and the data are centered and whitened to improve both their stationarity and signal-to-noise ratio (SNR). A convolutional neural network (CNN) is then trained in a supervised fashion to detect the residual signature of synthetic sources in the pre-processed images. Our method leads to a better trade-off between precision and recall than standard approaches in the field. It also outperforms a state-of-the-art algorithm based solely on a statistical framework. Besides, the exploitation of the spectral diversity improves the performance compared to a similar model built solely from spatio-temporal data.

The circumgalactic medium (CGM) is the most massive baryonic component of a spiral galaxy, shock heated to about $10^6$K for an $\rm L^{\star}$ galaxy. The CGM of the Milky Way has been well-characterized through X-ray absorption line spectroscopy. However, the paucity of bright background sources makes it challenging to probe the CGM of external galaxies. Previously, using broad OVI absorption as a signpost, we successfully detected the CGM of one galaxy in X-rays. Here we report on the detection of the OVII $K\alpha$ absorption line at the redshift of a spiral galaxy at $z\approx0.225$ using 1.2 Ms of Chandra observations. This is a robust detection, clearly showing the presence of the hot gas. The mass in the hot phase is at least an order of magnitude larger than that in the cooler phases detected in the UV. The presence of hot gas $116h^{-1}$kpc from the center of this galaxy provides credence to the existence of the extended CGM of the Milky Way. There has been a report of the detection of OVII absorption from the warm-hot intergalactic medium in this sightline using stacking analysis on an older dataset. We argue that the absorption line is from the CGM of the $z\approx0.225$ galaxy instead.

A negative correlation was found to exist between the low-energy spectral index and the redshift of gamma-ray bursts (GRBs) by Amati(2002). It was later confirmed by Geng(2013) and Gruber(2014), but the correlation was also found to be quite dispersive when the sample size was significantly expanded. In this study, we have established two even larger samples of gamma-ray bursts to further examine the correlation. One of our sample is consisted of 316 GRBs detected by the Swift satellite, and the other one is consisted of 80 GRBs detected by the Fermi satellite. It is found that there is no correlation between the two parameters for the Swift sample, but there does exist a weak negative correlation for the Fermi sample. The correlation becomes even more significant when the spectral index at the peak flux is considered. It is argued that the absence of the correlation in the Swift sample may be due to the fact that Swift has a very narrow energy response so that it could not measure the low-energy spectral index accurately enough. Further studies based on even larger GRB samples are solicited.

The Tayler-Spruit dynamo is one of the most promising mechanisms proposed to explain angular momentum transport during stellar evolution. Its development in proto-neutron stars spun-up by supernova fallback has also been put forward as a scenario to explain the formation of very magnetized neutron stars called magnetars. Using three-dimensional direct numerical simulations, we model the proto-neutron star interior as a stably stratified spherical Couette flow with the outer sphere that rotates faster than the inner one. We report the existence of two subcritical dynamo branches driven by the Tayler instability. They differ by their equatorial symmetry (dipolar or hemispherical) and the magnetic field scaling, which is in agreement with different theoretical predictions (by Fuller and Spruit, respectively). The magnetic dipole of the dipolar branch is found to reach intensities compatible with observational constraints on magnetars.

Once further confirmed in future analyses, the radius and mass measurement of HESS J1731-347 with $M=0.77^{+0.20}_{-0.17}~M_{\odot}$ and $R=10.4^{+0.86}_{-0.78}~\rm km$ will be among the lightest and smallest compact object ever detected. This raises many questions about its nature and opens up the window for different theories to explain such a measurement. In this article, we use the information from Doroshenko et al. (2022) on the mass, radius, and surface temperature together with the multi-messenger observations of neutron stars to investigate the possibility that HESS J1731-347 is one of the lightest observed neutron star (NS), a strange star (SS), a hybrid star (HS) with an early deconfinement phase transition, or a dark matter (DM) admixed neutron star. The nucleonic and quark matter are modeled within realistic equation of states (EoSs) with a self-consistent calculation of the pairing gaps in quark matter. By performing the joint analysis of the thermal evolution and mass-radius constraint, we find evidence that within a 1$\sigma$ confidence level, HESS J1731-347 is consistent with the neutron star scenario with the soft EoS as well as with a strange and hybrid star with the early deconfinement phase transition with a strong quark pairing and neutron star admixed with dark matter.

The bright disk of RX J1604 has a very simple axisymmetric structure and is well suited as a benchmark object for accurate photo-polarimetric measurements. We used archival data of RX J1604 from the ESO archive and carefully corrected the polarization signal for instrumental effects, also taking the interstellar polarization into account. We derive accurate radial disk profiles for the intrinsic polarized intensity, ${\hat{Q}}_{\varphi}(r)/I_{\star}$, and measure different profile peak radii for different bands because of the wavelength dependence of the dust opacity. The disk-integrated polarization is $\hat{Q}_{\varphi}/I_{\star} = 0.92 \pm 0.04\%$ for the R band and $1.51 \pm 0.11\%$ for the J band, indicating a red color for the polarized reflectivity of the disk. The intensity of the disk is $I_{\rm disk}/I_{\star} = 3.9 \pm 0.5 \%$ in the J band, and the fractional polarization is $\hat{p}_{\varphi} = 38 \pm 4\%$ for the J band and $42 \pm 2\%$ for the H band. The comparison with the IR excess for RX J1604 yields an apparent disk albedo of about $\Lambda_{I} \approx 0.16 \pm 0.08$. We also find that previously described shadows seen in the R band data are likely affected by calibration errors. Using dust scattering models for transition disks, We derive approximate J band values for the scattering albedo $\omega \approx 0.5$, scattering asymmetry $g \approx 0.5$, and scattering polarization $p_{\rm max} \approx 0.7$ for the dust. The positive R to J band color for the polarized reflectivity is mainly a result of the wavelength dependence of dust parameters because the scattering geometry is expected to be very similar for different colors. This work demonstrates the potential of accurate photo-polarimetric measurements of the circumstellar disk RX J1604 for the determination of dust scattering parameters that strongly constrain the physical properties of the dust.

We present a detailed study of the generation of large primordial non-Gaussianities during the slow-roll (SR) to ultra-slow roll (USR) transitions in the framework of Galileon inflation. We found out that due to having sharp transitions in the USR phase, which persist with a duration of $\Delta {\it N}_{\rm USR} \sim 2$ e-folds, we are able to generate the non-Gaussianity amplitude of the order: $|f_{\rm NL}| \sim {\it O}(10^{-2})$ in the SRI, $-5 < f_{\rm NL} < 5$ in the USR, and $-2 < f_{\rm NL} < 2$ in the SRII phases. As a result, we are able to achieve a cumulative average value of $|f_{\rm NL}| \sim {\it O}(1)$. This implies that our results strictly satisfy Maldacena's no-go theorem in the squeezed limit only for SRI, while they strictly violate the same condition in both the USR and SRII phases. The non-renormalization theorem in the Galileon theory helps to support our results regarding the generation of large mass primordial black holes along with large non-Gaussianities, which we show to be dependent on the specific positions of the transition wave numbers fixed at low scales.

This commentary addresses the anomalies in the results reported from the CRESU-SIS experiment at the Institute of Physics of Rennes, France. This experimental setup is dedicated to studying ion-molecule kinetic in the gas phase at very low temperatures using a uniform supersonic flow reactor. A reinterpretation of the latest study performed with this instrument highly suggests a dramatic decrease in flow density upon the injection of neutral reactants. In particular, these concerns can be related to the diffusion effects prevalent in the reported results on the vast majority of the kinetics experiments conducted with a uniform supersonic flow reactor. The scientific community in the field of low-temperature chemical kinetics in uniform supersonic flow would greatly benefit from being aware of and comprehending these highlighted anomalies because the evidence in this commentary calls into question many of the results published to date.

Understanding the discrepancy between the radii of observed hot Jupiters and standard 'radiative-convective' models remains a hotly debated topic in the exoplanet community. One mechanism which has been proposed to bridge this gap, and which has recently come under scrutiny, is the vertical advection of potential temperature from the irradiated outer atmosphere deep into the interior, heating the deep, unirradiated, atmosphere, warming the internal adiabat, and resulting in radius inflation. Specifically, a recent study which explored the atmosphere of WASP-76b using a 3D, non-grey, GCM suggested that their models lacked radius inflation, and hence any vertical enthalpy advection. Here we perform additional analysis of these, and related, models, focusing on an explicit analysis of vertical enthalpy transport and the resulting heating of the deep atmosphere compared with 1D models. Our results indicate that, after any evolution linked with initialisation, all the WASP-76b models considered here exhibit significant vertical enthalpy transport, heating the deep atmosphere significantly when compared with standard 1D models. Furthermore, comparison of a long time-scale (and hence near steady-state) model with a Jupiter-like internal-structure model suggests not only strong radius-inflation, but also that the model radius, $1.98 \mathrm{R_{J}}$, may be comparable with observations ($1.83\pm0.06 \mathrm{R_{J}}$). We thus conclude that the vertical advection of potential temperature alone is enough to explain the radius inflation of WASP-76b, and potentially other irradiated gas giants, albeit with the proviso that the exact strength of the vertical advection remains sensitive to model parameters, such as the inclusion of deep atmospheric drag.

The OH Megamaser emission in the merging galaxy Arp220 has been re-observed with the Multi-Element Radio Linked Interferometer Network (MERLIN) and the European VLBI Network (EVN). Imaging results of the OH line emission at the two nuclei are found to be consistent with earlier observations and confirm additional extended emission structures surrounding the nuclei. Detailed information about the distributed emission components around the two nuclei has been obtained using a concatenated MERLIN and EVN database with intermediate (40 mas) spatial resolution. Continuum imaging shows a relatively compact West nucleus and a more extended East nucleus in addition to an extended continuum ridge stretching below and beyond the two nuclei. Spectral line imaging show extended emission regions at both nuclei together with compact components and additional weaker components north and south of the West nucleus. Spectral line analysis indicates that the dominant OH line emission originates in foreground molecular material that is part of a large-scale molecular structure that engulfs the whole nuclear region. Compact OH components are representative of star formation regions within the two nearly edge-on nuclei and define the systemic velocities of East and West as 5425 km/s and 5360 km/s. The foreground material at East and West has a 100 km/s lower velocity at 5314 and 5254 km/s. These emission results confirm a maser amplification scenario where the background continuum and the line emission of the star formation regions are amplified by foreground masering material that is excited by the FIR radiation field originating in the two nuclear regions.

How many low-mass stars, brown dwarfs and free-floating planets are in the Milky Way? And how are they distributed in our Galaxy? Recent studies of Milky Way interlopers in high-redshift observations have revealed a 150-300 pc thick disk of these cool stars with 7% of the M-dwarfs in an oblate stellar halo. One can use the High Latitude Survey with the Roman Space Telescope to search for Galactic ultracool dwarfs (spectral classes M, L, T, and Y) to accurately model the 3D structure and the temperature and chemical evolution of the Milky Way disk in these low-mass (sub)stellar objects. Accurate typing has been shown to work on HST grism and photometric data using machine learning techniques. Such an approach can also be applied to Roman photometry, producing accurate photometric typing to within two subtypes. The High Latitude Survey provides enough statistical power to model the Milky Way structural components (thin and thick disks and halo) for M-, L- and T/Y-dwarfs. This approach has the benefit to allow us to constrain scale-lengths, scale-heights and densities, as well as the relative position of our Sun with respect to the disk of dwarf stars of our Milky Way. The total number of each brown dwarf type can be used to infer both the low-mass end of the Galaxy-wide Initial Mass Function (IMF) for the first time, the formation history of low-mass stellar and substellar objects, and the fraction of low-mass stars in the halo, a statistic that can test cold dark matter structure formation theories.

The future Euclid space satellite mission will offer an invaluable opportunity to constrain modifications to general relativity at cosmic scales. We focus on modified gravity models characterised, at linear scales, by a scale-independent growth of perturbations while featuring different testable types of derivative screening mechanisms at smaller nonlinear scales. We consider 3 specific models, namely Jordan-Brans-Dicke (JBD), the normal branch of Dvali-Gabadadze-Porrati (nDGP) gravity and $k$-mouflage (KM) gravity. We provide forecasts from spectroscopic and photometric primary probes by Euclid on the cosmological parameters and the extra parameters of the models, respectively, $\omega_{\rm BD}$, $\Omega_{\rm rc}$ and $\epsilon_{2,0}$, which quantify the deviations from general relativity. This analysis will improve our knowledge of the cosmology of these modified gravity models. The forecasts analysis employs the Fisher matrix method applied to weak lensing (WL); photometric galaxy clustering (GC$_{ph}$); spectroscopic galaxy clustering (GC$_{sp}$) and the cross-correlation (XC) between GC$_{ph}$ and WL. For the Euclid survey specifications we define three scenarios, characterised by different cuts in $\ell$ and $k$, to assess the constraining power of nonlinear scales. For each model we consider two fiducial values for the corresponding model parameter. In an optimistic setting at 68.3\% confidence interval, with Euclid alone we find the following percentage relative errors: for $\log_{10}{\omega_{\rm BD}}$, with a fiducial value of $\omega_{\rm BD}=800$, 35% using GC$_{sp}$ alone, 3.6% using GC$_{ph}$+WL+XC and 3.3% using GC$_{ph}$+WL+XC+GC$_{sp}$; for $\log_{10}{\Omega_{\rm rc}}$, with a fiducial value of $\Omega_{\rm rc}=0.25$, we find respectively 90%, 20% and 17%; finally, for $\epsilon_{2,0}=-0.04$ respectively 5%, 0.15% and 0.14%. (abridged)

Supernova (SN) H0pe was discovered as a new transient in James Webb Space Telescope (JWST) NIRCam images of the galaxy cluster PLCK G165.7+67.0 taken as part of the "Prime Extragalactic Areas for Reionization and Lensing Science" (PEARLS) JWST GTO program (# 1176) on 2023 March 30 (AstroNote 2023-96; Frye et al. 2023). The transient is a compact source associated with a background galaxy that is stretched and triply-imaged by the cluster's strong gravitational lensing. This paper reports spectra in the 950-1370 nm observer frame of two of the galaxy's images obtained with Large Binocular Telescope (LBT) Utility Camera in the Infrared (LUCI) in longslit mode two weeks after the \JWST\ observations. The individual average spectra show the [OII] doublet and the Balmer and 4000 Angstrom breaks at redshift z=1.783+/-0.002. The CIGALE best-fit model of the spectral energy distribution indicates that SN H0pe's host galaxy is massive (Mstar~6x10^10 Msun after correcting for a magnification factor ~7) with a predominant intermediate age (~2 Gyr) stellar population, moderate extinction, and a magnification-corrected star formation rate ~13 Msun/yr, consistent with being below the main sequence of star formation. These properties suggest that H0pe might be a type Ia SN. Additional observations of SN H0pe and its host recently carried out with JWST (JWST-DD-4446; PI: B. Frye) will be able to both determine the SN classification and confirm its association with the galaxy analyzed in this work.

The Milky Way stellar disk has both a thin and a thick component. The thin disk is composed mostly of younger stars ($\lesssim$8 Gyr) with a lower abundance of $\alpha$ elements, while the thick disk contains predominantly older stars ($\gtrsim$8--12 Gyr) with a higher $\alpha$ abundance, giving rise to an $\alpha$-bimodality most prominent at intermediate metallicities. A proposed explanation for the bimodality is an episode of clumpy star formation, where high-$\alpha$ stars form in massive clumps that appear in the first few Gyrs of the Milky Way's evolution, while low-$\alpha$ stars form throughout the disk and over a longer time span. To better understand the evolution of clumps, we track them and their constituent stars in two clumpy Milky Way simulations that reproduce the $\alpha$-abundance bimodality, one with 10% and the other with 20% supernova feedback efficiency. We investigate the paths that these clumps take in the chemical space ([O/Fe]--[Fe/H]) as well as their mass, star formation rate (SFR), formation location, lifetime, and merger history. The clumps in the simulation with lower feedback last longer on average, with several lasting hundreds of Myr. Some of the clumps do not reach high-$\alpha$, but the ones that do on average had a higher SFR, longer lifetime, greater mass, and form closer to the galactic center than the ones that do not. Most clumps that reach high-$\alpha$ merge with others and eventually spiral into the galactic center, but shed stars along the way to form most of the thick disk component.

We report the serendipitous discovery of Gaia17bpp/2MASS J19372316+1759029, a star with a deep single large-amplitude dimming event of $\sim$4.5 magnitudes that lasted over 6.5 years. Using the optical to IR spectral energy distribution (SED), we constrain the primary star to be a cool giant M0-III star with effective temperature $T_{\text{eff}}$=3,850 K and radius R=58 R$_{\odot}$. Based on the SED fitting, we obtained a bimodal posterior distribution of primary stellar masses at 1.5 M${\odot}$ and 3.7 M${\odot}$. Within the last 66 years of photometric coverage, no other significant dimming events of this depth and duration were identified in the optical light curves. Using a Gaussian Process, we fit a high-order Gaussian model to the optical and IR light curves and conclude the dimming event exhibits moderate asymmetries from optical to IR. At the minimum of the dimming event, the (W$_{1}$-W$_{2}$) color was bluer by $\sim$0.2 mag relative to the primary star outside the dimming event. The ingress and egress colors show a shallow reddening profile. We suggest that the main culprit of the dimming event is likely due to the presence of a large, optically thick disk transiting the primary giant star. By fitting a monochromatic transit model of an oblate disk transiting a star, we found good agreement with a slow-moving, 0.005 km sec$^{-1}$, disk with a $\sim$1.4 AU radius. We propose that Gaia17bpp belongs to a rare binary star population similar to the Epsilon Aurigae system, which consists of a secondary star enshrouded by an optically thick debris disk.

We introduce a numerical method specifically designed for investigating generic multifield models of inflation where a number of scalar fields $\phi^K$ are minimally coupled to gravity and live in a field space with a non-trivial metric $G_{IJ}(\phi^K)$. Our algorithm consists of three main parts. Firstly, we solve the field equations through the entire inflationary period, deriving predictions for observable quantities such as the spectrum of scalar perturbations, primordial gravitational waves, and isocurvature modes. We also incorporate the transfer matrix formalism to track the behavior of adiabatic and isocurvature modes on super-horizon scales and the transfer of entropy to scalar modes after the horizon crossing. Secondly, we interface our algorithm with Boltzmann integrator codes to compute the subsequent full cosmology, including the cosmic microwave background anisotropies and polarization angular power spectra. Finally, we develop a novel sampling algorithm able to efficiently explore a large volume of the parameter space and identify a sub-region where theoretical predictions agree with observations. In this way, sampling over the initial conditions of the fields and the free parameters of the models, we enable Monte Carlo analysis of multifield scenarios. We test all the features of our approach by analyzing a specific model and deriving constraints on its free parameters. Our methodology provides a robust framework for studying multifield inflation, opening new avenues for future research in the field.

We present an approximate treatment of the mixing between $\nu_{\tau}$ ($\bar\nu_{\tau}$) and a sterile species $\nu_s$ ($\bar\nu_s$) with a vacuum mass-squared difference of $\sim$ 10$^2$-10$^3$ keV$^2$ in protoneutron stars created in core-collapse supernovae. Including production of sterile neutrinos through both resonant flavor conversion and collisions, we track the evolution of the $\nu_{\tau}$ lepton number due to both escape of sterile neutrinos and diffusion. Our approach provides a reasonable treatment of the pertinent processes discussed in previous studies and serves a pedagogical purpose to elucidate the relevant physics. We also discuss refinements needed to study more accurately how flavor mixing with sterile neutrinos affects protoneutron star evolution.

A test fluid composed of relativistic collisionless neutral particles in the background of Kerr metric is expected to generate non-isotropic equilibrium configurations in which the corresponding stress-energy tensor exhibits pressure and temperature anisotropies. This arises as a consequence of the constraints placed on single-particle dynamics by Killing tensor symmetries, leading to a peculiar non-Maxwellian functional form of the kinetic distribution function describing the continuum system. Based on this outcome, in this paper the generation of Kerr-like metric by collisionless N-body systems of neutral matter orbiting in the field of a rotating black hole is reported. The result is obtained in the framework of covariant kinetic theory by solving the Einstein equations in terms of an analytical perturbative treatment whereby the gravitational field is decomposed as a prescribed background metric tensor described by the Kerr solution plus a self-field correction. The latter one is generated by the uncharged fluid at equilibrium and satisfies the linearized Einstein equations having the non-isotropic stress-energy tensor as source term. It is shown that the resulting self-metric is again of Kerr type, providing a mechanism of magnification of the background metric tensor and its qualitative features.

Non-ideal fluids are generally subject to the occurrence of non-isotropic pressure tensors, whose determination is fundamental in order to characterize their dynamical and thermodynamical properties. This requires the implementation of theoretical frameworks provided by appropriate microscopic and statistical kinetic approaches in terms of which continuum fluid fields are obtained. In this paper the case of non-relativistic magnetized fluids forming equilibrium toroidal structures in external gravitational fields is considered. Analytical solutions for the kinetic distribution function are explicitly constructed, to be represented by a Chapman-Enskog expansion around a Maxwellian equilibrium. In this way, different physical mechanisms responsible for the generation of non-isotropic pressures are identified and proved to be associated with the kinetic constraints imposed on single and collective particle dynamics by phase-space symmetries and magnetic field. As a major outcome, the validity of a polytropic representation for the kinetic pressure tensors corresponding to each source of anisotropy is established, whereby directional pressures exhibit a specific power-law functional dependence on fluid density. The astrophysical relevance of the solution for the understanding of fluid plasma properties in accretion-disc environments is discussed.

Astrophysical plasmas in the surrounding of compact objects and subject to intense gravitational and electromagnetic fields are believed to give rise to relativistic regimes. Theoretical and observational evidence suggest that magnetized plasmas of this type are collisionless and can persist for long times (e.g., with respect to a distant observer, coordinate, time), while exhibiting geometrical structures characterized by the absence of well-defined spatial symmetries. In this paper the problem is posed whether such configurations can correspond to some kind of kinetic equilibrium. The issue is addressed from a theoretical perspective in the framework of a covariant Vlasov statistical description, which relies on the method of invariants. For this purpose, a systematic covariant variational formulation of gyrokinetic theory is developed, which holds without requiring any symmetry condition on the background fields. As a result, an asymptotic representation of the relativistic particle magnetic moment is obtained from its formal exact solution, in terms of a suitably-defined invariant series expansion parameter (perturbative representation). On such a basis it is shown that spatially non-symmetric kinetic equilibria can actually be determined, an example being provided by Gaussian-like distributions. As an application, the physical mechanisms related to the occurrence of a non-vanishing equilibrium fluid 4-flow are investigated.

A review of the original thermodynamic formulation of the Tolman-Ehrenfest effect prescribing the temperature profile of uncharged fluid at thermal equilibrium forming stationary configurations in curved space-time is proposed. A statistical description based on relativistic kinetic theory is implemented. In this context the Tolman-Ehrenfest relation arises in the Schwarzschild space-time for collisionless uncharged particles at Maxwellian kinetic equilibrium. However, the result changes considerably when non-ideal fluids, i.e., non-Maxwellian distributions, are treated, whose statistical temperature becomes non-isotropic and gives rise to a tensor pressure. This is associated with phase-space anisotropies in the distribution function, occurring both for diagonal and non-diagonal metric tensors, exemplified by the Schwarzschild and Kerr metrics respectively. As a consequence, it is shown that for these systems it is not possible to define a Tolman-Ehrenfest relation in terms of an isotropic scalar temperature. Qualitative properties of the novel solution are discussed.

Non-ideal fluids are likely to be affected by the occurrence of pressure anisotropy effects, whose understanding for relativistic systems requires knowledge of the energy-momentum tensor. In this paper the case of magnetized jet plasmas at equilibrium is considered, in which both microscopic velocities of constituent particles as well as the continuum fluid flow are treated as relativistic ones. A theoretical framework based on covariant statistical kinetic approach is implemented, which permits the proper treatment of single-particle and phase-space kinetic constraints and, ultimately, the calculation of the system continuum fluid fields associated with physical observables. A Gaussian-like solution for the kinetic distribution function (KDF) is constructed, in which the physical mechanism responsible for the generation of temperature anisotropy is identified with magnetic moment conservation. A Chapman-Enskog representation of the same KDF is then obtained in terms of expansion around an equilibrium isotropic Juttner distribution. This permits the analytical calculation of the fluid 4-flow and stress-energy tensor and the consequent proof that the corresponding kinetic pressure tensor is non-isotropic. As a notable result, the validity of a polytropic representation for the perturbative non-isotropic pressure contributions is established, whereby directional pressures exhibit specific power-law functional dependences on fluid density.

In this study, we look into binaries undergoing gravitational radiation during a hyperbolic passage. Such hyperbolic events can be a credible source of gravitational waves in future detectors. We systematically calculate fluxes of gravitational radiation from such events in the presence of dark matter, also considering the effects of dynamical friction. We also investigate the binary dynamics through the changes in the orbital parameters by treating the potential due to dark matter spike and the dynamical friction effects as a perturbation term. An insight into the effects of such a medium on the binaries from the corresponding osculating elements opens up avenues to study binary dynamics for such events.

The ability of deep learning (DL) approaches to learn generalised signal and noise models, coupled with their fast inference on GPUs, holds great promise for enhancing gravitational-wave (GW) searches in terms of speed, parameter space coverage, and search sensitivity. However, the opaque nature of DL models severely harms their reliability. In this work, we meticulously develop a DL model stage-wise and work towards improving its robustness and reliability. First, we address the problems in maintaining the purity of training data by deriving a new metric that better reflects the visual strength of the "chirp" signal features in the data. Using a reduced, smooth representation obtained through a variational auto-encoder (VAE), we build a classifier to search for compact binary coalescence (CBC) signals. Our tests on real LIGO data show an impressive performance of the model. However, upon probing the robustness of the model through adversarial attacks, its simple failure modes were identified, underlining how such models can still be highly fragile. As a first step towards bringing robustness, we retrain the model in a novel framework involving a generative adversarial network (GAN). Over the course of training, the model learns to eliminate the primary modes of failure identified by the adversaries. Although absolute robustness is practically impossible to achieve, we demonstrate some fundamental improvements earned through such training, like sparseness and reduced degeneracy in the extracted features at different layers inside the model. Through comparative inference on real LIGO data, we show that the prescribed robustness is achieved at practically zero cost in terms of performance. Through a direct search on ~8.8 days of LIGO data, we recover two significant CBC events from GWTC-2.1, GW190519_153544 and GW190521_074359, and report the search sensitivity.

In this note we comment on a recent attempt by P. Burikham, T. Harko, K. Pimsamarn and S. Shahidi [Phys. Rev. D {\bf 107}, 064008 (2023)] to explain the galactic rotation curves as the result of the motion of time-like test particles in the Weyl geometric theory of gravity developed in [Eur. Phys. J. C {\bf 82}, 23 (2022)]. We show that the static, spherically symmetric solution found by the authors is not a solution of the assumed Weyl geometric theory with vectorial nonmetricity, but of the well-known conformally coupled scalar theory over Weyl integrable spacetime with gradient nonmetricity, instead. Besides, the solution found does not respect the gauge symmetry of the underlying theory.

We present the latest results from the search for light dark matter particle interactions with the DarkSide-50 dual-phase liquid argon time projection chamber. This analysis, based on the ionization signal only, improves the existing limits for spin-independent WIMP-nucleon interactions in the $[1.2, 3.6]$ GeV/c$^2$ mass range. The sensitivity is extended down to 40 MeV/c$^2$ by assuming the Migdal effect, responsible for an additional ionization signal from the recoiling atom. Finally, we set new constraints to interactions of dark matter particles with electrons in the final state, namely WIMPs, galactic axions, dark photons, and sterile neutrinos.

We present a strategy to obtain equations of general relativity for an irrotational dust continuum within a flow-orthogonal foliation of spacetime from the equations of Newtonian gravitation, and vice versa, without employing a weak field expansion or a limiting process on the speed of light. We argue that writing Newton's equations in a Lagrangian frame and relaxing integrability of vector gradients is sufficient to obtain equations that are identical to Einstein's equations in 3+1 form when respecting the Lorentzian signature of the time-parametrization. We discuss implications and provide an outlook on how to extend the obtained correspondence to more general spacetimes.

Theories of gravity with auxiliary fields are of particular interest since they are able to circumvent Lovelock's theorem while avoiding to introduce new degrees of freedom. This type of theories introduces derivatives of the stress-energy tensor in the modified Einstein equation. This peculiar structure of the field equations was shown to lead to spacetime singularities on the surface of stars. Here we focus on yet another problem afflicting gravity theories with auxiliary field. We show that such theories introduce deviations to the Standard Model unless one severely constrains the parameters of the theory, preventing them to produce significant phenomenology at large scales. We first consider the specific case of Palatini $f({\cal R})$ gravity, to clarify the results previously obtained in arXiv:astro-ph/0308111. We show that the matter fields satisfy the Standard Model field equations which reduce to those predicted by General Relativity in the local frame only at tree level, whereas at higher orders in perturbation theory they are affected by corrections that percolate from the gravity sector regardless of the specific $f({\cal R})$ model considered. Finally, we show that this is a more general issue affecting theories with auxiliary fields connected to the same terms responsible for the appearance of surface singularities.

A direct discovery of the cosmic neutrino background would bring to a closure the searches for relic left-over radiation predicted by the Hot Big Bang cosmology. Recently, the KATRIN experiment put a limit on the local relic neutrino overdensity with respect to the cosmological predicted average value at $\eta \lesssim 10^{11}$ [Phys. Rev. Lett. 129, 011806 (2022)]. In this work, we first examine to what extent such values of $\eta$ are conceivable. We show that even under cavalier assumptions, a cosmic origin of $\eta \gtrsim 10^4$ seems out of reach (with the caveat of forming bound objects under a new force,) but find that a hypothetical local source of low-energy neutrinos could achieve $\eta \sim 10^{11}$. Second, when such values are considered, we point out that the experimental signature in KATRIN and other neutrino-capture experiments changes, contrary to what has hitherto been assumed. Our results are model-independent and maximally accommodating as they only assume the Pauli exclusion principle. As intermittent physics target in the quest for C$\nu$B detection, we identify an experimental sensitivity to $\eta \sim 10^4$ for which conceivable sources exist; to resolve the effect of a degenerate Fermi gas for such overdensity an energy resolution of 10 meV is required.