Papers for Thursday, Sep 08 2022

We examine further the ability of the New Early Dark Energy model (NEDE) to resolve the current tension between the Cosmic Microwave Background (CMB) and local measurements of $H_0$ and the consequences for inflation. We perform new Bayesian analyses, including the current datasets from the ground-based CMB telescopes Atacama Cosmology Telescope (ACT), the South Pole Telescope (SPT), and the BICEP/Keck telescopes, employing an updated likelihood for the local measurements coming from the S$H_0$ES collaboration. Using the S$H_0$ES prior on $H_0$, the combined analysis with Baryonic Acoustic Oscillations (BAO), Pantheon, Planck and ACT improves the best-fit by $\Delta\chi^2 = -15.9$ with respect to $\Lambda$CDM, favors a non-zero fractional contribution of NEDE, $f_{\rm NEDE} > 0$, by $4.8\sigma$, and gives a best-fit value for the Hubble constant of $H_0 = 72.09$ km/s/Mpc (mean $71.48_{-0.81}^{+0.79}$ with $68\%$ C.L.). A similar analysis using SPT instead of ACT yields consistent results with a $\Delta \chi^2 = - 23.1$ over $\Lambda$CDM, a preference for non-zero $f_{\rm NEDE}$ of $4.7\sigma$ and a best-fit value of $H_0=71.77$ km/s/Mpc (mean $71.43_{-0.84}^{+0.84}$ with $68\%$ C.L.). We also provide the constraints on the inflation parameters $r$ and $n_s$ coming from NEDE, including the BICEP/Keck 2018 data, and show that the allowed upper value on the tensor-scalar ratio is consistent with the $\Lambda$CDM bound, but, as also originally found, with a more blue scalar spectrum implying that the simplest curvaton model is now favored over the Starobinsky inflation model.

A strong nuclear kilomaser, W1, has been found in the nearby galaxy NGC 253, associated with a forming super star cluster. Kilomasers could arise from the accretion disc around supermassive stars (>10^3 Msun), hypothetical objects that have been proposed as polluters responsible for the chemical peculiarities in globular clusters. The supermassive stars would form via runaway collisions, simultaneously with the cluster. Their discs are perturbed by stellar flybys, inspiralling and colliding stars. This raises the question if an accretion disc would at all be able to survive in such a dynamic environment and mase water lines. We investigated what the predicted maser spectrum of such a disc would look like using 2D hydrodynamic simulations and compared this to the W1 kilomaser. We derived model maser spectra from the simulations by using a general maser model for appropriate disc temperatures. All our model discs survived. The model maser spectra for the most destructive case for the simulations of M = 1000 Msun are a reasonable match with the W1 kilomaser spectrum in terms of scaling, flux values and some of the signal trends. Details in the spectrum suggest that a star of a few 1000 Msun might fit even better, with 10,000 Msun clearly giving too large velocities. Our investigations thus support the hypothesis that kilomasers could pinpoint supermassive stars.

We examine the spatial distribution and orbital pole correlations of satellites in a suite of zoom-in high-resolution dissipationless simulations of Miky Way (MW) sized haloes. We use the measured distribution to estimate the incidence of satellite configurations as flattened and as correlated in their orbital pole distribution as satellite system of the Milky Way. We confirm that this incidence is sensitive to the radial distribution of subhaloes and thereby to the processes that affect it, such as artificial disruption due to numerical effects and disruption due to the central disk. Controlling for the resolution effects and bracketing the effects of the disk, we find that the MW satellite system is somewhat unusual (at the $\approx 2-3\sigma$ level) but is statistically consistent with the $\Lambda$CDM model, in general agreement with results and conclusions of other recent studies.

Different physical processes in galaxy evolution, such as galaxy mergers that lead to coalescence of dual Active Galactic Nuclei (AGN) and outflows emanating from the narrow line region, can leave their imprint on the optical spectra of AGN in the form of double-peaked narrow emission lines. To investigate the neutral gas in the centres of such AGN, we have conducted a pilot survey of HI 21-cm absorption, using the upgraded Giant Metrewave Radio Telescope (uGMRT), in radio-loud AGN whose optical spectra show double-peaked [OIII] emission lines at z<=0.4 (median z~0.14). Among the eight sources for which we could obtain clean spectra, we detect HI 21-cm absorption in three sources (detection rate of 38(+36/-20)%) and find tentative indication of absorption in two other sources. The detection rate of HI 21-cm absorption is tentatively higher for the systems that show signatures of interaction or tidal disturbance (>~50%) in the ground-based optical images than that for the systems that appear single and undisturbed (~25%). This is consistent with the high incidence of HI 21-cm absorption observed in z<=0.2 galaxy mergers. Higher spatial resolution spectroscopy is required to confirm the origin of the HI absorbing gas, i.e. either gas infalling onto the radio-loud AGN, outflowing gas ejected by the AGN, or gas in rotation on the galactic-scale or circumnuclear discs.

We present a detailed strong lensing analysis of the massive and distant ($z=0.870$) galaxy cluster ACT-CL J0102$-$4915 (ACT0102, also known as El Gordo), taking advantage of new spectroscopic data from the Multi Unit Spectroscopic Explorer (MUSE) on the Very Large Telescope, and archival imaging from the Hubble Space Telescope. Thanks to the MUSE data, we measure secure redshifts for 374 single objects, including 23 multiply lensed galaxies, and 167 cluster members of ACT0102. The observed positions of 56 multiple images, along with their new spectroscopic redshift measurements, are used as constraints for our strong lensing model. Remarkably, some multiple images are detected out to a large projected distance of $\approx 1$ Mpc from the brightest cluster galaxy, allowing us to estimate a projected total mass value of $1.84_{-0.04}^{+0.03} \times 10^{15}\, \rm M_{\odot}$ within that radius. We find that we need two extended cluster mass components, the mass contributions from the cluster members and the additional lensing effect of a foreground ($z=0.633$) group of galaxies, to predict the positions of all multiple images with a root mean square offset of $0.75"$. The main cluster-scale mass component is centered very closely to the brightest cluster galaxy and the other extended mass component is located in the north-west region of the cluster. These two mass components have very similar values of mass projected within 300 kpc from their centers, namely $2.29_{-0.10}^{+0.09}\times10^{14}\,\rm M_{\odot}$ and $2.10_{-0.09}^{+0.08}\times10^{14}\,\rm M_{\odot}$, in agreement with the major merging scenario of ACT0102. We make publicly available the lens model, including the magnification maps and posterior distributions of the model parameter values, as well as the full spectroscopic catalogue containing all redshift measurements obtained with MUSE.

Massive disk galaxies like our Milky Way should host an ancient, metal-poor, and centrally concentrated stellar population. This population reflects the star formation and enrichment in the few most massive progenitor components that coalesced at high redshift to form the proto-Galaxy. While metal-poor stars are known to reside in the inner few kiloparsecs of our Galaxy, current data do not yet provide a comprehensive picture of such a metal-poor "heart" of the Milky Way. We use information from Gaia DR3, especially the XP spectra, to construct a sample of 2 million bright (BP $<15.5$ mag) giant stars within $30^\circ$ of the Galactic Center with robust [M/H] estimates, $\delta$ [M/H] $\lesssim 0.1$. For most sample members we can calculate orbits based on Gaia RVS velocities and astrometry. This sample reveals an extensive, ancient, and metal-poor population that includes $\sim 18,000$ stars with $-2.7<$ [M/H] $<-1.5$, representing a stellar mass of $\gtrsim 5\times 10^7$ M$_\odot$. The spatial distribution of these [M/H] $<-1.5$ stars has a Gaussian extent of only $\sigma_{\mathrm{R_{GC}}} \sim 2.7$ kpc around the Galactic center, with most of these orbits being confined to the inner Galaxy. At high orbital eccentricities, there is clear evidence for accreted halo stars in their pericentral orbit phase. Stars with [M/H] $< -2$ show no net rotation, whereas those with [M/H] $\sim -1$ are rotation dominated. Most of the tightly bound stars show $[\alpha/\text{Fe}]$-enhancement and [Al/Fe]-[Mn/Fe] abundance patterns expected for an origin in the more massive portions of the proto-Galaxy. These central, metal-poor stars most likely predate the oldest part of the disk ($\tau_{\text{age}}\approx 12.5$ Gyrs), which implies that they formed at $z\gtrsim 5$, forging the proto-Milky Way.

We build a multi-output generative model for quasar spectra and the properties of their black hole engines, based on a Gaussian process latent-variable model. This model treats every quasar as a vector of latent properties such that the spectrum and all physical properties of the quasar are associated with non-linear functions of those latent parameters; the Gaussian process kernel functions define priors on the function space. Our generative model is trained with a justifiable likelihood function that allows us to treat heteroscedastic noise and missing data correctly, which is crucial for all astrophysical applications. It can predict simultaneously unobserved spectral regions, as well as the physical properties of quasars in held-out test data. We apply the model to rest-frame ultraviolet and optical quasar spectra for which precise black hole masses (based on reverberation mapping measurements) are available. Unlike reverberation-mapping studies, which require multi-epoch data, our model predicts black hole masses from single-epoch spectra, even with limited spectral coverage. We demonstrate the capabilities of the model by predicting black hole masses and unobserved spectral regions. We find that we predict black hole masses at close to the best possible accuracy.

Baryonic feedback is expected to play a key role in regulating the star formation of low-mass galaxies by producing galaxy-scale winds associated with mass-loading factors $\beta\!\sim\!1\!-\!50$. We have tested this prediction using a sample of 19 nearby systems with stellar masses $10^7\!<\!M_\star/{\rm M}_{\odot}\!<\!10^{10}$, mostly lying above the main sequence of star-forming galaxies. We used MUSE@VLT optical integral field spectroscopy to study the warm ionised gas kinematics of these galaxies via a detailed modelling of their H$\alpha$ emission line. The ionised gas is characterised by irregular velocity fields, indicating the presence of non-circular motions of a few tens of km/s within galaxy discs, but with intrinsic velocity dispersion of $40$-$60$ km/s that are only marginally larger than those measured in main-sequence galaxies. Galactic winds, defined as gas at velocities larger than the galaxy escape speed, encompass only a few percent of the observed fluxes. Mass outflow rates and loading factors are strongly dependent on $M_\star$, star formation rate (SFR), SFR surface density and specific SFR. For $M_\star$ of $10^8$ M$_\odot$ we find $\beta\simeq0.02$, which is more than two orders of magnitude smaller than the values predicted by theoretical models of galaxy evolution. In our galaxy sample, baryonic feedback stimulates a gentle gas cycle rather than causing a large-scale blow out.

We investigate the connection between supermassive black holes (SMBHs) and their host dark matter halos in the local universe using the clustering statistics and luminosity function of AGN from the Swift/BAT AGN Spectroscopic survey (BASS DR2). By forward-modeling AGN activity into snapshot halo catalogs from N-body simulations, we test a scenario in which SMBH mass correlates with dark matter (sub)halo mass for fixed stellar mass. We compare this to a model absent of this correlation, where stellar mass alone determines the SMBH mass. We find that while both simple models are able to largely reproduce the abundance and overall clustering of AGN, the model in which black hole mass is tightly correlated with halo mass is preferred by the data by $1.8\sigma$. When including an independent measurement on the black hole mass-halo mass correlation, this model is preferred by $4.6\sigma$. We show that the clustering trends with black hole mass can further break the degeneracies between the two scenarios, and that our preferred model reproduces the measured clustering differences on 1-halo scales between large and small black hole masses. These results indicate that the halo binding energy is fundamentally connected to the growth of supermassive black holes.

Mapping the Galactic spiral structure is a difficult task since the Sun is located in the Galactic plane and because of dust extinction. For these reasons, molecular masers in radio wavelengths have been used with great success to trace the Milky Way spiral arms. Recently, Gaia parallaxes have helped in investigating the spiral structure in the Solar extended neighborhood. In this paper, we propose to determine the location of the spiral arms using Cepheids since they are bright, young supergiants with accurate distances (they are the first ladder of the extragalactic distance scale). They can be observed at very large distances; therefore, we need to take the Galactic warp into account. Thanks to updated mid-infrared photometry and to the most complete catalog of Galactic Cepheids, we derived the parameters of the warp using a robust regression method. Using a clustering algorithm, we identified groups of Cepheids after having corrected their Galactocentric distances from the (small) effects of the warp. We derived new parameters for the Galactic warp, and we show that the warp cannot be responsible for the increased dispersion of abundance gradients in the outer disk reported in previous studies. We show that Cepheids can be used to trace spiral arms, even at large distances from the Sun. The groups we identify are consistent with previous studies explicitly deriving the position of spiral arms using young tracers (masers, OB(A) stars) or mapping overdensities of upper main-sequence stars in the Solar neighborhood thanks to Gaia data.

If dark matter resides in a hidden sector minimally coupled to the Standard Model, another particle within the hidden sector might dominate the energy density of the early universe temporarily, causing an early matter-dominated era (EMDE). During an EMDE, matter perturbations grow more rapidly than they would in a period of radiation domination, which leads to the formation of microhalos much earlier than they would form in standard cosmological scenarios. These microhalos boost the dark matter annihilation signal, but this boost is highly sensitive to the small-scale cut-off in the matter power spectrum. If the dark matter is sufficiently cold, this cut-off is set by the relativistic pressure of the particle that dominates the hidden sector. We determine the evolution of dark matter density perturbations in this scenario, obtaining the power spectrum at the end of the EMDE. We analyze the suppression of perturbations due to the relativistic pressure of the dominant hidden sector particle and express the cut-off scale and peak scale for which the matter power spectrum is maximized in terms of the properties of this particle. We also supply transfer functions to relate the matter power spectrum with a small-scale cut-off resulting from the pressure of the dominant hidden sector particle to the matter power spectrum that results from a cold hidden sector. These transfer functions facilitate the quick computation of accurate matter power spectra in EMDE scenarios with initially hot hidden sectors and allow us to identify which models significantly enhance the microhalo abundance.

Metallicity distributions (MDs) of globular clusters (GCs) provide crucial clues for the assembly and star formation history of their host galaxies. GC colors, when GCs are old, have been used as a proxy of GC metallicities. Bimodal GC color distributions (CDs) observed in most early-type galaxies have been interpreted as bimodal MDs for decades, suggesting the presence of merely two GC subpopulations within single galaxies. However, the conventional view has been challenged by a new theory that nonlinear metallicity-to-color conversion can cause bimodal CDs from unimodal MDs. The unimodal MDs seem natural given that MDs involved many thousand protogalaxies. The new theory has been tested and corroborated by various observational and theoretical studies. Here we examine the nonlinear nature of GC color-metallicity relations (CMRs) using photometric and spectroscopic GC data of NGC 5128 (Centaurus A) and NGC 4594 (Sombrero), in comparison with stellar population simulations. We find that, with a slight offset in color, the overall shapes of observed and modeled CMRs agree well for all available colors. Diverse color-depending morphologies of GC CDs of the two galaxies are well reproduced based on their observed spectroscopic MDs via our CMR models. The results corroborate the nonlinear CMR interpretation of the GC color bimodality, shedding further light on theories of galaxy formation.

Pulsar timing arrays (PTAs) are sensitive to oscillations in the gravitational potential along the line-of-sight due to ultralight particle pressure. We calculate the probing power of PTAs for ultralight bosons across all frequencies, from those larger than the inverse observation time to those smaller than the inverse distance to the pulsar. We show that since the signal amplitude grows faster compared to the degradation in PTA sensitivity at frequencies smaller than inverse observation time, the discovery potential can be extended towards lower masses by over three decades, maintaining high precision. We demonstrate that existing 15-year PTA data can robustly detect or rule out an ultralight component in the mass range $10^{-26}\!-\!10^{-23}\,{\rm eV}$ with fraction down to percent level of the total dark matter. Non-detection, together with other bounds in different mass regimes, will imply that ultralight scalar/axion can comprise at most $1-10\%$ of dark matter in the $10^{-30}\!-\!10^{-17}$ eV regime. With 30 years of observation, current PTAs can extend the reach down to $0.1 \%$, while next-generation PTAs such as SKA can probe the $0.01\%$ level. We also generalize the analysis and derive predictions for probing ultralight vector (i.e. dark photon) and spin-2 dark matter components.

The narrow-line Seyfert 1 galaxy I Zwicky 1 shows a unique and complex system of ionised gas in outflow, which consists of an ultra-fast wind and a two-component warm absorber. In the last two decades, XMM-Newton monitored the source multiple times enabling the study of the long-term variability of the various outflows. Plasma in photoionisation equilibrium with the ionising source responds and varies accordingly to any change of the ionising luminosity. However, detailed modelling of the past RGS data has shown no correlation between the plasma ionisation state and the ionising continuum, revealing a complex long-term variability of the multi-phase warm absorber. Here, we present a new observation of I Zwicky 1 by XMM-Newton taken in early 2020 characterised by a lower X-ray flux state. The soft X-ray spectrum from the RGS reveals the two components of the warm absorber with $\log \xi \sim -1.0$ and $\log \xi \sim 1.7$. Comparing our results with the previous observations, the ionisation state of the two absorbing gas components is continuously changing, following the same unpredictable behaviour. The new results strengthen the scenario in which the ionisation state of the warm absorber is driven by the density of the gas rather than the ionising luminosity. In particular, the presence of a radiation driven, inhomogeneous clumpy outflow may explain both the variability in ionisation throughout the years and the line-locked N V system observed in the UV band. Finally, the EPIC-pn spectrum reveals an ultra-fast wind with an outflow velocity of $\sim 0.26c$ and ionisation parameter of $\log \xi \sim 3.8$.

We present the method behind HubPUG, a software tool built for recovering systemic proper motions (PMs) of Hubble Space Telescope (HST) fields with two epochs of observations by utilizing stars observed by Gaia as a foreground frame of reference. HST PM experiments have typically relied on the use of distant background galaxies or quasi-stellar objects (QSOs) as stationary sources against which to measure PMs. Without consistent profiles, background galaxies are more difficult to centroid, but benefit on-aggregate from their large numbers. QSOs, though they can be fit with stellar point-spread functions, are sparse, with most fields containing none. Historically, the use of stars as references against which to measure PMs would have been difficult because they have individual PMs of their own. However, Gaia has now provided positions and PMs for over 1.4 billion stars, which are much more likely to be well-imaged in the fields around targets versus background sources, have predictable stellar profiles, and require less observing time per-image for good signal-to-noise. This technique allows us to utilize the power of Gaia to measure the PM of targets too faint for Gaia to observe itself. We have recovered PMs for the Milky Way satellites Sculptor and Draco with comparable uncertainties over HST-only and Gaia-only measurements, limited primarily by the current capabilities of the Gaia data. We also show the promise of this method for satellites of M31 with a new PM measurement for Andromeda VII.

Comparisons of atmospheric retrievals can reveal powerful insights on the strengths and limitations of our data and modeling tools. In this paper, we examine a sample of 5 similar effective temperature (Teff) or spectral type L dwarfs to compare their pressure-temperature (P-T) profiles. Additionally, we explore the impact of an object's metallicity and the observations' signal-to-noise (SNR) on the parameters we can retrieve. We present the first atmospheric retrievals: 2MASS J15261405$+$2043414, 2MASS J05395200$-$0059019, 2MASS J15394189$-$0520428, and GD 165B increasing the small but growing number of L-dwarfs retrieved. When compared to atmospheric retrievals of SDSS J141624.08+134826.7, a low-metallicity d/sdL7 primary in a wide L+T binary, we find similar Teff sources have similar P-T profiles with metallicity differences impacting the relative offset between their P-T profiles in the photosphere. We also find that for near-infrared spectra, when the SNR is $\gtrsim80$ we are in a regime where model uncertainties dominate over data measurement uncertainties. As such, SNR does not play a role in the retrieval's ability to distinguish between a cloud-free and cloudless model, but may impact the confidence of the retrieved parameters. Lastly, we also discuss how to break cloud model degeneracies and the impact of extraneous gases in a retrieval model.

The Wide-field Imager for Solar Probe (WISPR) onboard Parker Solar Probe (PSP), observing in white light, has a fixed angular field of view, extending from 13.5 degree to 108 degree from the Sun and approximately 50 degree in the transverse direction. In January 2021, on its seventh orbit, PSP crossed the heliospheric current sheet (HCS) near perihelion at a distance of 20 solar radii. At this time, WISPR observed a broad band of highly variable solar wind and multiple coronal rays. For six days around perihelion, PSP was moving with an angular velocity exceeding that of the Sun. During this period, WISPR was able to image coronal rays as PSP approached and then passed under or over them. We have developed a technique for using the multiple viewpoints of the coronal rays to determine their location (longitude and latitude) in a heliocentric coordinate system and used the technique to determine the coordinates of three coronal rays. The technique was validated by comparing the results to observations of the coronal rays from Solar and Heliophysics Observatory (SOHO) / Large Angle and Spectrometric COronagraph (LASCO)/C3 and Solar Terrestrial Relations Observatory (STEREO)-A/COR2. Comparison of the rays' locations were also made with the HCS predicted by a 3D MHD model. In the future, results from this technique can be used to validate dynamic models of the corona.

Environment is one of the primary drivers of galaxy evolution, capable of transforming galaxies from star forming to quiescent via multiple mechanisms. Despite its importance, however, we still do not have a clear view of how environmental quenching proceeds even in the most extreme environments: galaxy clusters and their progenitor proto-clusters. Recent advances in infrared capabilities have enabled transformative progress not only in the identification of these structures but in detailed analyses of quiescence, obscured star formation, and molecular gas in (proto-)cluster galaxies across cosmic time. In this review, we will discuss the current state of the literature regarding the quenching of galaxies in (proto-)clusters from the observational, infrared perspective. Our improved understanding of environmental galaxy evolution comes from unique observables across the distinct regimes of the near-, mid-, and far-infrared, crucial in the push to high redshift where galaxy growth is dominated by highly extincted, infrared-bright galaxies.

The Legacy Survey of Space and Time (LSST) at the Vera C. Rubin Observatory will discover tens of thousands of extragalactic transients each night. The high volume of alerts demands immediate classification of transient types in order to prioritize observational follow-ups before events fade away. We use host galaxy features to classify transients, thereby providing classification upon discovery. In contrast to past work that focused on distinguishing Type Ia and core-collapse supernovae (SNe) using host galaxy features that are not always accessible (e.g., morphology), we determine the relative likelihood across $12$ transient classes based on only 19 host apparent magnitudes and colors from $10$ optical and IR photometric bands. We develop both binary and multiclass classifiers, using kernel density estimation to estimate the underlying distribution of host galaxy properties for each transient class. Even in this pilot study, and ignoring relative differences in transient class frequencies, we distinguish eight transient classes at purities significantly above the 8.3% baseline (based on a classifier that assigns labels uniformly and at random): tidal disruption events ($48\%\pm27\%$, where $\pm$ indicates the 95% confidence limit), SNe Ia-91bg ($32\%\pm18\%$), SNe Ia-91T ($23\%\pm11\%$), SNe Ib ($23\%\pm13\%$), SNe II ($17\%\pm2\%$), SNe IIn ($17\%\pm6\%$), SNe II P ($16\%\pm4\%$), and SNe Ia ($10\%\pm1\%$). We demonstrate that our model is applicable to LSST and estimate that our approach may accurately classify 59% of LSST alerts expected each year for SNe Ia, Ia-91bg, II, Ibc, SLSN-I, and tidal disruption events. Our code and dataset are publically available.

Quasi-Periodic Erupters (QPEs) are a remarkable class of objects exhibiting very large amplitude quasi-periodic X-ray flares. Although numerous dynamical models have been proposed to explain them, relatively little attention has been given to using the properties of their radiation to constrain their dynamics. Here we show that the observed luminosity, spectrum, repetition period, duty cycle, and fluctuations in the latter two quantities point toward a model in which: a main sequence star on a moderately eccentric orbit around a supermassive black hole periodically transfers mass to the Roche lobe of the black hole; orbital dynamics lead to mildly-relativistic shocks near the black hole; and thermal X-rays at the observed temperature are emitted by the gas as it flows away from the shock. Strong X-ray irradiation of the star by the flare itself augments the mass transfer, creates fluctuations in flare timing, and stirs turbulence in the stellar atmosphere that amplifies magnetic field to a level at which magnetic stresses can accelerate infall of the transferred mass toward the black hole.

Supermassive stars (SMSs), with masses $>10^5$ M$_\odot$, have been proposed as the possible progenitors of the most extreme supermassive black holes observed at redshifts $z>6-7$. In this scenario ('direct collapse'), a SMS accretes at rates $>0.1$ M$_\odot$ yr$^{-1}$ until it collapses to a black hole via the general-relativistic (GR) instability. Rotation plays a crucial role in the formation of such supermassive black hole seeds. The centrifugal barrier appears as particularly strong in this extreme case of star formation. Moreover, rotation impacts sensitively the stability of SMSs against GR, as well as the subsequent collapse. In particular, it might allow for gravitational wave emission and ultra-long gamma-ray bursts at black hole formation, which represents currently the main observational signatures proposed in the literature for the existence of such objects. Here, I present the latest models of SMSs accounting for accretion and rotation, and discuss some of the open questions and future prospects in this research line.

NASA sent the DART (Double Asteroid Redirection Test) mission to impact Dimorphos, the satellite of the asteroid binary system (65803) Didymos. DART will release LICIACube prior to impact to obtain high-resolution post-impact images. The impact will produce a crater and a large amount of material ejected at high speed (several tens of m/s), producing an ejecta cone that will quickly disperse. We analyzed an additional effect: the lofting of material at low velocity due to the generation of seismic waves that propagate inside Dimorphos, producing surface shaking far from the impact point. We divide the process into different stages: from the generation of impact-induced waves, the interaction of them with surface particles, the ejection of dust particles at velocities, and the prediction of the observability of the dust coma and trail. We anticipate the following observable effects: i) generation of a dust cloud that will produce a hazy appearance of Dimorphos' surface, detectable by LICIACube; ii) brightness increase of the binary system due to enhancement of the cross section produced by the dust cloud; iii) generation of a dust trail, similar to those observed in some Active Asteroids, which can last for several weeks after impact. Numerical prediction of the detectability of these effects depends on the amount and size distribution of ejected particles, which are largely unknown. In case these effects are observable, an inversion method can be applied to compute the amount of ejected material and its velocity distribution, and discuss the relevance of the shaking process.

Neutron stars provide a unique laboratory for studying matter at extreme pressures and densities. While there is no direct way to explore their interior structure, X-rays emitted from these stars can indirectly provide clues to the equation of state (EOS) of superdense nuclear matter through the inference of the star's mass and radius. However, inference of EOS directly from a star's X-ray spectra is extremely challenging and is complicated by systematic uncertainties. The current state of the art is to use simulation-based likelihoods in a piece-wise method, which first infer the star's mass and radius to reduce the dimensionality of the problem, and from those quantities infer the EOS. We demonstrate a series of enhancements to the state of the art, in terms of realistic uncertainty quantification and improved regression of physical properties with machine learning. We also demonstrate novel inference of the EOS directly from the high-dimensional spectra of observed stars, avoiding the intermediate mass-radius step. Our network is conditions on the sources of uncertainty of each star, allowing for natural and complete propagation of uncertainties to the EOS.

(Abridged) Physical processes such as accretion shocks are thought to be common in the protostellar phase, where the envelope component is still present, and they can release molecules from the dust to the gas phase, altering the original chemical composition of the disk. Consequently, the study of accretion shocks is essential for a better understanding of the physical processes at disk scales and their chemical output. The purpose of this work is to assess the characteristics of accretion shocks traced by sulfur-related species. We present ALMA high angular resolution observations (0.1") of the Class I protostar Oph-IRS 44. The continuum emission at 0.87 mm is observed, together with sulfur-related species such as SO, SO$_{2}$, and $^{34}$SO$_{2}$. Six lines of SO$_{2}$, two lines of $^{34}$SO$_{2}$, and one line of SO are detected toward IRS 44. The emission of all the detected lines peaks at ~0.1" (~14 au) from the continuum peak and we find infalling-rotating motions inside 30 au. However, only redshifted emission is seen between 50 and 30 au. Colder and more quiescent material is seen toward an offset region located at a distance of ~400 au from the protostar, and we do not find evidence of a Keplerian profile in these data. Accretion shocks are the most plausible explanation for the high temperatures, high densities, and velocities found for the SO$_{2}$ emission. When material enters the disk--envelope system, it generates accretion shocks that increase the dust temperature and desorb SO$_{2}$ molecules from dust grains. High-energy SO$_{2}$ transitions (~200 K) seem to be the best tracers of accretion shocks that can be followed up by future higher angular resolution ALMA observations and compared to other species to assess their importance in releasing molecules from the dust to the gas phase.

Much of coronal hole (CH) research is focused upon determining the boundary and calculating the open flux as accurately as possible. However, the observed boundary itself is worthy of investigation, and holds important clues to the physics transpiring at the interface between the open and closed fields. This Letter reports a powerful new method, an application of the correlation integral which we call correlation dimension mapping (CDM), by which the irregularity of a CH boundary can be objectively quantified. This method highlights the most important spatial scales involved in boundary dynamics, and also allows for easy temporal analysis of the boundary. We apply this method to an equatorial CH bounded on two sides by helmet streamers and on the third by a small pseudostreamer, which we observed at maximum cadence for an hour on 2015 June 4. We argue that the relevant spatial scales are in the range of $\sim 5-20$ Mm, and we find that boundary complexity depends measurably upon the nature of the neighboring closed structure. The boundary along the pseudostreamer shows signs of highly-localized, intermittent complexity variability, likely associated with abrupt changes in the magnetic topology, which would be elegantly explained by interchange reconnection. By contrast, the helmet streamer boundary supports long-lived high-complexity regions. These findings support the recent predictions of interchange reconnection occurring at very small scales in the corona.

White dwarfs are a class of stars with unique physical properties. They present many challenging problems whose solution requires the application of advanced theories of dense matter, state-of-the-art experimental techniques, and extensive computing efforts. New ground- and space-based observatories will soon provide an increasingly detailed view of white dwarf stars and reveal new phenomena that will challenge our models. This review is an introduction for researchers who are not in the field of white dwarf astrophysics with the intent to entice them to contribute their expertise to advance our knowledge of these exotic stars. We discuss a wide variety of currently unsolved or partially resolved problems that are broadly related to equations of state, transport processes and opacities.

Exoplanets orbiting in the habitable zone around M-dwarf stars have been prime targets in the search for life due to the long lifetimes of the host star, the prominence of such stars in the galaxy, and the apparent excess of terrestrial planets found around M-dwarfs. However, the heightened stellar activity of M-dwarfs and the often tidally locked planets in these systems have raised questions about the habitability of these planets. In this letter we examine another significant challenge that may exist: these systems seem to lack the architecture necessary to deliver asteroids to the habitable terrestrial planets, and asteroid impacts may play a crucial role in the origin of life. The most widely accepted mechanism for producing a stable asteroid belt and the late stage delivery of asteroids after gas disk dissipation requires a giant planet exterior to the snow line radius. We show that none of the observed systems with planets in the habitable zone of their star also contain a giant planet and therefore are unlikely to have stable asteroid belts. We consider the locations of observed giant planets relative to the snow line radius as a function of stellar mass and find that there is a population of giant planets outside of the snow line radius around M-dwarfs. Therefore, asteroid belt formation around M-dwarfs is generally possible. However, we find that multi-planetary system architectures around M-dwarfs can be quite different from those around more massive stars.

We present the results of a model-based search for continuous gravitational waves from the low-mass X-ray binary Scorpius X-1 using LIGO detector data from the third observing run of Advanced LIGO, Advanced Virgo and KAGRA. This is a semicoherent search which uses details of the signal model to coherently combine data separated by less than a specified coherence time, which can be adjusted to balance sensitivity with computing cost. The search covered a range of gravitational-wave frequencies from 25Hz to 1600Hz, as well as ranges in orbital speed, frequency and phase determined from observational constraints. No significant detection candidates were found, and upper limits were set as a function of frequency. The most stringent limits, between 100Hz and 200Hz, correspond to an amplitude h0 of about 1e-25 when marginalized isotropically over the unknown inclination angle of the neutron star's rotation axis, or less than 4e-26 assuming the optimal orientation. The sensitivity of this search is now probing amplitudes predicted by models of torque balance equilibrium. For the usual conservative model assuming accretion at the surface of the neutron star, our isotropically-marginalized upper limits are close to the predicted amplitude from about 70Hz to 100Hz; the limits assuming the neutron star spin is aligned with the most likely orbital angular momentum are below the conservative torque balance predictions from 40Hz to 200Hz. Assuming a broader range of accretion models, our direct limits on gravitational-wave amplitude delve into the relevant parameter space over a wide range of frequencies, to 500Hz or more.

In the hierarchical view of star formation, giant molecular gas clouds (GMCs) undergo fragmentation to form small-scale structures made up of stars and star clusters. Here we study the connection between young star clusters and cold gas across a range of extragalactic environments by combining the high resolution (1") PHANGS-ALMA catalogue of GMCs with the star cluster catalogues from PHANGS-HST. The star clusters are spatially matched with the GMCs across a sample of 11 nearby star-forming galaxies with a range of galactic environments (centres, bars, spiral arms, etc.). We find that after 4-6 Myr the star clusters are no longer associated with any gas clouds. Additionally, we measure the autocorrelation of the star clusters and GMCs as well as their cross-correlation to quantify the fractal nature of hierarchical star formation. Young ($\leq$ 10 Myr) star clusters are more strongly autocorrelated on kpc and smaller spatial scales than the >10 Myr stellar populations, indicating that the hierarchical structure dissolves over time.

Rocky exoplanets orbiting in the habitable zone (HZ) of nearby M dwarfs provide unique opportunities for characterizing their atmospheres and searching for biosignature gases. TRAPPIST-1e, a temperate Earth-sized exoplanet in the HZ of a nearby M dwarf, is arguably the most favorable target for ground- and space-based atmospheric characterization by the extremely large telescopes (ELTs) and the James Webb Space Telescope (JWST). To inform future observations in reflected and emitted lights using these upcoming telescopes, we simulate the high-resolution reflection and emission spectra for TRAPPIST-1e for both modern and prebiotic Earth-like atmospheric compositions. To demonstrate the effects of wavelength-dependent albedo on climate and spectra, we further simulate five albedo scenarios for each atmospheric composition: cloudy modern Earth-like, cloud-free modern Earth-like, cloudy ocean planet, 100 per cent cloudy planet, and wavelength-independent albedo of 0.31. We use the recent Mega-MUSCLES spectral energy distribution (SED) of TRAPPIST-1 for our models. We show that the O$_2$ + CH$_4$ and O$_3$ + CH$_4$ biosignature pairs as well as climate indicators (CO$_2$ and H$_2$O) show features in both high-resolution reflection and emission spectra of TRAPPIST-1e that the ELTs can search for. Our high-resolution database for modern and prebiotic Earth TRAPPIST-1e models with various surface compositions and cloud distributions provides a tool for observers to train retrieval algorithms and plan observation strategies to characterize this potentially habitable world.

The outer protoplanetary disks (PPDs) can be subject to the magnetorotational instability (MRI) and the vertical shear instability (VSI). While both processes can drive turbulence in the disk, existing numerical simulations have studied them separately. In this paper, we conduct global 3D non-ideal magnetohydrodynamic (MHD) simulations for outer PPDs with ambipolar diffusion and instantaneous cooling, and hence conductive to both instabilities. Given the range of ambipolar Els\"{a}sser numbers ($Am$) explored, it is found that the VSI turbulence dominates over the MRI when ambipolar diffusion is strong ($Am=0.1$); the VSI and MRI can co-exist for $Am=1$; and the VSI is overwhelmed by the MRI when ambipolar diffusion is weak ($Am=10$). Angular momentum transport process is primarily driven by MHD winds, while viscous accretion due to MRI and/or VSI turbulence makes a moderate contribution in most cases. Spontaneous magnetic flux concentration and formation of annular substructures remain robust in strong ambipolar diffusion dominated disks ($Am\leq1$) with the presence of the VSI. Ambipolar diffusion is the major contributor to the magnetic flux concentration phenomenon rather than advection.

Low wind and petaling effects, caused by the discontinuous apertures of telescopes, are poorly corrected -- if at all -- by commonly used workhorse wavefront sensors (WFSs). Wavefront petaling breaks the coherence of the point spread function core, splitting it into several side lobes, dramatically shutting off scientific throughput. We demonstrate the re-purposing of non-redundant sparse aperture masking (SAM) interferometers into low-order WFSs complementing the high-order pyramid WFS, on the SCExAO experimental platform at Subaru Telescope. The SAM far-field interferograms formed from a 7-hole mask are used for direct retrieval of petaling aberrations, which are almost invisible to the main AO loop. We implement a visible light dual-band SAM mode, using two disjoint 25 nm wide channels, that we recombine to overcome the one-lambda ambiguity of fringe-tracking techniques. This enables a control over petaling with sufficient capture range yet without conflicting with coronagraphic modes in the near-infrared. We present on-sky engineering results demonstrating that the design is able to measure petaling well beyond the range of a single-wavelength equivalent design.

The LHAASO Collaboration has observed ultrahigh-energy photons up to $1.4$PeV from $12$ $\gamma$-ray Galactic sources. In particular, the $\gamma$-ray spectra of the sources J2226+6057, J1908+0621, J1825-1326 have been published. We investigate the possibility of suggested origin pulsars near the sources as the PeVatrons. The pulsar is described by a rotating magnetic dipole. Assuming protons are uniform distributed out of the light cylinders, we obtain the Lorentz distribution of proton energy spectrum. It is found that the protons around pulsar could be accelerated to PeV at short times. The hadronic $\gamma$-ray spectra of the suggested origin pulsars are in good agreement with the LHAASO observed $\gamma$-ray spectra of the sources J2226+6057, J1908+0621, J1825-1326.

Current models of stellar evolution predict that stars more massive than $\sim6\:$M$_{\odot}$ should have completely depleted all lithium (Li) in their atmospheres by the time when they reach the He core burning phase. Against this, a non-negligible number of red giants with masses $>6\:$M$_{\odot}$ presenting strong Li lines have recently been reported. Motivated by this finding, we have carried out a spectroscopic survey of red supergiants (RSGs) in the Perseus Arm and a selection of young open clusters in the Magellanic Clouds to assess the presence of the Li I $\lambda$6708 doublet line. Based on a sample of >70 objects, close to one third of RSGs in the Perseus Arm display noticeable Li lines, with perhaps a trend towards a lower fraction among more luminous stars. The samples in the Magellanic Clouds are not as large, but hint at a metallicity dependence. Twenty one RSGs in 5 LMC clusters show a very high fraction of Li detection, around 40%. Conversely, 17 RSGs in 5 SMC clusters give only one secure detection. The interpretation of these observational results is not straightforward, but a mechanism for Li production seems most likely. Further characterisation work is ongoing, while theoretical studies into this matter are urgently needed.

This article provides detailed description on the fundamentals of aperture photometry analysis. The differential and all-sky aperture photometry techniques are described thoroughly to depict the difference between the two techniques and their selection for determining the stars' magnitudes and their respective magnitude errors. The crucial calibration parameters required for the all-sky photometry analysis such as atmospheric extinctioncoefficient, air-mass, zero point, color term and color index are discussed comprehensively with their extraction from the Sloan Digital Sky Survey (SDSS) archive. The all-sky aperture photometry technique is applied on the stars of an open cluster NGC 2420 to determine their calibrated magnitudes and magnitude errors in the g, r, and i bands. The images required for the analysis are extracted from data release DR12 of SDSS III archive. Herein, the photometry analysis is performed by the Makali'i: SUBARU Image Processor, a Windows-based software. This software has a simple yet effective GUI and it provides the starlight minus the background sky light value with a single click. This article would aid in providing the insight into the physics of aperture photometry by manually scanning the astronomical images. In addition, the g, r, and i magnitudes are transformed to B, V, and R band magnitudes of Johnson-Cousins UBVRI photometric system. The color magnitude diagram for both the standard photometry systems are also provided.

The fundamentals of estimating essential stellar parameters of an open cluster-NGC 2360 and globular clusters-NGC 5272 are presented extensively in this work. Here, the evaluation of stellar parameters, by manually fitting the appropriate isochrones on the color magnitude diagrams (CMDs), of the selected star clusters is discussed comprehensively. Aperture photometry and PSF fitting photometry are conducted on g, r, and i standard band filter images of Sloan Digital Sky Survey (SDSS) using the aperture photometry tool (APT) to obtain the respective CMDs. Further, to achieve the stellar parameters, isochrone fitting is described in detail. This work on stellar parameters evaluation has attained the following results: age of NGC 2360 is found to be 708 Myrs with metallicity, [Fe/H], of -0.15, whereas NGC 5272 is having age of 11.56 Gyrs with metallicity, [Fe/H], of -1.57. Additionally, the interstellar reddening, E(B-V), and distance modulus, DM, for NGC 2360 are obtained as 0.12 and 11.65, respectively. While, for NGC 5272, the interstellar reddening is attained as E(B-V)=0.015, and the distance modulus is DM=15.1. The values of these stellar parameters are found to be in close approximation with the results provided in the literature based on the IRAF analysis technique. The distribution of radii, masses, and temperatures are included along with the initial mass function (IMF) for both the start clusters. Thus, this article would aid in providing insight into the evaluation of stellar parameters by the astronomical photometry analysis which would successively upsurge the understanding of our universe. It should be noted that the cleaning of cluster population on the CMDs from the foreground/background stars, clearing of spurious objects, error estimations and the membership determination are not carried out in this work and are considered as separate project for analysis.

The recently-computed ExoMol line lists for isotopologues of AlH are used to analyse the blue spectrum (4000-4500 {\AA}) of Proxima Cen (M5.5 V). Comparison of the observed and computed spectra enables the identification of a large number of 27AlH lines of the A1{\Pi} - X1{\Sigma}+ band system: the spectral range covering 1-0, 0-0 and 1-1 bands is dominated by clearly resolved AlH lines. We reveal the diffuse nature of transitions close to the dissociation limit which appears in the form of increasingly wider(up to 5 {\AA}) and shallower (up to the continuum confusion limit) AlH line profiles. The predicted wavelengths of AlH diffuse lines are systematically displaced. The effect broadening by pre-dissociation states on the line profiles is included by increasing the radiative damping rate by up to 5 orders of magnitude. We determine empirical values of damping rates for a number of the clean 0-0 Q-branch transitions by comparing the observed and synthetic stellar spectra. We find excellent agreement between our damping rates and lifetimes available in the literature. A comparison of 27Al1H ExoMol and REALH spectra shows that the observed spectrum is better described by the ExoMol line list. A search for 26Al1H lines in the Proxima Cen spectrum does not reveal any notable features; giving an upper limit of 27Al1H/26Al1H {>} 100.

Observations of radiative cooling in a synchrotron source offer a possibility to further constrain its properties. Inverse Compton cooling is indicated in the radio spectra during the early phases of SN\,2020oi. It is shown that contrary to previous claims, observations are consistent with equipartition between relativistic electrons and magnetic field as well as a constant mass-loss rate of the progenitor star prior to the supernova explosion. The reason for this difference is the need to include cooling directly in the fitting procedure rather than estimating its effects afterward. It is emphasized that the inferred properties of the supernova ejecta are sensitive to the time evolution of the synchrotron self-absorption frequency; hence, great care should be taken when modeling spectra for which cooling and/or inhomogeneities are indicated. Furthermore, it is noted that the energies of the relativistic electrons in the radio emission regions in supernovae are likely too low for first-order Fermi acceleration to be effective.

New classes of astronomical objects are often discovered serendipitously. The enormous data volumes produced by recent high-time resolution, radio-telescope surveys imply that efficient algorithms are required for a discovery. Such algorithms are usually tuned to detect specific, known sources. Existing data sets therefore likely contain unknown astronomical sources, which will remain undetected unless algorithms are developed that can detect a more diverse range of signals. We present the Single-dish PARKES data challenge for finding the uneXpected (SPARKESX), a compilation of real and simulated high-time resolution observations. SPARKESX comprises three mock surveys from the Parkes "Murriyang" radio telescope. A broad selection of simulated and injected expected signals (such as pulsars, fast radio bursts), poorly characterised signals (plausible flare star signatures) and unknown unknowns are generated for each survey. The goal of this challenge is to aid in the development of new algorithms that can detect a wide-range of source types. We show how successful a typical pipeline based on the standard pulsar search software, PRESTO, is at finding the injected signals. The dataset is publicly available at https://doi.org/10.25919/fd4f-0g20.

Protoclusters are important for studying how halo mass and stellar mass assemble in the early universe. Finding signposts of such over-dense regions is a popular method to identify protocluster candidates. Hyperluminous infrared galaxies (HLIRGs), are expected to reside in overdense regions with massive halos. We study the Mpc-scale environment of the largest HLIRG sample to date and investigate whether they predominantly live in overdense regions. We first explore the surface density of Herschel 250 $\mu$m sources around HLIRGs and compare with that around random positions. Then, we compare the spatial distribution of neighbours around HLIRGs with that around randomly selected galaxies using a deep IRAC-selected catalogue with good-quality photometric redshifts. We also use a redshift-matched quasar sample and submillimeter galaxy (SMG) sample to validate our method, as previous clustering studies have measured the host halo masses of these populations. Finally, we adopt a Friends of Friends (FOF) algorithm to seek (proto)clusters that host HLIRGs. We find that HLIRGs tend to have more bright star-forming neighbours (with 250 $\mu$m flux density >10 mJy) within 100$\arcsec$ projected radius than a random galaxy at a 3.7$\sigma$ significance. In our 3D analysis, we find relatively weak excess of IRAC-selected sources within 3 Mpc around HLIRGs compared with random galaxy neighbours, mainly influenced by photometric redshift uncertainty and survey depth. We find a more significant difference (at a 4.7$\sigma$ significance) in the number of Low Frequency Array (LOFAR)-detected neighbours in the deepest EN1 field. HLIRGs at 3 < z < 4 show stronger excess compared to HLIRGs at 2 < z < 3, consistent with cosmic downsizing. Finally, we select and present a list of 30 most promising protocluster candidates for future follow-up observations.

To use strong gravitational lenses as an astrophysical or cosmological probe, models of their mass distributions are often needed. We present a new, time-efficient automation code for uniform modeling of strongly lensed quasars with GLEE, a lens modeling software, for high-resolution multi-band data. By using the observed positions of the lensed quasars and the spatially extended surface brightness distribution of the lensed quasar host galaxy, we obtain a model of the mass distribution of the lens galaxy. We apply this uniform modeling pipeline to a sample of nine strongly lensed quasars with HST WFC 3 images. The models show in most cases well reconstructed light components and a good alignment between mass and light centroids. We find that the automated modeling code significantly reduces the user input time during the modeling process. The preparation time of required input files is reduced significantly. This automated modeling pipeline can efficiently produce uniform models of extensive lens system samples which can be used for further cosmological analysis. A blind test through a comparison with the results of an independent automated modeling pipeline based on the modeling software Lenstronomy reveals important lessons. Quantities such as Einstein radius, astrometry, mass flattening and position angle are generally robustly determined. Other quantities depend crucially on the quality of the data and the accuracy of the PSF reconstruction. Better data and/or more detailed analysis will be necessary to elevate our automated models to cosmography grade. Nevertheless, our pipeline enables the quick selection of lenses for follow-up monitoring and further modeling, significantly speeding up the construction of cosmography-grade models. This is an important step forward to take advantage of the orders of magnitude increase in the number of lenses expected in the coming decade.

Venus and Earth provide an astonishingly different view of the evolution of a rocky planet, raising the question of why these two rocky worlds evolved so differently. The recently discovered transiting rocky planet LP 890-9c (TOI-4306c, SPECULOOS-2c) is a key to this question. SPECULOOS-2c (1.367 +0.055 -0.039 R_Earth) circles a relatively low-activity nearby (32 pc) M6V star in 8.46 days. SPECULOOS-2c receives 0.906 +/- 0.026 of the flux of modern Earth, putting it very close to the inner edge of the conservative Habitable Zone, where models differ strongly in their prediction of how long Earth-like planets can hold onto their water. Our atmosphere models show that the transmission spectra of the observable species can tell the difference between a hot, wet Young Earth, a steamy rocky planet caught in a runaway greenhouse at the brink of complete water loss, and a Venus-analog. Distinguishing these scenarios from the planet's spectra will provide critical new insights into when a hot terrestrial planet loses its water and becomes a Venus. SPECULOOS-2c is a prime target for observations with JWST. Observing it will also provide key insights to predict the long-term future of Earth.

Eclipsing binaries are important benchmark objects to test and calibrate stellar structure and evolution models. This is especially true for binaries with a fully convective M-dwarf component for which direct measurements of these stars' masses and radii are difficult using other techniques. Within the potential of M-dwarfs to be exoplanet host stars, the accuracy of theoretical predictions of their radius and effective temperature as a function of their mass is an active topic of discussion. Not only the parameters of transiting exoplanets but also the success of future atmospheric characterisation rely on accurate theoretical predictions. We present the analysis of five eclipsing binaries with low-mass stellar companions out of a sub-sample of 23, for which we obtained ultra high-precision light curves using the CHEOPS satellite. The observation of their primary and secondary eclipses are combined with spectroscopic measurements to precisely model the primary parameters and derive the M-dwarfs mass, radius, surface gravity, and effective temperature estimates using the PYCHEOPS data analysis software. Combining these results to the same set of parameters derived from TESS light curves, we find very good agreement (better than 1\% for radius and better than 0.2% for surface gravity). We also analyse the importance of precise orbits from radial velocity measurements and find them to be crucial to derive M-dwarf radii in a regime below 5% accuracy. These results add five valuable data points to the mass-radius diagram of fully-convective M-dwarfs.

The next-generation (3G/XG) ground-based gravitational-wave (GW) detectors such as Einstein Telescope (ET) and Cosmic Explorer (CE) will begin observing in the next decade. Due to the extremely high sensitivity of these detectors, the majority of stellar-mass compact-binary mergers in the entire Universe will be observed. It is also expected that 3G detectors will have significant sensitivity down to 2-7 Hz; the observed duration of binary neutron star signals could increase to several hours or days. The abundance and duration of signals will cause them to overlap in time, which may form a confusion noise that could affect the detection of individual GW sources when using naive matched filtering; Matched filtering is only optimal for stationary Gaussian noise. We create mock data for CE and ET using the latest population models informed by the GWTC-3 catalog and investigate the performance loss of matched filtering due to overlapping signals. We find the performance loss mainly comes from a deviation in the noise's measured amplitude spectral density. The redshift reach of CE (ET) can be reduced by 15-38 (8-21) % depending on the merger rate estimate. The direct contribution of confusion noise to the total SNR is generally negligible compared to the contribution from instrumental noise. We also find that correlated confusion noise has a negligible effect on the quadrature summation rule of network SNR for ET, but might reduce the network SNR of high detector-frame mass signals for detector networks including CE if no mitigation is applied. For ET, the null stream can mitigate the astrophysical foreground. For CE, we demonstrate that a computationally efficient, straightforward single-detector signal subtraction method suppresses the total noise to almost the instrument noise level; this will allow for near-optimal searches.

Vortices have long been speculated to play a role in planet formation, via the collection of dust in the pressure maxima that arise at the cores of vortices in protoplanetary discs. The question remains however: as dust collects in the core of a vortex, when does that vortex remain stable and able to collect further dust, and when and why does it break up? We study this question by running high resolution 2D simulations of dust laden vortices. By using the terminal velocity approximation in a local shearing box it was possible to efficiently run simulations of back-reacting dust in a gas at high resolution. Our results show how the stability of 2D dust laden vortices in protoplanetary discs depends on their size relative to the disc scale height, as well as the dust coupling. We find small vortices with semiminor axis much smaller than the scale height to be stable for the duration of the simulations ($t>2000$ orbits). Larger vortices, with semiminor axis smaller than but of order of the scale height, exhibit a drag instability after undergoing a long period of contraction where the core becomes progressively more dust rich. The lifetime of these vortices depends on the dust size, with larger dust grains causing the instability to occur sooner. For the size ranges tested in this paper, $\mu$m to mm sized grains, vortices survived for several hundreds of orbits. The result implies that the stability of vortices formed by vertical shear instability and zombie vortex instability, or the breakup of larger vortices through hydrodynamic instabilities, is affected by the presence of dust in the disc. The lifetimes observed in this paper, while shortened by the presence of dust for larger vortices, were still long enough to lead to considerable dust enrichment in the vortex cores.

MICADO is a workhorse instrument for the ESO ELT, allowing first light capability for diffraction limited imaging and long-slit spectroscopy at near-infrared wavelengths. The PSF Reconstruction (PSF-R) Team of MICADO is currently implementing, for the first time within all ESO telescopes, a software service devoted to the blind reconstruction of the PSF. This tool will work independently of the science data, using adaptive optics telemetry data, both for Single Conjugate (SCAO) and Multi-Conjugate Adaptive Optics (MCAO) allowed by the MORFEO module. The PSF-R service will support the state-of-the-art post-processing scientific analysis of the MICADO imaging and spectroscopic data. We provide here an update of the status of the PSF-R service tool of MICADO, after successfully fulfilling the Final Design Review phase, and discuss recent results obtained on simulated and real data gathered on instruments similar to MICADO.

Current state-of-the-art adaptive optics (AO) provides ground-based, diffraction-limited observations with high Strehl ratios (SR). However, a detailed knowledge of the point spread function (PSF) is required to fully exploit the scientific potential of these data. This is even more crucial for the next generation AO instruments that will equip 30-meter class telescopes, as the characterization of the PSF will be mandatory to fulfill the planned scientific requirements. For this reason, there is a growing interest in developing tools that accurately reconstruct the observed PSF of AO systems, the so-called PSF reconstruction. In this context, a PSF-R service is a planned deliverable for the MICADO@ELT instrument and our group is in charge of its development. In the case of MICADO, a blind PSF-R approach is being pursued to have the widest applicability to science cases. This means that the PSF is reconstructed without extracting information from the science data, relying only on telemetry and calibrations. While our PSF-R algorithm is currently being developed, its implementation is mature enough to test performances with actual observations. In this presentation we will discuss the reliability of our reconstructed PSFs and the uncertainties introduced in the measurements of scientific quantities for bright, on-axis observations taken with the SOUL+LUCI instrument of the LBT. This is the first application of our algorithm to real data. It demonstrates its readiness level and paves the way to further testing. Our PSF-R algorithm is able to reconstruct the SR and full-width at half maximum of the observed PSFs with errors smaller than 2% and 4.5%, respectively. We carried out the scientific evaluation of the obtained reconstructed PSFs thanks to a dedicated set of simulated observations of an ideal science case.

We investigate the internal kinematics of the young star-forming region NGC 346 in the Small Magellanic Cloud. We used two epochs of deep F555W and F814W Hubble Space Telescope ACS observations with an 11-year baseline to determine proper motions, and study the kinematics of different populations, as identified by their color-magnitude diagram and spatial distribution characteristics. The proper motion field of the young stars shows a complex structure with spatially coherent patterns. NGC 346 upper-main sequence and pre-main sequence stars follow very similar motion patterns, with the outer parts of the cluster being characterized both by outflows and inflows. The proper motion field in the inner ~10 pc shows a combination of rotation and inflow, indicative of inspiraling motion. The rotation velocity in this regions peaks at ~3 km/s, whereas the inflow velocity peaks at ~1 km/s. Sub-clusters and massive young stellar objects in NGC 346 are found at the interface of significant changes in the coherence of the proper motion field. This suggests that turbulence is the main star formation driver in this region. The similar kinematics observed in the metal-poor NGC 346 and the Milky Way star-forming regions suggest that the differences in the cooling conditions due to the different amounts of metallicity and dust density between the SMC and our Galaxy are too small to alter significantly the process of star clusters assembly and growth. The main characteristics of our findings are consistent with various proposed star cluster formation models.

The redshift space anisotropy of the bispectrum is generally quantified using multipole moments. The possibility of measuring these multipoles in any survey depends on the level of statistical fluctuations. We present a formalism to compute the statistical fluctuations in the measurement of bispectrum multipoles for galaxy surveys. We consider specifications of a {\it Euclid} like galaxy survey and present two quantities: the signal-to-noise ratio (SNR) which quantifies the detectability of a multipole, and the rank correlation which quantifies the correlation in measurement errors between any two multipoles. Based on SNR values, we find that {\it Euclid} can potentially measure the bispectrum multipoles up to $\ell=4$ across various triangle shapes, formed by the three {\bf k} vectors in Fourier space. In general, SNR is maximum for the linear triangles. SNR values also depend on the scales and redshifts of observation. While, $\ell \leq 2$ multipoles can be measured with ${\rm SNR}>5$ even at linear/quasi-linear ($k \lesssim 0.1 \,{\rm Mpc}^{-1}$) scales, for $\ell>2$ multipoles, we require to go to small scales or need to increase bin sizes. For most multipole pairs, the errors are only weakly correlated across much of the triangle shapes barring a few in the vicinity of squeezed and stretched triangles. This makes it possible to combine the measurements of different multipoles to increase the effective SNR.

We present the stellar radial velocity analysis of the central 1x1' of the young massive Small Magellanic Cloud star cluster NGC 346. Using VLT/MUSE integral field spectroscopy in combination with Hubble Space Telescope photometry we extract 103 spectra of cluster member stars suited to measure accurate line-of-sight kinematics. The cluster member stars show two distinct velocity groups at v1 = -3.3 (+0.3/-0.2) km/s and v2 = 2.6 (+0.1/-0.1) km/s, relative to the systemic velocity of 165.5+/-0.2 km/s, and hint for a third group at v3 = 9.4 (+0.1/-0.1 km/s. We show that there is neither a correlation between the velocity groups and the spatial location of the stars, nor their locus on optical color-magnitude diagrams, which makes the stellar velocity a key parameter to separate individual stellar components in such a young star cluster. Velocity group 2 shows clear rotation with Omega2 = -0.4 +/- 0.1 1/Myr, corresponding to -4.9+/-0.7 km/s at radial distance of 10 pc from the center, a possible remnant of the formation process of NGC 346 through the hierarchical collapse of the giant molecular cloud. The ionizing gas has lost any natal kinematic imprint and shows clear expansion, driven by far ultra violet fluxes and stellar winds of the numerous OB stars in the cluster center. The size of this expanding bubble and its expansion velocity of 7.9 km/s is in excellent agreement with the estimate that the latest star formation episode occurred about two million years ago.

Numerical models of solar flares typically focus on the behaviour of directly-heated flare models, adopting magnetic field- aligned, plane-parallel methodologies. With high spatial- and spectral-resolution ground-based optical observations of flares, it is essential also to understand the response of the plasma surrounding these strongly heated volumes. We investigate the effects of the extreme radiation field produced by a heated column of flare plasma on an adjacent slab of chromospheric plasma, using a two-dimensional radiative transfer model and considering the time-dependent solution to the atomic level populations and electron density throughout this model. The outgoing spectra of H$\alpha$ and Ca II 854.2 nm synthesised from our slab show significant spatial-, time-, and wavelength-dependent variations (both enhancements and reductions) in the line cores, extending on order 1 Mm into the non-flaring slab due to the incident transverse radiation field from the flaring boundary. This may lead to significant overestimates of the sizes of directly-heated flare kernels, if line-core observations are used. However, the radiation field alone is insufficient to drive any significant changes in continuum intensity, due to the typical photospheric depths at which they forms, so continuum sources will not have an apparent increase in size. We show that the line formation regions near the flaring boundary can be driven upwards in altitude by over 1 Mm despite the primary thermodynamic parameters (other than electron density) being held horizontally uniform. This work shows that in simple models these effects are significant and should be considered further in future flare modelling and interpretation.

Magnetic reconnection is ubiquitous in astrophysical systems, and in many such systems, the plasma suffers from significant cooling due to synchrotron radiation. We study relativistic magnetic reconnection in the presence of strong synchrotron cooling, where the ambient magnetization $\sigma$ is high and the magnetic compactness $\ell_{B}$ of the system is of order unity. In this regime, $e^{\pm}$ pair production from synchrotron photons is inevitable, and this process can regulate the magnetization $\sigma$ surrounding the current sheet. We investigate this self-regulation analytically and find a self-consistent steady state for a given magnetic compactness of the system and initial magnetization. This result helps estimate the self-consistent upstream magnetization in systems where plasma density is poorly constrained, and can be useful for a variety of astrophysical systems. As illustrative examples, we apply it to study the properties of reconnecting current sheets near the supermassive black hole of M87, as well as the equatorial current sheet outside the light cylinder of the Crab pulsar.

The astronomical applications greatly benefit from the knowledge of the instrument PSF. We describe the PSF Reconstruction algorithm developed for the LBT LUCI instrument assisted by the SOUL SCAO module. The reconstruction procedure considers only synchronous wavefront sensor telemetry data and a few asynchronous calibrations. We do not compute the Optical Transfer Function and corresponding filters. We compute instead a temporal series of wavefront maps and for each of these the corresponding instantaneous PSF. We tested the algorithm both in laboratory arrangement and in the nighttime for different SOUL configurations, adapting it to the guide star magnitudes and seeing conditions. We nick-named it "BRUTE", Blind Reconstruction Using TElemetry, also recalling the one-to-one approach, one slope-to one instantaneous PSF the algorithm applies.

We examine how the fraction $f$ of stars for which rotational modulation has been detected in Kepler light curves depends on the stellar mass $M_\star$ and age $t_\star$. Our sample consists of $\approx 850$ FGK stars hosting transiting planet candidates detected from the prime Kepler mission. For these stars, atmospheric parameters have been derived using high-resolution spectra from the California-Kepler survey, and rotational modulation has been searched in Kepler light curves homogeneously. We fit stellar models to the atmospheric parameters, Gaia parallax, and 2MASS magnitude of these stars and obtain samples drawn from the posterior probability distributions for their masses and ages under a given, uninformative prior. We combine them with the result of rotational modulation search to simultaneously infer the mass-age distribution of the sample as well as $f(M_\star, t_\star)$, in a manner that fully takes into account mass and age uncertainties of individual stars. We find that $f$ remains near unity up to $t_\star \sim 3\,\mathrm{Gyr}$ and drops to almost zero by $t_\star \sim 5\,\mathrm{Gyr}$, although the trend is less clearly detected for stars with $\lesssim 0.9\,M_\odot$ due to weaker age constraints. This finding is consistent with a view that the detection of rotational modulation is limited by photometric precision to younger stars that exhibit higher-amplitude modulation, and suggests that the detectability of rotational modulation in Kepler light curves is insensitive to metallicity and activity cycles for stars younger than the Sun.

We present the first detailed analysis of the radio halo in the merging galaxy cluster Abell 2256 using the LOFAR, uGMRT, and VLA. These observations combined with archival X-ray data allowed us to study the halo emission with unprecedented detail. The integrated radio emission from the entire halo is characterized by an ultra-steep spectrum, which can be described by a power law with $\alpha_{144 \rm MHz}^{1.5 \rm GHz}=-1.63\pm0.03$, and a radial steepening in the outer regions. The halo is significantly underluminous according to the scaling relations between radio power and mass at 1.4 GHz but not at 150 MHz; ultra-steep spectrum halos are predicted statistically underluminous. Despite the complex structure of this system, the radio halo morphology is remarkably similar to that of the X-ray emission. The radio surface brightness distribution across the halo is strongly correlated with the X-ray brightness of the intracluster medium. The derived correlations show sublinear slopes and there are distinct structures: the core is $\rm I_{R}\propto I_{X}^{1.51}$, the outermost region $\rm I_{R}\propto I_{X}^{0.41}$, and we find radio morphological connections with X-ray discontinuities. We also find a strong anti-correlation between the radio spectral index and the X-ray surface brightness, implying radial steepening. We suggests that the halo core is either related to old plasma from previous AGN activity, being advected, compressed and re-accelerated by mechanisms activated by the cold front or the turbulence is less strong and magnetic field is strong in the core. The change in the radio vs X-ray correlation slopes in the outer regions of the halo could be due to a radial decline of magnetic field, increase in the number density of seed particles or increasing turbulence. Our findings suggest that that the emitting volume is not homogenous according to turbulence re-acceleration models.

Gas accretion is an important process in the evolution of galaxies, but it has limited direct observational evidences. In this paper, we report the detection of a possible ongoing gas accretion event in a Blue Compact Dwarf (BCD) galaxy, MaNGA 8313-1901, observed by the Mapping Nearby Galaxies and Apache Point Observatory (MaNGA) program. This galaxy has a distinct off-centered blue clump to the northeast (the NE clump) that shows low metallicity and enhanced star-formation. The kinematics of the gas in the NE clump also seems to be detached from the host BCD galaxy. Together with the metallicity drop of the NE clump, it suggests that the NE clump likely has an external origin, such as the gas accretion or galaxy interaction, rather than an internal origin, such as an \hii~complex in the disk. After removing the underlying host component, we find that the spectrum of the "pure" clump can match very well with a modeled spectrum containing a stellar population of the young stars ($\le 7$ Myr) only. This may imply that the galaxy is experiencing an accretion of cold gas, instead of a merger event involving galaxies with significant pre-existing old stars. We also find signs of another clump (the SW clump) at the south-west corner of the host galaxy, and the two clumps may share the same origin of gas accretion.

Smalls scale challenges suggest some missing pieces in our current understandings of dark matter. A cascade theory for dark matter flow is proposed to provide extra insights, similar to the cascade in hydrodynamic turbulence. The energy cascade from small to large scales with a constant rate $\varepsilon_u$ ($\approx -4.6\times 10^{-7}m^2/s^3$) is a fundamental feature of dark matter flow. Energy cascade leads to a two-thirds law for kinetic energy $v_r^2$ on scale $r$ such that $v_r^2 \propto (\varepsilon_u r)^{2/3}$, as confirmed by N-body simulations. This is equivalent to a four-thirds law for mean halo density $\rho_s$ enclosed in the scale radius $r_s$ such that $\rho_s \propto \varepsilon_u^{2/3}G^{-1}r_s^{-4/3}$, as confirmed by data from galaxy rotation curves. By identifying relevant key constants, critical scales of dark matter might be obtained. The largest halo scale $r_l$ can be determined by $-u_0^3/\varepsilon_u$, where $u_0$ is the velocity dispersion. The smallest scale $r_{\eta}$ is dependent on the nature of dark matter. For collisionless dark matter, $r_{\eta} \propto (-{G\hbar/\varepsilon_{u}}) ^{1/3}\approx 10^{-13}m$, where $\hbar$ is the Planck constant. A uncertainty principle for momentum and acceleration fluctuations is also postulated. For self-interacting dark matter, $r_{\eta} \propto \varepsilon_{u}^2 G^{-3}(\sigma/m)^3$, where $\sigma/m$ is the cross-section. On halo scale, the energy cascade leads to an asymptotic slope $\gamma=-4/3$ for fully virialized halos with a vanishing radial flow, which might explain the nearly universal halo density. Based on the continuity equation, halo density is analytically shown to be closely dependent on the radial flow and mass accretion such that simulated halos can have different limiting slopes. A modified Einasto density profile is proposed accordingly.

We present a set of three-dimensional, global, general relativistic radiation magnetohydrodynamic simulations of thin, radiation-pressure-dominated accretion disks surrounding a non-rotating, stellar-mass black hole. The simulations are initialized using the Shakura-Sunyaev model with a mass accretion rate of $\dot{M} = 3 L_\mathrm{Edd}/c^2$ (corresponding to $L=0.17 L_\mathrm{Edd}$). Our previous work demonstrated that such disks are thermally unstable when accretion is driven by an $\alpha$-viscosity. In the present work, we test the hypothesis that strong magnetic fields can both drive accretion through the magneto-rotational instability and restore stability to such disks. We test four initial magnetic field configurations: 1) a zero-net-flux case with a single, radially extended set of magnetic field loops (dipole); 2) a zero-net-flux case with two radially extended sets of magnetic field loops of opposite polarity stacked vertically (quadrupole); 3) a zero-net-flux case with multiple radially concentric rings of alternating polarity (multi-loop); and 4) a net-flux, vertical magnetic field configuration (vertical). In all cases, the fields are initially weak, with the gas-to-magnetic pressure ratio $\gtrsim 100$. Based on the results of these simulations, we find that the dipole and multi-loop configurations remain thermally unstable like their $\alpha$-viscosity counterpart, in our case collapsing vertically on the local thermal timescale and never fully recovering. The vertical case, on the other hand, stabilizes and remains so for the duration of our tests (many thermal timescales). The quadrupole case is intermediate, showing signs of both stability and instability. The key stabilizing criteria is, $P_\mathrm{mag} \gtrsim 0.5P_\mathrm{tot}$ with strong toroidal fields near the disk midplane. We also report a comparison of our models to the standard Shakura-Sunyaev disk.

Observing the list of compatible second order equations of Absolute Parallelism (AP) found by Einstein and Mayer (they used D=4), we choose the one-parameter class of equations which take on a 3-linear form (when contra-frame density of some weight is in use). Spherically symmetric solutions to these equations are considered, and we try not to add any delta-sources (ie, $\delta(x)$-sources of unknown nature) during integrations allowed due to this high symmetry. Using two different ways to fix the radius and time, we have found that only non-static solutions (except for trivial one, of course) are possible. If D=5, such solutions, looking like a single wave moving along the radius, could serve as an expanding cosmological model (with a simple Hubble diagram). With one coordinate choice (gauge), a single second order equation remains and there exist spherically symmetric solutions with arising singularities. On the other hand, a more reasonable (covariant) choice of the radius and time reduces the problem to a system of two first order equations looking like Chaplygin gas dynamics, where solutions are seemingly free of emerging singularities and gradient catastrophe.

Recent sky surveys have discovered a large number of stellar substructures. It is highly likely that there are dark matter (DM) counterparts to these stellar substructures. We examine the implications of DM substructures for electron recoil (ER) direct detection (DD) rates in dual phase xenon experiments. We have utilized the results of the LAMOST survey and considered a few benchmark substructures in our analysis. Assuming that these substructures constitute $\sim 10\%$ of the local DM density, we study the discovery limits of DM-electron scattering cross sections considering one kg-year exposure and 1, 2, and 3 electron thresholds. With this exposure and threshold, it is possible to observe the effect of the considered DM substructure for the currently allowed parameter space. We also explore the sensitivity of these experiments in resolving the DM substructure fraction. For all the considered cases, we observe that DM having mass $\mathcal{O}(10)\,$MeV has a better prospect in resolving substructure fraction as compared to $\mathcal{O}(100)\,$MeV scale DM. We also find that within the currently allowed DM-electron scattering cross-section; these experiments can resolve the substructure fraction (provided it has a non-negligible contribution to the local DM density) with good accuracy for $\mathcal{O}(10)\,$MeV DM mass with one electron threshold.

A simple cooling model of white dwarf stars is re-analyzed in Palatini $f(R)$ gravity. Modified gravity affects the white dwarf structures and consequently their ages. We find that the resulting super-Chandrasekhar white dwarfs need more time to cool down than sub-Chandrasekhar ones, or when compared to the Newtonian models.

Invertible disformal transformations are a useful tool to investigate ghost-free scalar-tensor theories. By performing a higher-derivative generalization of the invertible disformal transformation on Horndeski theories, we construct a novel class of ghost-free scalar-tensor theories, which we dub generalized disformal Horndeski theories. Specifically, these theories lie beyond the quadratic/cubic DHOST class. We explore cosmological perturbations to identify a subclass where gravitational waves propagate at the speed of light and clarify the conditions for the absence of ghost/gradient instabilities for tensor and scalar perturbations. We also investigate the conditions under which a matter field can be consistently coupled to these theories without introducing unwanted extra degrees of freedom.

Neutrino emission in coincidence with gamma rays has been observed from the blazar TXS 0506+056 by the IceCube telescope. Neutrinos from the blazar had to pass through a dense spike of dark matter (DM) surrounding the central black hole. The observation of such a neutrino implies new upper bounds on the neutrino-DM scattering cross section as a function of DM mass. The constraint is stronger than existing ones for a range of DM masses, if the cross section rises linearly with energy. For constant cross sections, competitive bounds are also possible, depending on details of the DM spike.

Inference of the equation-of-state (EoS) of dense nuclear matter in neutron-star cores is a principal science goal of X-ray and gravitational-wave observations of neutron stars. In particular, gravitational-wave observations provide an independent probe of the properties of bulk matter in neutron star cores that can then be used to compare with theoretically derived equations of state. In this paper, we quantify the systematic errors arising from the application of EoS-independent \emph{quasi-universal relations} in the estimation of neutron star tidal deformabilities and radii from gravitational-wave measurements and introduce a strategy to correct for the systematic biases in the inferred radii. We apply this method to a simulated population of events expected to be observed by future upgrades of current detectors and the next-generation of ground-based observatories. We show that our approach can accurately correct for the systematic biases arising from approximate universal relations in the mass-radius curves of neutron stars. Using the posterior distributions of the mass and radius for the simulated population we infer the underlying EoS with a good degree of precision. Our method revives the possibility of using the universal relations for rapid Bayesian model selection of dense matter EoS in gravitational-wave observations.

For many students, introductory college science courses are often the only opportunity in their formal higher education to be exposed to science, shaping their view of the subject, their scientific literacy, and their attitudes towards their own ability in STEM. While science writing instruction has been demonstrated to impact attitudes and outlooks of STEM majors in their coursework, this instructional strategy has yet to be explored for non-majors. In this work, we investigate student attitudes towards STEM before and after taking a writing-intensive introductory astronomy course. We find that students cite writing about science as beneficial to their learning, deepening their understanding of science topics and their perspective on science as a field and finding writing to be a "bridge" between STEM content and their focus on humanities in their majors. Students also report increased perceptions of their own ability and confidence in engaging with STEM across multiple metrics, leaving the course more prepared to be informed, engaged, and science literate citizens.

Left-handed neutrinos interact attractively by Z-boson exchange. The Ginzburg-Landau mean field calculation and the Bogoliubov transformation suggest that the attractive force leads to neutrino pair condensate and neutrino super-fluidity. Neutrinos, as defined by quasi-particle in the super-fluid phase, behave as massless fermions. When the result of super-fluid formation is applied to the early universe, horizon scale pair condensate may become a component of dark energy. A further accretion of other fermions from thermal cosmic medium gives a seed of primordial neutron stars consisting of proton, neutron, electron, and neutrino in beta-equilibrium, surrounded by left-handed neutrino pair condensate. This possibility may provide a mechanism of giving a part or the whole of the dark matter in the present universe, if a properly chosen small fraction condenses to neutrino super-fluid and primordial neutron star not to over-close the universe. The proposal can be verified by measuring neutrino burst at primordial neutron star formation and by detecting super-fluid relic neutrinos in atomic experiments at laboratories.

It is known that axion haloscopes that operate to search for dark matter axions via the 2-photon anomaly are also sensitive to gravitational waves (GWs) through the inverse Gertsenshtein effect. Recently this way of searching for high frequency GWs has gained momentum as it has been shown that the strain sensitivity of such detectors, h_g, are of the same order of sensitivity as the axion-photon theta angle, \theta_a, which is related to the axion 2-photon coupling, g_{a\gamma\gamma}, by, \theta_a = g_{a\gamma\gamma}a, where, a, is the axion scalar field. This means after calculating the sensitivity of a haloscope to an axion signal, we also have calculated the order of magnitude sensitivity to a GW signal of the same spectral and temporal form. However, it is unlikely that a GW and an axion signal will be of the same form since physically the way the signals are generated are completely different. For GW detection, the spectral strain sensitivity in units strain per square root Hz, and is the natural way to compare the sensitivity of GW detectors due to its independence on the GW signal. Likewise, one can define a spectral axion-photon theta angle sensitivity in units of theta angle per square root Hz for axion detectors, which is independent of the axion signal. In this work we introduce a systematic way to calculate the spectral sensitivity of an axion haloscope so instrument comparison may be achieved independent of signal assumptions and only depends on the axion to signal transduction sensitivity and noise in the instrument. Thus, the calculation of the spectral sensitivity not only allows the comparison of dissimilar axion detectors independent of signal, but also allows comparison of the GW sensitivity in terms of spectral strain sensitivity, allowing comparisons to standard GW detectors based on optical interferometers and resonant-mass technology.

IceCube, a cubic-kilometer array of optical sensors built to detect atmospheric and astrophysical neutrinos between 1 GeV and 1 PeV, is deployed 1.45 km to 2.45 km below the surface of the ice sheet at the South Pole. The classification and reconstruction of events from the in-ice detectors play a central role in the analysis of data from IceCube. Reconstructing and classifying events is a challenge due to the irregular detector geometry, inhomogeneous scattering and absorption of light in the ice and, below 100 GeV, the relatively low number of signal photons produced per event. To address this challenge, it is possible to represent IceCube events as point cloud graphs and use a Graph Neural Network (GNN) as the classification and reconstruction method. The GNN is capable of distinguishing neutrino events from cosmic-ray backgrounds, classifying different neutrino event types, and reconstructing the deposited energy, direction and interaction vertex. Based on simulation, we provide a comparison in the 1-100 GeV energy range to the current state-of-the-art maximum likelihood techniques used in current IceCube analyses, including the effects of known systematic uncertainties. For neutrino event classification, the GNN increases the signal efficiency by 18% at a fixed false positive rate (FPR), compared to current IceCube methods. Alternatively, the GNN offers a reduction of the FPR by over a factor 8 (to below half a percent) at a fixed signal efficiency. For the reconstruction of energy, direction, and interaction vertex, the resolution improves by an average of 13%-20% compared to current maximum likelihood techniques in the energy range of 1-30 GeV. The GNN, when run on a GPU, is capable of processing IceCube events at a rate nearly double of the median IceCube trigger rate of 2.7 kHz, which opens the possibility of using low energy neutrinos in online searches for transient events.

We study the multi-fractional theory with $q$-derivatives, where the multi-fractional measure is considered to be in the time direction. The evolution of power spectra and also the expansion history of the universe are investigated in the $q$-derivatives theory. According to the matter power spectra diagrams, the structure growth would increase in the multi-fractional model, expressing incompatibility with low redshift measurements of large scale structures. Furthermore, concerning the diagrams of Hubble parameter evolution, there is a reduction in the value of Hubble constant which conflicts with local cosmological constraints. We also explore the multi-fractional model with current observational data, principally Planck 2018, weak lensing, supernovae, baryon acoustic oscillations (BAO), and redshift-space distortions (RSD) measurements. Numerical analysis reveals that the degeneracy between multi-fractional parameters makes them remain unconstrained under observations.

Non-axisymmetrical deformations of the crust on rapidly rotating neutron stars are one of the main targets of searches for continuous gravitational waves. The maximum ellipticity, or fractional difference in moments of inertia, that can be supported by deformations of the crust (known as "mountains") provides an important upper limit on the strength of these continuous gravitational wave sources. We use the formalism of Gittins et al 2021, along with a deforming force that acts mainly in the transverse direction, to obtain a maximum ellipticity of 7.4$\times$10$^{-6}$. This is larger than the original results of Gittins et al 2021 but consistent with earlier calculations by Ushomirsky et al 2000. This suggests that rotating neutron stars could be strong sources of continuous gravitational waves.

In this work, the current stability is discussed for cosmic strings with the bosonic superconductivity. A non-vanishing curvature of string generally induce the quantum instability of the current-carrying particle. Its decay rates are explored for various types of model parameters, curved string shapes, and decay processes. As a cosmological application, the stability is examined for superconducting strings in the string network and also for cosmic vortons by evaluating their cosmological evolution. The zero mode and hence the vorton cannot be stable in various cases, e.g., with a hierarchy between the current-carrying particle mass off the string and the string tension or with sizable couplings of the current-carrying particle to light species such as the Standard Model particles.

We address the question of the role of low-energy nuclear physics data in constraining neutron star global properties, e.g., masses, radii, angular momentum, and tidal deformability, in the absence of a phase transition in dense matter. To do so, we assess the capacity of 415 relativistic mean field and non-relativistic Skyrme-type interactions to reproduce the ground state binding energies, the charge radii and the giant monopole resonances of a set of spherical nuclei. The interactions are classified according to their ability to describe these characteristics and we show that a tight correlation between the symmetry energy and its slope is obtained providing $N=Z$ and $N\ne Z$ nuclei are described with the same accuracy (mainly driven by the charge radius data). By additionally imposing the constraints from isobaric analog states and neutron skin radius in $^{208}$Pb, we obtain the following estimates: $E_{sym,2}=31.8\pm 0.7$ MeV and $L_{sym,2}=58.1\pm 9.0$ MeV. We then analyze predictions of neutron star properties and we find that the 1.4$M_\odot$ neutron star (NS) radius lies between 12 and 14 km for the "better" nuclear interactions. We show that i) the better reproduction of low-energy nuclear physics data by the nuclear models only weakly impacts the global properties of canonical mass neutron stars and ii) the experimental constraint on the symmetry energy is the most effective one for reducing the uncertainties in NS matter. However, since the density region where constraints are required are well above densities in finite nuclei, the largest uncertainty originates from the density dependence of the EDF, which remains largely unknown.

Ultralight axion-like particles can contribute to the dark matter near the Sun, leading to a distinct, stochastic signature in terrestrial experiments. We search for such particles through their neutron-spin coupling by re-analyzing approximately 40 days of data from a K-$^3$He co-magnetometer with a new frequency-domain likelihood-based formalism that properly accounts for stochastic effects over all axion coherence times relative to the experimental time span. Assuming that axions make up all of the dark matter in the Sun's vicinity, we find a median 95% upper limit on the neutron-spin coupling of $2.4 \times 10^{-10}$ GeV$^{{-1}}$ for axion masses from 0.4 to 4 feV, which is about five orders of magnitude more stringent than previous laboratory bounds in that mass range. Although several peaks in the experiment's magnetic power spectrum suggest the rejection of a white-noise null hypothesis, further analysis of their lineshapes yields no positive evidence for a dark matter axion.