Papers for Tuesday, Apr 25 2023

AIMS: We present an analysis of an XMM-Newton observation of the dwarf novae V1504 Cyg during the decline from an outburst. Our goal is to search for evidence for an evaporated X-ray corona. Such a corona can be understood as an optically thin geometrically thick disc around a central part of an optically thick geometrically thin disc. METHODS: We study the X-ray spectra using a cooling flow model and the evolution of the amplitude of variability and power density spectra in UV and X-rays. RESULTS: The X-ray (pn) count rate increases from initially around 0.03 cps to 0.17 cps with a harder spectrum and a higher degree of variability. Meanwhile, the OM/UVW1 light curve follows a slow decline with decreasing amplitude of variability. For further study we split the X-ray data into two parts, and analysed them separately. Both parts are described by a cooling flow model, while the first low luminosity part requires an additional power law component suggesting presence of a wind. Spectral fitting revealed a higher temperature during the second brighter part. Timing analysis reveals a potential break frequency at log(f/Hz) = -3.02 during decline towards the quiescence. This detection agrees with optical data from Kepler observations. CONCLUSIONS: The X-ray nature of the break frequency supports the innermost parts of the disc as source of the variability. Moreover, a similar frequency was observed in several other cataclysmic variables and a sandwich model where a geometrically thick corona surrounds the geometrically thin disc is a possible accretion configuration.

The Apache Point Lunar Laser-ranging Operation (APOLLO) has been collecting lunar range measurements for 15 years at millimeter accuracy. The median nightly range uncertainty since 2006 is 1.7 mm. A recently added Absolute Calibration System (ACS), providing an independent assessment of APOLLO system accuracy and the capability to correct lunar range data, revealed a 0.4% systematic error in the calibration of one piece of hardware that has been present for the entire history of APOLLO. Application of ACS-based timing corrections suggests systematic errors are reduced to < 1 mm, such that overall data accuracy and precision are both 1 mm. This paper describes the processing of APOLLO/ACS data that converts photon-by-photon range measurements into the aggregated normal points that are used for physics analyses. Additionally we present methodologies to estimate timing corrections for range data lacking contemporaneous ACS photons, including range data collected prior to installation of the ACS. We also provide access to the full 15-year archive of APOLLO normal points (2006-04-06 to 2020-12-27).

We introduce an analytic surface density profile for dark matter halos that accurately reproduces the structure of simulated halos of mass $M_{\rm vir} = 10^{7-11}\ M_\odot$, making it useful for modeling line-of-sight perturbers in strong gravitational lensing models. The two-parameter function has an analytic deflection potential and is more accurate than the projected Navarro, Frenk & White (NFW) profile commonly adopted at this mass scale for perturbers, especially at the small radii of most relevant for lensing perturbations. Using a characteristic radius, $R_{-1}$, where the log slope of surface density is equal to $-1$, and an associated surface density, $\Sigma_{-1}$, we can represent the expected lensing signal from line-of-sight halos statistically, for an ensemble of halo orientations, using a distribution of {\em projected concentration} parameters, $\mathcal{C}_{\rm vir} := r_{\rm vir}/ R_{-1}$. Though an individual halo can have a projected concentration that varies with orientation with respect to the observer, the range of projected concentrations correlates with the usual three-dimensional halo concentration in a way that enables ease of use.

Attenuation curves in galaxies depend on dust chemical composition, content, and grain size distribution. Such parameters are related to intrinsic galaxy properties such as metallicity, star formation rate, and stellar age. Due to the lack of observational constraints at high redshift, dust empirical curves measured in the local Universe (e.g. Calzetti and SMC curves) have been employed to describe the dust attenuation at early epochs. We exploit the high sensitivity and spectral resolution of the JWST to constrain the dust attenuation curves in high-z galaxies. Our goals are to check whether dust attenuation curves evolve with redshift and quantify the dependence of the inferred galaxy properties on the assumed dust attenuation law. We develop a modified version of the SED fitting code BAGPIPES by including a detailed dust attenuation curve parametrization. Dust parameters are derived, along with galaxy properties, from the fit to the data from FUV to mm bands. Once applied to three star-forming galaxies at z = 7-8, we find that their attenuation curves differ from local templates. One out of three galaxies shows a characteristic MW bump, typically associated to the presence of small carbonaceous dust grains such as PAHs. This is one of the first evidences suggesting the presence of PAHs in early galaxies. Galaxy properties such as stellar mass and SFR inferred from SED fitting are strongly affected by the assumed attenuation curve, though the adopted star formation history also plays a major role. Our results highlight the importance of accounting for the potential diversity of dust attenuation laws when analyzing the properties of galaxies at the EoR, whose dust properties are still poorly understood. The application of our method to a larger sample of galaxies observed with JWST can provide us important insights into the properties of dust and galaxies in the early universe.

Fast radio bursts (FRBs) are millisecond-radio transients observed at cosmological distances. The nature of their progenitors is still under debate, although magnetars are invoked by most models. The FRB-magnetar connection was strengthened by the discovery of an FRB-like event from the Galactic magnetar SGR J1935+215. In this work, we aim to investigate how prevalent are magnetars like SGR~J1935+2154 within FRB progenitors. We carried out an FRB search in a sample of seven nearby (< 12 Mpc) galaxies with the Northern Cross radio telescope for a total of 692 h. We detected one 1.8~ms burst in the direction of M101 with fluence of $58 \pm 5$ Jy ms. Its dispersion measure of 303 pc cm$^{-3}$ places it most-likely beyond M101. Considering no significant detection coming indisputably from the selected galaxies, we place a 38 yr$^{-1}$ upper limit on the total - i.e., including the whole sample - burst rate at the 95\% confidence level. This upper limit constrains $\lambda_{\rm mag} < 0.42$~magnetar$^{-1}$ yr$^{-1}$ or, if combined with literature observations of a similar sample of nearby galaxies, it yields a joint constraint of $\lambda_{\rm mag} < 0.25$ magnetar$^{-1}$ yr$^{-1}$. We also provide the first constraints on the expected rate of FRBs hypothetically originating from ultraluminous X-ray (ULX) sources, since some of the galaxies observed during our observational campaign host confirmed ULXs. We obtain $< 13$ yr$^{-1}$ per ULX on the total sample of galaxies observed. Our results indicate that bursts with energies $E > 10^{34}$ erg from magnetars like SGR~J1935+2154 appear more rare compared to previous observations, and further disfavours them as unique progenitors for the cosmological FRB population, leaving more space open to the contribution from a population of more exotic magnetars, not born via core-collapsed supernovae.

We use JWST/NIRCam medium band photometry in a single pointing of the CAnadian NIRISS Unbiased Cluster Survey (CANUCS) to identify 118 Extreme Emission Line Galaxies (EELGs) over $1.7 \lesssim z \lesssim 6.7$, selected using a set of color cuts that target galaxies with extreme $\text{[OIII] + H}\beta$ and H$\alpha$ emission. We show that our medium band color selections are able to select galaxies based on emission line equivalent width (EW), which is advantageous to more commonly used selections since it does not require strong continuum emission, and can select galaxies with faint or red continuum fluxes. The median EWs of our sample is $EW(\text{H}\alpha) = 893 $ \AA\ and $ EW(\text{[OIII] + H}\beta) = 1255 $ \AA, and includes some objects with $EW(\text{[OIII] + H}\beta) \sim 3000$ \AA. These systems are mostly compact with low stellar mass (median $\log(M_\star/M_\odot) = 8.03$), low metallicity (median $Z = 0.14 Z_\odot$), little dust (median $A_V = 0.18$ mag) and high SSFR (median $SSFR = 1.18 \times 10^{-8}/yr$). Additionally, galaxies in our sample show increasing EW(\Ha) and EW(\OIIIHb) with redshift, an anti-correlation of EW(\Ha) with stellar mass, and no correlation between EW(\OIIIHb) and stellar mass. Finally, we present NIRSpec spectroscopy of 15 of the EELGs in our sample. These spectra confirm the redshifts and EWs of the EELGs calculated from the medium bands, which demonstrates the accuracy and efficiency of our color selections. Overall, we show that there are significant advantages to using medium band photometry to identify and study EELGs at a wide range of redshifts.

We explore the potential of the ongoing Zwicky-Transient-Facility (ZTF) survey for studying Intra-Night Optical Variability (INOV) of active galactic nuclei (AGN), in particular for picking rare events of large INOV amplitudes, whose detection may require extensive temporal coverage. For this, we have used the available high cadence subsets of the ZTF database to build a well-defined large sample of 53 blazars (BLs) and another sample of 132 radio-quiet quasars (RQQs), matched to the blazar sample in the redshift-magnitude plane. High-cadence ZTF monitoring of these two matched samples are available, respectively, for 156 and 418 intranight sessions. Median durations for both sets of sessions are 3.7 hours. The two classes of powerful AGN monitored in these sessions represent opposite extremes of jet activity. The present analysis of their ZTF light curves has revealed some strong INOV events which, although not exceptionally rare for blazars, are indeed so for RQQs, and their possible nature is briefly discussed.

The Event Horizon Telescope (EHT) has led to the first images of a supermassive black hole, revealing the central compact objects in the elliptical galaxy M87 and the Milky Way. Proposed upgrades to this array through the next-generation EHT (ngEHT) program would sharply improve the angular resolution, dynamic range, and temporal coverage of the existing EHT observations. These improvements will uniquely enable a wealth of transformative new discoveries related to black hole science, extending from event-horizon-scale studies of strong gravity to studies of explosive transients to the cosmological growth and influence of supermassive black holes. Here, we present the key science goals for the ngEHT and their associated instrument requirements, both of which have been formulated through a multi-year international effort involving hundreds of scientists worldwide.

Constraining when and how reionisation began is pivotal for understanding when the first galaxies formed. Lyman-alpha (Ly$\alpha$) emission from galaxies is currently our most promising probe of these early stages. At z>7 the majority of galaxies detected with Ly$\alpha$ are in candidate overdensities. Here we quantify the probability of these galaxies residing in large ionised bubbles. We create (1.6 Gpc)$^3$ reionising intergalactic medium (IGM) simulations, providing sufficient volume to robustly measure bubble size distributions around UV-bright galaxies and rare overdensities. We find $M_{\rm UV} \lesssim -16$ galaxies and overdensities are $\gtrsim$10-1000x more likely to trace ionised bubbles compared to randomly selected positions. The brightest galaxies and strongest overdensities have bubble size distributions with highest characteristic size and least scatter. We compare two models: gradual reionisation driven by numerous UV-faint galaxies versus more rapid reionisation by rarer brighter galaxies, producing larger bubbles at fixed neutral fraction. We demonstrate that recently observed z~7 overdensities are highly likely to trace large ionised bubbles, corroborated by their high Ly$\alpha$ detection rates. However, the z~8.7 association of Ly$\alpha$ emitters in EGS and GN-z11, with Ly$\alpha$ at z=10.6, are unlikely to trace large bubbles in our fiducial model -- 11% and 7% probability of >1 proper Mpc bubbles, respectively. Ly$\alpha$ detections at such high redshifts could be explained by: a less neutral IGM than previously expected; larger ionised regions at fixed neutral fraction; or if intrinsic Ly$\alpha$ flux is unusually strong in these galaxies. We discuss how to test these scenarios with JWST and the prospects for using upcoming wide-area surveys to distinguish between reionisation models.

We use Gemini integral field unit observations to map the stellar population properties in the inner region ($\sim680\times470$ pc$^2$) of the galaxy NGC 6868. In order to understand the physical and chemical properties of the stellar content of this galaxy, we performed stellar population synthesis using the starlight code with the MILES simple stellar population models. We measured the absorption line indices Fe4383, Mg$_2$, Mg$_b$, Fe5270, Fe5335 for the whole FoV, and used them to derive Fe3 and [MgFe]'. These indices were used to derive [$\alpha$/Fe]. This galaxy is dominated by old metal-rich populations (12.6 Gyr; 1.0 and 1.6 Z$_\odot$) with a negative metallicity gradient. We also found a recent ($\sim63$ Myr) metal-rich (1.6 Z$_{\odot}$) residual star formation in the centre of the galaxy. A dust lane with a peak extinction in the V band of 0.65 mag is seen. No signs of ordered stellar motion are found and the stellar kinematics is dispersion dominated. All indices show a spatial profile varying significantly along the FoV. Mg$_2$ shows a shallow gradient, compatible with the occurrence of mergers in the past. Mg$_b$ and Fe3 profiles suggest different enrichment processes for these elements. We observe three distinct regions: for $R<100$pc and $R>220$pc, Mg$_2$, Mg$_b$ anti correlate with respect to Fe3 and [MgFe]', and for $100 \text{pc}<R<220 \text{pc}$, they correlate, hinting at different enrichment histories. The [$\alpha$/Fe] profile is really complex and has a central value of $\sim 0.2$ dex. We interpret this as the result of a past merger with another galaxy with a different [$\alpha$/Fe] history, thus explaining the [$\alpha$/Fe] maps.

We show that the maximum shower depth ($X_{\rm max}$) distributions of Ultra-High Energy Cosmic Rays (UHECRs), as measured by fluorescence telescopes, can be augmented by building a mapping to observables collected by surface detectors. Using the publicly available data on "golden hybrid'' events from the Pierre Auger Observatory we demonstrate significant correlations between $X_{\rm max}$ and timing information from ground Cherenkov detectors. Using such a mapping we show how to incorporate a subset of ground data into the inference of the $X_{\rm max}$ distribution, where the size of this subset depends on the strength of the correlation found. With a simple linear fit model, we are able to effectively incorporate $\sim13\%$ of all ground data statistics. Finally, we use this augmented dataset to infer the composition of UHECRs and discriminate between hadronic models used in air shower development simulations, and show the results improve significantly due to the effectively larger statistics available.

The relation between the total mass contained in the globular clusters of a galaxy and the mass of its dark matter halo has been found observationally to be nearly linear over five decades of mass. However, the high-mass end of this relation is not well determined from previous data and shows large scatter. We analyze the globular cluster systems (GCSs) of a homogeneous sample of 11 brightest cluster galaxies (BCGs) through DOLPHOT photometry of their deep Hubble Space Telescope (HST) images in the F814W filter. We standardize the definition of $M_{GCS}$, the total GCS mass, by using the GC total population within a limiting radius of $0.1 R_{virial}$, while the dark-matter halo mass $M_h$ is determined from the weak-lensing calibration of $M_h$ versus $M_{bary}$. When these 11 BCGs are added to the previously studied homogeneous catalogue of Virgo member galaxies, a total value for $\eta = M_{GCS}/M_h$ is found to be $(3.0\pm1.8_{internal})\times10^{-5}$, slightly higher than previous estimates but with much reduced uncertainty. Perhaps more importantly, the results suggest that the relation continues to have a near-linear shape at the highest galaxy masses, strongly reinforcing the conclusion that accreted GCs make a major contribution to the GC populations at high galaxy mass.

Ultra-/hyperluminous X-ray sources (ULX/HLX) could be interesting laboratories to further improve our understanding of the supermassive black hole growth through super-Eddington accretion episodes and successive mergers of lighter holes. ULXs are thought to be powered by super-Eddington accretion onto stellar-mass compact objects, while HLXs may be accreting intermediate mass black holes (IMBH). However, a significant portion of the sample of ULX/HLX candidates derived from catalogue searches are background AGN. Here we build ULX and HLX samples from recent XMM-Newton, Swift-XRT and Chandra catalogues and the GLADE catalogue of galaxies. We aim to characterise the frequency, environment, hardness and variability of ULXs and HLXs to better assess their differences and understand their populations. After a thorough classification of these X-ray sources, we remove 42% of $S/N>3$ sources shown to be contaminants, to obtain the cleanest sample of ULX/HLX to date. From a sample of 1342 ULXs and 191 HLXs, we study the occupation fraction, hardness, variability, radial distribution and preferred environment of the sources. We build their Malmquist-corrected X-ray luminosity functions (XLF) and compare them with previous studies. We statistically compare ULXs and HLXs and assess the differences in their nature. The interpretation of HLXs as IMBHs is investigated. A significant break is seen in the XLF at $\sim 10^{40}$ erg/s. Our ULX sample, having $\leq 2$% of contaminants, confirms that ULXs are located preferentially in spiral galaxies and galaxies with higher star-formation rates. Unlike ULXs, HLXs seem to reside equally in spiral and lenticular/elliptical galaxies. 35% of the HLX candidates have an optical counterpart, and we estimate the mass of 120 of them in the range of $2000-10^5 M_\odot$. Most HLXs are consistent with an accreting massive black hole in a dwarf galaxy satellite.

What is the highest energy at which gravitons can be observed? We address this question by studying graviton-to-photon conversion - the inverse-Gertsenshtein effect - in the magnetic field of the Milky Way. We find that above $\sim 1~\mbox{PeV}$ the effective photon mass grows large enough to quench the conversion rate. The induced photon flux is comparable to the sensitivity of LHAASO to a diffuse $\gamma$-ray background, but only for graviton abundances of order $\Omega_{\text{gw}} h^2_0 \sim 1$. In the future, owing to a better understanding of $\gamma$-ray backgrounds, larger effective areas and longer observation times, sub-PeV shimmering gravitons with a realistic abundance of $\Omega_{\text{gw}} h^2_0 \sim 0.01$ could be detected. We show that this is achieved in a cosmologically-motivated scenario of post-recombination superheavy dark matter decay. Therefore, the sub-PeV range might be the ultimate energy frontier at which gravitons can be observed.

We present the first results study of the effects of the powerful Gamma Ray Burst GRB 221009A that occurred on October 9, 2022, and was serendipitously recorded by electron and proton detectors aboard the four spacecraft of the NASA THEMIS mission. Long-duration gamma-ray bursts (GRBs) are powerful cosmic explosions, signaling the death of massive stars, and, among them, GRB 221009A is so far the brightest burst ever observed due to its enormous energy ($E_{\gamma iso}\sim10^{55}$ erg) and proximity (the redshift is $z\sim 0.1505$). The THEMIS mission launched in 2008 was designed to study the plasma processes in the Earth's magnetosphere and the solar wind. The particle flux measurements from the two inner magnetosphere THEMIS probes THA and THE and ARTEMIS spacecraft THB and THC orbiting the Moon captured the dynamics of GRB 221009A with a high-time resolution of more than 20 measurements per second. This allowed us to resolve the fine structure of the gamma-ray burst and determine the temporal scales of the two main bursts spiky structure complementing the results from gamma-ray space telescopes and detectors.

We present a charge distribution-based emission model that calculates the infrared spectrum of fullerenes (C$_{60}$). Analysis of the modelled spectrum of C$_{60}$ in various charge states shows that the relative intensity of the features in the 5-10 $\mu$m versus 15-20 $\mu$m can be used to probe the C$_{60}$ charge state in interstellar spectra. We further used our model to simulate emission from polycyclic aromatic hydrocarbons (PAHs) and C$_{60}$ at five positions in the cavity of reflection nebula NGC~7023. Specifically, we modelled the 6.2/11.2 band ratio for circumcoronene and circumcircumcoronene and the 7.0/19.0 band ratio for C$_{60}$ as a function of the ionization parameter $\gamma$. A comparison of the model results with the observed band ratios shows that the $\gamma$ values in the cavity do not vary significantly, suggesting that the emission in the cavity does not originate from locations at the projected distances. Furthermore, we find that the C$_{60}$ derived $\gamma$ values are lower than the PAH-derived values by an order of magnitude. We discuss likely scenarios for this discrepancy. In one scenario, we attribute the differences in the derived $\gamma$ values to the uncertainties in the electron recombination rates of PAHs and C$_{60}$. In the other scenario, we suggest that PAHs and C$_{60}$ are not co-spatial resulting in different $\gamma$ values from their respective models. We highlight that experiments to determine necessary rates will be required in validating either one of the scenarios.

It remains an elusive goal to simultaneously model the astrophysics of Solar System accretion while reproducing the mantle chemistry of more than one inner terrestrial planet. Here, we used a multistage core-mantle differentiation model based on Rubie et al. (2011,2015) to track the formation and composition of Earth and Mars in various Grand Tack formation simulations. Prior studies showed that in order to recreate Earth's mantle composition, it must grow first from reduced (Fe-metal rich and O-poor) building blocks and then from increasingly oxidized (FeO rich) material. This accretion chemistry occurs when an oxidation gradient exists across the disk so that the innermost solids are reduced and increasingly oxidized material is found at greater heliocentric distances. For a suite of Grand Tack simulations, we investigated whether Earth and Mars can be simultaneously produced by the same oxidation gradient. Our model did not find an oxidation gradient that simultaneously reproduces the mantle composition of Earth and Mars. Due to its small mass and rapid formation, the formation history of Mars-like planets is very stochastic which decreases the likelihood of compatibility with an Earth-producing oxidation gradient in any given realization. To reconcile the accretion history and ideal chemistry of the Mars-like planet with the oxidation gradient of an Earth-producing disk, we determined where in the Earth-producing disk Mars must have formed. We find that the FeO-rich composition of the Martian mantle requires that Mars' building blocks must originate exterior to 1.0 astronomical units (AU).

We present a photometric analysis of globular cluster 47 Tuc (NGC\,104), using near-IR imaging data from the GeMS/GSAOI Galactic Globular Cluster Survey (G4CS) which is in operation at Gemini-South telescope.~Our survey is designed to obtain AO-assisted deep imaging with near diffraction-limited spatial resolution of the central fields of Milky Way globular clusters.~The G4CS near-IR photometry was combined with an optical photometry catalog obtained from Hubble Space Telescope survey data to produce a high-quality color-magnitude diagram that reaches down to K$_s\approx$ 21 Vega mag.~We used the software suite BASE-9, which uses an adaptive Metropolis sampling algorithm to perform a Markov chain Monte Carlo (MCMC) Bayesian analysis, and obtained probability distributions and precise estimates for the age, distance and extinction cluster parameters.~Our best estimate for the age of 47 Tuc is 12.42$^{+0.05}_{-0.05}$ $\pm$ 0.08 Gyr, and our true distance modulus estimate is (m$-$M)$_0$=13.250$^{+0.003}_{-0.003}$ $\pm$ 0.028 mag, in tight agreement with previous studies using Gaia DR2 parallax and detached eclipsing binaries.

The stellar dynamic-based black hole mass measurement of M87 is twice that determined via ionized gas kinematics; they disagree by more than 3$\sigma$, being the former closer to the mass estimated from the diameter of the gravitationally-lensed ring around the black hole. Using a deeper and more comprehensive ionized gas kinematic dataset, we aim to better constrain the complex morphology and kinematics of the nuclear ionized gas. We use both Narrow and Wide Field Mode integral field spectroscopic data from the Multi Unit Spectroscopic Explorer instrument to model the kinematics of multiple ionized gas emission lines. The dataset reveals complexities in the nuclear ionized gas kinematics not seen in earlier observations. Several ionized gas filaments, some with large flow velocities, can be traced down into the projected sphere of influence. We also find evidence of a partially-filled biconical outflow, aligned with the jet. The nuclear rotating ionized gas 'disk' is well resolved in our datacubes. The complexity of the nuclear morphology precludes the measurement of an accurate black hole mass from the ionized gas kinematics. The ionized gas kinematics can be explained with either a 6.6 $\times 10^{9}\rm~M_{\odot}$ black hole with a disk inclination of 25\deg, or a 3.5 $\times 10^{9}\rm~M_{\odot}$ black hole with an inclination of 42\deg, with more support for the former. The velocity isophotes of the sub-arcsecond ionized gas disk are twisted, and the position angle of the innermost gas disk tends towards a value perpendicular to the radio jet axis. An inclination close to 25\deg\ for the nuclear gas disk, and the warp in the sub-arcsec ionized gas disk, help to reconcile the contradictory nature of key earlier results: the discrepancy between stellar and ionized gas measurements, and the mis-orientation between the axes of the ionized gas disk and the jet.

Context. The second subclass of chemically peculiar stars, the CP2 stars, are early-type stars exhibiting anomalous abundances with variable line strengths possibly also accompanied by photometric variability that typically belong to the Galactic disk. However, a small fraction of these objects were recently found to be located far from the Galactic plane and are thought to belong to the Galactic halo, which is unexpected for this class of objects. Aims. Our study investigates the halo membership of the proposed halo CP2 star candidates based on their velocities and Galactic orbits, to determine their points of origin. In addition, we searched for further halo CP2 star candidates using an as-yet-unpublished sample of CP stars. Methods. Our analysis relied on the astrometry from the early third \textit{Gaia} data release and radial velocities based on the spectroscopy from LAMOST and SDSS, which was also employed in spectral classification. The light variability of the confirmed CP2 stars in our sample was analyzed using data from the ZTF and ATLAS surveys. Results. After filtering our initial sample using kinematic and spectroscopic criteria, we identified six CP2 stars with kinematical properties consistent with a halo membership. The orbits of these stars are in agreement with an origin in the Galactic disk where they were probably ejected through dynamical interactions or in the binary supernova scenario, making them the first runaway CP2 stars known.

We perform three-dimensional numerical simulations to study the hydrodynamic evolution of Sgr A East, the only known supernova remnant (SNR) in the center of our Galaxy, to infer its debated progenitor SN type and its potential impact on the Galactic center environment. Three sets of simulations are performed, each of which represents a represent a certain type of SN explosion (SN Iax, SN Ia or core-collapse SN) expanding against a nuclear outflow of hot gas driven by massive stars, whose thermodynamical properties have been well established by previous work and fixed in the simulations. All three simulations can simultaneously roughly reproduce the extent of Sgr A East and the position and morphology of an arc-shaped thermal X-ray feature, known as the "ridge". Confirming previous work, our simulations show that the ridge is the manifestation of a strong collision between the expanding SN ejecta and the nuclear outflow. The simulation of the core-collapse SN, with an assumed explosion energy of 5x10^50 erg and an ejecta mass of 10 M_sun, can well match the X-ray flux of the ridge, whereas the simulations of the SN Iax and SN Ia explosions underpredict its X-ray emission, due to a smaller ejecta mass. All three simulations constrain the age of Sgr A East to be <1500 yr and predict that the ridge should fade out over the next few hundred years. We address the implications of these results for our understanding of the Galactic center environment.

Supra-arcade downflows (SADs) are dark voids descending towards the post-reconnection flare loops and exhibit obvious variation in properties like width. However, due to the lack of further statistical studies, the mechanism behind such variations hitherto remains elusive. Here we statistically investigated widths of 81 SADs observed in one flare by the Solar Dynamics Observatory (SDO). For each of SADs, six moments were selected with equal time intervals to measure their widths at different stages of their evolution. It is found that most SADs show a roughly monotonous width decrease during their descents, while some SADs with small initial widths can have complex evolutions. 3D reconstruction results based on SDO and Solar Terrestrial Relations Observatory Ahead (STEREO-A) images and thermal properties analysis reveal that differences in magnetic and plasma environments may result in that SADs in the north are overall wider than those in the south. Additionally, correlation analysis between the width and other parameters of SADs was further conducted and revealed that: (1) SADs with different initial widths show no significant differences in their temperature and density evolution characteristics; (2) SADs with small initial widths usually appear in lower heights, where more frequent collisions between SADs could lead to their intermittent acceleration, width increment, and curved trajectories. These results indicate that SADs with different initial widths are produced the same way while different environments (magnetic field or plasma) could affect their subsequent width evolutions.

Core-collapse supernovae are a useful laboratory to probe the nature of exotic particles. If axionlike particles (ALPs) are produced in supernovae, they can affect the transfer of energy and leave traces in observational signatures. In this work, we develop two-dimensional supernova models including the effects of the production and the absorption of ALPs that couple with photons. It is found that the additional heating induced by ALPs can enhance the explosion energy E_exp; for moderate ALP-photon coupling, we find explosion energies ~0.6*10^51 erg compared to our reference model without ALPs of ~0.4*10^51 erg. Our findings also indicate that when the coupling constant is sufficiently high, the neutrino luminosities and mean energies are decreased because of the additional cooling of the proto-neutron star. The gravitational wave strain is also reduced because the mass accretion on the proto-neutron star is suppressed. Although the ALP-photon coupling can foster explodability, including enhancing the explosion energy closer to recent observations, more long-term simulations in spatially three-dimension are needed to draw robust conclusions.

Ongoing and upcoming galaxy surveys are providing precision measurements of galaxy clustering. However a major obstacle in its cosmological application is the stochasticity in the galaxy bias. We explore whether the principal component analysis (PCA) of galaxy correlation matrix in hyperspace of galaxy properties (e.g. magnitude and color) can reveal further information on mitigating this issue. Based on the hydrodynamic simulation TNG300-1, we analyze the cross power spectrum matrix of galaxies in the magnitude and color space of multiple photometric bands. (1) We find that the first principal component $E_i^{(1)}$ is an excellent proxy of the galaxy deterministic bias $b_{D}$, in that $E_i^{(1)}=\sqrt{\lambda^P(1)/P_{mm}}b_{D,i}$. Here $i$ denotes the $i$-th galaxy sub-sample. $\lambda^{(1)}$ is the largest eigenvalue and $P_{mm}$ is the matter power spectrum. We verify that this relation holds for all the galaxy samples investigated, down to $k\sim 2h/$Mpc. Since $E_i^{(1)}$ is a direct observable, we can utilize it to design a linear weighting scheme to suppress the stochasticity in the galaxy-matter relation. For an LSST-like magnitude limit galaxy sample, the stochasticity $\mathcal{S}\equiv 1-r^2$ can be suppressed by a factor of $\ga 2$ at $k=1h/$Mpc. This reduces the stochasticity-induced systematic error in the matter power spectrum reconstruction combining galaxy clustering and galaxy-galaxy lensing from $\sim 12\%$ to $\sim 5\%$ at $k=1h/$Mpc. (2) We also find that $\mathcal{S}$ increases monotonically with $f_\lambda$ and $f_{\lambda^2}$. $f_{\lambda,\lambda^2}$ quantify the fractional contribution of other eigenmodes to the galaxy clustering and are direct observables. Therefore the two provide extra information on mitigating galaxy stochasticity.

Massive black holes are fundamental constituents of our cosmos, from the Big Bang to today. Understanding their formation from cosmic dawn, their growth, and the emergence of the first, rare quasars in the early Universe remains one of our greatest theoretical and observational challenges. Hydrodynamic cosmological simulations self-consistently combine the processes of structure formation at cosmological scales with the physics of smaller, galaxy scales. They capture our most realistic understanding of massive black holes and their connection to galaxy formation and have become the primary avenue for theoretical research in this field. The space-based gravitational wave interferometer, LISA, will open up new investigations into the dynamical processes involving massive black holes. Multi-messenger astrophysics brings new exciting prospects for tracing the origin, growth and merger history of massive black holes across cosmic ages.

With thousands of exoplanets now identified, the characterization of habitable planets and the potential identification of inhabited ones is a major challenge for the coming decades. We review the current working definition of habitable planets, the upcoming observational prospects for their characterization and present an innovative approach to assess habitability and inhabitation. This integrated method couples for the first time the atmosphere and the interior modeling with the biological activity based on ecosystem modeling. We review here the first applications of the method to asses the likelihood and impact of methanogenesis for Enceladus, primitive Earth, and primitive Mars. Informed by these applications for solar system situations where habitability and inhabitation is questionned, we show how the method can be used to inform the design of future space observatories by considering habitability and inhabitation of Earth-like exoplanets around sun-like stars.

We investigate how the diverse star formation histories observed across galaxy masses emerged using models that evolve under gas accretion from host halos. They also include ejection of interstellar matter by supernova feedback, recycling of ejected matter and preventive feedback that partially hinders gas accretion. We consider three schemes of gas accretion: the fiducial scheme which includes the accretion of cold gas in low-mass halos and high-redshift massive halos as hinted by cosmological simulations; the flat scheme in which high-mass cold accretion is removed; and finally the shock-heating scheme which assumes radiative cooling of the shock-heated halo gas. The fiducial scheme reproduces dramatic diminishment in star formation rate (SFR) after its peak as observed for the present halo mass $M_{\rm vir}>10^{12.5}{\rm M}_\odot$ while other two schemes show reduced or negligible quenching. This scheme reproduces the high-mass slope in the SFR vs. stellar mass relation decreasing toward recent epochs whereas other two schemes show opposite trend which contradicts observation. Success in the fiducial scheme originates in the existence of high-mass cold-mode accretion which retards transition to the slow hot-mode accretion thereby inducing a larger drop in SFR. Aided by gas recycling, which creates monotonically increasing SFR in low-mass halos, this scheme can reproduce the downsizing galaxy formation. Several issues remain, suggesting non-negligible roles of missing physics. Feedback from active galactic nuclei could mitigate upturn of SFR in low-redshift massive halos whereas galaxy mergers could remedy early inefficient star formation.

Artificial intelligence methods are indispensable to identifying pulsars from large amounts of candidates. We develop a new pulsar identification system that utilizes the CoAtNet to score two-dimensional features of candidates, uses a multilayer perceptron to score one-dimensional features, and uses logistic regression to judge the scores above. In the data preprocessing stage, we performed two feature fusions separately, one for one-dimensional features and the other for two-dimensional features, which are used as inputs for the multilayer perceptron and the CoAtNet respectively. The newly developed system achieves 98.77\% recall, 1.07\% false positive rate and 98.85\% accuracy in our GPPS test set.

In this study, we employ and modify the Lorenz energy cycle (LEC) framework as another way to understand the atmospheric circulation on tidally locked terrestrial planets. It well describes the atmospheric general circulation in the perspective of energy transformation, involved with several dynamical processes. We find that on rapidly rotating, tidally locked terrestrial planets, mean potential energy (P$_{\rm M}$) and eddy potential energy (P$_{\rm E}$) are comparable to those on Earth, as they have similar steep meridional temperature gradients. Mean kinetic energy (K$_{\rm M}$) and eddy kinetic energy (K$_{\rm E}$) are larger than those on Earth, related to stronger winds. The two conversion paths, P$_{\rm M}\rightarrow$P$_{\rm E}\rightarrow$K$_{\rm E}$ and P$_{\rm M}\rightarrow$K$_{\rm M}\rightarrow$K$_{\rm E}$, are both efficient. The former is associated with strong baroclinic instabilities, and the latter is associated with Hadley cells. On slowly rotating, tidally locked terrestrial planets, weak temperature gradients in the free atmosphere and strong nightside temperature inversion make P$_{\rm M}$ and P$_{\rm E}$ are much smaller than those on Earth. Meanwhile, large day--night surface temperature contrast and small rotation rate make the overturning circulation extend to the globe, so that the main conversion path is P$_{\rm M}\rightarrow$K$_{\rm M}\rightarrow$K$_{\rm E}$. This study shows that the LEC analyses improve the understanding of the atmospheric circulation on tidally locked terrestrial planets.

This is the first paper in our series on the combined analysis of the Dynamics and stellar Population (DynPop) for the MaNGA survey in the final SDSS Data Release 17 (DR17). Here we present a catalogue of dynamically-determined quantities for over 10000 nearby galaxies based on integral-field stellar kinematics from the MaNGA survey. The dynamical properties are extracted using the axisymmetric Jeans Anisotropic Modelling (JAM) method, which was previously shown to be the most accurate for this kind of study. We assess systematic uncertainties using eight dynamical models with different assumptions. We use two orientations of the velocity ellipsoid: either cylindrically-aligned JAM$_{\rm cyl}$ or spherically-aligned JAM$_{\rm sph}$. We also make four assumptions for the models' dark vs. luminous matter distributions: (1) mass-follows-light, (2) free NFW dark halo, (3) cosmologically-constrained NFW halo, (4) generalized NFW dark halo, i.e. with free inner slope. In this catalogue, we provide the quantities related to the mass distributions (e.g. the density slopes and enclosed mass within a sphere of a given radius for total mass, stellar mass, and dark matter mass components). We also provide the complete models which can be used to compute the full luminous and mass distribution of each galaxy. Additionally, we visually assess the qualities of the models to help with model selections. We estimate the observed scatter in the measured quantities which decreases as expected with improvements in quality. For the best data quality, we find a remarkable consistency of measured quantities between different models, highlighting the robustness of the results.

We analyze the global stellar population, radial gradients and non-parametric star formation history of $\sim 10$K galaxies from the MaNGA Survey final data release 17 (DR17), based on stellar population synthesis and full-spectrum fitting, and relate them with dynamical properties of galaxies. We confirm that stellar population correlates with stellar velocity dispersion $\sigma_{\rm e}$ better than with stellar mass $M_{\ast}$, but also find that younger galaxies are more metal-poor at fixed $\sigma_{\rm e}$. Stellar age, metallicity, and mass-to-light ratio $M_{\ast}/L$ all decrease with galaxy rotation, while radial gradients become more negative (i.e., younger, more metal-poor, and lower $M_{\ast}/L$ in the outskirts). The trend between metallicity gradients and rotation reverses for slow rotators, which stand out for their more negative metallicity gradients than faster-rotating galaxies. We highlight a population of massive disk galaxies on the green valley, on the $(\sigma_{\rm e},\rm Age)$ plane, that show steep negative age and metallicity gradients, consistent with their old central bulges surrounded by young star-forming disks and metal-poor gas accretion. Galaxies with high $\sigma_{\rm e}$, steep total mass-density slope, low dark matter fraction, high $M_{\ast}/L$, and high metallicity have the highest star-formation rate at earlier times, and are currently quenched. We discover a population of low-mass star-forming galaxies with low rotation but physically distinct from the massive slow rotators. A catalogue of the population properties is provided publicly.

We present dynamical scaling relations, combined with the stellar population properties, for a subsample of about 6000 nearby galaxies with the most reliable dynamical models extracted from the full MaNGA sample of 10K galaxies. We show that the inclination-corrected mass plane (MP) for both early-type galaxies (ETGs) and late-type galaxies (LTGs), which links dynamical mass, projected half-light radius $R_{\rm e}$, and the second stellar velocity moment $\sigma_{\rm e}$ within $R_{\rm e}$, satisfies the virial theorem and is even tighter than the uncorrected one. We find a clear parabolic relation between $\lg(M/L)(<R_{\rm e})$, the total mass-to-light ratio within a sphere of radius $R_{\rm e}$, and $\lg\sigma_{\rm e}$, with the $M/L$ increasing with $\sigma_{\rm e}$ and for older stellar populations. However, the relation for ETGs is linear and the one for the youngest galaxies is constant. We confirm and improve the relation between average logarithmic total density slopes $\gamma_{_{\rm T}}$ and $\sigma_{\rm e}$: $\gamma_{_{\rm T}}$ become steeper with increasing $\sigma_{\rm e}$ until $\lg(\sigma_{\rm e}/{\rm km\,s^{-1}})\approx 2.2$ and then remain constant around $\gamma_{_{\rm T}}\approx -2.2$. The $\gamma_{_{\rm T}}-\sigma_{\rm e}$ variation is larger for LTGs than ETGs. At fixed $\sigma_{\rm e}$ the total density profiles steepen with galaxy age and for ETGs. We find generally low dark matter fractions, median $f_{\rm DM}(<R_{\rm e})=8$ per cent, within a sphere of radius $R_{\rm e}$. However, we find that $f_{\rm DM}(<R_{\rm e})$ depends on $\sigma_{\rm e}$ better than stellar mass: dark matter increases to a median $f_{\rm DM}=33$ percent for galaxies with $\sigma_{\rm e}\lesssim100{\rm km\,s^{-1}}$. The increased $f_{\rm DM}(<R_{\rm e})$ at low $\sigma_{\rm e}$ explains the parabolic $\lg(M/L)(<R_{\rm e})-\lg\sigma_{\rm e}$ relation.

We present the measurement of total and stellar/dark matter decomposed mass density profile around a sample of galaxy groups and clusters with dynamical masses derived from integral-field stellar kinematics from the MaNGA survey in Paper~I and weak lensing derived from the DECaLS imaging survey. Combining the two data sets enables accurate measurement of the radial density distribution from several kpc to Mpc scales. Intriguingly, we find that the excess surface density derived from stellar kinematics in the inner region cannot be explained by simply adding an NFW dark matter halo extrapolated from lensing measurement at a larger scale to a stellar mass component derived from the NASA-Sloan Atlas (NSA) catalog. We find that a good fit to both data sets requires a stellar mass normalization about 3 times higher than that derived from the NSA catalog, which would require an unrealistically too-heavy initial mass function for stellar mass estimation. If we keep the stellar mass normalization to that of the NSA catalog but allow a varying inner dark matter density profile, we obtain an asymptotic slope of $\gamma_{\rm gnfw}$= $1.82_{-0.25}^{+0.15}$, $\gamma_{\rm gnfw}$= $1.48_{-0.41}^{+0.20}$ for the group bin and the cluster bin respectively, significantly steeper than the NFW case. We also compare the total mass inner density slopes with those from Illustris-TNG300 and find that the values from the simulation are lower than the observation by at least $3\sigma$ level.

We report the discovery of a pair of exoplanets co-orbiting the red dwarf star GJ3470. The larger planet, GJ3470-d, was observed in a 14.9617-days orbit and the smaller planet, GJ3470-e, in a 14.9467-days orbit. GJ3470-d is sub-Jupiter size with a 1.4% depth and a duration of 3 hours, 4 minutes. The smaller planet, GJ3470-e, currently leads the larger planet by approximately 1.146-days and is extending that lead by about 7.5-minutes (JD 0.0052) per orbital cycle. It has an average depth of 0.5% and an average duration of 3 hours, 2 minutes. The larger planet, GJ3470-d, has been observed on seven separate occasions over a 3-year period, allowing for a very precise orbital period calculation. The last transit was observed by three separate observatories in Oklahoma and Arizona. The smaller planet, GJ3470-e, has been observed on five occasions over 2-years. Our data appears consistent with two exoplanets in a Horseshoe Exchange orbit. When confirmed, these will be the second and third exoplanets discovered and characterized by amateur astronomers without professional data or assistance. It will also be the first ever discovery of co-orbiting exoplanets in a Horseshoe Exchange orbit.

Previous models have shown that stochastic scattering of stars in a two-dimensional galaxy disc can generate a time-independent surface density distribution that is an exponential divided by radius when a constant inward scattering bias is present. Here we show, using a Markov chain model, that similar profiles result from an outward scattering bias, although the disc surface density decreases slowly with time because of a net stellar outflow. The trend towards a near-exponential surface profile is robust, as it exists even if the scattering intensity has moderate radial and time dependences, subject to some limitations on the scattering rates discussed in the text. The exponential scale length of the pseudo-equilibrium disc depends on the scattering bias, the scattering length, and the size of the disc where scattering is important.

We present the serendipitous discovery of a large double radio relic associated with the merging galaxy cluster PSZ2 G277.93+12.34 and a new odd radio circle, ORC J1027-4422, both found in deep MeerKAT 1.3 GHz wide-band data. The angular separation of the two arc-shaped cluster relics is 16 arcmin or 2.6 Mpc for a cluster redshift of z = 0.158. The thin southern relic, which shows a number of ridges/shocks including one possibly moving inwards, has a linear extent of 1.64 Mpc. In contrast, the northern relic is about twice as wide, twice as bright, but only has a largest linear size of 0.66 Mpc. Complementary SRG/eROSITA X-ray images reveal extended emission from hot intracluster gas between the two relics and around the narrow-angle tail (NAT) radio galaxy PMN J1033-4335 (z = 0.153) located just east of the northern relic. No radio halo associated with the PSZ2 cluster is detected. The radio morphologies of the NAT galaxy and the northern relic, which are also detected with the Australian Square Kilometer Array Pathfinder at 887.5 MHz, suggest both are moving in the same outward direction. The discovery of ORC J1027-4422 in a different part of the MeerKAT image makes it the 4th known single ORC. It has a diameter of 90" corresponding to 400 kpc at a tentative redshift of z = 0.3 and remains undetected in X-ray emission. We discuss similarities between galaxy and cluster mergers as the formation mechanisms for ORCs and radio relics, respectively.

Preflare activities contain critical information about the pre-cursors and causes of solar eruptions. Here we investigate the characteristics and origin of a group of broadband pulsations (BBPs) in the decimetric-metric wavelengths, taking place during the preflare stage of the M7.1 flare dated on 2011 September 24. The event was recorded by multiple solar instruments including the Nan\c{c}ay Radioheliograh that measure the properties of the radio source. The BBPs start $\sim$24 min before the flare onset, extending from $<$ 360 to above 800 MHz with no discernible spectral drift. The BBPs consist of two stages, during the first stage the main source remains stationary, during the second stage it moves outward along with a steepening extreme-ultraviolet (EUV) wave driven by the eruption of a high-temperature structure. In both stages, we observe frequent EUV brightenings and jets originating from the flare region. During the second stage, the BBPs become denser in number and stronger in general, with the level of the polarization increasing gradually from $<$ 20% to $>$ 60% in the right-handed sense. These observations indicate the steepening EUV wave is important to the BBPs during the second stage, while the preflare reconnections causing the jets and EUV brightenings are important in both stages. This is the first time such a strong association of an EUV wave with BBPs is reported. We suggest a reconnection plus shock-sweeping-across-loop scenario for the cause of the BBPs.

The origination and generation mechanisms of small magnetic flux ropes (SFRs), which are important structures in solar wind, are not clearly known. In present study, 1993 SFRs immersed in coronal holes, active regions, and quiet Sun solar wind are analyzed and compared. We find that the properties of SFRs immersed in three types of solar wind are signicantly different. The SFRs are further classifed into hot-SFRs, cold-SFRs, and normal-SFRs, according to whether the O7+/O6+ is 30% elevated or dropped inside SFRs as compared with background solar wind. Our studies show that the parameters of normal-SFRs are similar to background in all three types of solar wind. The properties of hot-SFRs and cold-SFRs seem to be lying in two extremes. Statistically, the hot-SFRs (cold-SFRs) are associated with longer (shorter) duration, lower (higher) speeds and proton temperatures, higher (lower) charge states, helium abundance, and FIP bias as compared with normal-SFRs and background solar wind. The anti-correlations between speed and O7+/O6+ inside hot-SFRs (normal-SFRs) are different from (similar to) those in background solar wind. Most of hot-SFRs and cold-SFRs should come from the Sun. Hot-SFRs may come from streamers associated with plasma blobs and/or small-scale activities on the Sun. Cold-SFRs may be accompanied by small-scale eruptions with lower-temperature materials. Both hot-SFRs and cold-SFRs could also be formed by magnetic erosions of ICMEs that do not contain or contain cold-filament materials. The characteristics of normal-SFRs can be explained reasonably by the two originations, from the Sun and generated in the heliosphere both.

Stellar winds blown out from massive stars ($\gtrsim 10M_{\odot}$) contain precious information on the progenitor itself, and in this context, the most important elements are carbon (C), nitrogen (N), and oxygen (O), which are produced by the CNO cycle in the H-burning layer. Although their X-ray fluorescence lines are expected to be detected in swept-up shock-heated circumstellar materials (CSMs) in supernova remnants (SNRs), particularly those of C and N have been difficult to detect so far. Here, we present a high-resolution spectroscopy of a young magnetar-hosting SNR RCW~103 with the Reflection Grating Spectrometer (RGS) onboard XMM-Newton and report on the detection of \ion{N}{7} Ly$\alpha$ (0.50~keV) line for the first time. By comparing the obtained abundance ratio of N to O (N/O$=3.8 \pm{0.1}$) with various stellar evolution models, we show that the progenitor of RCW~103 is likely to have a low-mass (10--12~$M_{\odot}$) and medium-rotation velocities ($\lesssim 100~\rm{km~s^{-1}}$). The results also rule out the possibility of dynamo effects in massive ($\geq35~M_{\odot}$) stars as a formation mechanism of the associated magnetar 1E~161348$-$5055. Our method is useful for estimating various progenitor parameters for future missions with microcalorimeters such as XRISM and Athena.

With high-sensitivity kiloparsec-scale radio polarimetry, we can examine the jet-medium interactions and get a better understanding of the blazar divide in radio-loud (RL) AGN. We are analyzing the radio polarimetric observations with the EVLA and GMRT of 24 quasars and BL Lacs belonging to the Palomar-Green (PG) sample. The RL quasars show extensive polarisation structures in their cores, jets, lobes, and hotspots, whereas preliminary results suggest that BL Lacs exhibit polarisation primarily in their cores and inner jet regions. These findings imply that both intrinsic (central engine-related) and extrinsic (environment-related) variables are important in the formation of the blazar subclasses. The Fanaroff-Riley (FR) dichotomy can also be studied assuming RL unification and looking through the lens of blazars. Due to the radio-unbiased nature of the optically/UV-selected PG sample, we find a large fraction of the PG quasars are restarted, distorted (S- or X-shaped), or have a hybrid FR morphology.

The Five-hundred-meter Aperture Spherical radio Telescope (FAST) has been running for several years. A new Ultra-Wide Bandwidth (UWB) receiver, simultaneously covering 500-3300 MHz, has been mounted in the FAST feed cabin and passed a series of observational tests. The whole UWB band is separated into four independent bands. Each band has 1048576 channels in total, resulted in a spectral resolution of 1 kHz. At 500-3300 MHz, the antenna gain is around 14.3-7.7 K/Jy, the aperture efficiency is around 0.56-0.30, the system temperature is around 88-130 K, and the HPBW is around 7.6-1.6 arcmin. The measured standard deviation of pointing accuracy is better than ~7.9 arcsec, when zenith angle (ZA) is within 26.4deg. The sensitivity and stability of the UWB receiver are confirmed to satisfy expectation by spectral observations, e.g., HI and OH. The FAST UWB receiver already has a good performance for taking sensitive observations in various scientific goals.

We study how stably stratified or semi-convective layers alter the tidal dissipation rates associated with the generation of internal waves in planetary interiors. We consider if these layers could contribute to the high rates of tidal dissipation observed for Jupiter and Saturn in our solar system. We use an idealised global spherical Boussinesq model to study the influence of stable stratification and semi-convective layers on tidal dissipation rates. We carry out analytical and numerical calculations considering realistic tidal forcing and measure how the viscous and thermal dissipation rates depend on the parameters relating to the internal stratification profile. We find that the strongly frequency-dependent tidal dissipation rate is highly dependent on the parameters relating to the stable stratification, with strong resonant peaks that align with the internal modes of the system. The locations and sizes of these resonances depend on the form and parameters of the stratification, which we explore both analytically and numerically. Our results suggest that stable stratification can significantly enhance the tidal dissipation in particular frequency ranges. Analytical calculations in the low frequency regime give us scaling laws for the key parameters, including the tidal quality factor $Q'$ due to internal gravity waves. Stably stratified layers can significantly contribute to tidal dissipation in solar and extrasolar giant planets, and we estimate substantial tidal evolution for hot Neptunes. Further investigation is needed to robustly quantify the significance of the contribution in realistic interior models, and to consider the contribution of inertial waves.

In their recent study, Labb\'e et al. used multi-band infrared images captured by the James Webb Space Telescope (JWST) to discover a population of red massive galaxies that formed approximately 600 million years after the Big Bang. The authors reported an extraordinarily large density of these galaxies, with stellar masses exceeding $10^{10}$ solar masses, which, if confirmed, challenges the standard cosmological model as suggested by recent studies. However, this conclusion is disputed. We contend that during the early epochs of the universe the stellar mass-to-light ratio could not have reached the values reported by Labb\'e et al. A model of galaxy formation based on standard cosmology provides support for this hypothesis, predicting the formation of massive galaxies with higher ultraviolet (UV) luminosity, which produce several hundred solar masses of stars per year and containing significant dust. These forecasts are consistent with the abundance of JWST/HST galaxies selected photometrically in the rest-frame UV wavelengths and with the properties of the recent spectroscopically-confirmed JWST/HST galaxies formed during that era. Discrepancies with Labb\'e et al. may arise from overestimation of the stellar masses, systematic uncertainties, absence of JWST/MIRI data, heavy dust extinction affecting UV luminosities, or misidentification of faint red AGN galaxies at closer redshifts. The current JWST/HST results, combined with a realistic galaxy formation model, provide strong confirmation of the standard cosmology.

This work concentrates on the effect of an irrotational forcing on a magnetized flow in the presence of rotation, baroclinicity, shear, or a combination of them. By including magnetic field in the model we can evaluate the occurrence of dynamo on both small and large scales. We aim at finding what are the minimum ingredients needed to trigger a dynamo instability and what is the relation between dynamo and the growth of vorticity. We use the Pencil code to run resistive MHD direct numerical simulations. We report no dynamo in all cases where only rotation is included, regardless on the equation of state. Conversely, the inclusion of a background sinusoidal shearing profile leads to an hydrodynamic instability that produces an exponential growth of the vorticity at all scales, starting from small ones. This is know as vorticity dynamo. The onset of this instability occurs after a rather long temporal evolution of several thousand turbulent turnover times. The vorticity dynamo in turn drives an exponential growth of the magnetic field, first at small scales, then also at large one. The instability then saturates and the magnetic field approximately reaches equipartition with the turbulent kinetic energy. During the saturation phase we can observe a winding of the magnetic field in the direction of the shearing flow. By varying the intensity of the shear we see that the growth rates of this instability change. The inclusion of the baroclinic term delays the onset of the vorticity dynamo but leads to a more rapid growth. We demonstrate how in the presence of shear, even a purely irrotational forcing amplifies the field to equipartition. At the same time, we confirm how this forcing alone does not lead to vorticity nor magnetic field growth, and this picture does not change in the presence of rotation or baroclinicity up to $256^3$ meshpoints.

After the prediction of meteor showers in the Earth's atmosphere caused by the particles originating in the nucleus of comet 21P/Giacobini-Zinner, we went on with the prediction of showers on the other three terrestrial planets. Based on our modeling of theoretical stream of the parent comet, we predicted several related meteorite (on Mercury) or meteor (on Venus and Mars) showers. There occurred the filaments, in the stream, with the particles coming to each planet from a similar direction. We found that this is a consequence of the specific distribution of argument of perihelion (peaked close to the value of $180^{\circ}$) and longitude of ascending node of the stream, and that the particles collide with each planet in an arc of their orbits being close to perihelion.

We present a timing and noise analysis of the Be/X-ray binary system Swift J0243.6+6124 during its 2017-2018 super-Eddington outburst using NICER/XTI observations. For the initial segments of the data that overlap with the Fermi/GBM pulse frequency history, we apply a synthetic pulse timing analysis to enrich the spin frequency history of the source. In addition, we employ phase-coherent timing analysis for NICER/XTI observations that extends beyond the Fermi/GBM frequency history. We show that the pulse profiles switch from double-peaked to single-peaked when the X-ray luminosity drops below $\sim$$7\times 10^{36}$ erg s$^{-1}$. We suggest that this transitional luminosity is associated with the transition from a pencil beam pattern to a hybrid beam pattern when the Coulomb interactions become ineffective to decelerate the accretion flow, which implies a dipolar magnetic field strength of $\sim$$5\times 10^{12}$ G. We also obtained the power density spectra (PDS) of the spin frequency derivative fluctuations. The red noise component of the PDS is found to be steeper ($\omega^{-3.36}$) than the other transient accreting sources. We find significantly high noise strength estimates above the super-Eddington luminosity levels, which may arise from the torque fluctuations due to interactions with the quadrupole fields at such levels.

We present a new homogeneous survey of VLT/FLAMES LR8 line-of-sight radial velocities (vlos) for 1604 resolved red giant branch stars in the Sculptor dwarf spheroidal galaxy. In addition, we provide reliable Ca II triplet metallicities, [Fe/H], for 1339 of these stars. From this combination of new observations (2257 individual spectra) with ESO archival data (2389 spectra), we obtain the largest and most complete sample of vlos and [Fe/H] measurements for individual stars in any dwarf galaxy. Our sample includes VLT/FLAMES LR8 spectra for 55% of the red giant branch stars at G $<20$ from Gaia DR3, and $>70$% of the brightest stars, G $<18.75$. Our spectroscopic velocities are combined with Gaia DR3 proper motions and parallax measurements for a new and more precise membership analysis. We look again at the global characteristics of Sculptor, deriving a mean metallicity of $\langle$[Fe/H]$\rangle = -1.82 \pm 0.45$ and a mean line-of-sight velocity of $\langle$vlos$\rangle = +111.2 \pm 0.25$km/s. There is a clear metallicity gradient in Sculptor, -0.7deg/dex, with the most metal-rich population being the most centrally concentrated. Furthermore, the most metal-poor population in Sculptor, [Fe/H]$<-2.5$, appears to show kinematic properties distinct from the rest of the stellar population. Finally, we combine our results with the exquisite Gaia DR3 multi-colour photometry to further investigate the colour-magnitude diagram of the resolved stellar population in Sculptor. Our detailed analysis shows a similar global picture as previous studies, but with much more precise detail, revealing that Sculptor has more complex properties than previously thought. This survey emphasises the role of the stellar spectroscopy technique and this galaxy as a benchmark system for modelling galaxy formation and evolution on small scales.

In this paper we introduce the Terzina telescope as a part of the NUSES space mission. This telescope aims to detect Ultra High Energy Cosmic Rays (UHECRs) through the Cherenkov light emission from the extensive air showers (EAS) that they create in the Earth's atmosphere. The Cherenkov photons are aligned along the shower axis inside about $\sim 0.2-1^{\circ}$, so that they become detectable by Terzina when it points towards the Earth's limb. A sun-synchronous orbit will allow the telescope to observe only the night side of the Earth's atmosphere. In this contribution, we focus on the description of the telescope detection goals, geometry, optical design and its photon detection camera composed of Silicon Photo-Multipliers (SiPMs). Moreover, we describe the full Monte Carlo simulation chain developed to estimate Terzina's performance for UHECR detection. The estimate of the radiation damage and light background rates, the readout electronics and trigger logic are briefly described. Terzina will be able to study the potential for future physics missions devoted to UHECR detection and to UHE neutrino astronomy. It is a pathfinder for missions like POEMMA or future constellations of similar satellites to NUSES.

Fast variability of the X-ray corona in black hole binaries can produce a soft lag by reverberation, where the reprocessed thermalized disc photons lag behind the illuminating hard X-rays. This lag is small, and systematically decreases with increasing mass accretion rate towards the hard-soft transition, consistent with a decreasing truncation radius between the thin disc and X-ray hot inner flow. However, the soft lag suddenly increases dramatically just before the spectrum becomes disc-dominated (hard-intermediate state). Interpreting this as reverberation requires that the X-ray source distance from the disc increases dramatically, potentially consistent with switching to X-rays produced in the radio jet. However, this change in lag behaviour occurs without any clear change in hard X-ray spectrum, and before the plasmoid ejection event which might produce such a source (soft-intermediate state). Instead, we show how the soft lag can be interpreted in the context of propagation lags from mass accretion rate fluctuations. These normally produce hard lags, as the model has radial stratification, with fluctuations from larger radii modulating the harder spectra produced at smaller radii. However, all that is required to switch the sign is that the hottest Comptonized emission has seed photons which allow it to extend down in energy below the softer emission from the slower variable turbulent region from the inner edge of the disc. Our model connects the timing change to the spectral change, and gives a smooth transition of the X-ray source properties from the bright hard state to the disc-dominated states.

We generalize previously derived analytic results for the one-loop power spectrum (PS) in scale-free models (with linear PS $P(k) \propto k^n$) to a broader class of such models in which part of the matter-like component driving the Einstein de Sitter expansion does not cluster. These models can be conveniently parametrized by $\alpha$, the constant logarithmic linear growth rate of fluctuations (with $\alpha=1$ in the usual case). For $-3< n<-1$, where the one-loop PS is both infrared and ultraviolet convergent and thus explicitly self-similar, it is characterized conveniently by a single numerical coefficient $c(n, \alpha)$. We compare the analytical predictions for $c(n=-2, \alpha)$ with results from a suite of $N$-body simulations with $\alpha \in [0.25, 1]$ performed with an appropriately modified version of the Gadget code. Although the simulations are of small ($256^3$) boxes, the constraint of self-similarity allows the identification of the converged PS at a level of accuracy sufficient to test the analytical predictions for the $\alpha$ dependence of the evolved PS. Good agreement for the predicted dependence on $\alpha$ of the PS is found. To treat the UV sensitivity of results which grows as one approaches $n =-1$, we derive exact results incorporating a regularisation $k_c$ and obtain expressions for $c(n, \alpha, k_c/k)$. Assuming that this regularisation is compatible with self-similarity allows us to infer a predicted functional form of the PS equivalent to that derived in effective field theory (EFT). The coefficient of the leading EFT correction at one loop has a strong dependence on $\alpha$, with a change in sign at $\alpha \approx 0.16$, providing a potentially stringent test of EFT.

Circumbinary discs around evolved post-asymptotic giant branch (post-AGB) binary systems show many similar properties to protoplanetary discs around young stars. Deficits of near-infrared (near-IR) flux in the spectral energy distributions (SEDs) of such systems hints towards large dust-free cavities, reminiscent of transition discs as commonly observed around young stars. We aim to assess the inner rim size of 6 post-AGB binary systems with such a lack in near-IR using resolved mid-IR high-angular resolution observations of VLTI/MATISSE and VLTI/MIDI. The inner rim of only one such system was previously resolved. We compare these inner rim sizes to 5 systems with available MATISSE data that were identified to host a disc starting at the dust sublimation radius. We used geometric ring models to estimate the inner rim sizes, the relative flux contributions of the star, the ring, and an over-resolved emission, the orientation of the ring, and the spectral dependencies of the components. We find that the dust inner rims of the targets with a lack of near-IR excess in their SEDs are 2.5 to 7.5 times larger than the theoretical dust sublimation radii while the systems that do not show such a deficit have inner rim sizes similar to their dust sublimation radii. Physical radii of the inner rims of these transition discs around post-AGB binaries are 3-25 au, which are larger than the disc sizes inferred for transition discs around young stars with VLTI/MIDI. This is due to the higher stellar luminosities of post-AGB systems compared to young stars, implying larger dust sublimation radii and thus larger physical transition disc inner radii. With mid-IR interferometric data we directly confirm the transition disc nature of six discs around post-AGB binary systems. Future observational and modelling efforts are needed to progress on the structure, origin, and evolution of these transition discs.

Small flares frequently occur in the quiet Sun. Previous studies have noted that they share many common characteristics with typical solar flares in active regions. However, their similarities and differences are not fully understood, especially their thermal properties. In this study, we performed imaging spectroscopic observations in the H$\alpha$ line taken with the Solar Dynamics Doppler Imager on the Solar Magnetic Activity Research Telescope (SMART/SDDI) at the Hida Observatory and imaging observations with the Atmospheric Imaging Assembly onboard Solar Dynamics Observatory (SDO/AIA). We analysed 25 cases of small flares in the quiet Sun over the thermal energy range of $10^{24}-10^{27}\,\mathrm{erg}$, paying particular attention to their thermal properties. Our main results are as follows: (1) We observe a redshift together with line centre brightening in the H$\alpha$ line associated with more than half of the small flares. (2) We employ differential emission measure analysis using AIA multi-temperature (channel) observations to obtain the emission measure and temperature of the small flares. The results are consistent with the Shibata & Yokoyama (1999, 2002) scaling law. From the scaling law, we estimated the coronal magnetic field strength of small flares to be 5 --15 G. (3) The temporal evolution of the temperature and the density shows that the temperature peaks precede the density peaks in more than half of the events. These results suggest that chromospheric evaporations/condensations play an essential role in the thermal properties of some of the small flares in the quiet Sun, as does for large flares.

We prepared a sample of oxirane doubly deuterated at one C atom and studied its rotational spectrum in the laboratory for the first time between 120~GHz and 1094~GHz. Accurate spectroscopic parameters up to eighth order were determined, and the calculated rest frequencies were used to identify $c$-CD$_2$CH$_2$O tentatively in the interstellar medium in the Atacama Large Millimeter/submillimeter Array Protostellar Interferometric Line Survey (PILS) of the Class 0 protostellar system IRAS 16293$-$2422. The $c$-CD$_2$CH$_2$O to $c$-C$_2$H$_4$O ratio was estimated to be $\sim$0.054 with $T_{\rm rot} = 125$ K. This value translates to a D-to-H ratio of $\sim$0.16 per H atom which is higher by a factor of 4.5 than the $\sim$0.036 per H atom obtained for $c$-C$_2$H$_3$DO. Such increase in the degree of deuteration referenced to one H atom in multiply deuterated isotopologs compared to their singly deuterated variants have been observed commonly in recent years.

Rocky planets and moons experiencing solar wind sputtering are continuously supplying their enveloping exosphere with ejected neutral atoms. To understand the quantity and properties of the ejecta, well established Binary Collision Approximation Monte Carlo codes like TRIM with default settings are used predominantly. Improved models such as SDTrimSP have come forward and together with new experimental data the underlying assumptions have been challenged. We introduce a hybrid model, combining the previous surface binding approach with a new bulk binding model akin to Hofs\"ass & Stegmaier (2023). In addition, we expand the model implementation by distinguishing between free and bound components sourced from mineral compounds such as oxides or sulfides. The use of oxides and sulfides also enables the correct setting of the mass densities of minerals, which was previously limited to the manual setting of individual atomic densities of elements. All of the energies and densities used are thereby based on tabulated data, so that only minimal user input and no fitting of parameters are required. We found unprecedented agreement between the newly implemented hybrid model and previously published sputter yields for incidence angles up to 45{\deg} from surface normal. Good agreement is found for the angular distribution of mass sputtered from enstatite MgSiO$_3$ compared to latest experimental data. Energy distributions recreate trends of experimental data of oxidized metals. Similar trends are to be expected from future mineral experimental data. The model thus serves its purpose of widespread applicability and ease of use for modelers of rocky body exospheres.

We investigate the role of the close environment on the nuclear activity of a sample of 436 nearby ($z<0.3$) QSO 2's -- selected from SDSS-III spectra, via comparison of their environment and interaction parameters with those of a control sample of 1308 galaxies. We have used the corresponding SDSS images to obtain the number of neighbour galaxies $N$, tidal strength parameter $Q$ and asymmetry parameters. We find a small excess of $N$ in the QSOs compared to its three controls, and no difference in $Q$. The main difference is an excess of asymmetry in the QSOs hosts, which is almost twice that of the control galaxies. This difference is not due to the hosts' morphology, since there is no difference in their Galaxy Zoo classifications. HST images of two highly asymmetric QSO 2 hosts of our sample show that both sources have a close companion (at projected separations $\sim$ 5 kpc), which we thus conclude is the cause of the observed asymmetry in the lower resolution SDSS images. The mean projected radius of the controls is $ \langle r \rangle = 8.53\pm$0.06 kpc, while that of the QSO hosts is $ \langle r \rangle = 9.39\pm$0.12 kpc, supporting the presence of interaction signatures in the outer regions of the QSO hosts. Our results favour a scenario in which nuclear activity in QSO 2's is triggered by close galaxy interactions -- when the distance between the host and companion is of the order of the galaxy radius, implying that they are already in the process of merger.

Imaging in radioastronomy is an ill-posed inverse problem. Particularly the Event Horizon Telescope (EHT) Collaboration faces two big limitations for the existing methods when imaging the active galactic nuclei (AGN): large and expensive surveys solving the problem with different optimization parameters must be done, and only one local minima for each instance is returned. With our novel nonconvex, multiobjective optimization modeling approach, we aim to overcome these limitations. To this end we used a multiobjective version of the genetic algorithm (GA): the Multiobjective Evolutionary Algorithm Based on Decomposition, or MOEA/D. GA strategies explore the objective function by evolutionary operations to find the different local minima, and to avoid getting trapped in saddle points. First, we have tested our algorithm (MOEA/D) using synthetic data based on the 2017 Event Horizon Telescope (EHT) array and a possible EHT + next-generation EHT (ngEHT) configuration. We successfully recover a fully evolved Pareto front of non-dominated solutions for these examples. The Pareto front divides into clusters of image morphologies representing the full set of locally optimal solutions. We discuss approaches to find the most natural guess among these solutions and demonstrate its performance on synthetic data. Finally, we apply MOEA/D to observations of the black hole shadow in Messier 87 (M87) with the EHT data in 2017. MOEA/D is very flexible, faster than any other Bayesian method and explores more solutions than Regularized Maximum Likelihood methods (RML). We have done two papers to present this new algorithm: the first explains the basic idea behind multi-objective optimization and MOEA/D and it is used to recover static images, while in the second paper we extend the algorithm to allow dynamic and (static and dynamic) polarimetric reconstructions.

We studied the properties of the young stellar populations in the NGC 299 cluster in the Small Magellanic Cloud using observations obtained with the Hubble Space Telescope in the $V, I$, and $H\alpha$ bands. We identified 252 stars with H$\alpha$ excess exceeding 5 $\sigma$ and an equivalent width of the H$\alpha$ emission line of at least 20 \r{A}, indicating that these stars are still undergoing accretion and therefore represent bona fide pre-main-sequence (PMS) objects. For all of them, we derived the mass, age, and mass accretion rate by comparing the observed photometry with theoretical models. We find evidence for the existence of two populations of PMS stars, with median ages of 25 and 50 Myr respectively. The average mass accretion rate for these PMS stars is $\sim 5 \times 10^{-9}$ M$_\odot$ yr$^{-1}$, which is comparable to the values found in other low-metallicity, low-density clusters in the Magellanic Clouds, but is about a factor of three lower than those measured for stars of similar mass and age in denser Magellanic Cloud stellar regions. Our findings support the hypothesis that both the metallicity and density of the forming environment can affect the mass accretion rate and thus the star formation process in a region. A study of the spatial distribution of both massive stars and (low-mass) PMS objects reveals that the former are clustered near the nominal centre of NGC 299, whereas the PMS stars are rather uniformly distributed over the field. To explore whether the stars formed in an initially more diffuse or compact structure, we studied the cluster's stellar density profile. We find a core radius $r_c\simeq 0.6$ pc and a tidal radius $r_t\simeq 5.5$ pc, with an implied concentration parameter $c \simeq 1$, suggesting that the cluster could be dispersing into the field.

The non-association of coronal mass ejections with high energetic flares is sparse. For this reason, the magnetic conditions required for the confinedness of major flares is a topic of active research. Using multi-instrument observations, we investigated the evolution and effects of confinedness in an X3.1 flare, which occurred in active region (AR) 12192. The decrease of net fluxes in the brightening regions, near the footpoints of the multi-sigmoidal AR in photosphere and chromosphere, indicative of flux cancellation favouring tether-cutting reconnection (TCR), is observed using the magnetic field observations of HMI/SDO and SOT/Hinode, respectively. The analysis of spectropolarimetric data obtained by the Interferometric Bidimensional Spectrometer over the brightening regions suggests untwisting of field lines, which further supports TCR. Filaments near polarity inversion line region, resulted from TCR of low lying sheared loops, undergo merging and form an elongated filament. The temperature and density differences between footpoints of the merged filament, revealed by DEM analysis, caused streaming and counter-streaming of plasma flow along the filament and unloads at its footpoints with an average velocity of $\approx$ 40 km s$^{-1}$. This results in decrease of mass of the filament (density decreased by $>50\%$), leading to its rise and expansion outwards. However, due to strong strapping flux, the filament separates itself instead of erupting. Further, the evolution of non-potential parameters describes the characteristics of confinedness of the flare. Our study suggests that the sigmoid-filament system exhibits upward catastrophe due to mass unloading, but gets suppressed by strong confinement of external poloidal field.

In this paper, we present observations of two high-resolution transit datasets obtained with ESPRESSO of the bloated sub-Saturn planet WASP-131~b. We have simultaneous photometric observations with NGTS and EulerCam. In addition, we utilised photometric lightcurves from {\tess}, WASP, EulerCam and TRAPPIST of multiple transits to fit for the planetary parameters and update the ephemeris. We spatially resolve the stellar surface of WASP-131 utilising the Reloaded Rossiter McLaughlin technique to search for centre-to-limb convective variations, stellar differential rotation, and to determine the star-planet obliquity for the first time. We find WASP-131 is misaligned on a nearly retrograde orbit with a projected obliquity of $\lambda = 162.4\substack{+1.3 \\ -1.2}^{\circ}$. In addition, we determined a stellar differential rotation shear of $\alpha = 0.61 \pm 0.06$ and disentangled the stellar inclination ($i_* = 40.9\substack{+13.3 \\ -8.5}^{\circ}$) from the projected rotational velocity, resulting in an equatorial velocity of $v_{\rm{eq}} = 7.7\substack{+1.5 \\ -1.3}$~km s$^{-1}$. In turn, we determined the true 3D obliquity of $\psi = 123.7\substack{+12.8 \\ -8.0}^{\circ}$, meaning the planet is on a perpendicular/polar orbit. Therefore, we explored possible mechanisms for the planetary system's formation and evolution. Finally, we searched for centre-to-limb convective variations where there was a null detection, indicating that centre-to-limb convective variations are not prominent in this star or are hidden within red noise.

We consider a scenario in which Sgr A* is in a massive black hole binary (MBHB) with an as-of-yet undetected supermassive or intermediate-mass black hole companion. Dynamical encounters between this MBHB and single stars in its immediate vicinity would eject hypervelocity stars (HVSs) with velocities beyond the Galactic escape velocity. In this work, we use existing HVS observations to constrain for the first time the existence of a companion to Sgr A*. We simulate the ejection of HVSs via the `MBHB slingshot' scenario and show that the population of HVSs detectable today depends strongly on the companion mass and the separation of the MBHB. We demonstrate that the lack of uncontroversial HVS candidates in \textit{Gaia} Data Release 3 places a firm upper limit on the mass of a possible Sgr A* companion. Within one milliparsec of Sgr A*, our results exclude a companion more massive than $1000 \, \mathrm{M_\odot}$. If Sgr A* recently merged with a companion black hole, our findings indicate that unless this companion was less massive than $500 \, \mathrm{M_\odot}$, this merger must have occurred at least $10$ Myr ago. These results complement and improve upon existing independent constraints on a companion to Sgr A* and show that large regions of its parameter space can now be ruled out.

We present an empirical model of age-dependent photospheric lithium depletion, calibrated using a large, homogeneously-analysed sample of 6200 stars in 52 open clusters, with ages from 2--6000 Myr and $-0.3<{\rm [Fe/H}]<0.2$, observed in the Gaia-ESO spectroscopic survey. The model is used to obtain age estimates and posterior age probability distributions from measurements of the Li I 6708A equivalent width for individual (pre) main sequence stars with $3000 < T_{\rm eff}/{\rm K} <6500$, a domain where age determination from the HR diagram is either insensitive or highly model-dependent. In the best cases, precisions of 0.1 dex in log age are achievable; even higher precision can be obtained for coeval groups and associations where the individual age probabilities of their members can be combined. The method is validated on a sample of exoplanet-hosting young stars, finding agreement with claimed young ages for some, but not others. We obtain better than 10 per cent precision in age, and excellent agreement with published ages, for seven well-studied young moving groups. The derived ages for young clusters ($<1$ Gyr) in our sample are also in good agreement with their training ages, and consistent with several published, model-insensitive lithium depletion boundary ages. For older clusters there remain systematic age errors that could be as large as a factor of two. There is no evidence to link these errors to any strong systematic metallicity dependence of (pre) main sequence lithium depletion, at least in the range $-0.29 < {\rm [Fe/H]} < 0.18$. Our methods and model are provided as software -- "Empirical AGes from Lithium Equivalent widthS" (EAGLES).

The central compact object within HESS J1731- 347 possesses unique mass and radius properties that renders it a compelling candidate for a self-bound star. In this research, we examine the capability of quark stars composed of colour superconducting quark matter to explain the latter object by using its marginalised posterior distribution and imposing it as a constraint on the relevant parameter space. Namely, we investigate quark matter for $N_f=2,3$ in the colour superconducting phase, incorporating perturbative QCD corrections, and we derive their properties accordingly. The utilised thermodynamic potential of this work possesses an MIT bag model formalism with the parameters being established as flavour-independent. In this instance, we conclude the favour of 3-flavour over 2-flavour colour superconducting quark matter, isolating our interest on the former. The parameter space is further confined due to the additional requirement for a high maximum mass ($M_{\text{TOV}} \geq 2.6 M_{\odot}$), accounting for GW$190814$'s secondary companion. We pay a significant attention on the speed of sound and the trace anomaly (proposed as a measure of conformality [\href{https://journals.aps.org/prl/abstract/10.1103/PhysRevLett.129.252702}{10.1103/PhysRevLett.129.252702}]). We conclude that it is possible for colour-flavour locked quark stars to reach high masses without violating the conformal bound or the $\langle \Theta \rangle _{\mu_B} \geq 0$ if the quartic coefficient value $\alpha_4$ does not exceed an upper limit which is solely dependent on the established $M_{\text{TOV}}$. For $M_{\text{TOV}}=2.6 M_{\odot}$, we find that the limit reads $\alpha_4 \leq 0.594$. Lastly, a further study takes place on the agreement of colour-flavour locked quark stars with additional astrophysical objects including the GW$170817$ and GW$190425$ events, followed by a relevant discussion.

Binary black holes (BHs) can be formed dynamically in the centers of star clusters. The high natal kicks for stellar-mass BHs used in previous works made it hard to retain BHs in star clusters. Recent studies of massive star evolution and supernovae (SN) propose kick velocities that are lower due to the fallback of the SN ejecta. We study the impact of these updates by performing $N$-body simulations following instantaneous gas expulsion. For comparison, we simulate two additional model sets with the previous treatment of stars: one with high kicks and another with artificial removal of the kicks. Our model clusters initially consist of about one hundred thousand stars, formed with centrally-peaked efficiency. We find that the updated treatment of stars, due to the fallback-scaled lower natal kicks, allows clusters to retain SN remnants after violent relaxation. The mass contribution of the retained remnants does not exceed a few percent of the total bound cluster mass during the early evolution. For this reason, the first giga year of evolution is not affected significantly by this effect. Nevertheless, during the subsequent long-term evolution, the retained BHs accelerate mass segregation, leading to the faster dissolution of the clusters.

The occurrence of a first-order hadron-quark matter phase transition at high baryon densities is investigated in astrophysical simulations of core-collapse supernovae, to decipher yet incompletely understood properties of the dense matter equation of state using neutrinos from such cosmic events. It is found that the emission of a non-standard second neutrino burst, dominated by electron-antineutrinos, is not only a measurable signal for the appearance of deconfined quark matter but also reveals information about the state of matter at extreme conditions encountered at the supernova interior. To this end, a large set of spherically symmetric supernova models is investigated, studying the dependence on the equation of state and on the stellar progenitor. General relativistic neutrino-radiation hydrodynamics is employed featuring three-flavor Boltzmann neutrino transport and a microscopic hadron-quark hybrid matter equation of state class, that covers a representative range of parameters. This facilitates the direct connection between intrinsic signatures of the neutrino signal and properties of the equation of state. In particular, a set of novel relations have been found empirically. These potentially provide a constraint for the onset density of a possible QCD phase transition, which is presently one of the largest uncertainties in modern investigations of the QCD phase diagram, from the future neutrino observation of the next galactic core-collapse supernova.

Axions and axion-like particles (ALPs) are a prominent dark matter candidate, drawing motivation in part from the axiverse of string theory. Axion-like particles can also arise as composite degrees of freedom of a dark sector, for example, as dark pions in dark Quantum Chromo-Dynamics. In a dark Standard Model (SM) wherein all 6 quark flavors are light while the photon is massive, one finds a rich low-energy spectrum of stable and ultralight particles, in the form of neutral and charged dark scalars, and complex neutral scalars analogous to the SM kaon, with mass splittings determined by the mass and charge of the dark quarks. The model finds a natural portal to the visible sector via kinetic coupling of the dark and visible photons, and consequent millicharges for dark matter. The dark matter can be a mixture of all these ultralight bosonic degrees of freedom, and exhibit both parity-even and parity-odd interactions, making the theory testable at a wide variety of experiments. In context of dark QCD with $N_f$ flavors of light quarks, this scenario predicts $N_f^2-1$ ultralight axion-like particles -- effectively an axiverse from dark QCD. This '$\pi$-axiverse' is consistent with but makes no recourse to string theory, and is complementary to the conventional string theory axiverse.

The inspiral phase of gravitational waves emitted by spinless compact binary systems is derived through the fourth-and-a-half post-Newtonian (4.5PN) order beyond quadrupole radiation, and the leading amplitude mode ($\ell$, m) = (2, 2) is obtained at 4PN order. We also provide the radiated flux, as well as the phase in the stationary phase approximation. Rough numerical estimates for the contribution of each PN order are provided for typical systems observed by current and future gravitational wave detectors.

In this work, we report an autoencoder-based 2D representation to classify a time-series as stochastic or non-stochastic, to understand the underlying physical process. Content-aware conversion of 1D time-series to 2D representation, that simultaneously utilizes time- and frequency-domain characteristics, is proposed. An autoencoder is trained with a loss function to learn latent space (using both time- and frequency domains) representation, that is designed to be, time-invariant. Every element of the time-series is represented as a tuple with two components, one each, from latent space representation in time- and frequency-domains, forming a binary image. In this binary image, those tuples that represent the points in the time-series, together form the ``Latent Space Signature" (LSS) of the input time-series. The obtained binary LSS images are fed to a classification network. The EfficientNetv2-S classifier is trained using 421 synthetic time-series, with fair representation from both categories. The proposed methodology is evaluated on publicly available astronomical data which are 12 distinct temporal classes of time-series pertaining to the black hole GRS 1915 + 105, obtained from RXTE satellite. Results obtained using the proposed methodology are compared with existing techniques. Concurrence in labels obtained across the classes, illustrates the efficacy of the proposed 2D representation using the latent space co-ordinates. The proposed methodology also outputs the confidence in the classification label.

Over the last few years, there has been an increasing interest in sub-solar mass black holes due to their potential to provide valuable information about cosmology or the black hole population. Motivated by this, we study observable phenomena connected to the merger of a sub-solar mass black hole with a neutron star. For this purpose, we perform new numerical-relativity simulations of a binary system composed of a black hole with mass $0.5M_\odot$ and a neutron star with mass $1.4 M_\odot$. We investigate the merger dynamics of this exotic system and provide information about the connected gravitational-wave and kilonova signals. Our study indicates that current gravitational-waveform models are unable to adequately describe such systems and that phenomenological relations connecting the binary parameters with the ejecta and remnant properties are not applicable to our system. Furthermore, we find a dependence of the kilonova signal on the azimuthal viewing angle due to the asymmetric mass ejection. This first-of-its-kind simulation opens the door for the study of sub-solar mass black hole - neutron star mergers and could serve as a testing ground for future model development.

We study the effect of an ultra-light primordial black hole (PBH) dominated phase on the gravitational wave (GW) spectrum generated by a cosmic string (CS) network formed as a result of a high-scale $U(1)$ symmetry breaking. A PBH-dominated phase leads to tilts in the spectrum via entropy dilution and generates a new GW spectrum from PBH density fluctuations, detectable at ongoing and planned near-future GW detectors. The combined spectrum has a unique shape with a plateau, a sharp tilted peak over the plateau, and a characteristic fall-off, which can be distinguished from the one generated in the combination of CS and any other matter domination or new exotic physics. We discuss how ongoing and planned future experiments can probe such a unique spectrum for different values of $U(1)$ breaking scale and PBH parameters such as initial mass and energy fraction.

Americium-beryllium (AmBe), a well-known tagged neutron source, is commonly used for evaluating the neutron detection efficiency of detectors used in ultralow background particle physics experiments, such as reactor neutrino and diffuse supernova neutrino background experiments. In particular, AmBe sources are used to calibrate neutron tagging by selecting the 4438-keV $\gamma$-ray signal, which is simultaneously emitted with a neutron signal. Therefore, analyzing the neutron and $\gamma$-ray emission properties of AmBe sources is crucial. In this study, we used the theoretical shape of a neutron energy spectrum, which was divided into three parts, to develop models of the energy spectrum and verify the results using experimental data. We used an AmBe source to measure the energy spectra of simultaneously emitted neutrons and $\gamma$-rays and determine the emission ratio of the neutrons with and without $\gamma$-ray emission. The measured spectrum was consistent with that obtained from the simulated result, whereas the measured emission ratio was significantly different from the corresponding simulated result. Here, we also discuss the feasibility of determining the neutron emission rates from the spectra divided into three parts.