Search

We present cosmological results from the measurement of clustering of galaxy, quasar and Lyman-$\alpha$ forest tracers from the first year of observations with the Dark Energy Spectroscopic Instrument (DESI Data Release 1). We adopt the full-shape (FS) modeling of the power spectrum, including the effects of redshift-space distortions, in an analysis which has been validated in a series of supporting papers. In the flat $\Lambda$CDM cosmological model, DESI (FS+BAO), combined with a baryon density prior from Big Bang Nucleosynthesis and a weak prior on the scalar spectral index, determines matter density to $\Omega_\mathrm{m}=0.2962\pm 0.0095$, and the amplitude of mass fluctuations to $\sigma_8=0.842\pm 0.034$. The addition of the cosmic microwave background (CMB) data tightens these constraints to $\Omega_\mathrm{m}=0.3056\pm 0.0049$ and $\sigma_8=0.8121\pm 0.0053$, while further addition of the the joint clustering and lensing analysis from the Dark Energy Survey Year-3 (DESY3) data leads to a 0.4% determination of the Hubble constant, $H_0 = (68.40\pm 0.27)\,{\rm km\,s^{-1}\,Mpc^{-1}}$. In models with a time-varying dark energy equation of state, combinations of DESI (FS+BAO) with CMB and type Ia supernovae continue to show the preference, previously found in the DESI DR1 BAO analysis, for $w_0>-1$ and $w_a<0$ with similar levels of significance. DESI data, in combination with the CMB, impose the upper limits on the sum of the neutrino masses of $\sum m_\nu < 0.071\,{\rm eV}$ at 95% confidence. DESI data alone measure the modified-gravity parameter that controls the clustering of massive particles, $\mu_0=0.11^{+0.45}_{-0.54}$, while the combination of DESI with the CMB and the clustering and lensing analysis from DESY3 constrains both modified-gravity parameters, giving $\mu_0 = 0.04\pm 0.22$ and $\Sigma_0 = 0.044\pm 0.047$, in agreement with general relativity. [Abridged.], Comment: This DESI Collaboration Key Publication is part of the 2024 publication series using the first year of observations (see https://data.desi.lbl.gov/doc/papers/). 55 pages, 10 figures

We present the samples of galaxies and quasars used for DESI 2024 cosmological analyses, drawn from the DESI Data Release 1 (DR1). We describe the construction of large-scale structure (LSS) catalogs from these samples, which include matched sets of synthetic reference `randoms' and weights that account for variations in the observed density of the samples due to experimental design and varying instrument performance. We detail how we correct for variations in observational completeness, the input `target' densities due to imaging systematics, and the ability to confidently measure redshifts from DESI spectra. We then summarize how remaining uncertainties in the corrections can be translated to systematic uncertainties for particular analyses. We describe the weights added to maximize the signal-to-noise of DESI DR1 2-point clustering measurements. We detail measurement pipelines applied to the LSS catalogs that obtain 2-point clustering measurements in configuration and Fourier space. The resulting 2-point measurements depend on window functions and normalization constraints particular to each sample, and we present the corrections required to match models to the data. We compare the configuration- and Fourier-space 2-point clustering of the data samples to that recovered from simulations of DESI DR1 and find they are, generally, in statistical agreement to within 2\% in the inferred real-space over-density field. The LSS catalogs, 2-point measurements, and their covariance matrices will be released publicly with DESI DR1., Comment: This DESI Collaboration Key Publication is part of the 2024 publication series using the first year of observations (see https://data.desi.lbl.gov/doc/papers/)

We present the measurements and cosmological implications of the galaxy two-point clustering using over 4.7 million unique galaxy and quasar redshifts in the range $0.1

The Dark Energy Spectroscopic Instrument (DESI) survey is mapping the large-scale distribution of millions of Emission Line Galaxies (ELGs) over vast cosmic volumes to measure the growth history of the Universe. However, compared to Luminous Red Galaxies (LRGs), very little is known about the connection of ELGs with the underlying matter field. In this paper, we employ a novel theoretical model, SHAMe-SF, to infer the connection between ELGs and their host dark matter subhaloes. SHAMe-SF is a version of subhalo abundance matching that incorporates prescriptions for multiple processes, including star formation, tidal stripping, environmental correlations, and quenching. We analyse the public measurements of the projected and redshift-space ELGs correlation functions at $z=1.0$ and $z=1.3$ from DESI One Percent data release, which we fit over a broad range of scales $r \in [0.1, 30]/h^{-1}$Mpc to within the statistical uncertainties of the data. We also validate the inference pipeline using two mock DESI ELG catalogues built from hydrodynamical (TNG300) and semi-analytical galaxy formation models (\texttt{L-Galaxies}). SHAMe-SF is able to reproduce the clustering of DESI-ELGs and the mock DESI samples within statistical uncertainties. We infer that DESI ELGs typically reside in haloes of $\sim 10^{11.8}h^{-1}$M$_{\odot}$ when they are central, and $\sim 10^{12.5}h^{-1}$M$_{\odot}$ when they are a satellite, which occurs in $\sim$30 \% of the cases. In addition, compared to the distribution of dark matter within halos, satellite ELGs preferentially reside both in the outskirts and inside haloes, and have a net infall velocity towards the centre. Finally, our results show evidence of assembly bias and conformity., Comment: 23 pages, 17 figures. To be submitted to A&A

The MUltiplexed Survey Telescope (MUST) is a 6.5-meter telescope under development. Dedicated to highly-multiplexed, wide-field spectroscopic surveys, MUST observes over 20,000 targets simultaneously using 6.2-mm pitch positioning robots within a ~5 deg2 field of view. MUST aims to carry out the first Stage-V spectroscopic survey in the 2030s to map the 3D Universe with over 100 million galaxies and quasars, spanning from the nearby Universe to redshift z~5.5, corresponding to around 1 billion years after the Big Bang. To cover this extensive redshift range, we present an initial conceptual target selection algorithm for different types of galaxies, from local bright galaxies, luminous red galaxies, and emission line galaxies to high-redshift (2 < z < 5.5) Lyman-break galaxies. Using Fisher forecasts, we demonstrate that MUST can address fundamental questions in cosmology, including the nature of dark energy, test of gravity theories, and investigations into primordial physics. This is the first paper in the series of science white papers for MUST, with subsequent developments focusing on additional scientific cases such as galaxy and quasar evolution, Milky Way physics, and dynamic phenomena in the time-domain Universe., Comment: To be submitted to SCPMA

The theory of interfacial properties in liquid-liquid or liquid-vapour systems is nearly 200 years old. The advent of computational tools has greatly advanced the field, mainly through the use of Molecular Dynamics simulations. Despite the successes and advances in the theory of interfacial phenomena for liquid-liquid systems, the study of solid-liquid interfaces remains a challenge both theoretically and experimentally. The main reason why the treatment of solid-liquid systems has fallen behind that of liquid-liquid systems is that there are complications that arise whenever an interface involving solid systems is considered involving both theory of the solid-liquid interface and the calculations using MD simulations. An example of the former is that, contrary to the liquid-liquid case, the interfacial properties of solids depend on the lattice orientation. The main complications in these calculations arise from the fact that for solids the ``mechanical route'' cannot be used. To overcome this problem, several numerical approaches were proposed. The main purpose of this review is to provide an overview of these different methodologies and to discuss their strengths and weaknesses. We classify these methodologies into two main groups: direct and indirect methods. Direct methods are those that can calculate directly the properties of interfaces, while in indirect approaches the properties of the interface are not the primary result of the simulations. We also included a discussion on the origin of the difficulties in considering solid interfaces from a thermodynamic point of view. In the second part of the review, we discuss two key related topics: nucleation theory and curved interfaces. They both represent an important problem in the study of interfaces and in the context of solid-liquid ones for which the research is still extremely active., Comment: 111 pages, 14 figures

Interoperability is the key to implementing a wide range of energy systems applications. It involves the seamless cooperation of different methods and components. With smart energy systems, interoperability faces challenges due to integrating differ-ent approaches and technologies. This includes dealing with heterogeneous approaches with various communication proto-cols and data formats. However, it is essential for smart energy systems to carry out thorough interoperability tests. They are usually diverse, and challenging, thus requiring careful consideration of compatibility issues and complex integration scenari-os. Overcoming these challenges requires a systematic approach that includes thorough test planning, rigorous testing, and continuous test monitoring. Although numerous testing approaches exist, most are more developed at the component/device level than at the system level. Consequently, there are few approaches and related facilities to test the interoperability of smart energy approaches and solutions at the system level. This work analyses existing interoperability test concepts, identi-fies enablers and the potential for harmonisation of procedures, and proposes further developments of these approaches., Comment: 14th Mediterranean Conference on Power Generation Transmission, Distribution and Energy Conversion (MED POWER 2024)

Wireless communications are critical in the constantly changing environment of IoT and RFID technologies, where thousands of devices can be deployed across a wide range of scenarios. Whether connecting to cloud servers or local fog/edge devices, maintaining seamless communications is difficult, especially in demanding contexts like industrial warehouses or remote rural areas. Opportunistic networks, when combined with edge devices, provide a possible solution to this challenge. These networks enable IoT devices, particularly mobile devices, to redirect information as it passes via other devices until it reaches an edge node. Using different communication protocols, this paper investigates their effects on response times and total messages received for a opportunistic assisted RFID system. Specifically, this article compares two communications technologies (Bluetooth 5 and Wize) when used for building a novel Opportunistic Edge Computing (OEC) identification system based on low-cost Single-Board Computers (SBCs). For such a comparison, measurements have been performed for quantifying packet loss and latency. The tests consisted in two experiments under identical conditions and scenarios, with a node located roadside, transmitting identification information, and a node located inside a moving vehicle that was driven at varying vehicle speeds. The obtained results show for Bluetooth 5 average latencies ranging between 700 and 950 ms with packet losses between 7% and 27%, whereas for Wize the average delay as between 150 and 370 ms with packet losses between 20% and 52%., Comment: Paper accepted in IEEE RFID 2024

Obtaining well-calibrated photometric redshift probability densities for galaxies without a spectroscopic measurement remains a challenge. Deep learning discriminative models, typically fed with multi-band galaxy images, can produce outputs that mimic probability densities and achieve state-of-the-art accuracy. However, such models may be affected by miscalibration that would result in discrepancies between the model outputs and the actual distributions of true redshifts. Our work develops a novel method called the Contrastive Learning and Adaptive KNN for Photometric Redshift (CLAP) that resolves this issue. It leverages supervised contrastive learning (SCL) and k-nearest neighbours (KNN) to construct and calibrate raw probability density estimates, and implements a refitting procedure to resume end-to-end discriminative models ready to produce final estimates for large-scale imaging data. The harmonic mean is adopted to combine an ensemble of estimates from multiple realisations for improving accuracy. Our experiments demonstrate that CLAP takes advantage of both deep learning and KNN, outperforming benchmark methods on the calibration of probability density estimates and retaining high accuracy and computational efficiency. With reference to CLAP, we point out that miscalibration is particularly sensitive to the method-induced excessive correlations among data instances in addition to the unaccounted-for epistemic uncertainties. Reducing the uncertainties may not guarantee the removal of miscalibration due to the presence of such excessive correlations, yet this is a problem for conventional deep learning methods rather than CLAP. These discussions underscore the robustness of CLAP for obtaining photometric redshift probability densities required by astrophysical and cosmological applications. This is the first paper in our series on CLAP., Comment: 22 + 6 pages, 9 + 5 figures

The large-scale linear halo bias encodes the relation between the clustering of dark-matter (DM) halos and that of the underlying matter density field. Although the primary dependence of bias on halo mass is well understood in the context of structure formation, the physical origins of the multiple additional relations at fixed halo mass, commonly known as secondary halo bias, have not been fully elucidated. Of particular relevance is the secondary dependence on halo assembly history, known as halo assembly bias. Our goal is to determine whether the properties of the initial regions from which $z=0$ halos originate produce any secondary bias at $z=0$. By analyzing these initial dependencies in connection with halo assembly bias, we intend to provide insight on the physical origins of the effect. To this end, we select halos at $z=0$ in the IllustrisTNG DM-only simulation and trace back the positions and velocities of their DM particles to $z=12$. The resulting initial regions are characterized according to several shape-related and kinematic properties. The secondary bias signal produced by these properties at $z=0$ is measured using an object-by-object bias estimator, which offers significant analytical advantages. We show that, when split by the properties of their initial DM clouds, $z=0$ halos display significant secondary bias, clearly exceeding the amplitude of the well-known halo assembly bias signal produced by concentration and age. The maximum bias segregation is measured for cloud velocity dispersion and radial velocity, followed by cloud concentration, sphericity, ellipticity and triaxiality. We further show that both velocity dispersion and radial velocity are also the properties of the initial clouds that most strongly correlate with halo age and concentration at fixed halo mass. Our results highlight the importance of linear effects in shaping halo assembly bias., Comment: 9 pages, 5 figures. Submitted to A&A, comments welcome

The connection between galaxies and dark matter halos encompasses a range of processes and play a pivotal role in our understanding of galaxy formation and evolution. Traditionally, this link has been established through physical or empirical models. Machine learning techniques are adaptable tools that handle high-dimensional data and grasp associations between numerous attributes. In particular, probabilistic models capture the stochasticity inherent to these complex relations. We compare different probabilistic machine learning methods to model the uncertainty in the halo-galaxy connection and efficiently generate galaxy catalogs that faithfully resemble the reference sample by predicting joint distributions of central galaxy properties conditioned to their host halo features. The analysis is based on the IllustrisTNG300 simulation. The methods model the distributions in different ways. We compare a multilayer perceptron that predicts the parameters of a multivariate Gaussian distribution, a multilayer perceptron classifier, and the method of normalizing flows. The classifier predicts the parameters of a Categorical distribution, which are defined in a high-dimensional parameter space through a Voronoi cell-based hierarchical scheme. We evaluate the model's performances under various sample selections based on halo properties. The three methods exhibit comparable results, with normalizing flows showing the best performance in most scenarios. The models reproduce the main features of galaxy properties distributions with high-fidelity and reproduce the results obtained with traditional, deterministic, estimators. Our results also indicate that different halos and galaxy populations are subject to varying degrees of stochasticity, which has relevant implications for studies of large-scale structure., Comment: 14 pages, 8 figures, 6 tables

Magnetite, a naturally abundant mineral, frequently interacts with water in both natural settings and various technical applications, making the study of its surface chemistry highly relevant. In this work, we investigate the hydrogen bonding dynamics and the presence of hydroxyl species at the magnetite-water interface using a combination of neural network potential-based molecular dynamics simulations and sum frequency generation vibrational spectroscopy. Our simulations, which involved large water systems, allowed us to identify distinct interfacial species, such as dissociated hydrogen and hydroxide ions formed by water dissociation. Notably, water molecules near the interface exhibited a preference for dipole orientation towards the surface, with bulk-like water behavior only re-emerging beyond 60 {\AA} from the surface. The vibrational spectroscopy results aligned well with the simulations, confirming the presence of a hydrogen bond network in the surface ad-layers. The analysis revealed that surface-adsorbed hydroxyl groups orient their hydrogen atoms towards the water bulk. In contrast, hydrogen-bonded water molecules align with their hydrogen atoms pointing towards the magnetite surface.

Soft robots have the potential to revolutionize the use of robotic systems with their capability of establishing safe, robust, and adaptable interactions with their environment, but their precise control remains challenging. In contrast, traditional rigid robots offer high accuracy and repeatability but lack the flexibility of soft robots. We argue that combining these characteristics in a hybrid robotic platform can significantly enhance overall capabilities. This work presents a novel hybrid robotic platform that integrates a rigid manipulator with a fully developed soft arm. This system is equipped with the intelligence necessary to perform flexible and generalizable tasks through imitation learning autonomously. The physical softness and machine learning enable our platform to achieve highly generalizable skills, while the rigid components ensure precision and repeatability., Comment: Corrected missing citation

Air-vented ionization chambers exposed to clinical radiation beams may suffer from recombination during the drift of the charge carriers towards the electrodes. Thus, dosimetry protocols recommend the use of a correction factor, usually denominated saturation factor ($k_{\rm sat}$), to correct the ionization chamber readout for the incomplete collection of charge. The two-voltage method is the recommended methodology for the calculation of the saturation factor, however, it is based on the early Boag model, which only takes into account the presence of positive and negative ions in the ionization chamber and does not account for the electric field screening or the free electron contribution to the signal. The numerical simulation shows a better agreement with the experimental data than the current analytical theories in terms of charge collection efficiency. The classical two-voltage method, systematically overestimates the saturation factor, with differences increasing with dose per pulse also present at low dose per pulse. These results may have implications for the dosimetry with ionization chambers in therapy modalities that use a dose per pulse higher than conventional radiotherapy such as intraoperative radiotherapy but also in conventional dose per pulse for ionization chambers that suffer from significant charge recombination.

Background and purpose: IMPT faces challenges in lung cancer treatment, like maintaining plan robustness for moving tumors against setup, range errors, and interplay effects. Proton Arc Therapy (PAT) is an alternative to maintain target coverage, potentially improving organ at risk (OAR) sparing, reducing beam delivery time (BDT), and enhancing patient experience. We aim to perform a systematic plan comparison study between IMPT and ELSA-PAT to assess its potential for lung cancer treatment. Material and Methods: 14 Lung ELSA-PAT plans were compared retrospectively with IMPT plans. 4D worst-case minimax robust optimization was performed, including 84 scenarios (3%,3 mm). Dosimetry assessment included target (CTV) and important OARs, on nominal and worst-case scenarios. Most relevant normal tissue complication probabilities (NTCP), target coverage robustness against interplay effect and beam delivery time (BDT) were evaluated. Results: CTV D95% and D98% showed no significant difference in comparison. PAT demonstrated better conformality by 66% (p = 0.00012) but delivered a higher heart mean dose (HMD,23%). There was a 2% increase in NTCP 2-year mortality risk with PAT. Total BDT was comparable among techniques. IMPT was more robust than PAT against interplay effect, considering both D1% (1,0 $\pm$ 0.8 Gy vs 1.1 $\pm$ 1.4 Gy) and D98% bandwidths (0.9$\pm$0.9 Gy vs 1.1 $\pm$ 1.3 Gy). Interpretation: both techniques provide a similar level of dose coverage to the target volume. Although PAT improved dose conformality, higher HMD translated into increased heart toxicity, presumably due to chosen planning methodology and OAR proximity to target. Increased energy layers and spots raised PAT beam delivery time, although it could improve daily treatment workflow.

The tumultuous effects of ultraviolet photons that source cosmic reionization, the subsequent compression and shock-heating of low-density regions, and the modulation of baryons in shallow potential wells induced by the passage of ionization fronts, collectively introduce perturbations to the evolution of the intergalactic medium in the post-reionization era. These enduring fluctuations persist deep into the post-reionization era, casting a challenge upon precision cosmology endeavors targeting tracers in this cosmic era. Simultaneously, these relics from reionization also present a unique opportunity to glean insights into the astrophysics that govern the epoch of reionization. In this work, we propose a first study of the cross-correlation of \lya forest and 21 cm intensity mapping, accounting for the repercussions of inhomogeneous reionization in the post-reionization era. We investigate the ability of SKA $\times$ DESI-like, SKA $\times$ MUST-like, and PUMA $\times$ MUST-like instrumental setups to achieve a high signal-to-noise ratio (SNR) in the redshift range $3.5 \leq z \leq 4$. Moreover, we assess how alterations in integration time, survey area, and reionization scenarios impact the SNR. Furthermore, we forecast the cross-correlation's potential to constrain cosmological parameters under varying assumptions: considering or disregarding reionization relics, marginalizing over reionization astrophysics, and assuming perfect knowledge of reionization. Notably, our findings underscore the remarkable capability of a futuristic PUMA $\times$ MUST-like setup, with a modest 100-hour integration time over a 100 sq. deg. survey, to constrain the ionization efficiency error to $\sigma_\zeta = 3.42 $., Comment: Comments welcome! (16 pages, 10 figures)

We compute baryon and electric charge deposition in high-energy heavy-ion collisions using the Color Glass Condensate (CGC) Effective Field Theory, where at leading order charge is deposited through multiple scatterings of valence quarks with a saturated gluon target. A simplified phenomenological formula is derived to describe charge deposition, from which the parametrical dependence with collisional energy and geometry can be extracted. We present an approximate analytical prediction of the so-called baryon stopping parameter $\alpha_B$, which shows excellent agreement with the state-of-the art extractions of $\alpha_B$ from experimental data. These results are further validated using the McDIPPER framework, by computing charge deposition at midrapidity across a range of collision energies ($\sqrt{s_{\rm NN}}= 62.4 - 5020$ GeV).

On large scales, measurements of the Lyman-$\alpha$ forest offer insights into the expansion history of the Universe, while on small scales, these impose strict constraints on the growth history, the nature of dark matter, and the sum of neutrino masses. This work introduces ForestFlow, a cosmological emulator designed to bridge the gap between large- and small-scale Lyman-$\alpha$ forest analyses. Using conditional normalizing flows, ForestFlow emulates the 2 Lyman-$\alpha$ linear biases ($b_\delta$ and $b_\eta$) and 6 parameters describing small-scale deviations of the 3D flux power spectrum ($P_\mathrm{3D}$) from linear theory. These 8 parameters are modeled as a function of cosmology $\unicode{x2013}$ the small-scale amplitude and slope of the linear power spectrum $\unicode{x2013}$ and the physics of the intergalactic medium. Thus, in combination with a Boltzmann solver, ForestFlow can predict $P_\mathrm{3D}$ on arbitrarily large (linear) scales and the 1D flux power spectrum ($P_\mathrm{1D}$) $\unicode{x2013}$ the primary observable for small-scale analyses $\unicode{x2013}$ without the need for interpolation or extrapolation. Consequently, ForestFlow enables for the first time multiscale analyses. Trained on a suite of 30 fixed-and-paired cosmological hydrodynamical simulations spanning redshifts from $z=2$ to $4.5$, ForestFlow achieves $3$ and $1.5\%$ precision in describing $P_\mathrm{3D}$ and $P_\mathrm{1D}$ from linear scales to $k=5\,\mathrm{Mpc}^{-1}$ and $k_\parallel=4\,\mathrm{Mpc}^{-1}$, respectively. Thanks to its parameterization, the precision of the emulator is also similar for both ionization histories and two extensions to the $\Lambda$CDM model $\unicode{x2013}$ massive neutrinos and curvature $\unicode{x2013}$ not included in the training set. ForestFlow will be crucial for the cosmological analysis of Lyman-$\alpha$ forest measurements from the DESI survey., Comment: 17 pages, 11 figures. Submitted to A&A

This contribution summarizes the main activities and objectives of the outreach project Astroaccesible, whose main aim is to carry the teaching and diffusion of astronomy among all kinds of collectives, focusing on blind and visually impaired (BVI) people. This project is led by a blind astronomer and aims to use a variety of resources based on different sensory channels, avoiding limiting the transmission of concepts to visual perception. This principle favors inclusion and benefits everyone, as the information is not presented using just one channel. This strategy is especially convenient for the nowadays typical data acquisition, where a variety of sources of information, not solely based on the collection of different spectral domains of electromagnetic radiation, is used. Moreover, the study of new multi-messenger astronomy could be much better understood using a multi-messenger teaching approach, favoring inclusion, motivation, and creativity., Comment: To appear as a proceeding of the 13th Cosmic-Ray International Studies and Multi-messenger Astroparticle Conference held on June 2024 in Trapani (Sicily, Italy). 8 pages and 4 figures

The standard Poisson-Boltzmann model for molecular electrostatics assumes a sharp variation of the permittivity and salt concentration along the solute-solvent interface. The discontinuous field parameters are not only difficult numerically, but also are not a realistic physical picture, as it forces the dielectric constant and ionic strength of bulk in the near-solute region. An alternative to alleviate some of these issues is to represent the molecular surface as a diffuse interface, however, this also presents challenges. In this work we analysed the impact of the shape of the interfacial variation of the field parameters in solvation and binding energy. However we used a hyperbolic tangent function ($\tanh(k_p x)$) to couple the internal and external regions, our analysis is valid for other definitions. Our methodology was based on a coupled finite element (FEM) and boundary element (BEM) scheme that allowed us to have a special treatment of the permittivity and ionic strength in a bounded FEM region near the interface, while maintaining BEM elsewhere. Our results suggest that the shape of the function (represented by $k_p$) has a large impact on solvation and binding energy. We saw that high values of $k_p$ induce a high gradient on the interface, to the limit of recovering the sharp jump when $k_p\to\infty$, presenting a numerical challenge where careful meshing is key. Using the FreeSolv database to compare with molecular dynamics, our calculations indicate that an optimal value of $k_p$ for solvation energies was around 3. However, more challenging binding free energy tests make this conclusion more difficult, as binding showed to be very sensitive to small variations of $k_p$. In that case, optimal values of $k_p$ ranged from 2 to 20.

The magnetite/water interface is commonly found in nature and plays a crucial role in various technological applications. However, our understanding of its structural and dynamical properties at the molecular scale remains still limited. In this study, we develop an efficient Behler-Parrinello neural network potential (NNP) for the magnetite/water system, paying particular attention to the accurate generation of reference data with density functional theory. Using this NNP, we performed extensive molecular dynamics simulations of the magnetite (001) surface across a wide range of water coverages, from the single molecule to bulk water. Our simulations revealed several new ground states of low coverage water on the Subsurface Cation Vacancy (SCV) model and yielded a density profile of water at the surface that exhibits marked layering. By calculating mean square displacements, we obtained quantitative information on the diffusion of water molecules on the SCV for different coverages, revealing significant anisotropy. Additionally, our simulations provided qualitative insights into the dissociation mechanisms of water molecules at the surface.

Inspired by Bourqui's work on anticanonical height zeta functions on Hirzebruch surfaces, we study height zeta functions of split toric varieties with Picard rank 2 over global function fields, with respect to height functions associated with big metrized line bundles. We show that these varieties can be naturally decomposed into a finite disjoint union of subvarieties, where precise analytic properties of the corresponding height zeta functions can be given. As application, we obtain asymptotic formulas for the number of rational points of large height on each subvariety, with explicit leading constants and controlled error terms., Comment: Updated version including asymptotic formulas for the number of rational points of large height. Comments are welcome!

The calculation of the interfacial free energy between two thermodynamic phases is crucial across different fields including materials science, chemistry, and condensed matter physics. This study delves into a thermodynamic approach to predict interfacial free energy along a coexistence line. For that purpose, we utilize Gibbs-Cahn integration method to accurately determine the interfacial free energy along various coexistence conditions, relying solely on data obtained from a single-state interfacial free energy calculation at given conditions of pressure and temperature. Unlike traditional methods that require individual measurements at each thermodynamic state, this approach significantly reduces computational demand and enhances efficiency. The integration scheme requires the computation of the excess interfacial energy (by means of unbiased NVT computer simulations in which the two phases involved coexist) at different points along the coexistence line, for later use this information as an input in our calculations. We apply this method for the Lennard-Jones and mW water model for the case of the liquid-solid interface, as well as Lennard-Jones and TIP4P/2005 models for liquid-vapor interfaces. Our results highlight the accuracy and effectiveness of this thermodynamic integration route for estimating the interfacial free energy along a coexistence line., Comment: Main text: 13 pages and 9 figures. Supporting material: 7 pages and 4 figures

Recent advancements in Earth system science have been marked by the exponential increase in the availability of diverse, multivariate datasets characterised by moderate to high spatio-temporal resolutions. Earth System Data Cubes (ESDCs) have emerged as one suitable solution for transforming this flood of data into a simple yet robust data structure. ESDCs achieve this by organising data into an analysis-ready format aligned with a spatio-temporal grid, facilitating user-friendly analysis and diminishing the need for extensive technical data processing knowledge. Despite these significant benefits, the completion of the entire ESDC life cycle remains a challenging task. Obstacles are not only of a technical nature but also relate to domain-specific problems in Earth system research. There exist barriers to realising the full potential of data collections in light of novel cloud-based technologies, particularly in curating data tailored for specific application domains. These include transforming data to conform to a spatio-temporal grid with minimum distortions and managing complexities such as spatio-temporal autocorrelation issues. Addressing these challenges is pivotal for the effective application of Artificial Intelligence (AI) approaches. Furthermore, adhering to open science principles for data dissemination, reproducibility, visualisation, and reuse is crucial for fostering sustainable research. Overcoming these challenges offers a substantial opportunity to advance data-driven Earth system research, unlocking the full potential of an integrated, multidimensional view of Earth system processes. This is particularly true when such research is coupled with innovative research paradigms and technological progress.

We study the asymptotic growth of the number of rational points of bounded height on smooth projective split toric varieties with Picard rank 2 over number fields, with respect to Arakelov height functions associated with big metrized line bundles. We show that these varieties can be naturally decomposed into a finite disjoint union of subvarieties, where explicit asymptotic formulas for the number of rational points of bounded height can be given. Additionally, we present various examples, including the case of Hirzebruch surfaces.

Context. Mock galaxy catalogues are essential for correctly interpreting current and future generations of galaxy surveys. Despite their significance in galaxy formation and cosmology, little to no work has been done to validate the predictions of these mocks for high-order clustering statistics. Aims. We compare the predicting power of the latest generation of empirical models used in the creation of mock galaxy catalogues: a 13-parameter Halo Occupation Distribution (HOD) and an extension of the SubHalo Abundance Matching technique (SHAMe). Methods. We build GalaxyEmu-Planck, an emulator that makes precise predictions for the two-point correlation function, galaxy-galaxy lensing (restricted to distances greater than 1 $h^{-1} {\rm Mpc}$ to avoid baryonic effects), and other high-order statistics resulting from the evaluation of SHAMe and HOD models. Results. We evaluate the precision of GalaxyEmu-Planck using two galaxy samples extracted from the FLAMINGO hydrodynamical simulation that mimic the properties of DESI-BGS and BOSS galaxies, finding that the emulator reproduces all the predicted statistics precisely. The HOD showed comparable performance when fitting galaxy clustering and galaxy-galaxy lensing. In contrast, the SHAMe model showed better predictions for higher-order statistics, especially regarding the galaxy assembly bias. We also tested the performance of the models after removing some of their extensions, finding that we can withdraw two of the HOD parameters without a loss of performance. Conclusions. The results of this paper validate the current generation of empirical models as a way to reproduce galaxy clustering, galaxy-galaxy lensing and other high-order statistics. The excellent performance of the SHAMe model with a small number of free parameters suggests that it is a valid method to extract cosmological constraints from galaxy clustering., Comment: 22 pages, 12 figures, submitted to A&A. Comments welcomed!

Our study focuses on characterizing the highly ionized gas within the Milky Way's (MW) Circumgalactic Medium (CGM) that gives rise to ionic transitions in the X-ray band 2 - 25 \AA. Utilizing stacked \Chandra/\ACISS\ \MEG\ and \LETG\ spectra toward QSO sightlines, we employ the self-consistent hybrid ionization code PHASE to model our data. The stacked spectra are optimally described by three distinct gas phase components: a \warm\ (\logT\ $\sim$ 5.5), \warmhot\ (\logT\ $\sim 6$), and \hot\ (\logT\ $\sim$ 7.5) components. These findings confirm the presence of the \hot\ component in the MW's CGM indicating its coexistence with a \warm\ and a \warmhot\ gas phases. We find this \hot\ component to be homogeneous in temperature but inhomogeneous in column density. The gas in the \hot\ component requires over-abundances relative to solar to be consistent with the Dispersion Measure (DM) from the Galactic halo reported in the literature. {For the hot phase we estimated a DM = $55.1^{+29.9}_{-23.7}$ pc cm$^{-3}$}. We conclude that this phase is either enriched in Oxygen, Silicon, and Sulfur, or has metallicity {over 6} times solar value, or a combination of both. We do not detect Fe L-shell absorption lines, implying O/Fe $\geq$ 4. The non-solar abundance ratios found in the super-virial gas component in the Galactic halo suggest that this phase arises from Galactic feedback.

When constructing galaxy mock catalogs based on suits of dark matter halo catalogs generated with approximated, calibrated or machine-learning approaches, the assignment of intrinsic properties for such tracers is a step of paramount relevance, given that these can shape the abundance of mock galaxy cluster and the spatial distribution of mock galaxies. We explore the possibility to assign properties of dark matter halos within the context of calibrated/learning approaches, explicitly using clustering information. The goal is to retrieve the correct signal of primary and secondary large-effective bias as a function of properties reconstructed solely based on phase-space properties of the halo distribution and dark matter density field. The algorithm reconstructs a set halo properties (such as virial mass, maximum circular velocity, concentration and spin) constraint to reproduce both primary and secondary (or assembly) bias. The key ingredients of the algorithm are the implementation of individually-assigned large-scale effective bias, a multi-scale approach to account for halo exclusion and a hierarchical assignment of halo properties. The method facilitates the assignment of halo properties aiming at replicating the large-scale effective bias, both primary and secondary. This improves over over previous methods found in the literature, especially at the high mass population. We have designed an strategy to reconstruct the main properties of dark matter halos obtained by calibrated/learning algorithms, in a way that the expected one and two-point statistics (on large scales) replicates the signal from detailed N-body simulations. We encourage the application of this strategy (or the implementation of our algorithm) for the generation of mock catalogs of dark matter halos based on approximated methods., Comment: 10 Pages, 6 Figures. Accepted for publication at A&A

In the study of crystal nucleation via computer simulations, hard spheres are arguably the most extensively explored model system. Nonetheless, even in this simple model system, the complex thermodynamics of crystal nuclei can sometimes give rise to counterintuitive results, such as the recent observation that the pressure inside a critical nucleus is lower than that of the surrounding fluid, seemingly clashing with the strictly positive Young--Laplace pressure we would expect in liquid droplets. Here, we re-derive many of the founding equations associated with crystal nucleation, and use the hard-sphere model to demonstrate how they give rise to this negative pressure difference. We exploit the fact that, in the canonical ensemble, a nucleus can be in a (meta)stable equilibrium with the fluid, and measure the surface stress for both flat and curved interfaces. Additionally, we explain the effect of defects on the chemical potential inside the crystal nucleus. Lastly, we present a simple, fitted thermodynamic model to capture the properties of the nucleus, including the work required to form critical nuclei.

Giant (>100 kpc) nebulae associated with active galaxies provide rich information about the circumgalactic medium (CGM) around galaxies, its link with the interstellar medium (ISM) of the hosts and the mechanisms involved in their evolution. We have studied the giant nebula associated with the Teacup (z=0.085) quasar based on VLT MUSE integral field spectroscopy to investigate whether the well known giant (~10 kpc) active galactic nucleus (AGN) induced outflow has an impact on the distribution of heavy elements in and outside the host galaxy. For this, we have mapped the O/H and N/O abundances in two spatial dimensions across the giant nebula and within the galaxy by means of comparing emission line ratios with photoionisation model predictions. We have found that the widely studied AGN driven outflow responsible for the 10 kpc ionised bubble is enhancing the gas metal abundance up to 10 kpc from the AGN. O/H is solar or slightly higher in the bubble edges, in comparison with the subsolar abundances across the rest of the nebula (median O/H~0.63 O/H_sun)). The main conclusion is that AGN feedback can produce metal enrichment at large extranuclear distances in galaxies (>~10 kpc)., Comment: Accepted for publicacion in A&A. 8 pages, 4 figures

In this note, we propose a generalized notion of distance between vacua in the theory of a scalar field $\phi$ with scalar potential $V(\phi)$ coupled to gravity. We propose the normalized tension of domain wall connecting different field values, with a varying normalization relative to a local energy scale, as the distance. We show this definition reproduces the usual moduli space distance for zero potential, as well as the $d\propto |\log \Lambda|$ behavior with the vacuum energy $\Lambda$ in the AdS case, previously proposed in the literature. In the case of large AdS we also obtain the expected exponent of mass versus distance in one particular case, when the mass of the light tower is $m\sim \sqrt \Lambda$ and there is a single extra dimension decompactifying. We also discuss the features and shortcomings of alternative but related proposals., Comment: 27 pages, 3 figures, typos corrected, references added

Being the most prominent HI line, Ly$\alpha$ permeates the cosmic web in emission. Despite its potential as a cosmological probe, its detection on large scales remains elusive. We present a new methodology to perform Ly$\alpha$ intensity mapping with broad-band optical images, by cross-correlating them with Ly$\alpha$ forest data using a custom one-parameter estimator. We also develop an analytical large-scale Ly$\alpha$ emission model with two parameters (average luminosity $\langle L_{\rm Ly\alpha} \rangle$ and bias $b_{\rm e}$) that respects observational constraints from QSO luminosity functions. We compute a forecast for DECaLS/BASS $g$-band images cross-correlated with DESI Ly$\alpha$ forest data, setting guidelines for reducing images into Ly$\alpha$ intensity maps. Given the transversal scales of our cross-correlation (26.4 arcmin, $\sim$33 cMpc/h), our study effectively integrates Ly$\alpha$ emission over all the cosmic volume inside the DESI footprint at $2.2 < z < 3.4$ (the $g$-band Ly$\alpha$ redshift range). Over the parameter space ($\langle L_{\rm Ly\alpha} \rangle$, $b_{\rm e}$) sampled by our forecast, we find a 3$\sigma$ of large-scale structure in Ly$\alpha$ likely, with a probability of detection of 23.95\% for DESI-DECaLS/BASS, and 54.93\% for a hypothetical DESI phase II with twice as much Ly$\alpha$ QSOs. Without a detection, we derive upper bounds on $\langle L_{\rm Ly\alpha} \rangle$ competitive with optimistic literature estimates ($2.3 \pm 1 \cdot 10^{\rm 41}$ erg/s/cMpc$^3$ for DESI, and $\sim$35\% lower for its hypothetical phase II). Extrapolation to the DESI-Rubin overlap shows that a detection of large-scale structure with Ly$\alpha$ intensity mapping using next-generation imaging surveys is certain. [abridged], Comment: 26 pages, 24 figures, submitted to MNRAS. Comments welcome!

With climate extremes' rising frequency and intensity, robust analytical tools are crucial to predict their impacts on terrestrial ecosystems. Machine learning techniques show promise but require well-structured, high-quality, and curated analysis-ready datasets. Earth observation datasets comprehensively monitor ecosystem dynamics and responses to climatic extremes, yet the data complexity can challenge the effectiveness of machine learning models. Despite recent progress in deep learning to ecosystem monitoring, there is a need for datasets specifically designed to analyse compound heatwave and drought extreme impact. Here, we introduce the DeepExtremeCubes database, tailored to map around these extremes, focusing on persistent natural vegetation. It comprises over 40,000 spatially sampled small data cubes (i.e. minicubes) globally, with a spatial coverage of 2.5 by 2.5 km. Each minicube includes (i) Sentinel-2 L2A images, (ii) ERA5-Land variables and generated extreme event cube covering 2016 to 2022, and (iii) ancillary land cover and topography maps. The paper aims to (1) streamline data accessibility, structuring, pre-processing, and enhance scientific reproducibility, and (2) facilitate biosphere dynamics forecasting in response to compound extremes.

Agent-based models describing social interactions among individuals can help to better understand emerging macroscopic patterns in societies. One of the topics which is worth tackling is the formation of different kinds of hierarchies that emerge in social spaces such as cities. Here we propose a Bonabeau-like model by adding a second class of agents. The fundamental particularity of our model is that only a pairwise interaction between agents of the opposite class is allowed. Agent fitness can thus only change by competition among the two classes, while the total fitness in the society remains constant. The main result is that for a broad range of values of the model parameters, the fitness of the agents of each class show a decay in time except for one or very few agents which capture almost all the fitness in the society. Numerical simulations also reveal a singular shift from egalitarian to hierarchical society for each class. This behaviour depends on the control parameter $\eta$, playing the role of the inverse of the temperature of the system. Results are invariant with regard to the system size, contingent solely on the quantity of agents within each class. Finally, a couple of scaling laws are provided thus showing a data collapse from different model parameters and they follow a shape which can be related to the presence of a phase transition in the model., Comment: 19 pages with 15 figures

We investigate low-mass central galaxies with Mstar = $10^{9.5}-10^{10}$ Msun/h, located near massive groups and galaxy clusters using the TNG300 and MDPL2-SAG simulations. We set out to study their evolution, aiming to find hints about the large-scale conformity signal they produce. We also use a control sample of low-mass central galaxies located far away from massive structures. For both samples, we find a sub-population of galaxies that were accreted by another halo in the past but are now considered central galaxies; we refer to these objects as former satellites. The fraction of former satellites is higher for quenched central galaxies near massive systems: 45% in TNG300 and 17% in MDPL2-SAG. Our results in TNG300 show that former satellites were typically hosted by massive dark matter halos (M200 $\geq 10^{13}$ Msun/h) at z$\sim$0.3, followed by a drop in halo mass at lower redshifts. In addition, we find a strong drop in the total gas mass at z$\leq$1 for quenched central galaxies near galaxy groups and clusters produced by these former satellites as well. By removing former satellites, the evolution of quenched central galaxies is fairly similar to those of the quenched control galaxies, showing small differences at low-z. For MDPL2-SAG, former satellites were hosted by less massive halos, with a mean halo mass around $10^{11}$ Msun/h, and the evolution remains equal before and after removing former satellites. We also measure the two-halo conformity, i.e., the correlation in the specific SFR between low-mass central galaxies and their neighbors at Mpc scales, and how former satellites contribute to the signal at z=0, 0.3, and 1. The conformity signal decreases from z=0 to z=1 in MDPL2-SAG but it increases in TNG300. However, after removing former satellites in TNG300, the signal is strongly reduced but almost does not change at z$\leq$0.3, and it disappears at z=1 (abridged)., Comment: 16 pages, 13 figures, submitted to A&A

Biological neural networks are characterized by their high degree of plasticity, a core property that enables the remarkable adaptability of natural organisms. Importantly, this ability affects both the synaptic strength and the topology of the nervous systems. Artificial neural networks, on the other hand, have been mainly designed as static, fully connected structures that can be notoriously brittle in the face of changing environments and novel inputs. Building on previous works on Neural Developmental Programs (NDPs), we propose a class of self-organizing neural networks capable of synaptic and structural plasticity in an activity and reward-dependent manner which we call Lifelong Neural Developmental Program (LNDP). We present an instance of such a network built on the graph transformer architecture and propose a mechanism for pre-experience plasticity based on the spontaneous activity of sensory neurons. Our results demonstrate the ability of the model to learn from experiences in different control tasks starting from randomly connected or empty networks. We further show that structural plasticity is advantageous in environments necessitating fast adaptation or with non-stationary rewards., Comment: 10 pages

Single-field models of accelerated expansion with nearly flat potentials, despite being able to provide observationally viable explanations for the early-time cosmic inflation and the late-time cosmic acceleration, are in strong tension with string theory evidence and the associated de Sitter swampland constraints. It has recently been argued that in an open universe, where the spatial curvature is negative (i.e., with $\Omega_k>0$), a new stable fixed point arises, which may lead to viable single-field-based accelerated expansion with an arbitrarily steep potential. Here we show, through a dynamical systems analysis and a Bayesian statistical inference of cosmological parameters, that the additional cosmological solutions based on the new fixed point do not render steep-potential, single-field, accelerated expansion observationally viable. We mainly focus on quintessence models of dark energy, but we also argue that a similar conclusion can be drawn for cosmic inflation., Comment: 13 pages, 6 figures. v2: expanded discussion and added references

Representations for black-box optimisation methods (such as evolutionary algorithms) are traditionally constructed using a delicate manual process. This is in contrast to the representation that maps DNAs to phenotypes in biological organisms, which is at the hear of biological complexity and evolvability. Additionally, the core of this process is fundamentally the same across nearly all forms of life, reflecting their shared evolutionary origin. Generative models have shown promise in being learnable representations for black-box optimisation but they are not per se designed to be easily searchable. Here we present a system that can meta-learn such representation by directly optimising for a representation's ability to generate quality-diversity. In more detail, we show our meta-learning approach can find one Neural Cellular Automata, in which cells can attend to different parts of a "DNA" string genome during development, enabling it to grow different solvable 2D maze structures. We show that the evolved genotype-to-phenotype mappings become more and more evolvable, not only resulting in a faster search but also increasing the quality and diversity of grown artefacts., Comment: 10 pages, 6 Figures, Proceedings of the 2024 Artificial Life Conference

We study a scenario with a transient phase of cosmological acceleration that could potentially be realized in asymptotic corners of String Theory moduli space. A very steep scalar potential is temporarily stabilized by the effect of a nonzero density of heavy states, leading to acceleration, in what amounts to a cosmological version of the Chameleon mechanism. The density of heavy states is diluted by cosmological expansion, weakening their effect. After roughly one $e$-fold their effect can no longer stabilize the potential, and the accelerating phase ends. We also study a scenario where there is no potential and the transient acceleration is achieved by the counterbalancing effects of light and heavy towers of states. In both cases we show that it is not possible to obtain more than $\mathcal{O}(1)$ $e$-folds without transplanckian field excursions. We also discuss the general EFT constraints on these models and explore a number of first attempts at concrete embeddings of the scenario in String Theory. These all turn out to face significant challenges., Comment: 27 pages, 6 figures

The Dark Energy Spectroscopic Instrument (DESI) will measure millions of quasar spectra by the end of its 5 year survey. Quasar redshift errors impact the shape of the Lyman-$\alpha$ forest correlation functions, which can affect cosmological analyses and therefore cosmological interpretations. Using data from the DESI Early Data Release and the first two months of the main survey, we measure the systematic redshift error from an offset in the cross-correlation of the Lyman-$\alpha$ forest with quasars. We find evidence for a redshift dependent bias causing redshifts to be underestimated with increasing redshift, stemming from improper modeling of the Lyman-$\alpha$ optical depth in the templates used for redshift estimation. New templates were derived for the DESI Year 1 quasar sample at $z > 1.6$ and we found the redshift dependent bias, $\Delta r_\parallel$, increased from $-1.94 \pm 0.15$ $h^{-1}$ Mpc to $-0.08 \pm 0.04$ $h^{-1}$ Mpc ($-205 \pm 15~\text{km s}^{-1}$ to $-9.0 \pm 4.0~\text{km s}^{-1}$). These new templates will be used to provide redshifts for the DESI Year 1 quasar sample., Comment: 29 pages, 9 figures, 5 tables

The standard cosmological model, with its six independent parameters, successfully describes our observable Universe. One of these parameters, the optical depth to reionization $\tau_\mathrm{reio}$, represents the scatterings that Cosmic Microwave Background (CMB) photons will experience after decoupling from the primordial plasma as the intergalactic medium transitions from neutral to ionized. $\tau_\mathrm{reio}$ depends on the neutral hydrogen fraction $x_\mathrm{HI}$, which, in turn, should theoretically depend on cosmology. We present a novel method to establish the missing link between cosmology and reionization timeline using symbolic regression. We discover the timeline has a universal shape well described by the Gompertz mortality law, applicable to any cosmology within our simulated data. Unlike the conventional tanh prescription, our model is asymmetric in time and a good fit to astrophysical constraints on $x_\mathrm{HI}$. By combining CMB with astrophysical data and marginalizing over astrophysics, we treat $\tau_\mathrm{reio}$ as a derived parameter, tightening its constraint to $<3\%$. This approach reduces the error on the amplitude of the primordial fluctuations by a factor of 2.3 compared to Planck's PR3 constraint and provides a commanding constraint on the ionization efficiency $\zeta_\mathrm{UV} = 26.9^{+2.1}_{-2.5}$. We expect further improvements in the near term as reionization constraints increase and our understanding of reionization advances., Comment: Correctly accounting for ionization efficiency enables us to place leading constraints by combining our universal asymmetric reionization model with CMB and damping wing constraints on $x_\mathrm{HI}$. By marginalizing over $\zeta_\mathrm{UV}$, we achieve powerful constraints on $A_\mathrm{s}$ and $\tau_\mathrm{reio}$, with uncertainties around 3% for the former

The motion of a graph is the minimum number of vertices that are moved by a non-trivial automorphism. Equivalently, it can be defined as the minimal degree of its automorphism group (as a permutation group on the vertices). In this paper we develop some results on permutation groups (primitive and imprimitive) with small minimal degree. As a consequence of such results we classify vertex-transitive graphs whose motion is $4$ or a prime number.

Broad absorption line quasars (BALs) exhibit blueshifted absorption relative to a number of their prominent broad emission features. These absorption features can contribute to quasar redshift errors and add absorption to the Lyman-alpha (LyA) forest that is unrelated to large-scale structure. We present a detailed analysis of the impact of BALs on the Baryon Acoustic Oscillation (BAO) results with the LyA forest from the first year of data from the Dark Energy Spectroscopic Instrument (DESI). The baseline strategy for the first year analysis is to mask all pixels associated with all BAL absorption features that fall within the wavelength region used to measure the forest. We explore a range of alternate masking strategies and demonstrate that these changes have minimal impact on the BAO measurements with both DESI data and synthetic data. This includes when we mask the BAL features associated with emission lines outside of the forest region to minimize their contribution to redshift errors. We identify differences in the properties of BALs in the synthetic datasets relative to the observational data, as well as use the synthetic observations to characterize the completeness of the BAL identification algorithm, and demonstrate that incompleteness and differences in the BALs between real and synthetic data also do not impact the BAO results for the LyA forest., Comment: 32 pages, 8 figures, accepted by JCAP

Abstract polytopes are combinatorial objects that generalise geometric objects such as convex polytopes, maps on surfaces and tilings of the space. Chiral polytopes are those abstract polytopes that admit full combinatorial rotational symmetry but do not admit reflections. In this paper we build chiral polytopes whose facets (maximal faces) are isomorphic to a prescribed regular cubic tessellation of the $n$-dimensional torus ($n \geq 2$). As a consequence, we prove that for every $d \geq 3$ there exist infinitely many chiral $d$-polytopes.