Source author record

Guglielmo Mazzola

Guglielmo Mazzola 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

18works
12topics
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

18 published item(s)

preprint2025arXiv

Quantum computing for genomics: conceptual challenges and practical perspectives

We assess the potential of quantum computing to accelerate computation of central tasks in genomics, focusing on often-neglected theoretical limitations. We discuss state-of-the-art challenges of quantum search, optimization, and machine learning algorithms. Examining database search with Grover's algorithm, we show that the expected speedup vanishes under realistic assumptions. For combinatorial optimization prevalent in genomics, we discuss the limitations of theoretical complexity in practice and suggest carefully identifying problems genuinely suited for quantum acceleration. Given the competition from excellent classical approximate solvers, quantum computing could offer a speedup in the near future only for a specific subset of hard enough tasks in assembly, gene selection, and inference. These tasks need to be characterized by core optimization problems that are particularly challenging for classical methods while requiring relatively limited variables. We emphasize rigorous empirical validation through runtime scaling analysis to avoid misleading claims of quantum advantage. Finally, we discuss the problem of trainability and data-loading in quantum machine learning. This work advocates for a balanced perspective on quantum computing in genomics, guiding future research toward targeted applications and robust validation.

preprint2022arXiv

Exponential challenges in unbiasing quantum Monte Carlo algorithms with quantum computers

Recently, Huggins et. al. [Nature, 603, 416-420 (2022)] devised a general projective Quantum Monte Carlo method suitable for implementation on quantum computers. This hybrid approach, however, relies on a subroutine -the computation of the local energy estimator on the quantum computer -that is intrinsically affected by an exponential scaling of the computational time with the number of qubits. By means of numerical experiments, we show that this exponential scaling manifests prominently already on systems below the point of "quantum advantage". For the prototypical transverse-field Ising model, we show that the required time resources to compete with classical simulations on around 40 qubits are already of the order of $10^{13}$ projective measurements, with an estimated running time of a few thousand years on superconducting hardware. These observations strongly suggest that the proposed hybrid method, in its present form, is unlikely to offer a sizeable advantage over conventional quantum Monte Carlo approaches.

preprint2022arXiv

Gauge invariant quantum circuits for $U(1)$ and Yang-Mills lattice gauge theories

Quantum computation represents an emerging framework to solve lattice gauge theories (LGT) with arbitrary gauge groups, a general and long-standing problem in computational physics. While quantum computers may encode LGT using only polynomially increasing resources, a major openissue concerns the violation of gauge-invariance during the dynamics and the search for groundstates. Here, we propose a new class of parametrized quantum circuits that can represent states belonging only to the physical sector of the total Hilbert space. This class of circuits is compact yet flexible enough to be used as a variational ansatz to study ground state properties, as well as representing states originating from a real-time dynamics. Concerning the first application, the structure of the wavefunction ansatz guarantees the preservation of physical constraints such as the Gauss law along the entire optimization process, enabling reliable variational calculations. As for the second application, this class of quantum circuits can be used in combination with timedependent variational quantum algorithms, thus drastically reducing the resource requirements to access dynamical properties.

preprint2022arXiv

Towards Quantum Advantage in Financial Market Risk using Quantum Gradient Algorithms

We introduce a quantum algorithm to compute the market risk of financial derivatives. Previous work has shown that quantum amplitude estimation can accelerate derivative pricing quadratically in the target error and we extend this to a quadratic error scaling advantage in market risk computation. We show that employing quantum gradient estimation algorithms can deliver a further quadratic advantage in the number of the associated market sensitivities, usually called greeks. By numerically simulating the quantum gradient estimation algorithms on financial derivatives of practical interest, we demonstrate that not only can we successfully estimate the greeks in the examples studied, but that the resource requirements can be significantly lower in practice than what is expected by theoretical complexity bounds. This additional advantage in the computation of financial market risk lowers the estimated logical clock rate required for financial quantum advantage from Chakrabarti et al. [Quantum 5, 463 (2021)] by a factor of ~7, from 50MHz to 7MHz, even for a modest number of greeks by industry standards (four). Moreover, we show that if we have access to enough resources, the quantum algorithm can be parallelized across 60 QPUs, in which case the logical clock rate of each device required to achieve the same overall runtime as the serial execution would be ~100kHz. Throughout this work, we summarize and compare several different combinations of quantum and classical approaches that could be used for computing the market risk of financial derivatives.

preprint2021arXiv

Improved accuracy on noisy devices by non-unitary Variational Quantum Eigensolver for chemistry applications

We propose a modification of the Variational Quantum Eigensolver algorithm for electronic structure optimization using quantum computers, named non-unitary Variational Quantum Eigensolver (nu-VQE), in which a non-unitary operator is combined with the original system Hamiltonian leading to a new variational problem with a simplified wavefunction Ansatz. In the present work, we use, as non-unitary operator, the Jastrow factor, inspired from classical Quantum Monte Carlo techniques for simulation of strongly correlated electrons. The method is applied to prototypical molecular Hamiltonians for which we obtain accurate ground state energies with shallower circuits, at the cost of an increased number of measurements. Finally, we also show that this method achieves an important error mitigation effect that drastically improves the quality of the results for VQE optimizations on today's noisy quantum computers. The absolute error in the calculated energy within our scheme is one order of magnitude smaller than the corresponding result using traditional VQE methods, with the same circuit depth.

preprint2020arXiv

Microcanonical and finite temperature ab initio molecular dynamics simulations on quantum computers

Ab initio molecular dynamics (AIMD) is a powerful tool to predict properties of molecular and condensed matter systems. The quality of this procedure is based on accurate electronic structure calculations. The development of quantum processors has shown great potential for the efficient evaluation of accurate ground and excited state energies of molecular systems, opening up new avenues for molecular dynamics simulations. In this work we address the use of variational quantum algorithms for the calculation of accurate atomic forces to be used in AIMD. In particular, we provide solutions for the alleviation of the statistical noise associated to the measurements of the expectation values of energies and forces, as well as schemes for the mitigation of the hardware noise sources (in particular, gate infidelities, qubit decoherence and readout errors). Despite the relative large error in the calculation of the potential energy, our results show that the proposed algorithms can provide reliable MD trajectories in the microcanonical (constant energy) ensemble. Further, exploiting the intrinsic noise arising from the quantum measurement process, we also propose a Langevin dynamics algorithm for the simulation of canonical, i.e., constant temperature, dynamics. Both algorithms (microcanonical and canonical) are applied to the simulation of simple molecular systems such as H2 and H3+. Finally, we also provide results for the dynamics of H2 obtained with IBM quantum computer ibmq_athens.

preprint2020arXiv

Non-adiabatic molecular quantum dynamics with quantum computers

The theoretical investigation of non-adiabatic processes is hampered by the complexity of the coupled electron-nuclear dynamics beyond the Born-Oppenheimer approximation. Classically, the simulation of such reactions is limited by the unfavourable scaling of the computational resources as a function of the system size. While quantum computing exhibits proven quantum advantage for the simulation of real-time dynamics, the study of quantum algorithms for the description of non-adiabatic phenomena is still unexplored. In this work, we propose a quantum algorithm for the simulation of fast non-adiabatic chemical processes together with an initialization scheme for quantum hardware calculations. In particular, we introduce a first-quantization method for the time evolution of a wavepacket on two coupled harmonic potential energy surfaces (Marcus model). In our approach, the computational resources scale polynomially in the system dimensions, opening up new avenues for the study of photophysical processes that are classically intractable.

preprint2020arXiv

TurboRVB: a many-body toolkit for {\it ab initio} electronic simulations by quantum Monte Carlo

TurboRVB is a computational package for {\it ab initio} Quantum Monte Carlo (QMC) simulations of both molecular and bulk electronic systems. The code implements two types of well established QMC algorithms: Variational Monte Carlo (VMC), and Diffusion Monte Carlo in its robust and efficient lattice regularized variant. A key feature of the code is the possibility of using strongly correlated many-body wave functions. The electronic wave function (WF) is obtained by applying a Jastrow factor, which takes into account dynamical correlations, to the most general mean-field ground state, written either as an antisymmetrized geminal product with spin-singlet pairing, or as a Pfaffian, including both singlet and triplet correlations. This wave function can be viewed as an efficient implementation of the so-called resonating valence bond (RVB) ansatz, first proposed by L. Pauling and P. W. Anderson in quantum chemistry and condensed matter physics, respectively. The RVB ansatz implemented in TurboRVB has a large variational freedom, including the Jastrow correlated Slater determinant as its simplest, but nontrivial case. Moreover, it has the remarkable advantage of remaining with an affordable computational cost, proportional to the one spent for the evaluation of a single Slater determinant. The code implements the adjoint algorithmic differentiation that enables a very efficient evaluation of energy derivatives, comprising the ionic forces. Thus, one can perform structural optimizations and molecular dynamics in the canonical NVT ensemble at the VMC level. For the electronic part, a full WF optimization is made possible thanks to state-of-the-art stochastic algorithms for energy minimization. The code has been efficiently parallelized by using a hybrid MPI-OpenMP protocol, that is also an ideal environment for exploiting the computational power of modern GPU accelerators.

preprint2020arXiv

Understanding dense hydrogen at planetary conditions

Materials at high pressures and temperatures are of great interest for planetary science and astrophysics, warm dense matter physics, and inertial confinement fusion research. Planetary structure models rely on our understanding of the behaviour of elements (and their mixtures) at exotic conditions that do not exist on Earth, and at the same time planets serve as natural laboratories for studying materials at extreme conditions. In this review we discuss the connection between modelling planetary interiors and high-pressure physics of hydrogen and helium. First, we summarise key experiments for determining the equation of state and phase diagram of hydrogen and helium as well as state-of-the-art theoretical approaches. We next briefly review our current knowledge of the internal structures of the giant planets in the Solar System, Jupiter and Saturn, and the importance of high pressure physics to their characterisation.

preprint2019arXiv

Evidence for supercritical behavior of high-pressure liquid hydrogen

Hydrogen exhibits unusual behaviors at megabar pressures, with consequences for planetary science, condensed matter physics and materials science. Experiments at such extreme conditions are challenging, often resulting in hard-to-interpret and controversial observations. We present a theoretical study of the phase diagram of dense hydrogen, using machine learning to overcome time and length scale limitations while describing accurately interatomic forces. We reproduce the re-entrant melting behavior and the polymorphism of the solid phase. In simulations based on the machine learning potential we find evidence for continuous metallization in the liquid, as a first-order liquid-liquid transition is pre-empted by freezing. This suggests a smooth transition between insulating and metallic layers in giant gas planets, and reconciles existing discrepancies between experiments as a manifestation of supercritical behavior.

preprint2019arXiv

Precise measurement of quantum observables with neural-network estimators

The measurement precision of modern quantum simulators is intrinsically constrained by the limited set of measurements that can be efficiently implemented on hardware. This fundamental limitation is particularly severe for quantum algorithms where complex quantum observables are to be precisely evaluated. To achieve precise estimates with current methods, prohibitively large amounts of sample statistics are required in experiments. Here, we propose to reduce the measurement overhead by integrating artificial neural networks with quantum simulation platforms. We show that unsupervised learning of single-qubit data allows the trained networks to accommodate measurements of complex observables, otherwise costly using traditional post-processing techniques. The effectiveness of this hybrid measurement protocol is demonstrated for quantum chemistry Hamiltonians using both synthetic and experimental data. Neural-network estimators attain high-precision measurements with a drastic reduction in the amount of sample statistics, without requiring additional quantum resources.

preprint2016arXiv

Accelerated nuclear quantum effects sampling with open path integrals

We numericaly demonstrate that, in double well models, the autocorrelation time of open path integral Monte Carlo simulations can be much smaller compared to standard ones using ring polymers. We also provide an intuitive explanation based on the role of instantons as transition states of the path integral pseudodynamics. Therefore we propose that, in all cases when the ground state approximation to the finite temperature partition function holds, open path integral simulations can be used to accelerate the sampling in realistic simulations aimed to explore nuclear quantum effects.

preprint2016arXiv

Geminal embedding scheme for optimal atomic basis set construction in correlated calculations

We introduce an efficient method to construct optimal and system adaptive basis sets for use in electronic structure and quantum Monte Carlo calculations. The method is based on an embedding scheme in which a reference atom is singled out from its environment, while the entire system (atom and environment) is described by a Slater determinant or its antisymmetrized geminal power (AGP) extension. The embedding procedure described here allows for the systematic and consistent contraction of the primitive basis set into geminal embedded orbitals (GEOs), with a dramatic reduction of the number of variational parameters necessary to represent the many-body wave function, for a chosen target accuracy. Within the variational Monte Carlo method, the Slater or AGP part is determined by a variational minimization of the energy of the whole system in presence of a flexible and accurate Jastrow factor, representing most of the dynamical electronic correlation. The resulting GEO basis set opens the way for a fully controlled optimization of many-body wave functions in electronic structure calculation of bulk materials, namely, containing a large number of electrons and atoms. We present applications on the water molecule, the volume collapse transition in cerium, and the high-pressure liquid hydrogen.

preprint2015arXiv

Ab-initio molecular dynamics simulation of liquid water by Quantum Monte Carlo

Although liquid water is ubiquitous in chemical reactions at roots of life and climate on the earth, the prediction of its properties by high-level ab initio molecular dynamics simulations still represents a formidable task for quantum chemistry. In this article we present a room temperature simulation of liquid water based on the potential energy surface obtained by a many-body wave function through quantum Monte Carlo (QMC) methods. The simulated properties are in good agreement with recent neutron scattering and X-ray experiments, particularly concerning the position of the oxygen-oxygen peak in the radial distribution function, at variance of previous Density Functional Theory attempts. Given the excellent performances of QMC on large scale supercomputers, this work opens new perspectives for predictive and reliable ab-initio simulations of complex chemical systems.

preprint2015arXiv

Understanding Quantum Tunneling through Quantum Monte Carlo Simulations

The tunneling between the two ground states of an Ising ferromagnet is a typical example of many-body tunneling processes between two local minima, as they occur during quantum annealing. Performing quantum Monte Carlo (QMC) simulations we find that the QMC tunneling rate displays the same scaling with system size, as the rate of incoherent tunneling. The scaling in both cases is $O(Δ^2)$, where $Δ$ is the tunneling splitting. An important consequence is that QMC simulations can be used to predict the performance of a quantum annealer for tunneling through a barrier. Furthermore, by using open instead of periodic boundary conditions in imaginary time, equivalent to a projector QMC algorithm, we obtain a quadratic speedup for QMC, and achieve linear scaling in $Δ$. We provide a physical understanding of these results and their range of applicability based on an instanton picture.

preprint2014arXiv

Distinct metallization and atomization transitions in dense liquid hydrogen

We perform molecular dynamics simulations driven by accurate Quantum Monte Carlo forces on dense liquid hydrogen. Recently it has been reported a complete atomization transition between a mixed-atomic liquid and a completely dissociated fluid in an almost unaccessible pressure range {[\emph{Nat. Commun.} {\bf 5}, 3487 (2014)]}. Here instead, in a much more interesting pressure range, we identify a different transition between the fully molecular liquid and the mixed-atomic fluid at $\sim$ 400 GPa, with numerical evidence supporting its metallic behavior. Therefore we predict that the metallization at finite temperature occurs in this partially dissociated molecular fluid, well before the complete atomization of the liquid. At high temperature this first-order transition becomes a crossover, in very good agreement with the experimental observation. Several systematic tests supporting the quality of our large scale calculations are also reported.

preprint2014arXiv

Unexpectedly high pressure for molecular dissociation in liquid hydrogen by a reliable electronic simulation

The study of the high pressure phase diagram of hydrogen has continued with renewed effort for about one century as it remains a fundamental challenge for experimental and theoretical techniques. Here we employ an efficient molecular dynamics based on the quantum Monte Carlo method, which can describe accurately the electronic correlation and treat a large number of hydrogen atoms, allowing a realistic and reliable prediction of thermodynamic roperties. We find that the molecular liquid phase is unexpectedly stable and the transition towards a fully atomic liquid phase occurs at much higher pressure than previously believed. The old standing problem of low temperature atomization is, therefore, still far from experimental reach.

preprint2012arXiv

Finite temperature electronic simulations beyond the Born-Oppenheimer approximation

We introduce a general technique to compute finite temperature electronic properties by a novel covariant formulation of the electronic partition function. By using a rigorous variational upper bound to the free energy we are led to the evaluation of a partition function that can be computed stochastically by sampling electronic wave functions and atomic positions (assumed classical). In order to achieve this target we show that it is extremely important to consider the non trivial geometry of the space defined by the wave function ansatz. The method can be extended to any technique capable to provide an energy value over a given wave function ansatz depending on several variational parameters and atomic positions. In particular we can take into account electronic correlation, by using the standard variational quantum Monte Carlo method, that has been so far limited to zero temperature ground state properties. We show that our approximation reduces correctly to the standard Born-Oppenheimer (BO) one at zero temperature and to the correct high temperature limit. At large enough temperatures this method allows to improve the BO, providing lower values of the electronic free energy, because within this method it is possible to take into account the electron entropy. We test this new method on the simple hydrogen molecule, where at low temperature we recover the correct BO low temperature limit. Moreover, we show that the dissociation of the molecule is possible at a temperature much smaller than the BO prediction. Several extension of the proposed technique are also discussed, as for instance the calculation of critical (magnetic, superconducting) temperatures, or transition rates in chemical reactions.