Galactic Systems Research Articles

A Radio-loud Semiregular Variable

As a byproduct of our search for Galactic stellar systems with gamma-ray emission, we have identified an unrelated cool and evolved star (IRC-10412) that attracted our attention due to its strong radio emission level with a spectral index matching, almost perfectly, the canonical +0.6 value expected from an ionized stellar wind. A follow-up observational analysis was undertaken given that these two properties are hard to reconcile as originating in the same stellar object. As a result, IRC-10412 has been classified as a new semiregular variable of SRb type in the asymptotic giant branch, and different but consistent estimates of its mass-loss parameter are reported. We propose that its unusually high radio emission arises from a ∼10−5 M ⊙ yr−1 stellar wind exposed to an external source of ionizing photons, possibly coming from nearby OB associations.

Galactic Stellar Black Hole Binaries: Spin Effects on Jet Emissions of High-Energy Gamma-Rays

In the last few decades, galactic stellar black hole X-ray binary systems (BHXRBs) have aroused intense observational and theoretical research efforts specifically focusing on their multi-messenger emissions (radio waves, X-rays, γ-rays, neutrinos, etc.). In this work, we investigate jet emissions of high-energy neutrinos and gamma-rays created through several hadronic and leptonic processes taking place within the jets. We pay special attention to the effect of the black hole’s spin (Kerr black holes) on the differential fluxes of photons originating from synchrotron emission and inverse Compton scattering and specifically on their absorption due to the accretion disk’s black-body radiation. The black hole’s spin (dimensionless spin parameter a*) enters into the calculations through the radius of the innermost circular orbit around the black hole, the RISCO parameter, assumed to be the inner radius of the accretion disk, which determines its optical depth τdisk. In our results, the differential photon fluxes after the absorption effect are depicted as a function of the photon energy in the range 1GeV ≤E≤103GeV. It is worth noting that when the black holes’ spin (α*) increases, the differential photon flux becomes significantly lower.

Study of the Effectiveness of Parallel Algorithms for Modeling the Dynamics of Collisionless Galactic Systems on GPUs

Study of the Effectiveness of Parallel Algorithms for Modeling the Dynamics of Collisionless Galactic Systems on GPUs

Investigation of the radial profile of galactic magnetic fields using rotation measure of background quasars

Probing magnetic fields in high-redshift galactic systems is crucial to investigate galactic dynamics and evolution. Utilizing the rotation measure of the background quasars, we have developed a radial profile of the magnetic field in a typical high-z galaxy. We have compiled a catalog of 59 confirmed quasar sightlines, having one intervening Mg II absorber in the redshift range 0.372 ≤ z abs ≤ 0.8. The presence of the foreground galaxy is ensured by comparing the photometric and spectroscopic redshifts within 3σ z-photo and visual checks. These quasar line-of-sights (LoS) pass through various impact parameters (D) up to 160 kpc, covering the circumgalactic medium of a typical Milky-Way type galaxy. Utilizing the residual rotation measure (RRM) of these sightlines, we estimated the excess in RRM dispersion, σ ex RRM. We found σ ex RRM decreases with increasing D. We translated σ ex RRM to average LoS magnetic field strength, 〈 B ∥ 〉 by considering a typical electron column density. Consequently, the decreasing trend is sustained in the magnetic field. In particular for sightlines with D ≤ 50 kpc and D> 50 kpc, 〈 B ∥ 〉 is found to be 2.39 ± 0.7 μG and 1.67 ± 0.38 μG, respectively. This suggests a clear indication of varying magnetic field from the disk to the circumgalactic medium. This work provides a methodology that, when applied to ongoing and future radio polarisation surveys such as LOFAR and SKA, promises to significantly enhance our understanding of magnetic field mapping in galactic systems.

Revisiting the Constraint on the Equation of State of Neutron Stars Based on Binary Neutron Star Mergers

The merger of a neutron star (NS)–NS binary can form different productions of compact remnants, among which a supramassive NS (SMNS) can create an internal plateau, and the following steep decay marks the collapse of the SMNS. The proportion of the SMNS and the corresponding collapse time are often used to constrain the NS equation of state (EOS). This paper revisits this topic by considering the effect of an accretion disk on a compact remnant, which is not considered in previous works. Compared with previous works, the collapse-time distribution (peaks ∼100 s) of SMNSs formed from an NS–NS merger is almost unaffected by the initial surface magnetic field (B s,i ) of the NS, but the total energy output of the magnetic dipole radiation from the SMNSs depends on B s,i significantly. Coupling the constraints from the SMNS fraction, we exclude some EOSs and obtain three candidate EOSs, i.e., DD2, ENG, and MPA1. By comparing the distributions of the collapse time and the luminosity of the internal plateau (in the short gamma-ray bursts) for observations obtained based on the three candidate EOSs, it is shown that only the EOS of ENG is favored. Our sample, based on the ENG EOS and a mass distribution motivated by Galactic systems, suggests that approximately 99% of NS–NS mergers collapse to form a black hole (BH) within 107s. This includes scenarios promptly forming a BH (36.5%), an SMNS (60.7%), or a stable NS that transitions into a BH or an SMNS following accretion (2.1%). It also indicates that the remnants for GW170817 and GW190425, and the second object of GW190814, are more likely to be BHs.

Discovery of a dormant 33 solar-mass black hole in pre-release Gaia astrometry

Context. Gravitational waves from black-hole (BH) merging events have revealed a population of extra-galactic BHs residing in short-period binaries with masses that are higher than expected based on most stellar evolution models – and also higher than known stellar-origin black holes in our Galaxy. It has been proposed that those high-mass BHs are the remnants of massive metal-poor stars. Aims. Gaia astrometry is expected to uncover many Galactic wide-binary systems containing dormant BHs, which may not have been detected before. The study of this population will provide new information on the BH-mass distribution in binaries and shed light on their formation mechanisms and progenitors. Methods. As part of the validation efforts in preparation for the fourth Gaia data release (DR4), we analysed the preliminary astrometric binary solutions, obtained by the Gaia Non-Single Star pipeline, to verify their significance and to minimise false-detection rates in high-mass-function orbital solutions. Results. The astrometric binary solution of one source, Gaia BH3, implies the presence of a 32.70 ± 0.82 M⊙ BH in a binary system with a period of 11.6 yr. Gaia radial velocities independently validate the astrometric orbit. Broad-band photometric and spectroscopic data show that the visible component is an old, very metal-poor giant of the Galactic halo, at a distance of 590 pc. Conclusions. The BH in the Gaia BH3 system is more massive than any other Galactic stellar-origin BH known thus far. The low metallicity of the star companion supports the scenario that metal-poor massive stars are progenitors of the high-mass BHs detected by gravitational-wave telescopes. The Galactic orbit of the system and its metallicity indicate that it might belong to the Sequoia halo substructure. Alternatively, and more plausibly, it could belong to the ED-2 stream, which likely originated from a globular cluster that had been disrupted by the Milky Way.

The oldest stars with low neutron-capture element abundances and origins in ancient dwarf galaxies

ABSTRACT We present a detailed chemical abundance and kinematic analysis of six extremely metal-poor (−4.2 ≤ [Fe/H] ≤−2.9) halo stars with very low neutron-capture abundances ([Sr/H] and [Ba/H]) based on high-resolution Magellan/MIKE spectra. Three of our stars have [Sr/Ba] and [Sr/H] ratios that resemble those of metal-poor stars in ultra-faint dwarf galaxies (UFDs). Since early UFDs may be the building blocks of the Milky Way, extremely metal-poor halo stars with low, UFD-like Sr and Ba abundances may thus be ancient stars from the earliest small galactic systems that were accreted by the proto-Milky Way. We label these objects as Small Accreted Stellar System (SASS) stars, and we find an additional 61 similar ones in the literature. A kinematic analysis of our sample and literature stars reveals them to be fast-moving halo objects, all with retrograde motion, indicating an accretion origin. Because SASS stars are much brighter than typical UFD stars, identifying them offers promising ways towards detailed studies of early star formation environments. From the chemical abundances of SASS stars, it appears that the earliest accreted systems were likely enriched by a few supernovae whose light element yields varied from system to system. Neutron-capture elements were sparsely produced and/or diluted, with r-process nucleosynthesis playing a role. These insights offer a glimpse into the early formation of the Galaxy. Using neutron-capture elements as a distinguishing criterion for early formation, we have access to a unique metal-poor population that consists of the oldest stars in the universe.

The Redshift Evolution of the M •–M ⋆ Relation for JWST’s Supermassive Black Holes at z > 4

JWST has detected many overmassive galactic systems at z > 4, where the mass of the black hole, M •, is 10–100 times larger than expected from local relations, given the host’s stellar mass, M ⋆. This paper presents a model to describe these overmassive systems in the high-z Universe. We suggest that the black hole mass is the main driver of high-z star formation quenching. Supermassive black holes globally impact their high-z galaxies because their hosts are physically small, and the black holes have duty cycles close to unity at z > 4. In this regime, we assume that black hole mass growth is regulated by the quasar’s output, while stellar mass growth is quenched by it and uncorrelated to the global properties of the host halo. We find that the ratio M •/M ⋆ controls the average star formation efficiency: if M •/M ⋆ > 8 × 1018(nΛ/ f Edd)[(Ω b M h )/(Ω m M ⋆) − 1], then the galaxy is unable to form stars efficiently. Once this ratio exceeds the threshold, a runaway process brings the originally overmassive system toward the local M •–M ⋆ relation. Furthermore, the M •–M ⋆ relation evolves with redshift as ∝(1 + z)5/2. At z ∼ 5, we find an overmassive factor of ∼55, in excellent agreement with current JWST data and the high-z relation inferred from those. Extending the black hole horizon farther in redshift and lower in mass will test this model and improve our understanding of the early coevolution of black holes and galaxies.

Analytic solution of chemical evolution models with Type Ia supernovae

Context. In recent years, a significant number of works have been focussed on finding analytic solutions for the chemical enrichment models of galactic systems, including the Milky Way. Some of these solutions, however, are not able to account for the enrichment produced by Type Ia supernovae (SNe) due to the presence of the delay time distributions (DTDs) in the models. Aims. We present a new analytic solution for the chemical evolution model of the Galaxy. This solution can be used with different prescriptions of the DTD, including the single- and double-degenerate scenarios, and allows for the inclusion of an arbitrary number of pristine gas infalls. Methods. We integrated the chemical evolution model by extending the instantaneous recycling approximation with the contribution of Type Ia SNe. This implies an extra term in the modelling that depends on the DTD. For DTDs that lead to non-analytic integrals, we describe them as a superposition of Gaussian, exponential, and 1/t functions using a restricted least-squares fitting method. Results. We obtained the exact solution for a chemical model with Type Ia SNe widely used in previous works, while managing to avoid numerical integration errors. This solution is able to reproduce the expected chemical evolution of the α and iron-peak elements in less computing time than numerical integration methods. We compare the pattern in the [Si/Fe] versus [Fe/H] plane observed by APOGEE DR17 with that predicted by the model. We find the low α sequence can be explained by a delayed gas infall. We exploit the applicability of our solution by modelling the chemical evolution of a simulated Milky Way-like galaxy from its star formation history. The implementation of our solution has been released as a PYTHON package. Conclusions. Our solution constitutes a promising tool for Galactic archaeology studies and it is able to model the observed trends in α element abundances versus [Fe/H] in the solar neighbourhood. We infer the chemical information of a simulated galaxy modelled without chemistry.

The impact of baryonic potentials on the gravothermal evolution of self-interacting dark matter haloes

ABSTRACT The presence of a central baryonic potential can have a significant impact on the gravothermal evolution of self-interacting dark matter (SIDM) haloes. We extend a semi-analytical fluid model to incorporate the influence of a static baryonic potential and calibrate it using controlled N-body simulations. We construct benchmark scenarios with varying baryon concentrations and different SIDM models, including constant and velocity-dependent self-interacting cross-sections. The presence of the baryonic potential induces changes in SIDM halo properties, including central density, core size, and velocity dispersion, and it accelerates the halo’s evolution in both expansion and collapse phases. Furthermore, we observe a quasi-universality in the gravothermal evolution of SIDM haloes with the baryonic potential, resembling a previously known feature in the absence of the baryons. By appropriately rescaling the physical quantities that characterize the SIDM haloes, the evolution of all our benchmark cases exhibits remarkable similarity. Our findings offer a framework for testing SIDM predictions using observations of galactic systems where baryons play a significant dynamical role.

Spectral and polarization properties of reflected X-ray emission from black-hole accretion discs for a distant observer: the lamp-post model

ABSTRACT Rebirth of X-ray polarimetric instruments will have a significant impact on our knowledge of compact accreting sources. The properties of inner-accreting regions of active galactic nuclei (AGNs) or X-ray binary systems, such as black-hole spin, their disc inclination and orientation, shape, and size of their corona, can be polarimetrically studied, parallelly to the well-known X-ray spectroscopic and timing techniques. In this work, we provide a new spectropolarimetric numerical estimate of X-rays in the lamp-post coronal model for a distant observer, including a polarized reflected radiation from the accretion disc. The local disc reflection was simulated using the codes titan and stokes and includes variable disc ionization as well as Monte Carlo treatment of Compton multiple scatterings. We introduce a relativistic code kynstokes based on our well-tested KY package that accounts for all relativistic effects on radiation near a black hole, apart from the returning radiation, and adds a possibility of polarized coronal emission. We study the spectrum, polarization degree and polarization angle at spatial infinity for various global system parameters and we demonstrate the difference at infinity, if analytical local reflection computations are used. We newly predict that in the hard X-rays the reflected component can be 25 per cent polarized and the total emission can be 9 per cent polarized in the most favourable, yet realistic configurations of radio-quiet AGNs. Thus, the relativistic disc reflection remains important for the interpretation of X-ray polarimetric observations.

Investigating galactic double white dwarfs for the sub-mHz gravitational wave mission ASTROD-GW

A large number of galactic binary systems emit gravitational waves (GW) continuously with frequencies below $\ensuremath{\sim}10\text{ }\text{ }\mathrm{mHz}$. The LISA mission could identify tens of thousands of binaries over years of observation and will be subject to the confusion noise around 1 mHz yielded by the unresolved sources. Beyond LISA, there are several missions have been proposed to observe GWs in the sub-mHz range where the galactic foreground is expected to be overwhelming the instrumental noises. In this study, we investigate the detectability of sub-mHz GW missions to detect the galactic double white dwarf (DWD) binaries and evaluate the confusion noise produced by the undistinguished DWDs. This confusion noise could also be viewed as a stochastic GW foreground and be effectively observed in the sub-mHz band. The parameter determinations for the modeled foreground are examined by employing different detector sensitivities and population models. By assuming the determined foregrounds could be subtracted from the data, we evaluate the residuals which are expected to have power spectral densities two orders of magnitude lower than the originals data.

GaiaData Release 3

Context.GaiaData Release 3 (DR3) provides a number of new data products that complement the early DR3 made available two years ago. Among these is the firstGaiacatalogue of eclipsing-binary candidates containing 2 184 477 sources with brightnesses from a few magnitudes to 20 mag in theGaiaG-band and covering the full sky.Aims.We present the catalogue, describe its content, provide tips for its use, estimate its quality, and show illustrative samples.Methods.Candidate selection is based on the results of variable object classification performed within theGaiaData Processing and Analysis Consortium. Candidates are then further filtered using eclipsing-binary-tailored criteria based on theG-band light curves. To find the orbital period, a large ensemble of trial periods is first acquired using three distinct period-search methods applied to the cleanedGlight curve of each source. TheGlight curve is then modelled with up to two Gaussians and a cosine for each trial period. The best combination of orbital period and geometric model is finally selected using Bayesian model comparison based on the BIC. A global ranking metric is provided to rank the quality of the chosen model between sources. The catalogue is restricted to orbital periods larger than 0.2 days.Results.Of ∼600 000 available crossmatches, about 530 000 of the candidates are classified as eclipsing binaries in the literature as well, and 93% of them have published periods compatible with theGaiaperiods. Catalogue completeness is estimated to be between 25% and 50%, depending on the sky region, relative to the OGLE4 catalogues of eclipsing binaries towards the Galactic Bulge and the Magellanic Clouds. The analysis of an illustrative sample of ∼400 000 candidates with significant parallaxes shows properties in the observational Hertzsprung-Russell diagram as expected for eclipsing binaries. The subsequent analysis of a subsample of detached bright candidates provides further hints for the exploitation of the catalogue. We also address the observed lack of short-period (less than a day) systems in the Magellanic Cloud in comparison to Galactic systems. The orbital periods, light-curve model parameters, and global rankings are all published in the catalogue with their related uncertainties where applicable.Conclusions.ThisGaiaDR3 catalogue of eclipsing-binary candidates constitutes the largest catalogue to date in terms of number of sources, sky coverage, and magnitude range.

Radiation Emission and Absorption by Astrophysical Jets from XRBs

High-energy particles and radiation such as protons, pions, muons, neutrinos and gamma-ray photons are known to emit from collimated outflows of magnetized astrophysical plasma known as jets. They are being ejected by Active Galactic Nuclei or X-ray binary systems consisting of a companion star accreting mass onto a black hole or a neutron star. Our work focuses on the calculation of the produced particle energy distributions and mainly on the intensity of gamma-rays. We apply this model to prominent examples of Black Hole XRBs such as Cygnus X-1, SS 433, GRS 1915+105, etc. We also study and work on calculating the absorption of emitted gamma-rays by soft and hard X-ray radiation fields originating from the system’s accretion disc, corona and companion star that could strongly affect the jet’s gamma-ray intensity that finally reaches the Earth.

From a Dinuclear System to Close Binary Cosmic Objects

Applying the ideas from microscopic objects to macroscopic stellar and galactic systems, the evolution of compact di-stars and di-galaxies is studied in the mass asymmetry coordinate. The formation of stable binary systems is analyzed. The role of symmetrization of an initially asymmetric binary system is revealed in the transformation of gravitational energy into internal energy of stars or galaxies accompanied by the release of a huge amount of energy. For the contact binary stars, the change of the orbital period is explained by evolution to symmetry in mass asymmetry coordinates. The matter transfer in binary black holes is studied. The conditions for the merger of black holes in a binary system are analyzed regarding the radiation of gravitational waves. Using the model based on the Regge-like laws, the Darwin instability effect in binary systems is discussed. New analytical formulas are derived for the period of orbital rotation and the relative distance between the components of a binary system. The impossibility of the appearance of a binary cosmic object from a single cosmic object is revealed.

A Special State Transition in the Blazar OT 081: Implication for the Unified State Transition Paradigm of Different-scale Black Hole Systems

Exploring the similar behavior of black hole systems with different scales will provide insight into the unified physical laws of black hole systems. Here, we report a special spectral state transition process in blazar OT 081, which is driven by a supermassive black hole. Based on the X-ray flux versus photon index distribution, stacked average X-ray spectra, L 2500Å − α ox parameters, and broadband spectral energy distributions, we conducted a pilot study on the spectral state transition for this blazar system. The analysis results suggest that the source transitions through a steady-state point from a strong corona, weak jet unsteady state to a weak corona, strong jet unsteady state. This state transition behavior is similar to that observed in the much smaller Galactic black hole binary systems, providing further observational evidence for the possible existence of a unified state transition paradigm for black hole systems with different scales.

A giant umbrella-like stellar stream around the tidal ring galaxy NGC 922

Context. Tidal ring galaxies (TRGs) are rarely observed in the local universe due to their intrinsically transient nature. The tidal ring structures are the result of strong interactions between gas-rich stellar disks and smaller galactic systems, and do not last longer than ∼500 Myr. Therefore, these are perfect scenarios in which to look for the debris of recently accreted dwarf galactic systems. Aims. Our goal is to study the low surface brightness stellar structures around the TRG NGC 922 and to revise the hypothesis of its formation in light of these new data. Methods. We present new deep images of the TRG NGC 922 and its surroundings from the DESI Legacy survey data and from our observations with an amateur telescope. These observations are compared with results from high-resolution N-body simulations that were designed to reproduce an alternative formation scenario for this peculiar galaxy. Results. Our new observations unveil that the low surface brightness stellar tidal structures around NGC 922 are much more complex than reported in previous works. In particular, the formerly detected tidal spike-like structure at the north-east of the central galaxy disk is not connected with the dwarf companion galaxy PGC 3080368, which has been suggested as the intruder triggering the ring formation of NGC 922. The deep images reveal that this tidal structure mainly has a fainter giant umbrella-like shape, and thus it was formed from the tidal disruption of a different satellite. Using the broad-band g, r, and z DESI LS images, we measured the photometric properties of this stellar stream, estimating a total absolute magnitude in the r band of Mr = −17.0 ± 0.03 mag and a total stellar mass for the stream of between 6.9 and 8.5×108 M⊙. We performed a set of N-body simulations to reproduce the observed NGC922-intruder interaction, suggesting a new scenario for the formation of its tidal ring from the infall of a gas-rich satellite around 150 Myr ago. Finally, our deep images also reveal a tidal shell around the dwarf galaxy PGC 3080368, a possible fossil of a recent merger with a smaller satellite, which may suggest it is in its first infall towards NGC 922.

A Semianalytical Line Transfer Model. III. Galactic Inflows

We present calculations of ultraviolet spectra resulting from the scattering of photons by gas infalling onto an isotropically emitting source of radiation. The model is based on an adaptation of the semianalytical line transfer (SALT) code of Scarlata & Panagia, and designed to interpret the inverse P Cygni profiles observed in the spectra of partially ionized galactic inflows. In addition to presenting the model, we explore the parameter space of the inflowing SALT model and recreate various physically motivated scenarios including spherical inflows, inflows with covering fractions less than unity, and galactic fountains (i.e., galactic systems with both an inflowing and outflowing component). The resulting spectra from inflowing gas show spectral features that could be misinterpreted as interstellar medium features in low resolution spectroscopy (σ ≈ 120 km s−1), suggesting that the total number of galactic systems with inflows is undercounted. Our models suggest that observations at medium resolution (R = 6000 or σ ≈ 50 km s−1) that can be obtained with 8 m class telescopes will be able to resolve the characteristic inverse P Cygni profiles necessary to identify inflows.

Galactic satellite systems in CDM, WDM and SIDM

ABSTRACT We investigate the population of bright satellites ($M_{*} \ge 10^{5} \, \mathrm{M}_{\odot }$) of haloes of mass comparable to that of the Milky Way in cosmological simulations in which the dark matter (DM) is either cold, warm, or self-interacting (CDM, WDM, and SIDM, respectively). The nature of the DM gives rise to differences in the abundance and structural properties of field haloes. In WDM, the main feature is a reduction in the total number of galaxies that form, reflecting a suppression of low-mass DM haloes and lower galaxy formation efficiency compared to CDM. For SIDM, the changes are structural, restricted to the central regions of haloes and dependent on the assumed self-interaction cross-section. We also consider different baryonic subgrid physics models for galaxy formation, in which supernova gas blowouts can or cannot induce the formation of a core in dwarf galaxies. Overall, the inclusion of baryons lessen the differences in the halo properties in the different DM models compared to DM-only simulations. This affects the satellite properties at infall and therefore their subsequent tidal stripping and survival rates. None the less, we find slightly less concentrated satellite radial distributions as the SIDM cross-section increases. Unfortunately, we also find that the satellite populations in simulations with baryon-induced cores in CDM and WDM can mimic the results found in SIDM, making the satellite stellar mass and maximum circular velocity functions heavily degenerate on the assumed nature of the DM and the adopted subgrid modelling. These degeneracies preclude using the brightest satellites of the Milky Way to constrain the nature of DM.

Broader view of bimetric MOND

All existing treatments of bimetric MOND (BIMOND) -- a class of relativistic versions of MOND -- have dealt with a rather restricted sub-class: The Lagrangian of the interaction between the gravitational degrees of freedom -- the two metrics -- is a function of a certain {\it single} scalar argument built from the difference in connections of the two metrics. I show that the scope of BIMOND is much richer: The two metrics can couple through several scalars to give theories that all have a "good" nonrelativistic (NR) limit -- one that accounts correctly, a-la MOND, for the dynamics of galactic systems, {\it including gravitational lensing}. This extended-BIMOND framework exhibits a qualitative departure from the way we think of MOND at present, as encapsulated, in all its aspects, by one "interpolating function" of one acceleration variable. After deriving the general field equations, I pinpoint the subclass of theories that satisfy the pivotal requirement of a good NR limit. These involve three independent, quadratic scalar variables. In the NR limit these scalars all reduce to the same acceleration scalar, and the NR theory then does hinge on one function of a {\it a single} acceleration variable -- representing the NR MOND "interpolating function", whose form is largely dictated by the observed NR galactic dynamics. However, these scalars behave differently, in different relativistic contexts. So, the full richness of the multi-variable Lagrangian, as it enters cosmology, for example, is hardly informed by what we learn from observations of galactic dynamics. In this paper, I present the formalism, with some generic examples. I also consider some cosmological solutions where the two metrics are small departures from one Friedman-Lemaitre-Robertson-Walker metric. This may offer a framework for describing cosmology within the extended BIMOND.