We reconstruct the cosmological background evolution under the scenario of dynamical dark energy through the Gaussian process approach, using the latest Dark Energy Spectroscopic Instrument (DESI) baryon acoustic oscillations (BAO) combined with other observations. Our results reveal that the reconstructed dark-energy equation-of-state (EoS) parameter $w(z)$ exhibits the so-called quintom-B behavior, crossing $-1$ from phantom to quintessence regime as the universe expands. We investigate under what situation this type of evolution could be achieved from the perspectives of field theories and modified gravity. In particular, we reconstruct the corresponding actions for $f(R)$, $f(T)$, and $f(Q)$ gravity, respectively. We explicitly show that, certain modified gravity can exhibit the quintom dynamics and fit the recent DESI data efficiently, and for all cases the quadratic deviation from the $\Lambda$CDM scenario is mildly favored.

In this article we use the latest cosmological observations, including SNe, BAO, CC and RSD, to reconstruct the cosmological evolution via the Gaussian process. At the background level, we find consistency with the quintom dynamics for different data combinations and divide the characteristics of dark energy into three different categories, which are negative-energy dark energy, late-dominated dark energy and oscillating dark energy, respectively. Considering the effect of modified gravity on the growth of matter perturbations, the reconstruction results at the perturbative level show that we only need minor corrections to general relativity. Furthermore, we provide theoretical interpretation for the three different types of dynamical dark-energy behavior, in the framework of modified gravity, scalar fields, and dark-energy equation-of-state parametrizations. Finally, we show that all of these models can be unified in the framework of effective field theory.

We apply a novel statistical analysis to measurements of 16 elemental abundances in 34,410 Milky Way disk stars from the final data release (DR17) of APOGEE-2. Building on recent work, we fit median abundance ratio trends [X/Mg] vs. [Mg/H] with a 2-process model, which decomposes abundance patterns into a "prompt" component tracing core collapse supernovae and a "delayed" component tracing Type Ia supernovae. For each sample star, we fit the amplitudes of these two components, then compute the residuals \Delta[X/H] from this two-parameter fit. The rms residuals range from ~0.01-0.03 dex for the most precisely measured APOGEE abundances to ~0.1 dex for Na, V, and Ce. The correlations of residuals reveal a complex underlying structure, including a correlated element group comprised of Ca, Na, Al, K, Cr, and Ce and a separate group comprised of Ni, V, Mn, and Co. Selecting stars poorly fit by the 2-process model reveals a rich variety of physical outliers and sometimes subtle measurement errors. Residual abundances allow comparison of populations controlled for differences in metallicity and [\alpha/Fe]. Relative to the main disk (R=3-13 kpc, |Z|<2 kpc), we find nearly identical abundance patterns in the outer disk (R=15-17 kpc), 0.05-0.2 dex depressions of multiple elements in LMC and Gaia Sausage/Enceladus stars, and wild deviations (0.4-1 dex) of multiple elements in \omega Cen. Residual abundance analysis opens new opportunities for discovering chemically distinctive stars and stellar populations, for empirically constraining nucleosynthetic yields, and for testing chemical evolution models that include stochasticity in the production and redistribution of elements.

We investigate the realization of quintom scenario for dynamical dark energy within modified gravity theories that can efficiently fit the recent observational datasets. Starting from a general effective field theory formulation of dark energy in metric-affine geometry, we derive the background action in unitary gauge and we demonstrate how both $f(T)$ and $f(Q)$ gravity can naturally realize quintom behavior through appropriate forms and parameter choices. Additionally, using the Gaussian process reconstruction of the latest DESI DR2 BAO data combined with SNe and CMB observations, we extract the reconstructed dark-energy equation-of-state parameter, showing that it exhibits quintom-type evolution, crossing the phantom divide from below. Moreover, through detailed parameter estimations and application of information criteria, we compare the model with the quadratic one and the $\Lambda$CDM model. Our results show that, due to its rich structure, modified gravity stands as one of the main candidates for the realization of the data-favoured dynamical dark energy.

We investigate the observational implications of Wald--Gauss--Bonnet (WGB) topological dark energy, a modified cosmological framework derived from the gravity-thermodynamics conjecture applied to the Universe's apparent horizon, with the Wald--Gauss--Bonnet entropy replacing the standard Bekenstein--Hawking one. Assuming a topological connection between the apparent horizon and interior black hole (BH) horizons, we derive modified Friedmann equations where the evolution of dark energy depends on BH formation and merger rates, which are approximated by the cosmic star formation rate. These equations introduce an additional, astrophysics--dependent contribution to the cosmological constant. We test two scenarios, one with a vanishing cosmological constant ($\Lambda = 0$) and another with a modified $\Lambda$ against late--Universe data (SNIa, BAO, Cosmic Chronometers) via a Bayesian analysis. Although the WGB framework is consistent with observations, information criteria statistically favor the standard $\Lambda$CDM model. An analysis of linear perturbations shows that the growth of cosmic structures is nearly indistinguishable from that of $\Lambda$CDM, with negligible dark energy clustering and minimal deviation in the effective Newton's constant. The standard thermal history is also preserved. In conclusion, WGB cosmology presents a phenomenologically rich alternative that connects dark energy to black hole astrophysics while remaining compatible with current cosmological data.

We investigate the relationship of kinematically misaligned galaxies with their large-scale environment, in the context of halo assembly bias. According to numerical simulations, halo age at fixed halo mass is intrinsically linked to the large-scale tidal environment created by the cosmic web. We investigate the relationship between distances to various cosmic web features and present-time gas accretion rate. We select a sub-sample of ~900 central galaxies from the MaNGA survey with defined global position angles (PA; angle at which velocity change is greatest) for their stellar and H$\alpha$ gas components up to a minimum of 1.5 effective radii ($R_e$). We split the sample by misalignment between the gas and stars as defined by the difference in their PA. For each central galaxy we find its distance to nodes and filaments within the cosmic web, and estimate the host halo's age using the central stellar mass to total halo mass ratio $M_{*}/M_{h}$. We also construct halo occupation distributions using a background subtraction technique for galaxy groups split using the central galaxy's kinematic misalignment. We find, at fixed halo mass, no statistical difference in these properties between our kinematically aligned and misaligned galaxies. We suggest that the lack of correlation could be indicative of cooling flows from the hot halo playing a far larger role than `cold mode' accretion from the cosmic web or a demonstration that the spatial extent of current large-scale integral field unit (IFU) surveys hold little information about large-scale environment extractable through this method.

We present how Tsallis cosmology can alleviate both $H_0$ and $\sigma_8$ tensions simultaneously. Such a modified cosmological scenario is obtained by the application of the gravity-thermodynamics conjecture, but using the non-additive Tsallis entropy, instead of the standard Bekenstein-Hawking one. Hence, one obtains modified Friedmann equations, with extra terms that depend on the new Tsallis exponent $\delta$ that quantifies the departure from standard entropy. We show that for particular $\delta$ choices we can obtain a phantom effective dark energy, which is known to be one of the sufficient mechanisms that can alleviate $H_0$ tension. Additionally, for the same parameter choice we obtain an increased friction term and an effective Newton's constant smaller than the usual one, and thus the $\sigma_8$ tension is also solved. These features act as a significant advantage of Tsallis modified cosmology.

We present a joint analysis of TTVs and Doppler data for the transiting exoplanet system TOI-4504. TOI-4504 c is a warm Jupiter-mass planet that exhibits the largest known transit timing variations (TTVs), with a peak-to-node amplitude of $\sim$ 2 days, the largest value ever observed, and a super-period of $\sim$ 930 d. TOI-4504 b and c were identified in public TESS data, while the TTVs observed in TOI-4504 c, together with radial velocity (RV) data collected with FEROS, allowed us to uncover a third, non-transiting planet in this system, TOI-4504 d. We were able to detect transits of TOI-4504 b in the TESS data with a period of 2.4261$\pm 0.0001$ days and derive a radius of 2.69$\pm 0.19$ R$_{\oplus}$. The RV scatter of TOI-4504 was too large to constrain the mass of TOI-4504 b, but the RV signals of TOI-4504 c \& d were sufficiently large to measure their masses. The TTV+RV dynamical model we apply confirms TOI-4504 c as a warm Jupiter planet with an osculating period of 82.54$\pm 0.02$ d, mass of 3.77$\pm 0.18$ M$_{\rm J}$ and a radius of 0.99$\pm 0.05$R$_{\rm J}$, while the non-transiting planet TOI-4504 d, has an orbital period of 40.56$\pm 0.04$ days and mass of 1.42$_{-0.06}^{+0.07}$ M$_{\rm J}$. We present the discovery of a system with three exoplanets: a hot sub-Neptune and two warm Jupiter planets. The gas giant pair is stable and likely locked in a first-order 2:1 mean-motion resonance (MMR). The TOI-4504 system is an important addition to MMR pairs, whose increasing occurrence supports a smooth migration into a resonant configuration during the protoplanetary disk phase.

We report the discovery and confirmation of two planets orbiting the metal-poor Sun-like star, HD 35843 (TOI 4189). HD 35843 c is a temperate sub-Neptune transiting planet with an orbital period of 46.96 days that was first identified by Planet Hunters TESS. We combine data from TESS and follow-up observations to rule out false-positive scenarios and validate the planet. We then use ESPRESSO radial velocities to confirm the planetary nature and characterize the planet's mass and orbit. Further analysis of these RVs reveals the presence of an additional planet, HD 35843 b, with a period of 9.90 days and a minimum mass of $5.84\pm0.84$ $M_{\oplus}$. For HD 35843 c, a joint photometric and spectroscopic analysis yields a radius of $2.54 \pm 0.08 R_{\oplus}$, a mass of$11.32 \pm 1.60 M_{\oplus}$, and an orbital eccentricity of $e = 0.15\pm0.07$. With a bulk density of $3.80 \pm 0.70$ g/cm$^3$, the planet might be rocky with a substantial H$_2$ atmosphere or it might be a ``water world". With an equilibrium temperature of $\sim$480 K, HD 35843 c is among the coolest $\sim 5\%$ of planets discovered by TESS. Combined with the host star's relative brightness (V= 9.4), HD 35843 c is a promising target for atmospheric characterization that will probe this sparse population of temperate sub-Neptunes.

We explore the origins of the young B-type stars found by Casetti-Dinescu et al.(2014) at the outskirts of the Milky-Way disk in the sky region of Leading Arm of the Magellanic Stream. High-resolution spectroscopic observations made with the MIKE instrument on the Magellan Clay 6.5m telescope for nine stars are added to the previous sample analyzed by Zhang et al. (2017). We compile a sample of fifteen young stars with well-determined stellar types, ages, abundances and kinematics. With proper motions from Gaia DR2 we also derive orbits in a realistic Milky-Way potential. We find that our previous radial-velocity selected LA candidates have substantial orbital angular momentum. The substantial amount of rotational component for these stars is in contrast with the near-polar Magellanic orbit, thus rendering these stars unlikely members of the LA. There are four large orbital-energy stars in our sample. The highest orbital-energy one has an age shorter than the time to disk crossing, with a birthplace $z=2.5$~kpc and $R_{\rm GC}\sim 28$~kpc. Therefore, the origin of this star is uncertain. The remaining three stars have disk runaway origin with birthplaces between 12 and 25 kpc from the Galactic center. Also, the most energetic stars are more metal poor ([Mg/H] =$-0.50\pm0.07$) and with larger He scatter ($\sigma_{\rm [He/H]} = 0.72$) than the inner disk ones ([Mg/H] $=0.12\pm0.36$, $\sigma_{\rm [He/H]} = 0.15$). While the former group's abundance is compatible with that of the Large Magellanic Cloud, it could also reflect the metallicity gradient of the MW disk and their runaway status via different runaway mechanisms.

Ultra-light primordial black holes (PBHs) with masses M_{\rm PBH}&lt;5\times 10^8\mathrm{g} can dominate transiently the energy budget of the Universe and reheat the Universe through their evaporation taking place before Big Bang Nucleosynthesis. The isocurvature energy density fluctuations associated to the inhomogeneous distribution of a population of such PBHs can induce an abundant production of GWs due to second-order gravitational effects. In this work, we discuss the effect of primordial non-Gaussianity on the clustering properties of PBHs and study the effect of a clustered PBH population on the spectral shape of the aforementioned induced GW signal. In particular, focusing on local-type non-Gaussianity we find a double-peaked GW signal with the amplitude of the low-frequency peak being proportional to the square of the non-Gaussian parameter $\tau_\mathrm{NL}$. Remarkably, depending on the PBH mass $M_{\rm PBH}$ and the initial abundance of PBHs at formation time, i.e. $\Omega_\mathrm{PBH,f}$, this double-peaked GW signal can lie well within the frequency bands of forthcoming GW detectors, namely LISA, ET, SKA and BBO, hence rendering this signal falsifiable by GW experiments and promoting it as a novel portal probing the primordial non-Gaussianity.

The APOGEE Open Cluster Chemical Abundances and Mapping (OCCAM) survey is used to probe the chemical evolution of the s-process element cerium in the Galactic disk. Cerium abundances were derived from measurements of Ce II lines in the APOGEE spectra using the Brussels Automatic Code for Characterizing High Accuracy Spectra (BACCHUS) in 218 stars belonging to 42 open clusters. Our results indicate that, in general, for Ages &lt; 4 Gyr, younger open clusters have higher [Ce/Fe] and [Ce/$\alpha$-element] ratios than older clusters. In addition, metallicity segregates open clusters in the [Ce/X]-Age plane (where X can be H, Fe, and the $\alpha$-elements O, Mg, Si, or Ca). These metallicity-dependant relations result in [Ce/Fe] and [Ce/$\alpha$] ratios with age that are not universal clocks. Radial gradients of [Ce/H] and [Ce/Fe] ratios in open clusters, binned by age, were derived for the first time, with d[Ce/H]dR$_{GC}$ being negative, while d[Ce/Fe]/dR$_{GC}$ is positive. [Ce/H] and [Ce/Fe] gradients are approximately constant over time, with the [Ce/Fe] gradient becoming slightly steeper, changing by $\sim$+0.009 dex-kpc$^{-1}$-Gyr$^{-1}$. Both the [Ce/H] and [Ce/Fe] gradients are shifted to lower values of [Ce/H] and [Ce/Fe] for older open clusters. The chemical pattern of Ce in open clusters across the Galactic disk is discussed within the context of s-process yields from AGB stars, $\sim$Gyr time delays in Ce enrichment of the interstellar medium, and the strong dependence of Ce nucleosynthesis on the metallicity of its AGB stellar sources.

We study the scenario of Kanadiakis horizon entropy cosmology which arises from the application of the gravity-thermodynamics conjecture using the Kaniadakis modified entropy. The resulting modified Friedmann equations contain extra terms that constitute an effective dark energy sector. We use data from Cosmic chronometers, Supernova Type Ia, HII galaxies, Strong lensing systems, and Baryon acoustic oscillations observations and we apply a Bayesian Markov Chain Monte Carlo analysis to construct the likelihood contours for the model parameters. We find that the Kaniadakis parameter is constrained around 0, namely, around the value where the standard Bekenstein-Hawking is recovered. Concerning the normalized Hubble parameter, we find $h=0.708^{+0.012}_{-0.011}$, a result that is independently verified by applying the $\mathbf{\mathbb{H}}0(z)$ diagnostic and, thus, we conclude that the scenario at hand can alleviate the $H_0$ tension problem. Regarding the transition redshift, the reconstruction of the cosmographic parameters gives $z_T=0.715^{+0.042}_{-0.041}$. Furthermore, we apply the AICc, BIC and DIC information criteria and we find that in most datasets the scenario is statistical equivalent to $\Lambda$CDM one. Moreover, we examine the Big Bang Nucleosynthesis (BBN) and we show that the scenario satisfies the corresponding requirements. Additionally, we perform a phase-space analysis, and we show that the Universe past attractor is the matter-dominated epoch, while at late times the Universe results in the dark-energy-dominated solution. Finally, we show that Kanadiakis horizon entropy cosmology accepts heteroclinic sequences, but it cannot exhibit bounce and turnaround solutions.

We investigate the strong coupling problem in modified teleparallel gravity theories using the effective field theory (EFT) approach, demonstrating that it is possible to shift the emergence of new degrees of freedom (DoFs) to lower orders in perturbation theory. We first focus on the case of $f(T)$ gravity, and we show that in its conformally equivalent form the scalar perturbations are non-dynamical up to the cubic action. We then propose a simple modification of the theory, which lowers the appearance of new DoFs to cubic order, compared to the quartic order in standard $f(T)$ gravity. Our work opens a new avenue to address the issue of strong coupling in modified teleparallel gravity, and suggests a new classification scheme of these theories based on the perturbative order at which new DoFs appear.

As part of our ongoing effort to investigate transit timing variations (TTVs) of known exoplanets, we monitored transits of the four exoplanets HAT-P-18b, HAT-P-19b, HAT-P-27b/WASP-40b and WASP-21b. All of them are suspected to show TTVs due to the known properties of their host systems based on the respective discovery papers. During the past three years 46 transit observations were carried out, mostly using telescopes of the Young Exoplanet Transit Initiative. The analyses are used to refine the systems orbital parameters. In all cases we found no hints for significant TTVs, or changes in the system parameters inclination, fractional stellar radius and planet to star radius ratio. However, comparing our results with those available in the literature shows that we can confirm the already published values.

We argue that the effect of cold dark matter in the cosmological setup can be explained by the coupling between the baryonic matter particles in terms of the long-range force having a graviton mass $m_g$ via the Yukawa gravitational potential. Such a quantum-corrected Yukawa-like gravitational potential is characterized by the coupling parameter $\alpha$, the wavelength parameter $\lambda$, which is related to the graviton mass via $m_g=\hbar/(\lambda c)$, that determines the range of the force and, finally, a Planck length quantity $l_0$ that makes the potential regular at the centre. The corrected Friedmann equations are obtained using Verlinde's entropic force interpretation of gravity based on the holographic scenario and the equipartition law of energy. The parameter $\alpha$ modifies Newton's constant as $G_{\rm eff}\rightarrow G(1+\alpha)$. We argue that dark matter is an apparent effect as no dark matter particle exists in this picture. Furthermore, the dark energy is also related to graviton mass and $\alpha$; in particular, we point out that the cosmological constant can be viewed as a self-interaction effect between gravitons. We further show that there exists a precise correspondence with Verlinde's emergent gravity theory, and due to the long-range force, the theory can be viewed as a non-local gravity theory. To this end, we performed the phase space analyses and estimated $\lambda\simeq 10^3$[Mpc] and $\alpha \in (0.0385,0.0450)$, respectively. Finally, from these values, for the graviton mass, we get $m_g\simeq 10^{-68}$ kg, and cosmological constant $\Lambda \simeq 10^{-52}\rm m^{-2}$. Further, we argue how this theory reproduces the MOND phenomenology on galactic scales via the acceleration of Milgrom $a_0 \simeq 10^{-10}\rm m/s^2$.

We report the discovery of the transiting exoplanet NGTS-12b by the Next Generation Transit Survey (NGTS). The host star, NGTS-12, is a V=12.38 mag star with an effective temperature of T$_{\rm eff}$=$5690\pm130$ K. NGTS-12b orbits with a period of $P=7.53$d, making it the longest period planet discovered to date by the main NGTS survey. We verify the NGTS transit signal with data extracted from the TESS full-frame images, and combining the photometry with radial velocity measurements from HARPS and FEROS we determine NGTS-12b to have a mass of $0.208\pm0.022$ M$_{J}$ and a radius of $1.048\pm0.032$ R$_{J}$. NGTS-12b sits on the edge of the Neptunian desert when we take the stellar properties into account, highlighting the importance of considering both the planet and star when studying the desert. The long period of NGTS-12b combined with its low density of just $0.223\pm0.029$ g cm$^{-3}$ make it an attractive target for atmospheric characterization through transmission spectroscopy with a Transmission Spectroscopy Metric of 89.4.

We study the Einstein-aether theory in Weyl integrable geometry. The scalar field which defines the Weyl affine connection is introduced in the gravitational field equation. We end up with an Einstein-aether scalar field model where the interaction between the scalar field and the aether field has a geometric origin. The scalar field plays a significant role in the evolution of the gravitational field equations. We focus our study on the case of homogeneous and isotropic background spacetimes and study their dynamical evolution for various cosmological models.

The unification of quantum mechanics and general relativity has long been elusive. Only recently have empirical predictions of various possible theories of quantum gravity been put to test, where a clear signal of quantum properties of gravity is still missing. The dawn of multi-messenger high-energy astrophysics has been tremendously beneficial, as it allows us to study particles with much higher energies and travelling much longer distances than possible in terrestrial experiments, but more progress is needed on several fronts. A thorough appraisal of current strategies and experimental frameworks, regarding quantum gravity phenomenology, is provided here. Our aim is twofold: a description of tentative multimessenger explorations, plus a focus on future detection experiments. As the outlook of the network of researchers that formed through the COST Action CA18108 ``Quantum gravity phenomenology in the multi-messenger approach (QG-MM)'', in this work we give an overview of the desiderata that future theoretical frameworks, observational facilities, and data-sharing policies should satisfy in order to advance the cause of quantum gravity phenomenology.