We present Denario, an AI multi-agent system designed to serve as a scientific research assistant. Denario can perform many different tasks, such as generating ideas, checking the literature, developing research plans, writing and executing code, making plots, and drafting and reviewing a scientific paper. The system has a modular architecture, allowing it to handle specific tasks, such as generating an idea, or carrying out end-to-end scientific analysis using Cmbagent as a deep-research backend. In this work, we describe in detail Denario and its modules, and illustrate its capabilities by presenting multiple AI-generated papers generated by it in many different scientific disciplines such as astrophysics, biology, biophysics, biomedical informatics, chemistry, material science, mathematical physics, medicine, neuroscience and planetary science. Denario also excels at combining ideas from different disciplines, and we illustrate this by showing a paper that applies methods from quantum physics and machine learning to astrophysical data. We report the evaluations performed on these papers by domain experts, who provided both numerical scores and review-like feedback. We then highlight the strengths, weaknesses, and limitations of the current system. Finally, we discuss the ethical implications of AI-driven research and reflect on how such technology relates to the philosophy of science. We publicly release the code at this https URL. A Denario demo can also be run directly on the web at this https URL, and the full app will be deployed on the cloud.

Massive black holes (MBHs) are typically hosted in the centres of massive galaxies but they appear to become rarer in lower mass galaxies, where nuclear star clusters (NSCs) frequently appear instead. The transition region, where both an MBH and NSC can co-exist, has been poorly studied to date and only a few dozen galaxies are known to host them. One avenue for detecting new galaxies with both an MBH and NSC is to look for accretion signatures of MBHs. Here, we use new SRG/eROSITA all-sky survey eRASS:4 data to search for X-ray signatures of accreting MBHs in NSCs, while also investigating their combined occupation fraction. We find significant detections for 18 galaxies (~8.3%), including one ultra-luminous X-ray source; however, only three galaxies (NGC2903, 4212, and 4639) have X-ray luminosities that are higher than the expected value from X-ray binaries, indicative of the presence of an MBH. In addition, the X-ray luminosity of six galaxies (NGC2903, 3384, 4321, 4365, 4639, and 4701) differs from previous studies and could indicate the presence of a variable active galactic nucleus. The combined occupation fraction of accreting MBHs and NSCs becomes non-zero for galaxy masses above ~10^7.5 M_sun and this result is slightly elevated as compared to the literature data. Our data extend, for the first time, towards the dwarf elliptical galaxy regime and identify promising MBH candidates for higher resolution follow-up observations. At most galaxy masses (and with the exception of three cases), the X-ray constraints are consistent with the expected emission from binary systems or an Eddington fraction of at most 0.01%, assuming a black holes mass of 10^6.5 M_sun. This work confirms the known complexities in similar-type of studies, while providing the appealing alternative of using X-ray survey data of in-depth observations of individual targets with higher resolution instruments.

Recent experiments observed fractional Chern insulators (FCI) in twisted bilayer MoTe$_2$ at zero magnetic field, yet even the single-particle model of this material is controversial, leading to unreliable predictions of the experimental phase diagram as discussed in [Yu et al., 2023]. In this light, we revisit the single-particle model of twisted bilayer MoTe$_2$. Utilizing large-scale density functional theory, we calculate the band structure of twisted AA-stacked bilayer MoTe$_2$ at various twist angles relevant to experiment. We find that a band inversion occurs near $4.41^\circ$ between the second and third bands. Our ab initio band structure is in qualitative agreement with [Wang et al., 2023], but shows important differences in the remote bands and in the $\Gamma$ valley. We incorporate two higher harmonic terms into the continuum model to capture the highest 3 valence bands per valley. We confirm that the two highest valence bands per valley have opposite Chern numbers with $|C|=1$ for small angles, and also use our model to predict a variety of Chern states in the remote bands accessible by displacement field. Finally, we perform DFT calculations and build models for the AB stacking configuration. Our work serves as a foundation for accurate determination of the correlated phases in twisted bilayer MoTe$_2$.

Fractionally filled Chern bands with strong interactions may give rise to fractional Chern insulator (FCI) states, the zero-field analogue of the fractional quantum Hall effect. Recent experiments have demonstrated the existence of FCIs in twisted bilayer MoTe$_2$ without external magnetic fields -- most robust at $\nu=-2/3$ -- as well as Chern insulators (CIs) at $\nu=-1$. Although the appearance of both of these states is theoretically natural in an interacting topological system, experiments repeatedly observe nonmagnetic states (lacking FCIs) at $\nu=-1/3$ and $-4/3$, a puzzling result which has not been fully theoretically explained. In this work, we perform Hartree-Fock and exact diagonalization calculations to test whether the standard MoTe$_2$ moiré model with the (greatly varying) parameter values available in the literature can reproduce the non-magnetic states at $\nu=-1/3$ and $-4/3$ in unison with the FCI at $\nu=-2/3$ and CI state at $\nu = -1$. We focus on the experimentally relevant twist angles and, crucially, include remote bands. We find that the parameters proposed in [Wang et al. (2023)] can nearly capture the experimental phenomena at $\nu=-1/3,-2/3,-1,-4/3$ simultaneously, though the predicted ground states at $\nu=-1/3$ are still mostly fully-spin-polarized and a larger dielectric constant $\epsilon>10$ than is typical of hexagonal boron nitride (h-BN) substrate $\epsilon\sim 6$ is required. Our results show the importance of remote bands in identifying the competing magnetic orders and lay the groundwork for further study of the realistic phase diagram.

In this work, using the time-dependent density functional theory, we address the electron tunneling triggered by short (single-cycle and several-cycle) optical pulses in narrow metallic gaps under conditions relevant for actual experiments. We identify photon-assisted tunneling with one-photon, two-photon, and higher-order photon absorption, and we discuss the effect of the tunneling barrier, applied bias, and strength of the optical field on transition from photon-assisted tunneling (weak optical fields) to the optical field emission at strong optical fields. The numerical single-electron calculations and an analytical strong-field theory model are used to gain deeper insights into the results of the time-dependent density functional theory calculations. Additionally, our parameter-free calculations allow us to retrieve and explain recent experimental results on optically induced transport in narrow metallic gaps.

Quantum technologies have the potential to solve certain computationally hard problems with polynomial or super-polynomial speedups when compared to classical methods. Unfortunately, the unstable nature of quantum information makes it prone to errors. For this reason, quantum error correction is an invaluable tool to make quantum information reliable and enable the ultimate goal of fault-tolerant quantum computing. Surface codes currently stand as the most promising candidates to build near term error corrected qubits given their two-dimensional architecture, the requirement of only local operations, and high tolerance to quantum noise. Decoding algorithms are an integral component of any error correction scheme, as they are tasked with producing accurate estimates of the errors that affect quantum information, so that they can subsequently be corrected. A critical aspect of decoding algorithms is their speed, since the quantum state will suffer additional errors with the passage of time. This poses a connundrum, where decoding performance is improved at the expense of complexity and viceversa. In this review, a thorough discussion of state-of-the-art decoding algorithms for surface codes is provided. The target audience of this work are both readers with an introductory understanding of the field as well as those seeking to further their knowledge of the decoding paradigm of surface codes. We describe the core principles of these decoding methods as well as existing variants that show promise for improved results. In addition, both the decoding performance, in terms of error correction capability, and decoding complexity, are compared. A review of the existing software tools regarding surface codes decoding is also provided.

Understanding the confinement mechanism in gauge theories and the universality of effective string-like descriptions of gauge flux tubes remains a fundamental challenge in modern physics. We probe string modes of motion with dynamical matter in a digital quantum simulation of a (2+1) dimensional gauge theory using a superconducting quantum processor with up to 144 qubits, stretching the hardware capabilities with quantum-circuit depths comprising up to 192 two-qubit layers. We realize the $Z_2$-Higgs model ($Z_2$HM) through an optimized embedding into a heavy-hex superconducting qubit architecture, directly mapping matter and gauge fields to vertex and link superconducting qubits, respectively. Using the structure of local gauge symmetries, we implement a comprehensive suite of error suppression, mitigation, and correction strategies to enable real-time observation and manipulation of electric strings connecting dynamical charges. Our results resolve a dynamical hierarchy of longitudinal oscillations and transverse bending at the end points of the string, which are precursors to hadronization and rotational spectra of mesons. We further explore multi-string processes, observing the fragmentation and recombination of strings. The experimental design supports 300,000 measurement shots per circuit, totaling 600,000 shots per time step, enabling high-fidelity statistics. We employ extensive tensor network simulations using the basis update and Galerkin method to predict large-scale real-time dynamics and validate our error-aware protocols. This work establishes a milestone for probing non-perturbative gauge dynamics via superconducting quantum simulation and elucidates the real-time behavior of confining strings.

The dynamics of quantum systems unfolds within a subspace of the state space or operator space, known as the Krylov space. This review presents the use of Krylov subspace methods to provide an efficient description of quantum evolution and quantum chaos, with emphasis on nonequilibrium phenomena of many-body systems with a large Hilbert space. It provides a comprehensive update of recent developments, focused on the quantum evolution of operators in the Heisenberg picture as well as pure and mixed states. It further explores the notion of Krylov complexity and associated metrics as tools for quantifying operator growth, their bounds by generalized quantum speed limits, the universal operator growth hypothesis, and its relation to quantum chaos, scrambling, and generalized coherent states. A comparison of several generalizations of the Krylov construction for open quantum systems is presented. A closing discussion addresses the application of Krylov subspace methods in quantum field theory, holography, integrability, quantum control, and quantum computing, as well as current open problems.

We show the presence of analytic $p_x + i p_y$ superconducting ground states in the Berry Trashcan -- a minimal model of rhombohedral graphene valid for $n \ge 4$ layers -- under short-range attractive interactions. We demonstrate that the model, whose dispersion consists of a flat bottom surrounded by steep walls of prohibitive kinetic energy, serves as a building block to understand superconductivity in the moiré-free limit. We find that the ground-state chirality has a ``ferromagnetic'' coupling to that of the uniform Berry curvature of the model, and compare the analytically obtained binding energies, excitation spectra and off-diagonal long-range order (ODLRO) with numerical exact diagonalization results. We show that the analytic structure of this model is that of a restricted spectrum generating algebra (RSGA), initially developed for quantum scars, and build a variety of other exact (but contrived) models with exact chiral superconductivity based on a method developed in Ref.[1]. However, under short range attraction, we show that the Berry Trashcan is the optimal and only realistic point in the class of GMP-like algebras to host a chiral superconductor state. A toy model in 1D and its related physics is also investigated. Our results reveal that chiral superconductivity is natural under attractive interactions in the Berry trashcan model of rhombohedral graphene in displacement field, although we make no claim about the origin of the attraction.

A large amount of effort has recently been put into understanding the barren plateau phenomenon. In this perspective article, we face the increasingly loud elephant in the room and ask a question that has been hinted at by many but not explicitly addressed: Can the structure that allows one to avoid barren plateaus also be leveraged to efficiently simulate the loss classically? We collect evidence-on a case-by-case basis-that many commonly used models whose loss landscapes avoid barren plateaus can also admit classical simulation, provided that one can collect some classical data from quantum devices during an initial data acquisition phase. This follows from the observation that barren plateaus result from a curse of dimensionality, and that current approaches for solving them end up encoding the problem into some small, classically simulable, subspaces. Thus, while stressing that quantum computers can be essential for collecting data, our analysis sheds doubt on the information processing capabilities of many parametrized quantum circuits with provably barren plateau-free landscapes. We end by discussing the (many) caveats in our arguments including the limitations of average case arguments, the role of smart initializations, models that fall outside our assumptions, the potential for provably superpolynomial advantages and the possibility that, once larger devices become available, parametrized quantum circuits could heuristically outperform our analytic expectations.

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

Chiral crystals offer an unique platform for controlling structural handedness through external stimuli. However, the ability to select between structural enantiomers remains challenging, both theoretically and experimentally. In this work, we demonstrate a two-step pathway for enantiomer selectivity in layered chiral NbOX$_2$ (X = Cl, Br, I) crystals based on photostriction-driven phase transitions. Ab-initio simulations reveal that optical excitation is capable of inducing a structural phase transition in NbOX$_2$ from the monoclinic ($C2$) ground state to the higher-symmetry ($C2/m$) structure. In the resulting transient high-symmetry state, an applied electric field breaks the residual inversion-symmetry degeneracy, selectively stabilizing one enantiomeric final state configuration over the other. Our results establish a combined optical-electrical control scheme for chiral materials, enabling reversible and non-contact enantiomer selection with potential applications in ultrafast switching, optoelectronics, and chiral information storage.

Intrinsic noise in pre-fault-tolerant quantum devices poses a major challenge to the reliable realization of unitary dynamics in quantum algorithms and simulations. To address this, we present a method for simulating open quantum system dynamics on a quantum computer, including negative dissipation rates in the Gorini-Kossakowski-Sudarshan-Lindblad (GKSL) master equation. Our approach lies beyond the standard Markovian approximation, enabling the controlled study of non-Markovian processes within a quantum simulation framework. Using this method, we develop a quantum algorithm for calculating ground-state properties that exploits feedback-controlled, noise-assisted dynamics. In this scheme, Lyapunov-based feedback steers the system toward a target virtual state under engineered noise conditions. This framework offers a promising strategy for harnessing current quantum hardware and advancing robust control protocols based on open system dynamics.