Search

One of the most challenging problems in the computational study of localization in quantum manybody systems is to capture the effects of rare events, which requires sampling over exponentially many disorder realizations. We implement an efficient procedure on a quantum processor, leveraging quantum parallelism, to efficiently sample over all disorder realizations. We observe localization without disorder in quantum many-body dynamics in one and two dimensions: perturbations do not diffuse even though both the generator of evolution and the initial states are fully translationally invariant. The disorder strength as well as its density can be readily tuned using the initial state. Furthermore, we demonstrate the versatility of our platform by measuring Renyi entropies. Our method could also be extended to higher moments of the physical observables and disorder learning.

We present a suite of six high-resolution chemo-dynamical simulations of isolated galaxies, spanning observed disk-dominated environments on the star-forming main sequence, as well as quenched, bulge-dominated environments. We compare and contrast the physics driving star formation and stellar feedback amongst the galaxies, with a view to modeling these processes in cosmological simulations. We find that the mass-loading of galactic outflows is coupled to the clustering of supernova explosions, which varies strongly with the rate of galactic rotation $\Omega = v_c/R$ via the Toomre length, leading to smoother gas disks in the bulge-dominated galaxies. This sets an equation of state in the star-forming gas that also varies strongly with $\Omega$, so that the bulge-dominated galaxies have higher mid-plane densities, lower velocity dispersions, and higher molecular gas fractions than their main sequence counterparts. The star formation rate in five out of six galaxies is independent of $\Omega$, and is consistent with regulation by the mid-plane gas pressure alone. In the sixth galaxy, which has the most centrally-concentrated bulge and thus the highest $\Omega$, we reproduce dynamical suppression of the star formation efficiency (SFE) in agreement with observations. This produces a transition away from pressure-regulated star formation., Comment: This is a Learning the Universe Publication, accepted for publication in ApJ. 29 pages, 14 figures, and analysis code at: https://github.com/sjeffreson/pressure_regulated_SF_analysis

Quantum error correction provides a path to reach practical quantum computing by combining multiple physical qubits into a logical qubit, where the logical error rate is suppressed exponentially as more qubits are added. However, this exponential suppression only occurs if the physical error rate is below a critical threshold. In this work, we present two surface code memories operating below this threshold: a distance-7 code and a distance-5 code integrated with a real-time decoder. The logical error rate of our larger quantum memory is suppressed by a factor of $\Lambda$ = 2.14 $\pm$ 0.02 when increasing the code distance by two, culminating in a 101-qubit distance-7 code with 0.143% $\pm$ 0.003% error per cycle of error correction. This logical memory is also beyond break-even, exceeding its best physical qubit's lifetime by a factor of 2.4 $\pm$ 0.3. We maintain below-threshold performance when decoding in real time, achieving an average decoder latency of 63 $\mu$s at distance-5 up to a million cycles, with a cycle time of 1.1 $\mu$s. To probe the limits of our error-correction performance, we run repetition codes up to distance-29 and find that logical performance is limited by rare correlated error events occurring approximately once every hour, or 3 $\times$ 10$^9$ cycles. Our results present device performance that, if scaled, could realize the operational requirements of large scale fault-tolerant quantum algorithms., Comment: 10 pages, 4 figures, Supplementary Information

Early-type galaxies (ETGs) are known to harbour dense spheroids of stars but scarce star formation (SF). Approximately a quarter of these galaxies have rich molecular gas reservoirs yet do not form stars efficiently. We study here the ETG NGC~524, with strong shear suspected to result in a smooth molecular gas disc and low star-formation efficiency (SFE). We present new spatially-resolved observations of the \textsuperscript{12}CO(2-1)-emitting cold molecular gas from the Atacama Large Millimeter/sub-millimeter Array (ALMA) and of the warm ionised-gas emission lines from SITELLE at the Canada-France-Hawaii Telescope. Although constrained by the resolution of the ALMA observations ($\approx37$~pc), we identify only $52$ GMCs with radii ranging from $30$ to $140$~pc, a low mean molecular gas mass surface density $\langle\Sigma_{\rm gas}\rangle\approx125$~M$_\odot$~pc$^{-2}$ and a high mean virial parameter $\langle\alpha_{\rm obs,vir}\rangle\approx5.3$. We measure spatially-resolved molecular gas depletion times ($\tau_{\rm dep}\equiv1/{\rm SFE}$) with a spatial resolution of $\approx100$~pc within a galactocentric distance of $1.5$~kpc. The global depletion time is $\approx2.0$~Gyr but $\tau_{\rm dep}$ increases toward the galaxy centre, with a maximum $\tau_{\rm dep,max}\approx5.2$~Gyr. However, no pure \ion{H}{II} region is identified in NGC~524 using ionised-gas emission-line ratio diagnostics, so the $\tau_{\rm dep}$ inferred are in fact lower limits. Measuring the GMC properties and dynamical states, we conclude that shear is the dominant mechanism shaping the molecular gas properties and regulating SF in NGC~524. This is supported by analogous analyses of the GMCs in a simulated ETG similar to NGC~524., Comment: 17 pages, 10 figures. To be published in MNRAS, accepted on May 27

NGC~613 is a nearby barred spiral galaxy with a nuclear ring. Exploiting high spatial resolution ($\approx20$ pc) Atacama Large Millimeter/sub-millimeter Array $^{12}$CO(1-0) observations, we study the giant molecular clouds (GMCs) in the nuclear ring and its vicinity, identifying $158$ spatially- and spectrally-resolved GMCs. The GMC sizes ($R_{\mathrm{c}}$) are comparable to those of the clouds in the Milky Way (MW) disc, but their gas masses, observed linewidths ($\sigma_{\mathrm{obs,los}}$) and gas mass surface densities are larger. The GMC size -- linewidth relation ($\sigma_{\mathrm{obs,los}}\propto R_{\mathrm{c}}^{0.77}$) is steeper than that of the clouds of the MW disc and centre, and the GMCs are on average only marginally gravitationally bound (with a mean virial parameter $\langle\alpha_{\mathrm{obs,vir}}\rangle\approx1.7$). We discuss the possible origins of the steep size -- linewidth relation and enhanced observed linewidths of the clouds and suggest that a combination of mechanisms such as stellar feedback, gas accretion and cloud-cloud collisions, as well as the gas inflows driven by the large-scale bar, may play a role., Comment: 15 pages, 8 figures, accepted for publication in MNRAS. arXiv admin note: text overlap with arXiv:2304.10471

Understanding how interacting particles approach thermal equilibrium is a major challenge of quantum simulators. Unlocking the full potential of such systems toward this goal requires flexible initial state preparation, precise time evolution, and extensive probes for final state characterization. We present a quantum simulator comprising 69 superconducting qubits which supports both universal quantum gates and high-fidelity analog evolution, with performance beyond the reach of classical simulation in cross-entropy benchmarking experiments. Emulating a two-dimensional (2D) XY quantum magnet, we leverage a wide range of measurement techniques to study quantum states after ramps from an antiferromagnetic initial state. We observe signatures of the classical Kosterlitz-Thouless phase transition, as well as strong deviations from Kibble-Zurek scaling predictions attributed to the interplay between quantum and classical coarsening of the correlated domains. This interpretation is corroborated by injecting variable energy density into the initial state, which enables studying the effects of the eigenstate thermalization hypothesis (ETH) in targeted parts of the eigenspectrum. Finally, we digitally prepare the system in pairwise-entangled dimer states and image the transport of energy and vorticity during thermalization. These results establish the efficacy of superconducting analog-digital quantum processors for preparing states across many-body spectra and unveiling their thermalization dynamics.

Observations show a tight correlation between the stellar mass of galaxies and their gas-phase metallicity (MZR). This relation evolves with redshift, with higher-redshift galaxies being characterized by lower metallicities. Understanding the physical origin of the slope and redshift evolution of the MZR may provide important insight into the physical processes underpinning it: star formation, feedback, and cosmological inflows. While theoretical models ascribe the shape of the MZR to the lower efficiency of galactic outflows in more massive galaxies, what drives its evolution remains an open question. In this letter, we analyze how the MZR evolves over $z=0-3$, combining results from the FIREbox cosmological volume simulation with analytical models. Contrary to a frequent assertion in the literature, we find that the evolution of the gas fraction does not contribute significantly to the redshift evolution of the MZR. Instead, we show that the latter is driven by the redshift-dependence of the inflow metallicity, outflow metallicity, and mass loading factor, whose relative importance depends on stellar mass. These findings also suggest that the evolution of the MZR is not explained by galaxies moving along a fixed surface in the space spanned by stellar mass, gas phase metallicity, and star formation rate., Comment: 6 pages, 4 figures. Submitted for publication on MNRAS

A kinematic misalignment of the stellar and gas components is a phenomenon observed in a significant fraction of galaxies. However, the underlying physical mechanisms are not well understood. A commonly proposed scenario for the formation of a misaligned component requires any pre-existing gas disc to be removed, via fly-bys or ejective feedback from an active galactic nucleus. In this Letter, we study the evolution of a Milky Way mass galaxy in the FIREbox cosmological volume that displays a thin, counter-rotating gas disc with respect to its stellar component at low redshift. In contrast to scenarios involving gas ejection, we find that pre-existing gas is mainly removed via the conversion into stars in a central starburst, triggered by a merging satellite galaxy. The newly-accreted, counter-rotating gas eventually settles into a kinematically misaligned disc. About 4.4 (8 out of 182) of FIREbox galaxies with stellar masses larger than 5e9 Msun at z=0 exhibit gas-star kinematic misalignment. In all cases, we identify central starburst-driven depletion as the main reason for the removal of the pre-existing co-rotating gas component, with no need for feedback from, e.g., a central active black hole. However, during the starburst, the gas is funneled towards the central regions, likely enhancing black hole activity. By comparing the fraction of misaligned discs between FIREbox and other simulations and observations, we conclude that this channel might have a non-negligible role in inducing kinematic misalignment in galaxies., Comment: 9 pages; 3 figures; submitted to ApJ Letters

We use high-resolution data from the millimetre-Wave Interferometric Survey of Dark Object Masses (WISDOM) project to investigate the connection between circumnuclear gas reservoirs and nuclear activity in a sample of nearby galaxies. Our sample spans a wide range of nuclear activity types including radio galaxies, Seyfert galaxies, low-luminosity active galactic nuclei (AGN) and inactive galaxies. We use measurements of nuclear millimetre continuum emission along with other archival tracers of AGN accretion/activity to investigate previous claims that at, circumnuclear scales (<100 pc), these should correlate with the mass of the cold molecular gas. We find that the molecular gas mass does not correlate with any tracer of nuclear activity. This suggests the level of nuclear activity cannot solely be regulated by the amount of cold gas around the supermassive black hole (SMBH). This indicates that AGN fuelling, that drives gas from the large scale galaxy to the nuclear regions, is not a ubiquitous process and may vary between AGN type, with timescale variations likely to be very important. By studying the structure of the central molecular gas reservoirs, we find our galaxies have a range of nuclear molecular gas concentrations. This could indicate that some of our galaxies may have had their circumnuclear regions impacted by AGN feedback, even though they currently have low nuclear activity. On the other hand, the nuclear molecular gas concentrations in our galaxies could instead be set by secular processes., Comment: 15 pages plus 3 in the appendix, 8 figures plus 1 in the appendix, 3 tables plus 4 in the appendix

Observations indicate that the central gas discs are smoother in early-type galaxies than their late-type counterparts, while recent simulations predict that the dynamical suppression of star formation in spheroid-dominated galaxies is preceded by the suppression of fragmentation of their interstellar media. The mass surface density power spectrum is a powerful tool to constrain the degree of structure within a gas reservoir. Specifically here, we focus on the power spectrum slope and aim to constrain whether the shear induced by a dominant spheroidal potential can induce sufficient turbulence to suppress fragmentation, resulting in the smooth central gas discs observed. We compute surface density power spectra for the nuclear gas reservoirs of fourteen simulated isolated galaxies and twelve galaxies observed as part of the mm-Wave Interferometric Survey of Dark Object Masses (WISDOM) project. Both simulated and observed galaxies range from disc-dominated galaxies to spheroids, with central stellar mass surface densities, a measure of bulge dominance, varying by more than an order of magnitude. For the simulations, the power spectra steepen with increasing central stellar mass surface density, thereby clearly linking the suppression of fragmentation to the shear-driven turbulence induced by the spheroid. The WISDOM observations show a different (but potentially consistent) picture: while there is no correlation between the power spectrum slopes and the central stellar mass surface densities, the slopes scatter around a value of 2.6. This is similar to the behaviour of the slopes of the simulated galaxies with high central stellar mass surface densities, and could indicate that high shear eventually drives incompressible turbulence., Comment: 14 pages, 7 figures + appendices (7 pages, 9 figures); accepted by MNRAS (October 09th, 2023)

Understanding what shapes the cold gas component of galaxies, which both provides the fuel for star formation and is strongly affected by the subsequent stellar feedback, is a crucial step towards a better understanding of galaxy evolution. Here, we analyse the HI properties of a sample of 46 Milky Way halo-mass galaxies, drawn from cosmological simulations (EMP-Pathfinder and FIREbox). This set of simulations comprises galaxies evolved self-consistently across cosmic time with different baryonic sub-grid physics: three different star formation models [constant star formation efficiency (SFE) with different star formation eligibility criteria, and an environmentally-dependent, turbulence-based SFE] and two different feedback prescriptions, where only one sub-sample includes early stellar feedback. We use these simulations to assess the impact of different baryonic physics on the HI content of galaxies. We find that the galaxy-wide HI properties agree with each other and with observations. However, differences appear for small-scale properties. The thin HI discs observed in the local Universe are only reproduced with a turbulence-dependent SFE and/or early stellar feedback. Furthermore, we find that the morphology of HI discs is particularly sensitive to the different physics models: galaxies simulated with a turbulence-based SFE have discs that are smoother and more rotationally symmetric, compared to those simulated with a constant SFE; galaxies simulated with early stellar feedback have more regular discs than supernova-feedback-only galaxies. We find that the rotational asymmetry of the HI discs depends most strongly on the underlying physics model, making this a promising observable for understanding the physics responsible for shaping the interstellar medium of galaxies., Comment: 13 pages, 10 figures + appendices (7 pages, 10 figures); updated to version accepted by MNRAS

Starburst (SB) galaxies are a rare population of galaxies with star formation rates (SFRs) greatly exceeding those of the majority of star-forming galaxies with similar stellar mass. It is unclear whether these bursts are the result of either especially large gas reservoirs or enhanced efficiencies in converting gas into stars. Tidal torques resulting from gas-rich galaxy mergers are known to enhance the SFR by funneling gas towards the centre. However, recent theoretical works show that mergers do not always trigger a SB and not all SB galaxies are interacting systems, raising the question of what drives a SB. We analyse a large sample of SB galaxies and a mass- and redshift-matched sample of control galaxies, drawn from the FIREbox cosmological volume at z=0-1. We find that SB galaxies have both larger molecular gas fractions and shorter molecular depletion times than control galaxies, but similar total gas masses. Control galaxies evolve towards the SB regime by gas compaction in their central regions, over timescales of about 70 Myr, accompanied by an increase in the fraction of ultra-dense and molecular gas. The driving mechanism behind the SB varies depending on the mass of the galaxy. Massive (Mstar > 1e10 Msun) galaxies undergoing intense, long-lasting SBs are mostly driven by galaxy interactions. Conversely, SBs in non-interacting galaxies are often triggered by a global gravitational instability, that can result in a breathing mode in low-mass galaxies., Comment: 21 pages, 18 figures, submitted to MNRAS

We present a study of the molecular gas of seven early-type galaxies with high angular resolution data obtained as part of the mm-Wave Interferometric Survey of Dark Object Masses (WISDOM) project with the Atacama Large Millimeter/submillimeter Array. Using a fixed spatial scale approach, we study the mass surface density ($\Sigma$) and velocity dispersion ($\sigma$) of the molecular gas on spatial scales ranging from $60$ to $120$pc. Given the spatial resolution of our data ($20$ - $70$pc), we characterise these properties across many thousands of individual sight lines ($\approx50,000$ at our highest physical resolution). The molecular gas along these sight lines has a large range ($\approx2$dex) of mass surface densities and velocity dispersions $\approx40\%$ higher than those of star-forming spiral galaxies. It has virial parameters $\alpha_\mathrm{vir}$ that depend weakly on the physical scale observed, likely due to beam smearing of the bulk galactic rotation, and is generally super-virial. Comparing the internal turbulent pressure ($P_\mathrm{turb}$) to the pressure required for dynamic equilibrium ($P_\mathrm{DE}$), the ratio $P_\mathrm{turb}$/$P_\mathrm{DE}$ is significantly less than unity in all galaxies, indicating that the gas is not in dynamic equilibrium and is strongly compressed, in apparent contradiction to the virial parameters. This may be due to our neglect of shear and tidal forces, and/or the combination of three-dimensional and vertical diagnostics. Both $\alpha_\mathrm{vir}$ and $P_\mathrm{turb}$ anti-correlate with the global star-formation rate of our galaxies. We therefore conclude that the molecular gas in early-type galaxies is likely unbound, and that large-scale dynamics likely plays a critical role in its regulation. This contrasts to the giant molecular clouds in the discs of late-type galaxies, that are much closer to dynamical equilibrium., Comment: 32 pages (16 of Appendices), 39 Figures (27 in Appendices). Accepted for publication in MNRAS

We report the discovery of the ``mm fundamental plane of black-hole accretion'', which is a tight correlation between the nuclear 1 mm luminosity ($L_{\rm \nu, mm}$), the intrinsic $2$ -- $10$~keV X-ray luminosity ($L_{\rm X,2-10}$) and the supermassive black hole (SMBH) mass ($M_{\rm BH}$) with an intrinsic scatter ($\sigma_{\rm int}$) of $0.40$ dex. The plane is found for a sample of 48 nearby galaxies, most of which are low-luminosity active galactic nuclei (LLAGN). Combining these sources with a sample of high-luminosity (quasar-like) nearby AGN, we find that the plane still holds. We also find that $M_{\rm BH}$ correlates with $L_{\rm \nu, mm}$ at a highly significant level, although such correlation is less tight than the mm fundamental plane ($\sigma_{\rm int}=0.51$ dex). Crucially, we show that spectral energy distribution (SED) models for both advection-dominated accretion flows (ADAFs) and compact jets can explain the existence of these relations, which are not reproduced by the standard torus-thin accretion disc models usually associated to quasar-like AGN. The ADAF models reproduces the observed relations somewhat better than those for compact jets, although neither provides a perfect prediction. Our findings thus suggest that radiatively-inefficient accretion processes such as those in ADAFs or compact (and thus possibly young) jets may play a key role in both low- and high-luminosity AGN. This mm fundamental plane also offers a new, rapid method to (indirectly) estimate SMBH masses., Comment: 9 pages, 5 figures

Understanding universal aspects of quantum dynamics is an unresolved problem in statistical mechanics. In particular, the spin dynamics of the 1D Heisenberg model were conjectured to belong to the Kardar-Parisi-Zhang (KPZ) universality class based on the scaling of the infinite-temperature spin-spin correlation function. In a chain of 46 superconducting qubits, we study the probability distribution, $P(\mathcal{M})$, of the magnetization transferred across the chain's center. The first two moments of $P(\mathcal{M})$ show superdiffusive behavior, a hallmark of KPZ universality. However, the third and fourth moments rule out the KPZ conjecture and allow for evaluating other theories. Our results highlight the importance of studying higher moments in determining dynamic universality classes and provide key insights into universal behavior in quantum systems.

We present high spatial resolution ($\approx24$ pc) Atacama Large Millimeter/sub-millimeter Array $^{12}$CO(2-1) observations of the central region of the nearby barred spiral galaxy NGC 5806. NGC 5806 has a highly structured molecular gas distribution with a clear nucleus, a nuclear ring and offset dust lanes. We identify $170$ spatially- and spectrally-resolved giant molecular clouds (GMCs). These clouds have comparable sizes ($R_{\mathrm{c}}$) and larger gas masses, observed linewidths ($\sigma_{\mathrm{obs,los}}$) and gas mass surface densities than those of clouds in the Milky Way disc. The size -- linewidth relation of the clouds is one of the steepest reported so far ($\sigma_{\mathrm{obs,los}}\propto R_{\mathrm{c}}^{1.20}$), the clouds are on average only marginally bound (with a mean virial parameter $\langle\alpha_{\mathrm{vir}}\rangle\approx2$), and high velocity dispersions are observed in the nuclear ring. These behaviours are likely due to bar-driven gas shocks and inflows along the offset dust lanes, and we infer an inflow velocity of $\approx120$ kms$^{-1}$ and a total molecular gas mass inflow rate of $\approx5$ M$_\odot$ yr$^{-1}$ into the nuclear ring. The observed internal velocity gradients of the clouds are consistent with internal turbulence. The number of clouds in the nuclear ring decreases with azimuthal angle downstream from the dust lanes without clear variation of cloud properties. This is likely due to the estimated short lifetime of the clouds ($\approx6$ Myr), which appears to be mainly regulated by cloud-cloud collision and/or shear processes. Overall, it thus seems that the presence of the large-scale bar and gas inflows to the centre of NGC 5806 affect cloud properties., Comment: Accepted for publication in MNRAS, 20 pages, 16 figures

Measurement has a special role in quantum theory: by collapsing the wavefunction it can enable phenomena such as teleportation and thereby alter the "arrow of time" that constrains unitary evolution. When integrated in many-body dynamics, measurements can lead to emergent patterns of quantum information in space-time that go beyond established paradigms for characterizing phases, either in or out of equilibrium. On present-day NISQ processors, the experimental realization of this physics is challenging due to noise, hardware limitations, and the stochastic nature of quantum measurement. Here we address each of these experimental challenges and investigate measurement-induced quantum information phases on up to 70 superconducting qubits. By leveraging the interchangeability of space and time, we use a duality mapping, to avoid mid-circuit measurement and access different manifestations of the underlying phases -- from entanglement scaling to measurement-induced teleportation -- in a unified way. We obtain finite-size signatures of a phase transition with a decoding protocol that correlates the experimental measurement record with classical simulation data. The phases display sharply different sensitivity to noise, which we exploit to turn an inherent hardware limitation into a useful diagnostic. Our work demonstrates an approach to realize measurement-induced physics at scales that are at the limits of current NISQ processors.

We measure empirical relationships between the local star formation rate (SFR) and properties of the star-forming molecular gas on 1.5 kpc scales across 80 nearby galaxies. These relationships, commonly referred to as "star formation laws," aim at predicting the local SFR surface density from various combinations of molecular gas surface density, galactic orbital time, molecular cloud free-fall time, and the interstellar medium dynamical equilibrium pressure. Leveraging a multiwavelength database built for the PHANGS survey, we measure these quantities consistently across all galaxies and quantify systematic uncertainties stemming from choices of SFR calibrations and the CO-to-H$_2$ conversion factors. The star formation laws we examine show 0.3-0.4 dex of intrinsic scatter, among which the molecular Kennicutt-Schmidt relation shows a $\sim$10% larger scatter than the other three. The slope of this relation ranges $\beta\approx0.9{-}1.2$, implying that the molecular gas depletion time remains roughly constant across the environments probed in our sample. The other relations have shallower slopes ($\beta\approx0.6{-}1.0$), suggesting that the star formation efficiency (SFE) per orbital time, the SFE per free-fall time, and the pressure-to-SFR surface density ratio (i.e., the feedback yield) may vary systematically with local molecular gas and SFR surface densities. Last but not least, the shapes of the star formation laws depend sensitively on methodological choices. Different choices of SFR calibrations can introduce systematic uncertainties of at least 10-15% in the star formation law slopes and 0.15-0.25 dex in their normalization, while the CO-to-H$_2$ conversion factors can additionally produce uncertainties of 20-25% for the slope and 0.10-0.20 dex for the normalization., Comment: 10 pages main text + 2 appendices. ApJL in press. Data products available at https://www.canfar.net/storage/list/phangs/RELEASES/Sun_etal_2023 . Slides summarizing key results can be found at https://www.dropbox.com/s/5gsegexeo9n0t05/Sun_et_PHANGS_2023.pptx?dl=0

Leakage of quantum information out of computational states into higher energy states represents a major challenge in the pursuit of quantum error correction (QEC). In a QEC circuit, leakage builds over time and spreads through multi-qubit interactions. This leads to correlated errors that degrade the exponential suppression of logical error with scale, challenging the feasibility of QEC as a path towards fault-tolerant quantum computation. Here, we demonstrate the execution of a distance-3 surface code and distance-21 bit-flip code on a Sycamore quantum processor where leakage is removed from all qubits in each cycle. This shortens the lifetime of leakage and curtails its ability to spread and induce correlated errors. We report a ten-fold reduction in steady-state leakage population on the data qubits encoding the logical state and an average leakage population of less than $1 \times 10^{-3}$ throughout the entire device. The leakage removal process itself efficiently returns leakage population back to the computational basis, and adding it to a code circuit prevents leakage from inducing correlated error across cycles, restoring a fundamental assumption of QEC. With this demonstration that leakage can be contained, we resolve a key challenge for practical QEC at scale., Comment: Main text: 7 pages, 5 figures

Indistinguishability of particles is a fundamental principle of quantum mechanics. For all elementary and quasiparticles observed to date - including fermions, bosons, and Abelian anyons - this principle guarantees that the braiding of identical particles leaves the system unchanged. However, in two spatial dimensions, an intriguing possibility exists: braiding of non-Abelian anyons causes rotations in a space of topologically degenerate wavefunctions. Hence, it can change the observables of the system without violating the principle of indistinguishability. Despite the well developed mathematical description of non-Abelian anyons and numerous theoretical proposals, the experimental observation of their exchange statistics has remained elusive for decades. Controllable many-body quantum states generated on quantum processors offer another path for exploring these fundamental phenomena. While efforts on conventional solid-state platforms typically involve Hamiltonian dynamics of quasi-particles, superconducting quantum processors allow for directly manipulating the many-body wavefunction via unitary gates. Building on predictions that stabilizer codes can host projective non-Abelian Ising anyons, we implement a generalized stabilizer code and unitary protocol to create and braid them. This allows us to experimentally verify the fusion rules of the anyons and braid them to realize their statistics. We then study the prospect of employing the anyons for quantum computation and utilize braiding to create an entangled state of anyons encoding three logical qubits. Our work provides new insights about non-Abelian braiding and - through the future inclusion of error correction to achieve topological protection - could open a path toward fault-tolerant quantum computing.

We demonstrate a high dynamic range Josephson parametric amplifier (JPA) in which the active nonlinear element is implemented using an array of rf-SQUIDs. The device is matched to the 50 $\Omega$ environment with a Klopfenstein-taper impedance transformer and achieves a bandwidth of 250-300 MHz, with input saturation powers up to -95 dBm at 20 dB gain. A 54-qubit Sycamore processor was used to benchmark these devices, providing a calibration for readout power, an estimate of amplifier added noise, and a platform for comparison against standard impedance matched parametric amplifiers with a single dc-SQUID. We find that the high power rf-SQUID array design has no adverse effect on system noise, readout fidelity, or qubit dephasing, and we estimate an upper bound on amplifier added noise at 1.6 times the quantum limit. Lastly, amplifiers with this design show no degradation in readout fidelity due to gain compression, which can occur in multi-tone multiplexed readout with traditional JPAs., Comment: 10 pages, 10 figures

Practical quantum computing will require error rates that are well below what is achievable with physical qubits. Quantum error correction offers a path to algorithmically-relevant error rates by encoding logical qubits within many physical qubits, where increasing the number of physical qubits enhances protection against physical errors. However, introducing more qubits also increases the number of error sources, so the density of errors must be sufficiently low in order for logical performance to improve with increasing code size. Here, we report the measurement of logical qubit performance scaling across multiple code sizes, and demonstrate that our system of superconducting qubits has sufficient performance to overcome the additional errors from increasing qubit number. We find our distance-5 surface code logical qubit modestly outperforms an ensemble of distance-3 logical qubits on average, both in terms of logical error probability over 25 cycles and logical error per cycle ($2.914\%\pm 0.016\%$ compared to $3.028\%\pm 0.023\%$). To investigate damaging, low-probability error sources, we run a distance-25 repetition code and observe a $1.7\times10^{-6}$ logical error per round floor set by a single high-energy event ($1.6\times10^{-7}$ when excluding this event). We are able to accurately model our experiment, and from this model we can extract error budgets that highlight the biggest challenges for future systems. These results mark the first experimental demonstration where quantum error correction begins to improve performance with increasing qubit number, illuminating the path to reaching the logical error rates required for computation., Comment: Main text: 6 pages, 4 figures. v2: Update author list, references, Fig. S12, Table IV

Accurately reproducing the thin cold gas discs observed in nearby spiral galaxies has been a long standing issue in cosmological simulations. Here, we present measurements of the radially resolved HI scale height in 22 non-interacting Milky Way-mass galaxies from the FIREbox cosmological volume. We measure the HI scale heights using five different approaches commonly used in the literature: fitting the vertical volume density distribution with a Gaussian, the distance between maximum and half-maximum of the vertical volume density distribution, a semi-empirical description using the velocity dispersion and the galactic gravitational potential, the analytic assumption of hydrostatic equilibrium, and the distance from the midplane which encloses $\gtrsim$60 per cent of the HI mass. We find median HI scale heights, measured using the vertical volume distribution, that range from ~100 pc in the galactic centres to ~800 pc in the outskirts and are in excellent agreement with recent observational results. We speculate that the presence of a realistic multiphase interstellar medium, including cold gas, and realistic stellar feedback are the drivers behind the realistic HI scale heights., Comment: 5 pages + appendix, 3 figures, updated to version accepted by MNRAS

Systems of correlated particles appear in many fields of science and represent some of the most intractable puzzles in nature. The computational challenge in these systems arises when interactions become comparable to other energy scales, which makes the state of each particle depend on all other particles. The lack of general solutions for the 3-body problem and acceptable theory for strongly correlated electrons shows that our understanding of correlated systems fades when the particle number or the interaction strength increases. One of the hallmarks of interacting systems is the formation of multi-particle bound states. In a ring of 24 superconducting qubits, we develop a high fidelity parameterizable fSim gate that we use to implement the periodic quantum circuit of the spin-1/2 XXZ model, an archetypal model of interaction. By placing microwave photons in adjacent qubit sites, we study the propagation of these excitations and observe their bound nature for up to 5 photons. We devise a phase sensitive method for constructing the few-body spectrum of the bound states and extract their pseudo-charge by introducing a synthetic flux. By introducing interactions between the ring and additional qubits, we observe an unexpected resilience of the bound states to integrability breaking. This finding goes against the common wisdom that bound states in non-integrable systems are unstable when their energies overlap with the continuum spectrum. Our work provides experimental evidence for bound states of interacting photons and discovers their stability beyond the integrability limit., Comment: 7 pages + 15 pages supplements

The star formation efficiency (SFE) has been shown to vary across different environments, particularly within galactic starbursts and deep within the bulges of galaxies. Various quenching mechanisms may be responsible, ranging from galactic dynamics to feedback from active galactic nuclei (AGN). Here, we use spatially-resolved observations of warm ionised gas emission lines from the imaging Fourier transform spectrograph SITELLE at the Canada-France-Hawaii Telescope (CFHT) and cold molecular gas (CO(2-1)) from the Atacama Large Millimeter/sub-millimeter Array (ALMA) to study the SFE in the bulge of the AGN-host galaxy NGC 3169. After distinguishing star-forming regions from AGN-ionised regions using emission-line ratio diagnostics, we measure spatially-resolved molecular gas depletion times (\tau_dep = 1/SFE) with a spatial resolution of \approx 100 pc within a galactocentric radius of 1.8 kpc. We identify a star-forming ring located at radii 1.25 \pm 0.6 kpc with an average \tau_dep of 0.3 Gyr. At radii < 0.9 kpc, however, the molecular gas surface densities and depletion times increase with decreasing radius, the latter reaching approximately 2.3 Gyr at a radius \approx 500 pc. Based on analyses of the gas kinematics and comparisons with simulations, we identify AGN feedback, bulge morphology and dynamics as the possible causes of the radial profile of SFE observed in the central region of NGC 3169., Comment: 25 pages, 15 figures, accepted for publication in MNRAS

We introduce a suite of cosmological volume simulations to study the evolution of galaxies as part of the Feedback in Realistic Environments project. FIREbox, the principal simulation of the present suite, provides a representative sample of galaxies (~1000 galaxies with Mstar > 10^8 Msun at z=0) at a resolution (~20 pc, m_b ~ 6x10^4 Msun) comparable to state-of-the-art galaxy zoom-in simulations. FIREbox captures the multiphase nature of the interstellar medium in a fully cosmological setting (L=22.1 Mpc) thanks to its exceptionally high dynamic range (~10^6) and the inclusion of multi-channel stellar feedback. Here, we focus on validating the simulation predictions by comparing to observational data. We find that simulated galaxies with Mstar < 10^{10.5-11} Msun have star formation rates, gas masses, and metallicities in broad agreement with observations. These galaxy scaling relations extend to low masses (Mstar ~ 10^7 Msun) and follow a (broken) power-law relationship. Also reproduced are the evolution of the cosmic HI density and the HI column density distribution at z~0-5. At low z, FIREbox predicts a peak in the stellar-mass--halo-mass relation, but also a higher abundance of massive galaxies and a higher cosmic star formation rate density than observed, showing that stellar feedback alone is insufficient to reproduce the properties of massive galaxies at late times. Given its high resolution and sample size, FIREbox offers a baseline prediction of galaxy formation theory in a $\Lambda$CDM Universe while also highlighting modeling challenges to be addressed in next-generation galaxy simulations., Comment: 31 pages, 17 figures, 3 tables, accepted for publication in MNRAS

We use high-resolution maps of the molecular interstellar medium (ISM) in the centres of eighty-six nearby galaxies from the millimetre-Wave Interferometric Survey of Dark Object Masses (WISDOM) and Physics at High Angular Resolution in Nearby GalaxieS (PHANGS) surveys to investigate the physical mechanisms setting the morphology of the ISM at molecular cloud scales. We show that early-type galaxies tend to have smooth, regular molecular gas morphologies, while the ISM in spiral galaxy bulges is much more asymmetric and clumpy when observed at the same spatial scales. We quantify these differences using non-parametric morphology measures (Asymmetry, Smoothness and Gini), and compare these measurements with those extracted from idealised galaxy simulations. We show that the morphology of the molecular ISM changes systematically as a function of various large-scale galaxy parameters, including galaxy morphological type, stellar mass, stellar velocity dispersion, effective stellar mass surface density, molecular gas surface density, star formation efficiency and the presence of a bar. We perform a statistical analysis to determine which of these correlated parameters best predicts the morphology of the ISM. We find the effective stellar mass surface (or volume) density to be the strongest predictor of the morphology of the molecular gas, while star formation and bars maybe be important secondary drivers. We find that gas self-gravity is not the dominant process shaping the morphology of the molecular gas in galaxy centres. Instead effects caused by the depth of the potential well such as shear, suppression of stellar spiral density waves and/or inflow affect the ability of the gas to fragment., Comment: 20 pages, 9 figures. Accepted for publication in MNRAS

The formation and evolution of stellar clusters is intimately linked to that of their host galaxies. To study this connection, we present the EMP-Pathfinder suite of cosmological zoom-in Milky Way-mass simulations. These simulations contain a sub-grid description for stellar cluster formation and evolution, allowing us to study the simultaneous formation and evolution of stellar clusters alongside their host galaxies across cosmic time. As a key ingredient in these simulations, we include the physics of the multi-phase nature of the interstellar medium (ISM), which enables studies of how the presence of a cold, dense ISM affects cluster formation and evolution. We consider two different star formation prescriptions: a constant star formation efficiency per free-fall time, as well as an environmentally-dependent, turbulence-based prescription. We identify two key results drawn from these simulations. Firstly, we find that tidal shock-driven disruption caused by the graininess of the cold ISM produces old ($\tau>10~$Gyr) stellar cluster populations with properties that are in excellent agreement with the observed populations in the Milky Way and M31. Importantly, the addition of the cold ISM addresses the areas of disagreement found in previous simulations that lacked the cold gas phase. Secondly, the formation of stellar clusters is extremely sensitive to the baryonic physics that govern the properties of the cold, dense gas reservoir in the galaxy. This implies that the demographics of stellar cluster populations represent an important diagnostic tool for constraining baryonic physics models in upcoming galaxy formation simulations that also include a description of the cold ISM., Comment: 30 pages + appendices, 12 figures. MNRAS submitted, comments and suggestions welcome!

This paper draws upon the evolutionary concepts of technological relatedness and knowledge complexity to enhance our understanding of the long-term evolution of Artificial Intelligence (AI). We reveal corresponding patterns in the emergence of AI - globally and in the context of specific geographies of the US, Japan, South Korea, and China. We argue that AI emergence is associated with increasing related variety due to knowledge commonalities as well as increasing complexity. We use patent-based indicators for the period between 1974-2018 to analyse the evolution of AI's global technological space, to identify its technological core as well as changes to its overall relatedness and knowledge complexity. At the national level, we also measure countries' overall specialisations against AI-specific ones. At the global level, we find increasing overall relatedness and complexity of AI. However, for the technological core of AI, which has been stable over time, we find decreasing related variety and increasing complexity. This evidence points out that AI innovations related to core technologies are becoming increasingly distinct from each other. At the country level, we find that the US and Japan have been increasing the overall relatedness of their innovations. The opposite is the case for China and South Korea, which we associate with the fact that these countries are overall less technologically developed than the US and Japan. Finally, we observe a stable increasing overall complexity for all countries apart from China, which we explain by the focus of this country in technologies not strongly linked to AI., Comment: Keywords: Artificial Intelligence; technological space; evolutionary economic geography; technological relatedness; knowledge complexity

In simple models of galaxy formation and evolution, star formation is solely regulated by the amount of gas present in the galaxy. However, it has recently been shown that star formation can be suppressed by galactic dynamics in galaxies that contain a dominant spheroidal component and a low gas fraction. This 'dynamical suppression' is hypothesised to also contribute to quenching gas-rich galaxies at high redshift, but its impact on the galaxy population at large remains unclear. In this paper, we assess the importance of dynamical suppression in the context of gas regulator models of galaxy evolution through hydrodynamic simulations of isolated galaxies, with gas-to-stellar mass ratios of 0.01-0.20 and a range of galactic gravitational potentials from disc-dominated to spheroidal. Star formation is modelled using a dynamics-dependent efficiency per free-fall time, which depends on the virial parameter of the gas. We find that dynamical suppression becomes more effective at lower gas fractions and quantify its impact on the star formation rate as a function of gas fraction and stellar spheroid mass surface density. We combine the results of our simulations with observed scaling relations that describe the change of galaxy properties across cosmic time, and determine the galaxy mass and redshift range where dynamical suppression may affect the baryon cycle. We predict that the physics of star formation can limit and regulate the baryon cycle at low redshifts ($z \lesssim 1.4$) and high galaxy masses ($M_{\ast} \gtrsim 3 \times 10^{10}~M_{\odot}$), where dynamical suppression can drive galaxies off the star formation main sequence., Comment: 12 pages, 6 figures, accepted by MNRAS (Nov. 03, 2020)

Quenched galaxies are often observed to contain a strong bulge component. The key question is whether this reflects a causal connection - can star formation be quenched dynamically by bulges or the spheroids of early-type galaxies? We systematically investigate the impact of these morphological components on star formation, by performing a suite of hydrodynamical simulations of isolated galaxies containing a spheroid. We vary the bulge mass and scale radius, while the total initial stellar, halo and gas mass are kept constant, with a gas fraction of 5 per cent. In addition, we consider two different sub-grid star formation prescriptions. The first follows most simulations in the literature by assuming a constant star formation efficiency per free-fall time, whereas in the second model it depends on the gas virial parameter, following high-resolution simulations of turbulent fragmentation. Across all simulations, central spheroids increase the gas velocity dispersion towards the galactic centre. This increases the gravitational stability of the gas disc, suppresses fragmentation and star formation, and results in galaxies hosting extremely smooth and quiescent gas discs that fall below the galaxy main sequence. These effects amplify when using the more sophisticated, dynamics-dependent star formation model. Finally, we discover a pronounced relation between the central stellar surface density and star formation rate (SFR), such that the most bulge-dominated galaxies show the strongest deviation from the main sequence. We conclude that the SFR of galaxies is not only set by the balance between accretion and feedback, but carries a (sometimes dominant) dependence on the gravitational potential., Comment: 25 pages, 22 figures, 2 tables; updated to version accepted by MNRAS (April 23, 2020)