Researcher profile

Jonathan E. Moussa

Jonathan E. Moussa contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
11works
0followers
6topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

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

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

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

Building this graph slice

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

Published work

11 published item(s)

preprint2022arXiv

Low-Depth Quantum Metropolis Algorithm

We construct a simple quantum version of the classical Metropolis algorithm to prepare and observe quantum thermal states. It induces both a quantum Markov chain that mixes the quantum thermal state and a classical Markov chain that mixes its observable measurements and enables a low-depth quantum circuit implementation. A Gaussian-filtered variant of quantum phase estimation enables thermalization times proportional to the reciprocal of temperature and the logarithm of biasing error. This matches the thermalization time of imaginary-time evolution, against which our algorithm performs favorably.

preprint2020arXiv

Engineered thermalization and cooling of quantum many-body systems

We develop a scheme for engineering genuine thermal states in analog quantum simulation platforms by coupling local degrees of freedom to driven, dissipative ancilla pseudospins. We demonstrate the scheme in a many-body quantum spin lattice simulation setting. A Born-Markov master equation describing the dynamics of the many-body system is developed, and we show that if the ancilla energies are periodically modulated, with a carefully chosen hierarchy of timescales, one can effectively thermalize the many-body system. Through analysis of the time-dependent dynamical generator, we determine the conditions under which the true thermal state is an approximate dynamical fixed point for general system Hamiltonians. Finally, we evaluate the thermalization protocol through numerical simulation and discuss prospects for implementation on current quantum simulation hardware.

preprint2020arXiv

Robust decompositions of quantum states

Classical-quantum computational complexity separations are an important motivation for the long-term development of digital quantum computers, but classical-quantum complexity equivalences are just as important in our present era of noisy intermediate-scale quantum devices for framing near-term progress towards quantum supremacy. We establish one such equivalence using a noisy quantum circuit model that can be simulated efficiently on classical computers. With respect to its noise model, quantum states have a robust decomposition into a sequence of operations that each extend the state by one qubit without spreading errors between qubits. This enables universal quantum sampling of states with an efficient representation in this robust form and observables with low quantum weight that can be sampled from general measurements on a few qubits and computational basis measurements on the remaining qubits. These robust decompositions are not unique, and we construct two distinct variants, both of which are compatible with machine-learning methodology. They both enable efficiently computable lower bounds on von Neumann entropy and thus can be used as finite-temperature variational quantum Monte Carlo methods.

preprint2014arXiv

Comment on "Self-Averaging Stochastic Kohn-Sham Density-Functional Theory"

In a recent Letter, Baer et al. present a stochastic method for Kohn-Sham density functional theory calculations. Their convergence criterion is the self-averaging total energy per electron, which requires a number of statistical samples that decreases with system size and enables a sublinear-scaling computational cost. However, the electron density, atomic forces, orbital energies, and many other physical quantities do not self-average, and this procedure causes their statistical errors to grow with system size. Convergence of non-self-averaging quantities requires that the number of samples be maintained with system size and is incompatible with a sublinear-scaling cost.

preprint2014arXiv

Cubic-scaling algorithm and self-consistent field for the random-phase approximation with second-order screened exchange

The random-phase approximation with second-order screened exchange (RPA+SOSEX) is a model of electron correlation energy with two caveats: its accuracy depends on an arbitrary choice of mean field, and it scales as $\mathcal{O}(n^5)$ operations and $\mathcal{O}(n^3)$ memory for $n$ electrons. We derive a new algorithm that reduces its scaling to $\mathcal{O}(n^3)$ operations and $\mathcal{O}(n^2)$ memory using controlled approximations and a new self-consistent field that approximates Brueckner coupled-cluster doubles (BCCD) theory with RPA+SOSEX, referred to as Brueckner RPA (BRPA) theory. The algorithm comparably reduces the scaling of second-order M$\mathrm{ø}$ller-Plesset (MP2) perturbation theory with smaller cost prefactors than RPA+SOSEX. Within a semiempirical model, we study H$_2$ dissociation to test accuracy and H$_n$ rings to verify scaling.

preprint2014arXiv

Quadratic constrained mixed discrete optimization with an adiabatic quantum optimizer

We extend the family of problems that may be implemented on an adiabatic quantum optimizer (AQO). When a quadratic optimization problem has at least one set of discrete controls and the constraints are linear, we call this a quadratic constrained mixed discrete optimization (QCMDO) problem. QCMDO problems are NP-hard, and no efficient classical algorithm for their solution is known. Included in the class of QCMDO problems are combinatorial optimization problems constrained by a linear partial differential equation (PDE) or system of linear PDEs. An essential complication commonly encountered in solving this type of problem is that the linear constraint may introduce many intermediate continuous variables into the optimization while the computational cost grows exponentially with problem size. We resolve this difficulty by developing a constructive mapping from QCMDO to quadratic unconstrained binary optimization (QUBO) such that the size of the QUBO problem depends only on the number of discrete control variables. With a suitable embedding, taking into account the physical constraints of the realizable coupling graph, the resulting QUBO problem can be implemented on an existing AQO. The mapping itself is efficient, scaling cubically with the number of continuous variables in the general case and linearly in the PDE case if an efficient preconditioner is available.

preprint2013arXiv

Comment on "Adiabatic Quantum Algorithm for Search Engine Ranking"

In their Letter, Garnerone et al. claim that an adiabatic quantum algorithm can extract information about a PageRank vector with either a polynomial or exponential reduction in time resources over the classical algorithm with comparable space resources. Here we argue that the quantum algorithm offers no obvious advantage over the classical algorithm and fails to preserve the fundamental stability property of the classical PageRank algorithm.

preprint2012arXiv

Accurate numerical integration of an electron exchange hole with a screened Coulomb interaction

The numerical implementation of an exchange-correlation functional is not always an accurate reproduction of its theoretical specification. For example, density functionals for exchange and correlation can be defined by an exchange or correlation hole function that is integrated with the Coulomb interaction to form an energy. This construction can be used to modify a density functional for use with any electron-electron interaction. Its most prominent use is in the Heyd-Scuseria-Ernzerhof (HSE) functional that generalizes the Perdew-Burke-Ernzerhof (PBE) model of exchange to a screened Coulomb interaction with an error function form. However, we find non-negligible numerical errors in the standard implementation of the HSE exchange hole integration. We formulate and implement a new method for evaluating the exchange hole integration that is simple, accurate, and efficient. Its numerical errors are bounded and minimized by applying basic elements of approximation theory.

preprint2012arXiv

Generalized unitary Bogoliubov transformation that breaks fermion number parity

The standard Bogoliubov transformation is generalized to enable fermion number parity breaking. The new transformation can diagonalize fermion Hamiltonians that are quadratic in fermion and number parity operators. This new variational freedom cannot lower the energy of Hamiltonians that conserve number parity. However, it has value as a minimal mechanism for continuously varying between states of different fermion number, regardless of number parity.

preprint2010arXiv

Approximate diagonalization method for many-fermion Hamiltonians

The limits of direct unitary transformation of many-fermion Hamiltonians are explored. Practical application of such transformations requires that effective many-body interactions be discarded over the course of a calculation. The truncation of the Hamiltonian leads to finite errors and in some cases divergences. A new formalism is proposed to manage errors and avoid divergences. Removing all interactions from a many-fermion Hamiltonian reduces it to fermion number operators allowing for direct calculation of eigenvalues. If the same transformations are applied to the bare fermions, eigenfermions are produced whose Slater determinants form eigenstates. This enables a hierarchy of diagonalization methods of increasing accuracy as fewer interactions are discarded from the Hamiltonian.