Source author record

Alessandro Roggero

Alessandro Roggero appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

20works
11topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

20 published item(s)

preprint2022arXiv

Entanglement and correlations in fast collective neutrino flavor oscillations

Collective neutrino oscillations play a crucial role in transporting lepton flavor in astrophysical settings like supernovae and neutron star binary merger remnants, which are characterized by large neutrino densities. In these settings, simulations in the mean-field approximation show that neutrino-neutrino interactions can overtake vacuum oscillations and give rise to fast collective flavor evolution on time-scales $t\proptoμ^{-1}$, with $μ$ proportional to the local neutrino density. In this work, we study the full out-of-equilibrium flavor dynamics in simple multi-angle geometries displaying fast oscillations in the mean field linear stability analysis. Focusing on simple initial conditions, we analyze the production of pair correlations and entanglement in the complete many-body-dynamics as a function of the number $N$ of neutrinos in the system, for up to thousands of neutrinos. Similarly to simpler geometries with only two neutrino beams, we identify three regimes: stable configurations with vanishing flavor oscillations, marginally unstable configurations with evolution occurring on long time scales $τ\approxμ^{-1}\sqrt{N}$, and unstable configurations showing flavor evolution on short time scales $τ\approxμ^{-1}\log(N)$. We present evidence that these fast collective modes are generated by the same dynamical phase transition which leads to the slow bipolar oscillations, establishing a connection between these two phenomena and explaining the difference in their time scales. We conclude by discussing a semi-classical approximation which reproduces the entanglement entropy at short to medium time scales and can be potentially useful in situations with more complicated geometries where classical simulation methods starts to become inefficient.

preprint2022arXiv

Hybridized Methods for Quantum Simulation in the Interaction Picture

Conventional methods of quantum simulation involve trade-offs that limit their applicability to specific contexts where their use is optimal. In particular, the interaction picture simulation has been found to provide substantial asymptotic advantages for some Hamiltonians, but incurs prohibitive constant factors and is incompatible with methods like qubitization. We provide a framework that allows different simulation methods to be hybridized and thereby improve performance for interaction picture simulations over known algorithms. These approaches show asymptotic improvements over the individual methods that comprise them and further make interaction picture simulation methods practical in the near term. Physical applications of these hybridized methods yield a gate complexity scaling as $\log^2 Λ$ in the electric cutoff $Λ$ for the Schwinger Model and independent of the electron density for collective neutrino oscillations, outperforming the scaling for all current algorithms with these parameters. For the general problem of Hamiltonian simulation subject to dynamical constraints, these methods yield a query complexity independent of the penalty parameter $λ$ used to impose an energy cost on time-evolution into an unphysical subspace.

preprint2022arXiv

Quantum Machine Learning with SQUID

In this work we present the Scaled QUantum IDentifier (SQUID), an open-source framework for exploring hybrid Quantum-Classical algorithms for classification problems. The classical infrastructure is based on PyTorch and we provide a standardized design to implement a variety of quantum models with the capability of back-propagation for efficient training. We present the structure of our framework and provide examples of using SQUID in a standard binary classification problem from the popular MNIST dataset. In particular, we highlight the implications for scalability for gradient-based optimization of quantum models on the choice of output for variational quantum models.

preprint2021arXiv

Dynamical Phase Transitions in models of Collective Neutrino Oscillations

Collective neutrino oscillations can potentially play an important role in transporting lepton flavor in astrophysical scenarios where the neutrino density is large, typical examples are the early universe and supernova explosions. It has been argued in the past that simple models of the neutrino Hamiltonian designed to describe forward scattering can support substantial flavor evolution on very short time scales $t\approx\log(N)/(G_Fρ_ν)$, with $N$ the number of neutrinos, $G_F$ the Fermi constant and $ρ_ν$ the neutrino density. This finding is in tension with results for similar but exactly solvable models for which $t\approx\sqrt{N}/(G_Fρ_ν)$ instead. In this work we provide a coherent explanation of this tension in terms of Dynamical Phase Transitions (DPT) and study the possible impact that a DPT could have in more realistic models of neutrino oscillations and their mean-field approximation.

preprint2021arXiv

Exact representations of many body interactions with RBM neural networks

Restricted Boltzmann Machines (RBM) are simple statistical models defined on a bipartite graph which have been successfully used in studying more complicated many-body systems, both classical and quantum. In this work, we exploit the representation power of RBMs to provide an exact decomposition of many-body contact interactions into one-body operators coupled to discrete auxiliary fields. This construction generalizes the well known Hirsch's transform used for the Hubbard model to more complicated theories such as Pionless EFT in nuclear physics, which we analyze in detail. We also discuss possible applications of our mapping for quantum annealing applications and conclude with some implications for RBM parameter optimization through machine learning.

preprint2021arXiv

Imaginary Time Propagation on a Quantum Chip

Evolution in imaginary time is a prominent technique for finding the ground state of quantum many-body systems, and the heart of a number of numerical methods that have been used with great success in quantum chemistry, condensed matter and nuclear physics. We propose an algorithm to implement imaginary time propagation on a quantum computer. Our algorithm is devised in the context of an efficient encoding into an optimized gate, drawing on the underlying characteristics of the quantum device, of a unitary operation in an extended Hilbert space. However, we proved that for simple problems it can be successfully applied to standard digital quantum machines. This work paves the way for porting quantum many-body methods based on imaginary-time propagation to near-term quantum devices, enabling the future quantum simulation of the ground states of a broad class of microscopic systems.

preprint2021arXiv

Nuclear two point correlation functions on a quantum-computer

The calculation of dynamic response functions is expected to be an early application benefiting from rapidly developing quantum hardware resources. The ability to calculate real-time quantities of strongly-correlated quantum systems is one of the most exciting applications that can easily reach beyond the capabilities of traditional classical hardware. Response functions of fermionic systems at moderate momenta and energies corresponding roughly to the Fermi energy of the system are a potential early application because the relevant operators are nearly local and the energies can be resolved in moderately short real time, reducing the spatial resolution and gate depth required. This is particularly the case in quasielastic electron and neutrino scattering from nuclei, a topic of great interest in the nuclear and particle physics communities and directly related to experiments designed to probe neutrino properties. In this work we use current quantum hardware and error mitigation protocols to calculate response functions for a highly simplified nuclear model through calculations of a 2-point real time correlation function for a modified Fermi-Hubbard model in two dimensions with three distinguishable nucleons on four lattice sites.

preprint2021arXiv

Preparation of excited states for nuclear dynamics on a quantum computer

We study two different methods to prepare excited states on a quantum computer, a key initial step to study dynamics within linear response theory. The first method uses unitary evolution for a short time $T=\mathcal{O}(\sqrt{1-F})$ to approximate the action of an excitation operator $\hat{O}$ with fidelity $F$ and success probability $P\approx1-F$. The second method probabilistically applies the excitation operator using the Linear Combination of Unitaries (LCU) algorithm. We benchmark these techniques on emulated and real quantum devices, using a toy model for thermal neutron-proton capture. Despite its larger memory footprint, the LCU-based method is efficient even on current generation noisy devices and can be implemented at a lower gate cost than a naive analysis would suggest. These findings show that quantum techniques designed to achieve good asymptotic scaling on fault tolerant quantum devices might also provide practical benefits on devices with limited connectivity and gate fidelity.

preprint2021arXiv

Spectral density reconstruction with Chebyshev polynomials

Accurate calculations of the spectral density in a strongly correlated quantum many-body system are of fundamental importance to study its dynamics in the linear response regime. Typical examples are the calculation of inclusive and semi-exclusive scattering cross sections in atomic nuclei and transport properties of nuclear and neutron star matter. Integral transform techniques play an important role in accessing the spectral density in a variety of nuclear systems. However, their accuracy is in practice limited by the need to perform a numerical inversion which is often ill-conditioned. In the present work we extend a recently proposed quantum algorithm which circumvents this problem. We show how to perform controllable reconstructions of the spectral density over a finite energy resolution with rigorous error estimates. An appropriate expansion in Chebyshev polynomials allows for efficient simulations also on classical computers. We apply our idea to reconstruct a simple model -- response function as a proof of principle. This paves the way for future applications in nuclear and condensed matter physics.

preprint2021arXiv

Standard Model Physics and the Digital Quantum Revolution: Thoughts about the Interface

Advances in isolating, controlling and entangling quantum systems are transforming what was once a curious feature of quantum mechanics into a vehicle for disruptive scientific and technological progress. Pursuing the vision articulated by Feynman, a concerted effort across many areas of research and development is introducing prototypical digital quantum devices into the computing ecosystem available to domain scientists. Through interactions with these early quantum devices, the abstract vision of exploring classically-intractable quantum systems is evolving toward becoming a tangible reality. Beyond catalyzing these technological advances, entanglement is enabling parallel progress as a diagnostic for quantum correlations and as an organizational tool, both guiding improved understanding of quantum many-body systems and quantum field theories defining and emerging from the Standard Model. From the perspective of three domain science theorists, this article compiles thoughts about the interface on entanglement, complexity, and quantum simulation in an effort to contextualize recent NISQ-era progress with the scientific objectives of nuclear and high-energy physics.

preprint2021arXiv

Stochastic Dynamics and Bound States of Heavy Impurities in a Fermi Bath

We investigate the dynamics of heavy impurities embedded in an ultra-cold Fermi gas by using a Generalized Langevin equation. The latter -- derived by means of influence functional theory -- describes the stochastic classical dynamics of the impurities and the quantum nature of the fermionic bath manifests in the emergent interaction between the impurities and in the viscosity tensor. By focusing on the two-impurity case, we predict the existence of bound states, in different conditions of coupling and temperature, and whose life-time can be analytically estimated. Our predictions should be testable using cold-gases platforms within current technology.

preprint2020arXiv

Spectral density estimation with the Gaussian Integral Transform

The spectral density operator $\hatρ(ω)=δ(ω-\hat{H})$ plays a central role in linear response theory as its expectation value, the dynamical response function, can be used to compute scattering cross-sections. In this work, we describe a near optimal quantum algorithm providing an approximation to the spectral density with energy resolution $Δ$ and error $ε$ using $\mathcal{O}\left(\sqrt{\log\left(1/ε\right)\left(\log\left(1/Δ\right)+\log\left(1/ε\right)\right)}/Δ\right)$ operations. This is achieved without using expensive approximations to the time-evolution operator but exploiting instead qubitization to implement an approximate Gaussian Integral Transform (GIT) of the spectral density. We also describe appropriate error metrics to assess the quality of spectral function approximations more generally.

preprint2019arXiv

Quantum Computing for Neutrino-nucleus Scattering

Neutrino-nucleus cross section uncertainties are expected to be a dominant systematic in future accelerator neutrino experiments. The cross sections are determined by the linear response of the nucleus to the weak interactions of the neutrino, and are dominated by energy and distance scales of the order of the separation between nucleons in the nucleus. These response functions are potentially an important early physics application of quantum computers. Here we present an analysis of the resources required and their expected scaling for scattering cross section calculations. We also examine simple small-scale neutrino-nucleus models on modern quantum hardware. In this paper, we use variational methods to obtain the ground state of a three nucleon system (the triton) and then implement the relevant time evolution. In order to tame the errors in present-day NISQ devices, we explore the use of different error-mitigation techniques to increase the fidelity of the calculations.

preprint2019arXiv

Short-depth circuits for efficient expectation value estimation

The evaluation of expectation values $Tr\left[ρO\right]$ for some pure state $ρ$ and Hermitian operator $O$ is of central importance in a variety of quantum algorithms. Near optimal techniques developed in the past require a number of measurements $N$ approaching the Heisenberg limit $N=\mathcal{O}\left(1/ε\right)$ as a function of target accuracy $ε$. The use of Quantum Phase Estimation requires however long circuit depths $C=\mathcal{O}\left(1/ε\right)$ making their implementation difficult on near term noisy devices. The more direct strategy of Operator Averaging is usually preferred as it can be performed using $N=\mathcal{O}\left(1/ε^2\right)$ measurements and no additional gates besides those needed for the state preparation. In this work we use a simple but realistic model to describe the bound state of a neutron and a proton (the deuteron) and show that the latter strategy can require an overly large number of measurement in order to achieve a reasonably small relative target accuracy $ε_r$. We propose to overcome this problem using a single step of QPE and classical post-processing. This approach leads to a circuit depth $C=\mathcal{O}\left(ε^μ\right)$ (with $μ\geq0$) and to a number of measurements $N=\mathcal{O}\left(1/ε^{2+ν}\right)$ for $0<ν\leq1$. We provide detailed descriptions of two implementations of our strategy for $ν=1$ and $ν\approx0.5$ and derive appropriate conditions that a particular problem instance has to satisfy in order for our method to provide an advantage.

preprint2016arXiv

Microscopically constrained mean field models from chiral nuclear thermodynamics

We explore the use of mean field models to approximate microscopic nuclear equations of state derived from chiral effective field theory across the densities and temperatures relevant for simu- lating astrophysical phenomena such as core-collapse supernovae and binary neutron star mergers. We consider both relativistic mean field theory with scalar and vector meson exchange as well as energy density functionals based on Skyrme phenomenology and compare to thermodynamic equa- tions of state derived from chiral two- and three-nucleon forces in many-body perturbation theory. Quantum Monte Carlo simulations of symmetric nuclear matter and pure neutron matter are used to determine the density regimes in which perturbation theory with chiral nuclear forces is valid. Within the theoretical uncertainties associated with the many-body methods, we find that select mean field models describe well microscopic nuclear thermodynamics. As an additional consistency requirement, we study as well the single-particle properties of nucleons in a hot/dense environment, which affect e.g., charged-current weak reactions in neutron-rich matter. The identified mean field models can be used across a larger range of densities and temperatures in astrophysical simulations than more computationally expensive microscopic models.

preprint2016arXiv

Thermal conductivity and impurity scattering in the accreting neutron star crust

We calculate the thermal conductivity of electrons for the strongly correlated multi-component ion plasma expected in the outer layers of neutron star's crust employing a Path Integral Monte Carlo (PIMC) approach. This allows us to isolate the low energy response of the ions and use it to calculate the electron scattering rate and the electron thermal conductivity. We find that the scattering rate is enhanced by a factor 2-4 compared to earlier calculations based on the simpler electron-impurity scattering formalism. This findings directly impacts the interpretation of thermal relaxation observed in transiently accreting neutron stars and has implications for the composition and nuclear reactions in the crust that occur during accretion.

preprint2015arXiv

Constraining the nuclear energy density functional with quantum Monte Carlo calculations

We study the problem of an impurity in fully polarized (spin-up) low density neutron matter with the help of an accurate quantum Monte Carlo method in conjunction with a realistic nucleon-nucleon interaction derived from chiral effective field theory at next-to-next-to-leading-order. Our calculations show that the behavior of the proton spin-down impurity is very similar to that of a polaron in a fully polarized unitary Fermi gas. We show that our results can be used to put tight constraints on the time-odd parts of the energy density functional, independent of the time-even parts, in the density regime relevant to neutron-rich nuclei and compact astrophysical objects such as neutron stars and supernovae.

preprint2014arXiv

Quantum Monte Carlo calculations of neutron matter with non-local chiral interactions

We present fully non-perturbative quantum Monte Carlo calculations with non-local chiral effective field theory (EFT) interactions for the ground state properties of neutron matter. The equation of state, the nucleon chemical potentials and the momentum distribution in pure neutron matter up to one and a half times the nuclear saturation density are computed with a newly optimized chiral EFT interaction at next-to-next-to-leading order. This work opens the way to systematic order by order benchmarking of chiral EFT interactions, and \emph{ab initio} prediction of nuclear properties while respecting the symmetries of quantum chromodynamics.

preprint2013arXiv

Quantum Monte Carlo with Coupled-Cluster wave functions

We introduce a novel many body method which combines two powerful many body techniques, viz., quantum Monte Carlo and coupled cluster theory. Coupled cluster wave functions are introduced as importance functions in a Monte Carlo method designed for the configuration interaction framework to provide rigorous upper bounds to the ground state energy. We benchmark our method on the homogeneous electron gas in momentum space. The importance function used is the coupled cluster doubles wave function. We show that the computational resources required in our method scale polynomially with system size. Our energy upper bounds are in very good agreement with previous calculations of similar accuracy, and they can be systematically improved by including higher order excitations in the coupled cluster wave function.

preprint2012arXiv

Dynamical Structure Factors in Quantum Many-Body Systems from Quantum Monte Carlo Calculations

An ab-initio method for determining the dynamical structure function of an interacting many--body quantum system has been devised by combining a generalized integral transform method with Quantum Monte Carlo methods. As a first application, the coherent and, separately, the incoherent excitation spectrum of bulk atomic 4He has been computed, both in the low and intermediate momentum range. The peculiar form of the kernel in the integral transform of the dynamical structure function allows to predict, without using any model, both position and width of the collective excitations in the maxon--roton region, as well as the second collective peak. A prediction of the dispersion of the single--particle modes described by the incoherent part is also presented.