Papers for Wednesday, May 24 2023

Solid states of strange-cluster matter called strangeon matter can form strangeon stars that are highly compact. We show that strangeon matter and strangeon stars can be recast into dimensionless forms by a simple reparametrization and rescaling, through which we manage to maximally reduce the number of degrees of freedom. With this dimensionless scheme, we find that strangeon stars are generally compact enough to feature a photon sphere that is essential to foster gravitational-wave (GW) echoes. Rescaling the dimension back, we illustrate its implications on the expanded dimensional parameter space, and calculate the GW echo frequencies associated with strangeon stars, showing that the minimum echo frequency is $\sim 8$ kHz for empirical parameter space that satisfies the GW170817 constraint, and can reduce to $\mathcal O(100)$ Hertz at the extended limit.

We compare the evolution of binary systems evolved in the MESA stellar evolution code to those in the COSMIC population synthesis code. Our aim is to convey the robustness of the equations that model binary evolution in the COSMIC code, particularly for the cases of high mass stars with closely orbiting compact object companions. Our larger goal is to accurately model the rates of these systems, as they are promising candidates for the progenitor systems behind energetic, longer lasting, radio bright GRB jets. These systems also may be key contributors to the rates of binary black hole mergers throughout our universe.

It is well known that the first structures that form from small fluctuations in a self-gravitating, collisionless and initially smooth cold dark matter (CDM) fluid are pancakes. We study the gravitational force generated by such pancakes just after shell-crossing, and find a simple analytical formula for the force along the collapse direction, which can be applied to both the single- and multi-stream regimes. The formula is tested on the early growth of CDM protohaloes seeded by two or three crossed sine waves. Adopting the high-order Lagrangian perturbation theory (LPT) solution as a proxy for the dynamics, we confirm that our analytical prediction agrees well with the exact solution computed by direct resolution of the Poisson equation, as long as the caustic structure remains locally sufficiently one-dimensional. These results are further confirmed by comparisons of the LPT predictions performed this way to measurements in Vlasov simulations performed with the public code ColDICE. We also show that the component of the force orthogonal to the collapse direction preserves its single stream nature by not changing qualitatively before and after the collapse, allowing sufficiently high-order LPT acceleration to be used to approximate it accurately as long as the LPT series converges. As expected, solving Poisson equation on the density field generated with LPT displacement provides a more accurate force than the LPT acceleration itself, as a direct consequence of the faster convergence of the LPT series for the positions than for the accelerations. This may provide a clue on improving standard LPT predictions. Our investigations represent a very needed first step to study analytically gravitational dynamics in the multi-stream regime, by estimating, at leading order in time and space the proper backreaction on the gravitational field inside the pancakes.

We present a novel, data-driven analysis of Galactic dynamics, using unsupervised machine learning -- in the form of density estimation with normalizing flows -- to learn the underlying phase space distribution of 6 million nearby stars from the Gaia DR3 catalog. Solving the collisionless Boltzmann equation with the assumption of approximate equilibrium, we calculate -- for the first time ever -- a model-free, unbinned, fully 3D map of the local acceleration and mass density fields within a 3 kpc sphere around the Sun. As our approach makes no assumptions about symmetries, we can test for signs of disequilibrium in our results. We find our results are consistent with equilibrium at the 10% level, limited by the current precision of the normalizing flows. After subtracting the known contribution of stars and gas from the calculated mass density, we find clear evidence for dark matter throughout the analyzed volume. Assuming spherical symmetry and averaging mass density measurements, we find a local dark matter density of $0.47\pm 0.05\;\mathrm{GeV/cm}^3$. We fit our results to a generalized NFW, and find a profile broadly consistent with other recent analyses.

We present the lifetime star formation histories (SFHs) for six ultra-faint dwarf (UFD; $M_V>-7.0$, $ 4.9<\log_{10}({M_*(z=0)}/{M_{\odot}})<5.5$) satellite galaxies of M31 based on deep color-magnitude diagrams constructed from Hubble Space Telescope imaging. These are the first SFHs obtained from the oldest main sequence turn-off of UFDs outside the halo of the Milky Way (MW). We find that five UFDs formed at least 50% of their stellar mass by $z=5$ (12.6 Gyr ago), similar to known UFDs around the MW, but that 10-40% of their stellar mass formed at later times. We uncover one remarkable UFD, And XIII, which formed only 10% of its stellar mass by $z=5$, and 75% in a rapid burst at $z\sim2-3$, a result that is robust to choices of underlying stellar model and is consistent with its predominantly red horizontal branch. This ''young'' UFD is the first of its kind and indicates that not all UFDs are necessarily quenched by reionization, which is consistent with predictions from several cosmological simulations of faint dwarf galaxies. SFHs of the combined MW and M31 samples suggest reionization did not homogeneously quench UFDs. We find that the least massive MW UFDs ($M_*(z=5) \lesssim 5\cdot10^4 M_{\odot}$) are likely quenched by reionization, whereas more massive M31 UFDs ($M_*(z=5) \gtrsim 10^5 M_{\odot}$) may only have their star formation suppressed by reionization and quench at a later time. We discuss these findings in the context of the evolution and quenching of UFDs.

We investigate the production efficiency of ionizing photons ($\xi_{ion}^*$) of 1174 galaxies with secure redshift at z=2-5 from the VANDELS survey to determine the relation between ionizing emission and physical properties of bright and massive sources. We constrain $\xi_{ion}^*$ and galaxy physical parameters by means of spectro-photometric fits performed with the BEAGLE code. The analysis exploits the multi-band photometry in the VANDELS fields, and the measurement of UV rest-frame emission lines (CIII]$\lambda 1909$, HeII$\lambda 1640$, OIII]$\lambda 1666$) from deep VIMOS spectra. We find no clear evolution of $\xi_{ion}^*$ with redshift within the probed range. The ionizing efficiency slightly increases at fainter $M_{UV}$, and bluer UV slopes, but these trends are less evident when restricting the analysis to a complete subsample at log(M$_{star}$/M$_{\odot}$)$>$9.5. We find a significant trend of increasing $\xi_{ion}^*$ with increasing EW(Ly$\alpha$), with an average log($\xi_{ion}^*$/Hz erg$^{-1}$)$>$25 at EW$>$50\AA, and a higher ionizing efficiency for high-EW CIII]$\lambda 1909$ and OIII]$\lambda 1666$ emitters. The most significant correlations are found with respect to stellar mass, specific star-formation rate (sSFR) and SFR surface density ($\Sigma_{SFR}$). The relation between $\xi_{ion}^*$ and sSFR shows a monotonic increase from log($\xi_{ion}^*$/Hz erg$^{-1}$) $\sim$24.5 at log(sSFR)$\sim$-9.5$yr^{-1}$ to $\sim$25.5 at log(sSFR)$\sim$-7.5$yr^{-1}$, a low scatter and little dependence on mass. The objects above the main-sequence of star-formation consistently have higher-than-average $\xi_{ion}^*$. A clear increase of $\xi_{ion}^*$ with $\Sigma_{SFR}$ is also found, with log($\xi_{ion}^*$/Hz erg$^{-1}$)$>$25 for objects at $\Sigma_{SFR}>$10 M$_{\odot}/yr/kpc^2$.(Abridged)

We define a physically-motivated measure for galactic bar length, called the dynamical length. The dynamical length of the bar corresponds to the radial extent of the orbits that are the backbone supporting the bar feature. We propose a direct observational technique using integral field unit spectroscopy to measure it. Identifying these orbits and using the dynamical length is a more faithful tracer of the secular evolution and influence of the bar. We demonstrate the success of the metric for recovering the maximal bar-parenting orbit in a range of simulations, and to show its promise we perform its measurement on a real galaxy. We also study the difference between traditionally used ellipse fit approaches to determine bar length and the dynamical length proposed here in a wide range of bar-forming N-body simulations of a stellar disc and dark matter halo. We find that ellipse fitting may severely overestimate measurements of the bar length by a factor of 1.5-2.5 relative to the extent of the orbits that are trapped and actually comprise the bar. This bias leads to overestimates of both bar mass and the ratio of corotation radius to bar length, i.e. the bar speed, affecting inferences about the evolution of bars in the real universe.

The 500ks Chandra ACIS-I observation of the field around the $z=6.31$ quasar SDSS J1030+0524 is currently the 5th deepest extragalactic X-ray survey. The rich multi-band coverage of the field allowed for an effective identification and redshift determination of the X-ray source counterparts: to date a catalog of 243 extragalactic X-ray sources with either a spectroscopic or photometric redshift estimate in the range $z\approx0-6$ is available over a 355 arcmin$^2$ area. Given its depth and the multi-band information, this catalog is an excellent resource to investigate X-ray spectral properties of distant Active Galactic Nuclei (AGN) and derive the redshift evolution of their obscuration. We performed a thorough X-ray spectral analysis for each object in the sample, measuring its nuclear column density $N_{\rm H}$ and intrinsic (de-absorbed) 2-10 keV rest-frame luminosity, $L_{2-10}$. Whenever possible, we also used the presence of the Fe K$_\alpha$ emission line to improve the photometric redshift estimates. We measured the fractions of AGN hidden by column densities in excess of $10^{22}$ and $10^{23}$cm$^{-2}$ ($f_{22}$ and $f_{23}$, respectively) as a function of $L_{2-10}$ and redshift, and corrected for selection effects to recover the intrinsic obscured fractions. At $z\sim 1.2$, we found $f_{22}\sim0.7-0.8$ and $f_{23}\sim0.5-0.6$, respectively, in broad agreement with the results from other X-ray surveys. No significant variations with X-ray luminosity were found within the limited luminosity range probed by our sample (log$L_{2-10}\sim 42.8-44.3$). When focusing on luminous AGN with log$L_{2-10}\sim44$ to maximize the sample completeness up to large cosmological distances, we did not observe any significant change in $f_{22}$ or $f_{23}$ over the redshift range $z\sim0.8-3$. Nonetheless, the obscured fractions we measure are significantly higher than ...

SDSS-V is carrying out a dedicated survey for white dwarfs, single and in binaries, and we report the analysis of the spectroscopy of cataclysmic variables (CVs) and CV candidates obtained during the final plug plate observations of SDSS. We identify eight new CVs, spectroscopically confirm 53 and refute eleven published CV candidates, and we report 21 new or improved orbital periods. Combined with previously published data, the orbital period distribution of the SDSS-V CVs does not clearly exhibit a period gap. This is consistent with previous findings that spectroscopically identified CVs have a larger proportion of short-period systems compared to samples identified from photometric variability. Remarkably, despite a systematic search, we find very few period bouncers. We estimate the space density of period bouncers to be $\simeq0.2\times10^{-6}\,\mathrm{pc}^{-3}$, i.e. they represent only a few per cent of the total CV population. This suggests that during their final phase of evolution, CVs either destroy the donor, e.g. via a merger, or that they become detached and cease mass transfer.

Chemical abundances are an essential tool in untangling the Milky Way's enrichment history. However, the evolution of the interstellar medium abundance gradient with cosmic time is lost as a result of radial mixing processes. For the first time, we quantify the evolution of many observational abundances across the Galactic disk as a function of lookback time and birth radius, $R_\text{birth}$. Using an empirical approach, we derive $R_\text{birth}$ estimates for 145,447 APOGEE DR17 red giant disk stars, based solely on their ages and [Fe/H]. We explore the detailed evolution of 6 abundances (Mg, Ca ($\alpha$), Mn (iron-peak), Al, C (light), Ce (s-process)) across the Milky Way disk using 87,426 APOGEE DR17 red giant stars. We discover that the interstellar medium had three fluctuations in the metallicity gradient $\sim 9$, $\sim 6$, and $\sim4$ Gyr ago. The first coincides with the end of high-$\alpha$ sequence formation around the time of the Gaia-Sausage-Enceladus disruption, while the others are likely related to passages of the Sagittarius dwarf galaxy. A clear distinction is found between present-day observed radial gradients with age and the evolution with lookback time for both [X/Fe] and [X/H], resulting from the significant flattening and inversion in old populations due to radial migration. We find the [Fe/H]--[$\alpha$/Fe] bimodality is also seen as a separation in the $R_\text{birth}$--[X/Fe] plane for the light and $\alpha$-elements. Our results recover the chemical enrichment of the Galactic disk over the past 12 Gyr, providing tight constraints on Galactic disk chemical evolution models.

Numerical simulations have become one of the key tools used by theorists in all the fields of astrophysics and cosmology. The development of modern tools that target the largest existing computing systems and exploit state-of-the-art numerical methods and algorithms is thus crucial. In this paper, we introduce the fully open-source highly-parallel, versatile, and modular coupled hydrodynamics, gravity, cosmology, and galaxy-formation code Swift. The software package exploits hybrid task-based parallelism, asynchronous communications, and domain-decomposition algorithms based on balancing the workload, rather than the data, to efficiently exploit modern high-performance computing cluster architectures. Gravity is solved for using a fast-multipole-method, optionally coupled to a particle mesh solver in Fourier space to handle periodic volumes. For gas evolution, multiple modern flavours of Smoothed Particle Hydrodynamics are implemented. Swift also evolves neutrinos using a state-of-the-art particle-based method. Two complementary networks of sub-grid models for galaxy formation as well as extensions to simulate planetary physics are also released as part of the code. An extensive set of output options, including snapshots, light-cones, power spectra, and a coupling to structure finders are also included. We describe the overall code architecture, summarize the consistency and accuracy tests that were performed, and demonstrate the excellent weak-scaling performance of the code using a representative cosmological hydrodynamical problem with $\approx$$300$ billion particles. The code is released to the community alongside extensive documentation for both users and developers, a large selection of example test problems, and a suite of tools to aid in the analysis of large simulations run with Swift.

The seed of dark matter can be generated from light spectator fields during inflation through a similar mechanism that the seed of observed large scale structures are produced from the inflaton field. The accumulated energy density of the corresponding excited modes, which is subdominant during inflation, dominates energy density of the universe later around the time of matter and radiation equality and plays the role of dark matter. For spin-2 spectator fields, Higuchi bound may seem to prevent excitation of such light modes since deviation of the inflationary background from the exact de Sitter spacetime is very small. However, sizable interactions with the inflaton field breaks (part of) isometries of the de Sitter space in the inflationary background and relaxes the Higuchi bound. Looking for this possibility in the context of effective field theory of inflation, we suggest a dark matter model consisting of spin-2 particles that produce during inflation.

We revise primordial black holes (PBHs) production in the axion-curvaton model, in light of recent developments in the computation of their abundance accounting for non-gaussianities (NGs) in the curvature perturbation up to all orders. We find that NGs intrinsically generated in such scenarios have a relevant impact on the phenomenology associated to PBHs and, in particular, on the relation between the abundance and the signal of second-order gravitational waves. We show that this model could explain both the totality of dark matter in the asteroid mass range and the tentative signal reported by the NANOGrav and IPTA collaborations in the nano-Hz frequency range. En route, we provide a new, explicit computation of the power spectrum of curvature perturbations going beyond the sudden-decay approximation.

Kinematic misalignment between gas and stellar components observed in a certain fraction of galaxies. It believed to be caused by acquisition of gas from the external reservoir by major or minor mergers, accretion from cosmological filaments or circumgalactic medium, etc. We aim to constrain possible sources of the gas that forms counter-rotating component. We derived the gas-phase oxygen abundance in 69 galaxies with kinematic misalignment between gas and stellar components from MaNGA DR17 survey and compared it with the metallicity expected according to the mass-metallicity relation. We found that the oxygen abundance of the counter-rotating gas in our sample is higher than 8.2 dex that excludes significant role of inflow of pristine gas. Meanwhile, there is a significant difference in the oxygen abundance of the counter-rotating gas between red and blue galaxies. In general, the oxygen abundance is lower than expected for their stellar mass in red galaxies, but is compatible with or even higher than typical values for their stellar mass in blue galaxies. We showed that the exchange of enriched gas between galaxies is the most plausible mechanism for explaining the metallicity of counter-rotating gas components in galaxies of all masses and colors. Meanwhile, minor mergers may play a significant role in the formation of counter-rotating gas components in red and quenched galaxies.

The mechanisms responsible for generating spin-orbit misalignments in exoplanetary systems are still not fully understood. It is unclear whether these misalignments are related to the migration of hot Jupiters or are a consequence of general star and planet formation processes. One promising method to address this question is to constrain the distribution of spin-orbit angle measurements for a broader range of planets beyond hot Jupiters. In this work, we present the sky-projected obliquity ($\lambda=-68.1_{-14.7}^{+21.2} \,^{\circ}$) for the warm sub-Saturn TOI-1842b, obtained through a measurement of the Rossiter-McLaughlin effect using WIYN/NEID. Using the projected obliquity, the stellar rotation period obtained from the TESS light curve, and the projected rotation velocity from spectral analysis, we infer the 3D spin-orbit angle ($\psi$) to be $\psi=73.3^{+16.3}_{-12.9} \,^{\circ}$. As the first spin-orbit angle determination made for a sub-Saturn-mass planet around a massive ($M_{\rm *}=1.45 \,{\rm M_\odot}$) star, our result presents an opportunity to examine the orbital geometries for new regimes of planetary systems. When combined with archival measurements, our observations of TOI-1842b support the hypothesis that the previously established prevalence of misaligned systems around hot, massive stars may be driven by planet-planet dynamical interactions. In massive stellar systems, multiple gas giants are more likely to form and can then dynamically interact with each other to excite spin-orbit misalignments.

Extreme, young stellar populations are considered the primary contributor to cosmic reionization. However, how Lyman-continuum (LyC) escapes these galaxies remains highly elusive because LyC escape can vary on sub-galactic scales that are technically challenging to observe in LyC emitters. We investigate the Sunburst Arc: a strongly lensed, LyC emitter at $z=2.37$. This galaxy reveals the exceptionally small scale (tens of parsecs) physics of LyC escape thanks to high magnification from strong lensing. Analyzing HST broadband and narrowband imaging, we find that the small ($<$100 pc) LyC leaking region shows distinctly extreme properties: a very blue UV slope ($\beta=-2.9\pm0.1$), high ionization state ([OIII]$\lambda 5007$/[OII]$\lambda 3727=11\pm3$ and [OIII]$\lambda 5007$/H$\beta=6.8\pm0.4$), strong oxygen emission (EW([OIII])$=1095\pm 40 \r{A}$), and high Lyman-$\alpha$ escape fraction ($0.3\pm 0.03$), none of which are found in any non-leaking regions of the galaxy. Moreover, a UV slope comparison with starburst population models indicates that the leaking region's UV emission consists of nearly ``pure'' stellar light with minimal contamination from surrounding nebular continuum emission and dust extinction. These results suggest a highly directional LyC escape such that LyC is produced and escapes from a small, extreme starburst region where the stellar feedback from an ionizing star cluster may create an anisotropic ``pencil beam'' viewing geometry in the surrounding gas. As a result, unabsorbed LyC directly escapes through these perforated hole(s). Importantly, such anisotropic escape processes imply that unfavorable sightline effects are a crucial contributor to the significant scatters between galaxy properties and LyC escape fraction in observations and that strong lensing uniquely reveals the small-scale physics that regulates the ionizing budget of galaxies for reionization.

Recent advances in cosmic-ray detectors have provided exceptional sensitivities of dark matter with high angular resolution. Motivated by this, we present a comprehensive study of cosmic-ray flux from dark matter decay in dwarf spheroidal galaxies (dSphs), with a focus on detectors possessing arcsecond-level field of view and/or angular resolution. We propose to use differential $D$-factors, which are estimated for various dSphs since such detectors are sensitive to their dark matter distributions. Our findings reveal that the resulting signal flux can experience a more than $O$(1-10) enhancement with different theoretical uncertainty compared to traditional estimations. Based on this analysis, we find that the Infrared Camera and Spectrograph (IRCS) installed on the 8.2m Subaru telescope can be a good dark matter detector for the mass in the eV range, particularly axion-like particles (ALPs). Observing the Draco or Ursa Major II galaxies with the IRCS for just a few nights will be sufficient to surpass the stellar cooling bounds for ALP dark matter with a mass in the range of $1\,{\rm eV} \lesssim m_a \lesssim 2\,\rm eV$.

A variety of energy sources, ranging from dynamic processes like magnetic reconnection and waves to quasi-steady terms like the plasma pressure, may contribute to the acceleration of the solar wind. We utilize a combination of charged particle and magnetic field observations from the Parker Solar Probe (PSP) to attempt to quantify the steady-state contribution of the proton pressure, the electric potential, and the wave energy to the solar wind proton acceleration observed by PSP between 13.3 and ~100 solar radii (RS). The proton pressure provides a natural kinematic driver of the outflow. The ambipolar electric potential acts to couple the electron pressure to the protons, providing another definite proton acceleration term. Fluctuations and waves, while inherently dynamic, can act as an additional effective steady-state pressure term. To analyze the contributions of these terms, we utilize radial binning of single-point PSP measurements, as well as repeated crossings of the same stream at different distances on individual PSP orbits (i.e. "fast radial scans"). In agreement with previous work, we find that the electric potential contains sufficient energy to fully explain the acceleration of the slower wind streams. On the other hand, we find that the wave pressure plays an increasingly important role in the faster wind streams. The combination of these terms can explain the continuing acceleration of both slow and fast wind streams beyond 13.3 RS.

Black holes can be produced in collapse of small-scale dark matter structures, which can happen at any time from the early to present-day universe. Microstructure black holes (MSBHs) can have a wide range of masses. Small MSBHs evaporate via Hawking radiation with lifetimes shorter than the age of the universe, but they are not subject to the usual early-universe bounds on the abundance of small primordial black holes. We investigate the possible signal of such a population of exploding, late-forming black holes, constraining their abundance with observations from diffuse extragalactic gamma- and x-ray sources, the galactic center, and dwarf spheroidal galaxies.

Using a semi-analytical model of the evolution of the Milky Way, we show how secular evolution can create distinct overdensities in the phase space of various properties (e.g. age vs metallicity or abundance ratios vs age) corresponding to the thin and thick discs. In particular, we show how key properties of the Solar vicinity can be obtained by secular evolution, with no need for external or special events, like galaxy mergers or paucity in star formation. This concerns the long established double-branch behaviour of [alpha/Fe] vs metallicity and the recently found non-monotonic evolution of the stellar abundance gradient, evaluated at the birth radii of stars. We extend the discussion to other abundance ratios and we suggest a classification scheme, based on the nature of the corresponding yields (primary vs secondary or odd elements) and on the lifetimes of their sources (short-lived vs long-lived ones). The latter property is critical in determining the single- or double- branch behavior of an elementary abundance ratio in the Solar neighborhood. We underline the high diagnostic potential of this finding, which can help to separate clearly elements with sources evolving on different timescales and help determining the site of e.g. the r-process(es). We define the "abundance distance" between the thin and thick disc sequences as an important element for such a separation. We also show how the inside-out evolution of the Milky Way disc leads rather to a single-branch behavior in other disc regions.

Black holes may form in present-day collapse of microscopic structures of dark matter. We show that, if microstructure black holes (MSBH) with mass $m\sim 10^{13}~g$ are produced, the spectrum of gamma rays from their evaporation agrees remarkably well with the GeV excess observed by Fermi Gamma-ray Space Telescope, while still avoiding all observational constraints. We also discuss the generic requirements for MSBHs to explain the GeV excess.

Star-forming galaxies populate a main sequence (MS), a well-defined relation between stellar mass (M*) and star-formation rate (SFR). Starburst (SB) galaxies lie significantly above the relation whereas quenched galaxies lie below the sequence. In order to study the evolution of galaxies on the SFR-M* plane and its connection to the gas content, we use the fact that recent episodes of star formation can be pinpointed by the existence of gamma-ray bursts (GRBs). Here we present sensitive [CI]-nondetections of z$\sim$2 ultra luminous infrared (ULIRG) GRB host galaxies. We find that our GRB hosts have similar molecular masses to those of other ULIRGs. However, unlike other ULIRGs, the GRB hosts are located at the MS or only a factor of a few above it. Hence, our GRB hosts are caught in the transition toward the SB phase. This is further supported by the estimated depletion times, which are similar to those of other transitioning galaxies. The GRB hosts are [CI]-dark galaxies, defined as having a [CI]/CO temperature brightness ratio of <0.1. Such a low [CI]/CO ratio has been found in high-density environments (nH > 10$^4$ cm$^{-3}$) where CO is shielded from photodissociation, leading to under-abundances of [CI]. This is consistent with the merger process that is indeed suggested for our GRB hosts by their morphologies.

Several asteroids are known to be shaped like toy tops. This paper models Top-Shaped Asteroids (TSAs) as Homogeneous Symmetric Lenses (HSLs), and derives their rotational, self-gravitational, and total energies as functions of their mass, density, and angular momentum. Then we raise, test, and ultimately reject the hypothesis that TSAs take the shape of lowest total energy, subject to the constraint that they keep the same mass, density, and angular momentum, while remaining HSLs. Other processes must control the shapes of TSAs. For completeness, we also describe a Core-Mantle Model for TSAs, as well as an Inverted Core-Mantle Model, and derive their self-gravitational energies, along with their rotational energies. The gravitational potential at the center of an HSL then is derived.

Thin galactic discs and nuclear stellar discs (NSDs) are fragile structures that can be easily disturbed by merger events. By studying the age of the stellar populations in present-day discs, we can learn about the assembly history of galaxies and place constraints on their past merger events. Following on the steps of our initial work, we explore the fragility of such disc structures in intermediate-mass-ratio dry encounters using the previously constructed $N$-body model of the Fornax galaxy NGC 1381 (FCC 170), which hosts both a thin galactic disc and a NSD. We dismiss major and minor encounters, as the former were previously shown to easily destroy thin-disc structures, whereas the latter take several Hubble times to complete in the specific case of FCC 170. The kinematics and structure of the thin galactic disc are dramatically altered by the mergers, whereas the NSD shows a remarkable resilience, exhibiting only a smooth increase of its size when compared to the model evolved in isolation. Our results suggest that thin galactic discs are better tracers for intermediate-mass-ratio mergers, while NSDs may be more useful for major encounters. Based on our simulations and previous analysis of the stellar populations, we concluded that FCC 170 has not experienced any intermediate-mass-ratio dry encounters for at least $\sim$10 Gyr, as indicated by the age of its thin-disc stellar populations.

Here we study the possible improvements of the existing constraints on the upper bound of graviton mass by the analysis of the stellar orbits around the supermassive black hole (SMBH) at the Galactic Center (GC) in the framework of Yukawa gravity. Main motivation for this study is recent detection of Schwarzschild precession in the orbit of S2 star around the SMBH at the GC by the GRAVITY Collaboration in 2020. They indicated that the orbital precession of the S2 star is close to the General Relativity (GR) prediction, but with possible small deviation from it, and parametrized this effect by introducing an ad hoc factor in the parameterized post-Newtonian (PPN) equations of motion. Here we use the value of this factor presented by GRAVITY in order to perform two-body simulations of the stellar orbits in massive gravity using equations of motion in the modified PPN formalism, as well as to constrain the range of massive interaction $\Lambda$. From the obtained values of $\Lambda$, and assuming that it corresponds to the Compton wavelength of graviton, we then calculated new estimates for the upper bound of graviton mass which are found to be independent, but consistent with the LIGO's estimate of graviton mass from the first gravitational wave (GW) signal GW150914. We also performed Markov chain Monte Carlo (MCMC) simulations in order to constrain the bounds on graviton mass in the case of a small deviation of the stellar orbits from the corresponding GR predictions and showed that our method could further improve previous estimates for upper bounds on the graviton mass. It is also demonstrated that such analysis of the observed orbits of S-stars around the GC in the frame of the Yukawa gravity represents a tool for constraining the upper bound for the graviton mass, as well as for probing the predictions of GR or other gravity theories.

Bow-shocks are produced in the local interstellar medium by the passage of fast stars from the Galactic thin-disk and thick-disk populations with velocities $V_* = $ 40-80 km/s. Stellar transits of local H I clouds occur every 3500-7000 yr on average and last between $10^4$ and $10^5$ yr. There could be 10-20 active bow shocks around low-mass stars inside clouds within 10-15 pc of the Sun. At local cloud distances of 3-10 pc, their turbulent wakes have transverse radial extents $R_{\rm wake} \approx$ 10-300 AU, angular sizes 10-100 arcsec, and Lyman-alpha surface brightnesses of 2-8 Rayleighs in gas with total hydrogen density $n_H \approx 0.1~{\rm cm}^{-3}$ and $V_* =$ 40-80 km/s. These transit wakes may cover an area fraction $f_A \approx (R_{\rm wake}/R_{\rm cl}) \approx 10^{-3}$ of local H I clouds and be detectable in IR (dust), UV (Lya, two-photon), or non-thermal radio emission. Turbulent heating in these wakes could produce the observed elevated rotational populations of H$_2$ ($J \geq 2$) and influence the endothermic formation of CH$^+$ in diffuse interstellar gas at $T > 10^3$ K.

The effects of the anisotropy on the fluid pulsation modes adopting the so-called Cowling approximation and tidal deformability of strange quark stars are investigated by using the numerical integration of the hydrostatic equilibrium, nonradial oscillations, and tidal deformability equations, being these equations modified from their standard form to include the anisotropic effects. The fluid matter inside the compact stars is described by the MIT bag model equation of state. For the anisotropy profile, we consider a local anisotropy that is both regular at the center and null at the star's surface. We find that the effect of the anisotropy is reflected in the fluid pulsation modes and tidal deformability. Finally, we analyze the correlation between the tidal deformability of the GW$170817$ event with the anisotropy.

Enceladus' plume consists mainly of a mixture of water vapor and solid ice particles that may originate from a subsurface ocean. The physical processes underlying Enceladus' plume particle dynamics are still being debated, and quantifying the particles' size distribution and launch velocities can help constrain these processes. Cassini's Visual and Infrared Mapping Spectrometer (VIMS) observed the Enceladus plume over a wavelength range of 0.9 micron to 5.0 microns for a significant fraction of Enceladus' orbital period on three dates in the summer of 2017. We find that the relative brightness of the plume on these different dates varies with wavelength, implying that the particle size distribution in the plume changes over time. These observations also enable us to study how the particles' launch velocities vary with time and observed wavelength. We find that the typical launch velocity of particles remains between 140 m/s and 148 m/s at wavelengths between 1.2 microns and 3.7 microns. This may not be consistent with prior models where particles are only accelerated by interactions with the vent walls and gas, and could imply that mutual particle collisions close to the vent are more important than previously recognized.

The magnetic field conditions in astrophysical relativistic jets can be probed by multiwavelength polarimetry, which has been recently extended to X-rays. For example, one can track how the magnetic field changes in the flow of the radiating particles by observing rotations of the electric vector position angle $\Psi$. Here we report the discovery of a $\Psi_{\mathrm x}$ rotation in the X-ray band in the blazar Mrk 421 at an average flux state. Across the 5 days of Imaging X-ray Polarimetry Explorer (IXPE) observations of 4-6 and 7-9 June 2022, $\Psi_{\mathrm x}$ rotated in total by $\geq360^\circ$. Over the two respective date ranges, we find constant, within uncertainties, rotation rates ($80 \pm 9$ and $91 \pm 8 ^\circ/\rm day$) and polarization degrees ($\Pi_{\mathrm x}=10\%\pm1\%$). Simulations of a random walk of the polarization vector indicate that it is unlikely that such rotation(s) are produced by a stochastic process. The X-ray emitting site does not completely overlap the radio/infrared/optical emission sites, as no similar rotation of $\Psi$ was observed in quasi-simultaneous data at longer wavelengths. We propose that the observed rotation was caused by a helical magnetic structure in the jet, illuminated in the X-rays by a localized shock propagating along this helix. The optically emitting region likely lies in a sheath surrounding an inner spine where the X-ray radiation is released.

Fluxes of the diffuse supernova neutrino background (DSNB) are calculated based on a new modeling of galactic chemical evolution, where a variable stellar initial mass function (IMF) depending on the galaxy type is introduced and black hole (BH) formation from the failed supernova is considered for progenitors heavier than 18$M_{\odot}$. The flux calculations are performed for different combinations of the star formation rate, nuclear equation of state, and neutrino mass hierarchy to examine the systematic effects from these factors. In any case, our new model predicts the enhanced DSNB $\bar{\nu}_{e}$ flux at $E_\nu \gtsim 30$~MeV and $E_\nu \ltsim 10$~MeV due to more frequent BH formation and a larger core collapse rate at high redshifts in the early-type galaxies, respectively. Event rate spectra of the DSNB $\bar{\nu}_{e}$ at a detector from the new model are shown and detectability at water-based Cherenkov detectors, SK-Gd and Hyper-Kamiokande, is discussed. In order to investigate impacts of the assumptions in the new model, we prepare alternative models based on the different IMF form and treatment of BH formation, and estimate discrimination capabilities between the new and alternative models at these detectors.

The physical conditions of the circumgalactic medium are probed by intervening absorption-line systems in the spectrum of background quasi-stellar objects out to the epoch of cosmic reionization. A correlation between the ionization state of the absorbing gas and the nature of the nearby galaxies has been suggested by the sources detected either in Lyalpha or [C ii] 158 m near to respectively highly-ionized and neutral absorbers. This is also likely linked to the global changes in the incidence of absorption systems of different types and the process of cosmic reionization. Here we report the detection of two [C ii]-emitting galaxies at redshift $z \sim 5.7$ that are associated with a complex high-ionization C iv absorption system. These objects are part of an overdensity of galaxies and have compact sizes (< 2.4 kpc) and narrow line widths (FWHM $\sim$ 62--64 km s-1). Hydrodynamic simulations predict that similar narrow [C ii] emission may arise from the heating of small ($\lesssim$ 3 kpc) clumps of cold neutral medium or a compact photodissociation region. The lack of counterparts in the rest-frame ultraviolet indicates severe obscuration of the sources that are exciting the [C ii] emission. These results may suggest a connection between the properties of the [C ii] emission, the rare overdensity of galaxies and the unusual high ionization state of the gas in this region.

The Cassini mission provided key measurements needed to determine the absolute age of Saturn's rings, including the extrinsic micrometeoroid flux at Saturn, the volume fraction of non-icy pollutants in the rings, and the total ring mass. These three factors constrain the ring age to be no more than a few 100 Myr (Kempf et al., 2023). Observations during the Cassini Grand Finale also showed that the rings are losing mass to the planet at a prodigious rate. Some of the mass flux falls as "ring rain" at high latitudes. However, the influx in ring rain is considerably less than the total measured mass influx of 4800 to 45000 kg/s at lower latitudes (Waite et al., 2018). In addition to polluting the rings, micrometeoroid impacts lead to ballistic transport, the mass and angular momentum transport due to net exchanges of meteoroid impact ejecta. Because the ejecta are predominantly prograde, they carry net angular momentum outward. As a result, ring material drifts inward toward the planet. Here, for the first time, we use a simple model to quantify this radial mass inflow rate for dense rings and find that, for plausible choices of parameters, ballistic transport and mass loading by meteoroids can produce a total inward flux of material in the inner B ring and in the C ring that is on the order of a few x 10^3 to a few x 10^4 kg/s, in agreement with measurements during the Cassini Grand Finale. From these mass inflow rates, we estimate that the remaining ring lifetime is ~15 to 400 Myr. Combining this with a revised pollution age of ~120 Myr, we conclude that Saturn's rings are not only young but ephemeral and probably started their evolution on a similar timescale to their pollution age with an initial mass of one to a few Mimas masses.

The young, compact, very high surface brightness but low excitation planetary nebula (PN) BD+303639 is one of the very few PNe that have been reported to exhibit the 21cm HI emission line. As part of a long-term programme to search for circumstellar atomic hydrogen, we observed the 21cm feature toward BD+303639 with the Five-hundred-meter Aperture Spherical radio Telescope (FAST). Assuming a direct association between the PN and the detected HI emission, these new observations show that this surrounding emission is significantly more spatially extended than indicated by previous interferometric observations, and can be resolved into two velocity components. The estimated HI mass is larger than 100M_sun, invalidating an origin from the host star itself or its ejecta for the emitting material. We discuss the possibility that the extended HI emission stems from the interstellar medium (ISM) swept out over time by the stellar wind. Moreover, we report tentative detections of HI absorption features lying near and blueward of the systemic velocity of this PN, which are probably from a stalled asterosphere at the outer boundary of the expanding ionized region. The mass of the gas producing the HI absorption is insufficient to solve the so-called `PN missing mass problem'. We demonstrate the capability of FAST to investigate the interaction process between a PN and the surrounding ISM.

The Cassini spacecraft provided key measurements during its more than twelve year mission that constrain the absolute age of Saturn's rings. These include the extrinsic micrometeoroid flux at Saturn, the volume fraction of non-icy pollutants in the rings, and a measurement of the ring mass. These observations taken together limit the ring exposure age to be < a few 100 Myr if the flux was persistent over that time (Kempf et al., 2023). In addition, Cassini observations during the Grand Finale further indicate the rings are losing mass (Hsu et al., 2018; Waite et al., 2018) suggesting the rings are ephemeral as well. In a companion paper (Durisen and Estrada, 2023), we show that the effects of micrometeoroid bombardment and ballistic transport of their impact ejecta can account for these loss rates for reasonable parameter choices. In this paper, we conduct numerical simulations of an evolving ring in a systematic way in order to determine initial conditions that are consistent with these observations.

We present timing solutions for 21 pulsars discovered in 350 MHz surveys using the Green Bank Telescope (GBT). All were discovered in the Green Bank North Celestial Cap pulsar survey, with the exception of PSR J0957-0619, which was found in the GBT 350 MHz Drift-scan pulsar survey. The majority of our timing observations were made with the GBT at 820 MHz. With a spin period of 37 ms and a 528-day orbit, PSR J0032+6946 joins a small group of five other mildly recycled wide binary pulsars, for which the duration of recycling through accretion is limited by the length of the companion's giant phase. PSRs J0141+6303 and J1327+3423 are new disrupted recycled pulsars. We incorporate Arecibo observations from the NANOGrav pulsar timing array into our analysis of the latter. We also observed PSR J1327+3423 with the Long Wavelength Array, and our data suggest a frequency-dependent dispersion measure. PSR J0957-0619 was discovered as a rotating radio transient, but is a nulling pulsar at 820 MHz. PSR J1239+3239 is a new millisecond pulsar (MSP) in a 4-day orbit with a low-mass companion. Four of our pulsars already have published timing solutions, which we update in this work: the recycled wide binary PSR J0214+5222, the non-eclipsing black widow PSR J0636+5128, the disrupted recycled pulsar J1434+7257, and the eclipsing binary MSP J1816+4510, which is in an 8.7 hr orbit with a redback-mass companion.

The galaxy-galaxy lensing technique allows us to measure the subhalo mass of satellite galaxies, studying their mass loss and evolution within galaxy clusters and providing direct observational validation for theories of galaxy formation. In this study, we use the weak gravitational lensing observations from DECaLS DR8, in combination with the redMaPPer galaxy cluster catalog from SDSS DR8 to accurately measure the dark matter halo mass of satellite galaxies. We confirm a significant increase in the stellar-to-halo mass ratio of satellite galaxies with their halo-centric radius, indicating clear evidence of mass loss due to tidal stripping. Additionally, we find that this mass loss is strongly dependent on the mass of the satellite galaxies, with satellite galaxies above $10^{11}~{\rm M_{\odot}/h}$ experiencing more pronounced mass loss compared to lower mass satellites, reaching 86\% at projected halo-centric radius $0.5R_{\rm 200c}$. The average mass loss rate, when not considering halo-centric radius, displays a U-shaped variation with stellar mass, with galaxies of approximately $4\times10^{10}~{\rm M_{\odot}/h}$ exhibiting the least mass loss, around 60\%. We compare our results with state-of-the-art hydrodynamical numerical simulations and find that the satellite galaxy stellar-to-halo mass ratio in the outskirts of galaxy clusters is higher compared to the predictions of the Illustris-TNG project about factor 5. Furthermore, the Illustris-TNG project's numerical simulations did not predict the observed dependence of satellite galaxy mass loss rate on satellite galaxy mass.

The non linear correlation between the UV and X-ray emission observed in Active Galactic Nuclei remains a puzzling question that challenged accretion models. While the UV emission originates from the cold disk, the X-ray emission is emitted by a hot corona whose physical characteristics and geometry are still highly debated. The Jet Emitting Disk - Standard Accretion Disk (JED-SAD) is a spectral model stemming from self similar accretion-ejection solutions. It is composed of an inner highly magnetized and hot accretion flow launching jets, the JED, and an outer SAD. The model has been successfully applied to X-ray binaries outbursts. The AGN UV X-ray correlation represent another essential test for the JED-SAD model. We use multiple AGN samples to explore the parameter space and identify the regions able to reproduce the observations. In this first paper, we show that JED-SAD model is able to reproduce the UV--X-ray correlation.

We report the discovery of the hyperluminous, highly obscured AGN WISE J190445.04+485308.9 (W1904+4853 hereafter, $L_{bol} = 1.1 \times 10^{13} \ L_{\odot}$) at z=0.415. Its well-sampled spectral energy distribution (SED) is dominated by infrared dust emission, though broad emission lines are detected in the optical spectra. These features suggest that W1904+4853 contains an actively accreting supermassive black hole hidden in its dusty cocoon, resembling the observed properties of Hot Dust-Obscured Galaxies (Hot DOGs), a population previously only identified at z>1.0. Using the broad component of the MgII emission line, we estimate a black hole mass of $log \ (M_{BH}/M_{\odot}) = 8.4 \pm 0.3$. The corresponding Eddington ratio of $1.4 \pm 0.2$ implies that the central black hole accretion is at the theoretical limit of isotropic accretion. The rest-frame UV-optical SED and [O II] emission line also indicate that the host galaxy of W1904+4853 harbors strong star formation activity at the rate of up to $\sim 45 \ M_{\odot} \ yr^{-1}$. With an estimated stellar mass of $3 \times 10^{10} \ M_{\odot}$, the host galaxy appears to be a starburst system with respect to the main sequence of the star-forming galaxies at the same redshift. Although blueshifted and asymmetric [O III] emission provides evidence of an outflow, we estimate it to be an order of magnitude smaller than the star formation rate, indicating that the current obscured AGN activity at the center has not yet produced significant feedback on the host galaxy star formation activity. W1904+4853 supports the interpretation that Hot DOGs are a rare transitional phase of AGN accretion in galaxy evolution, a phase that can persist into the present-day Universe.

Variability carries physical patterns and astronomical information of objects, and stellar light curve variations are essential to understand the stellar formation and evolution processes. The studies of variations in stellar photometry have the potential to expand the list of known stars, protostars, binary stars, and compact objects, which could shed more light on stages of stellar lifecycles. The progress in machine-learning techniques and applications has developed modern algorithms to detect and condense features from big data, which enables us to classify stellar light curves efficiently and effectively. We explore several deep-learning methods on variable star classifications. The sample of light curves is constructed with $\delta$ Scuti, $\gamma$ Doradus, RR Lyrae, eclipsing binaries, and hybrid variables from \textit{Kepler} observations. Several algorithms are applied to transform the light curves into images, continuous wavelet transform (CWT), Gramian angular fields, and recurrent plots. We also explore the representation ability of these algorithms. The processed images are fed to several deep-learning methods for image recognition, including VGG-19, GoogLeNet, Inception-v3, ResNet, SqueezeNet, and Xception architectures. The best transformation method is CWT, resulting in an average accuracy of 95.6\%. VGG-19 shows the highest average accuracy of 93.25\% among all architectures, while it shows the highest accuracy of 97.2\% under CWT transformation method. The prediction can reach $\sim1000$ light curves per second by using NVIDIA RTX 3090. Our results indicate that the combination of big data and deep learning opens a new path to classify light curves automatically.

A tau-Herculid meteor outburst or even a storm was predicted by several models to occur around 5~UT on 31~May, 2022 as a consequence of the break-up of comet 73P/Schwassmann-Wachmann 3 in 1995. The multi-instrument and multi-station experiment was carried-out within the Czech Republic to cover possible earlier activity of the shower between 21 and 1 UT on 30/31 May. Multi-station observations using video and photographic cameras were used for calculation of the atmospheric trajectories and heliocentric orbits of the meteors. Their arrival times are used for determination of the shower activity profile. Physical properties of the meteoroids are evaluated using various criteria based on meteor heights. Evolution of spectra of three meteors are studied as well. This annual but poor meteor shower was active for the whole night many hours before the predicted peak. A comparison with dynamical models shows that a mix of older material ejected after 1900 and fresh particles originating from the 1995 comet fragmentation event was observed. Radiant positions of both groups of meteors were identified and found to be in good agreement with simulated radiants. Meteoroids with masses between 10 mg and 10 kg were recorded. The mass distribution index was slightly higher than 2. A study of the physical properties shows that the tau-Herculid meteoroids belong to the most fragile particles observed ever, especially among higher masses of meteoroids. Exceptionally bright bolide observed during the dawn represents a challenge for the dynamical simulations as it is necessary to explain how to transfer a half metre body to the vicinity of the Earth at the same time as millimetre sized particles.

Using ASAS-SN Sky Patrol Photometic Database and Asteroid Terrestrial-impact Last Alert System (ATLAS) data, I found that the high-field polar AR UMa entered a long-lasting high state in 2022 October. This object is renowned for its small duty cycle, and short-lived high states have only been occasionally seen since the discovery. It appears that the present long-lasting high state is the first recorded one at least in the last 30 years and probably even more. Before entering the current long-lasting high state, this object showed three short-lived high states, which might have been precursors to the current state. Before these short-lived high states, the object had been in a low state for 8 years and probably more. I refined the orbital period to be 0.08050066(1) d. The object is still bright and current phenomenon provides a unique opportunity to study accretion processes onto a strongly magnetized white dwarf and to study the mechanism of maintaining the long-lasting high state.

Recent observations of the Milky Way (MW) found an unexpected steepening of the star-forming gas metallicity gradient around the time of the Gaia-Sausage-Enceladus (GSE) merger event. Here we investigate the influence of early ($t_{\mathrm{merger}}\lesssim5$ Gyr) massive ($M_{\mathrm{gas}}^{\mathrm{merger}}/M_{\mathrm{gas}}^{\mathrm{main}}(t_{\mathrm{merger}})\gtrsim10\%$) merger events such as the Gaia-Sausage Enceladus merger in the MW on the evolution of the cold gas metallicity gradient. We use the NIHAO-UHD suite of cosmological hydrodynamical simulations of MW-mass galaxies to study the frequency of massive early mergers and their detailed impact on the morphology and chemistry of the gaseous disks. We find a strong steepening of the metallicity gradient at early times for all four galaxies in our sample which is caused by a sudden increase in the cold gas disk size (up to a factor of 2) in combination with the supply of un-enriched gas ($\sim0.75$ dex lower compared to the main galaxy) by the merging dwarf galaxies. The mergers mostly affect the galaxy outskirts and lead to an increase in cold gas surface density of up to 200% outside of $\sim8$ kpc. The addition of un-enriched gas breaks the self-similar enrichment of the inter-stellar-medium and causes a dilution of the cold gas in the outskirts of the galaxies. The accreted stars and the ones formed later out of the accreted gas inhabit distinct tracks offset to lower [$\alpha$/Fe] and [Fe/H] values compared to the main galaxy's stars. We find that such mergers can contribute significantly to the formation of a second, low-$\alpha$ sequence as is observed in the MW.

Simulated galaxy distributions are suitable for developing filament detection algorithms. However, samples of observed galaxies, being of limited size, cause difficulties that lead to a discontinuous distribution of filaments. We created a new galaxy filament catalog composed of a continuous cosmic web with no lone filaments. The core of our approach is a ridge filter used within the framework of image analysis. We considered galaxies from the HyperLeda database with redshifts $0.02\leqslant z\leqslant 0.1$, and in the solid angle $120^\circ\leqslant {\rm RA}\leqslant 240^\circ$, $0^\circ\leqslant {\rm DEC}\leqslant 60^\circ$. We divided the sample into 16 two-dimensional celestial projections with redshift bin $\Delta z=0.005$, and compared our continuous filament network with a similar recent catalog covering the same region of the sky. We tested our catalog on two application scenarios. First, we compared the distributions of distance to nearest filament of various astrophysical sources (Seyfert galaxies and other active galactic nuclei, radio galaxies, low surface brightness galaxies, and dwarf galaxies), and found that all source types trace the filaments well, with no systematic differences. Next, among the HyperLeda galaxies, we investigated the dependence of $g-r$ color distribution on distance to nearest filament, and confirmed that early type galaxies are located on average further from the filaments than late type ones.

Machine learning has become an effective tool for processing the extensive data sets produced by large physics experiments. Gravitational-wave detectors are now listening to the universe with quantum-enhanced sensitivity, accomplished with the injection of squeezed vacuum states. Squeezed state preparation and injection is operationally complicated, as well as highly sensitive to environmental fluctuations and variations in the interferometer state. Achieving and maintaining optimal squeezing levels is a challenging problem and will require development of new techniques to reach the lofty targets set by design goals for future observing runs and next-generation detectors. We use machine learning techniques to predict the squeezing level during the third observing run of the Laser Interferometer Gravitational-Wave Observatory (LIGO) based on auxiliary data streams, and offer interpretations of our models to identify and quantify salient sources of squeezing degradation. The development of these techniques lays the groundwork for future efforts to optimize squeezed state injection in gravitational-wave detectors, with the goal of enabling closed-loop control of the squeezer subsystem by an agent based on machine learning.

Using the multiwavelength data sets, we studied the star formation activity in H II region Sh 2-61 (hereafter S61). We identified a clustering in the region and estimated the membership using the Gaia proper motion data. The physical environment of S61 is inspected using infrared to radio wavelength images. We also determined the Lyman continuum flux associated with the H II region and found that the H II region is formed by at least two massive stars (S1 and S2). We also analyzed the 12CO (J =3-2) JCMT data of S61, and a shell structure accompanying three molecular clumps are observed towards S61. We found that the ionized gas in S61 is surrounded by dust and a molecular shell. Many young stellar objects and three molecular clumps are observed at the interface of the ionized gas and the surrounding gas. The pressure at the interface is higher than in a typical cool molecular cloud.

The measurement of the cross-correlation function is crucial to assess magnification bias in galaxy surveys. Previous works used mini-tile subsampling, but accurately determining the integral constraint (IC) correction for unbiased estimation is challenging due to various factors. We present a new methodology for estimating the cross-correlation function, utilizing full field area and reducing statistical uncertainty. Covariance matrices were estimated by dividing each field into at least five patches using a k-mean clustering algorithm. Robustness was assessed by comparing spectroscopic and photometric lens samples, yielding compatible results. Cross-correlation and auto-correlation analyses in the GAMA fields revealed a stronger signal in GAMA15, likely due to rare large-scale structure combinations. Our findings highlight the robustness of the new methodology and suggest sample-specific effects. Subsequent papers in this series will explore other aspects of magnification bias and address potential biases from the GAMA15 signal on cosmological parameter constraints.

The main goal of this work, the second in a three-paper series, is to test the impact of a methodological improvement in measuring the magnification bias signal on a sample of submillimeter galaxies and its implications for constraining cosmological parameters. The analysis considers the angular cross-correlation function between a foreground sample of GAMA galaxies ($0.2<z<0.8$) and a background sample of H-ATLAS submillimeter galaxies ($1.2<z<4.0$). A refined methodology, discussed extensively in Paper I, is used. By interpreting the weak lensing signal within the halo model and employing an MCMC algorithm, the posterior distribution of the halo occupation distribution (HOD) and cosmological parameters is obtained for a flat $\Lambda$CDM model. The analysis incorporates the foreground angular auto-correlation function to account for galaxy clustering. The results demonstrate a remarkable improvement in uncertainties for both HOD and cosmological parameters compared to previous studies. However, when using the cross-correlation data alone, the estimation of $\sigma_8$ depends on prior knowledge of $\beta$, the logarithmic slope of the background number counts. Assuming a physically motivated prior distribution for $\beta$, mean values of $\Omega_m=0.18^{+0.03}{-0.03}$ and $\sigma_8=1.04^{+0.11}{-0.07}$ are obtained. These results may however be subject to an inherent bias in the data due to anomalous behavior observed in the G15 field. After excluding the G15 region, the mean values shift to $\Omega_m=0.30^{+0.05}{-0.06}$ and $\sigma_8=0.92^{+0.07}{-0.07}$. This becomes more apparent when adding the clustering of the foreground sample, but the dependence on $\beta$ information disappears, mitigating the aforementioned issue. Excluding the G15 region, the joint analysis yields mean values of $\Omega_m=0.36^{+0.03}{-0.07}$, $\sigma_8=0.90^{+0.03}{-0.03}$, and $h=0.76^{+0.14}_{-0.14}$.

This paper is the third in a series on submillimeter galaxy magnification bias, focusing on the tomographic scenario. It refines the methodology used to constrain the halo occupation distribution model and cosmological parameters within a flat $\Lambda$CDM model, using updated data. The study aims to optimize CPU time, explore strategies for analyzing different redshift bins, and assess the impact of excluding the GAMA15 field. The tomographic approach involves dividing the redshift range into bins and analyzing cross-correlation measurements between submillimeter and foreground galaxies. The results show good agreement between the mean-redshift and full model cases, with an increase in the minimum mass of lenses at higher redshifts. The inferred cosmological parameters have narrower posterior distributions, indicating reduced measurement uncertainties compared to previous studies. Excluding the GAMA15 field reduces the cross-correlation signal, suggesting sample variance within the large-scale structure. Extending the redshift range improves robustness against sample variance and produces similar but tighter constraints. The study highlights the importance of sample variance and redshift binning in tomographic analyses, and suggests using additional wide-area fields and updated foreground catalogues for more effective implementation.

We report the first $> 99\%$ confidence detection of X-ray polarization in BL Lacertae. During a recent X-ray/$\gamma$-ray outburst, a 287 ksec observation (2022 November 27-30) was taken using the Imaging X-ray Polarimetry Explorer ({\it IXPE}), together with contemporaneous multiwavelength observations from the Neil Gehrels {\it Swift} observatory and {\it XMM-Newton} in soft X-rays (0.3--10~keV), {\it NuSTAR} in hard X-rays (3--70~keV), and optical polarization from the Calar Alto, and Perkins Telescope observatories. Our contemporaneous X-ray data suggest that the {\it IXPE} energy band is at the crossover between the low- and high-frequency blazar emission humps. The source displays significant variability during the observation, and we measure polarization in three separate time bins. Contemporaneous X-ray spectra allow us to determine the relative contribution from each emission hump. We find $>99\%$ confidence X-ray polarization $\Pi_{2-4{\rm keV}} = 21.7^{+5.6}_{-7.9}\%$ and electric vector polarization angle $\psi_{2-4{\rm keV}} = -28.7 \pm 8.7^{\circ}$ in the time bin with highest estimated synchrotron flux contribution. We discuss possible implications of our observations, including previous {\it IXPE} BL Lacertae pointings, tentatively concluding that synchrotron self-Compton emission dominates over hadronic emission processes during the observed epochs.

Our knowledge of stellar evolution is driven by one-dimensional (1D) simulations. 1D models, however, are severely limited by uncertainties on the exact behaviour of many multi-dimensional phenomena occurring inside stars, affecting their structure and evolution. Recent advances in computing resources have allowed small sections of a star to be reproduced with multi-D hydrodynamic models, with an unprecedented degree of detail and realism. In this work, we present a set of 3D simulations of a convective neon-burning shell in a 20 M$_\odot$ star run for the first time continuously from its early development through to complete fuel exhaustion, using unaltered input conditions from a 321D-guided 1D stellar model. These simulations help answer some open questions in stellar physics. In particular, they show that convective regions do not grow indefinitely due to entrainment of fresh material, but fuel consumption prevails over entrainment, so when fuel is exhausted convection also starts decaying. Our results show convergence between the multi-D simulations and the new 321D-guided 1D model, concerning the amount of convective boundary mixing to include in stellar models. The size of the convective zones in a star strongly affects its structure and evolution, thus revising their modelling in 1D will have important implications for the life and fate of stars. This will thus affect theoretical predictions related to nucleosynthesis, supernova explosions and compact remnants.

Jupiter's deep abundances help to constrain the formation history of the planet and the environment of the protoplanetary nebula. Juno recently measured Jupiter's deep oxygen abundance near the equator to be 2.2$_{-2.1}^{+3.9}$ times the protosolar value (2$\sigma$ uncertainties). Even if the nominal value is supersolar, subsolar abundances cannot be ruled out. Here we use a state-of-the-art one-dimensional thermochemical and diffusion model with updated chemistry to constrain the deep oxygen abundance with upper tropospheric CO observations. We find a value of 0.3$_{-0.2}^{+0.5}$ times the protosolar value. This result suggests that Jupiter could have a carbon-rich envelope that accreted in a region where the protosolar nebula was depleted in water. However, our model can also reproduce a solar/supersolar water abundance if vertical mixing is reduced in a radiative layer where the deep oxygen abundance is obtained. More precise measurements of the deep water abundance are needed to discriminate between these two scenarios and understand Jupiter's internal structure and evolution.

Pisces VII/Triangulum III (Pisc~VII) was discovered in the DESI Legacy Imaging Survey and was shown to be a Local Group dwarf galaxy with follow-up imaging from the 4-m Telescopio Nazionale Galileo. However, this imaging was unable to reach the horizontal branch of Pisc VII, preventing a precision distance measurement. The distance bound from the red giant branch population placed Pisc VII as either an isolated ultra-faint dwarf galaxy or the second known satellite galaxy of Triangulum (M33). Using deep imaging from Gemini GMOS-N, we have resolved the horizontal branch of Pisc VII, and measure a distance of $D=962^{+32}_{-32}$~kpc, making Pisc VII a likely satellite of M33. We also remeasure its size and luminosity from this deeper data, finding $r_{\rm half}=186^{+58}_{-32}$ pc, $M_V=-5.7\pm0.3$ and $L=1.6^{+0.1}_{-0.2}\times10^4\,{\rm L}_\odot$. Given its position in the M33 halo, we argue that Pisc VII could support the theory that M33 is on its first infall to the Andromeda system. We also discuss the presence of blue stars in the colour-magnitude diagram of Pisc VII that are consistent with ages of 1.5 Gyr. If these are truly members of the galaxy, it would transform our understanding of how reionisation affects the faintest galaxies. However we cannot rule out a more ordinary explanation for these with current data. Future deep imaging and dynamics could allow significant insight into both the stellar populations of Pisc VII and the evolution of M33.

We look for possible evidence of a non-minimal coupling (NMC) between dark matter (DM) and gravity using data from the X-COP compilation of galaxy clusters. We consider a theoretically motivated NMC that may dynamically arise from the collective behavior of the coarse-grained DM field (e.g., via Bose-Einstein condensation) with averaging/coherence length $\L_{\mathrm{nmc}}$. In the Newtonian limit, the NMC modifies the Poisson equation by a term $\L_{\mathrm{nmc}}^2 \nabla^2 \rho$ proportional to the Laplacian of the DM density itself. We show that this term when acting as a perturbation over the standard Navarro-Frenk-White (NFW) profile of cold DM particles, can yield DM halo density profiles capable of correctly fitting galaxy clusters' pressure profiles with an accuracy comparable and in some cases even better than the standard cold DM NFW profile. We also show that the observed relation between the non-minimal coupling length scale and the virial mass found in Gandolfi et al., 2022 for Late Type Galaxies is consistent with the relation we find in the current work, suggesting that the previously determined power-law scaling law holds up to galaxy cluster mass scales.

We calculate the effect of gravitational lensing on the parity-odd power spectrum of the cosmic microwave background (CMB) polarization induced by axionlike particles (ALPs). Several recent works have reported a tantalizing hint of cosmic birefringence, a rotation of the linear polarization plane of CMB, which ALPs can explain. In future CMB observations, we can measure cosmic birefringence more precisely to get insight into ALPs. We find that the lensing effect is necessary to fit the observed EB power spectrum induced by cosmic birefringence in future CMB observations, including Simons Observatory and CMB-S4. We also show that the estimated ALPs parameters are biased if we ignore the lensing effect. Therefore, the lensing correction to the parity-odd power spectra must be included in future high-resolution CMB experiments.

Extragalactic jets are generated as bipolar outflows at the nuclei of active galaxies. Depending on their morphology, they are classified as Fanaroff-Riley type I (centre-brightened) and Fanaroff-Riley type II (edge-brightened) radio jets. However, this division is not sharp and observations of these sources at large scales often show intermediate jet morphologies or even hybrid jet morphologies with a FRI type jet on one side and a FRII type jet on the other. A good example of a radio galaxy that is difficult to classify as FRI or FRII is Hercules~A. This source shows jets with bright radio lobes (a common feature of FRII type jets) albeit without the hotspots indicative of the violent interaction between the jet and the ambient medium at the impact region, because the jets seem to be disrupted inside the lobes at a distance from the bow shocks surrounding the lobes. In this paper, we explore the jet physics that could trigger this peculiar morphology by means of three-dimensional relativisitic hydrodynamical simulations. Our results show that the large-scale morphological features of Hercules A jets and lobes can be reproduced by the propagation of a relativistically hot plasma outflow that is disrupted by helical instability modes, and generates a hot lobe that expands isotropically against the pressure-decreasing intergalactic medium. We also discuss the implications that this result may have for the host active nucleus in terms of a possible transition from high-excitation to low-excitation galaxy modes.

In close binary systems, tidal interactions and rotational effects can strongly influence stellar evolution as a result of mass-transfer, common envelope phases, ... All these aspects can only be treated following improvements of theoretical models, taking into account the breaking of spherical symmetry occurring in close binaries. Current models of binary stars are relying either on the so-called "Roche model" or the perturbative approach that in each case results on several assumptions concerning the gravitational, tidal and centrifugal potentials.We developed a new non-perturbative method to compute precise structural deformation of binary system in three dimensions that is valid even in the most distorted cases. We then compared our new method to the Roche and perturbative models for different orbital separations and binary components. We found that in the most distorted cases both Roche and perturbative models are significantly underestimating the deformation of binaries. The effective gravity and the overall structural deformations are also noticeably different in the most distorted cases leading, for the interpretation of observations, to modifications of the usual gravity darkening generally obtained through the Roche model. Moreover we found that the dipolar term of the gravitational potential, usually neglected by the perturbative theory, has the same order of magnitude than the leading tidal term in the most distorted cases. We developed a new method that is capable of precisely computing the deformations of binary system composed of any type of stars, even compact objects. For all stars studied the differences in deformation with respect to the Roche or perturbative models are significant in the most distorted cases impacting both the interpretation of observations and the theoretical structural depiction of these distorted bodies.

We present a new study of the Z~Cam-type eclipsing cataclysmic variable AY~Piscium with the aim of determining the fundamental parameters of the system and the structure of the accretion flow therein. We use time-resolved photometric observations supplemented by spectroscopy in the standstill, to which we applied our light-curve modeling techniques and the Doppler tomography method, to update system parameters. We found that the system has a massive white dwarf $M_{\rm WD}=0.90(4)$ \ms, a mass ratio $q=0.50(3)$, and the effective temperature of a secondary $T_2 = 4100(50)$~K. The system inclination is $i=74.^{\circ}8(7)$. The orbital period of the system $P_{\mathrm{orb}}=0.217320523(8)\;\mathrm{d}$ is continuously increasing with the rate of $\dot{P}_{\mathrm{orb}} = +7.6(5)\times10^{-9}$ d year$^{-1}$. The mass transfer rate varies between 2.4$\times$10$^{-10}$ M$_\odot$ year$^{-1}$ in quiescence up to 1.36$\times$10$^{-8}$ M$_\odot$ year$^{-1}$ in outburst. The accretion disk transitions from the cooler, flared, steady-state disk to a warmer state with a practically constant and relatively high disk height. The mass transfer rate is about 1.6$\times$10$^{-9}$ M$_\odot$ year$^{-1}$ in the standstill. The Balmer emission lines show a multi-component structure similar to that observed in long-orbital-period nova-like systems. Out of standstill, the system exhibits outburst bimodality, with long outbursts being more prominent. We conclude that the Balmer emission lines in AY~Psc are formed by the combination of radiation from the irradiated surface of the secondary, from the outflow zone, and from winds originating in the bright spot and the disk's inner part.

Interacting dark sectors may undergo changes in the number of their relativistic species during the early universe, due to a mass threshold $m$ (similar to changes in the Standard Model bath), and in doing so affect the cosmic history. When such changes occur close to recombination, i.e., for $m\sim (0.1-10)~\text{eV}$, the stringent bound on the effective number of neutrino species, $N_{\text{eff}}$, can be relaxed and the value of the Hubble expansion rate $H_0$ inferred from Cosmic Microwave Background (CMB) observations raised. We search for such sectors (with and without mass thresholds) in the latest cosmological datasets, including the full-shape (FS) of BOSS DR12 galaxy power spectrum. We perform a detailed analysis, accounting for the choice of prior boundaries and additionally exploring the possible effects of dark sector interactions with (a fraction of) the dark matter. We find $\Delta N_{\text{eff}}\leq 0.55\, (0.46)$ at 95% C.L. with (without) a mass threshold. While a significantly larger Hubble rate is achieved in this scenario, $H_0=69.01^{+0.66}_{-1.1}$, the overall fit to CMB+FS data does not provide a compelling advantage over the $\Lambda$CDM model. Furthermore, we find that dark matter interactions with the dark sector do not significantly improve the (matter fluctuations) $S_8$ tension with respect to the $\Lambda$CDM model. Our work provides model-independent constraints on (decoupled) dark sectors with mass thresholds around the eV scale.

Mean motion resonances (MMR) are a frequent phenomenon among extrasolar planetary systems. Current observations indicate that many systems have planets that are close to or inside the 2:1 MMR, when the orbital period of one of the planets is twice the other. Analytical models to describe this particular MMR can only be reduced to integrable approximations in a few specific cases. While there are successful approaches to the study of this MMR in the case of very elliptic and/or very inclined orbits using semi-analytical or semi-numerical methods, these may not be enough to completely understand the resonant dynamics. In this work, we propose to apply a well-established numerical method to assess the global portrait of the resonant dynamics, which consists in constructing dynamical maps. Combining these maps with the results from a semi-analytical method, helps to better understand the underlying dynamics of the 2:1 MMR, and to identify the behaviors that can be expected in different regions of the phase space and for different values of the model parameters. We verify that the family of stable resonant equilibria bifurcate from symmetric to asymmetric librations, depending on the mass ratio and eccentricities of the resonant planets pair. This introduces new structures in the phase space, that turns the classical V-shape of the MMR, in the semi-major axis vs. eccentricity space, into a sand clock shape. We construct dynamical maps for three extrasolar planetary systems, TOI-216, HD27894, and K2-24, and discuss their phase space structure and their stability in the light of the orbital fits available in the literature.

Understanding how the $\mathrm{H}_2$ molecule is formed under the chemical conditions of the interstellar media (ISM) is critical to the whole chemistry of it. Formation of $\mathrm{H}_2$ in the ISM requires a third body acting as a reservoir of energy. Polycyclic aromatic hydrocarbons (PAH's) are excellent candidates to play that role. In this work we simulated the collisions of hydrogen atoms with coronene to form $\mathrm{H}_2$ via the Eley-Rideal mechanism. To do so, we used Born-Oppenheimer (ab initio) Molecular Dynamics simulations. Our results show that that adsorption of H atoms and subsequent release of $\mathrm{H}_2$ readily happen on coronene for H atoms with kinetic energy as large as 1 eV. Special attention is paid to dissipation and partition of the energy released in the reactions. The capacity of coronene to dissipate collision and reaction energies depends varies with the reaction site. Inner sites dissipate energy easier and faster than edge sites, thus evidencing an interplay between the potential energy surface around the reaction center and its ability to cool the projectile. As for the the recombination of H atoms and the subsequent formation of $\mathrm{H}_{2}$, it is observed that $\sim 15\%$ of the energy is dissipated by the coronene molecule as vibrational energy and the remaining energy is carried by $\mathrm{H}_{2}$. The $\mathrm{H}_{2}$ molecules desorb from coronene with an excited vibrational state ($\upsilon \geq 3$), a large amount of translational kinetic energy ($\geq$ 0.4 eV) and with a small activation of the rotational degree of freedom.

Stellar obliquity, the angle between a planet's orbital axis and its host star's spin axis, traces the formation and evolution of a planetary system. In transiting exoplanet observations, only the sky-projected stellar obliquity can be measured, but this can be de-projected using an estimate of the stellar obliquity. In this paper, we introduce a flexible, hierarchical Bayesian framework that can be used to infer the stellar obliquity distribution solely from sky-projected stellar obliquities, including stellar inclination measurements when available. We demonstrate that while a constraint on the stellar inclination is crucial for measuring the obliquity of an individual system, it is not required for robust determination of the population-level stellar obliquity distribution. In practice, the constraints on the stellar obliquity distribution are mainly driven by the sky-projected stellar obliquities. When applying the framework to all systems with measured sky-projected stellar obliquity, which are mostly Hot Jupiter systems, we find that the inferred population-level obliquity distribution is unimodal and peaked at zero degrees. The misaligned systems have nearly isotropic stellar obliquities with no strong clustering near 90 degrees. The diverse range of stellar obliquities prefers dynamic mechanisms, such as planet-planet scattering after a convergent disk migration, which could produce both prograde and retrograde orbits of close-in planets with no strong inclination concentrations other than 0 degrees.

This paper investigates the effects of non-vanishing spatial curvature on the propagation of primordial gravitational waves produced during inflation. In particular, we consider tensor perturbations over a homogeneous and isotropic background, and describe the propagation of gravitational waves in the de Sitter phase with spatially curved geometries. We thus derive the expression of the primordial power spectrum at the horizon crossing, in the case of open and closed universes. Then, we analyze how tensor modes propagate in the post-inflationary era, showing the evolution of transfer functions in the radiation and matter epochs, as well as the matching conditions in the intermediate regime. To account for the intrinsic nature of different relativistic species, we also explore the corrections to the standard behavior of the radiation energy density. For this purpose, we introduce the effective number of degrees of freedom of relativistic particles contributing to the primordial energy and entropy densities. Under the subhorizon approximation, we obtain the spectral energy density of relic gravitational waves in terms of the curvature density parameter. Finally, we discuss the capability of present and future experiments to detect the primordial gravitational wave signal at different frequency regimes.

Standard binary evolutionary models predict a significant population of core helium-burning stars that lost their hydrogen-rich envelope after mass transfer via Roche-lobe overflow. However, there is a scarcity of observations of such stripped stars in the intermediate mass regime (~1.5 - 8$ M_{\odot}$), which are thought to be prominent progenitors of SN Ib/c. Especially at low metallicity, a significant fraction of these stars is expected to be only partially stripped, retaining a significant amount of hydrogen on their surfaces. For the first time, we discovered a partially stripped massive star in a binary with a Be-type companion located in the Small Magellanic Cloud (SMC) using a detailed spectroscopic analysis. The stripped-star nature of the primary is revealed by the extreme CNO abundance pattern and very high luminosity-to-mass ratio, which suggest that the primary is likely shell-hydrogen burning. Our target SMCSGS-FS 69 is the most luminous and most massive system among the known stripped star + Be binaries, with Mstripped ~3$ M_{\odot}$ and MBe ~17$ M_{\odot}$. Binary evolutionary tracks suggest an initial mass of Mini $\gtrsim 12 M_{\odot}$ for the stripped star and predict it to be in a transition phase towards a hot compact He star, which will eventually produce a stripped-envelope supernova. Our target marks the first representative of a so-far missing evolutionary stage in the formation pathway of Be X-ray binaries and double neutron star mergers.

This research explores a cosmological model proposed in arXiv:2304.11492 by Jusufi et al., which challenges the existence of dark matter and proposes a connection between dark matter, dark energy, and baryonic matter by altering Newton's law of gravity in large-scale structures. The model is described by a Yukawa-like gravitational potential with the coupling parameter $\alpha$ and wavelength parameter $\lambda$. Utilizing observational data from Supernovae Ia and $H(z)$ within a $1\sigma$ CL, the study determines the best-fit parameters, resulting in $\lambda=2693_{-1262}^{+1191}\, \rm Mpc$, $\alpha=0.416_{-0.326}^{+1.137}$, and a graviton mass of approximately $m_g \simeq 10^{-69}$ kg. The research also establishes a direct connection between the dark matter/dark energy density parameters and the angular radius of the black hole shadow of the SgrA and M87 black holes in the low-redshift limit, which is consistent with the findings of the Event Horizon Telescope. This discovery holds significant importance in the paper, as it challenges the conventional understanding of the universe's composition and gravitational forces.

We report the discovery of COOL J0335$-$1927, a quasar at $z$ = 3.27 lensed into three images with a maximum separation of 23.3" by a galaxy cluster at $z$ = 0.4178. We construct a parametric strong gravitational lens model using ground-based imaging, constrained by the redshift and positions of the quasar images as well as the positions of three other multiply-imaged background galaxies. Using our best-fit lens model, we calculate the predicted time delays between the three quasar images to be $\Delta$t$_{AB}=$ $241^{+41}_{-12}$ and $\Delta$t$_{AC}=$ $-64^{+3}_{-33}$ days. We also present g-band photometry from archival DECaLS imaging, and new multi-epoch observations obtained between September 18, 2022 UT and February 22, 2023 UT, which demonstrate significant variability in the quasar and which will eventually enable a measurement of the time delay between the three quasar images.

A fundamental property of coronal mass ejections (CMEs) is their radial expansion, which determines the increase in the CME radial size and the decrease in the CME magnetic field strength as the CME propagates. CME radial expansion can be investigated either by using remote observations or by in-situ measurements based on multiple spacecraft in radial conjunction. However, there have been only few case studies combining both remote and in-situ observations. It is therefore unknown if the radial expansion estimated remotely in the corona is consistent with that estimated locally in the heliosphere. To address this question, we first select 22 CME events between the years 2010 and 2013, which were well observed by coronagraphs and by two or three spacecraft in radial conjunction. We use the graduated cylindrical shell model to estimate the radial size, radial expansion speed, and a measure of the dimensionless expansion parameter of CMEs in the corona. The same parameters and two additional measures of the radial-size increase and magnetic-field-strength decrease with heliocentric distance of CMEs based on in-situ measurements are also calculated. For most of the events, the CME radial size estimated by remote observations is inconsistent with the in-situ estimates. We further statistically analyze the correlations of these expansion parameters estimated using remote and in-situ observations, and discuss the potential reasons for the inconsistencies and their implications for the CME space weather forecasting.

The origin and evolution of magnetic fields of neutron stars from birth has long been a source of debate. Here, motivated by recent simulations of the Hall cascade with magnetic helicity, we invoke a model where the large-scale magnetic field of neutron stars grows as a product of small-scale turbulence through an inverse cascade. We apply this model to a simulated population of neutron stars at birth and show how this model can account for the evolution of such objects across the $P\dot{P}$ diagram, explaining both pulsar and magnetar observations. Under the assumption that small-scale turbulence is responsible for large-scale magnetic fields, we place a lower limit on the spherical harmonic degree of the energy-carrying magnetic eddies of $\approx 40$. Our results favor the presence of a highly resistive pasta layer at the base of the neutron star crust. We further discuss the implications of this paradigm on direct observables, such as the nominal age and braking index of pulsars.

The chemical equilibrium distribution of 69 elements between gas and melt is modeled for bulk silicate Earth (BSE) material from 1000 - 4500 K and 1e-6 to 100 bar. The BSE melt is modeled as a non-ideal solution and the effects of different activity coefficients and ideal solution are studied. Results include 50% condensation temperatures, major gases of each element, and oxygen fugacity (fO2) of dry and wet BSE material. The dry BSE model excludes H, C, N, F, Cl, Br, I, S, Se, Te. The wet BSE model includes H and the other volatiles. Key conclusions are much higher condensation temperatures in silicate vapor than in solar composition gas at the same total P, a different condensation sequence in silicate vapor than in solar composition gas, good agreement between different activity coefficient models except for the alkalis, agreement, where overlap exists, with prior published work, condensation of Re, Mo, W, Ru, Os oxides instead of metals, a stability field for Ni-rich metal as reported by Lock et al. (2018), agreement between ideal solution (from this work and from Lock et al. 2018) and real solution condensation temperatures for elements with minor deviations from ideality in the oxide melt, similar 50% condensation temperatures, within a few degrees, in the dry and wet BSE models for the major elements Al, Ca, Fe, Mg, Si, and the minor elements Co, Cr, Li, Mn, Ti, V, and much lower 50 percent condensation temperatures for elements such as B, Cu, K, Na, Pb, Rb, which form halide, hydroxide, sulfide, selenide, telluride and oxyhalide gases. The latter results are preliminary because the poorly known solubilities and activities of volatile elements in silicate melts must be considered for the correct equilibrium distribution, condensation temperatures and mass balance of F, Cl, Br, I, H, S, Se and Te bearing species between melt and vapor (abridged).

Bioavailable nitrogen is thought to be a requirement for the origin and sustenance of life. Before the onset of biological nitrogen fixation, abiotic pathways to fix atmospheric N2 must have been prominent to provide bioavailable nitrogen to Earth's earliest ecosystems. Lightning has been shown to produce fixed nitrogen as nitrite and nitrate in both modern atmospheres dominated by N2 and O2 and atmospheres dominated by N2 and CO2 analogous to the Archaean Earth. However, a better understanding of the isotopic fingerprints of lightning-generated fixed nitrogen is needed to assess the role of this process on the early Earth. Here, we present results from spark discharge experiments in N2-CO2 and N2-O2 gas mixtures. Our experiments suggest that lightning-driven nitrogen fixation may have been similarly efficient in the Archaean atmosphere, compared to modern times. Measurements of the isotopic ratio {\delta}15N of the discharge-produced nitrite and nitrate in solution show very low values of -6 to -15 permil after equilibration with the gas phase with a calculated endmember composition of -17 permil. These results are much lower than most {\delta}15N values documented from the sedimentary rock record, which supports the development of biological nitrogen fixation earlier than 3.2 Ga. However, some Paleoarchean records (3.7 Ga) may be consistent with lightning-derived nitrogen input, highlighting the potential role of this process for the earliest ecosystems.

We propose a novel framework where baryon asymmetry can arise due to forbidden decay of dark matter (DM) enabled by finite temperature effects in the early universe. In order to implement it in a realistic setup, we consider the DM to be a singlet Dirac fermion which acquires a dark asymmetry from a scalar field $\Phi$ via Affleck-Dine mechanism. Due to finite-temperature effects, DM can decay in the early universe into leptons and a second Higgs doublet thereby transferring a part of the dark asymmetry into lepton asymmetry with the latter getting converted into baryon asymmetry subsequently via electroweak sphalerons. DM becomes stable below a critical temperature leading to a stable relic. While the scalar field $\Phi$ can play the role of inflaton with specific predictions for inflationary parameters, the setup also remains verifiable via astrophysical as well as laboratory based observations.

The density of relic neutrinos is expected to be enhanced due to clustering in our local neighbourhood at Earth. We introduce a novel analytical technique to calculate the neutrino overdensity, based on kinetic field theory. Kinetic field theory is a particle-based theory for cosmic structure formation and in this work we apply it for the first time to massive neutrinos. The gravitational interaction is expanded in a perturbation series and we take into account the first-order contribution to the local density of relic neutrinos. For neutrino masses that are consistent with cosmological neutrino mass bounds, our results are in excellent agreement with state-of-the-art calculations.

The recent studies of electroluminescence (EL) properties in two-phase argon detectors for dark matter searches have revealed the presence of unusual delayed pulses in the EL signal in the form of two slow components with time constants of about 5 and 50 $\mu$s. These components were shown to be present in the charge signal itself, which clearly indicates that drifting electrons are temporarily trapped on two states of metastable negative argon ions which have never been observed before. In this work, using the pressure dependence of the ratio of slow component contributions measured in experiment, it is shown that these states are those of two types of metastable negative molecular ions, $\mathrm{Ar}_2^{*-}(b \ ^4\Sigma_u^-)$ and $\mathrm{Ar}_2^{*-}(a \ ^4\Sigma_g^+)$ for the higher and lower energy level respectively.

We investigate the probability distribution of the spins of primordial black holes (PBHs) formed in the universe dominated by a perfect fluid with the linear equation of state $p=w\rho$, where $p$ and $\rho$ are the pressure and energy density of the fluid, respectively. We particularly focus on the parameter region $0<w\leq 1/3$ since the larger value of the spin is expected for the softer equation of state than that of the radiation fluid ($w=1/3$). The angular momentum inside the collapsing region is estimated based on the linear perturbation equation at the turn-around time which we define as the time when the linear velocity perturbation in the conformal Newtonian gauge takes the minimum value. The probability distribution is derived based on the peak theory with the Gaussian curvature perturbation. We find that the root mean square of the non-dimensional Kerr parameter $\sqrt{\langle a_{*}^2\rangle}$ is approximately proportional to $(M/M_{H})^{-1/3}(6w)^{-(1+2w)/(1+3w)}$, where $M$ and $M_{H}$ are the mass of the PBH and the horizon mass at the horizon entry, respectively. Therefore the typical value of the spin parameter decreases with the value of $w$. We also evaluate the mass and spin distribution $P(a_{*}, M)$, taking account of the critical phenomena. We find that, while the spin is mostly distributed in the range of $10^{-3.9}\leq a_{*}\leq 10^{-1.8}$ for the radiation-dominated universe, the peak of the spin distribution is shifted to the larger range $10^{-3.0}\leq a_{*}\leq 10^{-0.7}$ for $w=10^{-3}$.

The Cassini spacecraft's Grand Finale flybys through Saturn's ionosphere provided unprecedented insight into the composition and dynamics of the gas giant's upper atmosphere and a novel and complex spacecraft-plasma interaction. In this article, we further study Cassini's interaction with Saturn's ionosphere using three dimensional Particle-in-Cell simulations. We focus on understanding how electrons and ions, emitted from spacecraft surfaces due to the high-velocity impact of atmospheric water molecules, could have affected the spacecraft potential and low-energy plasma measurements. The simulations show emitted electrons extend upstream along the magnetic field and, for sufficiently high emission rates, charge the spacecraft to positive potentials. The lack of accurate emission rates and characteristics, however, makes differentiation between the prominence of secondary electron emission and ionospheric charged dust populations, which induce similar charging effects, difficult for Cassini. These results provide further context for Cassini's final measurements and highlight the need for future laboratory studies to support high-velocity flyby missions through planetary and cometary ionospheres.

We study graviton-photon conversion in magnetosphere of a pulsar and explore the possibility of detecting high frequency gravitational waves with pulsar observations. It is shown that conversion of one polarization mode of photons can be enhanced significantly due to strong magnetic fields around a pulsar. We also constrain stochastic gravitational waves in frequency range of $10^{8}-10^{9}\,$Hz and $10^{13}-10^{27}\,$Hz by using data of observations of the Crab pulsar and the Geminga pulsar. Our method widely fills the gap among existing high frequency gravitational wave experiments and boosts the frequency frontier in gravitational wave observations.

Rotations of axion fields in the early universe can produce dark matter and the matter-antimatter asymmetry of the universe. We point out that the rotation can generate an observable amount of a stochastic gravitational-wave (GW) background. It can be doubly enhanced in a class of models in which the equation of state of the rotations rapidly changes from a non-relativistic matter-like one to a kination-like one by 1) the so-called Poltergeist mechanism and 2) slower redshift of GWs compared to the axion kination fluid. In supersymmetric UV completion, future GW observations can probe the supersymmetry-breaking scale up to $10^7\,$GeV even if the axion does not directly couple to the Standard Model fields.

Post-inflationary evolution and (re)heating of the viable inflationary model, the $R^2$ one, is made more realistic by including the leptogenesis scenario into it. For this purpose, right-handed Majorana neutrinos with a large mass are added to the matter sector of the Standard Model to explain the neutrino oscillation experiments and the baryon asymmetry of the Universe. We have found parameters that characterize this model: non-minimal coupling of the Higgs field $\xi$ and the mass of the right-handed Majorana neutrino $M_{N_\alpha}$. We have analyzed the effect of these parameters on the reheating process and the resultant physical quantities: spectral indices and baryon asymmetry.