Source author record

Jonathan E. Moussa

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

ResearcherUnclaimed source record

quant-ph physics.comp-ph cond-mat.mtrl-sci physics.chem-ph cond-mat.str-el cond-mat.stat-mech cond-mat.mes-hall

Catalog footprint

What is connected

16works

7topics

4close collaborators

Actions

Connect this record

Open graph Browse works

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.

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.

preprint2016arXiv

Linear embedding of free energy minimization

Exact free energy minimization is a convex optimization problem that is usually approximated with stochastic sampling methods. Deterministic approximations have been less successful because many desirable properties have been difficult to attain. Such properties include the preservation of convexity, lower bounds on free energy, and applicability to systems without subsystem structure. We satisfy all of these properties by embedding free energy minimization into a linear program over energy-resolved expectation values. Numerical results on small systems are encouraging, but a lack of size consistency necessitates further development for large systems.

preprint2016arXiv

Minimax rational approximation of the Fermi-Dirac distribution

Accurate rational approximations of the Fermi-Dirac distribution are a useful component in many numerical algorithms for electronic structure calculations. The best known approximations use $O( \log (βΔ) \log (ε^{-1}))$ poles to achieve an error tolerance $ε$ at temperature $β^{-1}$ over an energy interval $Δ$. We apply minimax approximation to reduce the number of poles by a factor of four and replace $Δ$ with $Δ_{\mathrm{occ}}$, the occupied energy interval. This is particularly beneficial when $Δ\gg Δ_{\mathrm{occ}}$, such as in electronic structure calculations that use a large basis set.

preprint2016arXiv

Quantum circuits for qubit fusion

We consider four-dimensional qudits as qubit pairs and their qudit Pauli operators as qubit Clifford operators. This introduces a nesting, $C_1^2 \subset C_2^4 \subset C_3^2$, where $C_n^m$ is the $n$th level of the $m$-dimensional qudit Clifford hierarchy. If we can convert between logical qubits and qudits, then qudit Clifford operators are qubit non-Clifford operators. Conversion is achieved by qubit fusion and qudit fission using stabilizer circuits that consume a resource state. This resource is a fused qubit stabilizer state with a fault-tolerant state preparation using stabilizer circuits.

preprint2016arXiv

Transversal Clifford gates on folded surface codes

Surface and color codes are two forms of topological quantum error correction in two spatial dimensions with complementary properties. Surface codes have lower-depth error detection circuits and well-developed decoders to interpret and correct errors, while color codes have transversal Clifford gates and better code efficiency in the number of physical qubits needed to achieve a given code distance. A formal equivalence exists between color codes and folded surface codes, but it does not guarantee the transferability of any of these favorable properties. However, the equivalence does imply the existence of constant-depth circuit implementations of logical Clifford gates on folded surface codes. We achieve and improve this result by constructing two families of folded surface codes with transversal Clifford gates. This construction is presented generally for qudits of any dimension. The specific application of these codes to universal quantum computation based on qubit fusion is also discussed.

preprint2015arXiv

Multivalley effective mass theory simulation of donors in silicon

Last year, Salfi et al. made the first direct measurements of a donor wave function and found extremely good theoretical agreement with atomistic tight-binding [Salfi et al., Nat. Mater. 13, 605 (2014)]. Here, we show that multi-valley effective mass theory, applied properly, does achieve close agreement with tight-binding and hence gives reliable predictions. To demonstrate this, we variationally solve the coupled six-valley Shindo-Nara equations, including silicon's full Bloch functions. Surprisingly, we find that including the full Bloch functions necessitates a tetrahedral, rather than spherical, donor central cell correction to accurately reproduce the experimental energy spectrum of a phosphorus impurity in silicon. We cross-validate this method against atomistic tight-binding calculations, showing that the two theories agree well for the calculation of donor-donor tunnel coupling. Further, we benchmark our results by performing a statistical uncertainty analysis, confirming that derived quantities such as the wave function profile and tunnel couplings are robust with respect to variational energy fluctuations. Finally, we apply this method to exhaustively enumerate the tunnel coupling for all donor-donor configurations within a large search volume, demonstrating conclusively that the tunnel coupling has no spatially stable regions. Though this instability is problematic for reliably coupling donor pairs for two-qubit operations, we identify specific target locations where donor qubits can be placed with scanning tunneling microscopy technology to achieve reliably large tunnel couplings.

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

Comment on "Fast and Accurate Modeling of Molecular Atomization Energies with Machine Learning"

A Comment on the Letter by M. Rupp et al., Phys. Rev. Lett. 108 058301 (2012).

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.

Jonathan E. Moussa

What is connected

Connect this record

See the researcher in context

Building this map preview

16 published item(s)

Low-Depth Quantum Metropolis Algorithm

Engineered thermalization and cooling of quantum many-body systems

Robust decompositions of quantum states

Linear embedding of free energy minimization

Minimax rational approximation of the Fermi-Dirac distribution

Quantum circuits for qubit fusion

Transversal Clifford gates on folded surface codes

Multivalley effective mass theory simulation of donors in silicon

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

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

Quadratic constrained mixed discrete optimization with an adiabatic quantum optimizer

Comment on "Adiabatic Quantum Algorithm for Search Engine Ranking"

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

Comment on "Fast and Accurate Modeling of Molecular Atomization Energies with Machine Learning"

Generalized unitary Bogoliubov transformation that breaks fermion number parity

Approximate diagonalization method for many-fermion Hamiltonians