Papers for Friday, Jan 12 2024

Using images at multiple mid-infrared wavelengths, acquired in May 2018 using the VISIR instrument on ESO's Very Large Telescope (VLT), we study Jupiter's pole-to-pole thermal, chemical and aerosol structure in the troposphere and stratosphere. We confirm that the pattern of cool and cloudy anticyclonic zones and warm cloud-free cyclonic belts persists throughout the mid-latitudes, up to the polar boundaries, and evidence a strong correlation with the vertical maximum windshear and the locations of Jupiter's zonal jets. At high latitudes, VISIR images reveal a large region of mid-infrared cooling poleward $\sim$64$^{\circ}$N and $\sim$67$^{\circ}$S extending from the upper troposphere to the stratosphere, co-located with the reflective aerosols observed by JunoCam, and suggesting that aerosols play a key role in the radiative cooling at the poles. Comparison of zonal-mean thermal properties and high-resolution visible imaging from Juno allows us to study the variability of atmospheric properties as a function of altitude and jet boundaries, particularly in the cold southern polar vortex. However, the southern stratospheric polar vortex is partly masked by a warm mid-infrared signature of the aurora. Co-located with the southern main auroral oval, this warming results from the auroral precipitation and/or joule heating which heat the atmosphere and thus cause a significant stratospheric emission. This high emission results from a large enhancement of both ethane and acetylene in the polar region, reinforcing the evidence of enhanced ion-related chemistry in Jupiter's auroral regions.

The Galactic electron density model NE2001 describes the multicomponent ionized structure of the Milky Way interstellar medium. NE2001 forward models the dispersion and scattering of compact radio sources, including pulsars, fast radio bursts, AGNs, and masers, and the model is routinely used to predict the distances of radio sources lacking independent distance measures. Here we present the open-source package NE2001p, a fully Python implementation of NE2001. The model parameters are identical to NE2001 but the computational architecture is optimized for Python, yielding small (<1%) numerical differences between NE2001p and the Fortran code. NE2001p can be used on the command-line and through Python scripts available on PyPI. Future package releases will include modular extensions aimed at providing short-term improvements to model accuracy, including a modified thick disk scale height and additional clumps and voids. This implementation of NE2001 is a springboard to a next-generation Galactic electron density model now in development.

We report observations of CO($J$=9$\to$8) and OH$^+$($N$=1$\to$0) toward the four millimeter-selected lensed starburst galaxies SPT 2354-58, 0150-69, 0314-44, and 0452-50, using the Atacama Large Millimeter/submillimeter Array Compact Array (ALMA/ACA), as part of a larger study of OH$^+$ in the early universe. In this work, we use these observations for the main purpose of spectroscopic redshift measurements. For all sources except 0452-50, we confirm the previously reported most likely redshifts, and we find typical CO and OH$^+$ properties for massive starbursts. For 0452-50, we rule out the previously reported value of $z$=2.0105, measuring a firm redshift of $z$=5.0160 based on [OI], [CII], H$_2$O, and CO emission instead when adding in ancillary ALMA data. Previously, 0452-50 was considered an outlier in relations between dust temperature, far-infrared luminosity and redshift, which may have hinted at an unusually cold starburst with a dust temperature of only $T_{\rm dust}$=(21$\pm$2) K. Instead, our new measurements suggest it to be among highly luminous massive dusty starbursts at $z$$>$5, with rather typical properties within that population. We find a revised dust temperature of $T_{\rm dust}$=(76.2$\pm$2.5) K, and updated lensing-corrected far-infrared and infrared luminosities of (2.35$^{+0.09}_{-0.08}$) and (4.07$^{+0.25}_{-0.27}$)$\times$10$^{13}$ $L_\odot$, respectively - i.e., about an order of magnitude higher than previously reported. We thus do not find evidence for the existence of unusually cold starburst galaxies in the early universe that were missed by previous selection techniques.

We report the All-Sky Automated Survey for SuperNovae discovery of the tidal disruption event (TDE) ASASSN-23bd (AT 2023clx) in NGC 3799, a LINER galaxy with no evidence of strong AGN activity over the past decade. With a redshift of $z = 0.01107$ and a peak UV/optical luminosity of $(5.4\pm0.4)\times10^{42}$ erg s$^{-1}$, ASASSN-23bd is the lowest-redshift and least-luminous TDE discovered to date. Spectroscopically, ASASSN-23bd shows H$\alpha$ and He I emission throughout its spectral time series, and the UV spectrum shows nitrogen lines without the strong carbon and magnesium lines typically seen for AGN. Fits to the rising ASAS-SN light curve show that ASASSN-23bd started to brighten on MJD 59988$^{+1}_{-1}$, $\sim$9 days before discovery, with a nearly linear rise in flux, peaking in the $g$ band on MJD $60000^{+3}_{-3}$. Scaling relations and TDE light curve modelling find a black hole mass of $\sim$10$^6$ $M_\odot$, which is on the lower end of supermassive black hole masses. ASASSN-23bd is a dim X-ray source, with an upper limit of $L_{0.3-10\,\mathrm{keV}} < 1.0\times10^{40}$ erg s$^{-1}$ from stacking all \emph{Swift} observations prior to MJD 60061, but with soft ($\sim 0.1$ keV) thermal emission with a luminosity of $L_{0.3-2 \,\mathrm{keV}}\sim4\times10^{39}$ erg s$^{-1}$ in \emph{XMM-Newton} observations on MJD 60095. The rapid $(t < 15$ days) light curve rise, low UV/optical luminosity, and a luminosity decline over 40 days of $\Delta L_{40}\approx-0.7$ make ASASSN-23bd one of the dimmest TDEs to date and a member of the growing ``Low Luminosity and Fast'' class of TDEs.

Aims. We aim to expand the atmospheric exploration of AF Lep b by modeling all available observations obtained with SPHERE at VLT (between 0.95-1.65, at 2.105, and 2.253 $\mu$m, and NIRC2 at Keck (at 3.8 $\mu$m) with self-consistent atmospheric models. Methods. To understand the physical properties of this exoplanet, we used ForMoSA. This forward-modeling code compares observations with grids of pre-computed synthetic atmospheric spectra using Bayesian inference methods. We used Exo-REM, an atmospheric radiative-convective equilibrium model, including the effects of non-equilibrium processes and clouds. Results. From the atmospheric modeling we derive solutions at a low effective temperature of ~750 K. Our analysis also favors a metal-rich atmosphere (>0.4) and solar to super-solar carbon-to-oxygen ratio (~0.6). We tested the robustness of the estimated values for each parameter by cross-validating our models using the leave-one-out strategy, where all points are used iteratively as validation points. Our results indicate that the photometry point at 3.8 $\mu$m strongly drives the metal-rich and super-solar carbon-to-oxygen solutions. Conclusions. Our atmospheric forward-modeling analysis strongly supports the planetary nature of AF Lep b. Its spectral energy distribution is consistent with that of a young, cold, early-T super-Jovian planet. We recover physically consistent solutions for the surface gravity and radius, which allows us to reconcile atmospheric forward modeling with evolutionary models, in agreement with the previously published complementary analysis done by retrievals. Finally, we identified that future data at longer wavelengths are mandatory before concluding about the metal-rich nature of AF Lep b.

Dark matter may be detected in X-ray decay, including from the decay of the dark matter particles that make up the Milky Way (MW) halo. We use a range of density profiles to compute X-ray line intensity profiles, with a focus on the resonantly produced sterile neutrino dark matter candidate. Compared to the Navarro--Frenk--White density profile, we show that using an adiabatically contracted halo profile suppresses the line intensity in the halo outskirts and enhances it in the Galactic Centre (GC), although this enhancement is eliminated by the likely presence of a core within 3~kpc. Comparing our results to MW halo observations, other X-ray observations, and structure formation constraints implies a sterile neutrino mixing angle parameter $s_{11}\equiv\sin^{2}(2\theta)\times10^{11}\sim[3,4]$ (particle lifetime $\tau_{28}\equiv\tau/(10^{28}\mathrm{sec})\sim[1.0,1.3]$), which is nevertheless is strong tension with some reported non-detections. We make predictions for the likely decay flux that the XRISM satellite would measure in the GC, plus the Virgo and Perseus clusters, and outline further steps to determine whether the dark matter is indeed resonantly produced sterile neutrinos as detected in X-ray decay.

Rest-UV spectroscopy can constrain properties of the stellar populations, outflows, covering fractions, and can indirectly constrain the Lyman continuum escape fraction of galaxies. Many works have studied the rest-UV spectra of more massive star forming galaxies and low-mass galaxies selected via strong nebular line emission or via Ly$\alpha$ emission. However, studies of rest-UV spectroscopy have yet to be done on an unbiased sample at low mass during the epoch of peak star formation ($z\sim2$). We present a stacked rest-UV spectrum of a complete sample of 16 dwarf galaxies ($\rm \langle log(M^*/M_\odot)\rangle_{median} = 8.2$) at $z\sim2$. The rest-UV Keck/LRIS spectroscopy is complemented by rest-optical Keck/MOSFIRE spectroscopy and Hubble photometry. We find generally larger Ly$\alpha$ equivalent widths ($\rm EW_{Ly\alpha} = 11.2\;$\AA) when compared with higher mass ($\rm \langle log(M^*/M_\odot)\rangle_{median} = 10.3$) composites from KBSS ($\rm EW_{Ly\alpha} = -5\;$\AA). The average low- and high-ionization absorption line EWs ($\rm EW_{LIS}$ and $\rm EW_{HIS}$, respectively) are weaker ($\rm EW_{LIS}$=-1.18 \AA, $\rm EW_{HIS}=$-0.99 \AA) in dwarf galaxies than in higher mass galaxies ($\rm EW_{LIS}$=-2.04 \AA, $\rm EW_{HIS}=$-1.42 \AA). The LIS absorption lines are optically thick and is thus a good tracer of the neutral hydrogen covering fraction. Both higher $\rm EW_{Ly\alpha}$ and lower $\rm EW_{LIS}$ measurements imply that the escape fraction of ionizing radiation is larger in lower-mass galaxies at $z\sim2$.

We have identified 189 candidate $z > 1.3$ protoclusters and clusters in the LSST Deep Drilling Fields. This sample will enable the measurement of the metal enrichment and star formation history of clusters during their early assembly period through the direct measurement of the rate of supernovae identified through the LSST. The protocluster sample was selected from galaxy overdensities in a $Spitzer$/IRAC colour-selected sample using criteria that were optimised for protocluster purity using a realistic lightcone. Our tests reveal that $60-80\%$ of the identified candidates are likely to be genuine protoclusters or clusters, which is corroborated by a $\sim4\sigma$ stacked X-ray signal from these structures. We provide photometric redshift estimates for 47 candidates which exhibit strong peaks in the photo-$z$ distribution of their candidate members. However, the lack of a photo-$z$ peak does not mean a candidate is not genuine, since we find a stacked X-ray signal of similar significance from both the candidates that exhibit photo-$z$ peaks and those that do not. Tests on the lightcone reveal that our pursuit of a pure sample of protoclusters results in that sample being highly incomplete ($\sim4\%$) and heavily biased towards larger, richer, more massive, and more centrally concentrated protoclusters than the total protocluster population. Most ($\sim75\%$) of the selected protoclusters are likely to have a maximum collapsed halo mass of between $10^{13}-10^{14}$ M$_{\odot}$, with only $\sim25\%$ likely to be collapsed clusters above $10^{14}$ M$_{\odot}$. However, the aforementioned bias ensures our sample is $\sim50\%$ complete for structures that have already collapsed into clusters more massive than $10^{14}$ M$_{\odot}$.

Modeling of the stars in the red clump (RC) is challenging, due to the uncertainties associated with the physical processes happening in their core and during the helium flash. By probing the internal stellar structure, asteroseismology allows us to constrain the core properties of RC stars and eventually improve our understanding of this evolutionary phase. In this work, we aim to quantify the impact on the seismic properties of the RC stars of the two main core modeling uncertainties: core boundary mixing, and helium burning nuclear reaction rates. Using the MESA stellar evolution code, we compute models with different core boundary mixing as well as different $3\alpha$ and carbon alpha (C-alpha) nuclear reaction rates. We investigate the impact of those parameters on the period spacing $\Delta \Pi$, which is a probe of the region around the core. We find that different core boundary mixing schemes yield significantly different period spacings, with differences of 30s between the maximal $\Delta \Pi$ value computed with semiconvection and maximal overshoot. We show that increasing the rate of the C-alpha reaction lengthens the core helium burning phase, which extends the range of $\Delta \Pi$ covered by the models during their evolution. This results in a difference of 10s between the models computed with a nominal rate and a rate multiplied by 2, which is larger than the observational uncertainties. The effect of changing the 3$\alpha$ reaction rate is comparatively small. In conclusion, the core boundary mixing is the main source of uncertainty regarding the seismic modeling of RC stars. Moreover, the effect of the C-alpha reaction rate is non-negligible, even if difficult to disentangle from the effect of the mixing. Such degeneracy could be raised in the future, thanks to new seismic data from the PLATO mission and theoretical constraints from numerical simulations.

Although instruments for measuring the radial velocities (RVs) of stars now routinely reach sub-meter per second accuracy, the detection of low-mass planets is still very challenging. The rotational modulation and evolution of spots and/or faculae can induce variations in the RVs at the level of a few m/s in Sun-like stars. To overcome this, a multi-dimensional Gaussian Process framework has been developed to model the stellar activity signal using spectroscopic activity indicators together with the RVs. A recently published computationally efficient implementation of this framework, S+LEAF 2, enables the rapid analysis of large samples of targets with sizeable data sets. In this work, we apply this framework to HARPS observations of 268 well-observed targets with precisely determined stellar parameters. Our long-term goal is to quantify the effectiveness of this framework to model and mitigate activity signals for stars of different spectral types and activity levels. In this first paper in the series, we initially focus on the activity indicators (S-index and Bisector Inverse Slope), and use them to a) measure rotation periods for 49 slow rotators in our sample, b) explore the impact of these results on the spin-down of middle-aged late F, G & K stars, and c) explore indirectly how the spot to facular ratio varies across our sample. Our results should provide valuable clues for planning future RV planet surveys such as the Terra Hunting Experiment or the PLATO ground-based follow-up observations program, and help fine-tune current stellar structure and evolution models.

The ''twin birth'' of a positron and an electron by a photon in the presence of a nucleus, known as Bethe-Heitler pair production, is a key process in astroparticle physics. The Bethe-Heitler process offers a way of channeling energy stored in a population of relativistic protons (or nuclei) to relativistic pairs with extended distributions. Contrary to accelerated leptons, whose maximum energy is limited by radiative losses, the maximal energy of pairs is determined by the kinematics of the process and can be as high as the parent proton energy. We take a closer look at the features of the injected pair distribution, and provide a novel empirical function that describes the spectrum of pairs produced by interactions of single-energy protons with single-energy photons. The function is the kernel of the Bethe-Heitler pair production spectrum that replaces a double numerical integration involving the complex differential cross section of the process, and can be easily implemented in numerical codes. We further examine under which conditions Bethe-Heitler pairs produced in blazar jets can emit $\gamma$-ray photons via synchrotron radiation, thus providing an alternative to the inverse Compton scattering process for high-energy emission in jetted active galactic nuclei. After taking into consideration the broadband spectral characteristics of the source, the jet energetics, and the properties of radiation fields present in the blazar environment, we conclude that $\gamma$-rays in low- and intermediate-peaked blazars may arise from Bethe-Heitler pairs in regions of the jet with typical transverse size $\sim 10^{15}$ cm and co-moving magnetic field $5-500$ G.

Motivated by the recent surge in interest concerning white dwarf (WD) planets, this work presents the first numerical exploration of WD-driven atmospheric escape, whereby the high-energy radiation from a hot/young WD can trigger the outflow of the hydrogen-helium envelope for close-in planets. As a pilot investigation, we focus on two specific cases: a gas giant and a sub-Neptune-sized planet, both orbiting a rapidly cooling WD with mass $M_\ast$ = 0.6 \msun\ and separation $a$ = 0.02 AU. In both cases, the ensuing mass outflow rates exceed $10^{14}$ g sec$^{-1}$ for WD temperatures greater than $T_{\rm WD} \simeq$ 50,000 K. At $T_{\rm WD} \simeq$ 18,000 K [/22,000 K], the sub-Neptune [/gas giant] mass outflow rate approaches $10^{12}$ g sec$^{-1}$, i.e., comparable to the strongest outflows expected from close-in planets around late main-sequence stars. Whereas the gas giant remains virtually unaffected from an evolutionary standpoint, atmospheric escape may have sizable effects for the sub-Neptune, depending on its dynamical history, e.g., assuming that the hydrogen-helium envelope makes up 1 [/4] per cent of the planet mass, the entire envelope would be evaporated away so long as the planet reaches 0.02 AU within the first 230 [/130] Myr of the WD formation. We discuss how these results can be generalized to eccentric orbits with effective semi-major axis $a'=a/(1-e^2)^{1/4}$, which receive the same orbit-averaged irradiation. Extended to a much broader parameter space, this approach can be exploited to model the expected demographics of WD planets as a function of their initial mass, composition and migration history, as well as their potential for habitability.

We identified Active Galactic Nuclei (AGN) candidates as counterparts to unidentified gamma-ray sources (UGS) from the Fermi-LAT Fourth Source Catalogue at lower Galactic latitudes. Our methodology is based on the use of near- and mid-infrared photometric data from the VISTA Variables in the V\'ia L\'actea (VVV) and Wide-field Infrared Survey Explorer (WISE) surveys. The AGN candidates associated with the UGS occupy very different regions from the stars and extragalactic sources in the colour space defined by the VVV and WISE infrared colours. We found 27 near-infrared AGN candidates possibly associated with 14 Fermi-LAT sources using the VVV survey. We also found 2 blazar candidates in the regions of 2 Fermi-LAT sources using WISE data. There is no match between VVV and WISE candidates. We have also examined the K$_\mathrm{s}$ light curves of the VVV candidates and applied the fractional variability amplitude ($\mathrm{\sigma_{rms}}$) and the slope of variation in the K$_\mathrm{s}$ passband to characterise the near-infrared variability. This analysis shows that more than 85% of the candidates have slopes in the K$_\mathrm{s}$ passband $ > 10^{-4}$ mag/day and present $\mathrm{\sigma_{rms}}$ values consistent with a moderate variability. This is in good agreement with typical results seen from type-1 AGN. The combination of YJHK$_\mathrm{s}$ colours and K$_\mathrm{s}$ variability criteria was useful for AGN selection, including its use in identifying counterparts to Fermi $\gamma$-ray sources.

Recent asteroseismic studies have revealed that the convective core of $\gamma$ Doradus stars rotates faster than their radiative interior. We study the development of differential rotation near the convective core to test angular momentum transport processes that are typically adopted in stellar evolution models. Models that only include the advection of angular momentum by meridional circulation and shear instabilities cannot reproduce current rotational constraints, irrespective of the initial conditions. The latest formulation of internal magnetic fields based on the Tayler instability is indeed able to reproduce the internal rotation rate of post-main sequence stars, however, it appears too efficient during the main sequence and has thus been disfavoured. A less efficient version of the same transport process can simultaneously reproduce the rotation rate of the convective core, the rotation rate in radiative regions as probed by gravity-modes, and the surface rotational velocities of $\gamma$ Doradus stars. Our work suggests that there are additional physical processes apart from internal magnetic fields at work in the stellar interiors of post-main sequence stars.

The nova M31N 2023-11f (2023yoa) has been recently identified as the second eruption of a previously recognized nova, M31N 2013-10c, establishing the latter object as the 21st recurrent nova system thus far identified in M31. Here we present well sampled $R$-band lightcurves of both the 2013 and 2023 eruptions of this system. The photometric evolution of each eruption was quite similar as expected for the same progenitor system. The 2013 and 2023 eruptions each reached peak magnitudes just brighter than $R\sim16$, with fits to the declining branches of the eruptions yielding times to decline by two magnitudes of $t_2(R)=5.5\pm1.7$ and $t_2(R)=3.4\pm1.5$ days, respectively. M31N 2013-10c has an absolute magnitude at peak, $M_R=-8.8\pm0.2$, making it the most luminous known recurrent nova in M31.

SPectra Analysis and Retrievable Catalog Lab (SPARCL) at NOIRLab's Astro Data Lab was created to efficiently serve large optical and infrared spectroscopic datasets. It consists of services, tools, example workflows and currently contains spectra for over 7.5 million stars, galaxies and quasars from the Sloan Digital Sky Survey (SDSS) and the Dark Energy Spectroscopic Instrument (DESI) survey. We aim to eventually support the broad range of spectroscopic datasets that will be hosted at NOIRLab and beyond. Major elements of SPARCL include capabilities to discover and query for spectra based on parameters of interest, a fast web service that delivers desired spectra either individually or in bulk as well as documentation and example Jupyter Notebooks to empower users in their research. More information is available on the SPARCL website (https://astrosparcl.datalab.noirlab.edu).

We present a new model to calculate reflection spectra of thin accretion disks in Kerr spacetimes. Our model includes the effect of returning radiation, which is the radiation that is emitted by the disk and returns to the disk because of the strong light bending near a black hole. The major improvement with respect to the existing models is that it calculates the reflection spectrum at every point on the disk by using the actual spectrum of the incident radiation. Assuming a lamppost coronal geometry, we simulate simultaneous observations of NICER and NuSTAR of bright Galactic black holes and we fit the simulated data with the latest version of RELXILL (modified to read the table of REFLIONX, which is the non-relativistic reflection model used in our calculations). We find that RELXILL with returning radiation cannot fit well the simulated data when the black hole spin parameter is very high and the coronal height and disk's ionization parameter are low, and some parameters can be significantly overestimated or underestimated. We can find better fits and recover the correct input parameters as the value of the black hole spin parameter decreases and the values of the coronal height and of the disk's ionization parameter increase.

It has been claimed that the variability of field quasars resembles gravitational lensing by a large cosmological population of free-floating planets with mass of about 10 Earths. But Galactic photometric monitoring experiments, on the other hand, exclude a large population of such planetary-mass gravitational lenses. These apparently contradictory pieces of evidence can be reconciled if the objects under consideration have a mean column-density that lies between the critical column-densities for gravitational lensing in these two contexts. Dark matter in that form is known to be weakly collisional, so that a core develops in galaxy halo density profiles, and a preferred model has already been established. Here we consider what such a model implies for Q2237+0305, which is the best-studied example of a quasar that is strongly lensed by an intervening galaxy. We construct microlensing magnification maps appropriate to the four macro-images of the quasar -- all of which are seen through the bulge of the galaxy. Each of these maps exhibits a caustic network arising from the stars, plus many small, isolated caustics arising from the free-floating "planets" in the lens galaxy. The "planets" have little influence on the magnification histograms but a large effect on the statistics of the magnification gradients. We compare our predictions to the published OGLE photometry of Q2237+0305 and find that these data are consistent with the presence of the hypothetical "planets". However, the evidence is relatively weak because the OGLE dataset is not well suited to testing our predictions and requires low-pass filtering for this application. New data from a large, space-based telescope are desirable to address this issue.

Solar flare prediction studies have been recently conducted with the use of Space-Weather MDI (Michelson Doppler Imager onboard Solar and Heliospheric Observatory) Active Region Patches (SMARP) and Space-Weather HMI (Helioseismic and Magnetic Imager onboard Solar Dynamics Observatory) Active Region Patches (SHARP), which are two currently available data products containing magnetic field characteristics of solar active regions. The present work is an effort to combine them into one data product, and perform some initial statistical analyses in order to further expand their application in space weather forecasting. The combined data are derived by filtering, rescaling, and merging the SMARP with SHARP parameters, which can then be spatially reduced to create uniform multivariate time series. The resulting combined MDI-HMI dataset currently spans the period between April 4, 1996, and December 13, 2022, and may be extended to a more recent date. This provides an opportunity to correlate and compare it with other space weather time series, such as the daily solar flare index or the statistical properties of the soft X-ray flux measured by the Geostationary Operational Environmental Satellites (GOES). Time-lagged cross-correlation indicates that a relationship may exist, where some magnetic field properties of active regions lead the flare index in time. Applying the rolling window technique makes it possible to see how this leader-follower dynamic varies with time. Preliminary results indicate that areas of high correlation generally correspond to increased flare activity during the peak solar cycle.

The issue of radiation mechanisms had triggered in 1950-60s the first applications of plasma physics to understand the nature of radio galaxies. This interplay has steadily intensified during the past five decades, due to the premise of in-situ acceleration of relativistic electrons occurring in the lobes of radio galaxies. This article briefly traces the chain of these remarkable developments, largely from an observational perspective. We recount several observational and theoretical milestones established along the way and the lessons drawn from them. We also present a new observational clue about in-situ acceleration of the relativistic particles radiating in the lobes of radio galaxies, gleaned by us from the very recently published sensitive radio observations of a tailed radio source in the galaxy cluster Abell 1033.

The extreme temperatures and densities of many astrophysical environments tends to destabilize nuclear isomers by inducing transitions to higher energy states. Those states may then cascade to ground. However, not all environments destabilize all isomers. Nuclear isomers which retain their metastable character in pertinent astrophysical environments are known as astrophysically metastable nuclear isomers, or ``astromers''. Astromers can influence nucleosynthesis, altering abundances or even creating new pathways that would otherwise be inaccessible. Astromers may also release energy faster or slower relative to their associated ground state, acting as heating accelerants or batteries, respectively. In stable isotopes, they may even simply remain populated after a cataclysmic event and emit observable x- or $\gamma$-rays. The variety of behaviors of these nuclear species and the effects they can have merit careful consideration in nearly every possible astrophysical environment. Here we provide a brief overview of astromers past and present, and we outline future work that will help to illuminate their role in the cosmos.

Superfluid vortices pinned to nuclear lattice sites or magnetic flux tubes in a neutron star evolve abruptly through a sequence of metastable spatial configurations, punctuated by unpinning avalanches associated with rotational glitches, as the stellar crust spins down electromagnetically. The metastable configurations are approximately but not exactly axisymmetric, causing the emission of persistent, quasimonochromatic, current quadrupole gravitational radiation. The characteristic gravitational wave strain $h_0$ as a function of the spin frequency $f$ and distance $D$ from the Earth is bounded above by $h_0 = 1.2\substack{+1.3 \\ -0.9} \times 10^{-32} (f/30\;{\rm Hz})^{2.5} (D/1\;{\rm kpc})^{-1}$, corresponding to a Poissonian spatial configuration (equal probability per unit area, i.e. zero inter-vortex repulsion), and bounded below by $h_0 = 1.8\substack{+2.0 \\ -1.5} \times 10^{-50} (f/30\;{\rm Hz})^{1.5} (D/1\;{\rm kpc})^{-1}$, corresponding to a regular array (periodic separation, i.e.\ maximum inter-vortex repulsion). N-body point vortex simulations predict an intermediate scaling, $h_0 = 7.3\substack{+7.9 \\ -5.4} \times 10^{-42} (f/30\;{\rm Hz})^{1.9} (D/1\;{\rm kpc})^{-1}$, which reflects a balance between the randomizing but spatially correlated action of superfluid vortex avalanches and the regularizing action of inter-vortex repulsion. The scaling is calibrated by conducting simulations with ${N_{\rm v}} \leq 5\times10^3$ vortices and extrapolated to the astrophysical regime ${N_{\rm v}} \sim 10^{17} (f/30\;{\rm Hz})$. The scaling is provisional, pending future computational advances to raise ${N_{\rm v}}$ and include three-dimensional effects such as vortex tension and turbulence.

We present a novel, deep-learning based method -- dubbed Galactic-Seismology Substructures and Streams Hunter, or GS$^{3}$ Hunter for short, to search for substructures and streams in stellar kinematics data. GS$^{3}$ Hunter relies on a combined application of Siamese Neural Networks to transform the phase space information and the K-means algorithm for the clustering. As a validation test, we apply GS$^{3}$ Hunter to a subset of the Feedback in Realistic Environments (FIRE) cosmological simulations. The stellar streams and substructures thus identified are in good agreement with corresponding results reported earlier by the FIRE team. In the same vein, we apply our method to a subset of local halo stars from the Gaia Early Data Release 3 and GALAH DR3 datasets, and recover several, previously known dynamical groups, such as Thamnos 1+2, Hot Thick Disk, ED-1, L-RL3, Helmi 1+2, and Gaia-Sausage-Enceladus, Sequoia, VRM, Cronus, Nereus. Finally, we apply our method without fine-tuning to a subset of K-giant stars located in the inner halo region, obtained from the LAMOST Data Release 5 (DR5) dataset. We recover three, previously known structures (Sagittarius, Hercules-Aquila Cloud, and the Virgo Overdensity), but we also discover a number of new substructures. We anticipate that GS$^{3}$ Hunter will become a useful tool for the community dedicated to the search of stellar streams and structures in the Milky Way (MW) and the Local group, thus helping advance our understanding of the stellar inner and outer halos, and of the assembly and tidal stripping history in and around the MW.

Multispectral studies of nearby, forming stars provide insights into all classes of accreting systems. Objects which have magnetic fields, spin, and accrete produce jets and collimated outflows. Jets are seen in systems ranging from brown dwarf stars to supermassive black holes. Outflow speeds are typically a few times the escape speed from the launch region - 100s of \kms\ for young stars to nearly the speed of light for black-holes. Because many young stellar objects (YSOs) are nearby, we can see outflow evolution and measure proper motions on times scales of years. Because the shocks in YSO outflows emit in atoms, ions, and molecules in addition to the continuum, many physical properties such as temperatures, densities, and velocities can be measured. Momenta and kinetic energies can be computed. YSO outflows are a major source of feedback in the self-regulation of star formation. The lessons learned can be applied to much more distant and energetic cosmic sources such as AGN and galactic nuclear super winds - systems in which evolution occurs on time-scales of hundreds to millions of years. Some dense star-forming regions produce powerful explosions. The nearest massive star-forming region, Orion OMC1, powered a $\sim 10^{48}$ erg explosion about 550 years ago (that is when the light from the event would have reached the Solar System). The OMC1 explosion was likely powered by an N-body interaction which resulted in the formation of a compact, AU-scale binary or resulted in a protostellar merger. The binary or merger remnant, the $\sim$15 \Msol\ object known as radio source I (Src I) was ejected from the core with a speed of $\sim$10 \kms\ along with two other stars. The $\sim$10~\Msol\ BN object was ejected with $\sim$30~\kms\ and a $\sim$3~\Msol\ star was ejected with $\sim$55~\kms .

We search for high-velocity stars in the inner region of the Galactic bulge using a selected sample of red clump stars. Some of those stars might be considered hypervelocity stars (HVSs). Even though the HVSs ejection relies on an interaction with the supermassive black hole (SMBH) at the centre of the Galaxy, there are no confirmed detections of HVSs in the inner region of our Galaxy. With the detection of HVSs, ejection mechanism models can be constrained by exploring the stellar dynamics in the Galactic centre through a recent stellar interaction with the SMBH. Based on a previously developed methodology by our group, we searched with a sample of preliminary data from version 2 of the Vista Variables in the Via Lactea (VVV) Infrared Astrometric Catalogue (VIRAC2) and Gaia DR3 data, including accurate optical and NIR proper motions. This search resulted in a sample of 46 stars with transverse velocities larger than the local escape velocity within the Galactic bulge, of which 4 are prime candidate HVSs with high-proper motions consistent with being ejections from the Galactic centre. Adding to that, we studied a sample of reddened stars without a Gaia DR3 counterpart and found 481 stars with transverse velocities larger than the local escape velocity, from which 65 stars have proper motions pointing out of the Galactic centre and are candidate HVSs. In total, we found 69 candidate HVSs pointing away from the Galactic centre with transverse velocities larger than the local escape velocity.

We present low-frequency radio observations of the Galactic symbiotic recurrent nova RS Ophiuchi during its 2021 outburst. The observations were carried out with the upgraded Giant Metrewave Radio Telescope (uGMRT) spanning a frequency range of 0.15$-$1.4 GHz during 23$-$287 days post the outburst. The average value of the optically thin spectral index is $\alpha \sim$ $-$0.4 ($F_{\nu} \propto \nu^\alpha$), indicating a non-thermal origin of the radio emission at the observed frequencies. The radio light curves are best represented by shock-driven synchrotron emission, initially absorbed by a clumpy ionized circumbinary medium. We estimate the mass-loss rate of the red giant companion star to be $\dot{M} \sim$ 7.5 $\times$ 10$^{-8}$ $M_{\odot}$ yr$^{-1}$ for an assumed stellar wind velocity of 20 km/s. The 0.15--1.4 GHz radio light curves of the 2021 outburst are systematically brighter than those of the 2006 outburst. Considering similar shock properties between the two outbursts, this is indicative of a relatively higher particle number density in the synchrotron emitting plasma in the current outburst.

The decay of sunspot plays a key role in magnetic flux transportation in solar active regions (ARs). To better understand the physical mechanism of the entire decay process of a sunspot, an {\alpha}-configuration sunspot in AR NOAA 12411 was studied. Based on the continuum intensity images and vector magnetic field data with stray light correction from Solar Dynamics Observatory/Helioseismic and Magnetic Imager, the area, vector magnetic field and magnetic flux in the umbra and penumbra are calculated with time, respectively. Our main results are as follows: (1) The decay curves of the sunspot area in its umbra, penumbra, and whole sunspot take the appearance of Gaussian profiles. The area decay rates of the umbra, penumbra and whole sunspot are -1.56 MSH/day, -12.61 MSH/day and -14.04 MSH/day, respectively; (2) With the decay of the sunspot, the total magnetic field strength and the vertical component of the penumbra increase, and the magnetic field of the penumbra becomes more vertical. Meanwhile, the total magnetic field strength and vertical magnetic field strength for the umbra decrease, and the inclination angle changes slightly with an average value of about 20{\deg}; (3) The magnetic flux decay curves of the sunspot in its umbra, penumbra, and whole sunspot exhibit quadratic patterns, their magnetic flux decay rates of the umbra, penumbra and whole sunspot are -9.84 * 10^19 Mx/day, -1.59 * 10^20 Mx/day and -2.60 * 10^20 Mx/day , respectively. The observation suggests that the penumbra may be transformed into the umbra, resulting in the increase of the average vertical magnetic field strength and the reduction of the inclination angle in the penumbra during the decay of the sunspot.

There is compelling evidence that active galactic nuclei (AGNs) in high-density regions have undergone a different evolution than their counterparts in the field, indicating that they are strongly affected by their environment. To investigate the various factors that may affect the prevalence of AGNs in cluster galaxies, we selected a sample of 19 thoroughly studied X-ray-selected galaxy clusters from the LoCuSS survey. All these clusters are considered massive, with $M_{500}\gtrsim 2\times10^{14} M_{solar}$, and span a narrow redshift range between $z\sim$0.16 and 0.28. We divided the cluster surroundings into two concentric annuli with a width of $R_{500}$ radius. We further divided the cluster sample based on the presence of infalling X-ray-detected groups, cluster mass, or dynamical state. We found that the X-ray AGN fraction in the outskirts is consistent with the field, but it is significantly lower in cluster centres, in agreement with previous results for massive clusters. We show that these results do not depend on cluster mass. Furthermore, we did not find any evidence of a spatial correlation between infalling groups and AGNs. Nevertheless, a significant excess of X-ray AGNs is found in the outskirts of relaxed clusters at the 2$\sigma$ confidence level, compared both to non-relaxed clusters and to the field. Our results suggest that the mechanisms that trigger AGN activity may vary between cluster centres and the outskirts. Ram pressure can efficiently remove the gas from infalling galaxies, thereby triggering AGN activity in some cases. However, the reduced availability of gas globally diminishes the fraction of AGNs in cluster centers. The surplus of X-ray AGNs identified in the outskirts of relaxed clusters may be attributed to an increased frequency of galaxy mergers, a notion that is further supported by the disturbed morphology observed in several galaxies.

Self--interactions of dark matter particles impact the distribution of dark matter in halos. The exact nature of the self--interactions can lead to either expansion or collapse of the core within the halo lifetime, leaving distinctive signatures in the dark matter distributions not only at the halo center but throughout the virial region. Optical galaxy surveys, which precisely measure the weak lensing of background galaxies by massive foreground clusters, allow us to directly measure the matter distribution within clusters and probe subtle effects of self--interacting dark matter (SIDM) throughout the halo's full radial range. We compare the weak--lensing measurements reported by Shin et al. 2021, which use lens clusters identified by the Atacama Cosmology Telescope Survey and source galaxies from the Dark Energy Survey, with predictions from SIDM models having either elastic or dissipative self--interactions. To model the weak--lensing observables, we use cosmological N-body simulations for elastic self--interactions and semi-analytical fluid simulations for dissipative self--interactions. We find that current weak--lensing measurements already constrain the isotropic and elastic SIDM to a cross-section per mass of $\sigma/m<1~{\rm cm^2/g}$ at a $95\%$ confidence level. The same measurements also impose novel constraints on the energy loss per unit mass for dissipative SIDM. Upcoming surveys are anticipated to enhance the signal-to-noise of weak--lensing observables significantly making them effective tools for investigating the nature of dark matter, including self--interactions, through weak lensing.

The inner regions of protoplanetary disks are the locations where most of planets are thought to form and where processes that influence the global evolution of the disk, such as MHD-winds and photoevaporation, originate. Transition disks (TDs) with large inner dust cavities are the ideal targets to study the inner tens of au of disks, as the central emission can be fully disentangled from the outer disk emission. We present a homogeneous multi-wavelength analysis of the continuum emission in a sample of 11 TDs. We investigate the nature of the central emission close to the star, distinguishing between thermal dust and free-free emission. Spatially resolved measurements of continuum emission from archival ALMA data are combined with literature cm-wave observations to study the spectral indices of the inner and outer disks separately. While the emission from the outer disks is consistent with thermal dust emission, 10/11 of the spectral indices estimated for the central emission close to the star suggest that this emission is free-free emission, likely associated with an ionized jet or a disk wind. No correlation between the free-free luminosity and the accretion luminosity or the X-ray luminosity is found, arguing against the photoevaporative wind origin. A sub-linear correlation between the ionized mass loss rate and the accretion rate onto the star is observed, suggesting an origin in an ionized jet. The relative lack of mm-dust grains in the majority of inner disks in transition disks suggests that either such dust grains have drifted quickly towards the central star, grain growth is less efficient in the inner disk, or grains grow rapidly to planetesimal sizes in the inner disk. The observed correlation between the ionized mass loss rate and the accretion rate suggests the outflow is strictly connected with the stellar accretion and that accretion in these disks is driven by a jet.

We present photometric, spectroscopic and polarimetric observations of the intermediate-luminosity Type IIP supernova (SN) 2021gmj from 1 to 386 days after the explosion. The peak absolute V-band magnitude of SN 2021gmj is -15.5 mag, which is fainter than that of normal Type IIP SNe. The spectral evolution of SN 2021gmj resembles that of other sub-luminous supernovae: the optical spectra show narrow P-Cygni profiles, indicating a low expansion velocity. We estimate the progenitor mass to be about 12 Msun from the nebular spectrum and the 56Ni mass to be about 0.02 Msun from the bolometric light curve. We also derive the explosion energy to be about 3 x 10^{50} erg by comparing numerical light curve models with the observed light curves. Polarization in the plateau phase is not very large, suggesting nearly spherical outer envelope. The early photometric observations capture the rapid rise of the light curve, which is likely due to the interaction with a circumstellar material (CSM). The broad emission feature formed by highly-ionized lines on top of a blue continuum in the earliest spectrum gives further indication of the CSM at the vicinity of the progenitor. Our work suggests that a relatively low-mass progenitor of an intermediate-luminosity Type IIP SN can also experience an enhanced mass loss just before the explosion, as suggested for normal Type IIP SNe.

HESS J1832-093 is a member of the rare class of gamma-ray binaries, as recently confirmed by the detection of orbitally modulated X-ray and gamma-ray emission with a period of ~86 d. The spectral type of the massive companion star has been difficult to retrieve as there is no optical counterpart, but the system is coincident with a near-infrared source. Previous results have shown that the infrared counterpart is consistent with an O or B type star, but a clear classification is still lacking. We observed the counterpart twice, in 2019 and 2021, with the X-Shooter spectrograph operating on the VLT. The obtained spectra classify the counterpart as an O6 V type star. We estimate a distance to the source of $6.7 \pm 0.5$ kpc, although this estimate can be severely affected by the high extinction towards the source. This new O6 V classification for the companion star in HESS J1832-093 provides further support to an apparent grouping around a given spectral type for all discovered gamma-ray binaries that contain an O-type star. This may be due to the interplay between the initial mass function and the wind-momentum-luminosity relation.

We present a sample of 950 edge-on spiral galaxies found with the use of an artificial neural network in the Hubble Space Telescope COSMOS field. This is currently the largest sample of distant edge-on galaxies. For all galaxies we analyzed the 2D brightness distributions in the F814W filter and measured the radial and vertical exponential scales ($h$ and $h_z$ correspondingly) of the brightness distribution. By comparing the characteristics of distant galaxies with those of nearby objects, we conclude that thin stellar discs with $h/h_z \geq 10$ at $z \approx 0.5$ should be rarer than today. Both exponential scales of the stellar disc show evidence of luminosity-dependent evolution: in faint galaxies the $h$ and $h_z$ values do not change with $z$, in bright (and massive) spiral galaxies both scales, on average, grow towards our epoch.

The recent detection of significant neutrino flux from the inner Galactic plane by the IceCube detector has provided us valuable insights on the spectrum of cosmic rays in our Galaxy. This flux can be produced either by a population of Galactic point sources or by diffused emission from cosmic ray interactions with the interstellar medium or by a mixture of both. In this work, we compute diffused gamma-ray and neutrino fluxes produced by a population of giant molecular clouds (GMCs) in our Galaxy, assuming different parametrizations of the Galactic diffused cosmic ray distribution. In particular, we take into account two main cases: (I) constant cosmic ray luminosity in our Galaxy, and (II) space-dependent cosmic ray luminosity based on the supernovae distribution in our Galaxy. For Case-I, we found that the neutrino flux from GMCs is a factor of $\sim 10$ below compared to $\pi^0$ and KRA$_\gamma$ best-fitted models of IceCube observations at $10^5$ GeV. Instead, for Case-II the model can explain up to $\sim 90 \%$ of the neutrino flux at that energy. Moreover, for this scenario IceCube detector could be able to detect neutrino events from the Galactic centre regions. We then calculated the gamma-ray and neutrino fluxes from individual GMCs and noticed that several current and future Cherenkov telescopes and neutrino observatories have the right sensitivities to study these objects. In particular, very neutrino-bright region such as Aquila Rift is favourable for detection by the IceCube-Gen2 observatory.

Mildly irradiated mini-Neptunes have densities potentially consistent with them hosting substantial liquid water oceans (`Hycean' planets). The presence of CO2 and simultaneous absence of ammonia (NH3) in their atmospheres has been proposed as a fingerprint of such worlds. JWST observations of K2-18b, the archetypal Hycean, have found the presence of CO2 and the depletion of NH3 to <100 ppm; hence, it has been inferred that this planet may host liquid water oceans. In contrast, climate modelling suggests that many of these mini-Neptunes, including K2-18b, may likely be too hot to host liquid water. We propose a solution to this discrepancy between observation and climate modelling by investigating the effect of a magma ocean on the atmospheric chemistry of mini-Neptunes. We demonstrate that atmospheric NH3 depletion is a natural consequence of the high solubility of nitrogen species in magma at reducing conditions; precisely the conditions prevailing where a thick hydrogen envelope is in communication with a molten planetary surface. The magma ocean model reproduces the present JWST spectrum of K2-18b to < 3 sigma, suggesting this is as credible an explanation for current observations as the planet hosting a liquid water ocean. Spectral areas that could be used to rule out the magma ocean model include the >4um region, where CO2 and CO features dominate: Magma ocean models suggest a systematically lower CO2/CO ratio than estimated from free chemistry retrieval, indicating that deeper observations of this spectral region may be able to distinguish between oceans of liquid water and magma on mini-Neptunes.

The efficiency of radio emission is an important unknown parameter of early galaxies at cosmic dawn, as models with high efficiency have been shown to modify the cosmological 21-cm signal substantially, deepening the absorption trough and boosting the 21-cm power spectrum. Such models have been previously directly constrained by the overall extragalactic radio background as observed by ARCADE-2 and LWA-1. In this work, we constrain the clustering of high redshift radio sources by utilizing the observed upper limits on arcminute-scale anisotropy from the VLA at 4.9~GHz and ATCA at 8.7~GHz. Using a semi-numerical simulation of a plausible astrophysical model for illustration, we show that the clustering constraints on the radio efficiency are much stronger than those from the overall background intensity, by a factor that varies from 12 at redshift 7 to 30 at redshift 22. As a result, the predicted maximum depth of the global 21-cm signal is lowered by a factor of 5 (to 1700~mK), and the maximum 21-cm power spectrum peak at cosmic dawn is lowered by a factor of 24 (to $2\times 10^5$~mK$^2$). We conclude that the observed clustering is the strongest current direct constraint on such models, but strong early radio emission from galaxies remains viable for producing a strongly enhanced 21-cm signal from cosmic dawn.

We present 206 unpublished optical spectra of 104 type II supernovae obtained by the Xinglong 2.16m telescope and Lijiang 2.4m telescope during the period from 2011 to 2018, spanning the phases from about 1 to 200 days after the SN explosion. The spectral line identifications, evolution of line velocities and pseudo equivalent widths, as well as correlations between some important spectral parameters are presented. Our sample displays a large range in expansion velocities. For instance, the Fe~{\sc ii} $5169$ velocities measured from spectra at $t\sim 50$ days after the explosion vary from ${\rm 2000\ km\ s^{-1}}$ to ${\rm 5500\ km\ s^{-1}}$, with an average value of ${\rm 3872 \pm 949\ km\ s^{-1}}$. Power-law functions can be used to fit the velocity evolution, with the power-law exponent quantifying the velocity decline rate. We found an anticorrelation existing between H$\beta$ velocity at mid-plateau phase and its velocity decay exponent, SNe II with higher velocities tending to have smaller velocity decay rate. Moreover, we noticed that the velocity decay rate inferred from the Balmer lines (i.e., H$\alpha$ and H$\beta$) have moderate correlations with the ratio of absorption to emission for H$\alpha$ (a/e). In our sample, two objects show possibly flash-ionized features at early phases. Besides, we noticed that multiple high-velocity components may exist on the blue side of hydrogen lines of SN 2013ab, possibly suggesting that these features arise from complex line forming region. All our spectra can be found in WISeREP and Zenodo.

Type II supernovae (SNe II), which show abundant hydrogen in their spectra, belong to a class of SNe with diverse observed properties. It is commonly accepted that SNe II are produced by core collapse and explosion of massive stars. However, the large photometric and spectroscopic diversity of SNe II, and the mechanisms responsible for these diversities, have not been thoroughly understood. In this review, we first briefly introduce the optical characteristics and possible progenitors of each subtype of SNe II. We then highlight the role of the Chinese Space Station Telescope in future SN studies. With a deep limiting magnitude, the main survey project could detect SN IIP-like objects as distant as $z\sim 1.2$, and obtain UV-optical follow-up for peculiar transients, especially those long-lived events. With a high resolution and a large field of view, the main survey camera is powerful in linking a nearby SN with its progenitor, while the integral field spectrograph is powerful in revealing the SN environment. All this information has the potential to help enrich our understanding of supernova physics.

We present a study of a sample of 45 spectroscopically confirmed, UV luminous galaxies at $z\sim 6$. They were selected as bright Lyman-break galaxies (LBGs) using deep multi-band optical images in more than 2 deg$^2$ of the sky, and subsequently identified via their strong Ly$\alpha$ emission. The majority of these LBGs span an absolute UV magnitude range from $-22.0$ to $-20.5$ mag with Ly$\alpha$ equivalent width (EW) between $\sim$10 and $\sim$200 \AA, representing the most luminous galaxies at $z\sim 6$ in terms of both UV continuum emission and Ly$\alpha$ line emission. We model the SEDs of 10 LBGs that have deep infrared observations from HST, JWST, and/or Spitzer, and find that they have a wide range of stellar masses and ages. They also have high star-formation rates ranging from a few tens to a few hundreds of Solar mass per year. Five of the LBGs have JWST or HST images and four of them show compact morphology in these images, including one that is roughly consistent with a point source, suggesting that UV luminous galaxies at this redshift are generally compact. The fraction of our photometrically selected LBGs with strong Ly$\alpha$ emission ($\mathrm{EW}>25$ \AA) is about $0.2$, which is consistent with previous results and supports a moderate evolution of the IGM opacity at the end of cosmic reionization. Using deep X-ray images, we do not find evidence of strong AGN activity in these galaxies, but our constraint is loose and we are not able to rule out the possibility of any weak AGN activity.

Using the pair-count implementaion from the Corrfunc package we show that with a low discrepency sequence we can calculate the two-point correlation function more accurately than with random points at no extra computational cost.

Deriving atmospheric parameters of a large sample of stars is of vital importance to understand the formation and evolution of the Milky Way. Photometric surveys, especially those with near-ultraviolet filters, can offer accurate measurements of stellar parameters, with the precision comparable to that from low/medium resolution spectroscopy. In this study, we explore the capability of measuring stellar atmospheric parameters from CSST broad-band photometry (particularly the near-ultraviolet bands), based on synthetic colors derived from model spectra. We find that colors from the optical and near-ultraviolet filter systems adopted by CSST show significant sensitivities to the stellar atmospheric parameters, especially the metallicity. According to our mock data tests, the precision of the photometric metallicity is quite high, with typical values of 0.17 dex and 0.20 dex for dwarf and giant stars, respectively. The precision of the effective temperature estimated from broad-band colors are within 50 K.

The minimum variation timescale (MVT) and spectral lag of hundreds of X-ray bursts (XRBs) from soft gamma-ray repeater (SGR) J1935+2154 were analyzed in detail for the first time in our recent work, which are important probes for studying the physical mechanism and radiation region. In this work, we investigate their differential and cumulative distributions carefully and find that they follow power-law models. Besides, the distributions of fluctuations in both parameters follow the Tsallis $q$-Gaussian distributions and the $q$ values are consistent for different scale intervals. Therefore, these results indicate that both parameters are scale-invariant, which provides new parameters for the study of self-organized criticality systems. Interestingly, we find that the $q$ values for MVT and spectral lag are similar with duration and fluence, respectively.

We present ALMA observations of two dense gas tracers, HCN(1-0) and HCO$^{+}$(1-0), for three galaxies in the green valley and two galaxies on the star-forming main sequence with comparable molecular gas fractions as traced by the CO(1-0) emissions, selected from the ALMaQUEST survey. We investigate whether the deficit of molecular gas star formation efficiency (SFE$_{\rm mol}$) that leads to the low specific star formation rate in these green valley galaxies is due to a lack of dense gas (characterized by the dense gas fraction $f_{\rm dense}$) or the low star formation efficiency of dense gas (SFE$_{\rm dense}$). We find that SFE$_{\rm mol}$ as traced by the CO emissions, when considering both star-forming and retired spaxels together, is tightly correlated with SFE$_{\rm dense}$ and depends only weakly on $f_{\rm dense}$. The specific star formation rate (sSFR) on kpc scales is primarily driven by SFE$_{\rm mol}$ and SFE$_{\rm dense}$, followed by the dependence on $f_{\rm mol}$, and is least correlated with $f_{\rm dense}$ or the dense-to-stellar mass ratio ($R_{\rm dense}$). When compared with other works in the literature, we find that our green valley sample shows lower global SFE$_{\rm mol}$ as well as lower SFE$_{\rm dense}$ while exhibiting similar dense gas fractions when compared to star-forming and starburst galaxies. We conclude that the star formation of the 3 green valley galaxies with a normal abundance of molecular gas is suppressed mainly due to the reduced SFE$_{\rm dense}$ rather than the lack of dense gas.

Aims. We studied the first multi-spacecraft high-energy solar energetic particle (SEP) event of solar cycle 25, which triggered a ground level enhancement (GLE) on 28 October 2021, using data from multiple observers that were widely distributed throughout the heliosphere. Methods. We performed detail modelling of the shock wave and investigated the magnetic connectivity of each observer to the solar surface and examined the shock magnetic connection. We performed 3D SEP propagation simulations to investigate the role of particle transport in the distribution of SEPs to distant magnetically connected observers. Results. Observations and modelling show that a strong shock wave formed promptly in the low corona. At the SEP release time windows, we find a connection with the shock for all the observers. PSP, STA, and Solar Orbiter were connected to strong shock regions with high Mach numbers, whereas the Earth and other observers were connected to lower Mach numbers. The SEP spectral properties near Earth demonstrate two power laws, with a harder (softer) spectrum in the low-energy (high-energy) range. Composition observations from SIS (and near-Earth instruments) show no serious enhancement of flare-accelerated material. Conclusions. A possible scenario consistent with the observations and our analysis indicates that high-energy SEPs at PSP, STA, and Solar Orbiter were dominated by particle acceleration and injection by the shock, whereas high-energy SEPs that reached near-Earth space were associated with a weaker shock; it is likely that efficient transport of particles from a wide injection source contributed to the observed high-energy SEPs. Our study cannot exclude a contribution from a flare-related process; however, composition observations show no evidence of an impulsive composition of suprathermals during the event, suggestive of a non-dominant flare-related process.

We report results of optical identification and multi-wavelength study of a new polar-type magnetic cataclysmic variable (MCV), SRGA J213151.5+491400, discovered by Spectrum Roentgen-Gamma ($SRG$) observatory in the course of the all-sky survey. We present optical data from telescopes in Turkey (RTT-150 and T100 at the T\"UBITAK National Observatory), and in Russia (6-m and 1-m at SAO RAS), together with the X-ray data obtained with $ART-XC$ and $eROSITA$ telescopes aboard $SRG$ and the $NICER$ observatory. We detect SRGA J213151.5+491400 in a high state in 2020 (17.9 mag) that decreases about 3 mag into a low state (21 mag) in 2021. We find only one significant period using optical photometric time series analysis which reveals the white dwarf spin/orbital period to be 0.059710(1) days (85.982 min). The long slit spectroscopy in the high state yields a power law continuum increasing towards the blue with a prominent He II line along with the Balmer line emissions with no cyclotron humps; consistent with MCV nature. Doppler Tomography confirms the polar nature revealing ballistic stream accretion along with magnetic stream during the high state. These characteristics show that the new source is a polar-type MCV. $SRG$ $ART-XC$ detections yield an X-ray flux of (4.0-7.0)$\times$10$^{-12}$ erg cm$^2$ s$^{-1}$ in the high state. $eROSITA$ detects a dominating hot plasma component (kT$_{\rm{max}}$ $>$ 21 keV in the high state) declining to (4.0-6.0)$\times$10$^{-13}$ erg cm$^2$ s$^{-1}$ in 2021 (low state). The $NICER$ data obtained in the low state reveal a two-pole accretor showing a soft X-ray component at (6-7)$\sigma$ significance with a blackbody temperature of 15-18 eV. A soft X-ray component has never been detected for a polar in the low state before.

HESS J1825-137 is one of the most powerful and luminous TeV gamma-ray pulsar wind nebulae (PWNe), making it an excellent laboratory to study particle transportation around pulsars. We present a model of the (diffusive and advective) transport and radiative losses of electrons from the pulsar PSRJ1826-1334 powering HESSJ1825-137 using interstellar medium gas (ISM) data, soft photon fields and a spatially varying magnetic field. We find that for the characteristic age of 21 kyr, PSR J1826-1334 is unable to meet the energy requirements to match the observed X-ray and gamma-ray emission. An older age of 40 kyr, together with an electron conversion efficiency of 0.14 and advective flow of $v = 0.002c$, can reproduce the observed multi-wavelengh emission towards HESS J1825-137. A turbulent ISM with magnetic field of $B = 20\,{\mu}G$ to $60\,{\mu}G$ to the north of HESS J1825-137 (as suggested by ISM observations) is required to prevent significant gamma-ray contamination towards the northern TeV source HESS J1826-130.

The Generalized Uncertainty Principle (GUP) is motivated by the premise that spacetime fluctuations near the Planck scale impose a lower bound on the achievable resolution of distances, leading to a minimum length. Inspired by a semiclassical method that integrates the GUP into the partition function by deforming its phase space, we induce a modification on the thermodynamic quantities of the MIT bag model that we propose serves as an effective semiclassical description of deconfined quark matter in a space with minimal length. We investigate the consequences of this deformation on the zero-temperature limit, revealing a saturation limit for the energy density, pressure and baryon number density and an overall decrease of the thermodynamic quantities which suggests an enhanced stability against gravitational collapse. These findings extend existing research on GUP-deformed Fermi gases. Ultimately, our description introduces effects of quantum gravity in the equations of state for compact stars in a mathematically simple manner, suggesting potential for extension to more complex systems.

We apply methods of particle track reconstruction in High Energy Physics (HEP) to the search for distinct stellar populations in the Milky Way, using the Gaia EDR3 dataset. This was motivated by analogies between the 3D space points in HEP detectors and the positions of stars (which are also points in a coordinate space) and the way collections of space points correspond to particle trajectories in the HEP, while collections of stars from distinct populations (such as stellar streams) can resemble tracks. Track reconstruction consists of multiple steps, the first one being seeding. In this note we describe our implementation and results of the seeding step to the search for distinct stellar populations, and we indicate how the next step will proceed. Our seeding method uses machine learning tools from the FAISS library, such as the kNN nearest neighbor search.

To probe the star formation process, we present an observational investigation of the Pillar IV and an ionized knot HH 216 in the Eagle Nebula (M16). Pillar IV is known to host a Class I protostar that drives a bipolar outflow. The outflow has produced the bow shock, HH 216, which is associated with the red-shifted outflow lobe. The James Webb Space Telescope's near- and mid-infrared images (resolution $\sim$0.07 arcsec - 0.7 arcsec) reveal the protostar as a single, isolated object (below 1000 AU). The outer boundary of Pillar IV is depicted with the 3.3 $\mu$m Polycyclic aromatic hydrocarbon (PAH) emission. HH 216 is traced with the 4.05 $\mu$m Br$\alpha$ and the radio continuum emission, however it is undetected with 4.693 $\mu$m H$_{2}$ emission. HH 216 seems to be associated with both thermal and non-thermal radio emissions. High-resolution images reveal entangled ionized structures (below 3000 AU) of HH 216, which appear to be located toward termination shocks. New knots in 4.693 $\mu$m H$_{2}$ emission are detected, and are mainly found on Pillar IV's northern side. This particular result supports the previously proposed episodic accretion in the powering source of HH 216. One part of the ionized jet (extent $\sim$0.16 pc) is discovered on the southern side of the driving source. Using the $^{12}$CO($J$ = 1-0), $^{12}$CO($J$ = 3-2), and $^{13}$CO($J$ = 1-0) emission, observational signposts of Cloud-Cloud Collision (or interacting clouds) toward Pillar IV are investigated. Overall, our results suggest that the interaction of molecular cloud components around 23 and 26 km s$^{-1}$ might have influenced star formation activity in Pillar IV.

We use high-resolution cosmological simulations to compare the effect of bursty star formation histories on dwarf galaxy structure for two different subgrid supernovae (SNe) feedback models in dwarf galaxies with stellar masses from $5000 <$ M$_*$/M$_\odot$ $< 10^{9}$. Our simulations are run using two distinct supernova feedback models: superbubble and blastwave. We show that both models are capable of producing galaxies that are cored and reproduce observed scaling relations for metallicity, luminosity, mass, and size. We show that continuous bursty star formation and the resulting stellar feedback are able to sustain dark matter cores in the higher dwarf galaxy mass regime, while the majority of ultra-faint and classical dwarfs retain cuspy central dark matter density profiles. We find that both subgrid SN models are able to create bursty star formation histories. We find that effective core formation peaks at M$_*$/M$_\odot$ $\simeq 5 \times 10^{-3}$ for both feedback models. Galaxies simulated with superbubble feedback peak at lower mean burstiness values relative to blastwave feedback, indicating that core formation in the superbubble sample may be less motivated by the burstiness of star formation.

We present a JWST/NIRCam transmission spectrum from $3.9-5.0$ $\mu$m of the recently-validated sub-Earth GJ 341b ($\mathrm{R_P} = 0.92$ $\mathrm{R_{\oplus}}$, $\mathrm{T_{eq}} = 540$ K) orbiting a nearby bright M1 star ($\mathrm{d} = 10.4$ pc, $\mathrm{K_{mag}}=5.6$). We use three independent pipelines to reduce the data from the three JWST visits and perform several tests to check for the significance of an atmosphere. Overall, our analysis does not uncover evidence of an atmosphere. Our null hypothesis tests find that none of our pipelines' transmission spectra can rule out a flat line, although there is weak evidence for a Gaussian feature in two spectra from different pipelines (at 2.3 and $2.9\sigma$). However, the candidate features are seen at different wavelengths (4.3 $\mu$m vs 4.7 $\mu$m), and our retrieval analysis finds that different gas species can explain these features in the two reductions (CO$_2$ at $3.1\sigma$ compared to O$_3$ at $2.9\sigma$), suggesting that they are not real astrophysical signals. Our forward model analysis rules out a low mean molecular weight atmosphere ($< 350\times$ solar metallicity) to at least $3\sigma$, and disfavors CH$_4$-dominated atmospheres at $1-3\sigma$, depending on the reduction. Instead, the forward models find our transmission spectra are consistent with no atmosphere, a hazy atmosphere, or an atmosphere containing a species that does not have prominent molecular bands across the NIRCam/F444W bandpass, such as a water-dominated atmosphere. Our results demonstrate the unequivocal need for two or more transit observations analyzed with multiple reduction pipelines, alongside rigorous statistical tests, to determine the robustness of molecular detections for small exoplanet atmospheres.

We present a systematic study of cosmological parameter bias in weak lensing and large-scale structure analyses for upcoming imaging surveys induced by the interplay of intrinsic alignments (IA) and photometric redshift (photo-z) model mis-specification error. We first examine the degeneracies between the parameters of the Tidal Alignment - Tidal Torquing (TATT) model for IA and of a photo-z model including a mean shift ($\Delta \bar{z}$) and variance ($\sigma_{z}$) for each tomographic bin of lenses and sources, under a variety of underlying true IA behaviors. We identify strong degeneracies between: (1) the redshift scaling of the tidal alignment amplitude and the mean shift and variances of source bins, (2) the redshift scaling of the tidal torquing amplitude and the variance of the lowest-$z$ source bin, and (3) the IA source density weighting and the mean shift and variance of several source bins. We then use this information to guide our exploration of the level of cosmological parameter bias which can be induced given incorrect modelling of IA, photo-z, or both. We find that marginalizing over all the parameters of TATT is generally sufficient to preclude cosmological parameter bias in the scenarios we consider. However, this does not necessarily mean that IA and photo-z parameters are themselves unbiased, nor does it mean that the best-fit model is a good fit to the data. We also find scenarios where the inferred parameters produce $\chi^2_{\rm DOF}$ values indicative of a good fit but cosmological parameter bias is significant, particularly when the IA source density weighting parameter is not marginalized over.

Neural networks are being extensively used for modelling data, especially in the case where no likelihood can be formulated. Although in the case of X-ray spectral fitting, the likelihood is known, we aim to investigate the neural networks ability to recover the model parameters but also their associated uncertainties, and compare its performance with standard X-ray spectral fitting, whether following a frequentist or Bayesian approach. We apply Simulation-Based Inference with Neural Posterior Estimation (SBI-NPE) to X-ray spectra. We train a network with simulated spectra, and then it learns the mapping between the simulated spectra and their parameters and returns the posterior distribution. The model parameters are sampled from a predefined prior distribution. To maximize the efficiency of the training of the neural network, yet limiting the size of the training sample to speed up the inference, we introduce a way to reduce the range of the priors, either through a classifier or a coarse and quick inference of one or multiple observations. SBI-NPE is demonstrated to work equally well as standard X-ray spectral fitting, both in the Gaussian and Poisson regimes, both on simulated and real data, yielding fully consistent results in terms of best fit parameters and posterior distributions. The inference time is comparable to or smaller than the one needed for Bayesian inference. On the other hand, once properly trained, an amortized SBI-NPE network generates the posterior distributions in no time. We show that SBI-NPE is less sensitive to local minima trapping than standard fit statistic minimization techniques. We find that the neural network can be trained equally well on dimension-reduced spectra, via a Principal Component Decomposition, leading to a shorter inference time. Neural posterior estimation thus adds up as a complementary tool for X-ray spectral fitting. (abridged).

The ever-growing sample of observed supernovae enhances our capacity for comprehensive supernova population studies, providing a richer dataset for understanding the diverse characteristics of Type Ia supernovae and possibly that of their progenitors. Here, we present a data-driven analysis of observed Type Ia supernova photometric light curves collected in the Open Supernova Catalog. Where available, we add the environmental information from the host galaxy. We focus on identifying sub-classes of Type Ia supernovae without imposing the pre-defined sub-classes found in the literature to date. To do so, we employ an implicit-rank minimizing autoencoder neural network for developing low-dimensional data representations, providing a compact representation of the supernova light curve diversity. When we analyze light curves alone, we find that one of our resulting latent variables is strongly correlated with redshift, allowing us to approximately ``de-redshift'' the other latent variables describing each event. After doing so, we find that three of our latent variables account for $\sim$95\% of the variance in our sample, and provide a natural separation between 91T and 91bg thermonuclear supernovae. Of note, the 02cx subclass is not unambiguously delineated from the 91bg sample in our results, nor do either the over-luminous 91T or the under-luminous 91bg/02cx samples form a clearly distinct population from the broader sample of ``other'' SN Ia events. We identify the physical characteristics of supernova light curves which best distinguish SNe 91T from SNe 91bg \& 02cx, and discuss prospects for future refinements and applications to other classes of supernovae as well as other transients.

Spectral lightcurves consisting of time series single-pixel spectral measurements of spacecraft are used to infer the spacecraft's attitude and rotation. Two methods are used. One based on numerical optimisation of a regularised least squares cost function, and another based on machine learning with a neural network model. The aim is to work with minimal information, thus no prior is available on the attitude nor on the inertia tensor. The theoretical and practical aspects of this task are investigated, and the methodology is tested on synthetic data. Results are shown based on synthetic data.

Searches for signals at low signal-to-noise ratios frequently involve the Fast Fourier Transform (FFT). For high-throughput searches, we here consider FFT on the homogeneous mesh of Processing Elements (PEs) of a wafer-scale engine (WSE). To minimize memory overhead in the inherently non-local FFT algorithm, we introduce a new synchronous slide operation ({\em Slide}) exploiting the fast interconnect between adjacent PEs. Feasibility of compute-limited performance is demonstrated in linear scaling of Slide execution times with varying array size in preliminary benchmarks on the CS-2 WSE. The proposed implementation appears opportune to accelerate and open the full discovery potential of FFT-based signal processing in multi-messenger astronomy.

Differential equations of the form $\ddot R=-kR^\gamma$, with a positive constant $k$ and real parameter $\gamma$, are fundamental in describing phenomena such as the spherical gravitational collapse ($\gamma=-2$), the implosion of cavitation bubbles ($\gamma=-4$) and the orbital decay in binary black holes ($\gamma=-7$). While explicit elemental solutions exist for select integer values of $\gamma$, more comprehensive solutions encompassing larger subsets of $\gamma$ have been independently developed in hydrostatics (see Lane-Emden equation) and hydrodynamics (see Rayleigh-Plesset equation). This paper introduces a general explicit solution for all real $\gamma$, employing the quantile function of the beta distribution, readily available in most modern programming languages. This solution bridges between distinct fields and reveals insights, such as a critical branch point at $\gamma=-1$, thereby enhancing our understanding of these pervasive differential equations.

This study aims to investigate the strong gravitational lensing effects in $f(T)$ gravity. We present the theoretical analytic expressions for the lensing effects in $f(T)$ gravity, including deflection angle, magnification, and time delay. On this basis, we also take the plasma lensing effect into consideration. We compare the lensing effects between the General Relativity in a vacuum environment and the $f(T)$ gravity in a plasma environment. From a strongly lensed fast radio burst, the results indicate that in a plasma environment, General Relativity and $f(T)$ gravity can generate indistinguishable image positions, but the magnification and time delay on these positions are significantly different, which can be distinguished by current facilities in principle. Therefore, the discrepancies between observational results and theoretical expectations can serve as clues for a modified gravity theory and provide constraints on $f(T)$ gravity.

Gravitational wave observations of binary black hole mergers probe their astrophysical origins via the binary spin, namely the spin magnitudes and directions of each component black hole, together described by six degrees of freedom. However, the emitted signals primarily depend on two effective spin parameters that condense the spin degrees of freedom to those parallel and those perpendicular to the orbital plane. Given this reduction in dimensionality between the physically relevant problem and what is typically measurable, we revisit the question of whether information about the component spin magnitudes and directions can successfully be recovered via gravitational-wave observations, or if we simply extrapolate information about the distributions of effective spin parameters.To this end, we simulate three astrophysical populations with the same underlying effective-spin distribution but different spin magnitude and tilt distributions, on which we conduct full individual-event and population-level parameter estimation. We find that parameterized population models can indeed qualitatively distinguish between populations with different spin magnitude and tilt distributions at current sensitivity. However, it remains challenging to either accurately recover the true distribution or to diagnose biases due to model misspecification. We attribute the former to practical challenges of dealing with high-dimensional posterior distributions, and the latter to the fact that each individual event carries very little information about the full six spin degrees of freedom.

Neutrinos from a Galactic core collapse supernova can undergo elastic scattering with electrons in xenon atoms in liquid xenon based dark matter detectors giving rise to electrons of kinetic energy up to a few MeV. We calculate the scattered electron spectrum and the number of such elastic scattering events expected for a typical Galactic core collapse supernova in a xenon target. Although the expected number of events is small (compared to, for example, inelastic neutrino-nucleus charged current interaction with xenon nuclei, that also gives rise to final state electrons), the distinct spectral shape of the scattered electrons may allow identification of the elastic scattering events. Further, while the process is dominated by neutrinos and antineutrinos of electron flavor, it receives contributions from all the six neutrino species. Identification of the electron scattering events may, therefore, allow an estimation of the relative fractions of the total supernova explosion energy going into electron flavored and non-electron flavored neutrinos.

In this note, I describe a gravitational effect that generically limits the kinetic energy of a single massive elementary particle in the vicinity of a compact object. In the rest frame of a scattering trajectory, tidal accelerations have a quadratic dependence on the specific kinetic energy. As the kinetic energy is increased, the differences in the tidal potential over a Compton wavelength will at some point become large enough to create additional particles. A straightforward calculation reveals that neutrinos scattering off a $10 M_\odot$ black hole within three Schwarzschild radii are roughly limited to about $1~\text{GeV}$, so that an incident neutrino of significantly higher energy passing through such a region decays into a shower of neutrinos with individual energies below the threshold.

This study is devoted to the inference problem of extracting the nuclear matter properties directly from a set of mass-radius observations. We employ Bayesian neural networks (BNNs), which is a probabilistic model capable of estimating the uncertainties associated with its predictions. To simulate different noise levels on the $M(R)$ observations, we create three different sets of mock data. Our results show BNNs as an accurate and reliable tool for predicting the nuclear matter properties whenever the true values are not completely outside the training dataset statistics, i.e., if the model is not heavily dependent on its extrapolating capacities. Using real mass-radius pulsar data, the model predicted, for instance, $L_{\text{sym}}=39.80\pm17.52 $ MeV and $K_{\text{sym}}=-101.67\pm62.86 $ MeV ($2\sigma$ interval). Our study provides a valuable inference framework when new NS data becomes available.

Neutrinos from the cosmos have proven to be ideal for probing the nature of space-time. Previous studies on high-energy events of IceCube suggested that some of these events might be gamma-ray burst neutrinos, with their speeds varying linearly with their energy, implying also the coexistence of subluminal and superluminal propagation. However, a recent reanalysis of the data, incorporating revised directional information, reveals stronger signals that neutrinos are actually being slowed down compared to previous suggestion of neutrino speed variation. Thus, it is worth discussing its implications for the brane/string inspired framework of space-time foam, which has been used to explain previous observations. We revisit effects on neutrino propagation from specific foam models within the framework, indicating that the implied violation of Lorentz invariance could necessarily cause the neutrino to decelerate. We therefore argue that this sort of model is in agreement with the updated phenomenological indication just mentioned. An extended analysis of the revised IceCube data will further test these observations and stringy quantum gravity.

The Bodeker-Moore thermal friction is usually used to determine whether or not a bubble wall can run away. However, the friction on the wall is not necessarily a monotonous function of the wall velocity and could have a maximum before it reaches the Bodeker-Moore limit. In this paper, we compare the maximal hydrodynamic obstruction, i.e., a frictional force in local thermal equilibrium that originates from inhomogeneous temperature distribution across the wall, and the Bodeker-Moore thermal friction. We study the former in a fully analytical way, clarifying its physical origin and providing a simple expression for its corresponding critical phase transition strength above which the driving force cannot be balanced out by the maximal hydrodynamic obstruction. We find that for large parameter space, the maximal hydrodynamic obstruction is larger than the Bodeker-Moore thermal friction, indicating that the conventional criterion for the runaway behavior of the bubble wall must be modified. We also explain how to apply efficiently the modified criterion to particle physics models and discuss possible limitations of the analysis carried out in this paper.

Based on a recently proposed reinterpretation of gravitational wave memory that builds up on the definition of gravitational waves pioneered by Isaacson, we provide a unifying framework to derive both ordinary and null memory from a single well-defined equation at leading order in the asymptotic expansion. This allows us to formulate a memory equation that is valid for any unbound asymptotic energy-flux that preserves local Lorentz invariance. Using Horndeski gravity as a concrete example metric theory with an additional potentially massive scalar degree of freedom in the gravitational sector, the general memory formula is put into practice by presenting the first account of the memory correction sourced by the emission of massive field waves. Throughout the work, physical degrees of freedom are identified by constructing manifestly gauge invariant perturbation variables within an SVT decomposition on top of the asymptotic Minkowski background, which will in particular prove useful in future studies of gravitational wave memory within vector tensor theories.

It is shown that the evolution of an axially and reflection symmetric fluid distribution, satisfying the Tolman condition for thermal equilibrium, is not accompanied by the emission of gravitational radiation. This result, which was conjectured by Bondi many years ago, expresses the irreversibility associated to the emission of gravitational waves. The observational consequences emerging from this result are commented. The resulting models are not only non--dissipative and vorticity free, but also shear--free and geodesic, furthermore all their complexity factors vanish.

The equations of state (EoSs) governing neutron star (NS) matter obtained for both non-relativistic and relativistic mean-field models are systematically confronted with a diverse set of terrestrial data and astrophysical observations within the Bayesian framework. The terrestrial data spans from bulk properties of finite nuclei to the heavy-ion collisions which constrain the symmetric nuclear matter EoS and the symmetry energy up to twice the saturation density ($\rho_0$= 0.16 fm$^{-3}$). The astrophysical observations encompass the NS radius, the tidal deformability, and the lower bound on maximum mass. Three sets of EoSs distributions are generated by gradually updating them with different constraints: (i) we use only the maximum NS mass, (ii) then incorporate additional terrestrial data, (iii) and finally, we include further astrophysical observations. We compare these sets using the Kullback-Liebler divergence. Our results for the Kullback-Liebler divergence highlight the significant constraints imposed on the EoSs by the currently available lower bound of neutron star maximum mass and terrestrial data. The remaining astrophysical observations further refine the EoS within the density range $\sim$ 2-3$\rho_0$.

The symmetry-breaking first-order phase transition between superfluid phases $^3$He-A and $^3$He-B can be triggered extrinsically by ionising radiation or heterogeneous nucleation arising from the details of the sample cell construction. However, the role of potential homogeneous intrinsic nucleation mechanisms remains elusive. Discovering and resolving the intrinsic processes may have cosmological consequences, since an analogous first-order phase transition, and the production of gravitational waves, has been predicted for the very early stages of the expanding Universe in many extensions of the Standard Model of particle physics. Here we introduce a new approach for probing the phase transition in superfluid $^3$He. The setup consists of a novel stepped-height nanofluidic sample container with close to atomically smooth walls. The $^3$He is confined in five tiny nanofabricated volumes and assayed non-invasively by NMR. Tuning of the state of $^3$He by confinement is used to isolate each of these five volumes so that the phase transitions in them can occur independently and free from any obvious sources of heterogeneous nucleation. The small volumes also ensure that the transitions triggered by ionising radiation are strongly suppressed. Here we present the preliminary measurements using this setup, showing both strong supercooling of $^3$He-A and superheating of $^3$He-B, with stochastic processes dominating the phase transitions between the two. The objective is to study the nucleation as a function of temperature and pressure over the full phase diagram, to both better test the proposed extrinsic mechanisms and seek potential parallel intrinsic mechanisms.

We present a simple but powerful Lagrangian method that can be used to study the post-Newtonian evolution of a compact binary system with environment, including a dark matter spike, around it, and obtain the resulting gravitational wave emission. This formalism allows one to incorporate post-Newtonian effects up to any desired known order, as well as any other environmental effect around the binary, as long as their dissipation power or force formulae are known. In particular, in this work, we employ this method to study a black hole-black hole binary system of mass ratio $10^5$, by including post-Newtonian effects of order 1PN and 2.5PN as well as the effect of relativistic dynamical friction. We obtain the modified orbits and the corresponding modified gravitational waveform. Finally, we contrast these modifications against the LISA sensitivity curve in frequency space and show that this observatory can detect the associated signals.

Explore alien solar systems via local star power using interstellar photovoltaics, tailored for the particular target star for maximum power and low mass. We consider tailored organic thin-film photovoltaics. Key for sensing, sending more data back and powering A.I. to send back observational summaries and interesting events and observations. This plus other technology developments are necessary for exploring Alien Solar Systems in the not too distant future.