D-wave Solution Research Articles

ILP-based resource optimization realized by quantum annealing for optical wide-area communication networks—A framework for solving combinatorial problems of a real-world application by quantum annealing

Resource allocation of wide-area internet networks is inherently a combinatorial optimization problem that if solved quickly, could provide near real-time adaptive control of internet-protocol traffic ensuring increased network efficacy and robustness, while minimizing energy requirements coming from power-hungry transceivers. In recent works we demonstrated how such a problem could be cast as a quadratic unconstrained binary optimization (QUBO) problem that can be embedded onto the D-Wave Advantage™ quantum annealer system, demonstrating proof of principle. Our initial studies left open the possibility for improvement of D-Wave solutions via judicious choices of system run parameters. Here we report on our investigations for optimizing these system parameters, and how we incorporate machine learning (ML) techniques to further improve on the quality of solutions. In particular, we use the Hamming distance to investigate correlations between various system-run parameters and solution vectors. We then apply a decision tree neural network (NN) to learn these correlations, with the goal of using the neural network to provide further guesses to solution vectors. We successfully implement this NN in a simple integer linear programming (ILP) example, demonstrating how the NN can fully map out the solution space that was not captured by D-Wave. We find, however, for the 3-node network problem the NN is not able to enhance the quality of space of solutions.

Comparison between a quantum annealer and a classical approximation algorithm for computing the ground state of an Ising spin glass.

Finding the ground state of an Ising spin glass on general graphs belongs to the class of NP-hard problems, widely believed to have no efficient polynomial-time algorithms to solve them. An approach developed in computer science for dealing with such problems is to devise approximation algorithms; these are algorithms, whose run time scales polynomially with the input size, that provide solutions with provable guarantees on their quality in terms of the optimal unknown solution. Recently, several algorithms for the Ising spin-glass problem on a bounded degree graph that provide different approximation guarantees were introduced. D-Wave, a Canadian-based company, has constructed a physical realization of a quantum annealer and has enabled researchers and practitioners to access it via their cloud service. D-Wave is particularly suited for computing an approximation for the ground state of an Ising spin glass on its Chimera and Pegasus graphs-both with a bounded degree. To assess the quality of D-Wave's solution, it is natural to compare it to classical approximation algorithms specifically designed to solve the same problem. In this work, we compare the performance of a recently developed approximation algorithm to solve the Ising spin-glass problem on graphs of bounded degree against the performance of the D-Wave computer. We also compared the performance of D-Wave's computer in the Chimera architecture against the performance of a heuristic tailored specifically to handle the Chimera graph. We found that the D-Wave computer was able to find better approximations for all the random instances of the problem we studied-Gaussian weights, uniform weights, and discrete binary weights. Furthermore, the convergence times of D-Wave's computer were also significantly better. These results indicate the merit of D-Wave's computer under certain specific instances. More broadly, our method is relevant to a wider class of performance comparison studies, and we suggest that it is important to compare the performance of quantum computers not only against exact classical algorithms with exponential run-time scaling, but also against approximation algorithms with polynomial run-time scaling and a provable guarantee of performance.

Pairing symmetry of the one-band Hubbard model in the paramagnetic weak-coupling limit: A numerical RPA study

We study the spin-fluctuation-mediated superconducting pairing gap in a weak-coupling approach to the Hubbard model for a two dimensional square lattice in the paramagnetic state. Performing a comprehensive theoretical study of the phase diagram as a function of filling, we find that the superconducting gap exhibits transitions from p-wave at very low electron fillings to d_{x^2-y^2}-wave symmetry close to half filling in agreement with previous reports. At intermediate filling levels, different gap symmetries appear as a consequence of the changes in the Fermi surface topology and the associated structure of the spin susceptibility. In particular, the vicinity of a van Hove singularity in the electronic structure close to the Fermi level has important consequences for the gap structure in favoring the otherwise sub-dominant triplet solution over the singlet d-wave solution. By solving the full gap equation, we find that the energetically favorable triplet solutions are chiral and break time reversal symmetry. Finally, we also calculate the detailed angular gap structure of the quasi-particle spectrum, and show how spin-fluctuation-mediated pairing leads to significant deviations from the first harmonics both in the singlet d_{x^2-y^2} gap as well as the chiral triplet gap solution.

Chirald-wave superconductivity on the honeycomb lattice close to the Mott state

We study superconductivity on the honeycomb lattice close to the Mott state at half-filling. Due to the sixfold lattice symmetry and disjoint Fermi surfaces at opposite momenta, we show that several different fully gapped superconducting states naturally exist on the honeycomb lattice, of which the chiral $d+id'$-wave state has previously been shown to appear when superconductivity appears close to the Mott state. Using renormalized mean-field theory to study the t-J model and quantum Monte Carlo calculations of the Hubbard-U model we show that the $d+id'$-wave state is the favored superconducting state for a wide range of on-site repulsion U, from the intermediate to the strong coupling regime. We also investigate the possibility of a mixed chirality d-wave state, where the overall chirality cancels. We find that a state with $d+id'$-wave symmetry in one valley but $d-id'$-wave symmetry in the other valley is not possible in the t-J model without reducing the translational symmetry, due to the zero-momentum and spin-singlet nature of the superconducting order parameter. Moreover, any extended unit cells result either in disjoint Dirac points, which cannot harbor this mixed chirality state, or the two valleys are degenerate at the zone center, where valley hybridization prevents different superconducting condensates. We also investigate extended unit cells where the overall chirality cancels in real space. For supercells containing up to eight sites, including the Kekul\'{e} distortion, we find no energetically favorable d-wave solution with an overall zero chirality within the restriction of the t-J model.

Effect of nearest neighbor spin-singlet correlations in conventional graphene SNS Josephson junctions

Using the self-consistent tight-binding Bogoliubov-de Gennes formalism we have studied the effect of nearest neighbor spin-singlet bond (SB) correlations on Josephson coupling and proximity effect in graphene SNS Josephson junctions with conventional s-wave superconducting contacts. Despite the s-wave superconducting state in the contacts, the SB pairing state inside the junction has d-wave symmetry and clean, sharp interface junctions resemble a 'bulk-meets-bulk' situation with very little interaction between the two different superconducting states. In fact, due to a finite-size suppression of the superconducting state, a stronger SB coupling constant than in the bulk is needed in order to achieve SB pairing in a junction. For both short clean zigzag and armchair junctions a d-wave state that has a zero Josephson coupling to the s-wave state is chosen and therefore the Josephson current decreases when a SB pairing state develops in these junctions. In more realistic junctions, with smoother doping profiles and atomic scale disorder at the interfaces, it is possible to achieve some coupling between the contact s-wave state and the SB d-wave states. In addition, by breaking the appropriate lattice symmetry at the interface in order to induce another d-wave state, a non-zero Josephson coupling can be achieved which leads to a substantial increase in the Josephson current. We also report on the LDOS of the junctions and on a lack of zero energy states at interfaces despite the unconventional order parameters, which we attribute to the near degeneracy of the two d-wave solutions and their mixing at a general interface.

Resonating valence bonds and mean-fieldd-wave superconductivity in graphite

We investigate the possibility of inducing superconductivity in a graphite layer by electronic correlation effects. We use a phenomenological microscopic Hamiltonian which includes nearest neighbor hopping and an interaction term which explicitly favors nearest neighbor spin-singlets through the well-known resonance valence bond (RVB) character of planar organic molecules. Treating this Hamiltonian in mean-field theory, allowing for bond-dependent variation of the RVB order parameter, we show that both s- and d-wave superconducting states are possible. The d-wave solution belongs to a two-dimensional representation and breaks time reversal symmetry. At zero doping there exists a quantum critical point at the dimensionless coupling J/t = 1.91 and the s- and d-wave solutions are degenerate for low temperatures. At finite doping the d-wave solution has a significantly higher Tc than the s-wave solution. By using density functional theory we show that the doping induced from sulfur absorption on a graphite layer is enough to cause an electronically driven d-wave superconductivity at graphite-sulfur interfaces. We also discuss applying our results to the case of the intercalated graphites as well as the validity of a mean-field approach.

S- and d-wave solution of Eliashberg equations with finite bandwidth

In this work, we discuss the results of the direct solution of the Eliashberg equations with finite bandwidth, in the cases of s- and d-wave symmetry for the pair wave function and in the presence of scattering from impurities. We show that the reduction of the critical temperature T c due to the finite bandwidth depends on the value of the bandwidth itself, but is almost independent of the symmetry of the order parameter. The same happens for the shape of Z( ω) and Δ(ω). Moreover, we discuss the effect of the finite bandwidth on the shape of the quasiparticle density of states. The results clearly indicate that the infinite bandwidth approximation leads to an underestimation of the electron-boson coupling constant.

D-wave superconductivity in the hubbard model

The superconducting instabilities of the doped repulsive 2D Hubbard model are studied in the intermediate to strong coupling regime with the help of the dynamical cluster approximation. To solve the effective cluster problem we employ an extended noncrossing approximation, which allows for a transition to the broken symmetry state. At sufficiently low temperatures we find stable d-wave solutions with off-diagonal long-range order. The maximal T(c) approximately 150 K occurs for a doping delta approximately 20% and the doping dependence of the transition temperatures agrees well with the generic high- T(c) phase diagram.

Density-induced BCS to Bose-Einstein crossover

We investigate the zero-temperature BCS to Bose-Einstein crossover at the mean-field level, by driving it with the attractive potential and the particle density.We emphasize specifically the role played by the particle density in this crossover.Three different interparticle potentials are considered for the continuum model in three spatial dimensions, while both s- and d-wave solutions are analyzed for the attractive (extended) Hubbard model on a two-dimensional square lattice. For this model the peculiar behavior of the crossover for the d-wave solution is discussed.In particular, in the strong-coupling limit when approaching half filling we evidence the occurrence of strong correlations among antiparallel-spin fermions belonging to different composite bosons, which give rise to a quasi-long-range antiferromagnetic order in this limit.

S- and d-wave symmetries of the solutions of the Eliashberg equations

We examine the different possible symmetries of the superconducting gap obtained by solving the Eliashberg equations. We consider an electron-phonon interaction in a strong coupling scenario. The Coulomb pseudopotential plays the crucial role of providing the repulsion needed to favour the d-wave symmetry. But the key parameter that allows very anisotropic solutions even with very strong coupling is the small angular range of the interaction due to predominantly electron-phonon forward scattering that is found in the high- T c superconductors. We find both s- and d-wave solutions whose stability depends mainly on the angular range of the interaction.

Perturbed bifurcations in the BCS gap equation.

The anisotropic BCS gap equation is applied to a tight-binding model of layered high-temperature superconductors. The possible solutions have s-, d-, and mixed s- and d-wave symmetries using nearest-neighbor intralayer singlet pairing interaction of the same strength in the x and y directions. The transitions from d- or s- to mixed s- and d-wave solutions result from pitchfork bifurcations. In the case of slightly different pairing strength in the x and y directions, perturbed pitchfork bifurcations emerge, leading to a dramatic change in the physical properties of the superconducting state.