Papers for Friday, Sep 29 2023

The Nainital-Cape Survey was initiated more than two decades ago aiming to search for and study the pulsational variability in two subclasses of chemically peculiar (CP) stars, namely the Ap and Am stars. In this paper, we present the TESS photometry of 4 targets out of the 369 sample stars observed under the survey, which were not studied before using TESS data. Our results suggest that HD34060 is a rotational variable, HD25487 is of eclipsing nature, HD15550 exhibits pulsational variability while HD48953 is a non-variable star. The diverse variability detected in the studied sources places important constraints for the study of the internal structure and evolution of the CP stars in the presence of surface inhomogeneity, magnetic field, rotation and pulsation.

In theoretical models of galaxy evolution, black hole feedback is a necessary ingredient in order to explain the observed exponential decline in number density of massive galaxies. Most contemporary black hole feedback models in cosmological simulations rely on a constant radiative efficiency (usually $\eta \sim 0.1$) at all black hole accretion rates. We present a synthesis model for the spin-dependent radiative efficiencies of three physical accretion rate regimes, i.e. $\eta = \eta(j, \dot{M}_\mathrm{BH})$, for use in large-volume cosmological simulations. The three regimes include: an advection dominated accretion flow ($\dot{M}_\mathrm{BH} < 0.03\,\dot{M}_\mathrm{Edd}$), a quasar-like mode ($0.03 < \dot{M}_\mathrm{BH} / \dot{M}_\mathrm{Edd} < 0.3$), and a slim disc mode ($\dot{M}_\mathrm{BH} > 0.3\,\dot{M}_\mathrm{Edd}$). Additionally, we include a large-scale powerful jet at low accretion rates. The black hole feedback model we present is a kinetic model that prescribes mass loadings but could be used in thermal models directly using the radiative efficiency. We implemented our model into the \texttt{Simba} galaxy evolution model to determine if it is possible to reproduce galaxy populations successfully, and provide a first calibration for further study. Using a $2\times1024^3$ particle cosmological simulation in a $(150\,\mathrm{cMpc})^3$ volume, we found that the model is successful in reproducing the galaxy stellar mass function, black hole mass-stellar mass relationship, and stellar mass-halo mass relationship. Our model shines when we extrapolate to the galaxy group and cluster scale as it impressively predicts the observed baryon fraction within massive groups and low-mass clusters. Moving forward, this model opens new avenues for exploration of the impact of black hole feedback on galactic environments.

Modelling galaxy formation in hydrodynamic simulations has increasingly adopted various radiative transfer methods to account for photoionization feedback from young massive stars. However, the evolution of HII regions around stars begins in dense star-forming clouds and spans large dynamical ranges in both space and time, posing severe challenges for numerical simulations in terms of both spatial and temporal resolution that depends strongly on gas density ($\propto n^{-1}$). In this work, we perform a series of idealized HII region simulations using the moving-mesh radiation-hydrodynamic code Arepo-RT to study the effects of numerical resolution. The simulated results match the analytical solutions and the ionization feedback converges only if the Str\"omgren sphere is resolved by at least $10$--$100$ resolution elements and the size of each time integration step is smaller than $0.1$ times the recombination timescale. Insufficient spatial resolution leads to reduced ionization fraction but enhanced ionized gas mass and momentum feedback from the HII regions, as well as degrading the multi-phase interstellar medium into a diffuse, partially ionized, warm ($\sim8000$\,K) gas. On the other hand, insufficient temporal resolution strongly suppresses the effects of ionizing feedback. This is because longer timesteps are not able to resolve the rapid variation of the thermochemistry properties of the gas cells around massive stars, especially when the photon injection and thermochemistry are performed with different cadences. Finally, we provide novel numerical implementations to overcome the above issues when strict resolution requirements are not achievable in practice.

We present a new scheme for the classification of the in-situ and accreted globular clusters (GCs). The scheme uses total energy $E$ and $z$-component of the orbital angular momentum and is calibrated using [Al/Fe] abundance ratio. We demonstrate that such classification results in the GC populations with distinct spatial, kinematic, and chemical abundance distributions. The in-situ GCs are distributed within the central 10 kpc of the Galaxy in a flattened configuration aligned with the MW disc, while the accreted GCs have a wide distribution of distances and a spatial distribution close to spherical. In-situ and accreted GCs have different $\rm [Fe/H]$ distributions with the well-known bimodality present only in the metallicity distribution of the in-situ GCs. Furthermore, the accreted and in-situ GCs are well separated in the plane of $\rm [Al/Fe]-[Mg/Fe]$ abundance ratios and follow distinct sequences in the age--$\rm [Fe/H]$ plane. The in-situ GCs in our classification show a clear disc spin-up signature -- the increase of median $V_\phi$ at metallicities $\rm [Fe/H]\approx -1.3\div -1$ similar to the spin-up in the in-situ field stars. This signature signals the MW's disc formation, which occurred $\approx 11.7-12.7$ Gyrs ago (or at $z\approx 3.1-5.3$) according to GC ages. In-situ GCs with metallicities of $\rm [Fe/H]\gtrsim -1.3$ were thus born in the Milky Way disc, while lower metallicity in-situ GCs were born during early, turbulent, pre-disc stages of the evolution of the Galaxy and are part of its Aurora stellar component.

Ultra-luminous X-ray sources (ULXs) are sources observed to exceed the Eddington limit of a stellar-mass black hole (BH). A fraction of ULX sources show X-ray pulses which are evidence for accreting neutron stars (NSs). Theoretical studies have suggested that NSs dominate the compact objects of intrinsic ULXs, even though the majority of observed sample is non-pulsating, implying that X-ray pulses from many NS ULXs are unobservable. We use POSYDON to generate and study X-ray binary populations spanning starburst ages 5 to 1000Myr. Following theoretical predictions for the alignment of the NS spin axis with the accretion disc, we estimate the required accreted mass in ULXs so that the alignment suppresses observable X-ray pulses. While the properties of ULXs are sensitive to model assumptions, there are certain trends that the populations follow. Young and old stellar populations are dominated by BH and NS accretors, respectively. The donors go from massive H-rich main-sequence (MS) stars in young populations (<100Myr) to low-mass post-MS H-rich stars in older populations (>100Myr), with stripped He-rich giant stars dominating the populations at around 100Myr. In addition, we find that NS ULXs exhibit stronger geometrical beaming than BH ULXs, leading to an under-representation of NS accretors in observed populations. Coupled with our finding that X-ray pulses are suppressed in at least 60% of the NS ULXs, we suggest that the observed fraction of ULXs with detectable X-ray pulses is very small, in agreement with observations. This study investigates the effects of age on ULXs as well as the effects of different model assumptions on ULX demographics. We show that geometrical beaming and the mass-accretion phase are critical aspects of understanding ULX observations. Our results suggest that even though most ULXs have accreting NSs, those with observable X-ray pulses would be very few.

Post-common-envelope binaries (PCEBs) containing a white dwarf (WD) and a main-sequence (MS) star can constrain the physics of common envelope evolution and calibrate binary evolution models. Most PCEBs studied to date have short orbital periods ($P_{\rm orb} \lesssim 1$ d), implying relatively inefficient harnessing of binaries' orbital energy for envelope expulsion. Here, we present follow-up observations of five binaries from {\it Gaia} DR3 containing solar-type MS stars and probable ultramassive WDs ($M\gtrsim 1.2\,M_{\odot}$) with significantly wider orbits than previously known PCEBs, $P_{\rm orb} = 18-49$ d. The WD masses are much higher than expected for systems formed via stable mass transfer at these periods, and their near-circular orbits suggest partial tidal circularization when the WD progenitors were giants. These properties strongly suggest that the binaries are PCEBs. Forming PCEBs at such wide separations requires highly efficient envelope ejection, and we find that the observed periods can only be explained if a significant fraction of the energy released when the envelope recombines goes into ejecting it. Using 1D stellar evolution calculations, we show that the binding energy of massive AGB star envelopes is formally {\it positive} if recombination energy is included in the energy budget. This suggests that the star's envelope can be efficiently ejected if binary interaction causes it to recombine, and that a wide range of PCEB orbital periods can potentially result from Roche lobe overflow of an AGB star. This evolutionary scenario may also explain the formation of several wide WD+MS binaries discovered via self-lensing.

M82 is an archetypal starburst galaxy in the local Universe. The central burst of star formation, thought to be triggered by M82's interaction with other members in the M81 group, is driving a multiphase galaxy-scale wind away from the plane of the disk that has been studied across the electromagnetic spectrum. Here, we present new velocity-resolved observations of the [CII] 158$\mu$m line in the central disk and the southern outflow of M82 using the upGREAT instrument onboard SOFIA. We also report the first detections of velocity-resolved ($\Delta V = 10$ km s$^{-1}$) [CII] emission in the outflow of M82 at projected distances of $\approx1-2$ kpc south of the galaxy center. We compare the [CII] line profiles to observations of CO and HI and find that likely the majority ($>55$%) of the [CII] emission in the outflow is associated with the neutral atomic medium. We find that the fraction of [CII] actually outflowing from M82 is small compared to the bulk gas outside the midplane (which may be in a halo or tidal streamers), which has important implications for observations of [CII] outflows at higher redshift. Finally, by comparing the observed ratio of the [CII] and CO intensities to models of photodissociation regions, we estimate that the far-ultraviolet (FUV) radiation field in the disk is $\sim10^{3.5}~G_0$, in agreement with previous estimates. In the outflow, however, the FUV radiation field is 2-3 orders of magnitudes lower, which may explain the high fraction of [CII] arising from the neutral medium in the wind.

In this work we extend the non-ideal magnetohydrodynamics (MHD) solver in the moving mesh code AREPO to include the Hall effect. The core of our algorithm is based on an estimation of the magnetic field gradients by a least-square reconstruction on the unstructured mesh, which we also used in the companion paper for Ohmic and ambipolar diffusion. In an extensive study of simulations of a magnetic shock, we show that without additional magnetic diffusion our algorithm for the Hall effect becomes unstable at high resolution. We can however stabilise it by artificially increasing the Ohmic resistivity, $\eta_{\rm OR}$, so that it satisfies the condition $\eta_{\rm OR} \geq \eta_{\rm H} /5$, where $\eta_{\rm H}$ is the Hall diffusion coefficient. Adopting this solution we find second order convergence for the C-shock and are also able to accurately reproduce the dispersion relation of the whistler waves. As a first application of the new scheme, we simulate the collapse of a magnetised cloud with constant Hall parameter $\eta_{\rm H}$ and show that, depending on the sign of $\eta_{\rm H}$, the magnetic braking can either be weakened or strengthened by the Hall effect. The quasi-Lagrangian nature of the moving mesh method used here automatically increases the resolution in the forming core, making it well suited for more realistic studies with non-constant magnetic diffusivities in the future.

It is well known that supermassive black holes (SMBHs) and their host galaxies co-evolve. AGN feedback plays an important role on this symbiosis. To study the effect of the AGN feedback on the host galaxy, a popular method is to study the star-formation rate (SFR) as a function of the X-ray luminosity (L$_X$). However, hydrodynamical simulations suggest that the cumulative impact of AGN feedback on a galaxy is encapsulated in the mass of the SMBH, M$_{BH}$, rather than the L$_X$. In this study, we compare the SFR of AGN and non-AGN galaxies as a function of L$_X$, M$_{BH}$, Eddington ratio (n$_{Edd}$) and specific black hole accretion rate ($\lambda _{sBHAR}$). For that purpose, we use 122 X-ray AGN in the XMM-XXL field and 3371 galaxies from the VIPERS survey and calculate the SFR$_{norm}$ parameter, defined as the ratio of the SFR of AGN to the SFR of non-AGN galaxies with similar stellar mass, M$_*$, and redshift. Our datasets span a redshift range of $\rm 0.5\leq z\leq 1.2$. The results show that the correlation between SFR$_{norm}$ and M$_{BH}$ is stronger compared to that between SFR$_{norm}$ and L$_X$. A weaker correlation is found between SFR$_{norm}$ and $\lambda _{sBHAR}$. No correlation is detected between SFR$_{norm}$ and n$_{Edd}$. These results corroborate the idea that the M$_{BH}$ is a more robust tracer of the cumulative impact of the AGN feedback compared to the instantaneous accretion rate (L$_X$) and, thus, a better predictive parameter of the changes of the SFR of the host galaxy.

Evidence from different probes of the stellar initial mass function (IMF) of massive early-type galaxies (ETGs) has repeatedly converged on IMFs more bottom-heavy than in the Milky Way (MW). This consensus has come under scrutiny due to often contradictory results from different methods on the level of individual galaxies. In particular, a number of strong lensing probes are ostensibly incompatible with a non-MW IMF. Radial gradients of the IMF -- related to gradients of the stellar mass-to-light ratio $\Upsilon$ -- can potentially resolve this issue. We construct Schwarzschild models allowing for $\Upsilon$-gradients in seven massive ETGs with MUSE and SINFONI observations. We find dynamical evidence that $\Upsilon$ increases towards the center for all ETGs. The gradients are confined to sub-kpc scales. Our results suggest that constant-$\Upsilon$ models may overestimate the stellar mass of galaxies by up to a factor 1.5. For all except one galaxy, we find a radius where the total dynamical mass has a minimum. This minimum places the strongest constraints on the IMF outside the center and appears at roughly 1 kpc. We consider the IMF at this radius characteristic for the main body of each ETG. In terms of the IMF mass-normalization $\alpha$ relative to a Kroupa IMF, we find on average a MW-like IMF $<\alpha_{main}> = 1.03 \pm 0.19$. In the centers, we find concentrated regions with increased mass normalizations that are less extreme than previous studies suggested, but still point to a Salpeter-like IMF, $<\alpha_{cen}> = 1.54 \pm 0.15$

We conduct a long-timescale ($5000\,$d) 3-D simulation of a common-envelope event with a $2\,M_{\odot}$ red giant and a $1\,M_{\odot}$ main sequence companion, using the moving-mesh hydrodynamic solver MANGA. Starting with an orbital radius of $52\,R_{\odot}$, our binary shrinks to an orbital radius of $5\,R_{\odot}$ in $200\,$d. We show that over a timescale of about $1500\,$d, the envelope is completely ejected while $80$ per cent is ejected in about $400\,$d. The complete ejection of the envelope is solely powered by the orbital energy of the binary, without the need for late-time reheating from recombination or jets. Motivated by recent theoretical and observational results, we also find that the envelope enters a phase of homologous expansion about $550\,\rm d$ after the start of our simulation. We also run a simplified 1-D model to show that heating from the central binary in the envelope at late times does not influence the ejection. This homologous expansion of the envelope would likely simplify calculations of the observational implications such as light curves.

Black hole (BH) X-ray binaries cycle through different spectral states of accretion over the course of months to years. Although fluctuations in the BH mass accretion rate are generally recognized as the most important component of state transitions, it is becoming increasingly evident that magnetic fields play a similarly important role. In this article, we present the first radiative two-temperature (2T) general relativistic magnetohydrodynamics (GRMHD) simulations in which an accretion disk transitions from a quiescent state at an accretion rate of $\dot{M} \sim 10^{-10} \dot{M}_{\rm Edd}$ to a hard-intermediate state at an accretion rate of $\dot{M} \sim 10^{-2} \dot{M}_{\rm Edd}$. This huge parameter space in mass accretion rate is bridged by artificially rescaling the gas density scale of the simulations. We present two jetted BH models with varying degrees of magnetic flux saturation. We demonstrate that in `Standard and Normal Evolution' models, which are unsaturated with magnetic flux, the hot torus collapses into a thin and cold accretion disk when $\dot{M} \gtrsim 5\times 10^{-3} \dot{M}_{\rm Edd}$. On the other hand, in `Magnetically Arrested Disk' models, which are fully saturated with vertical magnetic flux, the plasma remains mostly hot with substructures that condense into cold clumps of gas when $\dot{M} \gtrsim 1 \times 10^{-2} \dot{M}_{\rm Edd}$. This suggests that the spectral signatures observed during state transitions are closely tied to the level of magnetic flux saturation.

We report the first detection of the X-ray polarization of the bright transient Swift J1727.8-1613 with the Imaging X-ray Polarimetry Explorer. The observation was performed at the beginning of the 2023 discovery outburst, when the source resided in the bright hard state. We find a time- and energy-averaged polarization degree of 4.1%+/-0.2% and a polarization angle of 2.2+/-1.3 degrees (errors at 68% confidence level; this translates to about 20-sigma significance of the polarization detection). This finding suggests that the hot corona emitting the bulk of the detected X-rays is elongated, rather than spherical. The X-ray polarization angle is consistent with that found in sub-mm wavelengths. Since the sub-mm polarization was found to be aligned with the jet direction in other X-ray binaries, this indicates that the corona is elongated orthogonal to the jet.

We present MESA-Web, a cloud resource with an online interface to the Modules for Experiments in Stellar Astrophysics (MESA) software instrument. MESA-Web allows learners to evolve stellar models without the need to download and install MESA. Since being released in 2015, MESA-Web has delivered over 17,000 calculations to over 2,200 unique learners and currently performs about 11 jobs per day. MESA-Web can be used as an educational tool for stars in the classroom or for scientific investigations. We report on new capabilities of MESA-Web introduced since its 2015 release including learner-supplied nuclear reaction rates, custom stopping conditions, and an expanded selection of input parameters. To foster collaboration we have created a Zenodo MESA-Web community hub where instructors can openly share examples of using MESA-Web in the classroom. We discuss two examples in the current community hub. The first example is a lesson module on Red Giant Branch stars that includes a suite of exercises designed to fit a range of learners and a Jupyter workbook for additional analysis. The second example is lesson materials for an upper-level Astronomy majors course in Stars and Radiation that includes an assignment verifying some of the expected trends that are presented in a popular stellar physics textbook.

Recent models for the inner structure of active galactic nuclei (AGN) aim at connecting the outer region of the accretion disk with the broad-line region and dusty torus through a radiatively accelerated, dusty outflow. Such an outflow not only requires the outer disk to be dusty and so predicts disk sizes beyond the self-gravity limit but requires the presence of nuclear dust with favourable properties. Here we investigate a large sample of type 1 AGN with near-infrared (near-IR) cross-dispersed spectroscopy with the aim to constrain the astrochemistry, location and geometry of the nuclear hot dust region. Assuming thermal equilibrium for optically thin dust, we derive the luminosity-based dust radius for different grain properties using our measurement of the temperature. We combine our results with independent dust radius measurements from reverberation mapping and interferometry and show that large dust grains that can provide the necessary opacity for the outflow are ubiquitous in AGN. Using our estimates of the dust covering factor, we investigate the dust geometry using the effects of the accretion disk anisotropy. A flared disk-like structure for the hot dust is favoured. Finally, we discuss the implication of our results for the dust radius-luminosity plane.

Theoretical and observational studies have suggested that ram-pressure stripping by the intracluster medium can be enhanced during cluster interactions, boosting the formation of the "jellyfish" galaxies. In this work, we study the incidence of galaxies undergoing ram-pressure stripping in 52 clusters of different dynamical states. We use optical data from the WINGS/OmegaWINGS surveys and archival X-ray data to characterise the dynamical state of our cluster sample, applying eight different proxies. We then compute the number of ram-pressure stripping candidates relative to the infalling population of blue late-type galaxies within a fixed circular aperture in each cluster. We find no clear correlation between the fractions of ram-pressure stripping candidates and the different cluster dynamical state proxies considered. These fractions also show no apparent correlation with cluster mass. To construct a dynamical state classification closer to a merging "sequence", we perform a visual classification of the dynamical states of the clusters, combining information available in optical, X-ray, and radio wavelengths. We find a mild increase in the RPS fraction in interacting clusters with respect to all other classes (including post-mergers). This mild enhancement could hint at a short-lived enhanced ram-pressure stripping in ongoing cluster mergers. However, our results are not statistically significant due to the low galaxy numbers. We note this is the first homogeneous attempt to quantify the effect of cluster dynamical state on ram-pressure stripping using a large cluster sample, but even larger (especially wider) multi-wavelength surveys are needed to confirm the results.

A key to understanding the formation of the first galaxies is to quantify the content of the molecular gas as the fuel for star formation activity through the epoch of reionization. In this paper, we use the 158$\mu$m [CII] fine-structure emission line as a tracer of the molecular gas in the interstellar medium (ISM) in a sample of $z=6.5-7.5$ galaxies recently unveiled by the Reionization Era Bright Line Emission Survey, REBELS, with the Atacama Large Millimeter/submillimeter Array. We find substantial amounts of molecular gas ($\sim10^{10.5}\ M_\odot$) comparable to those found in lower redshift galaxies for similar stellar masses ($\sim10^{10}\ M_\odot$). The REBELS galaxies appear to follow the standard scaling relations of molecular gas to stellar mass ratio ($\mu_{\rm mol}$) and gas depletion timescale ($t_{\rm dep}$) with distance to the star-forming main-sequence expected from extrapolations of $z\sim1-4$ observations. We find median values at $z\sim7$ of $\mu_{\rm mol}=2.6_{-1.4}^{4.1}$ and $t_{\rm dep}=0.5_{-0.14}^{+0.26}$ Gyr, indicating that the baryonic content of these galaxies is gas-phase dominated and little evolution from $z\sim7$ to 4. Our measurements of the cosmic density of molecular gas, log$(\rho_{\rm mol}/(M_\odot {\rm Mpc}^{-3}))=6.34^{+0.34}_{-0.31}$, indicate a steady increase by an order of magnitude from $z\sim7$ to 4.

Magnetic fields likely play an important role in star formation, but the number of directly measured magnetic field strengths remains scarce. We observed the 38.3 and 38.5 GHz Class II methanol (CH$_3$OH) maser lines toward the high mass star forming region NGC 6334F for the Zeeman effect. The observed spectral profiles have two prominent velocity features which can be further decomposed through Gaussian component fitting. In several of these fitted Gaussian components we find significant Zeeman detections, with $zB_{\rm los}$ in the range from 8 to 46 Hz. If the Zeeman splitting factor $z$ for the 38 GHz transitions is of the order of $\sim$1 Hz mG$^{-1}$, similar to that for several other CH$_3$OH maser lines, then magnetic fields in the regions traced by these masers would be in the range of 8-46 mG. Such magnetic field values in high mass star forming regions agree with those detected in the better-known 6.7 GHz Class II CH$_3$OH maser line. Since Class II CH$_3$OH masers are radiatively pumped close to the protostar and likely occur in the accretion disk or the interface between the disk and outflow regions, such fields likely have significant impact on the dynamics of these disks.

Intensity mapping observations measure galaxy clustering fluctuations from spectral-spatial maps, requiring stable noise properties on large angular scales. We have developed specialized readouts and analysis methods for achieving large-scale noise stability with Teledyne 2048$\times$2048 H2RG infrared detector arrays. We designed and fabricated a room-temperature low-noise ASIC Video8 amplifier to sample each of the 32 detector outputs continuously in sample-up-the-ramp mode with interleaved measurements of a stable reference voltage that remove current offsets and $1/f$ noise from the amplifier. The amplifier addresses rows in an order different from their physical arrangement on the array, modulating temporal $1/f$ noise in the H2RG to high spatial frequencies. Finally, we remove constant signal offsets in each of the 32 channels using reference pixels. These methods will be employed in the upcoming SPHEREx orbital mission that will carry out intensity mapping observations in near-infrared spectral maps in deep fields located near the ecliptic poles. We also developed a noise model for the H2RG and Video8 to optimize the choice of parameters. Our analysis indicates that these methods hold residual $1/f$ noise near the level of SPHEREx photon noise on angular scales smaller than $\sim30$ arcminutes.

The James Webb Space Telescope (JWST) has provided the first opportunity to study the atmospheres of terrestrial exoplanets and estimate their surface conditions. Earth-sized planets around Sun-like stars are currently inaccessible with JWST however, and will have to be observed using the next generation of telescopes with direct imaging capabilities. Detecting active volcanism on an Earth-like planet would be particularly valuable as it would provide insight into its interior, and provide context for the commonality of the interior states of Earth and Venus. In this work we used a climate model to simulate four exoEarths over eight years with ongoing large igneous province eruptions with outputs ranging from 1.8-60 Gt of sulfur dioxide. The atmospheric data from the simulations were used to model direct imaging observations between 0.2-2.0 $\mu$m, producing reflectance spectra for every month of each exoEarth simulation. We calculated the amount of observation time required to detect each of the major absorption features in the spectra, and identified the most prominent effects that volcanism had on the reflectance spectra. These effects include changes in the size of the O$_3$, O$_2$, and H$_2$O absorption features, and changes in the slope of the spectrum. Of these changes, we conclude that the most detectable and least ambiguous evidence of volcanism are changes in both O$_3$ absorption and the slope of the spectrum.

This work compares cosmological matching conditions used in approximating generic pre-inflationary phases of the universe. We show that the joining conditions for primordial scalar perturbations assumed by Contaldi et al. are inconsistent with the physically motivated Israel junction conditions, however, performing general relativistic matching with the aforementioned constraints results in unrealistic primordial power spectra. Eliminating the need for ambiguous matching, we look at an alternative semi-analytic model for producing the primordial power spectrum allowing for finite duration cosmological phase transitions.

Aims. We provide an in-depth analysis of the COSMOS-Web ring, an Einstein ring at z=2 that we serendipitously discovered in the COSMOS-Web survey and possibly the most distant lens discovered to date. Methods. We extract the visible and NIR photometry from more than 25 bands and we derive the photometric redshifts and physical properties of both the lens and the source with three different SED fitting codes. Using JWST/NIRCam images, we also produce two lens models to (i) recover the total mass of the lens, (ii) derive the magnification of the system, (iii) reconstruct the morphology of the lensed source, and (iv) measure the slope of the total mass density profile of the lens. Results. The lens is a very massive and quiescent (sSFR < 10^(-13) yr-1) elliptical galaxy at z = 2.02 \pm 0.02 with a total mass Mtot(<thetaE) = (3.66 \pm 0.36) x 10^11 Msun and a stellar mass M* = (1.37 \pm 0.14) x 10^11 Msun. Compared to SHMRs from the literature, we find that the total mass is consistent with the presence of a DM halo of mass Mh = 1.09^(+1.46)_(-0.57) x 10^13 Msun. In addition, the background source is a M* = (1.26 \pm 0.17) x 10^10 Msun star-forming galaxy (SFR=(78 \pm 15) Msun/yr) at z = 5.48 \pm 0.06. Its reconstructed morphology shows two components with different colors. Dust attenuation values from SED fitting and nearby detections in the FIR also suggest it could be partially dust-obscured. Conclusions. We find the lens at z=2. Its total, stellar, and DM halo masses are consistent within the Einstein ring, so we do not need any unexpected changes in our description of the lens (e.g. change its IMF or include a non-negligible gas contribution). The most likely solution for the lensed source is at z = 5.5. Its reconstructed morphology is complex and highly wavelength dependent, possibly because it is a merger or a main sequence galaxy with a heterogeneous dust distribution.

Dark Matter (DM) subhalos are predicted to perturb stellar streams; stream morphologies and dynamics can constrain the mass distribution of subhalos. Using FIRE-2 simulations of Milky Way-mass galaxies, we show that presence of a Large Magellanic Cloud (LMC)--analog significantly changes stream-subhalo encounter rates. Three key factors drive these changes. First, the LMC--analog brings in many subhalos, increasing encounter rates for streams near the massive satellite by up to 20--40%. Second, the LMC--analog displaces the host from its center-of-mass (inducing reflex motion), causing a north-south asymmetry in the density and radial velocity distribution of subhalos. This asymmetry results in encounter rates varying by 50--70% across the sky at the same distance. Finally, the LMC--mass satellite induces a density wake in the host's DM halo, further boosting the encounter rates near the LMC--analog. We also explore the influence of stream orbital properties, finding a 50% increase in encounters for streams moving retrograde to the LMC--analog's orbit in the opposite hemisphere. The dependence of encounter rates on stream location and orbit has important implications for where to search for new streams with spurs and gaps in the Milky Way.

The Roman Space Telescope will have the first advanced coronagraph in space, with deformable mirrors for wavefront control, low-order wavefront sensing and maintenance, and a photon-counting detector. It is expected to be able to detect and characterize mature, giant exoplanets in reflected visible light. Over the past decade the performance of the coronagraph in its flight environment has been simulated with increasingly detailed diffraction and structural/thermal finite element modeling. With the instrument now being integrated in preparation for launch within the next few years, the present state of the end-to-end modeling is described, including the measured flight components such as deformable mirrors. The coronagraphic modes are thoroughly described, including characteristics most readily derived from modeling. The methods for diffraction propagation, wavefront control, and structural and thermal finite-element modeling are detailed. The techniques and procedures developed for the instrument will serve as a foundation for future coronagraphic missions such as the Habitable Worlds Observatory.

The Baryon Acoustic Oscillations (BAO) are one of the most used probes to understand the accelerated expansion of the Universe. Traditional methods rely on fiducial model information within their statistical analysis, which may be a problem when constraining different families of models. The aim of this work is to provide a method that constrains $\theta_{BAO}$ through a model-independent and compare parameter estimation of the angular correlation function polynomial approach, using the covariance matrix from the galaxy sample from thin redshift bins, with the usual mock sample covariance matrix. We proposed a different approach to finding the BAO angular feature revisiting previous work in the literature, we take the bias between the correlation function between the bins and the whole sample. We used widths of $\delta z = 0.002$ separation for all samples as the basis for a sample covariance matrix weighted by the statistical importance of the redshift bin. We propose a different weighting scheme based only on random pair counting. We also propose an alternate shift parameter based only on the data. Each sample belongs to the Sloan Digital Sky Survey Luminous Red Galaxies (LRG): BOSS1, BOSS2, and eBOSS, with effective redshift $z_{eff}$: 0.35, 0.51, 0.71, respectively, and different numbers of bins with 50, 100, and 200 respectively. In addition, we correct the angular separation from the polynomial fit ($\theta_{fit}$) that encodes the BAO feature with a bias function obtained by comparing each bin correlation function with the correlation function of the whole set. We also tested the same correction choosing the bin at $z_{eff}$ and found that for eBOSS $\theta_{BAO}$ is in $1 \sigma$ agreement with the Planck 18 model. BOSS1 and BOSS2 $\theta_{BAO}$ agreed in $1\sigma$ with the Pantheon+ & S$H_0$ES Flat$\Lambda$CDM model, in tension with Planck 18.

A wealth of observations have long suggested that the vast majority of isolated classical dwarf galaxies ($M_*=10^7$-$10^9$ M$_\odot$) are currently star-forming. However, recent observations of the large abundance of "Ultra-Diffuse Galaxies" beyond the reach of previous large spectroscopic surveys suggest that our understanding of the dwarf galaxy population may be incomplete. Here we report the serendipitous discovery of an isolated quiescent dwarf galaxy in the nearby Universe, which was imaged as part of the PEARLS GTO program. Remarkably, individual red-giant branch stars are visible in this near-IR imaging, suggesting a distance of $31$ Mpc, and a wealth of archival photometry point to an sSFR of $2\times10^{-12}$ yr$^{-1}$. Spectra obtained with the Lowell Discovery Telescope find a recessional velocity consistent with the Hubble Flow and ${>}1500$ km/s separated from the nearest massive galaxy in SDSS, suggesting that this galaxy was either quenched from internal mechanisms or had a very high-velocity interaction with a nearby massive galaxy in the past. This analysis highlights the possibility that many nearby quiescent dwarf galaxies are waiting to be discovered and that JWST has the potential to identify them.

The low mean densities of sub-Neptunes imply that they formed within a few million years and accreted primordial envelopes. Because these planets receive a total X-ray and extreme ultra-violet flux that is comparable to the gravitational binding energy of their envelopes, their primordial hydrogen-helium atmospheres are susceptible to mass loss. Models of photoevaporating sub-Neptunes have so far assumed that envelope compositions remain constant over time. However, preferential loss of atmospheric hydrogen has the potential to change their compositions. Here, by modeling the thermal and compositional evolution of sub-Neptunes undergoing atmospheric escape with diffusive separation between hydrogen and helium, we show that planets with radii between 1.6 and 2.5 that of Earth can become helium-enhanced from billions of years of photoevaporation, obtaining helium mass fractions in excess of 40%. Atmospheric helium enhancement can be detected through transmission spectra, providing a novel observational test for whether atmospheric escape creates the radius valley.

We explore the use of a Neural Ratio Estimator (NRE) to determine the Hubble constant ($H_0$) in the context of time delay cosmography. Assuming a Singular Isothermal Ellipsoid (SIE) mass profile for the deflector, we simulate time delay measurements, image position measurements, and modeled lensing parameters. We train the NRE to output the posterior distribution of $H_0$ given the time delay measurements, the relative Fermat potentials (calculated from the modeled parameters and the measured image positions), the deflector redshift, and the source redshift. We compare the accuracy and precision of the NRE with traditional explicit likelihood methods in the limit where the latter is tractable and reliable, using Gaussian noise to emulate measurement uncertainties in the input parameters. The NRE posteriors track the ones from the conventional method and, while they show a slight tendency to overestimate uncertainties, they can be combined in a population inference without bias.

Much of a protoplanetary disk is thermally controlled by irradiation from the central star. Such a disk, long thought to have a smoothly flaring shape, is unstable to the so-called 'irradiation instability'. But what's the outcome of such an instability? In particular, is it possible that such a disk settles into a shape that is immune to the instability? We combine Athena++ with a simplified thermal treatment to show that passively heated disks settle into a 'staircase' shape. Here, the disk is punctuated by bright rings and dark gaps, with the bright rings intercepting the lion's share of stellar illumination, and the dark gaps hidden in their shadows. The optical surface of such a disk (height at which starlight is absorbed) resembles a staircase. Although our simulations do not have realistic radiative transfer, we use the RADMC3d code to show that this steady state is in good thermal equilibrium. It is possible that realistic disks reach such a state via ways not captured by our simulations. In contrast to our results here, two previous studies have claimed that irradiated disks stay smooth. We show here that they err on different issues. The staircase state, if confirmed by more sophisticated radiative hydrodynamic simulations, has a range of implications for disk evolution and planet formation.

We describe IXPE polarization observations of the Pulsar Wind Nebula (PWN) MSH15-52, the `Cosmic Hand'. We find X-ray polarization across the PWN, with B field vectors generally aligned with filamentary X-ray structures. High significance polarization is seen in arcs surrounding the pulsar and toward the end of the `jet', with polarization degree PD>70%, thus approaching the maximum allowed synchrotron value. In contrast, the base of the jet has lower polarization, indicating a complex magnetic field at significant angle to the jet axis. We also detect significant polarization from PSR B1509-58 itself. Although only the central pulse-phase bin of the pulse has high individual significance, flanking bins provide lower significance detections and, in conjunction with the X-ray image and radio polarization, can be used to constrain rotating vector model solutions for the pulsar geometry.

Vortex fiber nulling (VFN) is a new interferometric technique with the potential to unlock the ability to detect and spectroscopically characterize exoplanets at angular separations smaller than the conventional diffraction limit of $\lambda$/D. In early 2022, a VFN mode was added to the Keck Planet Imager and Characterizer (KPIC) instrument suite on Keck II. VFN operates by adding an azimuthal phase ramp to the incident wavefront so that light from the star at the center of the field is prevented from coupling into a single-mode fiber. One of the key performance goals of VFN is to minimize the ratio of on-axis starlight coupling to off-axis planet coupling, which requires minimizing the wavefront aberrations of light being injected into the fiber. Non-common path aberrations can be calibrated during the daytime and compensated for with the KPIC deformable mirror during nighttime observing. By applying different amplitudes of low-order Zernike modes, we determine which combinations maximize the system performance. Here we present our work developing and testing different procedures to estimate the incident aberrations, both in simulation and on the Keck bench. The current iteration of this calibration algorithm has been used successfully for VFN observing, and there are several avenues for improvement.

The fraction of the contribution from yet-unresolved gamma-ray sources in the Galactic diffuse gamma rays observed by the Tibet air shower array is an important key to interpreting recent multi-messenger observations. This paper shows a surprising fact: no Tibet diffuse events above 398TeV come from the gamma-ray sources newly detected above 100 TeV by LHAASO. Based on this observational fact, the contribution of sources unresolved by LHAASO to the Tibet diffuse events is estimated to be less than 31% above 398TeV with a 99% confidence level. Our result shows that unresolved sources make only a sub-dominant contribution to the Tibet diffuse events above 398 TeV and a large fraction of the events are truly a diffusive nature.

To understand the decaying phase of outbursts in the black hole (BH) X-ray binaries (BH-XRBs), we performed very long general relativistic magneto-hydrodynamic (GRMHD) simulations of a geometrically thin accretion disk around a Kerr BH with slowly rotating matter injected from outside. We thoroughly studied the flow properties, dynamical behavior of the accretion rate, magnetic flux rate, and jet properties during the temporal evolution. Due to the interaction between the thin disk and injected matter, the accretion flow near the BH goes through different phases. The sequence of phases is: soft state $\rightarrow$ soft-intermediate state $\rightarrow$ hard-intermediate state $\rightarrow$ hard state $\rightarrow$ quiescent state. For the accretion rate (and hence the luminosity) to decrease (as observed) in our model, the mass injection should not decay slower than angular momentum injection. We also observed quasi-periodic oscillations (QPOs) in the accretion flow. Throughout the evolution, we observed low-frequency QPOs (~10Hz) and high-frequency QPOs (\sim 200Hz). Our simple unified accretion flow model for state transitions is able to describe outbursts in BH-XRBs.

N6946-BH1 is the first plausible candidate for a failed supernova (SN), a peculiar event in which a massive star disappears without the expected bright SN, accompanied by collapse into a black hole (BH). Following a luminous outburst in 2009, the source experienced a significant decline in optical brightness, while maintaining a persistent infrared (IR) presence. While it was proposed to be a potential failed SN, such behavior has been observed in SN impostor events in nearby galaxies. Here, we present late-time observations of BH1, taken 14 years after disappearance, using JWST's NIRCam and MIRI instruments to probe a never-before-observed region of the object's spectral energy distribution. We show for the first time that all previous observations of BH1 (pre- and post-disappearance) are actually a blend of at least 3 sources. In the near-IR, BH1 is notably fainter than the progenitor but retains similar brightness to its state in 2017. In the mid-IR, the flux appears to have brightened compared to the inferred fluxes from the best-fitting progenitor model. The total luminosity of the source is between 13 - 25% that of the progenitor. We also show that the IR SED appears consistent with PAH features that arise when dust is illuminated by near-ultraviolet radiation. At present, the interpretation of N6946-BH1 remains uncertain. The observations match expectations for a stellar merger, but theoretical ambiguity in the failed SN hypothesis makes it hard to dismiss.

Gravity drives the collapse of molecular clouds through which stars form, yet the exact role of gravity in cloud collapse remains a complex issue. Studies point to a picture where star formation occurs in clusters. In a typical, pc-sized cluster-forming region, the collapse is hierarchical, and the stars should be born from regions of even smaller sizes ($\approx 0.1\;\rm pc$). The origin of this spatial arrangement remains under investigation. Based on a high-quality surface density map towards the Perseus region, we construct a 3D density structure, compute the gravitational potential, and derive eigenvalues of the tidal tensor ($\lambda_{\rm min}$, $\lambda_{\rm mid}$, $\lambda_{\rm max}$, $\lambda_{\rm min} < \lambda_{\rm mid} < \lambda_{\rm max}$), analyze the behavior of gravity at every location and reveal its multiple roles in cloud evolution. We find that fragmentation is limited to several isolated, high-density ``islands''. Surrounding them, is a vast amount of gas ($75 \%$ of the mass, $95 \%$ of the volume) that stays under the influence of extensive tides where fragmentation is suppressed. This gas will be transported towards these regions to fuel star formation. The spatial arrangement of regions under different tides explains the hierarchical and localized pattern of star formation inferred from the observations. Tides were first recognized by Newton, yet this is the first time its dominance in cloud evolution has been revealed. We expect this link between cloud density structure and role gravity to be strengthened by future studies, resulting in a clear view of the star formation process.

We report a detailed study of the magnetic-field structure of the Crab pulsar wind nebula, using the X-ray polarization data in 2--8~keV obtained with the Imaging X-ray Polarimetry Explorer. Contamination of the pulsar emission to the data of the nebula region was removed through application of a stringent pulsation phase-cut, extracting a phase range of 0.7--1.0 only. We found that the electric field vector polarization angle (PA) was about $130^{\circ}$ from north to east with the polarization degree (PD) of about 25\% at the pulsar position, indicating that the direction of the toroidal magnetic field is perpendicular to the pulsar spin axis in the region close to the termination shock. The PA gradually deviated from the angle as an increasing function of the distance from the pulsar. There was a region of a low PD to the west of the X-ray torus. Although such a region is expected to be located at the torus edge, where geometrical depolarization due to a steep spatial variation of the PA is expected, the observed low-PD region positionally deviated from the edge. We found that the region of low PD positionally coincided with a dense filament seen in the optical band, and conjecture that the low-PD region may be produced through deflection of the pulsar wind. By comparing the values of the PD at the pulsar position between the data and a model, in which toroidal and turbulent magnetic fields were considered, we estimated the fractional energy of the turbulent magnetic field to be about $2/3$ of the total. We also evaluated a potential polarization of the northern jet in the nebula and derived the PD and PA to be about $30\%$ and $120^{\circ}$, respectively.

The knowledge of the positions of the corotation resonance in spiral arms is a key way to estimate their pattern speed, which is a fundamental parameter determining the galaxy dynamics. Various methods for its estimation have been developed, but they all demonstrate certain limitations and a lack of agreement with each other. Here, we present a new method for estimating the corotation radius. This method takes into account the shape of the profile across the arm and its width and, thus, only photometric data is needed. The significance of the method is that it can potentially be used for the farthest galaxies with measurable spiral arms. We apply it to a sample of local galaxies from Savchenko et al. and compare the obtained corotation radii with those previously measured in the literature by other methods. Our results are in good agreement with the literature. We also apply the new method to distant galaxies from the COSMOS field. For the first time, corotation locations for galaxies with photometric redshifts up to $z\sim0.9$ are measured.

We present extended point spread function (PSF) models for the Hyper Suprime-Cam Subaru Strategic Program Public Data Release 3 (HSC-SSP PDR3) in all $\textit{g,r,i,Z}$ and $\textit{Y}$-bands. Due to its 8.2m primary mirror and long exposure periods, HSC combines deep images with wide-field coverage, making it one of the most suitable observing facilities for low surface brightness (LSB) studies. By applying a median stacking technique of point sources with different brightnesses, we show how to construct the HSC-SSP PDR3 PSF models to an extent of R $\sim$ 5.6 arcmin. These new PSFs provide the community with a crucial tool to characterise LSB properties at large angles. We apply our HSC PSFs and demonstrate that they behave reasonably in two cases: first, to generate a 2-D model of a bright star, and second, to remove the PSF-scattered light from an Ultra Deep image of the 400020 Galaxy And Mass Assembly (GAMA) group in the SXDS field. Our main focus in this second application is characterising the $\textit{r}$-band intra-halo light (IHL) component of 400020. Building on advanced source extraction techniques with careful consideration of PSF flux, we measure the IHL surface brightness (SB) group profile up to $\sim$ 31 mag arcsec$^{-2}$ and R = 300 kpc. We estimate the IHL fraction ($\mathrm{f_{IHL}}$) profile, with a mean of $\mathrm{f_{IHL}}$ $\sim$ 0.13. Our results show that not removing the PSF light can overestimate the IHL SB by $\sim$ 1.7 mag arcsec$^{-2}$ and the $\mathrm{f_{IHL}}$ by $\sim$ 30%.

Interest to astrophysical evidence for primordial black holes (PBHs) formed in the early Universe from initial cosmological perturbations has increased after the discovery of coalescing binary black holes with masses more than dozen solar ones by gravitational-wave (GW) observatories. We briefly discuss increasing evidence that PBHs can provide some fraction of detected merging binary BHs and can be related to an isotropic stochastic GW background recently discovered by pulsar timing arrays. We focus on PBHs with log-normal mass spectrum originated from isocurvature perturbations in the modified Affleck-Dine baryogenesis scenario by Dolgov and Silk (1993). We show that almost equal populations of astrophysical binary BHs from massive binary evolution and binary PBHs with log-normal mass spectrum can describe both the observed chirp mass distribution and effective spin -- mass ratio anti-correlation of the LVK binary BHs.

The Pierre Auger Collaboration has embraced the concept of open access to their research data since its foundation, with the aim of giving access to the widest possible community. A gradual process of release began as early as 2007 when 1% of the cosmic-ray data was made public, along with 100% of the space-weather information. In February 2021, a portal was released containing 10% of cosmic-ray data collected from 2004 to 2018, during Phase I of the Observatory. The Portal included detailed documentation about the detection and reconstruction procedures, analysis codes that can be easily used and modified and, additionally, visualization tools. Since then the Portal has been updated and extended. In 2023, a catalog of the 100 highest-energy cosmic-ray events examined in depth has been included. A specific section dedicated to educational use has been developed with the expectation that these data will be explored by a wide and diverse community including professional and citizen-scientists, and used for educational and outreach initiatives. This paper describes the context, the spirit and the technical implementation of the release of data by the largest cosmic-ray detector ever built, and anticipates its future developments.

We present a detailed study of cool-core systems in a sample of four galaxy clusters (RXCJ1504.1-0248, A3112, A4059, and A478) using archival X-ray data from the Chandra X-ray Observatory. Cool cores are frequently observed at the centers of galaxy clusters and are considered to be formed by radiative cooling of the intracluster medium (ICM). Cool cores are characterized by a significant drop in the ICM temperature toward the cluster center. We extract and analyze X-ray spectra of the ICM to measure the radial profiles of the ICM thermodynamic properties including temperature, density, pressure, entropy, and radiative cooling time. We define the cool-core radius as the turnover radius in the ICM temperature profile and investigate the relation between the cool-core radius and the properties of the host galaxy clusters. In our sample, we observe that the radiative cooling time of the ICM at the cool-core radius exceeds 10\,Gyr, with RXCJ1504.1-0248 exhibiting a radiative cooling time of $32^{+5}_{-11}$\,Gyr at its cool-core radius. These results indicate that not only radiative cooling but also additional mechanisms such as gas sloshing may play an important role in determining the size of cool cores. Additionally, we find that the best-fit relation between the cool-core radius and the cluster mass ($M_{500}$) is consistent with a linear relation. Our findings suggest that cool cores are linked to the evolution of their host galaxy clusters.

The Baikal-GVD is a deep-underwater neutrino telescope being constructed in Lake Baikal. After the winter 2023 deployment campaign the detector consists of 3456 optical modules installed on 96 vertical strings. The status of the detector and progress in data analysis are discussed in present report. The Baikal-GVD data collected in 2018-2022 indicate the presence of cosmic neutrino flux in high-energy cascade events consistent with observations by the IceCube neutrino telescope. Analysis of track-like events results in identification of first high-energy muon neutrino candidates. These and other results from 2018-2022 data samples are reviewed in this report.

Light curves of stars encapsulate a wealth of information about stellar oscillations and granulation, thereby offering key insights into the internal structure and evolutionary state of stars. Conventional asteroseismic techniques have been largely confined to power spectral analysis, neglecting the valuable phase information contained within light curves. While recent machine learning applications in asteroseismology utilizing Convolutional Neural Networks (CNNs) have successfully inferred stellar attributes from light curves, they are often limited by the local feature extraction inherent in convolutional operations. To circumvent these constraints, we present $\textit{Astroconformer}$, a Transformer-based deep learning framework designed to capture long-range dependencies in stellar light curves. Our empirical analysis, which focuses on estimating surface gravity ($\log g$), is grounded in a carefully curated dataset derived from $\textit{Kepler}$ light curves. These light curves feature asteroseismic $\log g$ values spanning from 0.2 to 4.4. Our results underscore that, in the regime where the training data is abundant, $\textit{Astroconformer}$ attains a root-mean-square-error (RMSE) of 0.017 dex around $\log g \approx 3 $. Even in regions where training data are sparse, the RMSE can reach 0.1 dex. It outperforms not only the K-nearest neighbor-based model ($\textit{The SWAN}$) but also state-of-the-art CNNs. Ablation studies confirm that the efficacy of the models in this particular task is strongly influenced by the size of their receptive fields, with larger receptive fields correlating with enhanced performance. Moreover, we find that the attention mechanisms within $\textit{Astroconformer}$ are well-aligned with the inherent characteristics of stellar oscillations and granulation present in the light curves.

The statistics of thermal gas pressure are a new and promising probe of cosmology and astrophysics. The large-scale cross-correlation between galaxies and the thermal Sunyaev-Zeldovich effect gives the bias-weighted mean electron pressure, $\langle b_\mathrm{h}P_e\rangle$. In this paper, we show that $\langle b_\mathrm{h}P_e\rangle$ is sensitive to the amplitude of fluctuations in matter density, for example $\langle b_\mathrm{h}P_e\rangle\propto \left(\sigma_8\Omega_\mathrm{m}^{0.81}h^{0.67}\right)^{3.14}$ at redshift $z=0$. We find that at $z<0.5$ the observed $\langle b_\mathrm{h}P_e\rangle$ is smaller than that predicted by the state-of-the-art hydrodynamical simulations of galaxy formation, MillenniumTNG, by a factor of $0.93$. This can be explained by a lower value of $\sigma_8$ and $\Omega_\mathrm{m}$, similar to the so-called "$S_8$ tension'' seen in the gravitational lensing effect, although the influence of astrophysics cannot be completely excluded. The difference between Magneticum and MillenniumTNG at $z<2$ is small, indicating that the difference in the galaxy formation models used by these simulations has little impact on $\langle b_\mathrm{h}P_e\rangle$ at this redshift range. At higher $z$, we find that both simulations are in a modest tension with the existing upper bounds on $\langle b_\mathrm{h}P_e\rangle$. We also find a significant difference between these simulations there, which we attribute to a larger sensitivity to the galaxy formation models in the high redshift regime. Therefore, more precise measurements of $\langle b_\mathrm{h}P_e\rangle$ at all redshifts will provide a new test of our understanding of cosmology and galaxy formation.

Emission in many astrophysical transients originates from a shocked fluid. A central engine typically produces an outflow with varying speeds, leading to internal collisions within the outflow at finite distances from the source. Each such collision produces a pair of forward and reverse shocks with the two shocked regions separated by a contact discontinuity (CD). As a useful approximation, we consider the head-on collision between two cold and uniform shells (a slower leading shell and a faster trailing shell) of finite radial width, and study the dynamics of shock propagation in planar geometry. We find significant differences between the forward and reverse shocks, in terms of their strength, internal energy production efficiency, and the time it takes for the shocks to sweep through the respective shells. We consider the subsequent propagation of rarefaction waves in the shocked regions and explore the cases where these waves can catch up with the shock fronts and thereby limit the internal energy dissipation. We demonstrate the importance of energy transfer from the trailing to leading shell through $pdV$ work across the CD. We outline the parameter space regions relevant for models of different transients,e.g., Gamma-ray burst (GRB) internal shock model, fast radio burst (FRB) blastwave model, Giant flare due to magnetars, and superluminous supernovae (SLSN) ejecta. We find that the reverse shock likely dominates the internal energy production for many astrophysical transients.

We report on the SRG/eROSITA detection of ultra-soft ($kT=47^{+5}_{-5}$ eV) X-ray emission ($L_{\mathrm{X}}=2.5^{+0.6}_{-0.5} \times 10^{43}$ erg s$^{-1}$) from the tidal disruption event (TDE) candidate AT 2022dsb $\sim$14 days before peak optical brightness. As the optical luminosity increases after the eROSITA detection, then the 0.2--2 keV observed flux decays, decreasing by a factor of $\sim 39$ over the 19 days after the initial X-ray detection. Multi-epoch optical spectroscopic follow-up observations reveal transient broad Balmer emission lines and a broad He II 4686A emission complex with respect to the pre-outburst spectrum. Despite the early drop in the observed X-ray flux, the He II 4686A complex is still detected for $\sim$40 days after the optical peak, suggesting the persistence of an obscured, hard ionising source in the system. Three outflow signatures are also detected at early times: i) blueshifted H$\alpha$ emission lines in a pre-peak optical spectrum, ii) transient radio emission, and iii) blueshifted Ly$\alpha$ absorption lines. The joint evolution of this early-time X-ray emission, the He II 4686A complex and these outflow signatures suggests that the X-ray emitting disc (formed promptly in this TDE) is still present after optical peak, but may have been enshrouded by optically thick debris, leading to the X-ray faintness in the months after the disruption. If the observed early-time properties in this TDE are not unique to this system, then other TDEs may also be X-ray bright at early times and become X-ray faint upon being veiled by debris launched shortly after the onset of circularisation.

Context. The heliosphere is formed by the interaction between the solar wind (SW) plasma emanating from the Sun and a magnetised component of local interstellar medium (LISM) inflowing on the Sun. A separation surface called the heliopause (HP) forms between the SW and the LISM. Aims. In this article, we define the nose of the HP and investigate the variations in its location. These result from a dependence on the intensity and direction of the interstellar magnetic field (ISMF), which is still not well known but has a significant impact on the movement of the HP nose, as we try to demonstrate in this paper. Methods. We used a parametric study method based on numerical simulations of various forms of the heliosphere using a time-dependent three-dimensional magnetohydrodynamic (3D MHD) model of the heliosphere. Results. The results confirm that the nose of the HP is always in a direction that is perpendicular to the maximum ISMF intensity directly behind the HP. The displacement of the HP nose depends on the direction and intensity of the ISMF, with the structure of the heliosphere and the shape of the HP depending on the 11-year cycle of solar activity. Conclusions. In the context of the planned space mission to send the Interstellar Probe (IP) to a distance of 1000 AU from the Sun, our study may shed light on the question as to which direction the IP should be sent. Further research is needed that introduces elements such as current sheet, reconnection, cosmic rays, instability, or turbulence into the models.

Solar irradiance variability has been monitored almost exclusively from the Earth's perspective. {We present a method to combine the unprecedented observations of the photospheric magnetic field and continuum intensity from outside the Sun-Earth line, which is being recorded by the Polarimetric and Helioseismic Imager on board the Solar Orbiter mission (SO/PHI), with solar observations recorded from the Earth's perspective to examine the solar irradiance variability from both perspectives simultaneously.} Taking SO/PHI magnetograms and continuum intensity images from the cruise phase of the Solar Orbiter mission and concurrent observations from the Helioseismic and Magnetic Imager onboard the Solar Dynamics Observatory (SDO/HMI) as input into the SATIRE-S model, we successfully reconstructed the total solar irradiance variability as apparent from both perspectives. In later stages of the SO mission, the orbital plane will tilt in such a way as to bring the spacecraft away from the ecliptic to heliographic latitudes of up to $33^{\circ}$. The current study sets the template for the reconstruction of solar irradiance variability as seen from outside the ecliptic from data that SO/PHI is expected to collect from such positions. {Such a reconstruction will be beneficial to factoring inclination into how the brightness variations of the Sun compare to those of other cool stars, whose rotation axes are randomly inclined.

In this study, we investigate the formation of primordial black holes (PBHs) in a scalar field inflationary model coupled to the Gauss-Bonnet (GB) term with fractional power-law potentials. The coupling function enhances the curvature perturbations, then results in the generation of PBHs and detectable secondary gravitational waves (GWs). % We identify three separate sets of parameters for the potential functions of the form $\phi^{1/3}$, $\phi^{2/5}$, and $\phi^{2/3}$. By adjusting the model parameters, we decelerate the inflaton during the ultra slow-roll (USR) phase and enhance curvature perturbations. % Our calculations predict the formation of PBHs with masses of ${\cal O}(10)M_{\odot}$, which are compatible with LIGO-Virgo observational data. Additionally, we find PBHs with masses around ${\cal O}(10^{-6})M_{\odot}$ and ${\cal O}(10^{-5})M_{\odot}$, which can explain ultrashort-timescale microlensing events in OGLE data. % Furthermore, our proposed mechanism could lead to the formation of PBHs in mass scales around ${\cal O}(10^{-14})M_{\odot}$ and ${\cal O}(10^{-13})M_{\odot}$, contributing to approximately 99\% of the dark matter in the universe. % We also study the production of secondary GWs in our model. In all cases of the model, the density parameter of secondary GWs $\Omega_{\rm GW_0}$ exhibits peaks that intersect the sensitivity curves of GWs detectors, providing a means to verify our findings using data of these detectors. % Our numerical results demonstrate a power-law behavior for the spectra of $\Omega_{\rm GW_0}$ with respect to frequency, given by $\Omega_{\rm GW_0} (f) \sim (f/f_c)^{n}$. Additionally, in the infrared regime where $f\ll f_{c}$, the power index takes a log-dependent form, specifically $n=3-2/\ln(f_c/f)$.

PSR J1953+1844 (i.e. M71E) is a millisecond pulsar (MSP) in a 53-min binary orbit discovered by the Five-hundred-meter Aperture Spherical radio Telescope (FAST). It has an orbital period of 53 minutes and a mass function of $2.3\times10^{-7}$ M$_\odot$. The possible redback origin of this system has been discussed by Pan et al. (2023). We discuss here an alternative evolution track for this binary system as being a descendant of ultra-compact X-ray binaries (UCXBs), which has a hydrogen-poor donor accreting onto a neutron-star with an orbital period of $\leq1~hr$. Noticed that some of UCXB systems hold an accreting millisecond X-ray pulsars (AMXPs) and a $\sim$0.01 M$_\odot$ donor. M71E has a very similar orbit as AMXPs, indicating that it might be evolved from an UCXB as PSR J1653--0158. The companion star should be significantly bloated and it most probably has a carbon and oxygen composition, otherwise a low inclination angle is required for a helium companion. The discovery of this binary system may shed lights on when and how a NS in UCXBs turns to be a radio pulsar.

The IceCube Neutrino Observatory is a cubic-kilometer high-energy neutrino detector deployed in the Antarctic ice. Two major event classes are charged-current electron and muon neutrino interactions. In this contribution, we discuss the inference of direction and energy for these classes using conditional normalizing flows. They allow to derive a posterior distribution for each individual event based on the raw data that can include systematic uncertainties, which makes them very promising for next-generation reconstructions. For each normalizing flow we use the differential entropy and the KL-divergence to its maximum entropy approximation to interpret the results. The normalizing flows correctly incorporate complex optical properties of the Antarctic ice and their relation to the embedded detector. For showers, the differential entropy increases in regions of high photon absorption and decreases in clear ice. For muons, the differential entropy strongly correlates with the contained track length. Coverage is maintained, even for low photon counts and highly asymmetrical contour shapes. For high-photon counts, the distributions get narrower and become more symmetrical, as expected from the asymptotic theorem of Bernstein-von-Mises. For shower directional reconstruction, we find the region between 1 TeV and 100 TeV to potentially benefit the most from normalizing flows because of azimuth-zenith asymmetries which have been neglected in previous analyses by assuming symmetrical contours. Events in this energy range play a vital role in the recent discovery of the galactic plane diffuse neutrino emission.

We present and evaluate the prospects for detecting coherent radio counterparts to gravitational wave (GW) events using Murchison Widefield Array (MWA) triggered observations. The MWA rapid-response system, combined with its buffering mode ($\sim4$ minutes negative latency), enables us to catch any radio signals produced from seconds prior to hours after a binary neutron star (BNS) merger. The large field of view of the MWA ($\sim1000\,\text{deg}^2$ at 120\,MHz) and its location under the high sensitivity sky region of the LIGO-Virgo-KAGRA (LVK) detector network, forecast a high chance of being on-target for a GW event. We consider three observing configurations for the MWA to follow up GW BNS merger events, including a single dipole per tile, the full array, and four sub-arrays. We then perform a population synthesis of BNS systems to predict the radio detectable fraction of GW events using these configurations. We find that the configuration with four sub-arrays is the best compromise between sky coverage and sensitivity as it is capable of placing meaningful constraints on the radio emission from 12.6\% of GW BNS detections. Based on the timescales of four BNS merger coherent radio emission models, we propose an observing strategy that involves triggering the buffering mode to target coherent signals emitted prior to, during or shortly following the merger, which is then followed by continued recording for up to three hours to target later time post-merger emission. We expect MWA to trigger on $\sim5\text{--}22$ BNS merger events during the LVK O4 observing run, which could potentially result in two detections of predicted coherent emission.

Earth likely acquired much of its inventory of volatile elements during the main stage of its formation. Some of Earth's proto-atmosphere must therefore have survived the giant impacts, collisions between planet-sized bodies, that dominate the latter phases of accretion. Here we use a suite of 1D hydrodynamic simulations and impedance match calculations to quantify the effect that pre-impact surface conditions (such as atmospheric pressure and presence of an ocean) have on the efficiency of atmospheric and ocean loss from proto-planets during giant impacts. We find that -- in the absence of an ocean -- lighter, hotter, and lower-pressure atmospheres are more easily lost. The presence of an ocean can significantly increase the efficiency of atmospheric loss compared to the no-ocean case, with a rapid transition between low and high loss regimes as the mass ratio of atmosphere to ocean decreases. However, contrary to previous thinking, the presence of an ocean can also reduce atmospheric loss if the ocean is not sufficiently massive, typically less than a few times the atmospheric mass. Volatile loss due to giant impacts is thus highly sensitive to the surface conditions on the colliding bodies. To allow our results to be combined with 3D impact simulations, we have developed scaling laws that relate loss to the ground velocity and surface conditions. Our results demonstrate that the final volatile budgets of planets are critically dependent on the exact timing and sequence of impacts experienced by their precursor planetary embryos, making atmospheric properties a highly stochastic outcome of accretion.

The Radio Neutrino Observatory in Greenland (RNO-G) is a radio-based ultra-high energy neutrino detector located at Summit Station, Greenland. It is still being constructed, with 7 stations currently operational. Neutrino detection works by measuring Askaryan radiation produced by neutrino-nucleon interactions. A neutrino candidate must be found amidst other backgrounds which are recorded at much higher rates -- including cosmic-rays and anthropogenic noise -- the origins of which are sometimes unknown. Here we describe a method to classify different noise classes using the latent space of a variational autoencoder. The latent space forms a compact representation that makes classification tractable. We analyze data from a noisy and a silent station. The method automatically detects and allows us to qualitatively separate multiple event classes, including physical wind-induced signals, for both the noisy and the quiet station.

The Milky Way is a barred spiral galaxy with "bar lanes" that bring gas towards the Galactic Center. Gas flowing along these bar lanes often overshoots, and instead of accreting onto the Central Molecular Zone, it collides with the bar lane on the opposite side of the Galaxy. We observed G5, a cloud which we believe is the site of one such collision, near the Galactic Center at (l,b) = (+5.4, -0.4) with the ALMA/ACA. We took measurements of the spectral lines $^{12}$CO J=2-1, $^{13}$CO J=2-1, C$^{18}$O J=2-1, H$_2$CO J=3$_{03}$-2$_{02}$, H$_{2}$CO J=3$_{22}$-2$_{21}$, CH$_{3}$OH J=4$_{22}$-3$_{12}$, OCS J=18-17 and SiO J=5-4. We observed a velocity bridge between two clouds at $\sim$50 km/s and $\sim$150 km/sin our position-velocity diagram, which is direct evidence of a cloud-cloud collision. We measured an average gas temperature of $\sim$60 K in G5 using H$_2$CO integrated intensity line ratios. We observed that the $^{12}$C/$^{13}$C ratio in G5 is consistent with optically thin, or at most marginally optically thick $^{12}$CO. We measured 1.5 x 10$^9$ cm$^{-2}$(K km/s)$^{-1}$ for the local X$_{CO}$, 10-20x less than the average Galactic value. G5 is strong direct observational evidence of gas overshooting the Central Molecular Zone (CMZ) and colliding with a bar lane on the opposite side of the Galactic center.

We describe enhancements to the digest2 software, a short-arc orbit classifier for heliocentric orbits. Digest2 is primarily used by the Near-Earth Object (NEO) community to flag newly discovered objects for a immediate follow-up and has been a part of NEO discovery process for more than 15 years. We have updated the solar system population model used to weight the digest2 score according to the 2023 catalog of known solar system orbits and extended the list of mean uncertainties for 140 observatory codes. Moreover, we have added Astrometry Data Exchange Standard (ADES) input format support to digest2, which provides additional information for the astrometry, such as positional uncertainties for each detection. The digest2 code was also extended to read the roving observer astrometric format as well as the ability to compute a new parameter from the provided astrometric uncertainties ($RMS'$) that can serve as an indicator of in-tracklet curvature when compared with tracklet's great-circle fit RMS. Comparison with the previous version of digest2 confirmed the improvement in accuracy of NEO identification and found that using ADES XML input significantly reduces the computation time of the digest2.

Supernova remnants are the nebular leftover of defunct stellar environments, resulting from the interaction between a supernova blastwave and the circumstellar medium shaped by the progenitor throughout its life. They display a large variety of non-spherical morphologies such as ears that shine non-thermally. % We have modelled the structure and the non-thermal emission of the supernova remnant G1.9+0.3 through 3D magnetohydrodynamic numerical simulations. We propose that the peculiar ear-shaped morphology of this supernova remnant results from the interaction of the its blast wave with a magnetized circumstellar medium, which was previously asymmetrically shaped by the past stellar wind emanating from the progenitor star or its stellar companion. We created synthetic non-thermal radio and x-ray maps from our simulated remnant structure, which are in qualitative agreement with observations, forming ears on the polar directions. Our synthetic map study explains the discrepancies between the measured non-thermal radio and X-ray surface brightness distributions assuming that the Inverse Compton process produces the observed X-ray emission.

The solar flares - which are the most prominent manifestation of the solar activity - typically manifest themselves as a single or a set of luminous arcs (magnetic flux tubes) rooted in the regions of opposite polarity in the photosphere. However, a careful analysis of the archival data by Hinode satellite sometimes reveals the surprising cases of the flaring arcs whose footpoints belong to the regions of the same polarity or to the areas without any appreciable magnetic field. Despite a counterintuitive nature of this phenomenon, it can be reasonably interpreted in the framework of the so-called 'topological model' of magnetic reconnection, where the magnetic null point is formed due to specific superposition of influences from the remote sources rather than by the local current systems. As a result, the energy release propagates along the separator of the flipping two-dome structure rather than along a fixed magnetic field line. Therefore, the luminous arc needs not to be associated anymore immediately with the magnetic sources. Here, we report both the observational cases of the above-mentioned type as well as provide their theoretical model and the numerical simulations.

About one-fifth of the Galactic globular clusters (GCs), dubbed Type II GCs, host distinct stellar populations with different heavy elements abundances. NGC 1851 is one of the most studied Type II GCs, surrounded by several controversies regarding the spatial distribution of its populations and the presence of star-to-star [Fe/H], C+N+O, and age differences. This paper provides a detailed characterization of its stellar populations through Hubble Space Telescope (HST), ground-based, and Gaia photometry. We identified two distinct populations with different abundances of s-process elements along the red-giant branch (RGB) and the sub-giant branch (SGB) and detected two sub-populations among both s-poor (canonical) and s-rich (anomalous) stars. To constrain the chemical composition of these stellar populations, we compared observed and simulated colors of stars with different abundances of He, C, N, and O. It results that the anomalous population has a higher CNO overall abundance compared to the canonical population and that both host stars with different light-element abundances. No significant differences in radial segregation between canonical and anomalous stars are detected, while we find that among their sub-populations, the two most chemical extremes are more centrally concentrated. Anomalous and canonical stars show different 2D spatial distributions outside ~3 arcmin, with the latter developing an elliptical shape and a stellar overdensity in the northeast direction. We confirm the presence of a stellar halo up to ~80 arcmin with Gaia photometry, tagging 14 and five of its stars as canonical and anomalous, respectively, finding a lack of the latter in the south/southeast field.

We present an analysis of the dense gas structures in the immediate surroundings of the massive young stellar object (MYSO) W42-MME, using the high-resolution (0$''$.31$\times$0$''$.25) ALMA dust continuum and molecular line data. We performed a dendrogram analysis of H$^{13}$CO$^{+}$ (4-3) line data to study multi-scale structures and their spatio-kinematic properties, and analyzed the fragmentation and dynamics of dense structures down to $\sim$2000 AU scale. Our results reveal 19 dense gas structures, out of which 12 are leaves and 7 are branches in dendrogram terminology. These structures exhibit transonic-supersonic gas motions (1$<\mathcal{M}<5$) with overvirial states ($\alpha_{\rm vir}\geq2$). The non-thermal velocity dispersion-size relation ($\sigma_{\rm nt}-L$) of dendrogram structures shows a weak negative correlation, while the velocity dispersion across the sky ($\delta\mathit{V_{\rm lsr}}$) correlates positively with structure size ($L$). Velocity structure function ($S_{2}(l)^{1/2}$) analysis of H$^{13}$CO$^{+}$ data reveals strong power-law dependencies with lag ($l$) up to a scale length of $\lesssim$ 6000 AU. The mass-size ($M-R$) relation of dendrogram structures shows a positive correlation with power-law index of 1.73$\pm$0.23, and the leaf L17 hosting W42-MME meets the mass-size conditions for massive star formation. Blue asymmetry is observed in the H$^{12}$CO$^{+}$ (4-3) line profiles of most of the leaves, indicating infall. Overall, our results observationally support the hierarchical and chaotic collapse scenario in the proximity of the MYSO W42-MME.

The Cherenkov Telescope Array (CTA) will be the next generation ground-based very-high-energy gamma-ray observatory, constituted by tens of Imaging Atmospheric Cherenkov Telescopes at two sites once its construction and commissioning are finished. Like its predecessors, CTA relies on Instrument Response Functions (IRFs) to relate the observed and reconstructed properties to the true ones of the primary gamma-ray photons. IRFs are needed for the proper reconstruction of spectral and spatial information of the observed sources and are thus among the data products issued to the observatory users. They are derived from Monte Carlo simulations, depend on observation conditions like the telescope pointing direction or the atmospheric transparency and can evolve with time as hardware ages or is replaced. Producing a complete set of IRFs from simulations for every observation taken is a time-consuming task and not feasible when releasing data products on short timescales. Consequently, interpolation techniques on simulated IRFs are investigated to quickly estimate IRFs for specific observation conditions. However, as some of the IRFs constituents are given as probability distributions, specialized methods are needed. This contribution summarizes and compares the feasibility of multiple approaches to interpolate IRF components in the context of the pyirf python software package and IRFs simulated for the Large-Sized Telescope prototype (LST-1). We will also give an overview of the current functionalities implemented in pyirf.

We use spectra and maps of the $J=1-0$ and $J=2-1$ DCO$^{+}$, DCN, DNC, $\rm N_2D^+$ lines and $1_{11}-1_{01}$ ortho- and para-NH$_{2}$D lines, obtained with the IRAM-30m telescope, as well as observations of their hydrogenated isotopologues to study deuteration processes in five high-mass star-forming regions. The temperature was estimated from CH$_3$CCH lines, also observed with the IRAM-30m telescope, and from NH$_3$ lines, observed with the 100-m radio telescope in Effelsberg, as well as using the integrated intensity ratios of the $J=1-0$ H$^{13}$CN and HN$^{13}$C lines and their main isotopologues. Applying a non-LTE radiative transfer model with RADEX, the gas density and the molecular column densities were estimated. D/H ratios are $0.001-0.05$ for DCO$^{+}$, $0.001-0.02$ for DCN, $0.001-0.05$ for DNC and $0.02-0.4$ for NH$_{2}$D. The D/H ratios decrease with increasing temperature in the range of $\rm 20-40 \,K$ and slightly vary at densities $n(\rm H_2) \sim 10^4-10^6\, cm^{-3}$. The deuterium fraction of $\rm N_2H^{+}$ is $0.008-0.1$ at temperatures in the range of $\rm 20-25\, K$ and at a density of $\sim 10^5\, \rm cm^{-3}$. We also estimate relative abundances and find $ \sim 10^{-11}-10^{-9}$ for DCO$^{+}$ and DNC, $ \sim 10^{-11}-10^{-10}$ for $\rm N_2D^+$ and $ \sim 10^{-10}-10^{-8}$ for NH$_{2}$D. The relative abundances of these species decrease with increasing temperature. However, the DCN/H$_2$ ratio is almost constant ($\sim 10^{-10}$). The observational results agree with the predictions of chemical models (although in some cases there are significant differences).

Previously, we demonstrated that MgII and CIV reverberation-mapped quasars (RM QSOs) are standardizable and that the cosmological parameters inferred using the broad-line region radius-luminosity (R-L) relation are consistent with those determined from better-established cosmological probes. With more data expected from ongoing and future spectroscopic and photometric surveys, it is imperative to examine how new QSO data sets of varied quality, with their own specific luminosity and time-delay distributions, can be best used to determine more restrictive cosmological parameter constraints. In this study, we test the effect of adding 25 OzDES MgII RM QSOs as well as 25 lower-quality SDSS RM CIV QSOs to our previous sample of RM QSOs. Although cosmological parameter constraints become tighter after adding these new QSOs, the new combined data sets have an increased maximum difference in R-L relation parameter values (with a > 2sigma difference in some R-L relation intercepts between cosmological models), and thus a lower standardizability for the larger MgII + CIV compilation. Different time-delay methodologies, particularly the ICCF and CREAM methods used for inferring time delays of SDSS RM QSOs, slightly affect cosmological and R-L relation parameter values, however, the effect is negligible for (smaller) compilations of robust time-delay detections.

We investigate the degree to which current ultrahigh energy cosmic ray observations above the ankle support a common maximum rigidity for all nuclei, often called a Peters cycle, over alternative scenarios for the cosmic ray spectra escaping sources. We show that a Peters cycle is not generally supported by the data when compared with these alternatives. We explore the observational signatures of non-Peters cycle scenarios, and the opportunities to explore both ultrahigh energy cosmic ray source conditions, as well as, physics beyond the Standard model they present.

At 327 MHz, the observed emission of PSR B1937+21 is greatly affected by scattering in the interstellar medium, on a timescale of order the pulse period. We use the bright impulsive giant pulses emitted by the pulsar to measure the impulse response of the interstellar medium and then recover the intrinsic emission of the pulsar by deconvolution -- revealing fine structure on timescales not normally observable. We find that the intrinsic widths of the main pulse and interpulse in the pulse profile are similar to those measured at higher frequencies. We detect 60,270 giant pulses which typically appear as narrow, ~100 ns bursts consisting of one to few nanoshots with widths $\lesssim \! 10$ ns. However, about 10% of the giant pulses exhibit multiple bursts which seem to be causally related to each other. We also report the first detection of giant micropulses in PSR B1937+21, primarily associated with the regular main pulse emission. These are distinct from giant pulses not only in the phases at which they occur, but also in their larger widths, of order a microsecond, and steeper energy distribution. These measurements place useful observational constraints on emission mechanisms for giant pulses as well as the regular radio emission of millisecond pulsars.

In the context of an evolutionary model, the outflow phase of an Active Galactic Nuclei (AGN) occurs at the peak of its activity, once the central SMBH is massive enough to generate sufficient power to counterbalance the potential well of the host galaxy. This phase plays a vital role in galaxy evolution. We aim to apply various selection methods to isolate powerful AGNs in the feedback phase, trace and characterise their outflows, and explore the link between AGN luminosity and outflow properties. We applied a combination of methods to the eROSITA Final Equatorial Depth survey (eFEDS) catalogue and isolated ~1400 candidates at z>0.5 out of ~11750 AGNs (~12\%). We tested the robustness of our selection on the small subsample of 50 sources with available good quality SDSS spectra at 0.5<z<1, for which we fitted the [OIII] emission line complex and searched for the presence of ionised gas outflows. We identified 23 quasars (~45\%) with evidence of ionised outflows based on the presence of significant broad and shifted components in the [OIII] line. They are on average more luminous and more obscured than the parent sample, although this may be ascribed to selection effects affecting the good quality SDSS spectra sample. By adding 118 outflowing quasars at 0.5<z<3.5 from the literature, we find a weak correlation between the maximum outflow velocity and AGN bolometric luminosity. On the contrary, we find strong correlations between mass outflow rate and outflow kinetic power with the AGN bolometric luminosity. About 30\% of our sample have kinetic coupling efficiencies >1\%. We find that the majority of the outflows have momentum flux ratios lower than 20 which rules out an energy-conserving nature. Our present work points to the unequivocal existence of a rather short AGN outflow phase, paving the way towards a new avenue to dissect AGN outflows in large samples within eROSITA and beyond.

Strong gravitational lensing offers constraints on the Hubble constant that are independent of other methods. However, those constraints are subject to uncertainties in lens models. Previous studies suggest that using an elliptical power law + external shear (EPL+XS) for the lensing galaxy can yield results that are precise but inaccurate. We examine such models by generating and fitting mock lenses which produces multiple images of a background quasar-like point source. Despite using the same model for input and output, we find statistical bias in the Hubble constant on the order of 3% to 5%, depending on whether the elliptical lenses have noise or not. The phase space distribution has a `flared' shape that causes the mass power law slope to be underestimated and the Hubble constant to be overestimated. The bias varies with image configuration, which we quantify through annulus length between images with the first and second time delays ($\Delta r_{1,2}$). The statistical bias is worse for configurations that have narrow annuli (e.g., symmetric cross configurations). Assuming a source at redshift 2.0 and an EPL+XS lens at redshift 0.3, we find that the bias can be reduced, but not eliminated, if we limit the sample to systems with annulus lengths $\Delta r_{1,2} \gtrsim 0.3$ arcsec. As lens samples grow, it may be helpful to prioritize this range of image configurations for follow-up observation and analysis.

We revisit the evolution of galaxy morphology in the COSMOS field over the redshift range $0.2\leq z \leq 1$, using a large and complete sample of 33,605 galaxies with a stellar mass of log($M_{\ast}$/M$_{\odot} )>9.5$ with significantly improved redshifts and comprehensive non-parametric morphological parameters. Our sample has 13,881 ($\sim41.3\%$) galaxies with reliable spectroscopic redshifts and has more accurate photometric redshifts with a $\sigma_{\rm NMAD} \sim 0.005$. This paper is the first in a series that investigates merging galaxies and their properties. We identify 3,594 major merging galaxies through visual inspection and find 1,737 massive galaxy pairs with log($M_\ast$/M$_\odot$)$>10.1$. Among the family of non-parametric morphological parameters including $C$, $A$, $S$, $Gini$, $M_{\rm 20}$, $A_{\rm O}$, and $D_{\rm O}$, we find that the outer asymmetry parameter $A_{\rm O}$ and the second-order momentum parameter $M_{\rm 20}$ are the best tracers of merging features than other combinations. Hence, we propose a criterion for selecting candidates of violently star-forming mergers: $M_{\rm 20}> -3A_{\rm O}+3$ at $0.2<z<0.6$ and $M_{\rm 20}> -6A_{\rm O}+3.7$ at $0.6<z<1.0$. Furthermore, we show that both the visual merger sample and the pair sample exhibit a similar evolution in the merger rate at $z<1$, with $\Re \sim(1+z)^{1.79 \pm 0.13}$ for the visual merger sample and $\Re \sim(1+z)^{2.02\pm 0.42}$ for the pair sample. The visual merger sample has a specific star formation rate that is about 0.16\,dex higher than that of non-merger galaxies, whereas no significant star formation excess is observed in the pair sample. This suggests that the effects of mergers on star formation differ at different merger stages.

Evolution of massive stars is dominated by interactions within binary systems. Therefore, it is necessary to investigate all forms of interaction in binary systems that may affect the evolution of the components. One of such `laboratories' is the massive eccentric binary system MACHO$\,$80.7443.1718 (ExtEV). We examine whether the light variability of the ExtEV can be explained by a wind-wind collision (WWC) binary system model. We conducted an analysis of broadband multi-color photometry of ExtEV, time-series space photometry from TESS, ground-based Johnson $UBV$ photometry, and time-series spectroscopy. We fitted an analytical model of light variations to the TESS light curve of ExtEV. We rule out the possibility of the presence of a disk around the primary component. We also argue that the non-linear wave-breaking scenario is not consistent with the observations of ExtEV. We refine the orbital parameters of ExtEV and find evidence for the presence of a tertiary component. Using evolutionary models we demonstrate that the primary component's mass is between 25 and 45$\,$M$_\odot$. We successfully reproduce light curve of ExtEV with our model, showing that the dominant processes shaping its light curve are atmospheric eclipse and light scattered in the WWC cone. We also estimate the primary's mass-loss rate due to stellar wind for $4.5\cdot 10^{-5}\,$M$_\odot\,{\rm yr}^{-1}$. ExtEV is not an extreme eccentric ellipsoidal variable, but an exceptional WWC binary system. The mass loss rate we derived exceeds theoretical predictions by up to two orders of magnitude. This implies that the wind in the system is likely enhanced by tidal interactions, rotation, and possibly also tidally excited oscillations. ExtEV represents a rare evolutionary phase of a binary system that may help to understand the role of a companion-driven enhanced mass loss in the evolution of massive binary systems.

We present the results of new time series photometric observations of 29 pre-white dwarf stars of PG 1159 spectral type, carried out in the years 2014-2022. For the majority of stars, a median noise level in Fourier amplitude spectra of 0.5-1.0 mmag was achieved. This allowed the detection of pulsations in the central star of planetary nebula Abell 72, consistent with g-modes excited in GW Vir stars, and variability in RX J0122.9-7521 that could be due to pulsations, binarity or rotation. For the remaining stars from the sample that were not observed to vary, we placed upper limits for variability. After combination with literature data, our results place the fraction of pulsating PG 1159 stars within the GW Vir instability strip at 36%. An updated list of all known PG 1159 stars is provided, containing astrometric measurements from the recent Gaia DR3 data, as well as information on physical parameters, variability, and nitrogen content. Those data are used to calculate luminosities for all PG 1159 stars to place the whole sample on the theoretical Hertzsprung-Russell diagram for the first time in that way. The pulsating stars are discussed as a group, and arguments are given that the traditional separation of GW Vir pulsators in "DOV" and "PNNV" stars is misleading and should not be used.

Combining galaxy clustering information from regions of different environmental densities can help break cosmological parameter degeneracies and access non-Gaussian information from the density field that is not readily captured by the standard two-point correlation function (2PCF) analyses. However, modelling these density-dependent statistics down to the non-linear regime has so far remained challenging. We present a simulation-based model that is able to capture the cosmological dependence of the full shape of the density-split clustering (DSC) statistics down to intra-halo scales. Our models are based on neural-network emulators that are trained on high-fidelity mock galaxy catalogues within an extended-$\Lambda$CDM framework, incorporating the effects of redshift-space, Alcock-Paczynski distortions and models of the halo-galaxy connection. Our models reach sub-percent level accuracy down to $1\,h^{-1}{\rm Mpc}$ and are robust against different choices of galaxy-halo connection modelling. When combined with the galaxy 2PCF, DSC can tighten the constraints on $\omega_{\rm cdm}$, $\sigma_8$, and $n_s$ by factors of 2.9, 1.9, and 2.1, respectively, compared to a 2PCF-only analysis. DSC additionally puts strong constraints on environment-based assembly bias parameters. Our code is made publicly available on Github.

We present a clustering analysis of the BOSS DR12 CMASS galaxy sample, combining measurements of the galaxy two-point correlation function and density-split clustering down to a scale of $1\,h^{-1}{\rm Mpc}$. Our theoretical framework is based on emulators trained on high-fidelity mock galaxy catalogues that forward model the cosmological dependence of the clustering statistics within an extended-$\Lambda$CDM framework, including redshift-space and Alcock-Paczynski distortions. Our base-$\Lambda$CDM analysis finds $\omega_{\rm cdm} = 0.1201\pm 0.0022$, $\sigma_8 = 0.792\pm 0.034$, and $n_s = 0.970\pm 0.018$, corresponding to $f\sigma_8 = 0.462\pm 0.020$ at $z \approx 0.525$, which is in agreement with Planck 2018 predictions and various clustering studies in the literature. We test single-parameter extensions to base-$\Lambda$CDM, varying the running of the spectral index, the dark energy equation of state, and the density of massless relic neutrinos, finding no compelling evidence for deviations from the base model. We model the galaxy-halo connection using a halo occupation distribution framework, finding signatures of environment-based assembly bias in the data. We validate our pipeline against mock catalogues that match the clustering and selection properties of CMASS, showing that we can recover unbiased cosmological constraints even with a volume 84 times larger than the one used in this study.

We present the first wide area (2.5 x 2.5 square degrees), deep (median noise of approximately 80 microJy per beam) LOFAR High Band Antenna image at a resolution of 1.2 arcseconds by 2 arcseconds. It was generated from an 8-hour International LOFAR Telescope (ILT) observation of the ELAIS-N1 field at frequencies ranging from 120 to 168 MHz with the most up-to-date ILT imaging strategy. This intermediate resolution falls between the highest possible resolution (0.3 arcseconds) achievable by using all International LOFAR Telescope (ILT) baselines and the standard 6-arcsecond resolution in the LoTSS (LOFAR Two-meter Sky Survey) image products utilizing the LOFAR Dutch baselines only. This is the first demonstration of the feasibility of approximately 1 arcsecond imaging using the ILT, providing unique information on source morphology at scales below the surface brightness limits of higher resolutions. The total calibration and imaging time is approximately 52,000 core hours, nearly five times more than producing a 6-arcsecond image. We also present a radio source catalog containing 2263 sources detected over the 2.5 x 2.5 square degrees image of the ELAIS-N1 field, with a peak intensity threshold of 5.5 sigma. The catalog has been cross-matched with the LoTSS deep ELAIS-N1 field radio catalog, and its flux density and positional accuracy have been investigated and corrected accordingly. We find that approximately 80% of sources that we expect to be detectable based on their peak brightness in the LoTSS 6-arcsecond image are detected in this image, which is approximately a factor of two higher than for 0.3 arcsecond imaging in the Lockman Hole, implying there is a wealth of information on these intermediate scales.

We use the stellar fossil record to constrain the stellar metallicity evolution and star-formation histories of the post-starburst regions within 45 local post-starburst galaxies from the MaNGA survey. The direct measurement of the regions' stellar metallicity evolution is achieved by a new two-step metallicity model that allows for stellar metallicity to change at the peak of the starburst. We also employ a Gaussian process noise model that accounts for correlated errors introduced by the observational data reduction or inaccuracies in the models. We find that a majority of post-starburst regions (69% at $>1\sigma$ significance) increased in stellar metallicity during the recent starburst, with an average increase of 0.8 dex and a standard deviation of 0.4 dex. A much smaller fraction of PSBs are found to have remained constant (22%) or declined in metallicity (9%, average decrease 0.4 dex, standard deviation 0.3 dex). The pre-burst metallicities of the post-starburst galaxies are in good agreement with the mass-metallicity relation of local star-forming galaxies. These results are consistent with hydrodynamic simulations, which suggest that mergers between gas-rich galaxies are the primary formation mechanism of local PSBs, and rapid metal recycling during the starburst outweighs the impact of dilution by any gas inflows. The final mass-weighted metallicities of the post-starburst galaxies are consistent with the mass-metallicity relation of local passive galaxies. Our results suggest that rapid quenching following a merger-driven starburst is entirely consistent with the observed gap between the stellar mass-metallicity relations of local star-forming and passive galaxies.

KIC 8462852 is a star in the Kepler field that exhibits almost unique behaviour. The deep, irregular and aperiodic dips in its light curve have been interpreted as the breakup of a large exocomet on a highly eccentric orbit whose post-disruption material obscures the star. It is hypothesised that a nearby M-dwarf, recently confirmed to be bound to the system, could be exciting planetesimals in a source belt to high eccentricities if its orbit is highly misaligned with the belt: an effect known as the 'Eccentric Kozai-Lidov Mechanism'. To quantify how often this effect is expected to occur, this paper presents a Monte Carlo model of wide binary stars with embedded, misaligned planetesimal belts. These belts collisionally erode over time until they are excited to high eccentricities on secular timescales by a companion star if its orbit is sufficiently misaligned. The large planetesimals then produce an observable dimming signature in the light curve for a set period of time which may or may not overlap with similar events. The model finds that, for dimming events that persist for 100 yr, the most likely companion stars are located at $10^2 - 10^4$ au, the most likely belts are at $10^2-10^3$ au and the system age is most likely to be $10^2 - 10^3$ Myr. However, the probability of observing one or more stars exhibiting this phenomenon in the Kepler field is $1.3 \times 10^{-3}$, such that it is unlikely this mechanism is driving the observations of KIC 8462852.

Astrometry holds the potential for testing fundamental physics through the effects of the Stochastic Gravitational Wave Background (SGWB) in the $\sim 1-100$ nHz frequency band on precision measurements of stellar positions. Such measurements are complementary to tests made possible by the detection of the SGWB using Pulsar Timing Arrays. Here, the feasibility of using astrometry for the identification of parity-violating signals within the SGWB is investigated. This is achieved by defining and quantifying a non-vanishing $EB$ correlation function within astrometric correlation functions, and investigating how one might estimate the detectability of such signals.

It is proved, using the curved line element of a spherically symmetric charged object in general relativity and the Schwinger discharge mechanism of quantum field theory, that the orbital periods $T_{\infty}$ of test particles around central compact objects as measured by flat-space asymptotic observers are fundamentally bounded from below. The lower bound on orbital periods becomes universal (independent of the mass $M$ of the central compact object) in the dimensionless $ME_{\text{c}}\gg1$ regime, in which case it can be expressed in terms of the electric charge $e$ and the proper mass $m_{e}$ of the lightest charged particle in nature: $T_{\infty}>{{2\pi e\hbar}\over{\sqrt{G}c^2 m^2_{e}}}$ (here $E_{\text{c}}=m^2_{e}/e\hbar$ is the critical electric field for pair production). The explicit dependence of the bound on the fundamental constants of nature $\{G,c,\hbar\}$ suggests that it may reflect a fundamental physical property of the elusive quantum theory of gravity.

We consider a scenario with axions/axion-like particles Chern-Simons gravity coupling, such that gravitational waves can be produced directly from axion wave tachyonic instability in the early universe after inflation. This axion gravity term is less constrained compared to the well-searched axion photon coupling and can provide a direct and efficient production channel for gravitational waves. Such stochastic gravitational waves can be detected by either space/ground-based gravitational wave detectors or pulsar timing arrays for a broad range of axion masses and decay constants.

The complex interplay between the Solar Wind and the lunar surface serves as a quintessential example of space weathering. However, uncertainties persist regarding the influence of plasma originating from Earth's ionosphere, necessitating a comprehensive understanding of its quantitative impact. Hitherto, the dearth of reliable models has impeded accurate computation of ion flux from Earth to the Moon under varying solar wind conditions.The objective of this study is to adapt a kinetic model for the challenging conditions of having both the Earth and the Moon in a single simulation box. IAPIC, the Particle-In-Cell Electromagnetic Relativistic Global Model was modified to handle the Sun-Earth-Moon system. It employs kinetic simulation techniques that have proven invaluable tools for exploring the intricate dynamics of physical systems across various scales while minimizing the loss of crucial physics information such as backscattering. The modeling allowed to derive the shape and size of the Earth's magnetosphere and allowed tracking the O$^+$ and H$^+$ ions escaping from the ionosphere to the Moon: $\mathrm{O^+}$ tends to escape towards the dayside magnetopause, while $\mathrm{H^+}$ travels deeper into the magnetotail, extending up to the Lunar surface. In addition, plasma temperature anisotropy and backstreaming ions were simulated, allowing for future comparison with the experiment. This study shows how a kinetic model can successfully be applied to study the transport of ions in the Earth-Moon environment. A second paper will detail the effect on the Lunar environment and the impact on the Lunar water.

In the pursuit of lunar exploration and the investigation of water presence on the lunar surface, a comprehensive understanding of plasma-surface interactions is crucial since the regolith's space weathering can create H$_2$O. However, the Moon is in the Earth's magnetotail for nearly 20\% of its orbit, which could affect this water creation on the side facing the Earth if this condition shields it from the solar wind. The objective of this study is to understand how the passage of the Moon in the Earth's magnetotail affects the plasma delivery near the lunar surface. The Particle-In-Cell Electromagnetic (EM) Relativistic Global Model, known as IAPIC, is employed to kinetically simulate the Solar Wind-Magnetosphere-Ionosphere-Moon Coupling. The Earth's magnetotail does not prevent the influx of solar wind ions and ionospheric ions into the solar environment; therefore the space weathering of the regolith is not stopped in these conditions. In addition, the charge separation of solar wind ions and electrons happens is modeled, leading to electric fields and charging of the lunar surface that can be validated by observations. The study of the Sun-Earth-Moon system provides insight into the lunar environment while in the magnetotail, which is essential to better interpret the results of future Lunar missions. It also provides insights in the Lunar charging in different conditions that could affect the human presence on the Moon.

AI accelerators have proliferated in data centers in recent years and are now almost ubiquitous. In addition, their computational power and, most importantly, their energy efficiency are up to orders of magnitude higher than that of traditional computing. Over the last years, various methods and optimizations have been tested to use these hybrid systems for simulations in the context of astroparticle physics utilizing CORSIKA. The main focus of this talk is the propagation of optical, i.e. fluorescence and Cherenkov, photons through low density inhomogeneous media in the context of the next generation CORSIKA8 simulation framework. Different techniques used and approximations, e.g. the atmospheric model, tested during the development will be presented. The trade-off between performance and precision allows the experiment to achieve its physical precision limited to the real resolution of the experiment and not invest power and time in vanishing precision gains. The additional comparison of classical CPU-based simulations with the new methods validates these methods and allows evaluation against a known baseline.

Tidal flat topography provides crucial insights for understanding tidal flats and their dynamic evolution. However, the wide-ranging and rapidly changing nature of tidal flats, which are periodically submerged in shallow water, pose challenges for many current monitoring methods in terms of both efficiency and precision. In this study, we considered the dynamic process of tidal flat submergence and utilized time-series Sentinel-2 images on Google Earth Engine (GEE) to calculate the tidal flat exposure frequency. This information was used to determine the spatial extent of the tidal flats, and subsequently, by employing ICESat-2 data, we established a 1D-linear regression model based on elevation and frequency values, which realizes the inversion of the tidal flat elevation within Cixi City. The study shows the following: (1) the tidal flat exposure frequency and ICESat-2 elevation data exhibit a strong positive correlation (R2=0.85); (2) the tidal flat area within Cixi City is 115.81 km2, and the overall accuracy is 95.36%; and (3) the elevation range of the tidal flats in the study area is between -0.42 and 2.73 m, and the mean absolute error (MAE) is 0.24 m. Additionally, we consider that the temporal resolution of remote sensing imagery plays a crucial role in determining the accuracy of the elevation inversion, and we found that higher tidal flats exhibit better inversion accuracy than lower tidal flats.

We perform an updated analysis on a long-lived axion domain wall (DW) network. By simulating the axion field on a 3D lattice and fitting an analytical model for the DW evolution, we identify the leading energy loss mechanisms of the DWs and compute the spectrum of axions emitted from the network. The contribution from the DWs to axion dark matter (DM) density is derived, with viable parameter space given. The application to both QCD axions and general axion-like particles (ALPs) are considered. Due to the new approaches taken, while our results bear consistency with earlier literature, notable discrepancies are also revealed, such as the prediction for DM abundance, which may have a profound impact on axion phenomenology at large.

Primordial black holes (PBHs) are plausible dark matter candidates that formed from the gravitational collapse of primordial density fluctuations. Current observational constraints allow asteroid-mass PBHs to account for all of the cosmological dark matter. We show that elastic scattering of a cosmological gravitational wave background, these black holes generate spectral distortions on the background of 0.3% for cosmologically relevant frequencies. Scattering from stellar objects induce much smaller distortions. Detectability of this signal depends on our ultimate understanding of the unperturbed background spectrum.

Recent years have witnessed a rise in interest in the geometrical trinity of General Relativity and its extensions. This interest has been fuelled by novel insights into the nature of gravity, the possibility to address computational and conceptual questions -- such as the determination of black hole entropy or the definition of gravitational energy-momentum -- from a new perspective. In particular, $f(Q)$ gravity has also inspired numerous works on black holes, wormholes, and cosmology. In the latter case, $f(Q)$ models have the potential to elucidate phenomena in both early and late-time cosmology without necessitating the inclusion of dark energy, the inflaton field, or dark matter. Particularly noteworthy is the role of $f(Q)$ theories in addressing cosmological tensions, presenting exciting possibilities for reshaping our understanding of gravity and its manifestations in cosmology. The emergence of intriguing new black hole solutions and the potential existence of wormhole solutions suggest the presence of novel physics within the realm of strong gravity. These phenomena have become increasingly measurable only in recent times, opening up exciting avenues for further exploration and discovery. This review is tailored to students and researchers alike. It offers a self-contained and pedagogical introduction to metric-affine geometry--The mathematical foundation and indispensable tool upon which the geometrical trinity of General Relativity as well as its various extensions are built.

A cosmological model based on holographic scenarios is formulated in a flat Friedmann-Robertson-Walker universe. To formulate this model, the cosmological horizon is assumed to have a general entropy and a general temperature (including Bekenstein-Hawking entropy and Gibbons-Hawking temperature, respectively). In addition, a holographic-like connection [Eur. Phys. J. C 83, 690 (2023) (arXiv:2212.05822)] and Padmanabhan's holographic equipartition law are assumed for the entropy and temperature, and the Friedmann and acceleration equations are derived from these. The derived Friedmann and acceleration equations include both the entropy and the temperature and are slightly complicated, but can be combined into a single simple equation, corresponding to a similar equation that describes the background evolution of the universe in time-varying $\Lambda (t)$ cosmologies. The simple equation depends on the entropy but not on the temperature because the temperatures in the Friedmann and acceleration equations cancel each other. These results imply that the holographic-like connection should be consistent with Padmanabhan's holographic equipartition law through the present model and that the entropy plays a more important role. When the Gibbons-Hawking temperature is used as the temperature, the Friedmann and acceleration equations are found to be equivalent to those for a $\Lambda(t)$ model. A particular case of the present model is also examined, applying a power-law corrected entropy.

Motivated by the recent heated debate on whether the masses of local objects, such as compact stars or black holes (BHs), may be affected by the large-scale, cosmological dynamics, we analyze the conditions under which, in a general relativity framework, such a coupling small/large scales is allowed. We shed light on some controversial arguments, which have been used to rule out the latter possibility. We argue that the actual observational quantity at play is the quasi-local Misner-Sharp mass (MS), and we find that the cosmological coupling occurs whenever the energy of the central object is quantified by it. Conversely, the decoupling occurs whenever the MS mass is fully equivalent to the (nonlocal) Arnowitt-Deser-Misner (ADM) mass. Consequently, for singular BHs embedded in cosmological backgrounds, like the Schwarzschild-de Sitter or McVittie solutions, we show that there is no cosmological coupling, confirming previous results in the literature. Furthermore, we show that nonsingular compact objects couple to the cosmological background, as quantified by their MS mass. We conclude that observational evidence of cosmological coupling of astrophysical BHs would be the smoking gun of their nonsingular nature.

Scalar fields interacting with the primordial curvature perturbation during inflation may communicate their statistics to the latter. This situation motivates the study of how the probability density function (PDF) of a light spectator field $\varphi$ in a pure de Sitter space-time, becomes non-Gaussian under the influence of a scalar potential ${\mathcal V(\varphi)}$. One approach to this problem is offered by the stochastic formalism introduced by Starobinsky and Yokoyama. It results in a Fokker-Planck equation for the time-dependent PDF $\rho (\varphi , t)$ describing the statistics of $\varphi$ which, in the limit of equilibrium gives one back the solution $\rho (\varphi) \propto \exp \big[ - \frac{8 \pi^2}{3 H^4} {\mathcal V(\varphi)} \big]$. We study the derivation of $\rho (\varphi , t)$ using quantum field theory tools. Our approach yields an almost Gaussian distribution function, distorted by minor corrections comprised of terms proportional to powers of $\Delta N \times \mathcal O(\partial_\varphi) {\mathcal V(\varphi)}$, where $\Delta N$ is the number of $e$-folds succeeding the Hubble-horizon crossing of $\varphi$'s wavelengths, and $\mathcal O(\partial_\varphi)$ stands for a derivative operator acting on ${\mathcal V(\varphi)}$. This general form is obtained perturbatively and remains valid even with loop corrections. Our solution satisfies a Fokker-Planck equation that receives corrections with respect to the one found within the stochastic approach, allowing us to comment on the validity of the standard equilibrium solution for generic potentials. We posit that higher order corrections to the Fokker-Planck equation may become important towards the equilibrium.

Experiments aimed at detecting ultralight dark matter typically rely on resonant effects, which are sensitive to the dark matter mass that matches the resonance frequency. In this study, we investigate the nucleon couplings of ultralight axion dark matter using a magnetometer operating in a nuclear magnetic resonance (NMR) mode. Our approach involves the use of a $^{21}$Ne spin-based sensor, which features the lowest nuclear magnetic moment among noble-gas spins. This configuration allows us to achieve an ultrahigh sensitivity of 0.73 fT/Hz$^{1/2}$ at around 5 Hz, corresponding to energy resolution of approximately 1.5$\times 10^{-23}\,\rm{eV/Hz^{1/2}}$. Our analysis reveals that under certain conditions it is beneficial to scan the frequency with steps significantly larger than the resonance width. The analytical results are in agreement with experimental data and the scan strategy is potentially applicable to other resonant searches. Further, our study establishes stringent constraints on axion-like particles (ALP) in the 4.5--15.5 Hz Compton-frequency range coupling to neutrons and protons, improving on prior work by several-fold. Within a band around 4.6--6.6 Hz and around 7.5 Hz, our laboratory findings surpass astrophysical limits derived from neutron-star cooling. Hence, we demonstrate an accelerated resonance search for ultralight dark matter, achieving an approximately 30-fold increase in scanning step while maintaining competitive sensitivity.

We study the origin of the real intermediate state subtraction problem and compare its different solutions. We show that the ambiguity in subtraction schemes arises from the on-shell approximation for the 2-point functions that reduces the Schwinger-Dyson equations to the Boltzmann limit. We also suggest a new subtraction scheme which, unlike the earlier definitions, never leads to negative scattering rates. This scheme also quantifies the validity of the on-shell limit in terms of an effective one-particle weight function $R(\Delta )$, where $\Delta$ measures the region around the resonance associated with the real state.

The Higgs-portal scalar dark matter (DM) model is a simple extension of the Standard Model (SM) to incorporate a DM particle to the SM, where a $Z_2$-odd real scalar field is introduced as a DM candidate. We consider this DM model in the context of 5-dimensional brane-world cosmology, where our 3-dimensional space is realized as a hyper-surface embedded in 4-dimensional space. In the setup, all the SM and DM fields reside on the hyper-surface while graviton lives in the bulk. We consider two well-known brane-world cosmologies, namely, the Randall-Sundrum (RS) and the Gauss-Bonnet (GB) brane-world cosmologies, in which the standard Big Bang cosmology is reproduced at low temperatures below the so-called ``transition temperature" while at high temperatures the expansion law of the universe is significantly modified. Such a non-standard expansion law directly impacts the prediction for the relic density of the Higgs-portal DM. We investigate the brane-world cosmological effects and identify the allowed model parameter region by combining the constraints from the observed DM relic density, and the direct and indirect DM detection experiments. It is well-known that only DM masses in the vicinity of half the Higgs boson mass are allowed in the Higgs-portal scalar DM model. We find that the allowed parameter region becomes more severely constrained and even disappears in the RS cosmology, while the GB cosmological effect significantly enlarges the allowed region. Upon discovering Higgs-portal DM, we can determine transition temperature in the GB brane-world cosmology.