Search

Cluster cosmology can benefit from combining multi-wavelength studies, which can benefit from characterising the correlation coefficients between different mass-observable relations. In this work, we aim to provide information on the scatter, the skewness, and the covariance of various mass-observable relations in galaxy clusters in cosmological hydrodynamic simulations. This information will help future analyses to better tackle accretion histories and projection effects and model mass observable relations for cosmology studies.We identify galaxy clusters in Magneticum Box2b simulations with mass $M_{\rm 200c}>10^{14} {\rm M}_\odot$ at redshift $z=0.24$ and $z=0.90$. Our analysis includes \Euclid-derived properties such as richness, stellar mass, lensing mass, and concentration. Additionally, we investigate complementary multi-wavelength data, including X-ray luminosity, integrated Compton-$y$ parameter, gas mass, and temperature. The impact of projection effects on mass-observable residuals and correlations is then examined. At intermediate redshift ($z=0.24$), projection effects impact lensing concentration, richness, and gas mass the most in terms of scatter and skewness of log-residuals of scaling relations. The contribution of projection effects can be significant enough to boost a spurious hot- vs. cold-baryons correlation and consequently hide underlying correlations due to halo accretion histories. At high redshift ($z=0.9$), the richness has a much lower scatter (of log-residuals), and the quantity that is most impacted by projection effects is the lensing mass. Lensing concentration reconstruction, in particular, is affected by deviations of the reduced-shear profile shape from the one derived by an NFW profile rather than interlopers in the line of sight., Comment: 26 pages, 22 figures, submitted to A&A, replied to referee comments

Beyond testing the current cosmological paradigm, cluster number counts can also be utilized to investigate the discrepancies currently affecting current cosmological measurements. In particular, cosmological studies based on cosmic shear and other large-scale structure probes routinely find a value of the amplitude of the fluctuations in the universe S8 smaller than the one inferred from the primary cosmic microwave background. In this work, we investigate this tension by measuring structure evolution across cosmic time as probed by the number counts of the massive halos with the first SRG/eROSITA All-Sky Survey cluster catalog in the Western Galactic Hemisphere complemented with the overlapping Dark Energy Survey Year-3, KiloDegree Survey, and Hyper Suprime-Cam data for weak lensing mass calibration, by implementing two different parameterizations and a model-agnostic method. In the first model, we measure the cosmic linear growth index as {\gamma} = 1.19 \pm 0.21, in tension with the standard value of {\gamma} = 0.55, but in good statistical agreement with other large-scale structures probes. The second model is a phenomenological scenario in which we rescale the linear matter power spectrum at low redshift to investigate a potential reduction of structure formation, providing similar results. Finally, in a third strategy, we consider a standard {\Lambda}CDM cosmology, but we separate the cluster catalog into five redshift bins, measuring the cosmological parameters in each and inferring the evolution of the structure formation, finding hints of a reduction. Interestingly, the S8 value inferred from eRASS1 cluster number counts, when we add a degree of freedom to the matter power spectrum, recovers the value inferred by cosmic shear studies.

In standard cosmology, the late Universe is assumed to be statistically homogeneous and isotropic. However, a recent study based on galaxy clusters by Migkas et al. (2021, arXiv:2103.13904) found an apparent spatial variation of approximately $9\%$ in the Hubble constant, $H_0$, across the sky. The authors utilised galaxy cluster scaling relations between various cosmology-dependent cluster properties and a cosmology-independent property, i.e., the temperature of the intracluster gas $(T)$. A position-dependent systematic bias of $T$ measurements can, in principle, result in an overestimation of apparent $H_0$ variations. In this study, we search for directional $T$ measurement biases by examining the scaling relation between the member galaxy velocity dispersion and the gas temperature $(\sigma_\mathrm{v}-T)$. Additionally, we search for apparent $H_0$ angular variations independently of $T$ by analysing the relations between the X-ray luminosity and Sunyaev-Zeldovich signal with the velocity dispersion, $L_\mathrm{X}-\sigma_\mathrm{v}$ and $Y_\mathrm{SZ}-\sigma_\mathrm{v}$. We utilise Monte Carlo simulations of isotropic cluster samples to quantify the statistical significance of any observed anisotropies. We find no significant directional $T$ measurement biases, and the probability that a directional $T$ bias causes the previously observed $H_0$ anisotropy is only $0.002\%$. On the other hand, from the joint analysis of the $L_\mathrm{X}-\sigma_\mathrm{v}$ and $Y_\mathrm{SZ}-\sigma_\mathrm{v}$ relations, the maximum variation of $H_0$ is found in the direction of $(295^\circ\pm71^\circ, -30^\circ\pm71^\circ)$ with a statistical significance of $3.64\sigma$, fully consistent with arXiv:2103.13904. Our findings strongly corroborate the previously detected spatial anisotropy of galaxy cluster scaling relations using a new independent cluster property, $\sigma_\mathrm{v}$., Comment: Submitted to the Astronomy & Astrophysics journal: 18 pages, 19 figures (main text), 7 figures (appendix)

Interactions between derivatives and fixpoints have many important applications in both computer science and mathematics. In this paper, we provide a categorical framework to combine fixpoints with derivatives by studying Cartesian differential categories with a fixpoint operator. We introduce an additional axiom relating the derivative of a fixpoint with the fixpoint of the derivative. We show how the standard examples of Cartesian differential categories where we can compute fixpoints provide canonical models of this notion. We also consider when the fixpoint operator is a Conway operator, or when the underlying category is closed. As an application, we show how this framework is a suitable setting to formalize the Newton-Raphson optimization for fast approximation of fixpoints and extend it to higher order languages., Comment: 16 pages, LICS 24 conference

For decades, Simultaneous Ascending Auction (SAA) has been the most widely used mechanism for spectrum auctions, and it has recently gained popularity for allocating 5G licenses in many countries. Despite its relatively simple rules, SAA introduces a complex strategic game with an unknown optimal bidding strategy. Given the high stakes involved, with billions of euros sometimes on the line, developing an efficient bidding strategy is of utmost importance. In this work, we extend our previous method, a Simultaneous Move Monte-Carlo Tree Search (SM-MCTS) based algorithm named $SMS^{\alpha}$ to incomplete information framework. For this purpose, we compare three determinization approaches which allow us to rely on complete information SM-MCTS. This algorithm addresses, in incomplete framework, the four key strategic issues of SAA: the exposure problem, the own price effect, budget constraints, and the eligibility management problem. Through extensive numerical experiments on instances of realistic size with an uncertain framework, we show that $SMS^{\alpha}$ largely outperforms state-of-the-art algorithms by achieving higher expected utility while taking less risks, no matter which determinization method is chosen., Comment: 15 pages, 9 figures, 2 pseudocodes

There is an abstract notion of connection in any tangent category. In this paper, we show that when applied to the tangent category of affine schemes, this recreates the classical notion of a connection on a module (and similarly, in the tangent category of schemes, this recreates the notion of connection on a quasi-coherent sheaf of modules). By contrast, we also show that in the tangent category of algebras, there are no non-trivial connections., Comment: Minor changes and typos fixes based on reviewer comments

We present radio observations of the galaxy cluster Abell S1136 at 888 MHz, using the Australian Square Kilometre Array Pathfinder radio telescope, as part of the Evolutionary Map of the Universe Early Science program. We compare these findings with data from the Murchison Widefield Array, XMM-Newton, the Wide-field Infrared Survey Explorer, the Digitised Sky Survey, and the Australia Telescope Compact Array. Our analysis shows the X-ray and radio emission in Abell S1136 are closely aligned and centered on the BCG, while the X-ray temperature profile shows a relaxed cluster with no evidence of a cool core. We find that the diffuse radio emission in the centre of the cluster shows more structure than seen in previous low-resolution observations of this source, which appeared formerly as an amorphous radio blob, similar in appearance to a radio halo; our observations show the diffuse emission in the Abell S1136 galaxy cluster contains three narrow filamentary structures visible at 888 MHz, between$\sim 80$ and 140 kpc in length; however the properties of the diffuse emission do not fully match that of a radio (mini-)halo or (fossil) tailed radio source., Comment: To appear in PASA

Cartesian differential categories provide a categorical framework for multivariable differential calculus and also the categorical semantics of the differential $\lambda$-calculus. Taylor series expansion is an important concept for both differential calculus and the differential $\lambda$-calculus. In differential calculus, a function is equal to its Taylor series if its sequence of Taylor polynomials converges to the function in the analytic sense. On the other hand, for the differential $\lambda$-calculus, one works in a setting with an appropriate notion of algebraic infinite sums to formalize Taylor series expansion. In this paper, we provide a formal theory of Taylor series in an arbitrary Cartesian differential category without the need for converging limits or infinite sums. We begin by developing the notion of Taylor polynomials of maps in a Cartesian differential category and then show how comparing Taylor polynomials of maps induces an ultrapseudometric on the homsets. We say that a Cartesian differential category is Taylor if maps are entirely determined by their Taylor polynomials. The main results of this paper are that in a Taylor Cartesian differential category, the induced ultrapseudometrics are ultrametrics and that for every map $f$, its Taylor series converges to $f$ with respect to this ultrametric. This framework recaptures both Taylor series expansion in differential calculus via analytic methods and in categorical models of the differential $\lambda$-calculus (or Differential Linear Logic) via infinite sums., Comment: To be published in the proceedings of MFPS2024

This paper explores two condensed-space interior-point methods to efficiently solve large-scale nonlinear programs on graphics processing units (GPUs). The interior-point method solves a sequence of symmetric indefinite linear systems, or Karush-Kuhn-Tucker (KKT) systems, which become increasingly ill-conditioned as we approach the solution. Solving a KKT system with traditional sparse factorization methods involve numerical pivoting, making parallelization difficult. A solution is to condense the KKT system into a symmetric positive-definite matrix and solve it with a Cholesky factorization, stable without pivoting. Although condensed KKT systems are more prone to ill-conditioning than the original ones, they exhibit structured ill-conditioning that mitigates the loss of accuracy. This paper compares the benefits of two recent condensed-space interior-point methods, HyKKT and LiftedKKT. We implement the two methods on GPUs using MadNLP.jl, an optimization solver interfaced with the NVIDIA sparse linear solver cuDSS and with the GPU-accelerated modeler ExaModels.jl. Our experiments on the PGLIB and the COPS benchmarks reveal that GPUs can attain up to a tenfold speed increase compared to CPUs when solving large-scale instances.

The current standard model of cosmology successfully describes a variety of measurements, but the nature of its main ingredients, dark matter and dark energy, remains unknown. Euclid is a medium-class mission in the Cosmic Vision 2015-2025 programme of the European Space Agency (ESA) that will provide high-resolution optical imaging, as well as near-infrared imaging and spectroscopy, over about 14,000 deg^2 of extragalactic sky. In addition to accurate weak lensing and clustering measurements that probe structure formation over half of the age of the Universe, its primary probes for cosmology, these exquisite data will enable a wide range of science. This paper provides a high-level overview of the mission, summarising the survey characteristics, the various data-processing steps, and data products. We also highlight the main science objectives and expected performance., Comment: Accepted for publication in the A&A special issue`Euclid on Sky'

Methods. We utilized the combined five SRG/eROSITA All-Sky Survey data (eRASS:5) to perform X-ray imaging and spectral analyses of the Centaurus cluster in various directions to large radii. Surface brightness (SB) profiles out to $2R_{200}$ were constructed. We acquired gas temperature, metallicity, and normalization per area profiles out to $R_{200}$. We compared our results with previous Centaurus studies, cluster outskirts measurements, and simulations. Comprehensive sky background analysis was done across the FoV, in particular, to assess the variation of the eROSITA Bubble emission that partially contaminates the field. Results. The processed X-ray images show the known sloshing-induced structures in the core. The core ($r\leq11~\mathrm{kpc}$) is better described with a 2T model than a 1T model. Here, we measured lower T from the cooler component (~1.0 keV) and higher Z ($\sim\!1.6Z_\odot$), signifying an iron bias. In the intermediate radial range, we observed prominent SB and normalization per area excesses in the eastern sector (Cen 45 location), reaching out to $R_{500}$. Temperature enhancements near the location of Cen 45 imply that the gas is shock-heated due to the interaction with Cen 30, the significant excess behind Cen 45 center might be the tail/ram-pressure-stripped gas. We found good agreement between the outskirt temperatures with the profile from simulations and fit from Suzaku outskirts measurements. We detected significant SB emission to the sky background level out to $R_{200}$ with a $3.5\sigma$ and followed by $2.9\sigma$ at $1.1R_{200}$. The metallicity at $R_{500}-R_{200}$ is low but within the ranges of other outskirts studies. Conclusions. We present the first measurement of ICM morphology and properties of Centaurus cluster sampling the whole azimuth beyond $30'$, increasing the probed volume by a factor of almost 30., Comment: Submitted to the Astronomy & Astrophysics journal: 19 pages, 10 figures (main text), 2 figure (appendix)

We investigate the potential of Graphics Processing Units (GPUs) to solve large-scale nonlinear programs with a dynamic structure. Using ExaModels, a GPU-accelerated automatic differentiation tool, and the interior-point solver MadNLP, we significantly reduce the time to solve dynamic nonlinear optimization problems. The sparse linear systems formulated in the interior-point method is solved on the GPU using a hybrid solver combining an iterative method with a sparse Cholesky factorization, which harness the newly released NVIDIA cuDSS solver. Our results on the classical distillation column instance show that despite a significant pre-processing time, the hybrid solver allows to reduce the time per iteration by a factor of 25 for the largest instance., Comment: 6 pages, 1 figure, 2 tables

Drazin inverses are a fundamental algebraic structure which have been extensively deployed in semigroup theory and ring theory. Drazin inverses can also be defined for endomorphisms in any category. However, beyond a paper by Puystjens and Robinson from 1987, not much has been done with Drazin inverses in category theory. As such, here we provide a survey of the theory of Drazin inverses from a categorical perspective. We introduce Drazin categories, in which every endomorphism has a Drazin inverse, and provide various examples including the category of matrices over a field, the category of finite length modules over a ring, and finite set enriched categories. We also introduce the notion of expressive rank and prove that a category with expressive rank is Drazin. Moreover, we study Drazin inverses in mere categories, in additive categories, and in dagger categories. In an arbitrary category, we show how a Drazin inverse corresponds to an isomorphism in the idempotent splitting, as well as explain how Drazin inverses relate to Leinster's notion of eventual image duality. In additive categories, we explore the consequences of the core-nilpotent decomposition and the image-kernel decomposition, which we relate back to Fitting's famous results. We then develop the notion of Drazin inverses for pairs of opposing maps, generalizing the usual notion of Drazin inverse for endomorphisms. As an application of this new kind of Drazin inverse, for dagger categories, we provide a novel characterization of the Moore-Penrose inverse in terms of being a Drazin inverse of the pair of a map and its adjoint.

The evolution of the cluster mass function traces the growth of the linear density perturbations and can be utilized for constraining the parameters of cosmological and alternative gravity models. In this context, we present new constraints on potential deviations from general relativity by investigating the Hu-Sawicki parametrization of the f(R) gravity with the first SRG-eROSITA All-Sky Survey (eRASS1) cluster catalog in the Western Galactic Hemisphere in combination with the overlapping Dark Energy Survey Year 3, KiloDegree Survey and Hyper Supreme Camera data for weak lensing mass calibration. For the first time, we present constraints obtained from cluster abundances only. When we consider massless neutrinos, we find a strict upper limit of log |fR0| < -4.31 at 95% confidence level. Massive neutrinos suppress structure growth at small scales, and thus have the opposite effect of f(R) gravity. We consequently investigate the joint fit of the mass of the neutrinos with the modified gravity parameter. We obtain log |fR0| < -4.12 jointly with \sum m_\nu < 0.44 e.V. at 95% confidence level, tighter than the limits in the literature utilizing cluster counts only. At log |fR0|= - 6, the number of clusters is not significantly changed by the theory. Consequently, we do not find any statistical deviation from general relativity from the study of eRASS1 cluster abundance. Deeper surveys with eROSITA, increasing the number of detected clusters, will further improve constraints on log |fR0| and investigate alternative gravity theories., Comment: 13 pages, 6 figures, submitted to A&A

Number counts of galaxy clusters across redshift are a powerful cosmological probe, if a precise and accurate reconstruction of the underlying mass distribution is performed -- a challenge called mass calibration. With the advent of wide and deep photometric surveys, weak gravitational lensing by clusters has become the method of choice to perform this measurement. We measure and validate the weak gravitational lensing (WL) signature in the shape of galaxies observed in the first 3 years of the DES Y3 caused by galaxy clusters selected in the first all-sky survey performed by SRG/eROSITA. These data are then used to determine the scaling between X-ray photon count rate of the clusters and their halo mass and redshift. We empirically determine the degree of cluster member contamination in our background source sample. The individual cluster shear profiles are then analysed with a Bayesian population model that self-consistently accounts for the lens sample selection and contamination, and includes marginalization over a host of instrumental and astrophysical systematics. To quantify the accuracy of the mass extraction of that model, we perform mass measurements on mock cluster catalogs with realistic synthetic shear profiles. This allows us to establish that hydro-dynamical modelling uncertainties at low lens redshifts ($z<0.6$) are the dominant systematic limitation. At high lens redshift the uncertainties of the sources' photometric redshift calibration dominate. With regard to the X-ray count rate to halo mass relation, we constrain all its parameters. This work sets the stage for a joint analysis with the number counts of eRASS1 clusters to constrain a host of cosmological parameters. We demonstrate that WL mass calibration of galaxy clusters can be performed successfully with source galaxies whose calibration was performed primarily for cosmic shear experiments., Comment: 27 pages, 18 figures, 2 appendices, submitted to A\&A

Characterising galaxy cluster populations from catalog of sources selected in astronomical surveys requires knowledge of sample incompleteness, known as selection function. The first All-Sky Survey (eRASS1) by eROSITA onboard Spectrum Roentgen Gamma (SRG) has enabled the collection of large samples of galaxy clusters detected in the soft X-ray band over the Western Galactic hemisphere. The driving goal consists in constraining cosmological parameters, which puts stringent requirements on accuracy, flexibility and explainability of the selection function models. We use a large set of mock observations of the eRASS1 survey and we process simulated data identically to the real eRASS1 events. We match detected sources to simulated clusters and we associate detections to intrinsic cluster properties. We train a series of models to build selection functions depending only on observable surface brightness data. We develop a second series of models relying on global cluster characteristics such as X-ray luminosity, flux, and expected instrumental count-rate as well as on morphological properties. We validate our models using our simulations and we rank them according to selected performance metrics. We validate the models with datasets of clusters detected in X-rays and via the Sunyaev-Zeldovich effect. We present the complete Bayesian population modelling framework developed for this purpose. Our results reveal the surface brightness characteristics most relevant to cluster selection in the eRASS1 sample, in particular the ambiguous role of central surface brightness at the scale of the instrument resolution. We have produced a series of user-friendly selection function models and demonstrated their validity and their limitations. Our selection function for bright sources reproduces well the catalog matches with external datasets. (abridged), Comment: 24 pages, 31 figures, 2 tables. Submitted to A&A

The cluster mass function traces the growth of linear density perturbations and provides valuable insights into the growth of structures, the nature of dark matter, and the cosmological parameters governing the Universe. The primary science goal of eROSITA, on board the {\it Spectrum Roentgen Gamma (SRG)} mission, launched in 2019, is to constrain cosmology through the evolution of cluster mass function. In this paper, we present the cosmological constraints obtained from 5259 clusters of galaxies detected over an area of 12791~deg$^2$ in the Western Galactic Hemisphere of the eROSITA's first All-Sky Survey (eRASS1). The common footprint region between the eROSITA Survey and DES, KiDS, and HSC surveys is used for calibration of the scaling between X-ray count rate and their total mass through measurements of their weak gravitational lensing signal. eRASS1 cluster abundances constrain the $\Lambda$CDM parameters, which are the energy density of the total matter to $\Omega_{\mathrm{m}}=0.29^{+0.01}_{-0.02}$, and the normalization of the density fluctuations to $\sigma_8=0.88\pm0.02$ and their combination yields $S_8=\sigma_8 (\Omega_\mathrm{m} / 0.3)^{0.5}=0.86\pm0.01$, consistent and at a similar precision with the state-of-the-art CMB measurements. eRASS1 cosmological experiment places a most stringent upper limit on the summed masses of left-handed light neutrinos to $\sum m_\nu< 0.22\mathrm{~eV}$ (95\% confidence interval). Combining eRASS1 cluster abundance measurements with CMB and ground-based neutrino oscillation experiments, we measure the summed neutrino masses to be $\sum m_\nu=0.08_{-0.02}^{+0.03}\mathrm{~eV}$ or $\sum m_\nu=0.12_{-0.01}^{+0.03}\mathrm{~eV}$ depending on the mass hierarchy scenario for neutrino eigenstates. eRASS1 cluster abundances significantly improve the constraints on the dark energy equation of state parameter to $w=-1.12\pm0.12$. (ABRIDGED), Comment: 41 pages, 22 figures, Accepted for publication by A&A

The first SRG/eROSITA All-Sky Survey (eRASS1) provides the largest intracluster medium-selected galaxy cluster and group catalog covering the western galactic hemisphere. Compared to samples selected purely on X-ray extent, the sample purity can be enhanced by identifying cluster candidates using optical and near-infrared data from the DESI Legacy Imaging Surveys. Using the red-sequence-based cluster finder eROMaPPer, we measured individual photometric properties (redshift $z_\lambda$, richness $\lambda$, optical center, and BCG position) for 12,000 eRASS1 clusters over a sky area of 13,116 deg$^2$, augmented by 247 cases identified by matching the candidates with known clusters from the literature. The median redshift of the identified eRASS1 sample is $z=0.31$, with 10% of the clusters at $z>0.72$. The photometric redshifts have an accuracy of $\delta z/(1+z)<0.005$ for $0.0595% for $z>0.05$. For these and further quality assessments of the eRASS1 identified catalog, we applied our identification method to a collection of galaxy cluster catalogs in the literature, as well as blindly on the full Legacy Surveys covering 24,069 deg$^2$. Using a combination of these cluster samples, we investigated the velocity dispersion-richness relation, finding $\log(\lambda)=2.401\times\log(\sigma)-5.074$ with an intrinsic scatter of $0.10\pm0.01$ dex. Our main result is the identified eRASS1 cluster catalog with a high purity and a well-defined X-ray selection process, enabling precise cosmological analyses presented in companion papers., Comment: 36 pages, 23 figures, 6 tables. Accepted for publication in A&A

Clusters of galaxies can be used as powerful probes to study astrophysical processes on large scales, test theories of the growth of structure, and constrain cosmological models. The driving science goal of the SRG/eROSITA All-Sky Survey (eRASS) is to assemble a large sample of X-ray-selected clusters with a well-defined selection function to determine the evolution of the mass function and, hence, the cosmological parameters. We present here a catalog of 12247 optically confirmed galaxy groups and clusters detected in the 0.2-2.3 keV as extended X-ray sources in a 13,116deg$^2$ region in the western Galactic hemisphere of the sky, which eROSITA surveyed in its first six months of operation. The clusters in the sample span the redshift range $0.003

Galaxy cluster gas temperatures ($T$) play a crucial role in many cosmological and astrophysical studies. However, it has been shown that $T$ measurements can vary between different X-ray telescopes. These $T$ biases can propagate to several cluster applications for which $T$ can be used. Thus, it is important to accurately cross-calibrate X-ray instruments to account for systematic biases. In this work, we present the cross-calibration between SRG/eROSITA and Chandra/ACIS, and between SRG/eROSITA and XMM-Newton/EPIC, using for the first time a large sample of galaxy cluster $T$. To do so, we use the first eROSITA All-Sky Survey data and a large X-ray flux-limited cluster catalog. We measure X-ray $T$ for 186 independent cluster regions with both SRG/eROSITA and Chandra/ACIS in a self-consistent way, for three energy bands; 0.7-7 keV (full), 0.5-4 keV (soft), and 1.5-7 keV (hard). We do the same with SRG/eROSITA and XMM-Newton/EPIC for 71 different cluster regions and all three bands. We find that SRG/eROSITA measures systematically lower $T$ than the other two instruments. For the full band, SRG/eROSITA returns 20$\%$ and 14$\%$ lower $T$ than Chandra/ACIS and XMM-Newton/EPIC respectively, when the two latter instruments measure $k_{\text{B}}T\approx 3$ keV each. The discrepancy increases to 38\% and 32\% when Chandra/ACIS and XMM-Newton/EPIC measure $k_{\text{B}}T\approx 10$ keV respectively. For low-$T$ galaxy groups, the discrepancy becomes milder. The soft band shows a marginally lower discrepancy than the full band. In the hard band, the cross-calibration of SRG/eROSITA and the other instruments show stronger differences. We could not identify any possible systematic biases that significantly alleviated the tension. Finally, we provide conversion factors between SRG/eROSITA, Chandra/ACIS, and XMM-Newton/EPIC $T$ which will be beneficial for future cluster studies., Comment: 24 pages, 15 figures, accepted in A&A (updated published version, replaced the submitted version v1)

As the closest galaxy cluster, the Virgo Cluster is an exemplary environment for the study of the large-scale filamentary structure and physical effects that are present in cluster outskirts but absent from the more easily studied inner regions. Here, we present an analysis of the SRG/eROSITA data from five all-sky surveys. eROSITA allows us to resolve the entire Virgo cluster and its outskirts on scales between 1 kpc and 3 Mpc, covering a total area on the sky of about 25$^\circ$ by 25$^\circ$. We utilized image manipulation techniques and surface brightness profiles to search for extended emission, surface brightness edges, and features in the outskirts. We employed a method of comparing mean and median profiles to measure gas clumping out to and beyond the virial radius. Surface brightness analysis of the cluster and individual sectors of the cluster reveal the full extent of previously identified cold fronts to the north and south. The emissivity bias due to gas clumping, which we quantify over three orders of magnitude in the radial range, is found to be mild, consistent with previous findings. We find uniform clumping measurements in all directions, with no enhancements along candidate filaments. We find an estimated virial gas mass of $M_{\mathrm{gas},r

The new Spectrum Roentgen Gamma (SRG) eROSITA X-ray telescope has a superior response to extended soft X-rays, which makes it ideal for studying low X-ray surface brightness emission of cosmic filaments. We search for extended X-ray emission between the two nearby galaxy clusters Abell~3667 and Abell~3651, which are separated by a projected transverse distance of ${\sim}13\,\mathrm{Mpc}$, using data from the SRG/eROSITA All-Sky Survey. We performed a detailed X-ray image analysis of the region between the two galaxy clusters and conducted a redshift analysis of the sources between them. We carried out a thorough surface brightness and spectral analysis between the clusters. The analysis was complemented with an X-ray pointed observation from XMM-Newton, infrared 2MASS data, and redshift information from NED. We discover an emission filament beyond the known radio relic northwest of A3667 and even beyond three times its virial radius. It is smoothly connected to A3651. The X-ray emission in the direction of the filament shows an enhancement of $30\pm3\,\%$ with a significance of $11\,\sigma$. The 2MASS map and redshift analysis show an alignment of galaxies along the filament and make a projection effect unlikely. Taking the redshift progression of galaxies within the filament into account, we estimate its three-dimensional length to be in the range of $25\,\mathrm{Mpc} - 32\,\mathrm{Mpc}$. The surface brightness analysis in combination with the temperature $T=({0.91}_{-0.11}^{+0.07})\,\mathrm{keV}$ and metallicity $Z=({0.10}_{-0.08}^{+0.05})\,\mathrm{Z_{\odot}}$ from the spectral analysis leads to estimates of a total flux, gas mass, and central baryon overdensity of $F_\mathrm{X}= (7.4\pm1.2)\times 10^{-12}\,\mathrm{erg\,s^{-1}\,cm^{-2}}$, $M_\mathrm{g}=(2.7^{+1.4}_{-0.8})\times 10^{14}\,\mathrm{M_\odot}$ and $\delta_0=215^{+86}_{-50}$., Comment: 12 pages and 13 figures in the main text plus 3 pages and 3 figures as appendix. Astronomy & Astrophysics, accepted on September 25, 2024. Now with included spectral analysis

The eROSITA telescope array aboard the Spektrum Roentgen Gamma (SRG) satellite began surveying the sky in December 2019, with the aim of producing all-sky X-ray source lists and sky maps of an unprecedented depth. Here we present catalogues of both point-like and extended sources using the data acquired in the first six months of survey operations (eRASS1; completed June 2020) over the half sky whose proprietary data rights lie with the German eROSITA Consortium. We describe the observation process, the data analysis pipelines, and the characteristics of the X-ray sources. With nearly 930000 entries detected in the most sensitive 0.2-2.3 keV energy range, the eRASS1 main catalogue presented here increases the number of known X-ray sources in the published literature by more than 60%, and provides a comprehensive inventory of all classes of X-ray celestial objects, covering a wide range of physical processes. A smaller catalogue of 5466 sources detected in the less sensitive but harder 2.3-5 keV band is the result of the first true imaging survey of the entire sky above 2 keV. We show that the number counts of X-ray sources in eRASS1 are consistent with those derived over narrower fields by past X-ray surveys of a similar depth, and we explore the number counts variation as a function of the location in the sky. Adopting a uniform all-sky flux limit (at 50% completeness) of F_{0.5-2 keV} > 5 \times 10^{-14}$ erg\,s$^{-1}$\,cm$^{-2}$, we estimate that the eROSITA all-sky survey resolves into individual sources about 20% of the cosmic X-ray background in the 1-2 keV range. The catalogues presented here form part of the first data release (DR1) of the SRG/eROSITA all-sky survey. Beyond the X-ray catalogues, DR1 contains all detected and calibrated event files, source products (light curves and spectra), and all-sky maps. Illustrative examples of these are provided., Comment: 39 pages, 23 figures. Accepted for publication in A&A. Accompanying eROSITA-DE Data Release 1

Cosmological studies have now entered Stage IV according to the Dark Energy Task Force prescription, thanks to new missions (Euclid, Rubin Observatory, SRG/eROSITA) that are expected to provide the required ultimate accuracy in the dark energy (DE) equation of state (EoS). However, none of these projects have the power to systematically unveil the galaxy cluster population at $z>1$. There therefore remains the need for an ATHENA-like mission to run independent cosmological investigations and scrutinise the consistency between the results from the $01$. Such samples will allow a detailed modelling of the evolution of cluster physics along with a standalone cosmological analysis. Our results suggest that survey B has the optimal design as it provides greater statistics. Remarkably, high-$z$ clusters, despite representing 15% or less of the full samples, allow a significant reduction of the uncertainty on the cosmological parameters: $\Delta w_a$ is reduced by a factor of 2.3 and $\Delta f_{NL}^{loc}$ by a factor of 3. Inventorying the high-$z$ X-ray cluster population can play a crucial role in ensuring overall cosmological consistency. This will be the major aim of future new-generation ATHENA-like missions., Comment: Accepted for publication in A&A, 17 pages, 18 figures

Aims. We aim to characterize the gas properties in the cluster outskirts ($R_{500}R_{200}$) of the A3391/95 system and to compare them to predictions. Methods. We performed X-ray image and spectral analyses using the eROSITA PV data to assess the gas morphology and properties in the outskirts and the filaments in the directions of the previously detected Northern and Southern Filament of the A3391/95 system. We took particular care of the foreground. Results. In the filament-facing outskirts of A3391 and the Northern Clump, we find higher temperatures than typical cluster outskirts profiles, with a significance of $1.6-2.8\sigma$, suggesting heating due to their connections with the filaments. We confirm SB excess in the profiles of the Northern, Eastern, and Southern Filaments. From spectral analysis, we detect hot gas of ~1 keV for the Northern and Southern Filaments. The filament metallicities are below 10% solar metallicity and the $n_e$ range between 2.6 and $6.3\times10^{-5}~\mathrm{cm^{-3}}$. The characteristic properties of the Little Southern Clump (~1.5$R_{200}$ from A3395S in the Southern Filament) suggest that it is a small galaxy group. Excluding the LSC from the analysis of the Southern Filament decreases the gas density by 30%. This shows the importance of taking into account any clumps to avoid overestimation of the gas measurement in the outskirts and filament regions. Conclusions. The $n_e$ of the filaments are consistent with the WHIM properties as predicted by cosmological simulations, but the temperatures are higher, close to the upper WHIM temperature limit. As both filaments are short and located in a denser environment, stronger gravitational heating may be responsible for this temperature enhancement. The metallicities are low, but still within the expected range from the simulations., Comment: Accepted for publication in Astronomy & Astrophysics: 17 pages, 6 figures (main text), 1 figure (appendix)

Three galaxy clusters selected from the XXL X-ray survey at high redshift and low mass ($z\sim1$ and $M_{500} \sim 1-2 \times 10^{14}$ M$_{\odot}$) were observed with NIKA2 to image their Sunyaev-Zel'dovich effect (SZ) signal. They all present an SZ morphology, together with the comparison with X-ray and optical data, that indicates dynamical activity related to merging events. Despite their disturbed intracluster medium, their high redshifts, and their low masses, the three clusters follow remarkably well the pressure profile and the SZ flux-mass relation expected from standard evolution. This suggests that the physics that drives cluster formation is already in place at $z \sim 1$ down to $M_{500} \sim 10^{14}$ M$_{\odot}$., Comment: to appear in Proc. of the mm Universe 2023 conference, Grenoble (France), June 2023, published by F. Mayet et al. (Eds), EPJ Web of conferences, EDP Sciences

The thermodynamical properties of the intracluster medium (ICM) are driven by scale-free gravitational collapse, but they also reflect the rich astrophysical processes at play in galaxy clusters. At low masses ($\sim 10^{14}$ M$_{\odot}$) and high redshift ($z \gtrsim 1$), these properties remain poorly constrained observationally, due to the difficulty in obtaining resolved and sensitive data. This paper aims at investigating the inner structure of the ICM as seen through the Sunyaev-Zel'dovich (SZ) effect in this regime of mass and redshift. Focus is set on the thermal pressure profile and the scaling relation between SZ flux and mass, namely the $Y_{\rm SZ} - M$ scaling relation. The three galaxy clusters XLSSC~072 ($z=1.002$), XLSSC~100 ($z=0.915$), and XLSSC~102 ($z=0.969$), with $M_{500} \sim 2 \times 10^{14}$ M$_{\odot}$, were selected from the XXL X-ray survey and observed with the NIKA2 millimeter camera to image their SZ signal. XMM-Newton X-ray data were used in complement to the NIKA2 data to derive masses based on the $Y_X - M$ relation and the hydrostatic equilibrium. The SZ images of the three clusters, along with the X-ray and optical data, indicate dynamical activity related to merging events. The pressure profile is consistent with that expected for morphologically disturbed systems, with a relatively flat core and a shallow outer slope. Despite significant disturbances in the ICM, the three high-redshift low-mass clusters follow remarkably well the $Y_{\rm SZ}-M$ relation expected from standard evolution. These results indicate that the dominant physics that drives cluster evolution is already in place by $z \sim 1$, at least for systems with masses above $M_{500} \sim 10^{14}$ M$_{\odot}$., Comment: 24 pages, published in A&A. Note that the title number had to be changed. This version matches the one from A&A

The notion of a Moore-Penrose inverse (M-P inverse) was introduced by Moore in 1920 and rediscovered by Penrose in 1955. The M-P inverse of a complex matrix is a special type of inverse which is unique, always exists, and can be computed using singular value decomposition. In a series of papers in the 1980s, Puystjens and Robinson studied M-P inverses more abstractly in the context of dagger categories. Despite the fact that dagger categories are now a fundamental notion in categorical quantum mechanics, the notion of a M-P inverse has not (to our knowledge) been revisited since their work. One purpose of this paper is, thus, to renew the study of M-P inverses in dagger categories. Here we introduce the notion of a Moore-Penrose dagger category and provide many examples including complex matrices, finite Hilbert spaces, dagger groupoids, and inverse categories. We also introduce generalized versions of singular value decomposition, compact singular value decomposition, and polar decomposition for maps in a dagger category, and show how, having such a decomposition is equivalent to having M-P inverses. This allows us to provide precise characterizations of which maps have M-P inverses in a dagger idempotent complete category, a dagger kernel category with dagger biproducts (and negatives), and a dagger category with unique square roots., Comment: In Proceedings QPL 2023, arXiv:2308.15489

Cartesian differential categories come equipped with a differential combinator which axiomatizes the fundamental properties of the total derivative from differential calculus. The objective of this paper is to understand when the Kleisli category of a monad is a Cartesian differential category. We introduce Cartesian differential monads, which are monads whose Kleisli category is a Cartesian differential category by way of lifting the differential combinator from the base category. Examples of Cartesian differential monads include tangent bundle monads and reader monads. We give a precise characterization of Cartesian differential categories which are Kleisli categories of Cartesian differential monads using abstract Kleisli categories. We also show that the Eilenberg-Moore category of a Cartesian differential monad is a tangent category., Comment: For the proceedings of MFPS2023

Previous work has shown that reverse differential categories give an abstract setting for gradient-based learning of functions between Euclidean spaces. However, reverse differential categories are not suited to handle gradient-based learning for functions between more general spaces such as smooth manifolds. In this paper, we propose a setting to handle this, which we call reverse tangent categories: tangent categories with an involution operation for their differential bundles., Comment: Accepted and published in the proceedings of CSL2024

This paper introduces a framework for solving alternating current optimal power flow (ACOPF) problems using graphics processing units (GPUs). While GPUs have demonstrated remarkable performance in various computing domains, their application in ACOPF has been limited due to challenges associated with porting sparse automatic differentiation (AD) and sparse linear solver routines to GPUs. We address these issues with two key strategies. First, we utilize a single-instruction, multiple-data abstraction of nonlinear programs. This approach enables the specification of model equations while preserving their parallelizable structure and, in turn, facilitates the parallel AD implementation. Second, we employ a condensed-space interior-point method (IPM) with an inequality relaxation. This technique involves condensing the Karush--Kuhn--Tucker (KKT) system into a positive definite system. This strategy offers the key advantage of being able to factorize the KKT matrix without numerical pivoting, which has hampered the parallelization of the IPM algorithm. By combining these strategies, we can perform the majority of operations on GPUs while keeping the data residing in the device memory only. Comprehensive numerical benchmark results showcase the advantage of our approach. Remarkably, our implementations -- MadNLP.jl and ExaModels.jl -- running on NVIDIA GPUs achieve an order of magnitude speedup compared with state-of-the-art tools running on contemporary CPUs., Comment: Accepted for publication in PSCC 2024

For decades, Simultaneous Ascending Auction (SAA) has been the most popular mechanism used for spectrum auctions. It has recently been employed by many countries for the allocation of 5G licences. Although SAA presents relatively simple rules, it induces a complex strategic game for which the optimal bidding strategy is unknown. Considering the fact that sometimes billions of euros are at stake in an SAA, establishing an efficient bidding strategy is crucial. In this work, we model the auction as a $n$-player simultaneous move game with complete information and propose the first efficient bidding algorithm that tackles simultaneously its four main strategic issues: the $\textit{exposure problem}$, the $\textit{own price effect}$, $\textit{budget constraints}$ and the $\textit{eligibility management problem}$. Our solution, called $SMS^\alpha$, is based on Simultaneous Move Monte Carlo Tree Search (SM-MCTS) and relies on a new method for the prediction of closing prices. By introducing a new reward function in $SMS^\alpha$, we give the possibility to bidders to define their own level of risk-aversion. Through extensive numerical experiments on instances of realistic size, we show that $SMS^\alpha$ largely outperforms state-of-the-art algorithms, notably by achieving higher expected utility while taking less risks., Comment: 20 pages, 19 figures, The paper has been submitted to IEEE journal for possible publication

The Virgo cluster is the nearest massive galaxy cluster and thus a prime target to study astrophysical processes in dense large-scale environments. In the radio band, we can probe the non-thermal components of the inter-stellar medium (ISM), intracluster medium (ICM) and of active galactic nuclei (AGN). With the ViCTORIA (Virgo Cluster multi-Telescope Observations in Radio of Interacting galaxies and AGN) project, we are carrying out multiple wide-field surveys of the Virgo cluster at different frequencies. We aim to investigate the impact of the environment on the evolution of galaxies and the contribution of AGN to the ICM-heating, from the inner cluster regions out to beyond the virial radius. We present a survey of the cluster at 120-168 MHz using LOFAR. We image a 132 deg$^2$ region of the cluster, reaching an order of magnitude greater sensitivity than existing wide-field radio surveys of this field at three times higher spatial resolution compared to other low-frequency observations. We developed a tailored data processing strategy to subtract the bright central radio galaxy M87 from the data. This allowed us to correct for the systematic effects due to ionospheric variation as a function of time and direction. In the final mosaic with a resolution of 9"x5", we reach a median noise level of 140 ${\mu}$Jy/beam inside the virial radius and 280 ${\mu}$Jy/beam for the full area. We detect 112 Virgo member galaxies and 114 background galaxies. In at least 18 cases, the radio morphology of the cluster member galaxies shows clear signs of ram-pressure stripping. This includes three previously unreported candidates. In addition, we reveal for the first time 150 kpc long tails from a previous epoch of AGN activity for NGC 4472 (M 49). While no cluster-scale diffuse radio sources are discovered, we find the presence of an extended radio signature of the W$'$-group., Comment: 18 pages, 12 figures, accepted for publication in A&A

A Hopf monad, in the sense of Brugui\`eres, Lack, and Virelizier, is a special kind of monad that can be defined for any monoidal category. In this note, we study Hopf monads in the case of a category with finite biproducts, seen as a symmetric monoidal category. We show that for biproducts, a Hopf monad is precisely characterized as a monad equipped with an extra natural transformation satisfying three axioms, which we call a fusion invertor. We will also consider three special cases: representable Hopf monads, idempotent Hopf monads, and when the category also has negatives. In these cases, the fusion invertor will always be of a specific form that can be defined for any monad. Thus in these cases, checking that a monad is a Hopf monad is reduced to checking one identity., Comment: Accepted in Theory and Applications of Categories (TAC)

In the category of sets and partial functions, $\mathsf{PAR}$, while the disjoint union $\sqcup$ is the usual categorical coproduct, the Cartesian product $\times$ becomes a restriction categorical analogue of the categorical product: a restriction product. Nevertheless, $\mathsf{PAR}$ does have a usual categorical product as well in the form $A \& B := A \sqcup B \sqcup (A \times B)$. Surprisingly, asking that a distributive restriction category (a restriction category with restriction products $\times$ and coproducts $\oplus$) has $A \& B$ a categorical product is enough to imply that the category is a classical restriction category. This is a restriction category which has joins and relative complements and, thus, supports classical Boolean reasoning. The first and main observation of the paper is that a distributive restriction category is classical if and only if $A \& B := A \oplus B \oplus (A \times B)$ is a categorical product in which case we call $\&$ the ''classical'' product. In fact, a distributive restriction category has a categorical product if and only if it is a classified restriction category. This is in the sense that every map $A \to B$ factors uniquely through a total map $A \to B \oplus \mathsf{1}$, where $\mathsf{1}$ is the restriction terminal object. This implies the second significant observation of the paper, namely, that a distributive restriction category has a classical product if and only if it is the Kleisli category of the exception monad $\_ \oplus \mathsf{1}$ for an ordinary distributive category. Thus having a classical product has a significant structural effect on a distributive restriction category. In particular, the classical product not only provides an alternative axiomatization for being classical but also for being the Kleisli category of the exception monad on an ordinary distributive category., Comment: We would like to dedicate this paper to Pieter Hofstra (1975-2022). Accepted in a special issue of Theory and Applications of Categories dedicated to Pieter Hofstra

In Linear Logic ($\mathsf{LL}$), the exponential modality $!$ brings forth a distinction between non-linear proofs and linear proofs, where linear means using an argument exactly once. Differential Linear Logic ($\mathsf{DiLL}$) is an extension of Linear Logic which includes additional rules for $!$ which encode differentiation and the ability of linearizing proofs. On the other hand, Graded Linear Logic ($\mathsf{GLL}$) is a variation of Linear Logic in such a way that $!$ is now indexed over a semiring $R$. This $R$-grading allows for non-linear proofs of degree $r \in R$, such that the linear proofs are of degree $1 \in R$. There has been recent interest in combining these two variations of $\mathsf{LL}$ together and developing Graded Differential Linear Logic ($\mathsf{GDiLL}$). In this paper we present a sequent calculus for $\mathsf{GDiLL}$, as well as introduce its categorical semantics, which we call graded differential categories, using both coderelictions and deriving transformations. We prove that symmetric powers always give graded differential categories, and provide other examples of graded differential categories. We also discuss graded versions of (monoidal) coalgebra modalities, additive bialgebra modalities, and the Seely isomorphisms, as well as their implementations in the sequent calculus of $\mathsf{GDiLL}$., Comment: In the proceedings of MFPS2023. Removed appendix from previous version to respect page limit. Minor corrections: the previous statement of one of our examples was incorrect, we thank Flavien Breuvart for explaining this to us. This has now been fixed. The rest of the paper remains unchanged

Tangent categories provide a categorical axiomatization of the tangent bundle. There are many interesting examples and applications of tangent categories in a variety of areas such as differential geometry, algebraic geometry, algebra, and even computer science. The purpose of this paper is to expand the theory of tangent categories in a new direction: the theory of operads. The main result of this paper is that both the category of algebras of an operad and its opposite category are tangent categories. The tangent bundle for the category of algebras is given by the semi-direct product, while the tangent bundle for the opposite category of algebras is constructed using the module of K\"ahler differentials, and these tangent bundles are in fact adjoints of one another. To prove these results, we first prove that the category of algebras of a coCartesian differential monad is a tangent category. We then show that the monad associated to any operad is a coCartesian differential monad. This also implies that we can construct Cartesian differential categories from operads. Therefore, operads provide a bountiful source of examples of tangent categories and Cartesian differential categories, which both recaptures previously known examples and also yield new interesting examples. We also discuss how certain basic tangent category notions recapture well-known concepts in the theory of operads., Comment: Fixed a minor typo, added table of contents

In this paper, we explain how the abstract notion of a differential bundle in a tangent category provides a new way of thinking about the category of modules over a commutative ring and its opposite category. MacAdam previously showed that differential bundles in the tangent category of smooth manifolds are precisely smooth vector bundles. Here we provide characterizations of differential bundles in the tangent categories of commutative rings and (affine) schemes. For commutative rings, the category of differential bundles over a commutative ring is equivalent to the category of modules over that ring. For affine schemes, the category of differential bundles over the Spec of a commutative ring is equivalent to the opposite category of modules over said ring. Finally, for schemes, the category of differential bundles over a scheme is equivalent to the opposite category of quasi-coherent sheaves of modules over that scheme.

We investigate how to port the standard interior-point method to new exascale architectures for block-structured nonlinear programs with state equations. Computationally, we decompose the interior-point algorithm into two successive operations: the evaluation of the derivatives and the solution of the associated Karush-Kuhn-Tucker (KKT) linear system. Our method accelerates both operations using two levels of parallelism. First, we distribute the computations on multiple processes using coarse parallelism. Second, each process uses a SIMD/GPU accelerator locally to accelerate the operations using fine-grained parallelism. The KKT system is reduced by eliminating the inequalities and the state variables from the corresponding equations, to a dense matrix encoding the sensitivities of the problem's degrees of freedom, drastically minimizing the memory exchange. We demonstrate the method's capability on the supercomputer Polaris, a testbed for the future exascale Aurora system. Each node is equipped with four GPUs, a setup amenable to our two-level approach. Our experiments on the stochastic optimal power flow problem show that the method can achieve a 50x speed-up compared to the state-of-the-art method., Comment: 23 pages, 8 figures

The operation of large-scale infrastructure networks requires scalable optimization schemes. To guarantee safe system operation, a high degree of feasibility in a small number of iterations is important. Decomposition schemes can help to achieve scalability. In terms of feasibility, however, classical approaches such as the alternating direction method of multipliers (ADMM) often converge slowly. In this work, we present primal decomposition schemes for hierarchically structured strongly convex QPs. These schemes offer high degrees of feasibility in a small number of iterations in combination with global convergence guarantees. We benchmark their performance against the centralized off-the-shelf interior-point solver Ipopt and ADMM on problems with up to 300,000 decision variables and constraints. We find that the proposed approaches solve problems as fast as Ipopt, but with reduced communication and without requiring a full model exchange. Moreover, the proposed schemes achieve a higher accuracy than ADMM.

Cartesian differential categories come equipped with a differential combinator that formalizes the directional derivative from multivariable calculus. Cartesian differential categories provide a categorical semantics of the differential lambda-calculus and have also found applications in causal computation, incremental computation, game theory, differentiable programming, and machine learning. There has recently been a desire to provide a (coordinate-free) characterization of Jacobians and gradients in Cartesian differential categories. One's first attempt might be to consider Cartesian differential categories which are Cartesian closed, such as models of the differential lambda-calculus, and then take the curry of the derivative. Unfortunately, this approach excludes numerous important examples of Cartesian differential categories such as the category of real smooth functions. In this paper, we introduce linearly closed Cartesian differential categories, which are Cartesian differential categories that have an internal hom of linear maps, a bilinear evaluation map, and the ability to curry maps which are linear in their second argument. As such, the Jacobian of a map is defined as the curry of its derivative. Many well-known examples of Cartesian differential categories are linearly closed, such as, in particular, the category of real smooth functions. We also explain how a Cartesian closed differential category is linearly closed if and only if a certain linear idempotent on the internal hom splits. To define the gradient of a map, one must be able to define the transpose of the Jacobian, which can be done in a Cartesian reverse differential category. Thus, we define the gradient of a map to be the curry of its reverse derivative and show this equals the transpose of its Jacobian. We also explain how a linearly closed Cartesian reverse differential category is precisely a linearly closed Cartesian differential category with an appropriate notion of transpose., Comment: In Proceedings ACT 2021, arXiv:2211.01102

Cartesian differential categories come equipped with a differential operator which formalises the total derivative from multivariable calculus. Cofree Cartesian differential categories always exist over a specified base category, where the general construction is based on Faa di Bruno's formula. A natural question to ask is, when given an arbitrary Cartesian differential category, how can one check if it is cofree without knowing the base category? In this paper, we provide characterisations of cofree Cartesian differential categories without specifying a base category. The key to these characterisations is, surprisingly, maps whose derivatives are zero, which we call differential constants. One characterisation is in terms of the homsets being complete ultrametric spaces, where the ultrametric is induced by differential constants, which is similar to the metric for power series. Another characterisation is as algebras of a monad. In either characterisation, the base category is the category of differential constants. We also discuss other basic properties of cofree Cartesian differential categories, such as the linear maps, and explain how many well-known Cartesian differential categories (such as polynomial or smooth functions) are not cofree.

This article presents a constrained policy optimization approach for the optimal control of systems under nonstationary uncertainties. We introduce an assumption that we call Markov embeddability that allows us to cast the stochastic optimal control problem as a policy optimization problem over the augmented state space. Then, the infinite-dimensional policy optimization problem is approximated as a finite-dimensional nonlinear program by applying function approximation, deterministic sampling, and temporal truncation. The approximated problem is solved by using automatic differentiation and condensed-space interior-point methods. We formulate several conceptual and practical open questions regarding the asymptotic exactness of the approximation and the solution strategies for the approximated problem. As a proof of concept, we provide a numerical example demonstrating the performance of the control policy obtained by the proposed method.

We report numerical results on solving constrained linear-quadratic model predictive control (MPC) problems by exploiting graphics processing units (GPUs). The presented method reduces the MPC problem by eliminating the state variables and applies a condensed-space interior-point method to remove the inequality constraints in the KKT system. The final condensed matrix is positive definite and can be efficiently factorized in parallel on GPU/SIMD architectures. In addition, the size of the condensed matrix depends only on the number of controls in the problem, rendering the method particularly effective when the problem has many states but few inputs and moderate horizon length. Our numerical results for PDE-constrained problems show that the approach is an order of magnitude faster than a standard CPU implementation. We also provide an open-source Julia framework that facilitates modeling (DynamicNLPModels.jl) and solution (MadNLP.jl) of MPC problems on GPUs.

A traced monad is a monad on a traced symmetric monoidal category that lifts the traced symmetric monoidal structure to its Eilenberg-Moore category. A long-standing question has been to provide a characterization of traced monads without explicitly mentioning the Eilenberg-Moore category. On the other hand, a symmetric Hopf monad is a symmetric bimonad whose fusion operators are invertible. For compact closed categories, symmetric Hopf monads are precisely the kind of monads that lift the compact closed structure to their Eilenberg-Moore categories. Since compact closed categories and traced symmetric monoidal categories are closely related, it is a natural question to ask what is the relationship between Hopf monads and traced monads. In this paper, we introduce trace-coherent Hopf monads on traced monoidal categories, which can be characterized without mentioning the Eilenberg-Moore category. The main theorem of this paper is that a symmetric Hopf monad is a traced monad if and only if it is a trace-coherent Hopf monad. We provide many examples of trace-coherent Hopf monads, such as those induced by cocommutative Hopf algebras or any symmetric Hopf monad on a compact closed category. We also explain how for traced Cartesian monoidal categories, trace-coherent Hopf monads can be expressed using the Conway operator, while for traced coCartesian monoidal categories, any trace-coherent Hopf monad is an idempotent monad. We also provide separating examples of traced monads that are not Hopf monads, as well as symmetric Hopf monads that are not trace-coherent., Comment: Final Version, Published in Compositionality