Papers for Friday, Mar 15 2024

A major focus of the planetary science and astrobiology community is the understanding of planetary habitability, including the myriad factors that control the evolution and sustainability of temperate surface environments such as that of Earth. The few substantial terrestrial planetary atmospheres within the Solar System serve as a critical resource in studying these habitability factors, from which models can be constructed for application to extrasolar planets. The recent Astronomy and Astrophysics and Planetary Science and Astrobiology Decadal Surveys both emphasise the need for an improved understanding of planetary habitability as an essential goal within the context of astrobiology. The divergence in climate evolution of Venus and Earth provides a major, accessible basis for understanding how the habitability of large rocky worlds evolves with time and what conditions limit the boundaries of habitability. Here, we argue that Venus can be considered an "anchor point" for understanding planetary habitability within the context of terrestrial planet evolution. We discuss the major factors that have influenced the respective evolutionary pathways of Venus and Earth, how these factors might be weighted in their overall influence, and the measurements that will shed further light on their impacts of these worlds' histories. We further discuss the importance of Venus with respect to both of the recent decadal surveys, and how these community consensus reports can help shape the exploration of Venus in the coming decades.

JAX is widely used in machine learning and scientific computing, the latter of which often relies on existing high-performance code that we would ideally like to incorporate into JAX. Reimplementing the existing code in JAX is often impractical and the existing interface in JAX for binding custom code requires deep knowledge of JAX and its C++ backend. The goal of JAXbind is to drastically reduce the effort required to bind custom functions implemented in other programming languages to JAX. Specifically, JAXbind provides an easy-to-use Python interface for defining custom so-called JAX primitives that support arbitrary JAX transformations.

This work analyses the present-day mass function (PDMF) of 93~star clusters utilizing Gaia DR3 data, with membership determined by the StarGo machine learning algorithm. The impact of unresolved binary systems on mass estimation is rigorously assessed, adopting three mass ratio profiles for correction. The PDMF is characterized by the power-law index, $\alpha$, derived through a robust maximum likelihood method that avoids biases associated with data binning. The value of $\alpha$ for stars between the completeness limited mass of Gaia with a mean 0.3 $M_\odot$ for our cluster samples and 2 $M_\odot$, exhibits stability for clusters younger than 200 Myr, decreasing for older clusters, particularly when considering stars within the half-mass radius. The PDMF of these star clusters is consistent with a dynamically evolved Kroupa IMF via the loss of low-mass stars. Cluster morphology shows a correlation with $\alpha$, as $\alpha$ values exhibit a decreasing trend from filamentary to tidal-tail clusters, mirroring the sequence of increasing cluster age. The dependence of $\alpha$ on total cluster mass is weak, with a subtle increase for higher-mass clusters, especially outside the half-mass radius. We do not observe a correlation between $\alpha$ and the mean metallicity of the clusters. Younger clusters have lower metallicity compared to their older counterparts, which indicates that the older clusters might have migrated to the solar neighbourhood from the inner disk. A comparison with numerical models incorporating a black hole population suggests the need for observations of distant, older, massive open clusters to determine whether or not they contain black holes.

We present PAPERCLIP (Proposal Abstracts Provide an Effective Representation for Contrastive Language-Image Pre-training), a method which associates astronomical observations imaged by telescopes with natural language using a neural network model. The model is fine-tuned from a pre-trained Contrastive Language-Image Pre-training (CLIP) model using successful observing proposal abstracts and corresponding downstream observations, with the abstracts optionally summarized via guided generation using large language models (LLMs). Using observations from the Hubble Space Telescope (HST) as an example, we show that the fine-tuned model embodies a meaningful joint representation between observations and natural language through tests targeting image retrieval (i.e., finding the most relevant observations using natural language queries) and description retrieval (i.e., querying for astrophysical object classes and use cases most relevant to a given observation). Our study demonstrates the potential for using generalist foundation models rather than task-specific models for interacting with astronomical data by leveraging text as an interface.

Recent observations have demonstrated that very-low mass stars and brown dwarfs are capable of sustaining strong magnetic fields despite their cool and neutral atmospheres. These kG field strengths are inferred based on strong highly circularly polarized GHz radio emission, a consequence of the electron cyclotron maser instability. Crucially, these observations imply the existence of energetic non-thermal electron populations, associated with strong current systems, as are found in the auroral regions of the magnetized planets of the Solar System. Intense auroral electron precipitation will lead to electron collisions with the H$_{2}$ gas that should ultimately generate the ion H$_{3}^{+}$. With this motivation, we targeted a sample of ultracool dwarfs, known to exhibit signatures associated with aurorae, in search of the K-band emission features of H$_{3}^{+}$ using the Keck telescopes on Mauna Kea. From our sample of 9 objects, we found no clear indication of H$_{3}^{+}$ emission features in our low-medium resolution spectra (R$\sim$3600). We also modeled the impact of an auroral electron beam on a brown dwarf atmosphere, determining the depth at which energetic beams deposit their energy and drive particle impact ionization. We find that the H$_{3}^{+}$ non-detections can be explained by electron beams of typical energies $\gtrsim$2-10~keV, which penetrate deeply enough that any H$_{3}^{+}$ produced is chemically destroyed before radiating energy through its infrared transitions. Strong electron beams could further explain the lack of UV detections, and suggest that most or nearly all of the precipitating auroral energy must ultimately emerge as thermal emissions deep in brown dwarf atmospheres.

The unprecedented infrared spectroscopic capabilities of JWST have provided high-quality interstellar medium (ISM) metallicity measurements and enabled characterization of the gas-phase mass-metallicity relation (MZR) for galaxies at $z \gtrsim 5$ for the first time. We analyze the gas-phase MZR and its evolution in a high-redshift suite of FIRE-2 cosmological zoom-in simulations at $z=5-12$ and for stellar masses $M_* \sim 10^6-10^{10} \rm{M}_\odot$. These simulations implement a multi-channel stellar feedback model and produce broadly realistic galaxy properties, including when evolved to $z=0$. The simulations predict very weak redshift evolution of the MZR over the redshift range studied, with the normalization of the MZR increasing by less than $0.01$ dex as redshift decreases from $z = 12$ to $z=5$. The median MZR in the simulations is well-approximated as a constant power-law relation across this redshift range given by $\log(Z/Z_\odot) = 0.37\log(M_*/\rm{M}_\odot) - 4.3$. We find good agreement between our best-fit model and recent observations made by JWST at high redshift. The weak evolution of the MZR at $z > 5$ contrasts with the evolution at $z \lesssim 3$, where increasing normalization of the MZR with decreasing redshift is observed and predicted by most models. The FIRE-2 simulations predict increasing scatter in the gas-phase MZR with decreasing stellar mass, in qualitative agreement with some observations.

The scatter about the mass-metallicity relation (MZR) has a correlation with the star formation rate (SFR) of galaxies. The lack of evidence of evolution in correlated scatter at $z\lesssim2.5$ leads many to refer to the relationship between mass, metallicity, and SFR as the Fundamental Metallicity Relation (FMR). Yet, recent high-redshift ($z>3$) JWST{} observations have challenged the fundamental (i.e., redshift-invariant) nature of the FMR. In this work, we show that the cosmological simulations Illustris, IllustrisTNG, and EAGLE all predict MZRs that exhibit scatter with a secondary dependence on SFR up to $z=8$. We introduce the concept of a ``strong'' FMR, where the strength of correlated scatter does not evolve with time, and a ``weak'' FMR, where there is some time evolution. We find that each simulation analysed has a weak FMR -- there is non-negligible evolution in the strength of the correlation with SFR. Furthermore, we show that the scatter is reduced an additional $\sim$10-40\% at $z\gtrsim3$ when using a weak FMR, compared to assuming a strong FMR. These results highlight the importance of avoiding coarse redshift binning when assessing the FMR.

Despite a burgeoning set of ultracool dwarf ($\leq$M7) radio detections, their radio emissions remain enigmatic. Open questions include the plasma source and acceleration mechanisms for the non-auroral "quiescent" component of these objects' radio emissions, which can trace Jovian synchrotron radiation belt analogs. Ultracool dwarf binary systems can provide test beds for examining the underlying physics for these plasma processes. We extend a recently developed occurrence rate calculation framework to compare the quiescent radio occurrence rate of binary systems to single objects. This generalized and semi-analytical framework can be applied to any set of astrophysical objects conceptualized as unresolved binary systems with approximately steady-state emission or absorption. We combine data available in the literature to create samples of 179 single ultracool dwarfs (82 M dwarfs, 74 L dwarfs, and 23 T/Y dwarfs) and 27 binary ultracool dwarf systems. Using these samples, we show that quiescent radio emissions occur in $54^{+11}_{-11}$ - $65^{+11}_{-12}$ per cent of binaries where both components are ultracool dwarfs, depending on priors. We also show that binarity enhances the ultracool dwarf quiescent radio occurrence rate relative to their single counterparts. Finally, we discuss potential implications for the underlying drivers of ultracool dwarf quiescent radio emissions, including possible plasma sources.

Future black hole (BH) imaging observations with better sensitivity and angular resolution are expected to resolve finer features corresponding to higher-order images of both hotspots that are produced in the accretion flow, as well as of the entire emitting region. In spherically symmetric spacetimes, the image order is determined by the maximum number of half-loops executed around the BH by the photons that form it. Due to the additional half-loop, consecutive-order images arrive after a delay time of approximately $\pi$ times the BH shadow radius. Furthermore, the deviation of the diameters from that of the shadow, widths, and flux-densities of consecutive-order images are exponentially demagnified by the lensing Lyapunov exponent, a characteristic of the spacetime. We compare the exact time delay between the appearance of the zeroth and first-order images of a hotspot to our best analytic estimate and find an error $\lesssim 50\%$ for hotspot locations within $\approx 10M$ from a Schwarzschild BH of mass $M$. We also explore the possibilities for inferring kinetic properties of hotspots, and also our inclination, from future BH movies. Furthermore, since the targets of such observations host geometrically-thick accretion flows (also jets), we obtain simple theoretical estimates of the variation in the diameters and widths of their first-order images, for varying disk scale-heights. We find realistically that the deviation of the former from that of the shadow is $\lesssim 30\%$ and that the latter remain $\lesssim 1.3M$. Finally, we estimate the error in recovering the lensing exponent, when using the first and second-order images, to be $\lesssim 20\%$. This provides further evidence that future observations could yield new and independent estimates of the shadow size and, in principle, of the lensing exponent, allowing us to robustly learn about the spacetimes of astrophysical BHs.

Exoplanets in the ultra-hot Jupiter regime provide an excellent laboratory for testing the impact of stellar irradiation on the dynamics and chemical composition of gas giant atmospheres. In this study, we observed two transits of the ultra-hot Jupiter WASP-189 b with MAROON-X/Gemini-North to probe its high-altitude atmospheric layers, using strong absorption lines. We derived posterior probability distributions for the planetary and stellar parameters by calculating the stellar spectrum behind the planet at every orbital phase during the transit. This was used to correct the Rossiter-McLaughlin imprint on the transmission spectra. Using differential transmission spectroscopy, we detect strong absorption lines of Ca+, Ba+, Na, H$\alpha$, Mg, Fe, and Fe+, providing an unprecedented and detailed view of the atmospheric chemical composition. Ca+ absorption is particularly well suited for analysis through time-resolved narrow-band spectroscopy, owing to its transition lines formed in high-altitude layers. The spectral absorption lines show no significant blueshifts that would indicate high-altitude day-to-night winds, and further analysis is needed to investigate the implications for atmospheric dynamics. These high signal-to-noise observations provide a benchmark data set for testing high-resolution retrievals and the assumptions of atmospheric models. We also simulate observations of WASP-189 b with ANDES/ELT, and show that ANDES will be highly sensitive to the individual absorption lines of a myriad of elements and molecules, including TiO and CO.

Located at less than 2pc away, Luhman16AB (WISE.J104915.57-531906.1) is the closest pair of brown dwarfs and third closest `stellar' system to Earth. An exoplanet candidate in the Luhman16 binary system was reported in 2017 based on a weak astrometric signature in the analysis of 12 HST epochs. An additional epoch collected in 2018 and re-analysis of the data with more advanced methods further increased the significance level of the candidate, consistent with a Neptune-mass exoplanet orbiting one of the Luhman16 brown dwarf components. We report the joint analysis of these previous data together with two new astrometric HST epochs we obtained to confirm or disprove this astrometric signature. Our new analysis rules out presence of a planet orbiting one component of the Luhman16AB system for masses M > 1.5 M_Nep (Neptune masses) and periods between 400 and 5000 days. However, the presence of third bodies with masses M < 3 M_Nep and periods between 2 and 400 days (~1.1yrs) can not be excluded. Our measurements make significant improvements to the characterization of this sub-stellar binary, including its mass-ratio 0.8305+/-0.0006, individual component masses 35.4+/-0.2 M_Jup and 29.4+/-0.2 M_Jup (Jupiter masses), and parallax distance 1.9960pc +/- 50AU. Comparison of the masses and luminosities of Luhman16AB to several evolutionary models shows persistent discrepancies in the ages of the two components, but strengthens the case that this system is a member of the 510+/-95 Myr Oceanus Moving Group.

The ultra luminous infrared galaxy (ULIRG) F13451+1232 is an excellent example of a galaxy merger in the early stages of active galactic nucleus (AGN) activity, a phase in which AGN-driven outflows are expected to be particularly important. However, previous observations have determined that the mass outflow rates of the warm ionised and neutral gas phases in F13451+1232 are relatively modest, and there has been no robust detection of molecular outflows. Using high spatial resolution ALMA CO(1-0) observations, we detect a kiloparsec-scale circumnuclear disk, as well as extended ($r\sim440$ pc), intermediate-velocity (300<|$v$|<400 km s$^{-1}$) cold molecular gas emission that cannot be explained by rotational disk motions. If interpreted as AGN-driven outflows, the mass outflow rates associated with this intermediate-velocity gas are relatively modest ($\dot{M}_\mathrm{out}=22$-$27$ M$_\odot$ yr$^{-1}$); however, we also detect a compact ($r_\mathrm{out}$<120 pc), high velocity (400<$v$<680 km s$^{-1}$) cold molecular outflow near the primary nucleus of F13451+1232, which carries an order of magnitude more mass ($\dot{M}_\mathrm{out}\sim230$ M$_\odot$ yr$^{-1}$) than (and several times the kinetic power of) the previously-detected warmer phases. Moreover, the similar spatial scales of this compact outflow and the radio structure indicate that it is likely accelerated by the small-scale ($r\sim130$ pc) AGN jet in the primary nucleus of F13451+1232. Considering the compactness of the nuclear outflow and intermediate-velocity non-rotating gas that we detect, we argue that high spatial-resolution observations are necessary to properly quantify the properties of AGN-driven outflows and their impacts on host galaxies.

The accretion of tidally disrupted planetary bodies is the current consensus model for the presence of photospheric metals commonly detected in white dwarfs. While most dynamical studies have considered a single star and associated planetary instabilities, several investigations have instead considered the influence of widely-bound stellar companions as potential drivers of white dwarf pollution. This study examines the prevalence of wide binaries among polluted white dwarfs using Gaia DR3 astrometry, where three samples are investigated: 71 DAZ stars with metals detected in the ultraviolet using Hubble, and two groups of DZ stars identified via SDSS spectroscopy, comprised of 116 warmer and 101 cooler sources. Each sample was searched for spatially-resolved, co-moving companions, and compared to the same analysis of thousands of field white dwarfs within overlapping regions of the Gaia HR diagram. The wide binary fraction of the DAZ sample is $10.6_{-3.2}^{+3.9}$ per cent, and within $1 \sigma$ of the corresponding field. However, the search yields wide binary fractions of less than 1.8 per cent for the two independent DZ star catalogues, which are each distinct from their fields by more than $3 \sigma$. Both sets of results support that pollution in white dwarfs is not the result of stellar companions, and the delivery of metals to white dwarf surfaces is caused by major planets. The discrepancy between the DAZ and DZ star wide binary fractions cannot be caused by white dwarf spectral evolution, suggesting these two populations may have distinct planetary architectures.

We investigate the effect of cutting-edge circumbinary disk (CBD) evolution models on massive black hole binary (MBHB) populations and the gravitational wave background (GWB). We show that CBD-driven evolution leaves a tell-tale signature in MBHB populations, by driving binaries towards an equilibrium eccentricity that depends on binary mass ratio. We find high orbital eccentricities ($e_{\rm b} \sim 0.5$) as MBHBs enter multi-messenger observable frequency bands. The CBD-induced eccentricity distribution of MBHB populations in observable bands is independent of the initial eccentricity distribution at binary formation, erasing any memory of eccentricities induced in the large-scale dynamics of merging galaxies. Our results suggest that eccentric MBHBs are the rule rather than the exception in upcoming transient surveys, provided that CBDs regularly form in MBHB systems. We show that the GWB amplitude is sensitive to CBD-driven preferential accretion onto the secondary, resulting in an increase in GWB amplitude $A_{\rm yr^{-1}}$ by over 100\% with just 10\% Eddington accretion. As we self consistently allow for binary hardening and softening, we show that CBD-driven orbital expansion does not diminish the GWB amplitude, and instead increases the amplitude by a small amount. We further present detection rates and population statistics of MBHBs with $M_{\rm b} \gtrsim 10^6 \, M_{\odot}$ in LISA, showing that most binaries have equal mass ratios and can retain residual eccentricities up to $e_{\rm b} \sim 10^{-3}$ due to CBD-driven evolution.

Combining the public JWST/NIRCam imaging programs CEERS, PRIMER and JADES, spanning a total area of $\sim$500 arcmin$^2$, we obtain a sample of $>$30,000 galaxies at $z\sim4-9$ that allows us to perform a complete, rest-optical selected census of the galaxy population at $z>3$. Comparing the stellar mass $M_*$ and the UV-slope $\beta$ distributions between JWST- and HST-selected samples, we generally find very good agreement and no significant biases. Nevertheless, JWST enables us to probe a small population of UV-red galaxies that was missing from previous HST-based LBG samples. We measure galaxy stellar mass functions (SMFs) at $z\sim4-9$ and show that they are broadly consistent with existing literature results. However, UV-red galaxies dominate the high-mass end of the SMF at least out to $z\sim6$. In particular the most massive galaxies typically show very red colors between $\lambda_{obs}\sim1.5\mu$m and $\sim4.5\mu$m, and thus JWST's unprecedented resolution and sensitivity at these wavelengths yields more accurate constraints on their abundance and masses. The implied redshift evolution of the high-mass end of the SMF suggests a rapid build-up of massive dust-obscured as well as quiescent galaxies from $z\sim6$ to $z\sim4$ as well as an enhanced efficiency of star formation towards earlier times ($z\gtrsim6$). We find the SMFs to be steep over the entire redshift range, and slightly steepening with redshift from $z\sim 4-6$, reaching values of $\approx-2$ at $z\gtrsim6$. Finally, we show that the galaxy mass density grows by a factor $\sim20\times$ in the $\sim1$ Gyr of cosmic time from $z\sim9$ to $z\sim4$. Our results emphasize the importance of rest-frame optically-selected samples in inferring accurate distributions of physical properties and studying the mass build-up of galaxies in the first 1.5 Gyr of cosmic history.

The connection between outer gas giants and inner super-Earths reflects their formation and evolutionary histories. Past work exploring this link has suggested a tentative positive correlation between these two populations, but these studies have been limited by small sample sizes and in some cases sample biases. Here we take a new look at this connection with a sample of 184 super-Earth systems with publicly available radial velocity data and fully resolved outer gas giants. We calculate the frequency of outer gas giants (GG) in super-Earth (SE) systems, dividing our sample into metal-rich ([Fe/H] $>$ 0) and metal-poor ([Fe/H]$\leq$0) hosts. We find P(GG$|$SE, [Fe/H]$>$0) = 28.0$^{+4.9}_{-4.6}\%$ and P(GG$|$SE, [Fe/H]$\leq$0) = 4.5$^{+2.6}_{-1.9}\%$. Comparing these conditional occurrence rates to field giant occurrence rates from Wittenmyer et al. 2020, we show that there is a distinct positive correlation between inner super-Earths and outer gas giants for metal-rich host stars at the 2.5$\sigma$ level, but this correlation disappears for metal-poor systems. We further find that, around metal-rich stars, the GG/SE correlation enhances for systems with giants that are more distant (beyond 3 AU) and/or more eccentric ($e > 0.2$), while gas giant multiplicity does not appear to affect the level of correlation. Such trends again disappear around metal-poor stars with the exception of systems of multiple giants in which we observe an anti-correlation. Our findings highlight the critical role metallicity (disk solid budget) plays in shaping the overall planetary architecture.

We investigate mass ejection from accretion disks formed during the collapse of rapidly-rotating Wolf-Rayet stars. The neutrino-cooled, black hole (BH) accretion disk system that forms at the center of the star -- and the ensuing outflows -- provide the conditions for these systems to be candidate $r$-process element production sites and potential progenitors of broad-lined Type Ic (Ic-BL) supernovae. Here we present global, long-term axisymmetric hydrodynamic simulations of collapsar disks that include angular momentum transport through shear viscosity, neutrino emission and absorption, a 19-isotope nuclear reaction network and nuclear statistical equilibrium solver, a pseudo-Newtonian BH with mass and spin modified by accreted matter, and self-gravity. Starting from a stellar profile collapsed in spherical symmetry, our models capture disk formation self-consistently, and are evolved until after the shock wave -- driven by disk winds -- reaches the surface of the star. None of our models achieve sufficient neutronization to eject significant amounts of $r$-process elements (detailed nucleosynthesis calculations will follow in a companion paper). Sufficient $^{56}$Ni is produced to power a typical type Ic-BL supernova light curve, but the average asymptotic velocity is a factor $\sim 2-3$ times too slow to account for the typical line widths in type Ic-BL supernova spectra. The gap in neutrino emission between BH formation and shocked disk formation, and the magnitude of the subsequent peak in emission, would be observable diagnostics of the internal conditions of the progenitor in a galactic collapsar. Periodic oscillations of the shocked disk prior to its expansion are also a potential observable through their impact on the the neutrino and gravitational wave signals.

Four years after the discovery of a unique DAQ white dwarf with a hydrogen-dominated and carbon-rich atmosphere, we report the discovery of four new DAQ white dwarfs, including two that were not recognized properly in the literature. We find all five DAQs in a relatively narrow mass and temperature range of $M=1.14-1.19~M_{\odot}$ and $T_{\rm eff}=13,000-17,000$ K. In addition, at least two show photometric variations due to rapid rotation with $\approx10$ min periods. All five are also kinematically old, but appear photometrically young with estimated cooling ages of about 1 Gyr based on standard cooling tracks, and their masses are roughly twice the mass of the most common white dwarfs in the solar neighborhood. These characteristics are smoking gun signatures of white dwarf merger remnants. Comparing the DAQ sample with warm DQ white dwarfs, we demonstrate that there is a range of hydrogen abundances among the warm DQ population, and the distinction between DAQ and warm DQ white dwarfs is superficial. We discuss the potential evolutionary channels for the emergence of the DAQ subclass, and suggest that DAQ white dwarfs are trapped on the crystallization sequence, and may remain there for a significant fraction of the Hubble time.

The flux ratios of strongly lensed quasars have previously been used to infer the properties of dark matter. In these analyses it is crucial to separate the effect of the main lensing galaxy and the low-mass dark matter halo population. In this work, we investigate flux-ratio perturbations resulting from general third- and fourth-order multipole perturbations to the main lensing galaxy's mass profile. We simulate four lens systems, each with a different lensing configuration, without multipoles. The simulated flux ratios are perturbed by 10-40 per cent by a population of low-mass haloes consistent with CDM and, in one case, also a satellite galaxy. This level of perturbation is comparable to the magnitude of flux-ratio anomalies in real data that has been previously analyzed. We then attempt to fit the simulated systems using multipoles instead of low-mass haloes. We find that multipoles with amplitudes of 0.01 or less can produce flux-ratio perturbations in excess of 40 per cent. In all cases, third- or fourth-order multipoles can individually reduce the magnitude of, if not eliminate, flux-ratio anomalies. When both multipole orders are jointly included, all simulated flux ratios can be fit to within the observational uncertainty. Our results indicate that low-mass haloes and multipoles are highly degenerate when modelling quadruply-imaged quasars based just on image positions and flux ratios. In the presence of this degeneracy, flux-ratio anomalies in lensed quasars alone cannot be used to place strong constraints on the properties of dark matter without additional information that can inform our priors.

Nitrogen is produced by CNO-cycling in massive stars, and can be ejected in significant amounts in supernova explosions. While in H-rich SNe, its [\ion{N}{II}] 6548, 6583 emission becomes obscured by strong H$\alpha$, in explosions of He stars, this nitrogen emission becomes more visible. We here explore the formation of this line, using the \texttt{SUMO} code to compute spectra for a grid of 1D models with parameterized mixing informed from new 2D simulations. Because the mass fraction of nitrogen in the ejecta decreases with larger He core masses, as more of the He/N zone gets processed by shell helium burning and is lost to winds, the [\ion{N}{II}] luminosity relative to the overall optical flux probes the He core mass. By comparing to large samples of data, we find that low-mass He cores ($M_{\rm preSN}\lesssim\ 3\ M_\odot$) are exclusively associated with Type IIb SNe, with the exception of Type Ib SN 2007Y. Seeing no strong nitrogen emission in other Type Ib SNe, the implication is either an origin from low-mass stars with the He/N layer (but not the He/C) layer peeled away, or from higher-mass He cores. We also see no clear nitrogen emission in Type Ic SNe. We discuss the diagnostic potential of this new line metric, and also dependencies on mass-loss-rate and metallicity.

We present early results from our program of ALMA Band 6 (1.3mm) molecular line mapping of a sample of nearby, well-studied examples of high-excitation, bipolar/pinched-waist and molecule-rich planetary nebulae (Hubble 5 and NGC 2440, 2818, 2899, 6302, and 6445). We have mapped these planetary nebulae (PNe) in isotopologues of CO as well as various molecular line tracers of high-energy irradiation, such as HCN, CN, HNC, and HCO+, with the complementary goals of establishing nebular kinematics as well as the zones of UV-heated and X-ray-ionized molecular gas within each nebula. The resulting high-resolution ALMA molecular emission-line maps reveal the regions of high-excitation bipolar PNe in which molecular gas, presumably ejected during asymptotic giant branch stages of the PN progenitor stars, survives and evolves chemically. We present a summary of molecular species detected to date in the sample nebulae, and we use example results for one PN (NGC 6455) to demonstrate the power of the ALMA data in revealing the structures, kinematics, and compositions of the equatorial molecular tori that are a common feature of the sample objects.

We derive the metallicity profile of the Milky Way low-$\alpha$ disc population from 2 to 20 kpc from the Galactic centre in 1 Gyr age bins using the astroNN catalogue, and show that it is highly structured, with a plateau between 4 and 7 kpc and a break at 10-12 kpc. We argue that these features result from the two main bar resonances, the corotation and the Outer Lindblad Resonance (OLR), respectively. We show that the break in the metallicity profile is most visible in stars having 7-8 Gyr, reaching an amplitude of about 0.4 dex, and is the signpost of the position of the bar OLR. The bar formation was accompanied by an episode of radial migration triggered by its slowing down and is responsible for spreading old metal-rich stars up to the OLR. The data show that the slowdown of the bar ended 6-7 Gyr ago. Based on numerical simulations that reproduce well the break observed in the metallicity profile, we argue that this implies that the bar formed in our Galaxy 8-10 Gyr ago. Analysis of the metallicity distribution as a function of radius shows no evidence of significant systematic outward radial migration after this first episode. We argue that the variation of the metallicity dispersion as a function of the guiding radius is dominated by the migration triggered by the bar, but also that the libration of orbits around the bar resonances induces a mixing that may have a significant impact on the observed metallicity dispersion. In contrast, the absence of a break in the metallicity profile of populations younger than about $\sim$6 Gyr and the flattening of the gradient at younger ages is interpreted as evidence that the strength of the bar has decreased, loosening its barrier effect and allowing the gas and metals on both sides of the OLR to mix, erasing the break. Beyond the OLR, stars younger than 7 Gyr show very small metallicity dispersion, suggesting no or limited mixing.

An exciting possibility to constrain dark matter (DM) scenarios is to search for their gravitational imprints on Black Hole (BH) observations. In this paper, we investigate the impact of self-interacting scalar field DM on the shadow radius of a Schwarzschild BH. We implement a self-consistent formulation, paying attention to the enhancement of the DM density due to the BH gravitational influence and the accretion flow onto the BH. First, we calculate the first-order correction to the shadow radius caused by a general DM environment. Then, we apply this perturbative method to the case of self-interacting scalar field DM and derive analytical expressions for the critical impact parameter. We find that self-consistency requirements, involving the lifetime and the mass of the central DM soliton, or the mass and the size of the extended virialized DM halo, ensure that the impact of the DM environment on the shadow radius is below the observational upper bound. This emphasizes the importance of taking into account the self-consistency constraints of the underlying DM scenario, which can strongly limit the range of possible DM density profiles and their impact on the shadow radius.

Identifying black holes is essential for comprehending the development of stars and uncovering novel principles of physics. Gravitational microlensing provides an exceptional opportunity to examine an undetectable population of black holes in the Milky Way. In particular, long-lasting events are likely to be associated with massive lenses, including black holes. We present an analysis of the Gaia18ajz microlensing event, reported by the Gaia Science Alerts system, which has exhibited a long timescale and features indicative of the annual microlensing parallax effect. Our objective is to estimate the parameters of the lens based on the best-fitting model. We utilized photometric data obtained from the Gaia satellite and terrestrial observatories to investigate a variety of microlensing models and calculate the most probable mass and distance to the lens, taking into consideration a Galactic model as a prior. Subsequently, we applied a mass-brightness relation to evaluate the likelihood that the lens is a main sequence star. We also describe DarkLensCode, an open-source routine which computes the distribution of probable lens mass, distance and luminosity employing the Galaxy priors on stellar density and velocity for microlensing events with detected microlensing parallax. We modelled Gaia18ajz event and found its two possible models with most likely Einstein timescale of $t_\mathrm{E}=316^{+36}_{-30}$ days and $t_\mathrm{E}=299^{+25}_{-22}$ days. Applying Galaxy priors for stellar density and motion, we calculated the most probable lens mass of $M_L = 5.6^{+7.5}_{-2.5} M_\odot$ located at $D_S = 1.05^{+0.78}_{-0.60}\,\text{kpc}$ or $M_L = 12.0^{+14.9}_{-5.4} M_\odot$ located at $D_S = 1.18^{+0.82}_{-0.63}\,\text{kpc}$. Our analysis of the blended light suggests that the lens is likely a dark remnant of stellar evolution, rather than a main sequence star.

We numerically simulate the motion of a black hole as it plunges radially through an ultralight dark matter soliton. We investigate the timescale in which dynamical friction reduces the kinetic energy of the black hole to a minimum, and consider the sensitivity of this timescale to changes in the ULDM particle mass, the total soliton mass, and the mass of the black hole. We contrast our numerical results with a semi-analytic treatment of dynamical friction, and find that the latter is poorly suited to this scenario. In particular, we find that the back-reaction of the soliton to the presence of the black hole is significant, resulting in oscillations in the coefficient of dynamical friction which cannot be described in the simple semi-analytical framework. Furthermore, we observe a late-time reheating effect, in which a significant amount of kinetic energy is transferred back to the black hole after an initial damping phase. This complicates the discussion of ULDM dynamical friction on the scales relevant to the final parsec problem.

We examine the fundamental plane of 91 Blazars which include FSRQs and BL Lacs with known X-ray luminosity ($L_{R}$), radio luminosity ($L_X$), and black hole mass measurements ($M$) to reflect the relationship between jet and accretion for blazars. The fundamental plane of Blazars are log$L_{R}$=${0.273}_{+0.059}^{-0.059}$log$L_X$+${0.695}_{+0.191}^{-0.191}$log$M$+${25.457}_{+2.728}^{-2.728}$ and log$L_{R}$=${0.190}_{+0.049}^{-0.049}$log$L_X$+${0.475}_{+0.157}^{-0.157}$log$M$+${28.568}_{+2.245}^{-2.245}$ after considering the effect of beam factor. Our results suggest that the jet of blazars has connection with accretion. We set the black hole spin energy as a new variable to correct the black hole mass and explore the effect of black hole spin on the fundamental relationship. We find that the fundamental plane of Blazars is effected by the black hole spin, which is similar to the previous work for AGNs. We additionally examine a new fundamental plane which is based on the black hole spin-mass energy ($M_{spin}$). The new fundamental plane (log$L_{R}$=${0.332}_{+0.081}^{-0.081}$log$L_X$+${0.502}_{+0.091}^{-0.091}$log$M_{spin}$+${22.606}_{+3.346}^{-3.346}$ with R-Square=0.575) shows that $M_{spin}$ has a better correlation coefficient comparing to the $M$ for fundamental plane of Blazars. These results support that the black hole spin should be considered as a important factor for the study of fundamental plane for Blazars. And these may further our understanding of the Blandford-Znajek process in blazars.

From the mid-19th century to the end of the 20th century, photographic plates served as the primary detectors for astronomical observations. Astronomical photographic observations in China began in 1901, and over a century, a total of approximately 30,000 astronomical photographic plates have been captured. These historical plates play an irreplaceable role in conducting long-term, time-domain astronomical research. To preserve and explore these valuable original astronomical observational data, Shanghai Astronomical Observatory has organized the transportation of plates taken at night from various stations across the country to the Sheshan Plate Archive for centralized preservation. For the first time, plate information statistics was performed. On this basis, the plates were cleaned and digitally scanned, and finally digitized images were acquired for 29,314 plates. In this study, using Gaia DR2 as the reference star catalog, astrometric processing has been carried out successfully on 15,696 single-exposure plates, including object extraction, stellar identification, and plate model computation. As a result, for long focal length telescopes, such as the 40cm double-tube refractor telescope and the 1.56m reflector telescope at the Shanghai Astronomical Observatory and the 1m reflector telescope at the Yunnan Astronomical Observatory, the astrometric accuracy obtained for their plates is approximately 0.1" to 0.3". The distribution of astrometric accuracy for medium and short focal length telescopes ranges from 0.3" to 1.0". The relevant data of this batch of plates, including digitized images and stellar catalog of the plates are archived and released by the National Astronomical Data Center. Users can access and download plate data based on keywords such as station, telescope, observation year, and observed celestial coordinates.

Galactic-scale outflows of molecular gas from star-forming galaxies constitute the most direct evidence for regulation of star formation. In the early universe ($ z > 4 $), such outflows have recently been inferred from gravitationally-lensed dusty star-forming galaxies (DSFGs) based on ubiquitous detections of OH absorption extending to more blueshifted velocities than [CII] or CO emission in spatially-integrated spectra. Because these lines are redshifted to sub-mm wavelengths, such measurements require careful corrections for atmospheric absorption lines, and a proper accounting of sometimes large variations in measurement uncertainties over these lines. Taking these factors into consideration, we re-analyze OH and [CII] data taken with ALMA for the five sources where such data is available, of which four were categorised as exhibiting outflows. Based on their spatially-integrated spectra alone, we find statistically significant ($ \geq 3 \sigma $) OH absorption more blueshifted than [CII] emission in only one source. By contrast, searching channel maps for signals diluted below the detection threshold in spatially-integrated spectra, we find evidence for a separate kinematic component in OH absorption in all five sources in the form of: (i) more blueshifted OH absorption than [CII] emission and/or (ii) a component in OH absorption exhibiting a different spatio-kinematic pattern than [CII] emission, the latter presumably tracing gas in a rotating disc. Providing a more complete and accurate assessment of molecular outflows in gravitationally-lensed DSFGs, we suggest methods to better assess the precision of corrections for atmospheric absorption and to more accurately measure the source continuum in future observations.

We identify a fading AGN SDSS J220141.64+115124.3 from the internal Product Launch-11 (MPL-11) in Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey. The central region with a projected radius of $\sim$2.4 kpc is characterized as LINER-like line ratios while the outskirts extended to $\sim$15 kpc show Seyfert-like line ratios. The [OIII]$\lambda$5007 luminosity of the Seyfert regions is a factor of 37 (2) higher than the LINER regions without (with) dust attenuation correction, suggesting that the AGN activity decreases at least $\sim$8 $\times$ 10$^3$ yrs ($\sim$2.4 kpc/light-speed) ago. We model the emission line spectra in the central region with double Gaussian components (a narrow core and a broad wing) and analyze the properties of each component. The narrow core component mostly co-rotates with the stellar disc, whereas the broad wing component with a median of the velocity dispersion $\sim$300 km s$^{-1}$ is related to a wind outflow. The kinematic position angle (PA) of the ionized gas shows a $\sim$20{\deg} twist from the galaxy center to 1.5 effective radius. The median of the PA difference between the gas and stellar components is as large as $\sim$50{\deg} within 0.4 effective radius. The tidal feature in DESI image and star-gas misalignment suggest this galaxy is a merger remnant. Combining all these observational results as well as public available X-ray and MIR luminosities, we confirm this is a fading AGN, the merger process kick-started the central engine to quasar phase which ionized gas composed of tidal debris, and now the activity of the central black hole decreases. The discontinuity in [OIII]$\lambda$5007 flux and EQW maps is due to multiple AGN outbursts triggered by merger remnant gas inflows.

Self-interacting dark matter (SIDM) has the potential to significantly influence galaxy formation in comparison to the cold, collisionless dark matter paradigm (CDM), resulting in observable effects. This study aims to elucidate this influence and to demonstrate that the stellar mass Tully-Fisher relation imposes robust constraints on the parameter space of velocity-dependent SIDM models. We present a new set of cosmological hydrodynamical simulations that include the SIDM scheme from the TangoSIDM project and the SWIFT-EAGLE galaxy formation model. Two cosmological simulations suites were generated: one (Reference model) which yields good agreement with the observed $z=0$ galaxy stellar mass function, galaxy mass-size relation, and stellar-to-halo mass relation; and another (WeakStellarFB model) in which the stellar feedback is less efficient, particularly for Milky Way-like systems. Both galaxy formation models were simulated under four dark matter cosmologies: CDM, SIDM with two different velocity-dependent cross sections, and SIDM with a constant cross section. While SIDM does not modify global galaxy properties such as stellar masses and star formation rates, it does make the galaxies more extended. In Milky Way-like galaxies, where baryons dominate the central gravitational potential, SIDM thermalises, causing dark matter to accumulate in the central regions. This accumulation results in density profiles that are steeper than those produced in CDM from adiabatic contraction. The enhanced dark matter density in the central regions of galaxies causes a deviation in the slope of the Tully-Fisher relation, which significantly diverges from the observational data. In contrast, the Tully-Fisher relation derived from CDM models aligns well with observations.

The study of protoplanetary disks has become increasingly important with the Kepler satellite finding that exoplanets are ubiquitous around stars in our galaxy and the discovery of enormous diversity in planetary system architectures and planet properties. High-resolution near-IR and ALMA images show strong evidence for ongoing planet formation in young disks. The JWST MIRI mid-INfrared Disk Survey (MINDS) aims to (1) investigate the chemical inventory in the terrestrial planet-forming zone across stellar spectral type, (2) follow the gas evolution into the disk dispersal stage, and (3) study the structure of protoplanetary and debris disks in the thermal mid-IR. The MINDS survey will thus build a bridge between the chemical inventory of disks and the properties of exoplanets. The survey comprises 52 targets (Herbig Ae stars, T Tauri stars, very low-mass stars and young debris disks). We primarily obtain MIRI/MRS spectra with high S/N (~100-500) covering the complete wavelength range from 4.9 to 27.9 {\mu}m. For a handful of selected targets we also obtain NIRSpec IFU high resolution spectroscopy (2.87-5.27 {\mu}m). We will search for signposts of planet formation in thermal emission of micron-sized dust - information complementary to near-IR scattered light emission from small dust grains and emission from large dust in the submillimeter wavelength domain. We will also study the spatial structure of disks in three key systems that have shown signposts for planet formation, TW Hya and HD 169142 using the MIRI coronagraph at 15.5 {\mu}m and 10.65 {\mu}m respectively and PDS70 using NIRCam imaging in the 1.87 {\mu}m narrow and the 4.8 {\mu}m medium band filter. ...

The question of the global topology of the Universe (cosmic topology) is still open. In the $\Lambda$CDM concordance model it is assumed that the space of the Universe possesses the trivial topology of $\mathbb{R}^3$ and thus that the Universe has an infinite volume. As an alternative, we study in this paper one of the simplest non-trivial topologies given by a cubic 3-torus describing a universe with a finite volume. To probe cosmic topology, we analyse certain structure properties in the cosmic microwave background (CMB) using Betti Functionals and the Euler Characteristic evaluated on excursions sets, which possess a simple geometrical interpretation. Since the CMB temperature fluctuations $\delta T$ are observed on the sphere $\mathbb{S}^2$ surrounding the observer, there are only three Betti functionals $\beta_k(\nu)$, $k=1,2,3$. Here $\nu=\delta T/\sigma_0$ denotes the temperature threshold normalized by the standard deviation $\sigma_0$ of $\delta T$. Analytic approximations of the Gaussian expectations for the Betti functionals and an exact formula for the Euler characteristic are given. It is shown that the amplitudes of $\beta_0(\nu)$ and $\beta_1(\nu)$ decrease with increasing volume $V=L^3$ of the cubic 3-torus universe. Since the computation of the $\beta_k$'s from observational sky maps is hindered due to the presence of masks, we suggest a method yielding lower and upper bounds for them and apply it to four Planck 2018 sky maps. It is found that the $\beta_k$'s of the Planck maps lie between those of the torus universes with side-lengths $L=2.0$ and $L=3.0$ in units of the Hubble length and above the infinite $\Lambda$CDM case. These results give a further hint that the Universe has a non-trivial topology.

One of the intriguing mechanisms of the Sun is the formation of the bipolar magnetic regions (BMRs) in the solar convection zone which are observed as regions of concentrated magnetic fields of opposite polarity on photosphere. These BMRs are tilted with respect to the equatorial line, which statistically increases with latitude. The thin flux tube model, employing the rise of magnetically buoyant flux loops and their twist by Coriolis force, is a popular paradigm for explaining the formation of tilted BMRs. In this study, we assess the validity of the thin flux tube model by analyzing the tracked BMR data obtained through the Automatic Tracking Algorithm for BMRs (AutoTAB). Our observations reveal that the tracked BMRs exhibit the expected collective behaviors. We find that the polarity separation of BMRs increases over their normalized lifetime, supporting the assumption of a rising flux tube from the CZ. Moreover, we observe an increasing trend of the tilt with the flux of the BMR, suggesting that rising flux tubes associated with lower flux regions are primarily influenced by drag force and Coriolis force, while in higher flux regions, magnetic buoyancy dominates. Furthermore, we observe Joy's law dependence for emerging BMRs from their first detection, indicating that at least a portion of the tilt observed in BMRs can be attributed to the Coriolis force. Notably, lower flux regions exhibit a higher amount of fluctuations associated with their tilt measurement compared to stronger flux regions, suggesting that lower flux regions are more susceptible to turbulent convection.

We study the metallicity distribution and evolution in the Galactic disk based on the largest sample of open star clusters in the Galaxy. From the catalogue of 1879 open clusters in the range of Galactocentric distance (R_GC) from 4 to 20 kpc, we investigate the variation of metallicity in the Galactic disk as functions of R_GC, vertical distance (Z), and ages of the clusters. In the direction perpendicular to the Galactic plane, variation in metallicity is found to follow a stepped linear relation. We estimate a vertical metallicity gradient d[Fe/H]/dZ of -0.545+/-0.046 dex/kpc for |Z| < 0.487 kpc, and -0.075+/-0.093 dex/kpc for 0.487 < |Z| < 1.8 kpc. On average, metallicity variations above and below the Galactic plane are found to change at similar rates. The change in metallicity in the radial direction is also found to follow a two-function linear relation. We obtain a radial metallicity gradient d[Fe/H]/d[R_GC] of -0.070+/-0.002 dex/kpc for 4.0<R_GC<12.8 kpc, and -0.005+/-0.018 dex/kpc for 12.8< R_GC < 20.5 kpc which clearly shows a strong variation in the metallicity gradient when moving from the inner to the outer Galactic disk. Age-metallicity relation (AMR) is found to follow a steeper negative slope of -0.031+/-0.006 dex/Gyr for clusters older than 240 Myr, however, there is some hint of positive metallicity age gradient for younger clusters.

We present a novel approach to identify galaxy clusters that are undergoing a merger using a deep learning approach. This paper uses massive galaxy clusters spanning $0 \leq z \leq 2$ from \textsc{The Three Hundred} project, a suite of hydrodynamic re-simulations of 324 large galaxy clusters. Mock, idealised Compton-{\it y} and X-ray maps were constructed for the sample, capturing them out to a radius of $2R_{200}$. The idealised nature of these maps mean they do not consider observational effects such as foreground or background astrophysical objects, any spatial resolution limits or restriction on X-ray energy bands. Half of the maps belong to a merging population as defined by a mass increase $\Delta${\it M/M} $\geq$ 0.75, and the other half serve as a control, relaxed population. We employ a convolutional neural network architecture and train the model to classify clusters into one of the groups. A best-performing model was able to correctly distinguish between the two populations with a balanced accuracy (BA) and recall of 0.77, ROC-AUC of 0.85, PR-AUC of 0.55 and $F_{1}$ score of 0.53. Using a multichannel model relative to a single channel model, we obtain a 3\% improvement in BA score, and a 6\% improvement in $F_{1}$ score. We use a saliency interpretation approach to discern the regions most important to each classification decision. By analysing radially binned saliency values we find a preference to utilise regions out to larger distances for mergers with respect to non-mergers, greater than $\sim1.2 R_{200}$ and $\sim0.7 R_{200}$ for SZ and X-ray respectively.

We examine various physical processes that may explain the shallow high-mass slope of the IMF as well as the low SFR in star-forming molecular clouds (MCs) in the Central Molecular Zone (CMZ). We show that the strong tidal field and the tidal shear experienced by the CMZ have opposite effects on the collapse of density fluctuations and nearly compensate, but in any case have a negligible impact and can not explain these unusual properties. Similarly, we show that the intense magnetic field in the CMZ provides a negligible pressure support and, for the high densities at play should not modify the probability density function (PDF) of the turbulent gas flow in the clouds, thus affecting negligibly the slope of the IMF. However, we show that, in contrast to MCs in the Galactic disk, the ones in the CMZ experience only one single episode of turbulence injection at large scale, most likely due dominantly to bar gas inflow. Indeed, their rather short lifetime, due to their high mean densities, is similar to one typical turbulence crossing time. Consequently, according to the Hennebelle-Chabrier theory of star formation, within this 'single turbulence episode' scenario, the cloud experiences one single field of turbulence induced density fluctuations, leading eventually to gravitationally unstable prestellar cores. As shown in Hennebelle & Chabrier (2013}, this yields a flatter IMF than usual and leads to the correct observed slope for the CMZ star-forming clouds. Similarly, this single large scale turbulence event within the cloud lifetime yields a 5 to 6 lower SFR than under usual MW cloud conditions, again in agreement with the observed values. Therefore, we suggest that this 'single large scale turbulence injection' episode can explain both the shallow IMF high-mass slope and low SFR of clouds in the CMZ.

The planets' gravitational interaction causes rhythmic changes in Earth's orbital parameters (also called Milankovi\'c cycles), which have powerful applications in geology and astrochronology. For instance, the primary astronomical eccentricity cycle due to the secular frequency term (g2-g5) (~405 kyr in the recent past) utilized in deep-time analyses is dominated by Venus' and Jupiter's orbits, aka long eccentricity cycle. The widely accepted and long-held view is that (g2-g5) was practically stable in the past and may hence be used as a "metronome" to reconstruct accurate ages and chronologies. However, using state-of-the-art integrations of the solar system, we show here that (g2-g5) can become unstable over long time scales, without major changes in, or destabilization of, planetary orbits. The (g2-g5) disruption is due to the secular resonance $\sigma_{12}$ = (g1 - g2) + (s1 - s2), a major contributor to solar system chaos. We demonstrate that entering/exiting the $\sigma_{12}$ resonance is a common phenomenon on long time scales, occurring in ~40% of our solutions. During $\sigma_{12}$-resonance episodes, (g2-g5) is very weak or absent and Earth's orbital eccentricity and climate-forcing spectrum are unrecognizable compared to the recent past. Our results have fundamental implications for geology and astrochronology, as well as climate forcing because the paradigm that the longest Milankovi\'c cycle dominates Earth's astronomical forcing, is stable, and has a period of ~405 kyr requires revision.

Utilizing NICER observations, we present an analysis of the soft X-ray re-brightening event of GRS 1915+105 observed in 2021. During this event, we observed the emergence of a stable, long-lasting low-frequency quasi-periodic oscillation (LFQPO) with frequencies ranging from 0.17 to 0.21 Hz. Through a careful spectral analysis, we demonstrate that a low-temperature Compton-thick gas model well characterizes the emitted radiation. By examining the spectrum and identifying numerous absorption lines, we discerned a transition in the wind properties. This transition was marked by a shift from a state characterized by low speed, high column density, and high ionization degree to one featuring still low speed but low column density and ionization degree. Intriguingly, the presence or absence of the QPO signal is perfectly correlated with these distinct wind characteristics. The low-speed wind observed could be indicative of a 'failed wind', while the observed shift implies a transition from a magnetically to a thermally driven wind. Notably, this QPO signal exclusively manifested itself during the magnetically driven phase, suggesting the possibility of a novel perturbation associated with magnetic effects.

Merging of stellar-mass binary black holes (BBH) could take place within the accretion disk of active galactic nuclei (AGN). The resulting BH remnant is likely to accrete the disk gas at a super-Eddington rate, launching a fast, quasi-spherical outflow (wind). Particles will be accelerated by shocks driven by the wind, subsequently interacting with the shocked disk gas or radiation field through hadronic processes and resulting in the production of high-energy neutrinos and potential electromagnetic (EM) emissions. This study delves into the intricate evolution of the shock driven by the remnant BH wind within AGN disks. Subsequently, we calculated the production of neutrinos and the expected detection numbers for a single event, along with their contributions to the overall diffuse neutrino background. Our analysis, considering various scenarios, reveals considerable neutrino production and possible detection by IceCube for nearby events. The contribution of the remnant BH winds on the diffuse neutrino background is minor due to the low event rate density, but it can be improved to some extent for some optimistic parameters. We also propose that there could be two neutrino/EM bursts, one originating from the premerger BBH wind and the other from the remnant BH wind, with the latter typically having a time gap to the GW event of around tens of days. When combined with the anticipated gravitational waves (GW) emitted during the BBH merger, such a system emerges as a promising candidate for joint observations involving neutrinos, GWs, and EM signals.

Expanding X-ray halo or rings appear when short pulses of X-ray radiation from a background source are scattered by clouds of dust in the Milky Way. We study the X-ray rings of the brightest gamma-ray burst (GRB) 221009A, detected by the {\it Swift} X-Ray Telescope. The rings center on the GRB position and their angular radii increase with time. We identify five major expanding rings, and our modeling of their expansion history suggests that they are scattered off, respectively, from five dusty clouds at distances of 0.4-13 kpc from the observer. Given an assumed prompt X-ray fluence of this GRB, the fluxes of those rings suggest that these clouds have dust grain column densities of $10^{7\sim8}~\mathrm{cm^{-2}}$. More interestingly, our time-dependent spectral analysis of these rings show that they all experience spectral softening, i.e., getting softer as they expand, with spectral indices ranging from 2.2 to 5, consistent with what the dust scattering model predicts.

We analyze measurements of the thermal Sunyaev-Zeldovich (tSZ) effect arising in the circumgalactic medium (CGM) of $L^*$ galaxies, reported by Bregman et al. 2022 and Das et al. 2023. In our analysis we use the Faerman et al. 2017 and Faerman et al. 2020 CGM models, a new power-law model (PLM), and the TNG100 simulation. For a given $M_{\rm vir}$, our PLM has four parameters; the fraction, $f_{\rm hCGM}$, of the halo baryon mass in hot CGM gas, the ratio, $\phi_T$, of the actual gas temperature at the virial radius to the virial temperature, and the power-law indicies, $a_{P,{\rm th}}$ and $a_n$ for the thermal electron pressure and the hydrogen nucleon density. The B+22 Compton-$y$ profile implies steep electron pressure slopes ($a_{P,{\rm th}}\simeq 2$). For isothermal conditions the temperature is at least $1.1\times 10^6$ K, with a hot CGM gas mass of up to $3.5\times 10^{11}$ M$_\odot$ for a virial mass of $2.75\times 10^{12}$ M$_\odot$. However, if isothermal the gas must be expanding out of the halos. An isentropic equation of state is favored for which hydrostatic equilibrium is possible. The B+22 and D+23 results are consistent with each other and with recent (0.5-2 keV) CGM X-ray observations by Zhang et al. 2024 of Milky Way mass systems. For $M_{\rm vir}\simeq 3\times 10^{12}$ M$_\odot$, the scaled Compton pressure integrals, $E(z)^{-2/3}Y_{500}/M_{\rm vir,12}^{5/3}$, lie in the narrow range, $2.5\times 10^{-4}$ to $5.0\times 10^{-4}$ kpc$^2$, for all three sets of observations. TNG100 underpredicts the tSZ parameters by factors $\sim 0.5$ dex for the $L^*$ galaxies, suggesting that the feedback strengths and CGM gas losses are overestimated in the simulated halos at these mass scales.

The globular cluster ultraluminous X-ray source, RZ2109, is a complex and unique system which has been detected at X-ray, ultra-violet, and optical wavelengths. Based on almost 20 years of Chandra and XMM-Newton observations, the X-ray luminosity exhibits order-of-magnitude variability, with the peak flux lasting on the order of a few hours. We perform robust time series analysis on the archival X-ray observations and find that this variability is periodic on a timescale of 1.3 $\pm 0.04$ days. The source also demonstrates broad [OIII] 5007 Angstrom emission, which has been observed since 2004, suggesting a white dwarf donor and therefore an ultra-compact X-ray binary. We present new spectra from 2020 and 2022, marking eighteen years of observed [OIII] emission from this source. Meanwhile, we find that the globular cluster counterpart is unusually bright in the NUV/UVW2 band. Finally, we discuss RZ2109 in the context of the eccentric Kozai Lidov mechanism and show that the observed 1.3 day periodicity can be used to place constraints on the tertiary configuration, ranging from 20 minutes (for a 0.1 ${\rm M}_\odot$ companion) to approximately 95 minutes (for a 1 ${\rm M}_\odot$ companion), if the eccentric Kozai Lidov mechanism is at the origin of the periodic variability.

The detection of blue-shifted absorption lines likely associated with ionized Iron K-shell transitions in the X-ray spectra of many Active Galactic Nuclei (AGN) suggests the presence of a highly ionized gas outflowing with mildly relativistic velocities (0.03c-0.6c), named Ultra-Fast Outflow (UFO). Within the SUBWAYS project we characterized these winds starting from a sample of 22 radio-quiet quasars at 0.1 < z < 0.4, and compared the results with similar studies in the literature on samples of 42 local radio-quiet Seyfert galaxies and 14 high redshift radio-quiet quasars. The scope of our work is a statistical study of UFO parameters and incidence, considering key physical properties of the sources, e.g. supermassive black hole (SMBH) mass, bolometric luminosity, accretion rates and Spectral Energy Distribution, with the aim of gaining new insights into the UFO launching mechanisms. We find indications that highly luminous AGN with steeper X-ray/UV ratio, are more likely to host UFO. The presence of UFO is not significantly related to any other AGN property in our sample. These findings suggest that the UFO phenomenon may be transient. Focusing on AGN with UFO, other important results are: (1) faster UFO have larger ionization parameters and column densities; (2) X-ray radiation plays a more crucial role in driving highly ionized winds compared to UV; (3) the correlation between outflow velocity and luminosity is significantly flatter than what expected for radiatively driven winds; (4) more massive BH experience higher wind mass-losses, suppressing accretion of matter onto the BH; (5) the UFO launching radius is positively correlated with the Eddington ratio. Furthermore, our analysis suggest the involvement of multiple launching mechanisms, including radiation pressure and magneto-hydrodynamic processes, rather than pointing to a single, universally applicable mechanism.

Redbacks are millisecond pulsar binaries with low mass, irradiated companions. These systems have a rich phenomenology that can be used to probe binary evolution models, pulsar wind physics, and the neutron star mass distribution. A number of high-confidence redback candidates have been identified through searches for variable optical and X-ray sources within the localisation regions of unidentified but pulsar-like Fermi-LAT gamma-ray sources. However, these candidates remain unconfirmed until pulsations are detected. As part of the TRAPUM project, we searched for radio pulsations from six of these redback candidates with MeerKAT. We discovered three new radio millisecond pulsars, PSRs J0838$-$2527, J0955$-$3947 and J2333$-$5526, confirming their redback nature. PSR J0838$-$2827 remained undetected for two years after our discovery despite repeated observations, likely due to evaporated material absorbing the radio emission for long periods of time. While, to our knowledge, this system has not undergone a transition to an accreting state, the disappearance, likely caused by extreme eclipses, illustrates the transient nature of spider pulsars and the heavy selection bias in uncovering their radio population. Radio timing enabled the detection of gamma-ray pulsations from all three pulsars, from which we obtained 15-year timing solutions. All of these sources exhibit complex orbital period variations consistent with gravitational quadrupole moment variations in the companion stars. These timing solutions also constrain the binary mass ratios, allowing us to narrow down the pulsar masses. We find that PSR J2333$-$5526 may have a neutron star mass in excess of 2 M$_{\odot}$.

In this letter we present the first systematic spectroscopic measurements of the Near-InfraRed (NIR) hydrogen recombination lines $Pa_{\alpha}$ ($\lambda = 1.875 \mu m$) and $Br_{\beta}$ ($\lambda = 2.626 \mu m$), produced by pre-main sequence (PMS) stars. These stars, consisting of T Tauri and Herbig AeBe stars are located in the massive Galactic star-formation region NGC 3603. The measurements were obtained from JWST NIRSpec, using multi-object spectroscopy (MOS) mode. Utilising existing empirical relations between $L_{acc}$ and $L_{Br_{\gamma}}$, we have used our new measurements to formulate, for the first time, an empirical relationship between the accretion luminosity $L_{acc}$ of the stars, and the line luminosities $L_{line}$ of both $Pa_{\alpha}$ and $Br_{\beta}$. These relationships are: $\log_{10}( \frac{L_{acc}}{L_{\odot}}) = 1.42 (\pm 0.18) \times \log_{10}( \frac{L_{Pa_{\alpha}}}{L_{\odot}}) + 3.33 (\pm 0.42)$ and $\log_{10}( \frac{L_{acc}}{L_{\odot}}) = 1.47 (\pm 0.18) \times \log_{10}( \frac{L_{Br_{\beta}}}{L_{\odot}}) + 4.60 (\pm 0.57)$. These new relationships are key for roughly estimating the accretion rates for large samples of PMS stars with JWST.

We introduce cosmocnc, a Python package for computing the number count likelihood of galaxy cluster catalogues in a fast, flexible and accurate way. cosmocnc offers three types of likelihoods: an unbinned, a binned, and an extreme value likelihood. It also supports the addition of stacked cluster data, which is modelled consistently with the cluster catalogue. The unbinned likelihood, which is the main focus of the code, can take an arbitrary number of mass observables as input and deal with several complexities in the data, such as variations in the properties of the cluster observable across the survey footprint, the possibility of different clusters having measurements for different combinations of mass observables, redshift measurement uncertainties, and the presence on unconfirmed detections in the catalogue. If there are more than one mass observables, the unbinned likelihood is computed with the backward convolutional approach, a novel approach that is first implemented in cosmocnc. After developing the likelihood formalism and describing its implementation, we validate the code with synthetic Simons-Observatory-like catalogues, finding excellent agreement between their properties and cosmocnc's predictions and obtaining constraints on cosmological and scaling relation parameters featuring negligible biases. cosmocnc is publicly available at github.com/inigozubeldia/cosmocnc.

To investigate the effects of stellar feedback on the gravitational state of giant molecular clouds (GMCs), we study $^{12}$CO and $^{13}$CO ALMA maps of nine GMCs distributed throughout the Large Magellanic Cloud (LMC), the nearest star-forming galaxy to our own. We perform noise and resolution matching on the sample, working at a common resolution of 3.5 arcseconds (0.85 pc at the LMC distance of 50 kpc), and use the \textit{SCIMES} clustering algorithm to identify discrete substructure, or "clumps." We supplement these data with three tracers of recent star formation: $8\mu$m surface brightness, continuum-subtracted H$\alpha$ flux, and interstellar radiation field energy density inferred from dust emission. The $^{12}$CO clumps identified cover a range of 3.6 dex in luminosity-based mass and 2.4 dex in average $8\mu$m surface brightness, representative of the wide range of conditions of the interstellar medium in the LMC. Our observations suggest evidence for increased turbulence in these clouds. While the turbulent linewidths are correlated with clump surface density, in agreement with previous observations, we find even better correlation with the three star formation activity tracers considered, suggesting stellar energy injection plays a significant role in the dynamical state of the clumps. The excess linewidths we measure do not appear to result from opacity broadening. $^{12}$CO clumps are found to be typically less gravitationally bound than $^{13}$CO clumps, with some evidence of the kinetic-to-gravitational potential energy ratio increasing with star-formation tracers. Further multi-line analysis may better constrain the assumptions made in these calculations.

We report here results from pulse arrival time delay analysis of the eclipsing high mass X-ray binary pulsar LMC X-4 using observations made with the Rossi X-ray Timing Explorer, XMM-Newton, NuSTAR and AstroSat. Combining the orbital parameters determined from these observations with the historical measurements dating back to 1998, we have extended the $T_{\pi/2}$ epoch history of LMC X-4 by about 4600 binary orbits spanning about 18 years. We also report mid-eclipse time measurements ($T_{ecl}$) using data obtained from wide-field X-ray monitors of MAXI-GSC and Swift-BAT. Combining the new $T_{\pi/2}$ and $T_{ecl}$ estimates with all the previously reported values, we have significantly improved the orbital evolution measurement, which indicates that the orbital period is evolving at a time scale ($P_{\rm orb}/\dot{P}_{\rm orb}$ ) of about 0.8 Myr. For the first time in an accreting X-ray pulsar system, we confirm the existence of a second derivative of the orbital period, having an evolution time scale ($\dot{P}_{orb}/\ddot{P}_{orb}$) of about 55 yr. Detection of a second derivative of the orbital period in LMC X-4 makes its orbital evolution timescale more uncertain, which may also be true for other HMXBs. Independent solutions for the orbital evolution measurement using the mid-eclipse data and the pulse timing data are consistent with each other, and help us put an upper limit of 0.009 on the eccentricity of the binary system.

The internal structure and abundance of dark matter halos and subhalos are powerful probes of the nature of dark matter. In order to compare observations with dark matter models, accurate theoretical predictions of these quantities are needed. We present a fast and accurate method to describe the tidal evolution of subhalos within their parent halo, based on a semi-analytic approach. We first consider idealized N-body simulations of subhalos within their host halo, using a generalized mass density profile that describes their properties in a variety of dark matter models at infall, including popular warm, cold, and self-interacting ones. Using these simulations we construct tidal "tracks" for the evolution of subhalos based on their conditions at infall. Second, we use the results of these simulations to build semi-analytic models (SAMs) for tidal effects, including stripping and heating and implement them within the code GALACTICUS. Our SAMs can accurately predict the tidal evolution of both cored and cuspy subhalos, including the bound mass and density profiles, providing a powerful and efficient tool for studying the post-infall properties of subhalos in different dark matter models.

The Interstellar Boundary Explorer (IBEX) satellite collects data on energetic neutral atoms (ENAs) that provide insight into the heliosphere, the region surrounding our solar system and separating it from interstellar space. IBEX collects information on these particles and on extraneous ``background'' particles. While IBEX records how and when the different particles are observed, it does not distinguish between heliospheric ENA particles and incidental background particles. To address this issue, all IBEX data has historically been manually labeled as ``good'' ENA data, or ``bad'' background data. This manual culling process is incredibly time-intensive and contingent on subjective, manually-induced decision thresholds. In this paper, we develop a three-stage automated culling process, called LOTUS, that uses random forests to expedite and standardize the labelling process. In Stage 1, LOTUS uses random forests to obtain probabilities of observing true ENA particles on a per-observation basis. In Stage 2, LOTUS aggregates these probabilities to obtain predictions within small windows of time. In Stage 3, LOTUS refines these predictions. We compare the labels generated by LOTUS to those manually generated by the subject matter expert. We use various metrics to demonstrate that LOTUS is a useful automated process for supplementing and standardizing the manual culling process.

As the gravitational wave detector network is upgraded and the sensitivity of the detectors improves, novel scientific avenues open for exploration. For example, tests of general relativity will become more accurate as smaller deviations can be probed. Additionally, the detection of lensed gravitational waves becomes more likely. However, these new avenues could also interact with each other, and a gravitational wave event presenting deviations from general relativity could be mistaken for a lensed one. Here, we explore how phenomenological deviations from general relativity or binaries of exotic compact objects could impact those lensing searches focusing on a single event. We consider strong lensing, millilensing, and microlensing and find that certain phenomenological deviations from general relativity may be mistaken for all of these types of lensing. Therefore, our study shows that future candidate lensing events would need to be carefully examined to avoid a false claim of lensing where instead a deviation from general relativity has been seen.

The study of the Two-Body and Circular Restricted Three-Body Problems in the field of aerospace engineering and sciences is deeply important because they help describe the motion of both celestial and artificial satellites. With the growing demand for satellites and satellite formation flying, fast and efficient control of these systems is becoming ever more important. Global linearization of these systems allows engineers to employ methods of control in order to achieve these desired results. We propose a data-driven framework for simultaneous system identification and global linearization of both the Two-Body Problem and Circular Restricted Three-Body Problem via deep learning-based Koopman Theory, i.e., a framework that can identify the underlying dynamics and globally linearize it into a linear time-invariant (LTI) system. The linear Koopman operator is discovered through purely data-driven training of a Deep Neural Network with a custom architecture. This paper displays the ability of the Koopman operator to generalize to various other Two-Body systems without the need for retraining. We also demonstrate the capability of the same architecture to be utilized to accurately learn a Koopman operator that approximates the Circular Restricted Three-Body Problem.

Following our previous studies on gravitational dark matter (GDM) production in the early Universe of preheating,we forecast high-frequency gravitational wave (GW) spectra as the indirect probe of such GDM. We use proper lattice simulations to handle the resonance,and to solve the GW equation of motion with resonance induced scalar field excitations as the source term.Our numerical results show that the Higgs scalar excitations in the Higgs preheating model give rise to magnitudes of GW energy density spectra of order $10^{-10}$ at frequencies $10-10^{3}$ MHz, whereas the inflaton fluctuation excitations in the inflaton self-resonant preheating model yield magnitudes of GW energy density spectrum up to $10^{-9}~(10^{-11})$ at frequencies near $30~(2)$ MHz for the index $n=4~(6)$.

We analyze the recent EMPRESS data on the $^4$He abundance in Big Bang nucleosynthesis as a function of the Higgs vacuum expectation value $v$ and compare our calculation to the recently published work of Burns et al.~[1]. The EMPRESS result for the $^4$He abundance can be explained within $2\sigma$ by $0.03 \leq \delta v / v \leq 0.07$. However, as first noted in~[1], a Higgs VEV in this range would worsen the discrepancy between theory and experiment for the deuterium abundance significantly.

It is known that the four-dimensional effective field theory arising from heterotic string theory is general relativity with both a Chern-Simons and Gauss-Bonnet term. We study the propagation of gravitational waves in this combination of Chern-Simons and Gauss-Bonnet gravity, both of which have an associated scalar field, the axion and the dilaton respectively, that are kinetically coupled. We review how the combination of dynamical Chern-Simons and Gauss-Bonnet gravities can arise from string theory as corrections to general relativity and show how the gravitational wave waveform is modified in such a theory. We compare our results to a novel framework recently introduced for parametrizing the parity-violating sector (Chern-Simons), and use that to guide our construction of a similar parametrization for the parity-conserving (Gauss-Bonnet) sector. In general, we find that the contributions from the parity-violating and parity-conserving sectors are similar. Moreover, the kinetic coupling between the axion and dilaton introduces an extra contribution to the parity-violating sector of the gravitational waves. Using our parametrization, we are able to comment on initial constraints for the theory parameters, including the time variations of the axion and dilaton.

We consider F-term hybrid inflation (FHI) and SUSY breaking in the context of a B-L extension of MSSM which largely respects a global U(1) R symmetry. The hidden sector Kaehler manifold enjoys an enhanced SU(1,1)/U(1) symmetry with the scalar curvature determined by the achievement of a SUSY-breaking de Sitter vacuum without ugly tuning. FHI turns out to be consistent with data, provided that the magnitude of the emergent soft tadpole term is confined in the range (1.2-100) TeV, and it is accompanied with the production of B-L cosmic strings. If these are metastable they interpret the present observations from PTA experiments on the stochastic background of gravitational waves with dimensionless tension Gmu~(1-9.2)x10^-8. The mu parameter of MSSM arises by appropriately adapting the Giudice-Masiero mechanism and facilitates the out-of-equilibrium decay of the R saxion at a reheat temperature lower than about 71 GeV. Due to the prolonged matter dominated era the gravitational wave signal is suppressed at high frequencies. The SUSY mass scale turns out to lie in the PeV region.

The recent advancements in black hole imaging have opened a new era of probing horizon-scale physics with electromagnetic radiation. However, a feature of the observed images, a bright ring encircling a relatively dark region, has not sufficiently proved the existence of event horizons. It thus requires extreme care when studying the possibility of using such image features to examine quantum effects that may change the classical picture of black holes slightly or drastically. In this work, we investigate the image of a horizonless compact object, whose interior metric satisfies the 4D semi-classical Einstein equation non-perturbatively for the Planck constant, and whose entropy agrees with the Bekenstein-Hawking formula. Although the absence of an event horizon allows light rays to pass through the dense interior, the extremely strong redshift significantly darkens the image, making it almost identical to the classical black-hole image. In particular, if there is light emission a bit inside the surface of the object, the intensity around the inner shadow is slightly enhanced, which could be a future observable prediction to characterize the object. We also find through a phenomenological parameter that the image is further darkened due to interactions inside. Thus, the image is consistent with current observations, and the object could be a candidate for black holes in quantum theory.

We explore the anomaly detection framework based on Normalizing Flows (NF) models introduced in \cite{PhysRevC.106.065802} to detect the presence of a large (destabilising) dense matter phase transition in neutron star (NS) observations of masses and radii, and relate the feasibility of detection with parameters of the underlying mass-radius sequence, which is a functional of the dense matter equation of state. Once trained on simulated data featuring continuous $M(R)$ solutions (i.e., no phase transitions), NF is used to determine the likelihood of a first-order phase transition in a given set of $M(R)$ observations featuring a discontinuity, i.e., perform the anomaly detection. Different mock test sets, featuring two branch solutions in the $M(R)$ diagram, were parameterized by the NS mass at which the phase transition occurs, $M_c$, and the radius difference between the heaviest hadronic star and lightest hybrid star, $\Delta R$. We analyze the impact of these parameters on the NF performance in detecting the presence of a first-order phase transition. Among the results, we report that given a set of 15 stars with radius uncertainty of $0.2$ km, a detection of a two-branch solution is possible with 95\% accuracy if $\Delta R > 0.4$ km.

Matter in compact stars is dense enough that transient events within the star could have sufficiently high energies to produce detectable gravitational waves (GWs). These GWs could be used to constrain the equation of state (EoS) for matter in the star and could reveal that there is more than one type of EoS at play in the population, implying that multiple types of compact stars exist. One of these types could be quark stars, composed almost entirely of stable quark matter, and observing GWs is a way to test for the strange matter EoS. Here we explore the possibility that, if fundamental (f-) mode oscillations in pulsars are induced by a pulsar glitch, then these oscillations might produce detectable GWs. We use the existing population of pulsars and their glitches, as well as a much larger synthesized population, along with 15 EoSs (8 for neutron stars and 7 for quark stars) to generate frequencies, damping times, and GW strengths for each. We find that of the EoSs examined, all quark star EoSs produce narrower distributions of f-mode frequency than neutron star EoSs. This result, along with other elements of the data, could be used to differentiate between GWs (or other signals from f-modes) originating from neutron stars and quark stars and thus could confirm the existence of quark stars. We also find that GW astronomy is a potentially viable method for detecting a larger population of pulsars which are not observable electromagnetically and that future GW observatories have the possibility to greatly expand this capability.