Researcher profile

Alessandro Lupi

Alessandro Lupi contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
20works
0followers
7topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

20 published item(s)

preprint2025arXiv

Super-Eddington accretion in high-redshift quasar hosts: Black-hole driven outflows, galaxy quenching, and the nature of little red dots

The advent of the James Webb Space Telescope has revolutionised our understanding of the high-redshift Universe through its detection of bright, massive galaxies up to $z\gtrsim 10$ and its identification of peculiar sources called `little red dots' (LRDs). The origin of both classes of objects remains uncertain but is likely linked to the formation and early growth of the first massive black holes (MBHs), which may be more easily explained by invoking phases of super-Eddington accretion. In this study, we used a state-of-the-art zoom-in cosmological simulation of a quasar host to investigate whether these objects could resemble any of the peculiar sources observed with JWST during their assembly. We find that the impact of MBH feedback on star formation is typically moderate, with outflows preferentially escaping perpendicular to the galactic disc. However, for approximately ten percent of the galaxy's lifetime, the system enters a distinct quenched phase following rapid MBH growth driven by super-Eddington accretion. This phase culminates in a powerful feedback event, during which the MBH jet and disc-driven winds interact directly with the galactic disc and carve out a central cavity. We also find that, during the history of the quasar host progenitor, the spectral properties of the system can resemble both LRDs and quenched galaxies, depending on the specific evolutionary stage considered. These findings suggest that both conditions may represent transient phases in the life cycle of high-redshift galaxies.

preprint2023arXiv

A general relativistic extension to mesh-free methods for hydrodynamics

The detection of gravitational waves has opened a new era for astronomy, allowing for the combined use of gravitational wave and electromagnetic emissions to directly probe the physics of compact objects, still poorly understood. So far, the theoretical modelling of these sources has mainly relied on standard numerical techniques as grid-based methods or smoothed particle hydrodynamics, with only a few recent attempts at using new techniques as moving-mesh schemes. Here, we introduce a general relativistic extension to the mesh-less hydrodynamic schemes in the code GIZMO, which benefits from the use of Riemann solvers and at the same time perfectly conserves angular momentum thanks to a generalised leap-frog integration scheme. We benchmark our implementation against many standard tests for relativistic hydrodynamics, either in one or three dimensions, and also test the ability to preserve the equilibrium solution of a Tolman-Oppenheimer-Volkoff compact star. In all the presented tests, the code performs extremely well, at a level at least comparable to other numerical techniques.

preprint2022arXiv

Black hole spin evolution in warped accretion discs

Massive black holes (BHs) inhabiting galactic nuclei can be described by two parameters only, i.e. mass and spin, that change through cosmic time in response to accretion and merger events. While most numerical simulations accurately track the BH mass, spin evolution is rarely taken into account. In this work, we implement and validate a self-consistent sub-grid model for the evolution of the BH mass and spin via gas accretion in the hydrodynamics code GIZMO. The model assumes that accretion from resolved scales does not occur instantaneously, but is mediated by a sub-grid geometrically thin $α$-disc. After validating our model semi-analytically, we test it in an idealized environment consisting of a circumnuclear disc, where gas accretion onto the accretion disc is consistently determined by GIZMO. In the absence of any accretion-related feedback, the spin evolution closely traces that observed in the semi-analytical models and depends on the free parameters of our implementation, such as the initial BH spin, angular momentum of the accretion disc, and the radius at which the gas inflow circularises. In GIZMO, we also couple our model with the biconical-outflow model presented in a companion paper, wherein the feedback axis is always aligned with the BH spin. In this last case, the evolution of the central BH differs significantly from the previous cases, since the feedback process modifies the gas dynamics and its inflow rates from resolved scales. Such an interaction cannot be modeled by simple semi-analytical models and should be treated using full N-body hydrodynamical simulations.

preprint2022arXiv

Disc instability and bar formation: view from the IllustrisTNG simulations

We make use of z = 0 samples of strongly barred and unbarred disc galaxies from the TNG100 and TNG50 cosmological hydrodynamical simulations to assess the performance of the simple disc instability criterion proposed by Efstathiou, Lake & Negroponte (1982) (ELN-criterion). We find that strongly barred galaxies generally assemble earlier, are more star-dominated in their central regions, and have more massive and more compact discs than unbarred galaxies. The ELN-criterion successfully identifies ~75% and ~80% of the strongly barred and the unbarred galaxies, respectively. Strongly barred galaxies that the criterion fails to identify tend to have more extended discs, higher spin values and bars that assembled later than is typical for the bulk of the barred population. The bars in many of these cases appear to be produced by an interaction with a close neighbour (i.e. to be externally triggered) rather than to result from secular growth in the disc. On the other hand, we find that unbarred galaxies misclassified as barred by the ELN-criterion typically have stellar discs similar to those of barred galaxies, although more extended in the vertical direction and less star-dominated in their central regions, possibly reflecting later formation times. In addition, the bulge component of these galaxies is significantly more prominent at early times than in the strongly barred sample. Thus, the ELN-criterion robustly identifies secular bar instabilities in most simulated disc galaxies, but additional environmental criteria are needed to account for interaction-induced bar formation.

preprint2022arXiv

Enhanced star formation in $z\sim6$ quasar companions

Quasars powered by supermassive black holes (MBH, $>10^8~M_{\odot}$) at $z\sim 6$ are predicted to reside in cosmic over-dense regions. However, observations so far could not confirm this expectation due to limited statistics. The picture is further complicated by the possible effects of quasar outflows (i.e. feedback) that could either suppress or stimulate the star formation rate (SFR) of companion galaxies, thus modifying the expected bias. Here we quantify feedback effects on the properties and detectability of companions by comparing cosmological zoom-in simulations of a quasar in which feedback is either included or turned-off. With respect to the no-feedback case, companions (a) directly impacted by the outflow have their SFR increased by a factor $2-3$, and (b) tend to be more massive. Both effects shift the [CII]158$μ$m and UV luminosity functions toward brighter magnitudes. This leads us to conclude that quasar feedback slightly increases the effective quasar bias, boosting the number density of observable quasar companions, in agreement with what has been found around the brightest quasars of recent ALMA [CII] surveys. Deeper observations performed with JWST and/or ALMA will improve the statistical significance of this result by detecting a larger number of fainter quasar companions.

preprint2022arXiv

Forming massive seed black holes in high-redshift quasar host progenitors

The presence of massive black holes (BHs) with masses of order $10^9\rm\, M_\odot$, powering bright quasars when the Universe was less than 1 Gyr old, poses strong constraints on their formation mechanism. Several scenarios have been proposed to date to explain massive BH formation, from the low-mass seed BH remnants of the first generation of stars to the massive seed BHs resulting from the rapid collapse of massive gas clouds. However, the plausibility of some of these scenarios to occur within the progenitors of high-z quasars has not yet been thoroughly explored. In this work, we investigate, by combining dark-matter only N-body simulations with a semi-analytic framework, whether the conditions for the formation of massive seed BHs from synchronised atomic-cooling halo pairs and/or dynamically-heated mini-haloes are fulfilled in the overdense regions where the progenitors of a typical high-redshift quasar host form and evolve. Our analysis shows that the peculiar conditions in such regions, i.e. strong halo clustering and high star formation rates, are crucial to produce a non-negligible number of massive seed BH host candidates: we find $\approx1400$ dynamically heated metal-free mini-haloes, including one of these which evolves to a synchronised pair and ends up in the massive quasar-host halo by $z=6$. This demonstrates that the progenitors of high-redshift quasar host haloes can harbour early massive seed BHs. Our results further suggest that multiple massive seed BHs may form in or near the quasar host's progenitors, potentially merging at lower redshifts and yielding gravitational wave events.

preprint2022arXiv

High-redshift quasars and their host galaxies II: multiphase gas and stellar kinematics

Observations of $z \gtrsim 6$ quasars provide information on the early phases of the most massive black holes (MBHs) and galaxies. Current observations at sub-mm wavelengths trace cold and warm gas, and future observations will extend information to other gas phases and the stellar properties. The goal of this study is to examine the gas life cycle in a $z \gtrsim 6$ quasar: from accretion from the halo to the galaxy and all the way into the MBH, to how star formation and the MBH itself affect the gas properties. Using a very-high resolution cosmological zoom-in simulation of a $z=7$ quasar including state-of-the-art non-equilibrium chemistry, MBH formation, growth and feedback, we investigate the distribution of the different gas phases in the interstellar medium across cosmic time. We assess the morphological evolution of the quasar host using different tracers (star- or gas-based) and the thermodynamic distribution of the MBH accretion-driven outflows, finding that obscuration in the disc is mainly due to molecular gas, with the atomic component contributing at larger scales and/or above/below the disc plane. Moreover, our results also show that molecular outflows, if present, are more likely the result of gas being lifted near the MBH than production within the wind because of thermal instabilities. Finally, we also discuss how different gas phases can be employed to dynamically constrain the MBH mass, and argue that resolutions below $\sim 100$ pc yield unreliable estimates because of the strong contribution of the nuclear stellar component to the potential at larger scales.

preprint2022arXiv

Maximally informed Bayesian modelling of disc galaxies

Dissecting the underlying structure of galaxies is of main importance in the framework of galaxy formation and evolution theories. While a classical bulge+disc decomposition of disc galaxies is usually taken as granted, this is only rarely solidly founded upon the full exploitation of the richness of data arising from spectroscopic studies with integral field units. In this work we describe a fully Bayesian estimation method of the global structure of disc galaxies which makes use of the wealth of photometric, kinematic, and mass-to-light ratio data, and that can be seen as a first step towards a machine-learning approach, certainly needed when dealing with larger samples of galaxies. Ours is a novel, hybrid line of action in tackling the problem of galactic parameter estimation, neither purely photometric nor orbit-based. Being rooted on a nested sampler, our code, which is available publicly as an online repository, allows for a statistical assessment of the need for multiple components in the dissecting process. As a first case-study the GPU-optimized code is applied to the S0 galaxy NGC-7683, finding that in this galaxy a pseudo-bulge, possibly the remnant of a bar-like structure, does exist in the center of the system. These results are then tested against the publicly available, orbit-based code DYNAMITE, finding substantial agreement.

preprint2022arXiv

Morphological decomposition of TNG50 galaxies: methodology and catalogue

We present MORDOR (MORphological DecOmposeR, a new algorithm for structural decomposition of simulated galaxies based on stellar kinematics. The code measures the properties of up to five structural components (a thin/cold and a thick/warm disc, a classical and a secular bulge, and a spherical stellar halo), and determines the properties of a stellar bar (if present). A comparison with other algorithms presented in the literature yields overall good agreement, with MORDOR displaying a higher flexibility in correctly decomposing systems and identifying bars in crowded environments (e.g. with ongoing fly-bys, often observable in cosmological simulations). We use MORDOR to analyse galaxies in the TNG50 simulation and find the following: ($i$) the thick disc component undergoes the strongest evolution in the binding energy-circularity plane, as expected when disc galaxies decrease their turbulent-rotational support with cosmic time; ($ii$) smaller galaxies (with stellar mass, $10^{9} \lesssim M_{*} / {\rm M_{\odot}} \leq 5 \times 10^{9}$) undergo a major growth in their disc components after $z\sim 1$, whereas ($iii$) the most massive galaxies ($5 \times 10^{10} < M_{*} / {\rm M_{\odot}} \leq 5\times10^{11}$) evolve toward more spheroidal dominated objects down to $z=0$ due to frequent gravitational interactions with satellites; ($iv$) the fraction of barred galaxies grows rapidly at high redshift and stabilizes below $z\sim 2$, except for the most massive galaxies that show a decrease in the bar occupation fraction at low redshift; ($v$) galaxies with $M_{*} \sim 10^{11}~{\rm M_{\odot}}$ exhibit the highest relative occurrence of bars at $z=0$, in agreement with observational studies. We publicly release MORDOR and the morphological catalogue of TNG50 galaxies.

preprint2022arXiv

Non-isotropic feedback from accreting spinning black holes

Active galactic nuclei (AGN) are massive black holes (BHs) caught in the act of accreting gas at the centre of their host galaxies. Part of the accreting mass is converted to energy and released into the surrounding medium, in a process loosely referred to as AGN feedback. Most numerical simulations include AGN feedback as a sub-grid model, wherein energy or momentum (or both) is coupled to the nearby gas. In this work, we implement a new momentum-driven model in the hydrodynamics code GIZMO, in which accretion from large scales is mediated by a sub-grid accretion disc model, and gas particles are stochastically kicked over a bi-conical region, to mimic observed kinetic winds. The feedback cone&#39;s axis can be set either parallel to the angular momentum of the gas surrounding the BH or to the BH spin direction, which is self-consistently evolved within the accretion-disc model. Using a circumnuclear disc (CND) as a test bed, we find that (i) the conical shape of the outflow is always visible and is weakly dependent on the launching orientation and aperture, resulting in comparable mass inflows and outflows; (ii) the cone&#39;s orientation is also similar amongst our tests, and it is not always the same as the initial value, due to the interaction with the CND playing a crucial role in shaping the outflow; and (iii) the velocity of the outflow, instead, differs and strongly depends on the interplay with the CND.

preprint2022arXiv

On the low ortho-to-para H2 ratio in star-forming filaments

The formation of stars and planetary systems is a complex phenomenon, which relies on the interplay of multiple physical processes. Nonetheless, it represents a crucial stage for our understanding of the Universe, and in particular of the conditions leading to the formation of key molecules (e.g. water) on comets and planets. Herschel observations demonstrated that stars form out of gaseous filamentary structures in which the main constituent is molecular hydrogen (H$_2$). Depending on its nuclear spin H$_2$ can be found in two forms: `ortho&#39; with parallel spins and `para&#39; where the spins are anti-parallel. The relative ratio among these isomers, i.e. the ortho-to-para ratio (OPR), plays a crucial role in a variety of processes related to the thermodynamics of star-forming gas and to the fundamental chemistry affecting the deuteration of water in molecular clouds, commonly used to determine the origin of water in Solar System&#39;s bodies. Here, for the first time, we assess the evolution of the OPR starting from the warm neutral medium, by means of state-of-the-art three-dimensional magneto-hydrodynamic simulations of turbulent molecular clouds. Our results show that star-forming clouds exhibit a low OPR ($\ll 0.1$) already at moderate densities ($\sim$1000 cm$^{-3}$). We also constrain the cosmic rays ionisation rate, finding that $10^{-16}\,\rm s^{-1}$ is the lower limit required to explain the observations of diffuse clouds. Our results represent a step forward in the understanding of the star and planet formation process providing a robust determination of the chemical initial conditions for both theoretical and observational studies.

preprint2022arXiv

The evolution of the barred galaxy population in the TNG50 simulation

We use the magnetic-hydrodynamical simulation TNG50 to study the evolution of barred massive disc galaxies. Massive spiral galaxies are already present as early as $z=4$, and bar formation takes place already at those early times. The bars grow longer and stronger as the host galaxies evolve, with the bar sizes increasing at a pace similar to that of the disc scale lengths. The bar fraction mildly evolves with redshift for galaxies with $M_{*}\geq10^{10}M\odot$, being greater than $\sim40\%$ at $0.5<z<3$ and $\sim30\%$ at $z=0$. When bars larger than a given physical size ($\geq 2\,\rm kpc$) or the angular resolution limit of twice the I-band angular PSF FWHM of the HST are considered, the bar fraction dramatically decreases with increasing redshift, reconciling the theoretical predictions with observational data. We find that barred galaxies have an older stellar population, lower gas fractions and star formation rates than unbarred galaxies. In most cases, the discs of barred galaxies assembled earlier and faster than the discs of unbarred galaxies. We also find that barred galaxies are typical in haloes with larger concentrations and smaller spin parameters than unbarred galaxies. Furthermore, the inner regions of barred galaxies are more baryon-dominated than those of unbarred galaxies but have comparable global stellar mass fractions. Our findings suggest that the bar population could be used as a potential tracer of the buildup of disc galaxies and their host haloes. With this paper, we release a catalogue of barred galaxies in TNG50 at $6$ redshifts between $z=4$ and $z=0$.

preprint2022arXiv

The role of bars on the dynamical-friction driven inspiral of massive perturbers

In this paper, we systematically explore the impact of a galactic bar on the inspiral time-scale of a massive object (MO) within a Milky Way-like galaxy. We integrate the orbit of MOs in a multi-component galaxy model via a semi-analytical approach that accounts for dynamical friction generalized to rotationally supported backgrounds. We compare the MO evolution in a galaxy featuring a Milky Way-like rotating bar to the evolution within an analogous axisymmetric galaxy without the bar. In agreement with previous studies, we find that the bar presence may significantly affect the inspiral, sometimes making it shorter by a factor of a few, sometimes hindering it for a Hubble time. The erratic behaviour is mainly impacted by the relative phase at which the MO encounters the stronger bar-induced resonances. In particular, the effect of the bar is more prominent for initially in-plane, prograde MOs, especially those crossing the bar co-rotation radius or outer Lindblad resonance. In the barred galaxy, we find the sinking of the most massive MOs (>~10^7.5 Msun) approaching the galaxy from large separations (>~8 kpc) to be most efficiently hampered. Neglecting the effect of global torques associated with the non-symmetric mass distribution is thus not advisable even within an idealized, smooth galaxy model; we further note that spiral patterns are unlikely to affect the inspiral due to their transient and fluctuating nature. We speculate that the sinking efficiency of massive black holes involved in minor galaxy mergers may be hampered in barred galaxies, making them less likely to host a gravitational wave signal accessible to low-frequency detectors.

preprint2021arXiv

Dynamical evolution of massive perturbers in realistic multi-component galaxy models I: implementation and validation

Galaxies are self-gravitating structures composed by several components encompassing spherical, axial and triaxial symmetry. Although real systems feature heterogeneous components whose properties are intimately connected, semi-analytical approaches often exploit the linearity of the Poisson&#39;s equation to represent the potential and mass distribution of a multi-component galaxy as the sum of the individual components. In this work, we expand the semi-analytical framework developed in Bonetti et al. (2020) by including both a detailed implementation of the gravitational potential of exponential disc (modelled with a ${\rm sech}^2$ and an exponential vertical profile) and an accurate prescription for the dynamical friction experienced by massive perturbers in composite galaxy models featuring rotating disc structures. Such improvements allow us to evolve arbitrary orbits either within or outside the galactic disc plane. We validate the results obtained by our numerical model against public semi-analytical codes as well as full N-body simulations, finding that our model is in excellent agreement to the codes it is compared with. The ability to reproduce the relevant physical processes responsible for the evolution of massive perturber orbits and its computational efficiency make our framework perfectly suited for large parameter-space exploration studies.

preprint2020arXiv

A new proxy to estimate the cosmic-ray ionisation rate in dense cores

Cosmic rays are a global source of ionisation, and the ionisation fraction represents a fundamental parameter in the interstellar medium. Ions couple to magnetic fields, affect the chemistry, and the dynamics of star-forming regions as well as planetary atmospheres. However, the cosmic-ray ionisation rate represents one of the bottlenecks for astrochemical models, and its determination is one of the most puzzling problems in astrophysics. While for diffuse clouds reasonable values have been provided from H$_3^+$ observations, for dense clouds, due to the lack of rotational transitions, this is not possible, and estimates are strongly biased by the employed model. We present here an analytical expression, obtained from first principles, to estimate the cosmic-ray ionisation rate from observational quantities. The theoretical predictions are validated with high-resolution three-dimensional numerical simulations and applied to the well known core L1544; we obtained an estimate of $ζ_2 \sim 2-3 \times 10^{-17}$ s$^{-1}$. Our results and the analytical formulae provided represent the first model-independent, robust tool to probe the cosmic-ray ionisation rate in the densest part of star-forming regions (on spatial scales of $R \leq 0.05$ pc). An error analysis is presented to give statistical relevance to our study.

preprint2020arXiv

Difficulties in Mid-Infrared selection of AGN in dwarf galaxies

While massive black holes (MBHs) are known to inhabit all massive galaxies, their ubiquitous presence in dwarf galaxies has not been confirmed yet, with only a limited number of sources detected so far. Recently, some studies proposed infrared emission as an alternative way to identify MBHs in dwarfs, based on a similar approach usually applied to quasars. In this study, by accurately combining optical and infrared data taking into account resolution effects and source overlapping, we investigate in detail the possible limitations of this approach with current ground-based facilities, finding a quite low ($\sim$0.4 per cent) fraction of active MBH in dwarfs that are luminous in mid-infrared, consistent with several previous results. Our results suggest that the infrared selection is strongly affected by several limitations that make the identification of MBHs in dwarf galaxies currently prohibitive, especially because of the very poor resolution compared to optical surveys, and the likely contamination by nearby sources, although we find a few good candidates worth further follow-ups. Optical, X-ray and radio observations, therefore, still represent the most secure way to search for MBH in dwarfs.

preprint2020arXiv

Dynamical friction-driven orbital circularisation in rotating discs: a semi-analytical description

We present and validate a novel semi-analytical approach to study the effect of dynamical friction on the orbits of massive perturbers in rotating stellar discs. We find that dynamical friction efficiently circularises the orbit of co-rotating perturbers, while it constantly increases the eccentricity of counter-rotating ones until their angular momenta reverse, then once again promoting circularisation. Such &#34;drag toward circular corotation&#34; could shape the distribution of orientations of kinematically decoupled cores in disc galaxies, naturally leading to the observed larger fraction of co-rotating cores.

preprint2020arXiv

Globular Cluster Formation from Colliding Substructure

We investigate a scenario where the formation of Globular Clusters (GCs) is triggered by high-speed collisions between infalling atomic-cooling subhalos during the assembly of the main galaxy host, a special dynamical mode of star formation that operates at high gas pressures and is intimately tied to LCDM hierarchical galaxy assembly. The proposed mechanism would give origin to &#34;naked&#34; globulars, as colliding dark matter subhalos and their stars will simply pass through one another while the warm gas within them clashes at highly supersonic speed and decouples from the collisionless component, in a process reminiscent of the Bullet galaxy cluster. We find that the resulting shock-compressed layer cools on a timescale that is typically shorter than the crossing time, first by atomic line emission and then via fine-structure metal-line emission, and is subject to gravitational instability and fragmentation. Through a combination of kinetic theory approximation and high-resolution $N$-body simulations, we show that this model may produce: (a) a GC number-halo mass relation that is linear down to dwarf galaxy scales and agrees with the trend observed over five orders of magnitude in galaxy mass; (b) a population of old globulars with a median age of 12 Gyr and an age spread similar to that observed; (c) a spatial distribution that is biased relative to the overall mass profile of the host; and (d) a bimodal metallicity distribution with a spread similar to that observed in massive galaxies.

preprint2020arXiv

The [CII]--SFR correlation in dwarf galaxies across cosmic time

Current galaxy observations suggest that a roughly linear correlation exists between the [CII] emission and the star formation rate, either as spatially-resolved or integrated quantities. Observationally, this correlation seems to be independent of metallicity, but the very large scatter does not allow to properly assess whether this is true. On the other hand, theoretical models tend to suggest a metallicity dependence of the correlation. In this study, we investigate the metallicity evolution of the correlation via a high-resolution zoom-in cosmological simulation of a dwarf galaxy employing state-of-the-art sub-grid modelling for gas cooling, star formation, and stellar feedback, and that self-consistently evolves the abundances of metal elements out of equilibrium. Our results suggest that the correlation should evolve with metallicity, in agreement with theoretical predictions, but also that this evolution can be hardly detected in observations, because of the large scatter. We also find that most of the [CII] emission is associated with neutral gas at low-intermediate densities, whereas the highest emissivity is produced by the densest regions around star-forming regions.

preprint2020arXiv

The AGORA high-resolution galaxy simulations comparison project: Public data release

As part of the AGORA High-resolution Galaxy Simulations Comparison Project (Kim et al. 2014, 2016) we have generated a suite of isolated Milky Way-mass galaxy simulations using 9 state-of-the-art gravito-hydrodynamics codes widely used in the numerical galaxy formation community. In these simulations we adopted identical galactic disk initial conditions, and common physics models (e.g., radiative cooling and ultraviolet background by a standardized package). Subgrid physics models such as Jeans pressure floor, star formation, supernova feedback energy, and metal production were carefully constrained. Here we release the simulation data to be freely used by the community. In this release we include the disk snapshots at 0 and 500Myr of evolution per each code as used in Kim et al. (2016), from simulations with and without star formation and feedback. We encourage any member of the numerical galaxy formation community to make use of these resources for their research - for example, compare their own simulations with the AGORA galaxies, with the common analysis yt scripts used to obtain the plots shown in our papers, also available in this release.