Papers for Tuesday, Jan 30 2024

Cosmological observational programs often compare their data not only with $\Lambda$CDM, but also with dark energy (DE) models, whose time-dependent equations of state (EoS) differ from a cosmological constant. We identified a generic issue in the standard procedure of computing the expansion history for models with time-dependent EoS, which leads to bias in the interpretation of the results. In order to compute the evolution of models with time-dependent EoS parameter $w$ in a consistent manner, we introduce an enhanced computational procedure, which accounts for the correct choice of initial conditions in the respective backward-in-time and forward-in-time evolution of the equations of motion. We implement our enhanced procedure in an amended version of the code CLASS, where we focus on exemplary DE models which are based on the CPL parameterization, studying cases with monotonically increasing and decreasing $w$ over cosmic time. Our results reveal that a cosmological DE model with a decreasing $w$ of the form $w(a)=-0.9 + 0.1(1-a)$ could provide a resolution to the Hubble tension problem. Moreover, we find characteristic signatures in the late expansion histories of models, allowing a phenomenological discrimination of DE candidates. Finally, we argue that our enhanced scheme should be implemented as a novel consistency check for cosmological models within current Monte-Carlo-Markov-Chain (MCMC) methods. Our enhanced computational procedure avoids the interpretational bias to which the standard procedure is unwittingly exposed. As a result, DE models can be better constrained. If implemented into MCMC codes, we expect that our scheme will contribute to providing a significant improvement in the determined accuracy of cosmological model parameters.

The exquisite angular resolution and sensitivity of JWST is opening a new window for our understanding of the Universe. In nearby galaxies, JWST observations are revolutionizing our understanding of the first phases of star formation and the dusty interstellar medium. Nineteen local galaxies spanning a range of properties and morphologies across the star-forming main sequence have been observed as part of the PHANGS-JWST Cycle 1 Treasury program at spatial scales of $\sim$5-50pc. Here, we describe pjpipe, an image processing pipeline developed for the PHANGS-JWST program that wraps around and extends the official JWST pipeline. We release this pipeline to the community as it contains a number of tools generally useful for JWST NIRCam and MIRI observations. Particularly for extended sources, pjpipe products provide significant improvements over mosaics from the MAST archive in terms of removing instrumental noise in NIRCam data, background flux matching, and calibration of relative and absolute astrometry. We show that slightly smoothing F2100W MIRI data to 0.9" (degrading the resolution by about 30 percent) reduces the noise by a factor of $\approx$3. We also present the first public release (DR1.0.1) of the pjpipe processed eight-band 2-21 $\mu$m imaging for all nineteen galaxies in the PHANGS-JWST Cycle 1 Treasury program. An additional 55 galaxies will soon follow from a new PHANGS-JWST Cycle 2 Treasury program.

Cosmological probes constructed in large-scale surveys are independent of the underlying theory of gravity, and their relativistic descriptions are indeed applicable to any theory of gravity. It was shown that the presence of fluctuations with wavelength much larger than the characteristic scales of the surveys has no impact on cosmological probes, if the matter content is adiabatic and the Einstein equations are used. In this paper we study the sensitivity of cosmological probes to infrared fluctuations in Horndeski theory. We find that the extra degree of freedom in the Horndeski scalar field can induce sensitivity to infrared fluctuations in the cosmological probes, even when the matter components are adiabatic on large scales. A generalized adiabatic condition including the extra dof, in contrast, guarantees that cosmological probes are devoid of infrared sensitivity, and this solution corresponds to the adiabatic modes \`a la Weinberg in Horndeski theory, which can be removed by a coordinate transformation in the infrared limit. We discuss the implications of our findings and the connections to the initial conditions.

Obtaining spectroscopic observations of the progenitors of core-collapse supernovae is often unfeasible due to an inherent lack of knowledge as to which stars will go supernova and when they will explode. In this letter, we present photometric and spectroscopic observations of the progenitor activity of SN 2023fyq in the preceding 150 days before the He-rich progenitor exploded as a Type Ibn supernova. The progenitor of SN 2023fyq shows an exponential rise in flux prior to core-collapse. Complex He I emission line features are observed, with a P-Cygni like profile, as well as an evolving broad base with velocities on the order of 10,000 km/s, possibly due to electron scattering. The luminosity and evolution of SN 2023fyq are consistent with a faint Type Ibn, reaching a peak r-band magnitude of 18.1 mag, although there is some uncertainty in the distance to the host, NGC 4388, located in the Virgo cluster. We present additional evidence of asymmetric He-rich material being present prior to the explosion of SN 2023fyq, as well as after, suggesting this material has survived the ejecta-CSM interaction. Broad [O I] and the Ca II triplet lines are observed at late phases, confirming that SN 2023fyq was a genuine supernova rather than a non-terminal interacting transient. SN 2023fyq provides insight into the final moments of a massive star's life, highlighting that the progenitor is likely highly unstable before core-collapse.

$\omega$ Centauri is considered the most massive globular cluster of the Milky Way and likely the former nuclear star cluster of a galaxy accreted by the Milky Way. It is speculated to contain an intermediate-mass black hole (IMBH) from several dynamical models. However, uncertainties regarding the location of the cluster center or the retention of stellar remnants limit the robustness of the IMBH detections reported so far. In this paper, we derive and study the stellar kinematics from the highest-resolution spectroscopic data yet, using the Multi Unit Spectroscopic Explorer (MUSE) in the narrow field mode (NFM) and wide field mode (WFM). Our exceptional data near the center reveal for the first time that stars within the inner 20" ($\sim$0.5 pc) counter-rotate relative to the bulk rotation of the cluster. Using this dataset, we measure the rotation and line-of-sight velocity dispersion (LOSVD) profile out to 120$''$ with different centers proposed in the literature. We find that the velocity dispersion profiles using different centers match well with those previously published. Based on the counter--rotation, we determine a kinematic center and look for any signs of an IMBH using the high-velocity stars close to the center. We do not find any significant outliers $>$60 km/s within the central 20$''$, consistent with no IMBH being present at the center of $\omega$ Centauri. A detailed analysis of Jeans' modeling of the putative IMBH will be presented in the next paper of the series.

The finding of active galactic nuclei (AGN) in dwarf galaxies has important implications for galaxy evolution and supermassive black hole formation models. Yet, how AGN in dwarf galaxies form is still debated, in part due to scant demographics. We make use of the MaNGA survey, comprising $\sim$10,000 galaxies at z $<$ 0.15, to identify AGN dwarf galaxies using a spaxel by spaxel classification in three spatially-resolved emission line diagnostic diagrams (the [NII-, [SII]- and [OI]-BPT) and the WHAN diagram. This yields a sample of 664 AGN dwarf galaxies, the largest to date, and an AGN fraction of $\sim20\%$ that is significantly larger than that of single-fiber-spectroscopy studies (i.e. $\sim1\%$). This can be explained by the lower bolometric luminosity ($< 10^{42}$ erg s$^{-1}$) and accretion rate (sub-Eddington) of the MaNGA AGN dwarf galaxies. We additionally identify 1,176 SF-AGN (classified as star-forming in the [NII]-BPT but as AGN in the [SII]- and [OI]-BPT), 122 Composite, and 173 LINER sources. The offset between the optical center of the galaxy and the median position of the AGN spaxels is more than 3 arcsec for $\sim$62\% of the AGN, suggesting that some could be off-nuclear. We also identify seven new broad-line AGN with log $M_\mathrm{BH}$ = 5.0 - 5.9 $M_\mathrm{\odot}$. Our results show how integral-field spectroscopy is a powerful tool for uncovering faint and low-accretion AGN and better constraining the demographics of AGN in dwarf galaxies.

Isotopologue ratios are anticipated to be one of the most promising signs of life that can be observed remotely. On Earth, carbon isotopes have been used for decades as evidence of modern and early metabolic processes. In fact, carbon isotopes may be the oldest evidence for life on Earth, though there are alternative geological processes that can lead to the same magnitude of fractionation. However, using isotopologues as biosignature gases in exoplanet atmospheres presents several challenges. Most significantly, we will only have limited knowledge of the underlying abiotic carbon reservoir of an exoplanet. Atmospheric carbon isotope ratios will thus have to be compared against the local interstellar medium or, better yet, their host star. A further substantial complication is the limited precision of remote atmospheric measurements using spectroscopy. The various metabolic processes which cause isotope fractionation cause less fractionation than anticipated measurement precision (biological fractionation is typically 2 to 7%). While this level of precision is easily reachable in the laboratory or with special in situ instruments, it is out of reach of current telescope technology to measure isotope ratios for terrestrial exoplanet atmospheres. Thus, gas isotopologues are poor biosignatures for exoplanets given our current and foreseeable technological limitations.

Recent surveys of close white dwarf binaries as well as single white dwarfs have provided evidence for the late appearance of magnetic fields in white dwarfs, and a possible generation mechanism a crystallization and rotation-driven dynamo has been suggested. A key prediction of this dynamo is that magnetic white dwarfs rotate, at least on average, faster than their non-magnetic counterparts and/or that the magnetic field strength increases with rotation. Here we present rotation periods of ten white dwarfs within 40 pc measured using photometric variations. Eight of the light curves come from TESS observations and are thus not biased towards short periods, in contrast to most period estimates that have been reported previously in the literature. These TESS spin periods are indeed systematically shorter than those of non-magnetic white dwarfs. This means that the crystallization and rotation-driven dynamo could be responsible for a fraction of the magnetic fields in white dwarfs. However, the full sample of magnetic white dwarfs also contains slowly rotating strongly magnetic white dwarfs which indicates that another mechanism that leads to the late appearance of magnetic white dwarfs might be at work, either in addition to or instead of the dynamo. The fast-spinning and massive magnetic white dwarfs that appear in the literature form a small fraction of magnetic white dwarfs, and probably result from a channel related to white dwarf mergers.

High mass-loss rates of red supergiants (RSGs) drastically affect their evolution and final fate, yet their mass-loss mechanism remains poorly understood. Various empirical prescriptions scaled with luminosity have been derived in the literature, yielding results with a dispersion of 2-3 orders of magnitude. We aim to determine an accurate mass-loss rate relation with luminosity and other parameters using a large, clean sample of RSGs and explain the discrepancy between previous works. We assembled a sample of 2,219 RSG candidates in the Large Magellanic Cloud, with ultraviolet to mid-infrared photometry in up to 49 filters. We determined the luminosity of each RSG by integrating the spectral energy distribution and the mass-loss rate using the radiative transfer code DUSTY. Our derived RSG mass-loss rates range from $10^{-9} M_\odot$ yr$^{-1}$ to $10^{-5} M_\odot$ yr$^{-1}$, depending mainly on the luminosity. The average mass-loss rate is $9.3\times 10^{-7} M_\odot$ yr$^{-1}$ for $\log{(L/L_\odot)}>4$. We established a mass-loss rate relation as a function of the luminosity and the effective temperature. Furthermore, we found a turning point in the mass-loss rate versus luminosity relation at approximately $\log{(L/L_\odot)} = 4.4$, indicating enhanced rates beyond this limit. We show that this enhancement correlates with photometric variability. Moreover, we compare our results with prescriptions from the literature, finding an agreement with those assuming steady-state winds. Additionally, we examine the effect of different assumptions on our models and found that radiatively driven winds result in higher mass loss rates by 2-3 orders of magnitude, which are unrealistically high for RSGs. Finally, we found 21% of our sample to constitute current binary candidates with a minor effect on our mass-loss relation.

The Bright Transient Survey (BTS) aims to obtain a classification spectrum for all bright ($m_\mathrm{peak}\,\leq\,18.5\,$mag) extragalactic transients found in the Zwicky Transient Facility (ZTF) public survey. BTS critically relies on visual inspection ("scanning") to select targets for spectroscopic follow-up, which, while effective, has required a significant time investment over the past $\sim5$ yr of ZTF operations. We present $\texttt{BTSbot}$, a multi-modal convolutional neural network, which provides a bright transient score to individual ZTF detections using their image data and 25 extracted features. $\texttt{BTSbot}$ is able to eliminate the need for daily human scanning by automatically identifying and requesting spectroscopic follow-up observations of new bright transient candidates. $\texttt{BTSbot}$ recovers all bright transients in our test split and performs on par with scanners in terms of identification speed (on average, $\sim$1 hour quicker than scanners). We also find that $\texttt{BTSbot}$ is not significantly impacted by any data shift by comparing performance across a concealed test split and a sample of very recent BTS candidates. $\texttt{BTSbot}$ has been integrated into Fritz and $\texttt{Kowalski}$, ZTF's first-party marshal and alert broker, and now sends automatic spectroscopic follow-up requests for the new transients it identifies. During the month of October 2023, $\texttt{BTSbot}$ selected 296 sources in real-time, 93% of which were real extragalactic transients. With $\texttt{BTSbot}$ and other automation tools, the BTS workflow has produced the first fully automatic end-to-end discovery and classification of a transient, representing a significant reduction in the human-time needed to scan. Future development has tremendous potential for creating similar models to identify and request follow-up observations for specific types of transients.

The modeling of convection is a long standing problem in stellar physics. Up-to-now, all ad hoc models rely on a free parameter alpha (among others) which has no real physical justification and is therefore poorly constrained. However, a link exists between this free parameter and the entropy of the stellar adiabat. Prescriptions, derived from 3D stellar atmospheric models, are available that provide entropy as a function of stellar atmospheric parameters (effective temperature, surface gravity, chemical composition). This can provide constraints on alpha through the development of entropy-calibrated models. Several questions arise as these models are increasingly used. Which prescription should be used? How do uncertainties impact entropy-calibrated models? We aim to study the three existing prescriptions and determine which one should be used, and how. We implemented the entropy-calibration method into the stellar evolution code Cesam2k20 and performed comparisons with the Sun and the alpha Cen system. In addition, we used data from the CIFIST grid of 3D atmosphere models to evaluate the accuracy of the prescriptions. Of the three entropy prescriptions available, we determine which one best reproduces the entropies of the 3D models. We also demonstrate that the entropy obtained from this prescription should be corrected for the evolving chemical composition and for an entropy offset different between various EoS tables, following a precise procedure, otherwise classical parameters obtained from the models will be strongly biased. Finally, we also provide table with entropy of the adiabat of the CIFIST grid, as well as fits of these entropies. We performed a precise examination of entropy-calibrated modelling, and gave recommendations on which adiabatic entropy prescription to use, how to correct it and to implement the method into a stellar evolution code.

Although the scientific principles of anthropogenic climate change are well-established, existing calculations of the warming effect of carbon dioxide rely on spectral absorption databases, which obscures the physical foundations of the climate problem. Here we show how CO2 radiative forcing can be expressed via a first-principles description of the molecule's key vibrational-rotational transitions. Our analysis elucidates the dependence of carbon dioxide's effectiveness as a greenhouse gas on the Fermi resonance between the symmetric stretch mode $\nu_1$ and bending mode $\nu_2$. It is remarkable that an apparently accidental quantum resonance in an otherwise ordinary three-atom molecule has had such a large impact on our planet's climate over geologic time, and will also help determine its future warming due to human activity. In addition to providing a simple explanation of CO2 radiative forcing on Earth, our results may have implications for understanding radiation and climate on other planets.

The hot intracluster medium (ICM) provides a unique laboratory to test multi-scale physics in numerical simulations and probe plasma physics. Utilizing archival Chandra observations, we measure density fluctuations in the ICM in a sample of 80 nearby (z<1) galaxy clusters and infer scale-dependent velocities within regions affected by mergers (r<R2500c), excluding cool-cores. Systematic uncertainties (e.g., substructures, cluster asymmetries) are carefully explored to ensure robust measurements within the bulk ICM. We find typical velocities ~220 (300) km/s in relaxed (unrelaxed) clusters, which translate to non-thermal pressure fractions ~4 (8) per cent, and clumping factors ~1.03 (1.06). We show that density fluctuation amplitudes could distinguish relaxed from unrelaxed clusters in these regions. Comparison with density fluctuations in cosmological simulations shows good agreement in merging clusters. Simulations underpredict the amplitude of fluctuations in relaxed clusters on length scales <0.75 R2500c, suggesting these systems are most sensitive to missing physics in the simulations. In clusters hosting radio halos, we examine correlations between gas velocities, turbulent dissipation rate, and radio emission strength/efficiency to test turbulent re-acceleration of cosmic ray electrons. We measure a weak correlation, driven by a few outlier clusters, in contrast to some previous studies. Finally, we present upper limits on effective viscosity in the bulk ICM of 16 clusters, showing it is systematically suppressed by at least a factor of 8, and the suppression is a general property of the ICM. Confirmation of our results with direct velocity measurements will be possible soon with XRISM.

Multi-dimensional nature of core-collapse supernova (CCSN) leads to asymmetric matter ejection and neutrino emission, that potentially accounts for the origin of neutron star (NS) kick. Asymmetric neutrino radiation fields are, in general, accompanied by large-scale inhomogeneous fluid distributions, in particular for electron-fraction ($Y_e$) distributions. Recently, it has also been revealed that lower $Y_e$ environments in proto-neutron star envelope can offer preferable conditions for collective neutrino oscillations. In this paper, we show that a dipole asymmetry of fast neutrino-flavor conversion (FFC), one of the collective neutrino oscillation modes, can power a NS kick, and that it would generate a characteristic correlation between asymmetric distributions of heavy elements in the ejecta and the direction of NS kick. We strengthen our argument for the FFC-driven NS kick mechanism by performing axisymmetric neutrino transport simulations with full Boltzmann neutrino transport. We show that this mechanism can generate sufficient linear momentum of neutrinos to account for typical proper motions of NS. Although more detailed studies are necessary, the present study opens a new channel to give a natal NS kick.

We investigate the influence of the environment on the evolution of barred and unbarred disc galaxies with a mass $>10^{10}\Msun$ from z=1 down to z=0, employing the TNG50 magnetic-hydrodynamical simulation. We find that 49% of z=1 disc galaxies undergoes a morphological transformation, transitioning into either a lenticular or spheroidal, while the other 51% retains the massive disc. The morphological alteration is mostly influenced by the environment. Lenticular and spheroidal galaxies tend to exist in denser environments and have more frequent mergers compared to disc galaxies. We find that over half of the barred galaxies (60.2%) retain the bar structure and have experienced fewer mergers compared to those galaxies that lose their bars (5.6%). These latter ones start with weaker and shorter bars at z=1 influenced by tidal interactions and are frequently observed in more populated areas. Additionally, our study reveals that less than 20% of unbarred galaxies will never develop a bar and exhibit the quietest merger history. Unbarred galaxies that experience bar formation after z=1 exhibit more frequent instances of merging events. Furthermore, tidal interactions with a close companion may account for bar formation in at least one-third of the cases. Our findings highlight that stable bars are prevalent in disc galaxies. Bar evolution may nonetheless be affected by the environment. Interactions with nearby companions or tidal forces caused by mergers have the capacity to disrupt the disc. This perturbance may materialise as the dissolution of the bar, the formation of a bar, or, in its most severe form, the complete destruction of the disc, resulting in morphological transformation. Bars that are weak and short at z=1 and undergo major or minor mergers may eventually dissolve, whereas unbarred galaxies that enter crowded environments or experience a merger may develop a bar.

Approximating neutrino oscillations as subgrid physics is an appealing prospect for simulators of core-collapse supernovae and neutron-star mergers. Because flavor instabilities quickly lead to quasisteady states in oscillation calculations, it is widely believed that flavor mixing can be approximated in astrophysical simulations by mapping unstable states onto the appropriate asymptotic ones. Subgrid models of this kind, however, are not self-consistent. The miscidynamic theory of quantum-coherent gases furnishes a subgrid model that is.

The lightcurves of asteroids are essential for determining their physical characteristics, including shape, spin, size, and surface composition. However, most asteroids are missing some of these basic physical parameters due to lack of photometric data. Although a few telescopes or surveys are specially designed for photometric lightcurve observations of asteroids, many ground-based and space-based sky surveys for hunting new exoplanets, transient events, etc., should capture numerous small Solar System objects. This will benefit the physical studies of these objects. In order to extract data of these moving objects from time-domain photometric surveys, we have developed a new method using the model tree algorithm in the field of machine learning. A dedicated module is built to automatically identify moving objects in dataset, and extract their photometric and astrometric data. As the first application of this novel method, we have analyzed data in five fields of the Yunnan-Hong Kong wide field photometric (YNHK) survey, from which 538 lightcurves of 211 asteroids are successfully extracted. Meanwhile, we also tested the method based on the data from NASA's Transiting Exoplanet Survey Satellite, and the result proves the reliability of our method. With derived lightcurves of 13 asteroids from the YNHK survey, we have determined their synodic spin periods, among which the periods of 4 asteroids are estimated for the first time. In future, we are going to apply this method to search for small objects in the outer part of the Solar System from the Chinese Space Station Telescope survey.

The timing properties of the Z-type low-mass X-ray binaries provide insights into the emission components involved in producing the unique Z-shaped track in the hardness-intensity diagrams of these sources. In this work, we investigate the AstroSat and NICER observations of the GX 340+0 covering the complete 'Z'-track from the horizontal branch (HB) to the extended flaring branch (EFB). For the first time, we present the Z-track as seen in soft X-rays using the AstroSat/SXT and NICER (the soft colour is defined as a ratio of 3-6 keV to 0.5-3 keV). The shape of the track is distinctly different in soft X-rays, strongly suggesting the presence of additional components active in soft X-rays. The detailed timing analysis revealed significant quasi-periodic oscillation throughout the HB and the normal branch (NB) using LAXPC and the first NICER detection of 33.1 +/- 1.1 Hz horizontal branch oscillation (HBO) in 3-6 keV. The oscillations at the HB/NB vertex are observed to have higher frequencies (41-52 Hz) than the HB oscillations (16-31 Hz) and NB oscillations (6.2-8 Hz) but significantly lower rms (~1.6%). The HB oscillation is also limited to the energy range of 3-20 keV, indicating an association of HBO origin with the non-thermal component. It is also supported by earlier studies that found the strongest X-ray polarisation during HB.

The FAST All Sky HI survey (FASHI) is broader in frequency band and deeper in detection sensitivity than the Arecibo Legacy Fast ALFA survey (ALFALFA). FASHI is designed to cover the entire sky observable by the Five-hundred-meter Aperture Spherical radio Telescope (FAST). To efficiently expand the OH megamaser (OHM) sample at the rest frequency of 1667.35903 MHz, we directly matched the IRAS Point Source Catalog Redshift (PSCz) catalog with the corresponding FASHI data cube. We then obtained 27 OHMs, including 9 previously known and 18 new ones, within a redshift range of $0.14314\lesssim z_{\rm OH} \lesssim0.27656$. We also measured the hyperfine ratio of nine OHMs between the 1667 and 1665 MHz lines. The ratio ranges from 1.32 to 15.22, with an average of $R_{1667:1665}=4.74$. To fit the $L_{\rm OH}$ and $L_{\rm FIR}$ relationship, we have ${\rm log}L_{\rm OH}= (1.58\pm0.10){\rm log}L_{\rm FIR}-(15.93\pm1.20)$, which is almost the same as the previous observations. All 27 OHMs have relatively high FIR luminosities, suggesting that the host system has undergone a very recent and violent starburst. The multi-wavelength SED fitting results also show that the OH host galaxies are massive ($M_* \simeq 10^{11}M_\odot$) galaxies with the SFR above the star-forming main sequence. The high $A_V$ value of the OH galaxies from the SED fit is consistent with the high IR luminosities. On the other hand, not all ULIRGs have detectable OH emission, suggesting that the OH emission may be triggered within a specific stage of the merger. In general, FAST, with its 19-beam array and UWB receiver, will be a powerful tool for observing more OHMs and unraveling their mystery in the future.

We present NICER observations of accreting X-ray pulsar 1A 0535+262 during its faint state ($\lesssim 6\times10^{36}$ erg/s) observed in several type-I and II outbursts. We discovered a transition of temporal and spectral properties around the luminosity $L_{\rm t}=3.3\times10^{35}$ erg/s, below which spectra are relatively soft and the pulse profiles have only a narrow peak. The spectra are harder and a secondary hump gradually appears in pulse profiles when $L \gtrsim L_{\rm t}$. We discuss possible physical mechanisms for this transition, including different Comptonization seed photons, the disappearance of gas shocks on the neutron star surface, and the combination of plasma and vacuum polarization effects.

The empirical relations between rotation period, chromospheric activity, and age can be used to estimate stellar age. To calibrate these relations, we present a catalog, including the masses and ages of 52,321 FGK dwarfs, 47,489 chromospheric activity index $logR^{+}_{HK}$, 6,077 rotation period $P_{rot}$ and variability amplitude $S_{ph}$, based on data from LAMOST DR7, Kepler and Gaia DR3. We find a pronounced correlation among $P_{rot}$, age, and [Fe/H] throughout the main-sequence phase for F dwarfs. However, the decrease of $logR^{+}_{HK}$ over time is not significant except for those with [Fe/H] $<$ $-$0.1. For G dwarfs, both $P_{rot}$ and $logR^{+}_{HK}$ are reliable age probes in the ranges $\sim$ 2-11 Gyr and $\sim$ 2-13 Gyr, respectively. K dwarfs exhibit a prominent decrease in $logR^{+}_{HK}$ within the age range of $\sim$ 3-13 Gyr when the relation of $P_{rot}-\tau$ is invalid. These relations are very important for promptly estimating the age of a vast number of stars, thus serving as a powerful tool in advancing the fields of exoplanet properties, stellar evolution, and Galactic-archaeology.

Jupiter Trojans (JTs) librate about the Lagrangian stationary centers L4 and L5 associated with this planet on a typically small-eccentricity and moderate-inclination heliocentric orbits. The physical and orbital properties of JTs provide important clues about the dynamical evolution of the giant planets in the early Solar System, as well as populations of planetesimals in their source regions. Here we use decade long observations from the Catalina Sky Survey (station G96) to determine the bias-corrected orbital and magnitude distributions of JTs. We distinguish the background JT population, filling smoothly the long-term stable orbital zone about L4 and L5 points, and collisional families. We find that the cumulative magnitude distribution of JTs (the background population in our case) has a steep slope for $H\leq 9$, followed with a moderately shallow slope till $H\simeq 14.5$, beyond which the distribution becomes even shallower. At $H=15$ we find a local power-law exponent $0.38\pm 0.01$. We confirm the asymmetry between the magnitude limited background populations in L4 and L5 clouds characterized by a ratio $1.45\pm 0.05$ for $H<15$. Our analysis suggests an asymmetry in the inclination distribution of JTs, with the L4 population being tighter and the L5 population being broader. We also provide a new catalog of the synthetic proper elements for JTs with an updated identification of statistically robust families (9 at L4, and 4 at L5). The previously known Ennomos family is found to consist of two, overlapping Deiphobus and Ennomos families.

The Fermi Large Area Telescope (Fermi-LAT) has been widely used to search for Weakly Interacting Massive Particle (WIMP) dark matter signals due to its unparalleled sensitivity in the GeV energy band. The leading constraints for WIMP by Fermi-LAT are obtained from the analyses of dwarf spheroidal galaxies within the Local Group, which are compelling targets for dark matter searches due to their relatively low astrophysical backgrounds and high dark matter content. In the meantime, the search for heavy dark matter with masses above TeV remains a compelling and relatively unexplored frontier. In this study, we utilize 14-year Fermi-LAT data to search for dark matter annihilation and decay signals in 8 classical dwarf spheroidal galaxies within the Local Group. We consider secondary emission caused by electromagnetic cascades of prompt gamma rays and electrons/positrons from dark matter, which enables us to extend the search with Fermi-LAT to heavier dark matter cases. We also update the dark matter subhalo model with informative priors respecting the fact that they reside in subhalos of our Milky Way halo aiming to enhance the robustness of our results. We place constraints on dark matter annihilation cross section and decay lifetime for dark matter masses ranging from $10^3$ GeV to $10^{11}$ GeV, where our limits are more stringent than those obtained by many other high-energy gamma-ray instruments.

The Nuclear Stellar Disk has been a highly active star-forming region in the Milky Way for approximately the last 30 million years. Despite hosting prominent clusters like Arches, Quintuplet, and Nuclear Stellar, their combined mass is less than 10% of the expected stellar mass, leading to the "missing cluster problem." Various factors, including high stellar density and tidal forces, contribute to this absence of clusters. Traces of dissolving clusters may exist as co-moving groups of stars, shedding light on the region's star formation history. Our analysis, utilizing integral field spectroscopy and astrophotometric data, reveals a group of young stellar objects in the NSD sharing velocities and positions, potentially indicating remnants of dissolved clusters or stellar associations. This finding contributes valuable insights into the understanding of the missing clusters problem in the Galactic center.

This thesis comprises the first three chapters dedicated to providing an overview of Gamma Ray-Bursts (GRBs), their properties, the instrumentation used to detect them, and Artificial Intelligence (AI) applications in the context of GRBs, including a literature review and future prospects. Considering both the current and the next generation of high X-ray monitors, such as Fermi-GBM and HERMES Pathfinder (an in-orbit demonstration of six 3U nano-satellites), the research question revolves around the detection of long and faint high-energy transients, potentially GRBs, that might have been missed by previous detection algorithms. To address this, two chapters introduce a new data-driven framework, DeepGRB. In Chapter 4, a Neural Network (NN) is described for background count rate estimation for X/gamma-ray detectors, providing a performance evaluation in different periods, including both solar maxima, solar minima periods, and one containing an ultra-long GRB. The application of eXplainable Artificial Intelligence (XAI) is performed for global and local feature importance analysis to better understand the behavior of the NN. Chapter 5 employs FOCuS-Poisson for anomaly detection in count rate observations and estimation from the NN. DeepGRB demonstrates its capability to process Fermi-GBM data, confirming cataloged events and identifying new ones, providing further analysis with estimates for localization, duration, and classification. The chapter concludes with an automated classification method using Machine Learning techniques that incorporates XAI for eventual bias identification.

Future observations of exoplanets will hopefully reveal detailed constraints on planetary compositions. Recently, we have developed and introduced chemcomp (Schneider & Bitsch 2021a), which simulates the formation of planets in viscously evolving protoplanetary disks by the accretion of pebbles and gas. The chemical composition of planetary building blocks (pebbles and gas) is traced by including a physical approach of the evaporation and condensation of volatiles at evaporation lines. We have now open-sourced the chemcomp code to enable comparisons between planet formation models and observational constraints by the community. The code can be found at https://github.com/AaronDavidSchneider/chemcomp, is easy to use (using configuration files) and comes with a detailed documentation and examples.

Owing to the exceptional sensitivity of the Vera C. Rubin Observatory, we predict that its upcoming LSST images will be contaminated by numerous flares from centimeter-scale space debris in Low Earth Orbits (LEO). Millisecond-duration flares from these LEO objects are expected to produce detectable image streaks of a few arcseconds with AB magnitudes brighter than 14.

We present smoothed-particle-hydrodynamics (SPH) simulations of the binary-driven hypernova (BdHN) scenario of long gamma-ray bursts (GRBs), focusing on the binary stability during the supernova (SN) explosion. The BdHN progenitor is a binary comprised of a carbon-oxygen (CO) star and a neutron star (NS) companion. The core collapse of the CO leads to an SN explosion and a newborn NS ($\nu$NS) at its center. Ejected material accretes onto the NS and the $\nu$NS. BdHNe of type I have compact orbits of a few minutes, the NS reaches the critical mass, forming a black hole (BH), and the energy release is $\gtrsim 10^{52}$ erg. BdHNe II have longer periods of tens of minutes to hours; the NS becomes more massive, remains stable, and the system releases $\sim 10^{50}$-$10^{52}$ erg. BdHN III have longer periods, even days, where the accretion is negligible, and the energy released is $\lesssim 10^{50}$ erg. We assess whether the system remains gravitationally bound after the SN explosion, leading to an NS-BH in BdHN I, an NS-NS in BdHN II and III, or if the SN explosion disrupts the system. The existence of bound systems predicts an evolutionary connection between the long and short GRB populations. We determine the binary parameters for which the binary remains bound after the BdHN event. For these binaries, we derive fitting formulas of the numerical results for the main parameters, e.g., the mass loss, the SN explosion energy, orbital period, eccentricity, center-of-mass velocity, and the relation between the initial and final binary parameters, which are useful for outlined astrophysical applications.

The thermal evolution of isolated neutron stars is a key element in unraveling their internal structure and composition and establishing evolutionary connections among different observational subclasses. Previous studies have predominantly focused on one-dimensional or axisymmetric two-dimensional models. In this study, we present the thermal evolution component of the novel three-dimensional magnetothermal code MATINS (MAgneto-Thermal evolution of Isolated Neutron Star). MATINS employs a finite volume scheme and integrates a realistic background structure, along with state-of-the-art microphysical calculations for the conductivities, neutrino emissivities, heat capacity, and superfluid gap models. This paper outlines the methodology employed to solve the thermal evolution equations in MATINS, along with the microphysical implementation which is essential for the thermal component. We test the accuracy of the code and present simulations with non-evolving magnetic fields of different topologies to produce temperature maps of the neutron star surface. Additionally, for a specific magnetic field topology, we show one fully coupled evolution of magnetic field and temperature. Subsequently, we use a ray-tracing code to link the neutron star surface temperature maps obtained by MATINS with the phase-resolved spectra and pulsed profiles that would be detected by distant observers. This study, together with our previous article focused on the magnetic formalism, presents in detail the most advanced evolutionary code for isolated neutron stars, with the aim of comparison with their timing properties, thermal luminosities and the associated X-ray light curves.

The INSPIRE project has built the largest sample of ultra-compact massive galaxies (UCMGs) at 0.1<z<0.4 and obtained their star formation histories (SFHs). Due to their preserved very old stellar populations, relics are the perfect systems to constrain the earliest epochs of mass assembly in the Universe and the formation of massive early-type galaxies. The goal of this work is to investigate whether a correlation exists between the degree of relicness (DoR), quantifying the fraction of stellar mass formed at z>2, and the other stellar population parameters.We use the Full-Index-Fitting method to fit the INSPIRE spectra to single stellar population (SSP) models. This allows us to measure, for the first time, the low-mass end slope of the IMF, as well as stellar metallicity [M/H], [Mg/Fe], [Ti/Fe] and [Na/Fe] ratios, and study correlations between them and the DoR. Similarly to normal-sized galaxies, UCMGs with larger stellar masses have overall higher metallicities. We found a correlation between the low-mass end of the IMF slope and the DoR, that, however, breaks down for systems with a more extended SFH. An even stronger dependency is found between the IMF and the fraction of mass formed at high-z. At equal velocity dispersion and metallicity, galaxies with a higher DoR have a dwarf-richer IMF than that of low-DoR counterparts. This might indicate that the cosmic epoch and formation mechanisms influence the fragmentation of the star formation cloud and hence might be the explanation for IMF variations detected in massive ETGs.

Achieving a comprehensive understanding of the star and planet formation process is one of the fundamental tasks of astrophysics, requiring detailed knowledge of the physical conditions during the different phases of this process. During the earliest stages, i.e., concerning physical processes in molecular clouds and filaments, the column density N(H2), dust temperature T and dust emissivity index \b{eta} of these objects can be derived by adopting a modified blackbody fit of the far-infrared to (sub-)millimeter spectral energy distributions. However, this often applied method is based on various assumptions. In addition, the observational basis and required, but only assumed cloud properties, such as a limited wavelength-coverage of the spectral energy distribution and dust properties, respectively, may differ between different studies. We review the basic limitations of this method and evaluate their impact on the derived physical properties of the objects of interest, i.e., molecular clouds and filaments. We find that the highest uncertainty when applying this method is introduced by the often poorly constrained dust properties. Therefore, we propose to first derive the optical depth and subsequently the column density with the help of a suitable dust model as the optical depth can be obtained with high accuracy, especially at longer wavelengths. The method provides reliable results up to the high densities and corresponding optical depths observed in molecular clouds. Considering typically used observational data, i.e., measurements obtained with far-infrared instruments like Herschel/PACS, JCMT/SCUBA-2 and SOFIA/HAWC+, data at four wavelengths are sufficient to obtain accurate results. Furthermore, we find that the dust emissivity index \b{eta} derived with this method is not suitable as an indicator of dust grain size.

One of many ways for the James-Webb Space Telescope (JWST) to capture astronomical signals is the Mid-Infrared Instrument (MIRI) Imaging mode. To make this data ready for analysis, the JWST standard reduction pipeline has three stages and many mandatory and optional steps to produce analysis-ready data. At the end of stage 3, there is a resampled 2-dimensional image for each band/wavelength, an estimated source catalog, and a source mask (segmentation image) locating these sources. This study focuses on enhancing this source mask part so that it can detect more point sources, previously cataloged after older missions, without spuriously "detecting" false positives. Combined use of the fraction of a resampled image and a derivative image seemed to improve the capability to detect unWISE catalog-located sources better than original segmentation images in 7 different real cases with the MIRI F770W filter. A few approaches are recommended to make better use of these value-sliced and derivative images.

We analyze the possibility of global anisotropy of the universe. We consider an altered FLRW metric in which there are different scale factors along the 3 different axes of space. We construct the corresponding altered Friedmann equations. We show that any initial anisotropies decrease into the future. At late times, the difference in Hubble parameters changes as $1/\sqrt{t}$ in a radiation dominated era and as $1/t$ in a matter dominated era. We use constraints from Big Bang Nucleosynthesis and the Cosmic Microwave Background to constrain the level of anisotropies at early times. We also examine how the approach back in time to the singularity is radically altered; happening much more abruptly, as a function of density, in an anisotropic universe. We also mention improved bounds that can arise from measurements of primordial gravitons, WIMPs, and neutrinos.

Neutron star low-mass X-ray binaries accrete via Roche-lobe overflow from a stellar companion that is $\lesssim$ 1 M$_{\odot}$. The accretion disk in these systems can be externally illuminated by X-rays that are reprocessed by the accreting material into an emergent reflection spectrum comprised of emission lines superimposed onto the reprocessed continuum. Due to proximity to the compact object, strong gravity effects are imparted to the reflection spectrum that can be modeled to infer properties of the NS itself and other aspects of the accreting system. This short review discusses the field of reflection modeling in neutron star low-mass X-ray binary systems with the intention to highlight the work that was awarded the 2023 AAS Newton Lacy Pierce Prize, but also to consolidate key information as a reference for those entering this subfield.

Here we present the current state of knowledge on the long-term evolution of Saturn's moon system due to tides within Saturn. First we provide some background on tidal evolution, orbital resonances and satellite tides. Then we address in detail some of the present and past orbital resonances between Saturn's moons (including the Enceladus-Dione and Titan-Hyperion resonances) and what they can tell us about the evolution of the system. We also present the current state of knowledge on the spin-axis dynamics of Saturn: we discuss arguments for a (past or current) secular resonance of Saturn's spin precession with planetary orbits, and explain the links of this resonance to the tidal evolution of Titan and a possible recent cataclysm in the Saturnian system. We also address how the moons' orbital evolution, including resonances, affects the evolution of their interiors. Finally, we summarize the state of knowledge about the Saturnian system's long-term evolution and discuss prospects for future progress.

We present a catalog of $\sim$52,000 extragalactic HII regions and their spectroscopic properties obtained using Integral Field Spectroscopy (IFS) from MUSE observations. The sample analyzed in this study contains 678 galaxies within the nearby Universe (0.004 < z < 0.06) covering different morphological types and a wide range of stellar masses (6 < log(M$_{*}$/M$_{\odot}$) < 13). Each galaxy was analyzed using the Pipe3D and pyHIIextractor codes to obtain information of the ionized gas and underlying stellar populations. Specifically, the fluxes, equivalent widths, velocities and velocity dispersions of 30 emission lines covering the wavelength range between $\lambda$4750A to $\lambda$9300A, were extracted and were used to estimate luminosity weighted ages and metallicities of the underlying stellar populations from each HII region (of the original sample we detect HII regions in 539 galaxies). In addition, we introduce and apply a novel method and independent of any intrinsic physical property to estimate and decontaminate the contribution of the diffuse ionized gas. Using the final catalog, we explore the dependence of properties of the HII regions on different local and global galaxy parameters: (i) Hubble type, (ii) stellar mass, (iii) galactocentric distance, and (iv) the age and metallicity of the underlying/neighbour stellar populations. We confirm known relations between properties of the HII regions and the underlying stellar populations (in particular with the age) uncovered using data of lower spatial and spectral resolution. Furthermore, we describe the existence of two main families of diffuse ionized gas different for galaxies host or not of HII region

We measure the line-of-sight accelerations of 26 binary pulsars due to the Milky Way's gravitational potential, and produce a 3-dimensional map of the acceleration field of the Galaxy. Acceleration measurements directly give us the change in the line-of-sight velocity at present day, without requiring any assumptions inherent to kinematic modeling. We measure the Oort limit ($\rho_0=0.062\pm0.017$ \msun/pc$^3$) and the dark matter density in the midplane ($\rho_{0,\textrm{DM}}=-0.010\pm0.018$ \msun/pc$^3$); these values are similar to, but have smaller uncertainties than previous pulsar timing measurements of these quantities. Here, we provide for the first time, values for the Oort constants and the slope of the rotation curve from direct acceleration measurements. We find that $A=15.4\pm2.6$ km/s/kpc and $B=-13.1\pm2.6$ km/s/kpc (consistent with results from \textit{Gaia}), and the slope of the rotation curve near the Sun is $-2\pm5$ km/s/kpc. We show that the Galactic acceleration field is clearly asymmetric, but due to data limitations it is not yet clear which physical processes drive this asymmetry. We provide updated models of the Galactic potential that account for various sources of disequilibrium; these models are incompatible with commonly used kinematic potentials. This indicates that use of kinematically derived Galactic potentials in precision tests (e.g., in tests of general relativity with pulsar timing) may be subject to larger uncertainties than reported. The acceleration data indicates that the mass of the Galaxy within the Solar circle is $2.3 \times 10^{11}$ M$_\odot$, roughly twice as large as currently accepted models. Additionally, the residuals of the acceleration data compared to existing Galactic models have a dependence on radial position; this trend can be explained if the Sun has an additional acceleration away from the Galactic center.

The precise definition of the lower mass limit of red supergiant stars (RSGs) is an open question in astrophysics and does not attract too much attention. Here we assemble a spectroscopic evolved cool star sample with 6,602 targets, including RSGs, asymptotic giant branch stars, and red giant branch stars, in the Large Magellanic Cloud based on \textit{Gaia} DR3 and SDSS-IV/APOGEE-2. The reference spectrum of each stellar population is built according to the quantile range of relative intensity ($1\%\sim99\%$). Five different methods, e.g., chi-square ($\chi^2$), cosine similarity (CS), machine learning (ML), equivalent width (EW), and line ratio (LR), are used in order to separate different stellar populations. The ML and $\chi^2$ provide the best and relatively consistent prediction of certain population. The derived lower limit of the RSG population is able to reach to the $\rm K_S$-band tip of red giant branch ($\rm K_S~$$\approx12.0$ mag), indicating a luminosity as low as about $10^{3.5}~L_{\sun}$, which corresponds to a stellar radius only about $100~R_{\sun}$. Given the mass-luminosity relation of $L/L_\sun =f(M/M_\sun)^3$ with $f\approx15.5\pm3$ and taking into account of the mass loss of faint RSGs up to now, the minimal initial mass of the RSG population would be about $6.1\pm0.4~M_\sun$, which is much lower than the traditional threshold of $8~M_\sun$ for the massive stars. This is the first spectroscopic evidence, indicating that the lower mass limit of RSG population is around $6~M_\sun$. However, the destinies of such faint RSGs are still elusive and may have large impact on the stellar evolutionary and supernova models.

In this paper, we report our multi-angle observations of the transverse oscillation of a prominence and a filament induced by an EUV wave originating from the farside of the Sun on 2014 September 1. The prominence oscillation was simultaneously observed by both Atmospheric Imaging Assembly (AIA) onboard the Solar Dynamics Observatory (SDO) spacecraft and Extreme-UltraViolet Imager (EUVI) onboard the Behind Solar Terrestrial Relations Observatory (STEREO) spacecraft. The speed of the shock travelling in the interplanetary space exceeds that of the EUV wave, and the coronal dimming area experiences minimal growth. This indicates that the shock wave is driven by the CME, while the EUV wave freely propagates after the lateral motion of the CME flanks has stopped. The observed oscillation direction of the prominence, determined through three-dimensional reconstruction, further supports this point. Moreover, The detailed investigation of the oscillations in the prominence and filament induced by the EUV wave reveals initial amplitudes of 16.08 and 2.15 Mm, periods of 1769 and 1863 s, damping time scales of 2640 and 1259 s, and damping ratios of 1.49 and 0.68, respectively. The radial component of magnetic field, as derived from the prominence and filament oscillation measurements, was estimated to be 5.4 G and 4.1 G, respectively. In turn, utilizing the onset times of both the prominence and filament oscillation, the average speeds of the EUV wave are determined to be 498 km s$^{-1}$ and 451 km s$^{-1}$, respectively.

Occupying the intermediate-mass regime of the accretion--jet parameter space, radio continuum emission from active galactic nuclei with black hole mass M_BH <~ 10^6 Msun (low-mass AGNs) is a valuable probe to the physics of relativistic jets. Yet the number of low-mass AGNs with radio detection is rather limited so far (~ 40 in total). In this work we make two efforts to search for radio counterparts for the largest sample of optically selected low-mass AGNs. First, we collect counterparts from the recent data releases of SKA pathfinders such as LOFAR Two-metre Sky Survey (LoTSS). Additionally, we deeply mine in Faint Images of the Radio Sky at Twenty-Centimeters (FIRST), fitting the FIRST images of the optical AGNs with an elaborate procedure optimized to detect faint radio sources. We have obtained 151 radio sources (mainly from the SKA pathfinders), including 102 new reliable sources (S/N >= 5) and 23 new candidates (3.5 <= S/N < 5). The majority of these new sources (119 of 125) have flux densities lower than the threshold of the official FIRST catalog. The new sources have rest-frame 20 cm power (P_20cm) from 1.98 x 10^20 to 1.29 x 10^23 W/Hz. For low-z Seyfert galaxies P_20cm correlates with M_BH intrinsically and positively, yet only marginally with Eddington ratio L/L_EDD. In terms of the logN--logS relation for the expanding Universe, the limiting flux density for the completeness of our LoTSS sources turns out to be 0.45 mJy at 1.4 GHz; i.e., complete to such a flux-density level that is four times deeper than the official FIRST catalog.

We have measured the spherically averaged bispectrum of the SDSS DR17 main galaxy sample, considering a volume-limited $[273\, \rm Mpc]^3$ data cube with mean galaxy number density $1.76 \times 10^{-3} \, {\rm Mpc}^{-3}$ and median redshift $0.093$. Our analysis considers $\sim 1.37 \times 10^{8}$ triangles, for which we have measured the binned bispectrum and analyzed its dependence on the size and shape of the triangle. It spans wavenumbers $k_1=(0.082-0.472)\,{\rm Mpc}^{-1}$ for equilateral triangles, and a smaller range of $k_1$ (the largest side) for triangles of other shapes. For all shapes, we find that the measured bispectrum is well modelled by a power law $A\big(k_1/1 Mpc^{-1}\big)^{n}$, where the best-fit values of $A$ and $n$ vary with the shape. The parameter $A$ is the minimum for equilateral triangles and increases as the shape is deformed to linear triangles where the two largest sides are nearly aligned, reaching its maximum value for squeezed triangles. The values of $n$ are all negative, $|n|$ is minimum $(3.31 \pm 0.17)$ for squeezed triangles, and $4.12 \pm 0.16$ for equilateral. We have also analyzed mock galaxy samples constructed from $\Lambda$CDM N-body simulations by applying a simple Eulerian bias prescription where the galaxies reside in regions where the smoothed density field exceeds a threshold. We find that the bispectrum from the mock samples with bias $b_1=1.2$ is in good agreement with the SDSS results.

A key feature of the Trappist-1 system is its monotonic decrease in bulk density with growing distance from the central star, which indicates an ice mass fraction that is zero in the innermost planets, b and c, and about 10\% in planets d through h. Previous studies suggest that the density gradient of this system could be due to the growth of planets from icy planetesimals that progressively lost their volatile content during their inward drift through the protoplanetary disk. Here we investigate the alternative possibility that the planets formed in a dry protoplanetary disk populated with pebbles made of phyllosilicates, a class of hydrated minerals with a water fraction possibly exceeding 10 wt\%. We show that the dehydration of these minerals in the inner regions of the disk and the outward diffusion of the released vapor up to the ice-line location allow the condensation of ice onto grains. Pebbles with water mass fractions consistent with those of planets d--h would have formed at the snow-line location. In contrast, planets b and c would have been accreted from drier material in regions closer to the star than the phyllosilicate dehydration line.

Starting from more than 11,200 short-period (less than 0.5 days) EW-type eclipsing binary candidates with the All-Sky Automated Survey for Supernovae (ASAS-SN) V-band light curves, we use MCMC and neural networks (NNs) to obtain the mass ratio ($q$), orbital inclination ($incl$), fill-out factor ($f$) and temperature ratio ($T_s/T_p$). After cross-matching with the Gaia DR3 database, the final sample contains parameters of 2,399 A-type and 8,712 W-type contact binaries (CBs). We present the distributions of parameters of these 11,111 short-period CBs. The mass ratio ($q$) and fill-out factor ($f$) are found to obey log-normal distributions, and the remaining parameters obey normal distributions. There is a significant period-temperature correlation of these CBs. Additionally, the temperature ratio (${T_s}$/${T_p}$) tends to increase as the orbital period decreases for W-type CBs. There is no significant correlation between them for A-type CBs. The mass ratio and fill-out factor ($q-f$) diagram suggest there is no significant correlation between these two parameters. A clear correlation exists between the mass ratio and radius ratio. The radius ratio increases with the mass ratio. Moreover, the deep fill-out CBs tend to fall on the upper boundary of the $q$$-$${R_s}$/${R_p}$ distribution, while the shallow fill-out CBs fall on the lower boundary.

We report the phase-resolved spectral results of the first Galactic Pulsating Ultra-Luminous X-ray source (PULX) Swift J0243.6+6124, modeling at its 2017-2018 outburst peak using data collected by the Hard X-ray Modulation Telescope (Insight-HXMT). The broad energy coverage of Insight-HXMT allows us to obtain more accurate spectral continuum to reduce the coupling of broad iron line profiles with other components. We use three different continuum spectrum models but obtain similar iron line results. For the first time, we detected the pulse characteristics of the broad iron line in a PULX. The variation in width and intensity of this iron line with $\sigma \sim 1.2-1.5$\,keV has a phase offset of about 0.25 from the pulse phase. We suggest that the uneven irradiation of the thick inner disk by the accretion column produces the modulated variation of the broad iron line. In addition, the non-pulsed narrow line is suggested to come from the outer disk region.

We present the detection of the three 13C isotopologs of HCCNC and HNCCC toward TMC-1 using the QUIJOTE line survey. In addition, the D species has also been detected for these two isomers of HCCCN, whereas the 15N isotopolog was only detected for HCCNC. Using high-J lines of HCCNC and HNCCC, we were able to derive very precise rotational temperatures, column densities, and subsequently the isotopic abundance ratios. We found that 12C/13C is around 90 for the three possible substitutions in both isomers. These results are slightly different from what has been found for the most abundant isomer HCCCN, for which abundances of 105, 95, and 66 were found for each one of the three possible positions of 13C. The H/D abundance ratio was found to be 31+/-4 for HCCNC and of 53+/-6 for HNCCC. The latter is similar to the H/D abundace ratio derived for HCCCN (59). The 14N/15N isotopic abundance ratio in HCCNC is 243+/-24.

The Sun has been observed through a telescope for four centuries. However, its study made a prodigious leap at the end of the nineteenth century with the appearance of photography and spectroscopy, then at the beginning of the following century with the invention of the coronagraph and monochromatic filters, and finally in the second half of the twentieth century with the advent of space exploration (satellites, probes). This makes it possible to observe the radiations hidden by the Earth's atmosphere (Ultra Violet, X-rays, $\gamma$) and to carry out ''in situ'' measurements in the solar environment. This article retraces the major stages of this fantastic epic in which renowned scientists such as Janssen, Deslandres, d'Azambuja, Lyot and Dollfus entered the scene, giving the Paris-Meudon Observatory a pioneering role in the history of solar physics until 1960. After this golden age, space exploration required large resources shared between nations, which could no longer be implemented within teams or even individual institutes. The development of numerical simulation, a new research tool, also required the pooling of supercomputers.

The spatial distribution of dwarf galaxies around their host galaxies is a critical test for the standard model of cosmology because it probes the dynamics of dark matter halos and is independent of the internal baryonic processes of galaxies. Co-moving planes-of-satellites have been found around the Milky Way, the Andromeda galaxy, and the nearby CenA galaxy, which seem to be at odds with the standard model of galaxy formation. Another nearby galaxy group, with a putative flattened distribution of dwarfs, is the M81 group. We present a quantitative analysis of the distribution of the M81 satellites using a Hough transform to detect linear structures. We confirm a flattened distribution of dwarf galaxies. Depending on the morphological type, we find a minor-to-major axis ratio of the satellite distribution to be 0.5 (all types) or 0.3 (dSph), which is in line with previous results for the M81 group. Comparing the orientation of this flattened structure in 3D with the surrounding large-scale matter distribution, we find a strong alignment with the local sheet and the planes-of-satellites around the Andromeda galaxy and Cen A. Employing velocities for a sub-sample of the dwarfs, we find no signal of co-rotation. Comparing the flattening and motion of the M81 dwarf galaxy system with IllustrisTNG50 we find good agreement between observations and simulations, but caution that i) velocity information of half of the satellites is missing, ii) velocities are coming mainly from dwarf irregulars clustered around NGC3077, which may hint towards an infall of a dwarf galaxy group and iii) some of the dwarfs may actually be tidal dwarf galaxies. From the missing velocities, we predict that the observed frequency within TNG may range between 2 to 29 per cent. Any conclusions about the agreement/disagreement with cosmological models needs to wait for a more complete picture of the dwarf galaxy system.

Alpha Centauri is a triple stellar system, and it contains the closest star to Earth (Proxima Centauri). Over the last decades, the stars in Alpha Cen and their orbits have been investigated in great detail. However, the possible scenarios for planet formation and evolution in this triple stellar system remain to be explored further. First, we present a 3D hydrodynamical simulation of the circumstellar discs in the binary Alpha Cen AB. Then, we compute stability maps for the planets within Alpha Cen obtained through N-body integrations. Last, we estimate the radial velocity (RV) signals of such planets. We find that the circumstellar discs within the binary cannot exceed 3 au in radius and that the available dust mass to form planets is about 30 $M_\oplus$. Planets around A and B are stable if their semimajor axes are below 3 au, while those around C are stable and remain unperturbed by the binary AB. For rocky planets, the planetary mass has only a mild effect on the stability. Therefore, Alpha Cen could have formed and hosted rocky planets around each star, which may be detected with RV methods in the future. The exoplanetary hunt in this triple stellar system must continue.

SS 433 is a microquasar, a stellar binary system with collimated relativistic jets. We observed SS 433 in gamma rays using the High Energy Stereoscopic System (H.E.S.S.), finding an energy-dependent shift in the apparent position of the gamma-ray emission of the parsec-scale jets. These observations trace the energetic electron population and indicate the gamma rays are produced by inverse-Compton scattering. Modelling of the energy-dependent gamma-ray morphology constrains the location of particle acceleration and requires an abrupt deceleration of the jet flow. We infer the presence of shocks on either side of the binary system at distances of 25 to 30 parsecs and conclude that self-collimation of the precessing jets forms the shocks, which then efficiently accelerate electrons.

The Spectrometer/Telescope for Imaging X-rays (STIX) on board Solar Orbiter (SolO) provides a unique opportunity to systematically perform stereoscopic X-ray observations of solar flares with current and upcoming X-ray missions at Earth. These observations will produce the first reliable measurements of hard X-ray (HXR) directivity in decades, providing a new diagnostic of the flare-accelerated electron angular distribution and helping to constrain the processes that accelerate electrons in flares. However, such observations must be compared to modelling, taking into account electron and X-ray transport effects and realistic plasma conditions, all of which can change the properties of the measured HXR directivity. Here, we show how HXR directivity, defined as the ratio of X-ray spectra at different spacecraft viewing angles, varies with different electron and flare properties (e.g., electron angular distribution, highest energy electrons, and magnetic configuration), and how modelling can be used to extract these typically unknown properties from the data. Lastly, we present a preliminary HXR directivity analysis of two flares, observed by the Fermi Gamma-ray Burst Monitor (GBM) and SolO/STIX, demonstrating the feasibility and challenges associated with such observations, and how HXR directivity can be extracted by comparison with the modelling presented here.

The equilibrium resulting in a recombining plasma with arbitrary metallicity Z, heated by a mean radiation field E as well as by sound waves dissipation due to thermal conduction, dynamic and bulk viscosity is analyzed. Generally, the heating by acoustic waves dissipation induces drastic changes in the range of temperature where the thermochemical equilibrium may exist. An additional equilibrium state appears which is characterized by a lower ionization and higher gas pressure than the equilibrium resulting when the wave dissipation is neglected. The above effects are sensibly to the values of the gas parameters as well as the wavelength and intensity of the acoustic waves. Implications in the interstellar gas, in particular, in the high velocity clouds are outlined.

The Jeans criterion is one cornerstone in our understanding of gravitational fragmentation. A critical limitation of the Jeans criterion is that the background density is assumed to be a constant, which is often not true in dynamic conditions such as star-forming regions. For example, during the formation phase of the high-density gas filaments in a molecular cloud, a density increase rate $\dot \rho$ implies a mass accumulation time of $t_{\rm acc}= \rho / \dot \rho= - \rho (\nabla \cdot (\rho \vec{v}))^{-1}$. The system is non-stationary when the mass accumulation time becomes comparable to the free-fall time $t_{\rm ff} = 1 / \sqrt{G \rho}$. We study fragmentation in non-stationary settings, and find that accretion can significantly increase in the characteristic mass of gravitational fragmentation ( $\lambda_{\rm Jeans,\; aac}= \lambda_{\rm Jeans} (1 + t_{\rm ff} / t_{\rm acc})^{1/3}$, $m_{\rm Jeans,\, acc} = m_{\rm Jeans} (1 + t_{\rm ff} / t_{\rm acc})$). In massive star-forming regions, this mechanism of transport-driven super-jeans fragmentation can contribute to the formation of massive stars by causing order-of-magnitude increases in the mass of the fragments.

We present an unbiased $\lambda$ 3 mm spectral line survey (between 84.5 and 115.8 GHz), conducted by the Purple Mountain Observatory 13.7 meter radio telescope, together with updated modeling results, towards the carbon-rich Asymptotic Giant Branch star, IRC+10216 (CW Leo). A total of 75 spectral lines (96 transitions) are detected, and identified to arise from 19 molecules: C$_2$H, $l$-C$_3$H, C$_4$H, CN, C$_3$N, HC$_3$N, HC$_5$N, HCN, HNC, CH$_3$CN, MgNC, CO, $c$-C$_3$H$_2$, SiC$_2$, SiO, SiS, CS, C$_2$S, C$_3$S, and their isotopologues. Among them, one molecular emission line (H$^{13}$CCCN $J=13-12$) is discovered in IRC+10216 for the first time. The excitation temperature, column density, and fractional abundance of the detected species are deduced by assuming they are in local thermodynamic equilibrium. In addition, the isotopic ratios of [$^{12}$C]/[$^{13}$C], [$^{32}$S]/[$^{34}$S], [$^{28}$Si]/[$^{29}$Si], and [$^{12}$C$^{34}$S]/[$^{13}$C$^{32}$S] are obtained and found to be consistent with previous studies. Finally, we summarize all of the 106 species detected in IRC+10216 to date with their observed and modeled column densities for the convenience of future studies.

In this paper, I will present the recent progress on circular-ribbon flares (CRFs) and their related activities, including coronal jets, filaments, CMEs, radio bursts, coronal dimmings, and coronal loop oscillations. Owing to the prevalence of 3D magnetic null points and the corresponding fan-spine topology in the solar atmosphere, CRFs are regularly observed in UV, EUV, and H$\alpha$ passbands. Spine reconnection and fan reconnection around the null points are predominantly responsible for the energy release and subsequent particle acceleration. Slipping reconnection at QSLs may explain the sequential brightening or rapid degradation of the circular ribbons. Periodic or quasi-periodic acceleration and precipitation of nonthermal particles in the chromosphere produce observed QPPs of CRFs in multiple altitudes as well as wavelengths. Like two-ribbon flares, the injected high-energy particles result in explosive evaporation in circular and inner ribbons, which is characterized by simultaneous blueshifts in the coronal emission lines and redshifts in the chromospheric emission lines. Homologous CRFs residing in the same active region present similar morphology, evolution, and energy partition. The peculiar topology of CRFs with closed outer spines facilitates remote brightenings and EUV late phases, which are uncommon in two-ribbon flares. Besides, CRFs are often accompanied by coronal jets, type III radio bursts, CMEs, shock waves, coronal dimmings, and kink oscillations in coronal loops and filaments. Magnetohydrodynamics numerical simulations are very helpful to understand the key problems that are still unclear up to now. Multiwavelength and multipoint observations with state-of-the-art instruments are enormously desired to make a breakthrough.

We introduce an alternative method for the calculation of sky maps from data taken with gamma-ray telescopes. In contrast to the established method of smoothing the 2D histogram of reconstructed event directions with a static kernel, we apply a Kernel Density Estimation (KDE) where the kernel size of each gamma-ray candidate is related to its estimated direction uncertainty. Exploiting this additional information implies a gain in resulting image quality, which is validated using both simulations and data. For the tested simulation and analysis configuration, the achieved improvement can only be matched with the classical approach by removing events with lower reconstruction quality, reducing the data set by a considerable amount.

M 87 is the best target for studying black hole accretion and jet formation. Reanalysis of the EHT public data at 230 GHz shows a core-knots structure at the center and jet features (Miyoshi et al. 2022a). We here compare this with the new results of GMVA at 86 GHz showing a spatially resolved central core (Lu et al. 2023a). There are similarities and differences between the two. At 86 GHz,"two bright regions" are seen on the ring in the core."Core-Knot-Westknot", triple structure in the 230 GHz image shows apparent appearance of two peaks similar to the "two bright regions" when convolved with the GMVA beam. This similarity suggests that both frequencies reveal the same objects in the core area. Protrusions are observed on both the south and north sides of the core at both frequencies, becoming prominent and wing-like at 230 GHz. The 86 GHz image shows a triple ridge jet structure, while the 230 GHz image shows only a bright central ridge with two roots. Both frequencies show a shade between the core and the central ridge. To detect the faint features from the EHT2017 data, we found that the use of all baseline data is essential. Using all including the ultrashort baseline data, revealed the jet and faint structures. Without the ultrashort baselines, these structures were not detectable. The lack of detection of any faint structures other than the ring in the M 87 data by the EHTC is presumably due to the exclusion of ultrashort baselines from their analysis.

We present a comprehensive characterization of the evolved thermally pulsing asymptotic giant branch (TP-AGB) star R Hydrae, building on the techniques applied in Stellar Evolution in Real Time I (Moln\'ar et al. 2019) to T Ursae Minoris. We compute over 3000 theoretical TP-AGB pulse spectra using MESA and GYRE and combine these with classical observational constraints and nearly 400 years of measurements of R Hya's period evolution to fit R Hya's evolutionary and asteroseismic features. Two hypotheses for the mode driving R Hya's period are considered. Solutions that identify this as the fundamental mode (FM) as well as the first overtone (O1) are consistent with observations. Using a variety of statistical tests, we find that R Hya is most likely driven by the FM and currently occupies the ``power down'' phase of an intermediate pulse (TP ~ 9-16). We predict that its pulsation period will continue to shorten for millennia. Using supplementary calculations from the Monash stellar evolution code, we also find that R Hya is likely to have undergone third dredge-up in its most recent pulse. The MESA+GYRE model grid used in this analysis includes exact solutions to the adiabatic equations of stellar oscillation for the first 10 radial-order pressure modes for every time step in every evolutionary track. The grid is fully open-source and packaged with a data visualization application. This is the first publicly available grid of TP-AGB models with seismology produced with MESA.

Detecting cosmic gamma rays at high rates is the key to time-resolve the acceleration of particles within some of the most powerful events in the universe. Time-resolving the emission of gamma rays from merging celestial bodies, apparently random bursts of gamma rays, recurring novas in binary systems, flaring jets from active galactic nuclei, clocking pulsars, and many more became a critical contribution to astronomy. For good timing on account of high rates, we would ideally collect the naturally more abundant, low energetic gamma rays in the domain of one giga electronvolt in large areas. Satellites detect low energetic gamma rays but only in small collecting areas. Cherenkov telescopes have large collecting areas but can only detect the rare, high energetic gamma rays. To detect gamma rays with lower energies, Cherenkov-telescopes need to increase in precision and size. But when we push the concept of the --far/tele-- seeing Cherenkov telescope accordingly, the telescope's physical limits show more clearly. The narrower depth-of-field of larger mirrors, the aberrations of mirrors, and the deformations of mirrors and mechanics all blur the telescope's image. To overcome these limits, we propose to record the --full/plenum-- Cherenkov-light field of an atmospheric shower, i.e. recording the directions and impacts of each individual Cherenkov photon simultaneously, with a novel class of instrument. This novel Cherenkov plenoscope can turn a narrow depth-of-field into the perception of depth, can compensate aberrations, and can tolerate deformations. We design a Cherenkov plenoscope to explore timing by detecting low energetic gamma rays in large areas.

The exoplanet diversity has been linked to the disc environment in which they form, where the host star metallicity and the formation pathways play a crucial role. In the context of the core accretion paradigm, the initial stages of planet formation require the growth of dust material from micrometre size to planetesimal size bodies before core accretion can kick in. Although numerous studies have been conducted on planetesimal formation, it is still poorly understood how this process takes place in low metallicity stellar environments. We explore planetesimals formation in stellar environments primarily with low metallicity. We performed global 1D viscous disc evolution simulations including grain growth, evaporation and condensation of chemical species at ice lines. We followed the formation of planetesimals during disc evolution and tested different metallicities, disc sizes and turbulent viscosity strengths. We find that at solar and sub-solar metallicities, there is a significant enhancement in the midplane dust-to-gas mass ratios at the ice lines but this leads to planetesimal formation only at the water ice line. In our simulations, [Fe/H] = -0.6 is the lowest limit of metallicity for planetesimal formation where a few Earth masses of planetesimals could form. For such extreme disc environments, large discs are more conducive than small discs for forming large amounts of planetesimals at a fixed metallicity, because the pebble flux can be maintained for a longer time resulting in a longer and more efficient planetesimal formation phase. At lower metallicities, planetesimal formation is less supported in quiescent discs compared to turbulent discs, because the pebble flux can be maintained for a longer time. The amount of planetesimals formed at sub-solar metallicities in our simulations places a limit on core sizes that could possibly only result in the formation of super-Earths.

We present a comprehensive analysis of the presence of Very Massive Stars (VMS > $100 M_{\odot}$) in the integrated spectra of 13 UV-bright star-forming galaxies at $2.2 \lesssim z \lesssim 3.6$ taken with the Gran Telescopio Canarias (GTC). These galaxies have very high UV absolute magnitudes ($M_{\rm UV} \simeq -24$), intense star-formation ($SFR \simeq 100-1000$ $M_{\odot}$ yr$^{-1}$), and metallicities in the range of 12+log(O/H) $\simeq8.10-8.50$ inferred from strong rest-optical lines. The GTC rest-UV spectra reveal spectral features indicative of very young stellar populations with VMS, such as strong P-Cygni line profiles in the wind lines N V $\lambda 1240$ and C IV $\lambda 1550$ along with intense and broad He II $\lambda 1640$ emission ($EW_{0}$ (HeII) $\simeq 1.40-4.60$ \r{A}). Comparison with known VMS-dominated sources and typical galaxies without VMS reveals that some UV-bright galaxies closely resemble VMS-dominated clusters (e.g., R136 cluster). The presence of VMS is further supported by a quantitative comparison of the observed strength of the He II emission with population synthesis models with and without VMS, where models with VMS are clearly preferred. Employing an empirical threshold for $EW_{0}$ (HeII) $\geq 3.0$ \r{A}, along with the detection of other VMS-related spectral profiles N IV $\lambda 1486, 1719$), we classify 9 out of 13 UV-bright galaxies as VMS-dominated sources. This high incidence of VMS-dominated sources in the UV-bright galaxy population ($\approx 70\%$) contrasts significantly with the negligible presence of VMS in typical $L_{\rm UV}^{*}$ LBGs at similar redshifts ($<1\%$). Our results thus indicate that VMS are common in UV-bright galaxies, suggesting a different, top-heavy IMF with upper mass limits between $175 M_{\odot}$ and $475 M_{\odot}$.

We compute the Gaseous Dynamical Friction (GDF) force experienced by massive perturbers on elliptical Keplerian orbits. In this paper, we investigate the density wake morphology, dynamical friction force, and secular orbital evolution for massive single perturbers as well as equal mass binaries embedded in an homogenous, static background flow. In all cases, the rate-of-change in semi-major axis is found to be negative (as expected), whereas the rate-of-change in eccentricity is negative for strictly-subsonic trajectories and positive for strictly-supersonic trajectories. Transonic orbits can experience both positive and negative torques during the course of an orbit, with some growing in eccentricity and others circularising. We observe all initial orbits becoming highly supersonic and eccentric (over sufficiently long timescales) due to a relentless semi-major axis decay increasing the Mach number and subsequent eccentricity driving. We compare our findings to previous studies for rectilinear and circular motion, while also making our data for orbital decay available.

The cosmic metals are believed to originate from stellar and supernovae (SNe) nucleosynthesis, dispersed into the interstellar medium (ISM) through stellar winds and supernova explosions. In this paper, we present the clear evidence of metal enrichment by a type Ic SN 2011jm in the galaxy NGC 4809, utilizing high spatial resolution Integral Field Units (IFU) observations obtained from the Very Large Telescope (VLT)/Multi Unit Spectroscopic Explorer (MUSE). Despite SN 2011jm being surrounded by metal-deficient ISM ($\sim 0.25 \ Z_\odot$) at a scale about 100 pc, we clearly detect enriched oxygen abundance ($\sim 0.35 \ Z_\odot$) and a noteworthy nitrogen-to-oxygen ratio at the SN site. Remarkably, the metal pollution is confined to a smaller scale ( $\leq$ 13 pc). We posit that the enhanced ionized metal stems from stellar winds emitted by massive stars or previous SNe explosions. This observation may represent the first direct detection of chemical pollution by stellar feedback in star-forming galaxies beyond the Local Volume.

Detailed timing analyses of three cataclysmic variables, namely 2MASS J09213414$-$5939068, IGR J16167$-$4957, and V667 Pup are carried out using the long-baseline and high-cadence optical photometric data from the Transiting Exoplanet Survey Satellite (TESS). Periods of 908.12$\pm$0.05 s and 990.10$\pm$0.06 s are observed in the optical variation of 2MASS J09213414$-$5939068 that were not found in earlier studies and appear to be probable spin and beat periods of the system, respectively. The presence of multiple periods at spin, beat, and other sidebands indicates that 2MASS J09213414$-$5939068 likely belongs to an intermediate polar class of magnetic cataclysmic variables that seems to be accreted via a disc-overflow mechanism. Clear evidence of a period of 582.45$\pm$0.04 s is found during the TESS observations of IGR J16167$-$4957, which can be interpreted as the spin period of the system. Strong modulation at this frequency supports its classification as an intermediate polar, where accretion may primarily be governed by a disc. The dominance of the spin pulse unveils the disc-fed dominance accretion in V667 Pup, but the detection of the previously unknown beat period of 525.77$\pm$0.03 s suggests that a portion of the material is also accreted through a stream. Moreover, the double-peaked structure observed in the optical spin pulse profile of V667 Pup suggests the possibility of a two-pole accretion geometry, where each pole accretes at a different rate and is separated by 180$^\circ$.

Context: Planet-disk interactions play a crucial role in the understanding of planet formation and disk evolution. There are multiple numerical tools available to simulate these interactions, including the often-used FARGO code and its variants. Many of the codes were extended over time to include additional physical processes with a focus on their accurate modeling. Aims: We introduce FargoCPT, an updated version of FARGO incorporating other previous enhancements to the code, to provide a simulation environment tailored to study interactions between stars, planets, and disks, ensuring accurate representation of planet systems, hydrodynamics, and dust dynamics with a focus on usability. Methods: The radiation-hydrodynamics part of FargoCPT uses a second-order upwind scheme in 2D polar coordinates supporting multiple equations of state, radiation transport, heating and cooling, and self-gravity. Shocks are considered using artificial viscosity. Integration of the N-body system is achieved by leveraging the REBOUND code. The dust module utilizes massless tracer particles, adapting to drag laws for the Stokes and Epstein regimes. Moreover, FargoCPT provides mechanisms to simulate accretion onto the stars and planets. Results: The code has been tested in practice by its use in various publications. Additionally, it comes with an automated test suite to test the physics modules. Conclusions: FargoCPT offers a unique set of simulation capabilities within the current landscape of publicly available planet-disk interaction simulation tools. Its structured interface and underlying technical updates are intended to assist researchers in the ongoing exploration of planet formation.

Large-scale characterization of exoplanetary atmospheres is on the horizon, thereby making it possible in the future to extract their statistical properties. In this context, by using a well validated model in the solar system, we carry out three-dimensional magnetohydrodynamic simulations to compute nonthermal atmospheric ion escape rates of unmagnetized rocky exoplanets as a function of their radius based on fixed stellar radiation and wind conditions. We find that the atmospheric escape rate is, unexpectedly and strikingly, a nonmonotonic function of the planetary radius $R$ and that it evinces a maximum at $R \sim 0.7\,R_\oplus$. This novel nonmonotonic behavior may arise from an intricate tradeoff between the cross-sectional area of a planet (which increases with size, boosting escape rates) and its associated escape velocity (which also increases with size, but diminishes escape rates). Our results could guide forthcoming observations because worlds with certain values of $R$ (such as $R \sim 0.7\,R_\oplus$) might exhibit comparatively higher escape rates when all other factors are constant.

We present a new catalog of stars for which detected solar-like oscillations and magnetic activity measurements are both available from chromospheric spectroscopic observations. Our results were obtained by exploiting NASA TESS mission light curves for active stars observed within the Mount Wilson Observatory HK project and the HK survey of the Hamburg Robotic Telescope TIGRE. We analyzed the light curves for a total of 191 stars by adopting recent techniques based on Bayesian analysis and model comparison to assess the detection of a power excess originating from solar-like oscillations. We characterized the oscillations in a total of 34 targets, for which we provide estimates for the global asteroseismic parameters of $\nu_\mathrm{max}$ (the frequency of maximum oscillation power), $\Delta\nu$ (the large frequency separation), and for the amplitude of the solar-like oscillation envelope $A_\mathrm{max}$. We provide strong statistical evidence for the detection of solar-like oscillations in 15 stars of our sample, identify six further stars where a detection is likely, and 13 stars for which oscillations cannot be ruled out. The key parameters extracted in this work will be exploited for a detailed stellar modeling of the targets and to calibrate relations that connect the level of the measured magnetic activity to the suppression induced on the global oscillation amplitudes. This opens the possibility of shedding light on the interplay between magnetic fields and oscillations. Because of their relatively high brightness, the targets may also be of interest for future dedicated follow-up observations using both photometry and spectropolarimetry.

X-ray binaries hosting a compact object have been among the main targets of the Imaging X-ray Polarimetry Explorer (IXPE) since its launch, due to their high brightness in the 2-8 keV energy band. The spectropolarimetric analysis performed so far has proved to be of great importance in providing constraints on the accretion geometry of these systems. However, the data statistics is not enough to unambiguously disentangle the contribution of the single components to the net observed polarimetric signal. In this work, we aim to present a model for computing the polarization degree and polarization angle of the boundary layer around weakly magnetized neutron stars in low-mass X-ray binaries in the soft state. The main motivation is to provide strong theoretical support to data interpretation of observations performed by IXPE or future satellites for X-ray polarimetry. The results were obtained by modeling the boundary layer as an equatorial belt around the compact object and locally approximating it as a plane-parallel scattering atmosphere, for which the associated radiative transfer equation for polarized radiation in the Thomson limit was solved. The polarimetric quantities were then transformed from the comoving frame to the observer frame using the numerical methods formerly developed for X-ray pulsars. For typical values of the optical depth and electron temperature of the boundary layer of these systems in a soft state, the polarization degree was less then 0.5\%, while the polarization angle was rotated by $\protect \la 5^{\circ}$ with respect to the neutron star spin axis due to special and general relativistic effects for fast rotation, the amount progressively decreasing for lower spin frequencies. The derived quantities can be used to remove degeneracy when multicomponent spectropolarimetry is performed.

We demonstrate that astrophysical constraints on the dense-matter equation of state place an upper bound on the color-superconducting gap in dense matter above the transition from nuclear matter to quark matter. Pairing effects in the color-flavor locked (CFL) quark matter phase increase the pressure at high density, and if this effect is sufficiently large then the requirements of causality and mechanical stability make it impossible to reach such a pressure in a way that is consistent with what is known at lower densities. The intermediate-density equation of state is inferred by considering extensions of chiral effective field theory (CEFT) to neutron star densities, and conditioning these using current astrophysical observations of neutron star radius, maximum mass, and tidal deformability (PSR J0348+0432, PSR J1624-2230, PSR J0740+6620, GW170817). At baryon number chemical potential $\mu = 2.6~\text{GeV}$ we find a 95% upper limit on the CFL pairing gap $\Delta$ of $457~\text{MeV}$ using overly conservative assumptions and $216~\text{MeV}$ with more reasonable assumptions. This constraint may be strengthened by future astrophysical measurements as well as by future advances in high density QCD calculations.

We seek to understand the evolution of Wolf-Rayet central stars by comparing the diffuse X-ray emission from their wind-blown bubbles with that from their hydrogen-rich counterparts with predictions from hydrodynamical models. We simulate the dynamical evolution of heat-conducting wind-blown bubbles using a post-AGB-model of 0.595~Msun, allowing for variations of its evolutionary timescale and wind power. For Wolf-Rayet central stars, the wind is hydrogen-poor, more dense, and slower compared to O-type central stars. We use the CHIANTI software to compute the X-ray properties of bubble models along the evolutionary paths and explicitly allow for non-equilibrium ionisation of key chemical elements. A sample of 12 planetary nebulae with diffuse X-ray emission -- seven harbouring an O-type and five a Wolf-Rayet nucleus -- is used to test the bubble models. The properties of most hydrogen-rich bubbles (X-ray temperature, X-ray luminosity, size) and their central stars (photon and wind luminosity) are fairly well represented by bubble models of our 0.595~Msun AGB remnant. The bubble evolution of Wolf-Rayet objects is different, thanks to the high radiation cooling of their carbon- and oxygen-rich winds. The bubble formation is delayed, and eventually, evaporation begins, leading to chemically stratified bubbles. The bubbles of the youngest Wolf-Rayet objects appear chemically uniform, the chemically stratified bubbles of the evolved Wolf-rayet objects have excessively low characteristic temperatures that cannot be explained by our modelling. The formation of nebulae with O-type nuclei follows mainly a single path, but the formation pathways leading to the Wolf-Rayet-type objects appear diverse.

Space-based observations show that the solar atmosphere from the solar chromosphere to the solar corona is filled with small-scale jets and is linked with small-scale explosions. These jets may be produced by mechanisms similar to that of large-scale flares and such jets may be related to the heating of corona and chromosphere as well as the acceleration of solar wind. The chromospheric anemone jets on the Sun remain puzzling because their footpoints (or bright knots) have not been well resolved and the formation process of such enigmatic small-scale jets remains unclear. We propose a new model for chromospheric jets using the three-dimensional magnetohydrodynamic (MHD) simulations, which show that the continuous, upward rising of small-scale twisted magnetic flux ropes in a magnetized solar chromosphere drive small-scale magnetic reconnection and the launching of several small-scale jets during the evolution of the chromospheric anemone jets. Our new, self-consistent, three-dimensional computer modeling of small-scale, but ever-changing flux rope emergence in the magnetized solar atmosphere is fully consistent with observations and provides a universal mechanism for nanoflare and jet formation.

Almost all current and future high-contrast imaging instruments will use a Pyramid wavefront sensor (PWFS) as a primary or secondary wavefront sensor. The main issue with the PWFS is its nonlinear response to large phase aberrations, especially under strong atmospheric turbulence. Most instruments try to increase its linearity range by using dynamic modulation, but this leads to decreased sensitivity, most prominently for low-order modes, and makes it blind to petal-piston modes. In the push toward high-contrast imaging of fainter stars and deeper contrasts, there is a strong interest in using the PWFS in its unmodulated form. Here, we present closed-loop lab results of a nonlinear reconstructor for the unmodulated PWFS of the Magellan Adaptive Optics eXtreme (MagAO-X) system based on convolutional neural networks (CNNs). We show that our nonlinear reconstructor has a dynamic range of >600 nm root-mean-square (RMS), significantly outperforming the linear reconstructor that only has a 50 nm RMS dynamic range. The reconstructor behaves well in closed loop and can obtain >80% Strehl at 875 nm under a large variety of conditions and reaches higher Strehl ratios than the linear reconstructor under all simulated conditions. The CNN reconstructor also achieves the theoretical sensitivity limit of a PWFS, showing that it does not lose its sensitivity in exchange for dynamic range. The current CNN's computational time is 690 microseconds, which enables loop speeds of >1 kHz. On-sky tests are foreseen soon and will be important for pushing future high-contrast imaging instruments toward their limits.

Relativistically energetic, non-thermal population of electrons can give rise to unique CMB spectral distortion signatures which can be significantly different from thermal Sunyaev-Zeldovich signal or $y$-distortion. These signatures depend upon the spectrum of non-thermal electrons, therefore, a detection can inform us about the existence and abundance of non-thermal electrons in our universe. Using public CMB maps and data, we derive upper limits on non-thermal $y$-parameter for a relativistic, power-law electron distribution. With future CMB experiments, we may be in a position to detect or put significantly tighter constraints on these signals which can affect our understanding of non-thermal electron distributions in our universe.

We present JWST MIRI Medium Resolution Spectrograph (MRS) observations of the $\beta$ Pictoris system. We detect an infrared excess from the central unresolved point source from 5 to 7.5 $\mu$m which is indicative of dust within the inner $\sim$7 au of the system. We perform PSF subtraction on the MRS data cubes and detect a spatially resolved dust population emitting at 5 $\mu$m. This spatially resolved hot dust population is best explained if the dust grains are in the small grain limit (2$\pi$a$\ll$$\lambda$). The combination of unresolved and resolved dust at 5 $\mu$m could suggest that dust grains are being produced in the inner few au of the system and are then radiatively driven outwards, where the particles could accrete onto the known planets in the system $\beta$ Pic b and c. We also report the detection of an emission line at 6.986 $\mu$m that we attribute to be [Ar II]. We find that the [Ar II] emission is spatially resolved with JWST and appears to be aligned with the dust disk. Through PSF subtraction techniques, we detect $\beta$ Pic b at the 5$\sigma$ level in our MRS data cubes and present the first mid-IR spectrum of the planet from 5 to 7 $\mu$m. The planet's spectrum is consistent with having absorption from water vapor between 5 and 6.5 $\mu$m. We perform atmosphere model grid fitting on spectra and photometry of $\beta$ Pic b and find that the planet's atmosphere likely has a sub-stellar C/O ratio.

We have developed a suite of novel infrared-blocking filters made by embedding diamond scattering particles in a polyimide aerogel substrate. Our developments allow us to tune the spectral performance of the filters based on both the composition of the base aerogel material and the properties of the scattering particles. Our filters are targeted for use in a variety of applications, from ground-based cryogenic telescope experiments to space-based planetary science probes. We summarize the design, fabrication, and characterization of these filters. We investigate several polyimide base aerogel formulations and the effects of loading them with diamond scattering particles of varying sizes and relative densities.

Context. Jets are commonly detected in post-asymptotic giant branch (post-AGB) binaries and originate from an accretion process onto the companion of the post-AGB primary. These jets are revealed by high-resolution spectral time series. Aims. This paper is part of a series. In this work, we move away from our previous parametric modelling and include a self-similar wind model that allows the physical properties of post-AGB binaries to be characterised. This model describes magnetically driven jets from a thin accretion disk threaded by a large-scale, near equipartition vertical field. Methods. We expanded our methodology in order to simulate the high-resolution dynamic spectra coming from the obscuration of the primary by the jets launched by the companion. We present the framework to exploit the self-similar jet models for post-AGB binaries. We performed a parameter study to investigate the impact of different parameters (inclination, accretion rate, inner and outer launching radius) on the synthetic spectra. Results. We successfully included the physical jet models into our framework. The synthetic spectra have a very similar orbital phase coverage and absorption strengths as the observational data. The magnetohydrodynamic (MHD) jet models provide a good representation of the actual jet creation process in these evolved binaries. Challenges remain, however, as the needed high-accretion rate would induce accretion disks that are too hot in comparison to the data. Moreover, the rotational signature of the models is not detected in the observations. In future research, we will explore models with a higher disk ejection efficiency and even lower magnetisation in order to solve some of the remaining discrepancies between the observed and synthetic dynamic spectra.

To understand its evolution and the effects of its eruptive events, the Sun is permanently monitored by multiple satellite missions. The optically-thin emission of the solar plasma and the limited number of viewpoints make it challenging to reconstruct the geometry and structure of the solar atmosphere; however, this information is the missing link to understand the Sun as it is: a three-dimensional evolving star. We present a method that enables a complete 3D representation of the uppermost solar layer (corona) observed in extreme ultraviolet (EUV) light. We use a deep learning approach for 3D scene representation that accounts for radiative transfer, to map the entire solar atmosphere from three simultaneous observations. We demonstrate that our approach provides unprecedented reconstructions of the solar poles, and directly enables height estimates of coronal structures, solar filaments, coronal hole profiles, and coronal mass ejections. We validate the approach using model-generated synthetic EUV images, finding that our method accurately captures the 3D geometry of the Sun even from a limited number of 32 ecliptic viewpoints ($|\text{latitude}| \leq 7^\circ$). We quantify uncertainties of our model using an ensemble approach that allows us to estimate the model performance in absence of a ground-truth. Our method enables a novel view of our closest star, and is a breakthrough technology for the efficient use of multi-instrument datasets, which paves the way for future cluster missions.

Water is an important component of exoplanets, with its distribution, i.e., whether at the surface or deep inside, fundamentally influencing the planetary properties. The distribution of water in most exoplanets is determined by yet-unknown partitioning coefficients at extreme conditions. Our new first-principles molecular dynamics simulations reveal that water strongly partitions into iron over silicate at high pressures and thus would preferentially stay in a planet's core. Furthermore, we model planet interiors by considering the effect of water on density, melting temperature, and water partitioning. The results shatter the notion of water worlds as imagined before: the majority of the bulk water budget (even more than 95%) can be stored deep within the core and the mantle, and not at the surface. For planets more massive than ~6 Earth's mass and Earth-size planets (of lower mass and small water budgets), the majority of water resides deep in the cores of planets, Whether water is assumed to be at the surface or at depth can affect the radius by up to 25% for a given mass. This has drastic consequences for the inferred water distribution in exoplanets from mass-radius data.

The Cadmium Zinc Telluride Imager (CZTI) aboard AstroSat has good sensitivity to Gamma Ray Bursts (GRBs), with close to 600 detections including about 50 discoveries undetected by other missions. However, CZTI was not designed to be a GRB monitor and lacks localisation capabilities. We introduce a new method of localising GRBs using "shadows" cast on the CZTI detector plane due to absorption and scattering by satellite components and instruments. Comparing the observed distribution of counts on the detector plane with simulated distributions with the AstroSat Mass Model, we can localise GRBs in the sky. Our localisation uncertainty is defined by a two-component model, with a narrow Gaussian component that has close to 50% probability of containing the source, and the remaining spread over a broader Gaussian component with an 11.3 times higher $\sigma$. The width ($\sigma$) of the Gaussian components scales inversely with source counts. We test this model by applying the method to GRBs with known positions and find good agreement between the model and observations. This new ability expands the utility of CZTI in the study of GRBs and other rapid high-energy transients.

In this work we investigate the merger rate of primordial black hole-neutron star (PBH-NS) binaries in two widely-studied modified gravity (MG) models: Hu-Sawicki $f(R)$ gravity and the normal branch of Dvali-Gabadadze-Porrati (nDGP) gravity. In our analysis, we take into account the effects of MG on the halo properties including halo mass function, halo concentration parameter, halo density profile, and velocity dispersion of dark matter particles. We find that these MG models, due to their stronger gravitational field induced by an effective fifth force, predict enhanced merger rates compared to general relativity. This enhancement is found to be redshift-dependent and sensitive to model parameters, PBH mass and fraction. Assuming PBH mass range of $5-50 M_{\odot}$, we compare the predicted merger rate of PBH-NS binaries with those inferred from LIGO-Virgo-KAGRA observations of gravitational waves (GWs). We find that the merger rates obtained from MG models will be consistent with the GW observations, if the abundance of PBHs is relatively large, with the exact amount depending on the MG model and its parameter values, as well as PBH mass. We also establish upper limits on the abundance of PBHs in these MG frameworks while comparing with the existing non-GW constraints, which can potentially impose even more stringent constraints.

We study the effect of box size on the matter power spectrum obtained via cosmological N-body simulations. Within the framework of the cosmic screening approach, we show that the relative deviation between the spectra for our largest comoving box with L = 5632 Mpc/h and those for L = 280, 560, 1680, 4480, 5120 Mpc/h boxes consistently increases with decreasing box size in the latter set in the redshift range $0\leq z\leq 80$ for the considered values. As an additional demonstrative example, at redshift zero, we determine the values $k_{1\%}$ corresponding to the modes at which relative deviations reach 1\%.

Transient signals arising from instrumental or environmental factors, commonly referred to as glitches, constitute the predominant background of false alarms in the detection of gravitational waves in data collected from ground-based detectors. Therefore, effective data analysis methods for vetoing glitch-induced false alarms are crucial to enhancing the sensitivity of a search for gravitational waves. We present a veto method for glitches that impact matched filtering-based searches for binary inspiral signals. The veto uses unphysical sectors in the space of chirp time parameters as well as an unphysical extension including negative chirp times to efficiently segregate glitches from gravitational wave signals in data from a single detector. Inhabited predominantly by glitches but nearly depopulated of genuine gravitational wave signals, these unphysical sectors can be efficiently explored using Particle Swarm Optimization. In a test carried out on data taken from both LIGO detectors spanning multiple observation runs, the veto was able to reject $99.9\%$ of glitches with no loss of injected signals detected with a signal-to-noise ratio $\geq 9.0$. Our results show that extending a matched filter search to unphysical parts of a signal parameter space promises to be an effective strategy for mitigating glitches.

As with the laser interferometer gravitational-wave observatory (LIGO), the matched filtering technique will be critical to the data analysis of gravitational wave detection by space-based detectors, including LISA, Taiji and Tianqin. Waveform templates are the basis for such matched filtering techniques. To construct ready-to-use waveform templates, numerical relativity waveforms are a starting point. Therefore, the accuracy issue of numerical relativity waveforms is critically important. There are many investigations regarding this issue with respect to LIGO. But unfortunately there are few results on this issue with respect to space-based detectors. The current paper investigates this problem. Our results indicate that the existing numerical relativity waveforms are as accurate as 99% with respect to space-based detectors, including LISA, Taiji and Tianqin. Such an accuracy level is comparable to that with respect to LIGO.

Antonie (Anton) Pannekoek (1873-1960) is remembered as one of the initiators of the field of stellar atmospheres. A second part of his research concerned Galactic astronomy. He was convinced that the sidereal system was built up of clouds of stars in a smooth, low-density stratum. In addition there were dark clouds together with streaks with little or no extinction in between. Pannekoek looked at bright star clouds and estimated their distance from their contribution to star counts. He found values of tens of kpc, which would mean their distribution was similar in extent to that of Shapleys globular cluster system. Later he had to reduce his distance by a factor over two, and later still retract the method. He developed a rigorous method of estimating distances of dark clouds from modeling star counts off and on the cloud, preceding Wolf's quick and dirty method. He should have received more credit for this. He started isophotal maps of the northern and southern Milky Way, first from visual observations, later from photographic surface photometry using out-of-focus exposures. I compare Pannekoeks maps with detailed photographic surface photometry of the south by the group in Bochum and to the almost all-sky mapping by the Pioneer 10 spacecraft, free of zodiacal light, from beyond the asteroid belt. This shows Panneloeks maps to be surprisingly accurate. The legacy of Pannekoek in the area of Galactic research consists of his mapping of the structure of the nearby part of the Galaxy, the distances of dark clouds, and isophotal maps of the Milky Way. His other contributions turned out inconclusive or wrong as a result of his conviction, resulting from his many years of observing and mapping the Milky Way, that the nearby distribution is characterized primarily by more or less isolated clouds of stars and by dust restricted to isolated dark clouds and streaks.

Considering accretion onto a charged dilaton black hole, the fundamental equations governing accretion, general analytic expressions for critical points, critical velocity, critical speed of sound, and ultimately the mass accretion rate are obtained. A new constraint on the dilation parameter coming from string theory is found and the case for polytropic gas is delved into a detailed discussion. It is found that the dialtion and the adiabatic index of accreted material have deep effects on the accretion process.

We investigate graviton-photon oscillations sourced by cosmological magnetic fields from Gertsenshtein effect. We adopt a robust perturbative approach and we find that the conversion probability from graviton to photon can be resonantly enhanced in monochromatic, multi-chromatic and scale invariant spectrum models of stochastic magnetic field fluctuations. In addition, the expansion of the Universe acts as a decoherence factor, which demands a natural discretization scheme along the line of sight. Including also decoherence from cosmic acceleration, we find that conversion probabilities for stochastic magnetic fields are completely different than results predicted from existing magnetic domain-like models in a wide range of magnetic strengths and coherence lengths. Resonances can be tested by radio telescopes as a probe of high frequency gravitational wave sources and primordial magnetogenesis mechanisms.

Stimulated radiation and gravitational waves (GWs) are two of the most important predictions made by Albert Einstein. In this work, we demonstrate that stimulated GW radiation can occur within gravitational atoms, which consist of Kerr black holes and the surrounding boson clouds formed through superradiance. The presence of GWs induces mixing between different states of the gravitational atoms, leading to resonant transitions between two states when the GW wavenumber closely matches the energy difference. Consequently, the energy and angular momentum released from these transitions lead to the amplification of GWs, resulting in an exponential increase in the transition rate. Remarkably, the transitions complete within a much shorter time compared to the lifetime of the cloud. These stimulated transitions give rise to a novel GW signal that is strong and directed, distinguished from the previously predicted continuous GWs originating from clouds of ultralight bosons.

Numerous tabletop experiments have been dedicated to exploring the manifestations of screened scalar field dark energy, such as symmetron or chameleon fields. Precise theoretical predictions require simulating field configurations within the respective experiments. This paper focuses onto the less-explored environment-dependent dilaton field, which emerges in the strong coupling limit of string theory. Due to its exponential self-coupling, this field can exhibit significantly steeper slopes compared to symmetron and chameleon fields, and the equations of motion can be challenging to solve with standard machine precision. We present the first exact solution for the geometry of a vacuum region between two infinitely extended parallel plates. This solution serves as a benchmark for testing the accuracy of numerical solvers. By reparametrizing the model and transforming the equations of motion, we show how to make the model computable across the entire experimentally accessible parameter space. To simulate the dilaton field in one- and two-mirror geometries, as well as spherical configurations, we introduce a non-uniform finite difference method. Additionally, we provide an algorithm for solving the stationary Schr\"odinger equation for a fermion in one dimension in the presence of a dilaton field. The algorithms developed here are not limited to the dilaton field, but can be applied to similar scalar-tensor theories as well. We demonstrate such applications at hand of the chameleon and symmetron field. Our computational tools have practical applications in a variety of experimental contexts, including gravity resonance spectroscopy (qBounce), Lunar Laser Ranging (LLR), and the upcoming Casimir and Non-Newtonian Force Experiment (CANNEX). A Mathematica implementation of all algorithms is provided.

We examine the justification for taking the Event Horizon Telescope's famous 2019 image to be a reliable representation of the region surrounding a black hole. We argue that it takes the form of a robustness argument, with the resulting image being robust across variation in a range of data-analysis pipelines. We clarify the sense of "robustness" operating here and show how it can account for the reliability of astrophysical inferences, even in cases -- like the EHT -- where these inferences are based on experiments that are (for all practical purposes) unique. This has consequences far beyond the 2019 image.