Papers for Friday, Jan 19 2024

The primary purpose of this paper is to see how well a recently proposed new model fits (a) the position of the baryon acoustic oscillations (BAO) features observed in the large-scale distribution of galaxies and (b) the angular size measured for the sound horizon due to BAO imprinted in the cosmic microwave background (CMB) anisotropy. The new model is a hybrid model that combines the tired light (TL) theory with a variant of the ${\Lambda}CDM$ model in which the cosmological constant is replaced with a covarying coupling constants' (CCC) parameter ${\alpha}$. This model, dubbed the CCC+TL model, can fit the supernovae type 1a Pantheon+ data as accurately as the ${\Lambda}CDM$ model, and also fit the angular size of cosmic dawn galaxies observed by the James Webb Space Telescope, which is in tension with the ${\Lambda}CDM$ model. The results we obtained are $151.0 (\pm5.1)$ Mpc for the absolute BAO scale at the current epoch, and the angular size of the sound horizon ${\theta}_{sh}=0.60{\deg}$ matching Planck's observations at the surface of the last scattering when the baryon density is set to 100% of the matter density and |${\alpha}$| is increased by 5.6%. It remains to be seen if the new model is consistent with the CMB power spectrum, the big-bang nucleosynthesis of light elements, and other critical observations.

V1741 Sgr (= SPICY 71482/Gaia22dtk) is a Classical T Tauri star on the outskirts of the Lagoon Nebula. After at least a decade of stability, in mid-2022, the optical source brightened by ~3 mag over two months, remained bright until early 2023, then dimmed erratically over the next four months. This event was monitored with optical and infrared spectroscopy and photometry. Spectra from the peak (October 2022) indicate an EX Lup-type (EXor) accretion outburst, with strong emission from H I, He I, and Ca II lines and CO bands. At this stage, spectroscopic absorption features indicated a temperature of T ~ 4750 K with low-gravity lines (e.g., Ba II and Sr II). By April 2023, with the outburst beginning to dim, strong TiO absorption appeared, indicating a cooler T ~ 3600 K temperature. However, once the source had returned to its pre-outburst flux in August 2023, the TiO absorption and the CO emission disappeared. When the star went into outburst, the source's spectral energy distribution became flatter, leading to bluer colours at wavelengths shorter than ~1.6 microns and redder colours at longer wavelengths. The brightening requires a continuum emitting area larger than the stellar surface, likely from optically thick circumstellar gas with cooler surface layers producing the absorption features. Additional contributions to the outburst spectrum may include blue excess from hotspots on the stellar surface, emission lines from diffuse gas, and reprocessed emission from the dust disc. Cooling of the circumstellar gas would explain the appearance of TiO, which subsequently disappeared once this gas had faded and the stellar spectrum reemerged.

Recent works have uncovered an excess signal in the parity-odd four-point correlation function measured from the BOSS spectroscopic galaxy survey. If physical in origin, this could indicate evidence for new parity-breaking processes in the scalar sector, most likely from inflation. At heart, these studies compare the observed four-point correlator to the distribution obtained from parity-conserving mock galaxy surveys; if the simulations underestimate the covariance of the data, noise fluctuations may be misinterpreted as a signal. To test this, we reanalyse the BOSS CMASS + LOWZ parity-odd dataset with the noise distribution modeled using the newly developed GLAM-Uchuu suite of mocks. These comprise full N-body simulations that follow the evolution of $2000^3$ dark matter particles in a $\Lambda$CDM universe, and represent a significant upgrade compared to the formerly MultiDark-Patchy mocks, which were based on an alternative (non N-body) gravity solver. We find no significant evidence for parity-violation in the BOSS dataset (with a baseline detection significance of $1.4\sigma$), suggesting that the former signal ($>3.5\sigma$ with our data cuts) could be caused by an underestimation of the covariance in MultiDark-Patchy. The significant differences between results obtained with the two sets of BOSS-calibrated galaxy catalogs showcases the heightened sensitivity of beyond-two-point analyses to the treatment of non-linear effects and indicates that previous constraints may suffer from large systematic uncertainties.

While the shape of the Ly$\alpha$ profile is viewed as one of the best tracers of ionizing-photon escape fraction ($f_{esc}$) within low redshift (z~0.3) surveys of the Lyman continuum, this connection remains untested at high redshift. Here, we combine deep, rest-UV Keck/LRIS spectra of 80 objects from the Keck Lyman Continuum Spectroscopic Survey with rest-optical Keck/MOSFIRE spectroscopy in order to examine potential correlations between Ly$\alpha$ profile shape and the escape of ionizing radiation within z~3 star-forming galaxies. We measure the velocity separation between double-peaked Ly$\alpha$ emission structure (v$_{\rm sep}$), between red-side Ly$\alpha$ emission peaks and systemic (v$_{\rm Ly\alpha,red}$), and between red-side emission peaks and low-ionization interstellar absorption lines (v$_{\rm Ly\alpha-LIS}$). We find that the IGM-corrected ratio of ionizing to non-ionizing flux density is significantly higher in KLCS objects with lower v$_{\rm Ly\alpha,red}$. We find no significant trend between measures of ionizing-photon escape and v$_{\rm Ly\alpha-LIS}$. We compare our results to measurements of z~0.3 "Green Peas" from the literature and find that KLCS objects have larger v$_{\rm sep}$ at fixed v$_{\rm Ly\alpha,red}$, larger $f_{esc}$ at fixed v$_{\rm Ly\alpha,red}$, and higher v$_{\rm Ly\alpha,red}$ overall than z~0.3 analogs. We conclude that the Ly$\alpha$ profile shapes of our high-redshift sources are fundamentally different, and that measurements of profile shape such as v$_{\rm Ly\alpha,red}$ map on to $f_{esc}$ in different ways. We caution against building reionization-era $f_{esc}$ diagnostics based purely on Ly$\alpha$ profiles of low-redshift dwarf galaxies. Tracing v$_{\rm sep}$, v$_{\rm Ly\alpha,red}$, and $f_{esc}$ in a larger sample of z~3 galaxies will reveal how these variables may be connected for galaxies at the epoch of reionization.

The interacting binary system NGC 4490/85 = Arp 269 is intermediate in mass between the Milky Way/Large Magellanic Cloud and the Large/Small Magellanic Cloud binary systems. It is a system of 14 known galaxies. We estimate the total Newtonian gravitating mass of the NGC\,4490/85 group to be M_T = (1.37 +- 0.43) times 10^{12} M_Sun using radial velocities and projected separations of its 13~candidate members. The system of dwarf satellites in the group demonstrates signs of coherent rotation in the same direction as that of the extended HI-shell surrounding the central interacting galaxy pair. The origin of this phase-space correlated population of star-forming late-type satellite galaxies raises questions in view of the planes-of-satellites observed around more massive galaxy pairs that are, however, made up of old early-type dwarf galaxies. We also report the detection of a candidate stellar Plume near the binary. This elongated structure of low surface brightness is a likely optical counterpart to the HI-tail north of NGC 4490/85, recently discovered by the FAST radio telescope.

The recent Gaia Focused Product Release contains radial velocity time-series for more than 9,000 Gaia long-period photometric variables. Here we search for binary systems with large radial velocity amplitudes to identify candidates with massive, unseen companions. Eight targets have binary mass function $f(M)>1\ M_\odot$, three of which are eclipsing binaries. The remaining five show evidence of ellipsoidal modulations. We fit spectroscopic orbit models to the Gaia radial velocities, and fit the spectral energy distributions of three targets. For the two systems most likely to host dark companions, J0946 and J1640, we use PHOEBE to fit the ASAS-SN light curves and Gaia radial velocities. The derived companion masses are $>3 M_\odot$, but the high Galactic dust extinctions towards these objects limit our ability to rule out main sequence companions or subgiants hotter than the photometric primaries. These systems are similar to other stellar-mass black hole impostors, notably the Unicorn (V723 Mon) and the Giraffe (2M04123153$+$6738486). While it is possible that J1640 and J0946 are similar examples of stripped giant star binaries, high-resolution spectra can be used to determine the nature of their companions.

JWST provides an unprecedented opportunity for unbiased surveys of H$\alpha$-emitting galaxies at $z>4$ with the NIRCam wide-field slitless spectroscopy (WFSS). In this work, we present a census of Ly$\alpha$ escape fraction ($f_{esc, Ly\alpha}$) of 165 star-forming galaxies at $z=4.9-6.3$ using their H$\alpha$ emission directly measured from FRESCO NIRCam/WFSS data. We search for Ly$\alpha$ emission of each H$\alpha$-emitting galaxy in VLT/MUSE data. The overall $f_{esc, Ly\alpha}$ measured by stacking is $f_{esc, Ly\alpha}$ is $0.090\pm0.006$. We find that $f_{esc, Ly\alpha}$ displays a strong dependence on the observed UV slope ($\beta_{\rm obs}$) and E(B-V), such that the bluest galaxies ($\beta_{\rm obs}\sim-2.5$) have the largest escape fractions ($f_{\rm esc, Ly\alpha}\approx0.6$), indicative of the crucial role of dust and gas in modulating the escape of Ly$\alpha$ photons. $f_{esc, Ly\alpha}$ is less well related to other parameters, including the UV luminosity and stellar mass, and the variation in $f_{esc, Ly\alpha}$ with them can be explained by their underlying coupling with E(B-V) or $\beta_{\rm obs}$. Our results suggest a tentative decline in $f_{esc, Ly\alpha}$ at $z\gtrsim 5$, implying increasing intergalactic medium attenuation towards higher redshift. Furthermore, the dependence of $f_{esc, Ly\alpha}$ on $\beta_{\rm obs}$ is proportional to that of the ionizing photon escape fraction ($f_{\rm esc, LyC}$), indicating the escape of Ly$\alpha$ and ionizing photon may be regulated by similar physical processes. With $f_{esc, Ly\alpha}$ as a proxy to $f_{\rm esc, LyC}$, we infer that UV-faint ($M_{\rm UV}>-16$) galaxies contribute $>70\%$ of the total ionizing emissivity at $z=5-6$. If these relations hold during the epoch of reionization, UV-faint galaxies can contribute the majority of UV photon budget to reionize the Universe.

A number of stellar astrophysical phenomena, such as tidal novae and planetary engulfment, involve sudden injection of sub-binding energy in a thin layer within the star, leading to mass ejection of the stellar envelope. We use a 1D hydrodynamical model to survey the stellar response and mass loss for various amounts ($E_{\mathrm{dep}}$) and locations of the energy deposition. We find that the total mass ejection has a nontrivial dependence on $E_{\mathrm{dep}}$ due to the varying strengths of mass ejection events, which are associated with density/pressure waves breaking out from the stellar surface. The rapid occurrence of multiple breakouts may present a unique observational signature for sudden envelope heating events in stars.

We present the most precise constraints to date for the mass and age distributions of single ultracool dwarfs in the solar neighborhood, based on an updated volume-limited sample of 504 L, T, and Y dwarfs within 25 pc. We develop a Monte Carlo approach using the $\langle V/V_{\rm max}\rangle$ statistic to correct for incompleteness and obtain a space density of $(1.83_{-0.15}^{+0.16})\times10^{-2}$ pc$^{-3}$ for spectral types L0-Y2. We calculate bolometric luminosities for our sample, using an updated "super-magnitude" method for the faintest objects. We use our resulting luminosity function and a likelihood-based population synthesis approach to simultaneously constrain the mass and age distributions. We employ the fraction of young L0-L7 dwarfs as a novel input for this analysis that is crucial for constraining the age distribution. For a power-law mass function $\frac{dN}{dM} \propto M^{-\alpha}$ we find $\alpha=0.58_{-0.20}^{+0.16}$, indicating an increase in numbers toward lower masses, consistent with measurements in nearby star-forming regions. For an exponential age distribution $b(t) \propto e^{-\beta t}$ we find $\beta=-0.44\pm0.14$, i.e., a population with fewer old objects than often assumed, which may reflect dynamical heating of the Galactic plane as much as the historical brown dwarf birthrate. We compare our analysis to Kirkpatrick et al. (2021), who used a similar volume-limited sample. Although our mass function measurements are numerically consistent, their assumption of a flat age distribution is disfavored by our analysis, and we identify several important methodological differences between our two studies. Our calculation of the age distribution of solar neighborhood brown dwarfs is the first based on a volume-limited sample.

Galactic X-ray sources are diverse, ranging from active M dwarfs to compact object binaries, and everything in between. The X-ray landscape of today is rich, with point source catalogs such as those from XMM-Newton, Chandra, and Swift, each with $\gtrsim10^5$ sources and growing. Furthermore, X-ray astronomy is on the verge of being transformed through data releases from the all-sky SRG/eROSITA survey. Many X-ray sources can be associated with an optical counterpart, which in the era of Gaia, can be determined to be Galactic or extragalactic through parallax and proper motion information. Here, I present a simple diagram -- the ``X-ray Main Sequence", which distinguishes between compact objects and active stars based on their optical color and X-ray-to-optical flux ratio ($F_X/F_\textrm{opt}$). As a proof of concept, I present optical spectroscopy of six exotic accreting WDs discovered using the X-ray Main Sequence as applied to the XMM-Newton catalog. Looking ahead to surveys of the near future, I additionally present SDSS-V optical spectroscopy of new systems discovered using the X-ray Main Sequence as applied to the SRG/eROSITA eFEDS catalog.

The dynamics of the Local Group (LG), especially concerning the contributions of the Milky Way (MW) and Andromeda (M31) galaxies, is sensitive to the presence of dark energy. This work compares the evolution of the LG by considering it as a two-body problem in a homogeneous and isotropic expanding spacetime, i.e. the McVitte spacetime (McV) versus the spherically symmetric metric for LG dynamics with the Cosmological Constant, i.e. the De Sitter-Schwarzschild spacetime (DsS). Using the Timing Argument (which links LG dynamics to LG mass), calibrated by the IllustrisTNG simulations, we find that the McV spacetime predicts a lower mass for the LG: $\left(4.20 \pm 0.61\right) \cdot 10^{12} M_{\odot}$ for McV spacetime vs. $\left(4.65 \pm 0.75\right) \cdot 10^{12} M_{\odot}$ for DsS spacetime ($68 \% ,$ CL). Due to uncertainties in tangential velocity measurements, the masses are indistinguishable. However, with future astrometric measurements, we demonstrate that the predicted masses will be distinguishable, indicating different LG histories. By independently estimating the total mass of MW and M31, we compare the possible upper bounds for the Cosmological Constant in these scenarios. We find a tighter upper bound for the DsS spacetime model, $\Lambda < 3.3 \,\Lambda_{\text{CMB}}$, compared to $\Lambda < 8.4\, \Lambda_{\text{CMB}}$ for the McV spacetime (where $\Lambda_{\text{CMB}}$ is the mean value from Planck). Future astrometric measurements, such as those from JWST, hold the potential to independently detect dark energy for both spacetime models independent from Planck's value.

Abridged: We measured photometric and spectroscopic $P_{\rm rot}$ for a large sample of nearby bright M dwarfs with spectral types from M0 to M9, as part of our continual effort to fully characterize the Guaranteed Time Observation programme stars of the CARMENES survey. We determine $P_{\rm rot}$ for 129 stars. Combined with the literature, we tabulate $P_{\rm rot}$ for 261 stars, or 75% of our sample. We evaluate the plausibility of all periods available for this sample by comparing them with activity signatures and checking for consistency between multiple measurements. We find that 166 of these stars have independent evidence that confirmed their $P_{\rm rot}$. There are inconsistencies in 27 periods, which we classify as debated. A further 68 periods are identified as provisional detections that could benefit from independent verification. We provide an empirical relation for the $P_{\rm rot}$ uncertainty as a function of the $P_{\rm rot}$ value, based on the dispersion of the measurements. We show that published formal errors seem to be often underestimated for periods $\gtrsim 10$ d. We highlight the importance of independent verification on $P_{\rm rot}$ measurements, especially for inactive M dwarfs. We examine rotation-activity relations with emission in X-rays, H$\alpha$, Ca II H & K, and surface magnetic field strengths. We find overall agreement with previous works, as well as tentative differences in the partially versus fully convective subsamples. We show $P_{\rm rot}$ as a function of stellar mass, age, and galactic kinematics. With the notable exception of three transiting planet systems and TZ Ari, all known planet hosts in this sample have $P_{\rm rot} \gtrsim 15$ d. This indicates that important limitations need to be overcome before the radial velocity technique can be routinely used to detect and study planets around young and active stars.

We consider small-scale jet-like events that might make the solar wind, as has been suggested in recent studies. We show that the events referred to as "coronal jets" and as "jetlets" both fall on a power-law distribution that also includes large-scale eruptions and spicule-sized features; all of the jet-like events could contribute to the solar wind. Based on imaging and magnetic field data, it is plausible that many or most of these events might form by the same mechanism: Magnetic flux cancelation produces small-scale flux ropes, often containing a cool-material minifilament. This minifilament/flux rope erupts and reconnects with adjacent open coronal field, along which "plasma jets" flow and contribute to the solar wind. The erupting flux ropes can contain twist that is transferred to the open field, and these become Alfv\'enic pulses that form magnetic switchbacks, providing an intrinsic connection between switchbacks and the production of the solar wind.

We present new analytical solutions for the evolution of protoplanetary discs (PPDs) where magnetohydrodynamic (MHD) wind-driven processes dominate. Our study uses a 1D model which incorporates equations detailing angular momentum extraction by MHD winds and mass-loss rates. Our solutions demonstrate that the disc retains its initial state during the early phases; however, it rapidly evolves towards a self-similar state in the later stages of disc evolution. The total disc mass undergoes a continuous decline over time, with a particularly rapid reduction occurring beyond a certain critical time threshold. This gradual decrease in mass is influenced by the wind parameters and the initial surface density of the disc. In the MHD wind-dominated regime, we show that the disc's lifespan correlates positively with the magnetic lever arm up to a certain threshold, irrespective of the initial disc size. PPDs with a larger magnetic lever arm are found to maintain significantly higher total disc mass over extended periods compared to their counterparts. The mass ejection-to-accretion ratio increases in efficient wind scenarios and is further amplified by a steeper initial surface density profile. Our analysis also reveals varied evolutionary trajectories in the plane of accretion rate and total disc mass, influenced by magnetic parameters and initial disc size. In scenarios with efficient MHD winds, discs with bigger sizes have extended operation time for mechanisms governing planet formation.

Accreting main-sequence stars expand significantly when the mass accretion timescale is much shorter than their thermal timescales. This occurs during mass transfer from an evolved giant star onto a main-sequence companion in a binary system, and is an important phase in the formation of compact binaries including X-ray binaries, cataclysmic variables, and gravitational-wave sources. In this study, we compute 1D stellar models of main-sequence accretors with different initial masses and accretion rates. The calculations are used to derive semi-analytical approximations to the maximum expansion radius. We assume that mass transfer remains fully conservative as long as the inflated accretor fits within its Roche lobe, leading stars to behave like hamsters, stuffing excess material behind their expanding cheeks. We suggest a physically motivated prescription for the mass growth of such "hamstars", which can be used to determine mass-transfer efficiency in rapid binary population synthesis models. With this prescription, we estimate that progenitors of high-mass X-ray binaries and gravitational-wave sources may have experienced highly non-conservative mass transfer. In contrast, for low-mass accretors, the accretion timescale can exceed the thermal timescale by a larger factor without causing significant radial expansion.

Post-starburst galaxies are believed to be in a rapid transition between major merger starbursts and quiescent ellipticals, where AGN feedback is suggested as one of the processes responsible for the quenching. To study the role of AGN feedback, we constructed a sample of post-starburst candidates with AGN and indications of ionized outflows. We use MUSE/VLT observations to resolve the properties of the stars and multi-phased gas in five of them. All the galaxies show signatures of interaction/merger in their stellar or gas properties, with some galaxies at an early stage of interaction with companions at distances $\sim$50 kpc, suggesting that optical post-starburst signatures may be present well before the final starburst and coalescence. We detect narrow and broad kinematic components in multiple transitions in all the galaxies. Our detailed analysis of their kinematics and morphology suggests that, contrary to our expectation, the properties of the broad kinematic components are inconsistent with AGN-driven winds in 3 out of 5 galaxies. The two exceptions are also the only galaxies in which spatially-resolved NaID P-Cygni profiles are detected. In some cases, the observations are more consistent with interaction-induced galactic-scale flows, an often overlooked process. These observations raise the question of how to interpret broad kinematic components in interacting and perhaps also in active galaxies, in particular when spatially-resolved observations are not available or cannot rule out merger-induced galactic-scale motions. We suggest that NaID P-Cygni profiles are more effective outflow tracers, and use them to estimate the energy that is carried by the outflow.

Here we provide a review of Mira variables, their basic properties, and Period-Luminosity Relations with an emphasis on their role in measuring the Hubble Constant. The usage of multiple independent distance indicators and methods is crucial to cross-checking systematic uncertainties in distance measurements and in reinforcing previous findings of the Hubble tension. To this end, Mira variables serve as an alternative Type Ia Supernova calibrator to the more commonly-used Cepheid variables or Tip of the Red Giant Branch method. They also have the potential to expand the number of local SN Ia calibrators by calibrating previously-inaccessible SNe Ia. Short-period ($P \lesssim 400$ d) O-rich Miras are a ubiquitous older population that can reach galaxies not hosting the younger Cepheids variables or out of reach to the old but fainter Tip of the Red Giant Branch. With the current and upcoming focus on infrared observations, Miras, which can be discovered and characterized using exclusively near-infrared and infrared observations, will be particularly useful in obtaining distances to astrophysical objects. Long-period Miras ($P \gtrsim 400$ d) are highly luminous variables that have the potential to measure $H_0$ directly, excluding Type Ia SNe altogether in the distance ladder.

Proxima b is a rocky exoplanet in the habitable zone of the nearest star system and a key test case in the search for extraterrestrial life. Here, we investigate the characterization of a potential Earth-like atmosphere around Proxima b in reflected light via molecule mapping, combining high resolution spectroscopy (HRS) and high contrast imaging, using the first-generation integral field spectrograph HARMONI on the $39$-m Extremely Large Telescope. We simulate comprehensive observations of Proxima b at an assumed $45^{\circ}$ inclination using HARMONI's High Contrast Adaptive Optics mode, with spatial resolution $\sim 8$mas ($3.88$mas/spaxel) and spectral resolving power $R\simeq17,000$ between $1.538$--$1.678 \mu m$, containing the spectral features of water, carbon dioxide and methane. Tellurics, stellar features, and additional noise sources are included, and removed using established molecule mapping techniques. We find that HARMONI's current focal plane mask (FPM) is too large and obscures the orbit of Proxima b and thus explore smaller and offset FPMs to yield a detection. A $\rm{S/N}=5$ detection of Proxima b's reflected light, suitable for atmospheric characterisation, is possible with such modifications, requiring a minimum of $20$ hours, but ideally at least $30$ hours of integration time. We highlight that such detections do not scale with the photon noise, hence suitably detailed simulations of future instruments for the ELTs are needed to fully understand their ability to perform HRS observations of exoplanet atmospheres. Alterations to the HARMONI FPM design are feasible at this stage, but must be considered in context of other science cases.

It has recently been shown that local primordial non-Gaussianities (PNG) with significant amplitude ($|f_{\rm NL}| \sim 1000$), at small (Mpc) scales, can help in forming simulated galaxies with more disky baryonic kinematics than in the Gaussian case, while generating matter power spectra that can differ by up to 20% from the Gaussian case at non-linear scales. Here, we explore in detail the consequences of such small-scale PNG on the dark matter halo profiles. We show in particular that, for negative $f_{\rm NL}$, dark matter halos formed in collisionless simulations are not always well described by the traditional Navarro-Frenk-White (NFW) profiles, as supported by their sparsity distribution. We conclude that NFW profiles are not as clear attractors for the density profiles of dark matter halos in the presence of PNG than in the case of a Gaussian contrast density field. We show how a minimal extension of the NFW profile can describe halos both in the Gaussian and non-Gaussian cases. From the combination of our sparsity analysis and the quality of the adjustments of the density profiles with a minimal extension to NFW, we conclude that $z=1$ halos carry the most interesting information about PNG.

Supermassive black holes require a reservoir of cold gas at the centre of their host galaxy in order to accrete and shine as active galactic nuclei (AGN). Major mergers have the ability to drive gas rapidly inwards, but observations trying to link mergers with AGN have found mixed results due to the difficulty of consistently identifying galaxy mergers in surveys. This study applies deep learning to this problem, using convolutional neural networks trained to identify simulated post-merger galaxies from survey-realistic imaging. This provides a fast and repeatable alternative to human visual inspection. Using this tool, we examine a sample of ~8500 Seyfert 2 galaxies (L[OIII] ~ $10^{38.5 - 42}$ erg/s) at z < 0.3 in the Sloan Digital Sky Survey and find a merger fraction of $2.19_{-0.17}^{+0.21}$% compared with inactive control galaxies, in which we find a merger fraction of $2.96_{-0.20}^{+0.26}$%, indicating an overall lack of mergers among AGN hosts compared with controls. However, matching the controls to the AGN hosts in stellar mass and star formation rate reveals that AGN hosts in the star-forming blue cloud exhibit a ~$2\times$ merger enhancement over controls, while those in the quiescent red sequence have significantly lower relative merger fractions, leading to the observed overall deficit due to the differing $M_{\ast} - $SFR distributions. We conclude that while mergers are not the dominant trigger of all low-luminosity, obscured AGN activity in the nearby Universe, they are more important to AGN fuelling in galaxies with higher cold gas mass fractions as traced through star formation.

Context: Asteroseismology, a powerful approach for obtaining internal structure and stellar properties, requires surface temperature and chemical composition information to determine mass and age. High-resolution spectroscopy is a valuable technique for precise stellar parameters (including surface temperature) and chemical composition analysis. Aim: We aim to combine spectroscopic parameters with asteroseismology to test stellar models. Method: Using high-resolution optical and near-IR spectra from GIARPS at the Telescopio Nazionale Galileo, we conducted a detailed spectroscopic analysis of 16 stars photometrically selected to be on the red giant and red clump. Stellar parameters and chemical abundances for light elements (Li, C, N, F), Fe peak, $\alpha$ and n-capture elements were derived using a combination of equivalent widths and spectral synthesis techniques, based on atomic and molecular features. Ages were determined through asteroseismic scaling relations and compared with ages based on chemical clocks, [Y/Mg] and [C/N]. Results: Spectroscopic parameters confirmed the stars as part of the red giant branch and red clump. Two objects, HD 22045 and HD 24680 exhibited relatively high Li abundances, with HD 24680 potentially being a Li-rich giant resulting from mass transfer with an intermediate-mass companion, which already underwent its AGB phase. Stellar parameters derived from scaling different sets of relations were consistent with each other. For what concerns ages, the values based on asteroseismology were in excellent agreement with those derived from theoretical evolutionary tracks, but did not align with ages derived from the chemical clocks [Y/Mg] and [C/N].

We present the spectral analysis of bright steady states in an outburst of the transient neutron star low-mass X-ray binary (NS-LMXB) Swift J1858.6-0814 observed with NICER. We detected an ionized iron K absorption line (H-like Fe) at 6.97keV in the spectrum. We estimated the photoionization parameter using the ratio of the equivalent widths (EWs) of the FeXXVI (H-like) (17+\-5eV) and FeXXV (He-like) (<3eV) and discuss the origin of the iron absorption line. The irradiated gas producing the absorption line would locate within (3-6)*1E9cm from the X-ray source. We suggest that the observed H-like Fe absorption line originates from the highly-ionized gas in the inner accretion disk in Swift J1858.6-0814.

The NASA Astrophysics Data System (ADS) is the primary Digital Library portal for researchers in astronomy and astrophysics. Over the past 30 years, the ADS has gone from being an astronomy-focused bibliographic database to an open digital library system supporting research in space and (soon) earth sciences. This paper describes the evolution of the ADS system, its capabilities, and the technological infrastructure underpinning it. We give an overview of the ADS's original architecture, constructed primarily around simple database models. This bespoke system allowed for the efficient indexing of metadata and citations, the digitization and archival of full-text articles, and the rapid development of discipline-specific capabilities running on commodity hardware. The move towards a cloud-based microservices architecture and an open-source search engine in the late 2010s marked a significant shift, bringing full-text search capabilities, a modern API, higher uptime, more reliable data retrieval, and integration of advanced visualizations and analytics. Another crucial evolution came with the gradual and ongoing incorporation of Machine Learning and Natural Language Processing algorithms in our data pipelines. Originally used for information extraction and classification tasks, NLP and ML techniques are now being developed to improve metadata enrichment, search, notifications, and recommendations. we describe how these computational techniques are being embedded into our software infrastructure, the challenges faced, and the benefits reaped. Finally, we conclude by describing the future prospects of ADS and its ongoing expansion, discussing the challenges of managing an interdisciplinary information system in the era of AI and Open Science, where information is abundant, technology is transformative, but their trustworthiness can be elusive.

We present the results of an ultra-deep radio continuum survey, containing $\sim480$ hours of observations, of the Galactic globular cluster 47 Tucanae with the Australia Telescope Compact Array. This comprehensive coverage of the cluster allows us to reach RMS noise levels of 1.19 $\mu Jy~\textrm{beam}^{-1}$ at 5.5 GHz, 940 $nJy~\textrm{beam}^{-1}$ at 9 GHz, and 790 $nJy~\textrm{beam}^{-1}$ in a stacked 7.25 GHz image. This is the deepest radio image of a globular cluster, and the deepest image ever made with the Australia Telescope Compact Array. We identify ATCA J002405.702-720452.361, a faint ($6.3\pm1.2$ $\mu Jy$ at 5.5 GHz, $5.4\pm0.9$ $\mu Jy$ at 9 GHz), flat-spectrum ($\alpha=-0.31\pm0.54$) radio source that is positionally coincident with the cluster centre and potentially associated with a faint X-ray source. No convincing optical counterpart was identified. We use radio, X-ray, optical, and UV data to show that explanations involving a background active galactic nucleus, a chromospherically active binary, or a binary involving a white dwarf are unlikely. The most plausible explanations are that the source is an undiscovered millisecond pulsar or a weakly accreting black hole. If the X-ray source is associated with the radio source, the fundamental plane of black hole activity suggests a black hole mass of $\sim54-6000$ M$_{\odot}$, indicating an intermediate-mass black hole or a heavy stellar-mass black hole.

A new mechanism is proposed to account for the formation of retrograde hot Jupiter in coplanar star-planet system via close encounter between a Jupiter mass planet and a brown dwarf mass planet. After long timescale scattering between several Jupiter mass planets with inner orbits, the remaining planets still rotating around the star could have large semimajor with large eccentricity. If there exists a brown dwarf mass planet in distant orbit around the star, planetary encounter may happen. After encounter, the Jupiter mass planet may rotate around the star in a retrograde orbit with extremely large eccentricity and the periastron can reach 0.01 AU, which means that, within the first several orbits around the star, tidal interaction from the star can shrink the semimajor axis of the planet quickly. Thus, the Jupiter mass planet is isolated from the brown dwarf mass planet due to the quick decrease of its apastron distance and eventually evolves into a retrograde hot Jupiter.

HD 54879 is the most recently discovered magnetic O-type star. Previous studies ruled out a rotation period shorter than 7 years, implying that HD 54879 is the second most slowly-rotating known magnetic O-type star. We report new high-resolution spectropolarimetric measurements of HD 54879, which confirm that a full stellar rotation cycle has been observed, and derive a rotation period of $P = 2562^{+63}_{-58}$ d (about 7.02 yr). The radial velocity of HD 54879 has been stable over the last decade of observations. We explore equivalent widths and longitudinal magnetic fields calculated from lines of different elements, and conclude the atmosphere of HD 54879 is likely chemically homogeneous, with no strong evidence for chemical stratification or lateral abundance nonuniformities. We present the first detailed magnetic map of the star, with an average surface magnetic field strength of 2954 G, and a strength for the dipole component of 3939 G. There is a significant amount of magnetic energy in the quadrupole components of the field (23 percent). Thus, we find HD 54879 has a strong magnetic field with a significantly complex topology.

The Galactic diffuse X-ray emission (GDXE) can be spatially segmented into the Galactic center X-ray Emission (GCXE), the Galactic ridge X-ray emission (GRXE), and the Galactic bulge X-ray emission (GBXE). The X-ray spectra of the GDXE are expressed by the assembly of compact X-ray sources, which are either the white dwarfs (WDs), or the X-ray active stars, consisting of binaries with late type stars. The WDs have either strong magnetic field (mCV), or weak magnetic field (non-mCV). The WDs and X-ray active stars are collectively called as compact X-ray stars. However, spectral fittings by the assembly of all compact X-ray stars for the GCXE, GRXE, and GBXE are rejected leaving significant excess near the energies of K$\alpha$, He$\alpha$, Ly$\alpha$ lines. These excesses are found in the collisional ionization equilibrium (CIE) plasma. Thus the spectra of the GRXE and GBXE are improved by adding CIE-SNRs. However the GCXE spectrum is still unacceptable with significant data excess due to the radiative recombination emission (RP-plasma). Then the GCXE fit is significantly improved by adding aged RP-SNRs. The aged RP-SNRs would be made by a past big flare of Sgr~A$^*$ emitting either hard X-rays or low-energy cosmic-rays. The big flares may excite Fe and Ni atoms in cold diffuse gas (CM), and emit fluorescent X-ray lines. The CIE-SNRs, RP-SNRs and CM are called as diffuse X-ray sources. This paper presents the spectral fits by the assembly of all the compact and diffuse X-ray sources together with high quality spectra and combined fit among all GDXE of GCXE, GRXE, and GBXE. This provides the first scenario to quantitatively and comprehensively predict the origin of the GDXE spectra.

We investigate an SNR W28 with the Suzaku archive data and report the results of spatial resolved analyses. We carry out spectral analysis using a recombining plasma (RP) model with an element-dependent initial ionization temperature, and obtain the ionization temperatures to be $\sim0.5$~keV for Ne, $\sim0.7$~keV for Mg, $\sim1.0$~keV for Si, $\sim1.2$~keV for S, $\sim1.4$~keV for Ar, $\sim1.7$~keV for Ca, and $\sim0.7$~keV for Fe in the RP-initial phase. In addition to northeast regions where RP have been reported, we find that the ionization temperature in the southeast and southwest regions show a similar trend to the central region, in the RP-initial phase. Furthermore, the elapsed time from the RP-initial phase to present is shorter, $\sim300$~yr in the central region and longer, $\sim10^3$-$10^4$~yr in the outside regions. Our results cannot be explained by simple scenarios of thermal conduction due to molecular clouds or adiabatic cooling (rarefaction), and indicate that more complex mechanism or other scenarios are required. Also, we estimate the ejecta mass $\gtrsim14M_{\odot}$, which indicates a SNR derived a massive star.

The low-energy $\gamma$-ray (0.1-30 MeV) sky has been relatively unexplored since the decommissioning of the COMPTEL instrument on the Compton Gamma-Ray Observatory (CGRO) satellite in 2000. However, the study of this part of the energy spectrum (the ``MeV gap") is crucial for addressing numerous unresolved questions in high-energy and multi-messenger astrophysics. Although several large MeV $\gamma$-ray missions like AMEGO and e-ASTROGAM are being proposed, they are predominantly in the developmental phase, with launches not anticipated until the next decade at the earliest. In recent times, there has been a surge in proposed CubeSat missions as cost-effective and rapidly implementable ``pathfinder" alternatives. A MeV CubeSat dedicated to $\gamma$-ray astronomy has the potential to serve as a demonstrator for future, larger-scale MeV payloads. This paper presents a $\gamma$-ray payload design featuring a CdZnTe crystal calorimeter module developed by IDEAS. We report the detailed results of simulations to assess the performance of this proposed payload and compare it with those of previous $\gamma$-ray instruments.

We explore the properties of an 'almost' dark cloud of neutral hydrogen (HI) using data from the Widefield ASKAP L-band Legacy All-sky Survey (WALLABY). Until recently, WALLABY J103508-283427 (also known as H1032-2819 or LEDA 2793457) was not known to have an optical counterpart, but we have identified an extremely faint optical counterpart in the DESI Legacy Imaging Survey Data Release 10. We measured the mean g-band surface brightness to be $27.0\pm0.3$ mag arcsec$^{-2}$. The WALLABY data revealed the cloud to be closely associated with the interacting group Klemola 13 (also known as HIPASS J1034-28 and the Tol 9 group), which itself is associated with the Hydra cluster. In addition to WALLABY J103508-283427/H1032-2819, Klemola 13 contains ten known significant galaxies and almost half of the total HI gas is beyond the optical limits of the galaxies. By combining the new WALLABY data with archival data from the Australia Telescope Compact Array (ATCA), we investigate the HI distribution and kinematics of the system. We discuss the relative role of tidal interactions and ram pressure stripping in the formation of the cloud and the evolution of the system. The ease of detection of this cloud and intragroup gas is due to the sensitivity, resolution and wide field of view of WALLABY, and showcases the potential of the full WALLABY survey to detect many more examples.

Recent numerical studies suggest that magnetic fields play an important role in primordial star formation in the early universe. However, the detailed evolution of the magnetic field in the collapse phase still has uncertainties because of the complicated physics associated with turbulence in a collapsing magnetized system. Here, we perform a suite of numerical MHD simulations that follow the collapse of magnetized, turbulent primordial gas clouds to investigate the evolution of the magnetic field associated with the turbulence, assuming a polytropic equation of state with exponent $\gamma_{\rm eff}$ and with various numerical resolutions. In addition, we generalize the analytic theory of magnetic field growth/saturation so that it can deal with various exponents $\gamma_{\rm eff}$ and turbulence energy spectra. We find that the numerical results are well reproduced by the theory for various $\gamma_{\rm eff}$ through the collapse phase during the formation of the first stars. The magnetic field is eventually amplified by a factor of $10^{12}$ -- $10^{15}$ due to kinematic and non-linear turbulent dynamo effects and reaches 3% -- 100% of the equipartition level, depending on $\gamma_{\rm eff}$. We also find that the transition between the kinematic and non-linear stages can be analytically estimated. These results indicate that the strong magnetic field accompanied by supersonic turbulence is a general property and suggest that it can play a crucial role in the formation of the first stars.

We report the discovery of 9 new hot molecular cores in the Deep South (DS) region of Sagittarius B2 using Atacama Large Millimeter/submillimeter Array Band 6 observations. We measure the rotational temperature of CH$_3$OH and derive the physical conditions present within these cores and the hot core Sgr B2(S). The cores show heterogeneous temperature structure, with peak temperatures between 252 and 662 K. We find that the cores span a range of masses (203-4842 M$_\odot$) and radii (3587-9436 AU). CH$_3$OH abundances consistently increase with temperature across the sample. Our measurements show the DS hot cores are structurally similar to Galactic Disk hot cores, with radii and temperature gradients that are comparable to sources in the Disk. They also show shallower density gradients than Disk hot cores, which may arise from the Central Molecular Zone's higher density threshold for star formation. The hot cores have properties which are consistent with those of Sgr B2(N), with 3 associated with Class II CH$_3$OH masers and one associated with an UCHII region. Our sample nearly doubles the high-mass star forming gas mass near Sgr B2(S) and suggest the region may be a younger, comparably massive counterpart to Sgr B2(N) and (M). The relationship between peak CH$_3$OH abundance and rotational temperature traced by our sample and a selection of comparable hot cores is qualitatively consistent with predictions from chemical modeling. However, we observe constant peak abundances at higher temperatures ($T \gtrsim 250$ K), which may indicate mechanisms for methanol survival that are not yet accounted for in models.

{We explored the {softness parameter} in the infrared, whose main purpose is the characterisation of the hardness of the incident ionising radiation in emission-line nebulae. This parameter is obtained from the combination of mid-infrared wavelength range transitions corresponding to consecutive ionisation stages in star-forming regions. We compiled observational data from a sample of star-forming galaxies (SFGs), including luminous and ultraluminous infrared galaxies (LIRGs and ULIRGs), to study the softness parameter and its equivalent expression in two dimensions, the softness diagram. We compared them with predictions from photoionisation models to determine the shape of the ionising continuum energy distribution in each case. We also used the measured emission-line ratios as input for HCmistry-Teff-IR, a code that performs a Bayesian-like comparison with photoionisation model predictions in order to quantify the equivalent effective temperature (T*) and the ionisation parameter. We found similar average values within the errors of the softness parameter in (U)LIRGs (-0.57) in the rest of the SFGs (-0.51), which could be interpreted as indicative of a similar incident radiation field. This result is confirmed from the analysis using HCm-Teff-IR, which simultaneously points to a slightly lower, although similar within the errors, T* scale for (U)LIRGs, even when a higher dust-to-gas mass ratio is considered in the models for these objects. These derived T* values are compatible with the ionisation from massive stars, without any need of harder ionising sources, both for (U)LIRGs and the rest of the SFGs. However, the derived T* in (U)LIRGs do not show any correlation with metallicity. This could be interpreted as a sign that their similar average T* values are due to the attenuation of the energetic incident flux from massive stars by the heated dust mixed with the gas.

Context. Slow waves in solar coronal loops are strongly damped. The current theory of damping by thermal conduction cannot explain some observational features.\n Aims. We investigate the propagation of slow waves in a coronal loop built up from strands of different temperatures. \n Methods. We consider the loop to have a multithermal, Gaussian temperature distribution. The different propagation speeds in different strands lead to an multithermal apparent damping of the wave, similar to observational phase mixing. We use an analytical model to predict the damping length and propagation speed for the slow waves, including in imaging with filter telescopes. \n Results. We compare the damping length due to this multithermal apparent damping with damping due to thermal conduction and find that the multithermal apparent damping is more important for shorter period slow waves. We have found the influence of instrument filters on the wave's propagation speed and damping. This allows us to compare our analytical theory to forward models of numerical simulations. \n Conclusions. We find that our analytical model matches the numerical simulations very well. Moreover, we offer an outlook for using the slow wave properties to infer the loop's thermal properties.

Optical broad-band spectral shape measurements of gamma-ray bursts (GRBs) are typically made starting an hour or more after the trigger event. With our automated, rapid-response system, the Burst Simultaneous Three-channel Imager (BSTI) on the Nazarbayev University Transient Telescope at Assy-Turgen Astrophysical Observatory (NUTTelA-TAO), we began measurements of GRB200925B 129 s after the Swift BAT trigger. The temporal decay log slopes in the g', r', and i' bands in the time interval 129 s to 1029 s are -0.43 \pm 0.31, -0.43 \pm 0.15, and -0.72 \pm 0.14, respectively. During the decay phase, a shift in color from red to blue, a change in log slope of \{beta} from -2.73 to -1.52 was measured. The evolution in the optical spectral slope is consistent with a decrease in extinction caused by dust destruction.

Future research studies of cosmic microwave background polarization seems likely to provide a more improved upper bound of $r \le 0.03$ on the tensor-to-scalar ratio(r). In our work, we have done the reheating study of mutated hilltop inflation(MHI), a model falling in the broad category of small field inflation. We have parameterized reheating in terms of various parameters like reheating duration $N_{\text{rh}}$, reheating temperature $T_{\text{rh}}$ and effective equation of state $\overline{\omega }_{\text{rh}}$ using observationally viable values of scalar power spectrum amplitude $A_{\text{s}}$ and scalar spectral index $n_{\text{s}}$. In our study, working over a range of $\overline{\omega }_{\text{rh}}$, we found that the MHI potential is well consistent with combined Planck and BK18 observations for $\overline{\omega }_{\text{rh}} > 0$ within a particular range of model's parameter space and the lower values of the model parameter in MHI generate considerably smaller r compared to normal hilltop potential without any incompatibility of $n_s$ with observational data, making MHI a better choice in accordance to recent and future studies.

The initial-final mass relation (IFMR) plays a crucial role in understanding stellar structure and evolution by linking a star's initial mass to the mass of the resulting white dwarf. This study explores the IFMR in the initial mass range $0.8 \leq M_\mathrm{ini} / M_\odot \leq 4$ using full PARSEC evolutionary calculations supplemented with COLIBRI computations to complete the ejection of the envelope and obtain the final core mass. Recent works have shown that the supposed monotonicity of the IFMR is interrupted by a kink in the initial mass range $M_\mathrm{ini} \approx 1.65-2.10 M_\odot$, due to the interaction between recurrent dredge-up episodes and stellar winds in carbon stars evolving on the thermally-pulsing asymptotic giant branch phase. To reproduce the IFMR non-monotonic behavior we investigate the role of convective overshooting efficiency applied to the base of the convective envelope ($f_\mathrm{env}$) and to the borders of the pulse-driven convective zone ($f_\mathrm{pdcz}$), as well as its interplay with mass loss. We compare our models to observational data and find that $f_\mathrm{env}$ must vary with initial mass in order to accurately reproduce the IFMR's observed kink and slopes. We find some degeneracy between the overshooting parameters when only the IFMR information is used. Nonetheless, this analysis provides valuable insights into the internal mixing processes during the TP-AGB phase.

The temperature of the interstellar medium (ISM) is governed by several physical process, among which radiative cooling, external UV/cosmic ray heating, and the mechanical work by compression and expansion. In regimes where the dynamical effect is important, the temperature deviates from that derived by simply balancing the heating and cooling functions. This renders the expression of the gas energy evolution with a simple equation of state (EOS) less straightforward. Given a cooling function, the behavior of the gas is subject to the combined effect of dynamical compression and radiative cooling. The goal of the present work is to derive the effective EOS of a collapsing gas within a full fluid solution. We solve the Navier-Stokes equations with a parametric cooling term in spherical coordinate and look for a self-similar collapse solution. We present a solution which describes a cloud that is contracting while losing energy through radiation. This yields an effective EOS that can be generally applied to various ISM context, where the cooling function is available from first principles, and expressed as powerlaw product of the density and temperature. Our findings suggest that a radiatively cooling gas under self-gravitating collapse can easily manifest an effective polytropic EOS, even isothermal in many scenarios. The present model provides theoretical justification for the simplifying isothermal assumptions of simulations at various scales, and can also provide a more realistic thermal recipe without additional computation cost.

The formation of $e^\pm$ plasmas within pulsar magnetospheres through quantum electrodynamics (QED) cascades in vacuum gaps is widely acknowledged. This paper aims to investigate the effect of photon polarization during the QED cascade occurring over the polar cap of a pulsar. We employ a Monte Carlo-based QED algorithm that accurately accounts for both spin and polarization effects during photon emission and pair production in both single-particle and particle-in-cell (PIC) simulations. Our findings reveal distinctive properties in the photon polarization of curvature radiation (CR) and synchrotron radiation (SR). CR photons exhibit high linear polarization parallel to the plane of the curved magnetic field lines, whereas SR photons, on average, demonstrate weak polarization. As the QED cascade progresses, SR photons gradually dominate over CR photons, thus reducing the average degree of photon polarization. Additionally, our study highlights an intriguing observation: the polarization of CR photons enhances $e^\pm$ pair production by approximately 5%, in contrast to the inhibition observed in laser-plasma interactions. Our self-consistent QED PIC simulations in the corotating frame reproduce the essential results obtained from single-particle simulations.

In a born-again planetary nebula (PN), processed H-deficient material has been injected inside the old, H-rich nebula as a result of a very late thermal pulse (VLTP) event. Long-slit spectra have been used to unveil the chemical and physical differences between these two structures, but the ejection and shaping processes remain still unclear. In order to peer into the morpho-kinematics of the H-deficient ejecta in the born-again PN A 58, we present the first integral field spectroscopic observations of a born-again PN as obtained with GTC MEGARA. We detect emission from the H$\alpha$, He I, [O III], [N II] and [S II] emission lines, which help us unveil the expansion patterns of the different structures. In combination with ALMA and Hubble Space Telescope data we are able to produce a complete view of the H-deficient ionized and molecular ejecta in A 58. We propose an hourglass structure for the ionized material that embraces molecular high-velocity polar components, while bisected by an expanding toroidal molecular and dusty structure. Our results leverage the role of a companion in shaping the VLTP ejecta in this born-again PN.

Interstellar aromatic molecules such as polycyclic aromatic hydrocarbons and polycyclic nitrogen and oxygen bearing molecules are thought to be abundant in the interstellar medium. In this class of molecules, benzonitrile ($c$-C$_6$H$_5$CN) plays an important role as a proxy for benzene. It has been detected through rotational emission in several astrophysical sources and is one of the simplest N-bearing polar aromatic molecules. Even in the cold ISM, the population of the rotational levels of benzonitrile might not be at equilibrium. Consequently, modeling its detected emission lines requires a prior computation of its quenching rate coefficients by the most abundant species in the ISM (He or H$_2$). In this paper, we focus on the excitation of c-C$_6$H$_5$CN by collision with He. We compute the first potential energy surface (PES) using the explicitly correlated coupled cluster method in conjunction with large basis sets. The PES obtained is characterized by a potential well depth of -97.2 cm$^{-1}$ and an important anisotropy. Scattering computations of the rotational (de-)excitation of c-C$_6$H$_5$CN by He atoms are performed by means of the coupled states approximation that allow to obtain collisional rates for rotational states up to $j$ = 9 and temperatures up to 40 K. These rate coefficients are then used to examine the effect of C$_6$H$_5$CN excitation induced by collisions with para-H$_2$ in molecular clouds by carrying out simple radiative transfer calculations of the excitation temperatures and show that non-equilibrium effects can be expected for H$_2$ densities up to 10$^5$-10$^6$ cm$^{-3}$.

Several planetary nebulae (PNe) have been found to contain both polycyclic aromatic hydrocarbon (PAH-like) species and fullerenes (C$_{60}$) distinguished by their mid-infrared emission. Previous laboratory and astronomical studies suggest that the formation of both species could be related to the decomposition, by photochemical processing, of hydrogenated amorphous carbon (HAC) grains. Then, HACs and, seemingly, big-fullerene related species (e.g., carbon onions) have been suggested as potential carriers of the UV bump at 2175{\AA} and the far-UV rise common to interstellar extinction curves. Our goal is to investigate the UV bump with the possible presence of a HAC extinction component in the International Ultraviolet Explorer (IUE) spectra of C-rich PNe; both with detected and non-detected fullerenes. The considered sample includes three C$_{60}$-PNe (Tc 1, IC 418, and IC 2501) and the non-C$_{60}$-PN Hen 2-5. Independently of the presence of C$_{60}$ in their circumstellar envelopes, we found that the UV bump in all sample PNe is well explained by interstellar extinction, suggesting that species different from those of the foreground insterstellar medium, e.g., large fullerene-related species like carbon onions, are not the carrier. Interestingly, we found that PNe Tc 1 and Hen 2-5 show an absorption in the FUV rise. Their IUE continuum spectra may be very well reproduced by including the extinction curve of HAC-like very small grains (VSG). The possible presence of both species, HAC-like grains and fullerenes (C$_{60}$), in Tc 1 could support the HAC photochemical processing scenario for the formation of fullerenes in the complex circumstellar envelopes of PNe.

Among the compact objects observed in gravitational wave merger events a few have masses in the gap between the most massive neutron stars (NSs) and least massive black holes (BHs) known. Their nature and the formation of their merging binaries are not well understood. We report on pulsar timing observations using the Karoo Array Telescope (MeerKAT) of PSR J0514-4002E, an eccentric binary millisecond pulsar in the globular cluster NGC 1851 with a total binary mass of $3.887 \pm 0.004$ solar masses. The companion to the pulsar is a compact object and its mass (between $2.09$ and $2.71$ solar masses, 95% confidence interval) is in the mass gap, so it either is a very massive NS or a low-mass BH. We propose the companion was formed by a merger between two earlier NSs.

The notion of liquid water beneath the ice layer at the south polar layered deposits of Mars is an interesting possibility given the implications for astrobiology, and possible human habitation. A body of liquid water located at a depth of 1.5 km has been inferred from radar data in the South Polar Cap. However, the high temperatures that would facilitate the existence of liquid water or brine at that depth are not consistent with estimations of heat flow that are based on the lithosphere's flexure. Attempts to reconcile both issues have been inconclusive or otherwise unsuccessful. Here, we analyse the possible role of subsurface ammonia and methanol in maintaining water in a liquid state at subsurface temperatures that are compatible with the lithosphere strength. Our results indicate that the presence of these compounds at the base of the south polar layered deposits can reconcile the existence of liquid water with previous estimations of surface heat flow.

Multiband photometric investigations for eight binary systems of the W Ursae Majoris (W UMa)-type are presented. Six systems are presented for the first time to analyze their light curves. All the analyzed systems have a temperature below 5000 K and an orbital period of less than 0.28 days. We extracted primary and secondary minima from the ground-based observations of these systems. According to a few observations reported in the literature, linear fits were considered in the O-C diagrams, and new ephemerides were presented. Light curve solutions were performed using the PHysics Of Eclipsing BinariEs (PHOEBE) code. The results of the mass ratio and fillout factor indicate that the systems are contact binary stars. Six of them showed the O'Connell effect, and a cold starspot on each companion was required for light curve solutions. Their absolute parameters were estimated and evaluated by two other methods. In this study, the empirical relationship between the orbital period and semi-major axis was updated using a sample consisting of 414 contact binary systems and the Monte Carlo Markov Chain (MCMC) approach. Also, using Machine Learning (ML) and the Artificial Neural Network (ANN) model, the relationship between $P-T_1-M_1$ was updated for a better estimation of the mass of the primary star.

Spider pulsars continue to provide promising candidates for neutron star mass measurements. Here we present the discovery of PSR~J1910$-$5320, a new millisecond pulsar discovered in a MeerKAT observation of an unidentified \textit{Fermi}-LAT gamma-ray source. This pulsar is coincident with a recently identified candidate redback binary, independently discovered through its periodic optical flux and radial velocity. New multi-color optical light curves obtained with ULTRACAM/NTT in combination with MeerKAT timing and updated SOAR/Goodman spectroscopic radial velocity measurements allow a mass constraint for PSR~J1910$-$5320. \texttt{Icarus} optical light curve modelling, with streamlined radial velocity fitting, constrains the orbital inclination and companion velocity, unlocking the binary mass function given the precise radio ephemeris. Our modelling aims to unite the photometric and spectroscopic measurements available by fitting each simultaneously to the same underlying physical model, ensuring self-consistency. This targets centre-of-light radial velocity corrections necessitated by the irradiation endemic to spider systems. Depending on the gravity darkening prescription used, we find a moderate neutron star mass of either $1.6\pm0.2$ or $1.4\pm0.2$ $M_\odot$. The companion mass of either $0.45\pm0.04$ or $0.43^{+0.04}_{-0.03}$ $M_\odot$ also further confirms PSR~J1910$-$5320 as an irradiated redback spider pulsar.radiated redback spider pulsar.

We present forecasts for constraints on the Hu \& Sawicki $f(R)$ modified gravity model using realistic mock data representative of future cluster and weak lensing surveys. We create mock thermal Sunyaev-Zel'dovich effect selected cluster samples for SPT-3G and CMB-S4 and the corresponding weak gravitational lensing data from next-generation weak-lensing (ngWL) surveys like Euclid and Rubin. We employ a state-of-the-art Bayesian likelihood approach that includes all observational effects and systematic uncertainties to obtain constraints on the $f(R)$ gravity parameter $\log_{10}|f_{R0}|$. In this analysis we vary the cosmological parameters $[\Omega_{\rm m}, \Omega_\nu h^2, h^2, A_s, n_s, \log_{10}|f_{R0}|]$, which allows us to account for possible degeneracies between cosmological parameters and $f(R)$ modified gravity. The analysis accounts for $f(R)$ gravity via its effect on the halo mass function which is enhanced on cluster mass scales compared to the expectations within general relativity (GR). Assuming a fiducial GR model, the upcoming cluster dataset SPT-3G$\times$ngWL is expected to obtain an upper limit of $\log_{10}|f_{R0}| < -5.95$ at $95\,\%$ credibility, which significantly improves upon the current best bounds. The CMB-S4$\times$ngWL dataset is expected to improve this even further to $\log_{10}|f_{R0}| < -6.23$. Furthermore, $f(R)$ gravity models with $\log_{10}|f_{R0}| \geq -6$, which have larger numbers of clusters, would be distinguishable from GR with both datasets. We also report degeneracies between $\log_{10}|f_{R0}|$ and $\Omega_{\mathrm{m}}$ as well as $\sigma_8$ for $\log_{10}|f_{R0}| > -6$ and $\log_{10}|f_{R0}| > -5$ respectively. Our forecasts indicate that future cluster abundance studies of $f(R)$ gravity will enable substantially improved constraints that are competitive with other cosmological probes.

We report the first interstellar identification of protonated acetylene, C2H3+, a fundamental hydrocarbon, in the z=0.89 molecular absorber toward the gravitationally lensed quasar PKS1830-211. The molecular species is identified from clear absorption features corresponding to the 2_12-1_01 (rest frequency 494.034 GHz) and 1_11-0_00 (431.316 GHz) ground-state transitions of ortho and para forms of C2H3+, respectively, in ALMA spectra toward the southwestern image of PKS1830-211, where numerous molecules, including other hydrocarbons, have already been detected. From the simple assumption of local thermodynamic equilibrium (LTE) with cosmic microwave background photons and an ortho-to-para ratio of three, we estimate a total C2H3+ column density of 2 x 10^12 cm^-2 and an abundance of 10^-10 compared to H_2. However, formation pumping could affect the population of metastable states, yielding a C2H3+ column density higher than the LTE value by a factor of a few. We explore possible routes to the formation of C2H3+, mainly connected to acetylene and methane, and find that the methane route is more likely in PDR environment. As one of the initial hydrocarbon building blocks, C2H3+ is thought to play an important role in astrochemistry, in particular in the formation of more complex organic molecules.

The number of satellites on low orbit has dramatically increased over the past years, raising concerns among the astronomical community about their impact on observations. Spectroscopic observations represent a large fraction of professional data, and spectrographs lack spatial information that can reveal the presence of a satellite. We simulated how often satellites contaminate spectrograph observations by using realistic constellations with over 400,000 objects. We also measured how a spectrum is affected by using real data from different scientific targets and a scaled solar analogue as the satellite, and using standard tools to measure astrophysical parameters and compare them with the clean spectrum. The fraction of affected spectra varies dramatically with the elevation of the sun, with a maximum of 10% at twilight and a nightly average of about 2%. Because of the fast motion of the satellites and the limiting magnitude of the spectrographs, high-resolution instruments are essentially blind to most satellites. For lower resolution spectrographs, the effect on the measured astrophysical parameters depends strongly on the signal-to-noise of the exposure, longer exposures on brighter targets being the least affected at <=1%. Satellites that are brighter and/or higher than the constellation satellites, while less numerous, can also contaminate spectra. While the fraction of affected spectra is likely to remain low, some of these contaminated spectra will be difficult to identify, as it is already the case with existing satellites and asteroids. The best mitigation is to ensure that their brightness is fainter than V=7, that their absolute magnitude V1000km is also fainter than 7, and, whenever possible, to shoot multiple exposures.

Broad sets of spectroscopic observations comprising multiple lines represent an excellent opportunity for diagnostics of the properties of the prominence plasma and the dynamics of their fine structures. However, they also bring significant challenges when they are compared with synthetic spectra provided by radiative transfer modeling. In this work, we provide a statistical spectroscopic analysis of a unique dataset of coordinated prominence observations in the Lyman lines (Ly_alpha to Ly_delta) and the Mg II k and h lines. The observed data were obtained by the Solar Ultraviolet Measurements of Emitted Radiation (SUMER) spectrograph on board of the Solar and Heliospheric Observatory (SoHO) satellite and the Interface Region Imaging Spectrograph (IRIS) on 22 October 2013. We focus on the following profile characteristics: the shape of the observed line profiles based on the number of distinct peaks, the integrated line intensity, the center-to-peak ratio describing the depth of the reversal of two-peaked profiles, and the asymmetry of these peaks. We show that the presence of noise has a negligible effect on the integrated intensity of all observed lines, but it significantly affects the classification of spectral profiles using the number of distinct peaks, the reversal depth, and also the peak asymmetry. We also demonstrate that by taking the influence of noise into account, we can assess which profile characteristics in which spectral lines are suitable for diagnostics of different properties of the observed prominence.

The measure of the abundance of galaxy clusters in the Universe is a sensitive probe of cosmology, sensitive to both the expansion history and the growth of structure. Density fluctuations across the finite survey volume induce noise to this measure, often referred to as Super-Sample Covariance (SSC). In the past for unbinned cluster analysis such noise has not been included in the cluster likelihood. In this paper, we present a derivation of the unbinned likelihood accounting for Super-Sample Covariance by using a Gauss-Poisson Compound (GPC) likelihood. We show that deriving the unbinned likelihood with SSC from the expansion of the GPC formalism to the second order in density perturbation is not sufficient, preventing us from using analytical methods already explored in the literature. In order to solve this issue, we still used the GPC model to derive an alternative "hybrid" likelihood, by using standard redshift bins. Using simulated dark matter halo catalog obtained by the PINOCCHIO algorithm, we found that the hybrid likelihood, accounting for both Poisson noise and SSC, increases the dispersion of the parameter posteriors by 25 per cent using 2,500 clusters, compared to the standard Poisson likelihood.

Ejection of large boulder-like debris is a vigorous form of cometary activity that is unlikely induced by water ice out-gassing alone but rather associated with the sublimation of super-volatile ices. Though perceived on several comets, actual pattern and mechanism of such activity are still unclear. Here we report on a specialized observation of ejections of decimeter- to meter-sized boulders on comet 67P/Churyumov-Gerasimenko outbound between 2.5 and 3.3 AU from the Sun. With a common source region, these events recurred in local morning. The boulders of elongated shapes were ejected in clusters at low inclinations comparable to the solar elevation below 40 degrees at the time. We show that these chunks could be propelled by the surrounding, asymmetric gas field that produced a distinct lateral acceleration. Possibly both water and carbon dioxide have contributed to their mobilization, while the season and local topography are among deciding factors. The mechanisms for sustaining regular activity of comets at large heliocentric distances are likely more diverse and intricate than previously thought.

Aims. Measuring dynamical masses of substellar companions is a powerful tool to test models of mass-luminosity-age relations, as well as determining observational features that constrain the boundary between stellar and substellar companions. In order to dynamically constrain the mass of such companions, we use multiple exoplanet measurement techniques to remove degeneracies in the orbital fits of these objects and place tight constraints on their model-independent masses. Methods. We combine long-period radial-velocity data from the CORALIE survey with relative astrometry from direct imaging with VLT/SPHERE, along with astrometric accelerations from Hipparcos-Gaia eDR3 to perform a combined orbital fit and measure precise dynamical masses of two newly discovered benchmark brown dwarfs. Results. We report the discovery of HD112863B and HD206505B, which are two new benchmark likely brown dwarfs that sit at the substellar-stellar boundary, with precise dynamical masses. We perform an orbital fit which yields dynamical masses for HD112863B and HD206505B to be $77.1^{+2.9}_{-2.8}~M_{\rm{Jup}}$ and $79.8\pm1.8~M_{\rm{Jup}}$ respectively. The orbital period for HD112863B is determined to be $21.59\pm0.05$ years and the orbital period of HD206505B is determined to be ${50.9}_{-1.5}^{+1.7}$ years. From the $H$ and $K$ band photometry from IRDIS data taken with VLT/SPHERE, we estimate the spectral types of both HD112863B and HD206505B to be early-mid L-types.

We show how to significantly improve difference image analysis (DIA) of gravitationally lensed quasars over long periods of time using Gaia proper motions. DIA requires the subtraction of a reference image from the individual images of a monitoring campaign, using stars in the field to align the images. Since the proper motion of the stars can be of the same order as the pixel size during a several-year campaign, we use Gaia DR3 proper motions to enable a correct image alignment. The proper motion corrected star positions can be aligned by the ISIS package. DIA is carried out using the HOTPAnTS package. We apply point spread function (PSF) photometry to obtain light curves and add a proper motion correction of the PSF star to GALFIT. We apply our method to the light curves of the three gravitationally lensed quasars HE1104-1805, HE2149-2745 and Q2237+0305 in the R and V band, respectively, obtained using 1 m telescopes of the Las Cumbres Observatory from 2014 to 2022. We show that the image alignment and the determination of the lensed quasar positions is significantly improved by this method. The light curves of individual quasar images display intrinsic quasar variations and are affected by chromatic microlensing.

The Transiting Exoplanet Survey Satellite (TESS) and the upcoming PLATO mission (PLAnetary Transits and Oscillations of stars) represent two space-based missions with complementary objectives in the field of exoplanet science. While TESS aims at detecting and characterizing exoplanets around bright and nearby stars on a relative short-period orbit, PLATO will discover a wide range of exoplanets including rocky planets within the habitable zones of their stars. We analyze mono-transit events in TESS data around stars that will or could be monitored by the PLATO mission, offering a unique opportunity to bridge the knowledge gap between the two missions and gain deeper insights into exoplanet demographics and system architectures. We found $48$ TESS mono-transit events around stars contained in the all-sky PLATO Input Catalog; of these, at least four will be imaged on the first long-pointing PLATO field, LOPS2. We uniformly vetted this sample to rule out possible false positive detections thus removing $10$ signals from the original sample. We developed an analytic method which allows us to estimate both the orbital period and inclination of a mono-transit planet candidate using only the shape of the transit. We derived the orbital period and inclination estimates for $30$ TESS mono-transit planet candidates. Finally, we investigated whether these candidates are amenable targets for a CHEOPS observing campaign.

Short Gamma-Ray Bursts (GRBs) are often associated with NSNS or BHNS mergers. The discovery of GW/GRB 170817A has enhanced our understanding, revealing that the interaction between relativistic jets and post-merger outflows influences the observed radiation. However, the nature of compact binary merger event suggests that the system can be more complex than the uniform jet interacting with a homologously expanding wind. We consider here an improved scenario by performing a set of two-dimensional, large scale numerical simulations, and we investigate the interaction between short GRB jets and post-merger disk wind outflows. We focus on two types of configurations, arising from NSNS and BHNS mergers. The simulations consider the effects of the r-process nucleosynthesis in the accretion disk wind on its pressure profile. The main properties of the jet, such as its energy distribution and collimation degree, are estimated from our simulations. We found that a) the impact of the r-process on initial wind pressure leads to significant changes in the jet collimation and cocoon expansion; b) the angular structure of thermal and kinetic energy components in the jets, cocoons, and winds differ with respect to simple homologous models, hence it would affect the predictions of GRB afterglow emission; c) the temporal evolution of the structure reveals conversion of thermal to kinetic energy being different for each component in the system (jet, cocoon, and wind); d) post-merger environments influence energy structure and material dispersion, altering the interaction between jets and disk winds. %Our study underscores the importance of post-merger disk wind in the jet propagation, emphasizing the need for careful parameter selection to avoid interpretation degeneracy in the electromagnetic counterparts.

The rotational energy of a fluid parcel changes during isotropic expansion or compression. In solar convection, rotation absorbs energy from convection and inhibits it, causing the motion of fluid parcels larger than a critical size to become vibration. Turbulence and inertial oscillations can cause the deformation of fluid parcels to deviate from isotropic, altering the equilibrium position of the vibration and forming motion larger than the critical size, respectively, the large granules within the granules and probably the mesogranulation. The change in rotational energy of granules during convection causes their rotation speed to differ from the local speed, forming a statistically significant solar radial differential rotation. The meridional circulation driven by radial differential rotation transports angular momentum towards the equator, forming the latitudinal differential rotation. A model constructed by combining mixing length theory explains why granule size and temperature distribution are independent of latitude, and the structure produced by this mechanism is similar to the characteristics of supergranules.

Short GRBs are produced by relativistic jets arising from binary NS-NS or NS-BH mergers. Since the detection of the first unambiguous off-axis GRB 170817A, we learned that energy distribution in the jet plays an important role in explaining the GRB emission. The structure and dynamics are modified during the first seconds of the jet interaction with a post-merger environment. Conventional studies often assume this environment as a simple homologous and symmetrically expanding wind. However, post-merger outflows exhibit complex dynamics influenced by the accretion disc evolution. Moreover, the r-process nucleosynthesis influences the thermodynamics and properties of the post-merger neutron-rich environment. In this work, we study the impact of realistic post-merger disc outflow over the jet dynamics at large scales. We find the results are substantially different from the typical model with symmetric homologous wind.

The in-situ characterization of moon-magnetosphere interactions at Jupiter and the mapping of moon auroral footpaths require accurate global models of the magnetospheric magnetic field. In this study, we compare the ability of two widely-used current sheet models, Khurana-2005 (KK2005) and Connerney-2020 (CON2020) combined with the most recent measurements acquired at low, medium, and high latitudes. With the adjustments of the KK2005 model to JRM33, we show that in the outer and middle magnetosphere (R>15RJ), JRM33+KK2005 is found to be the best model to reproduce the magnetic field observations of Galileo and Juno as it accounts for local time effects. JRM33+CON2020 gives the most accurate representation of the inner magnetosphere. This finding is drawn from comparisons with Juno in-situ magnetic field measurements and confirmed by contrasting the timing of the crossings of the Io, Europa, and Ganymede flux tubes identified in the Juno particles data with the two model estimates. JRM33+CON2020 also maps more accurately the UV auroral footpath of Io, Europa, and Ganymede observed by Juno than JRM33+KK2005. The JRM33+KK2005 model predicts a local time asymmetry in position of the moons' footprints, which is however not detected in Juno's UV measurements.This could indicate that local time effects on the magnetic field are marginal at the orbital locations of Io, Europa, and Ganymede. Finally, the accuracy of the models and their predictions as a function of hemisphere, local time, and longitude is explored.

The Cosmic microwave background (CMB) anisotropies predicted by two cosmological models are compared, one of them is the standard model of general relativity with cold dark matter and cosmological constant, whereas the second model is based on a consistent vector-tensor theory of gravitation explaining solar system and cosmological observations. It is proved that the resulting differences -- between the anisotropies of both models -- are due to the so-called late integrated Sachs Wolfe effect and, consequently, cross correlations between maps of CMB temperatures and tracers of the dark matter distribution could be used in future to select one of the above models. The role of reionization is analysed in detail.

Young protostars that undergo episodic accretion can provide insight into the impact on their circumstellar environments while matter is accreted from the disk onto the protostar. IRAS 22343+7501 is a four component protostar system with one of those being a fading outbursting protostar referred to as L1251 VLA 6. Given the rarity of YSOs undergoing this type of accretion, L1251 VLA 6 can elucidate the fading phase of the post-outburst process. Here we examine structure in the disk around L1251 VLA 6 at frequencies of 10 GHz and 33 GHz with the Karl G. Jansky Very Large Array (VLA). We model the disk structure using Markov chain Monte Carlo (MCMC). This method is then combined with a parametric ray-tracing code to generate synthetic model images of an axisymmetric disk, allowing us to characterize the radial distribution of dust in the system. The results of our MCMC fit show that the most probable values for the mass and radius are consistent with values typical of Class I objects. We find that the total mass of the disk is $0.070^{+0.031}_{-0.2} \rm ~ M_{\sun}$ and investigate the conditions that could cause the accretion outburst. We conclude that the eruption is not caused by gravitational instability and consider alternative explanations and trigger mechanisms.

Fully convective M dwarfs typically remain rapidly rotating and magnetically active for billions of years, followed by an abrupt and mass-dependent transition to slow rotation and quiescence. A robust understanding of this process is complicated by difficulties in estimating M-dwarf ages and potential dependencies on other variables such as birth environment or metallicity. To isolate the effect of mass, we consider M dwarfs in wide binaries. We identify 67 widely separated, fully convective (0.08-0.35M$_\odot$) M-dwarf binary systems using Gaia and measure the H$\alpha$ feature for each component. We classify the pairs into three categories: systems where both components are active, systems where both are inactive, and candidate transition systems, where one component is active and the other inactive. We gather higher-resolution spectra of the candidate transition systems to verify that their behavior does not result from an unresolved third component, yielding one new triple with surprising activity levels. Neglecting this triple, we find 22 active, 36 inactive, and 8 transition pairs. Our results are consistent with the epoch of spindown for these binaries being primarily determined by mass, with mild second-order effects; we place a 1$\sigma$ upper limit of 0.5Gyr or 25% on the dispersion in the mass-dependent spindown relation. Our findings suggest that the large dispersion in spindown epoch previously observed for field stars of a given mass may stem from differences in birth environment, in addition to modest intrinsic stochasticity. We also see evidence that the wide binary population is dispersed over time due to dynamical processing.

We present a strong lensing analysis of COOL J1241+2219, the brightest known gravitationally lensed galaxy at $z \geq 5$, based on new multi-band Hubble Space Telescope (HST) imaging data. The lensed galaxy has a redshift of z=5.043, placing it shortly after the end of the Epoch of Reionization, and an AB magnitude z_AB=20.47 mag (Khullar et al. 2021). As such, it serves as a touchstone for future research of that epoch. The high spatial resolution of HST reveals internal structure in the giant arc, from which we identify 15 constraints and construct a robust lens model. We use the lens model to extract cluster mass and lensing magnification. We find that the mass enclosed within the Einstein radius of the z=1.001 cluster lens is M(<5.77'')=$1.079^{+0.023}_{-0.007}$, significantly lower than other known strong lensing clusters at its redshift. The average magnification of the giant arc is $<\mu_{arc}>=76^{+40}_{-20}$, a factor of $2.4^{+1.4}_{-0.7}$ greater than previously estimated from ground-based data; the flux-weighted average magnification is $<\mu_{arc}>=92^{+37}_{-31}$ We update the current measurements of the stellar mass and star formation rate (SFR) of the source for the revised magnification, $\log(M_\star/M_{\odot})=9.7\pm0.3$ and ${\rm SFR} = 10.3^{+7.0}_{-4.4}$ $ M_{\odot} $yr$^{-1}$. The powerful lensing magnification acting upon COOL J1241+2219 resolves the source and enables future studies of the properties of its star formation on a clump-by-clump basis. The lensing analysis presented here will support upcoming multiwavelength characterization with HST and JWST data of the stellar mass assembly and physical properties of this high-redshift lensed galaxy.

Imaging interferometric data in radio astronomy requires the use of non-linear algorithms that rely on different assumptions on the source structure and may produce non-unique results. This is especially true for Very Long Baseline Interferometry (VLBI) observations, where the sampling of Fourier space is very sparse. A basic tenet in standard VLBI imaging techniques is to assume that the observed source structure does not evolve during the observation. However, the recent VLBI results of the supermassive black hole (SMBH) at our Galactic Center (Sagittarius A$^*$, SgrA*), recently reported by the Event Horizon Telescope Collaboration (EHTC), require the development of dynamic imaging algorithms, since it exhibits variability at minute timescales. In this paper, we introduce a new non-convex optimization problem that extends the standard Maximum Entropy Method (MEM), for reconstructing intra-observation dynamical images from interferometric data that evolves in every integration time. We present a rigorous mathematical formalism to solve the problem via the primal-dual approach. We build a Newton strategy and we give its numerical complexity. We also give a strategy to iteratively improve the obtained solution and finally, we define a novel figure of merit to evaluate the quality of the recovered solution. Then, we test the algorithm, called ngMEM, in different synthetic datasets, with increasing difficulty. Finally, we compare it with another well-established dynamical imaging method. Within this comparison we identified a significant improvement of the ngMEM reconstructions. Moreover, the evaluation of the integration time evolution scheme and the time contribution showed to play a crucial role for obtaining good dynamic reconstructions.

While the standard X-ray variability of black hole X-ray binaries (BHXBs) is stochastic and noisy, there are two known BHXBs that exhibit exotic `heartbeat'-like variability in their light curves: GRS 1915+105 and IGR J17091-3624. In 2022, IGR J17091-3624 went into outburst for the first time in the NICER/NuSTAR era. These exquisite data allow us to simultaneously track the exotic variability and the corresponding spectral features with unprecedented detail. We find that as in typical BHXBs, the outburst began in the hard state, then the intermediate state, but then transitioned to an exotic soft state where we identify two types of heartbeat-like variability (Class V and a new Class X). The flux-energy spectra show a broad iron emission line due to relativistic reflection when there is no exotic variability, and absorption features from highly ionized iron when the source exhibits exotic variability. Whether absorption lines from highly ionized iron are detected in IGR J17091-3624 is not determined by the spectral state alone, but rather is determined by the presence of exotic variability; in a soft spectral state, absorption lines are only detected along with exotic variability. Our finding indicates that IGR J17091-3624 can be seen as a bridge between the most peculiar BHXB GRS 1915+105 and `normal' BHXBs because it alternates between the conventional and exotic behavior of BHXBs. We discuss the physical nature of the absorbing material and exotic variability in light of this new legacy dataset.

IGR J17091-3624 is a black hole X-ray binary (BHXB), often referred to as the 'twin' of GRS 1915+105 because it is the only other known BHXB that can show exotic 'heartbeat'-like variability that is highly structured and repeated. Here we report on observations of IGR J17091-3624 from its 2022 outburst, where we detect an unusually coherent quasi-periodic oscillation (QPO) when the broadband variability is low (total fractional rms $\lesssim$ 6%) and the spectrum is dominated by the accretion disk. Such spectral and variability behavior is characteristic of the soft state of typical BHXBs (i.e., those that do not show heartbeats), but we also find that this QPO is strongest when there is some exotic heartbeat-like variability (so-called Class V variability). This QPO is detected at frequencies between 5 and 8 Hz and has Q-factors (defined as the QPO frequency divided by the width) $\gtrsim$ 50, making it one of the most highly coherent low-frequency QPO ever seen in a BHXB. The extremely high Q factor makes this QPO distinct from typical low-frequency QPOs that are conventionally classified into Type-A/B/C QPOs. Instead, we find evidence that archival observations of GRS 1915+105 also showed a similarly high-coherence QPO in the same frequency range, suggesting that this unusually coherent and strong QPO may be unique to BHXBs that can exhibit 'heartbeat'-like variability.

We currently have two different hypotheses to solve the missing mass problem: dark matter (DM) and modified Newtonian dynamics (MOND). In this work, we use Bayesian inference applied to the Spitzer Photometry and Accurate Rotation Curves (SPARC) galaxies' rotation curves to see which hypothesis fares better. For this, we represent DM by two widely used cusped and cored profiles, Navarro-Frenk-White (NFW) and Burkert. We parameterize MOND by a widely used radial-acceleration relation (RAR). Our results show a preference for the cored DM profile with high Bayes factors in a substantial fraction of galaxies. Interestingly enough, MOND is typically preferred by those galaxies which lack precise rotation curve data. Our study also confirms that the choice of prior has a significant impact on the credible interval of the characteristic MOND acceleration. Overall, our analysis comes out in favor of dark matter.

By measuring photoelectron tracks, the gas pixel detectors of the Imaging X-ray Polarimetry Explorer satellite provide estimates of the photon detection location and its electric vector position angle (EVPA). However, imperfections in reconstructing event positions blur the image and EVPA-position correlations result in artificial polarized halos around bright sources. We introduce a new model describing this "polarization leakage" and use it to recover the on-orbit telescope point-spread functions, useful for faint source detection and image reconstruction. These point spread functions are more accurate than previous approximations or ground-calibrated products ($\Delta \chi^2\approx 3\times 10^{4}$ and $4 \times 10^4$ respectively for a bright $10^6$-count source). We also define an algorithm for polarization leakage correction substantially more accurate than existing prescriptions ($\Delta \chi^2\approx 1\times 10^{3}$). These corrections depend on the reconstruction method, and we supply prescriptions for the mission-standard "Moments" methods as well as for "Neural Net" event reconstruction. Finally, we present a method to isolate leakage contributions to polarization observations of extended sources and show that an accurate PSF allows the extraction of sub-PSF-scale polarization patterns.

The effect of gravity in Maxwell's equations is often treated as a medium property. The commonly used formulation is based on managing Maxwell's equations in exactly the same form as in Minkowski spacetime and expressing the effect of gravity as a set of constitutive relations. We show that such a set of Maxwell's equations is, in fact, a combination of the electric and magnetic fields defined in two different non-covariant ways, both of which fail to identify the associated observer's four-vectors. The suggested constitutive relations are also ad hoc and unjustified. To an observer with a proper four-vector, the effect of gravity can be arranged as effective polarizations and magnetizations appearing in both the homogeneous and inhomogeneous parts. Modifying the homogeneous part by gravity is inevitable to any observer, and the result cannot be interpreted as the medium property. For optical properties one should directly handle Maxwell's equations in curved spacetime.

We consider a charged BTZ black hole in asymptotically AdS space-time of massive gravity to study the effect of the thermal fluctuations on the black hole thermodynamics. We consider the Einstein-Born-Infeld solution and investigate critical points and stability. We also compare the results with the case of Einstein-Maxwell solutions. Besides, we find that thermal fluctuations, which appear as a logarithmic term in the entropy, affect the stability of the black hole and change the phase transition point. Moreover, we study the geometrical thermodynamics and find that the behaviour of the linear Maxwell solution is the same as the nonlinear one.

The species scale is a field-dependent UV cut-off for any effective field theory weakly coupled to gravity. In this letter we show that, in the context of inflationary cosmology, a detection of primordial gravitational waves will set an upper bound on the decay rate $|\Lambda'_s/\Lambda_s|$ of the species scale. Specifically, we derive this in terms of the tensor-to-scalar ratio $r$ of power spectra of primordial perturbations. Given the targets of current and next generation experiments, we show that any successful detection would signify that this upper limit is of the order of unity, which is consistent with recent discussions in the literature.

Adiabatic binary inspiral in the small mass ratio limit treats the small body as moving along a geodesic of a large Kerr black hole, with the geodesic slowly evolving due to radiative backreaction. Up to initial conditions, geodesics are typically parameterized in two ways: using the integrals of motion energy $E$, axial angular momentum $L_z$, and Carter constant $Q$; or, using orbit geometry parameters semi-latus rectum $p$, eccentricity $e$, and (cosine of ) inclination $x_I \equiv \cos I$. The community has long known how to compute orbit integrals as functions of the orbit geometry parameters, i.e., as functions expressing $E(p, e, x_I)$, and likewise for $L_z$ and $Q$. Mappings in the other direction -- functions $p(E, L_z, Q)$, and likewise for $e$ and $x_I$ -- have not yet been developed in general. In this note, we develop generic mappings from ($E$, $L_z$, $Q$) to ($p$, $e$, $x_I$). The mappings are particularly simple for equatorial orbits ($Q = 0$, $x_I = \pm1$), and can be evaluated efficiently for generic cases. These results make it possible to more accurately compute adiabatic inspirals by eliminating the need to use a Jacobian which becomes singular as inspiral approaches the last stable orbit.

This dissertation aims to deepen the understanding of the primordial composition of the Universe in the temperature range 300 MeV>T>0.02 MeV. I exploit known properties of elementary particles and apply methods of kinetic theory and statistical physics to advance the understanding of the cosmic plasma. Within the Big Bang model, we begin by considering the Universe being a highly energetic fireball, an ultra-relativistic plasma exhibiting distinct properties. Fundamental particles such as quarks, leptons, and even heavier gauge bosons play a crucial role in the understanding of the early Universe. Our research focuses on the investigation of these fundamental particles as constituents of the dense Universe plasma during the epoch which transits from primordial quark-gluon plasma to the era of normal hadron matter, passing through the decoupling of neutrinos and addressing in detail the electron-positron antimatter plasma.

The effective action in renormalizable quantum theory of gravity provides entropy because the total Hamiltonian vanishes. Since it is a renormalization group invariant that is constant in the process of cosmic evolution, we can show conservation of entropy, that is an ansatz in the standard cosmology. Here we study renormalizable quantum gravity that exhibits conformal dominance at high energy beyond the Planck scale. The current entropy of the universe is derived by calculating the effective action under the scenario of quantum gravity inflation caused by its dynamics. We then argue that ghost modes must be unphysical, but necessary for the Hamiltonian to vanish and for entropy to exist in gravitational systems.

In this lecture we discuss the properties of dense hadronic matter inside neutron stars. In particular, we pay attention to the role of strangeness in the core of neutron stars, by analysing the presence of baryons and mesons with strangeness. We consider two interesting possible scenarios in their interior, that is, the existence of hyperons leading to the so-called hyperon puzzle and the presence of a kaon condensed phase inside neutron stars.

We have studied the stability of C$_{59}$ anions as a function of time, from their formation on femtosecond timescales to their stabilization on second timescales and beyond, using a combination of theory and experiments. The C$_{59}^-$ fragments were produced in collisions between C$_{60}$ fullerene anions and neutral helium gas at a velocity of 90 km/s (corresponding to a collision energy of 166\,eV in the center-of-mass frame). The fragments were then stored in a cryogenic ion-beam storage ring at the DESIREE facility where they were followed for up to one minute. Classical molecular dynamics simulations were used to determine the reaction cross section and the excitation energy distributions of the products formed in these collisions. We found that about 15 percent of the C$_{59}^-$ ions initially stored in the ring are intact after about 100 ms, and that this population then remains intact indefinitely. This means that C$_{60}$ fullerenes exposed to energetic atoms and ions, such as stellar winds and shock waves, will produce stable, highly reactive products, like C$_{59}$, that are fed into interstellar chemical reaction networks.

We investigate symmetric Metric-Affine Theories of Gravity with a Lagrangian containing all operators of dimension up to four that are relevant to free propagation in flat space. Complementing recent work in the antisymmetric case, we derive the conditions for the existence of a single massive particle with good properties, in addition to the graviton.

We study the statistics of scalar perturbations in models of inflation with small and rapid oscillations in the inflaton potential (resonant non-Gaussianity). We do so by deriving the wavefunction $\Psi[\zeta(\boldsymbol{x})]$ non-perturbatively in $\zeta$, but at first order in the amplitude of the oscillations. The expression of the wavefunction of the universe (WFU) is explicit and does not require solving partial differential equations. One finds qualitative deviations from perturbation theory for $ |\zeta| \gtrsim \alpha^{-2}$, where $\alpha \gg 1$ is the number of oscillations per Hubble time. Notably, the WFU exhibits distinct behaviours for negative and positive values of $\zeta$ (troughs and peaks respectively). While corrections for $\zeta <0$ remain relatively small, of the order of the oscillation amplitude, positive $\zeta$ yields substantial effects, growing exponentially as $e^{\pi\alpha/2}$ in the limit of large $\zeta$. This indicates that even minute oscillations give large effects on the tail of the distribution.