Source author record

Masayuki Ohzeki

Masayuki Ohzeki 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

50works
14topics
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

50 published item(s)

preprint2026arXiv

Interpretation of Unfair Sampling in Quantum Annealing by Node Centrality

In applications where multiple optimal solutions are needed, transverse-field quantum annealing (QA) is known to sample degenerate ground states in a strongly biased manner. Despite extensive empirical observations, it remains unclear which features of degenerate ground states are preferentially sampled and why by QA. Here we analyze the final states using degenerate perturbation theory to characterize the preference among them. In this analysis, the adjacency matrix of the graph composed by the ground states naturally emerges, and we can predict the eigenvector centralities (one of the node centralities) are related to the probabilities of these states. We verify this prediction on toy models where degeneracy is lifted at first and second order, and we show that second-order weights encode local barrier information, relating sampling fairness to the flatness of the local energy landscape. Finally, this perspective suggests two practical routes toward fair sampling -- promoting connectivity of the graph and reducing heterogeneity of centralities -- and we illustrate consistency with higher-order drivers and minor-embedding transformations.

preprint2026arXiv

Kernel Learning for Regression via Quantum Annealing Based Spectral Sampling

While quantum annealing (QA) has been developed for combinatorial optimization, practical QA devices operate at finite temperature and under noise, and their outputs can be regarded as stochastic samples close to a Gibbs--Boltzmann distribution. In this study, we propose a QA-in-the-loop kernel learning framework that integrates QA not merely as a substitute for Markov-chain Monte Carlo sampling but as a component that directly determines the learned kernel for regression. Based on Bochner's theorem, a shift-invariant kernel is represented as an expectation over a spectral distribution, and random Fourier features (RFF) approximate the kernel by sampling frequencies. We model the spectral distribution with a (multi-layer) restricted Boltzmann machine (RBM), generate discrete RBM samples using QA, and map them to continuous frequencies via a Gaussian--Bernoulli transformation. Using the resulting RFF, we construct a data-adaptive kernel and perform Nadaraya--Watson (NW) regression. Because the RFF approximation based on $\cos(\bmω^{\top}Δ\bm{x})$ can yield small negative values and cancellation across neighbors, the Nadaraya--Watson denominator $\sum_j k_{ij}$ may become close to zero. We therefore employ nonnegative squared-kernel weights $w_{ij}=k(\bm{x}_i,\bm{x}_j)^2$, which also enhances the contrast of kernel weights. The kernel parameters are trained by minimizing the leave-one-out NW mean squared error, and we additionally evaluate local linear regression with the same squared-kernel weights at inference. Experiments on multiple benchmark regression datasets demonstrate a decrease in training loss, accompanied by structural changes in the kernel matrix, and show that the learned kernel tends to improve $R^2$ and RMSE over the baseline Gaussian-kernel NW. Increasing the number of random features at inference further enhances accuracy.

preprint2025arXiv

Graph minor embedding can affect sampling degenerate ground states using quantum annealing

Quantum annealing, as currently implemented in hardware, cannot fairly sample all ground states. Graph minor embedding, which maps a problem to the hardware graph of quantum annealers, affects sampling all states. In this study, we demonstrate the influence of graph minor embedding on fair sampling of degenerate ground states. For two embedded models that introduce auxiliary variables, numerical simulations of Schrödinger evolution revealed that fairness varies significantly depending on the embedding, and the chain strength is related to ground-state fairness. Using perturbation theory, we found that chain strength determines the energy landscape around ground states, with flatter landscapes having higher probabilities of being obtained.

preprint2022arXiv

Sparse Signal Reconstruction with QUBO Formulation in l0-regularized Linear Regression

An l0-regularized linear regression for a sparse signal reconstruction is implemented based on the quadratic unconstrained binary optimization (QUBO) formulation. In this method, the signal values are quantized and expressed as bit sequences. By transforming l0-norm to a quadratic form of these bits, the fully quadratic objective function is provided and optimized by the solver specialized for QUBO, such as the quantum annealer. Numerical experiments with a commercial quantum annealer show that the proposed method performs slightly better than conventional methods based on orthogonal matching pursuit (OMP) and the least absolute shrinkage and selection operator (LASSO) under several limited conditions.

preprint2021arXiv

Benchmark of quantum-inspired heuristic solvers for quadratic unconstrained binary optimization

Recently, inspired by quantum annealing, many solvers specialized for unconstrained binary quadratic programming problems have been developed. For further improvement and application of these solvers, it is important to clarify the differences in their performance for various types of problems. In this study, the performance of four quadratic unconstrained binary optimization problem solvers, namely D-Wave Hybrid Solver Service (HSS), Toshiba Simulated Bifurcation Machine (SBM), Fujitsu DigitalAnnealer (DA), and simulated annealing on a personal computer, was benchmarked. The problems used for benchmarking were instances of real problems in MQLib, instances of the SAT-UNSAT phase transition point of random not-all-equal 3-SAT(NAE 3-SAT), and the Ising spin glass Sherrington-Kirkpatrick (SK) model. Concerning MQLib instances, the HSS performance ranked first; for NAE 3-SAT, DA performance ranked first; and regarding the SK model, SBM performance ranked first. These results may help understand the strengths and weaknesses of these solvers.

preprint2020arXiv

A simple relation between frustration and transition points in diluted spin glasses

We investigate a possible relation between frustration and phase-transition points in spin glasses. The relation is represented as a condition of the number of frustrated plaquettes in the lattice at phase-transition points at zero temperature and was reported to provide very close points to the phase-transition points for several lattices. Although there has been no proof of the relation, the good correspondence in several lattices suggests the validity of the relation and some important role of frustration in the phase transitions. To examine the relation further, we present a natural extension of the relation to diluted lattices and verify its effectiveness for bond-diluted square lattices. We then confirm that the resulting points are in good agreement with the phase-transition points in a wide range of dilution rate. Our result supports the suggestion from the previous work for non-diluted lattices on the importance of frustration to the phase transition of spin glasses.

preprint2020arXiv

Breaking limitation of quantum annealer in solving optimization problems under constraints

Quantum annealing is a generic solver for optimization problems that uses fictitious quantum fluctuation. The most groundbreaking progress in the research field of quantum annealing is its hardware implementation, i.e., the so-called quantum annealer, using artificial spins. However, the connectivity between the artificial spins is sparse and limited on a special network known as the chimera graph. Several embedding techniques have been proposed, but the number of logical spins, which represents the optimization problems to be solved, is drastically reduced. In particular, an optimization problem including fully or even partly connected spins suffers from low embeddable size on the chimera graph. In the present study, we propose an alternative approach to solve a large-scale optimization problem on the chimera graph via a well-known method in statistical mechanics called the Hubbard-Stratonovich transformation or its variants. The proposed method can be used to deal with a fully connected Ising model without embedding on the chimera graph and leads to nontrivial results of the optimization problem. We tested the proposed method with a number of partition problems involving solving linear equations and the traffic flow optimization problem in Sendai and Kyoto cities in Japan.

preprint2020arXiv

Inequality for local energy of Ising models with quenched randomness and its application

In this study, we extend the lower bound on the average of the local energy of the Ising model with quenched randomness [J. Phys. Soc. Jpn. 76, 074711 (2007)] obtained for a symmetric distribution to an asymmetric one. Compared with the case of symmetric distribution, our bound has a non-trivial term. By applying the acquired bound to a Gaussian distribution, we obtain the lower bounds on the expectation of the square of the correlation function. Thus, we demonstrate that in the Ising model in a Gaussian random field, the spin-glass order parameter generally has a finite value at any temperature, regardless of the forms of the other interactions.

preprint2020arXiv

Probing the Universality of Topological Defect Formation in a Quantum Annealer: Kibble-Zurek Mechanism and Beyond

The number of topological defects created in a system driven through a quantum phase transition exhibits a power-law scaling with the driving time. This universal scaling law is the key prediction of the Kibble-Zurek mechanism (KZM), and testing it using a hardware-based quantum simulator is a coveted goal of quantum information science. Here we provide such a test using quantum annealing. Specifically, we report on extensive experimental tests of topological defect formation via the one-dimensional transverse-field Ising model on two different D-Wave quantum annealing devices. We find that the quantum simulator results can indeed be explained by the KZM for open-system quantum dynamics with phase-flip errors, with certain quantitative deviations from the theory likely caused by factors such as random control errors and transient effects. In addition, we probe physics beyond the KZM by identifying signatures of universality in the distribution and cumulants of the number of kinks and their decay, and again find agreement with the quantum simulator results. This implies that the theoretical predictions of the generalized KZM theory, which assumes isolation from the environment, applies beyond its original scope to an open system. We support this result by extensive numerical computations. To check whether an alternative, classical interpretation of these results is possible, we used the spin-vector Monte Carlo model, a candidate classical description of the D-Wave device. We find that the degree of agreement with the experimental data from the D-Wave annealing devices is better for the KZM, a quantum theory, than for the classical spin-vector Monte Carlo model, thus favoring a quantum description of the device. Our work provides an experimental test of quantum critical dynamics in an open quantum system, and paves the way to new directions in quantum simulation experiments.

preprint2020arXiv

Some inequalities for correlation functions of Ising models with quenched randomness

Correlation inequalities have played an essential role in the analysis of ferromagnetic models but have not been established in spin glass models. In this study, we obtain some correlation inequalities for the Ising models with quenched randomness, where the distribution of the interactions is symmetric. The acquired inequalities can be regarded as an extension of the previous results, which were limited to the local energy for a spin set, to the local energy for a pair of spin sets. Besides, we also obtain some correlation inequalities for asymmetric distribution.

preprint2020arXiv

Upper bound on the second derivative of the quenched pressure in spin-glass models: weak Griffiths second inequality

The Griffiths first and second inequalities have played an important role in the analysis of ferromagnetic models. In spin-glass models, although the counterpart of the Griffiths first inequality has been obtained, the counterpart of the Griffiths second inequality has not been established. In this study, we generalize the method in the previous work [J. Phys. Soc. Jpn. 76, 074711 (2007)] to the case with multi variables for both symmetric and non-symmetric distributions of the interactions, and derive some correlation inequalities for spin-glass models. Furthermore, by combining the acquired equalities in symmetric distributions, we show that there is a non-trivial positive upper bound on the second derivative of the quenched pressure with respect to the strength of the randomness, which is a weak result of the counterpart of the Griffiths second inequality in spin-glass models for general symmetric distributions.

preprint2019arXiv

Fair Sampling by Simulated Annealing on Quantum Annealer

Conventional quantum annealing does not sample all ground states fairly. We demonstrate that fair sampling can be achieved by performing simulated annealing on a quantum annealer. We discuss the problems that occur when implementing this method and propose an alternative way to overcome them. We numerically verify the fair sampling ability of our method in a small-scale toy model.

preprint2019arXiv

Sparse Modeling in Quantum Many-Body Problems

This review paper describes the basic concept and technical details of sparse modeling and its applications to quantum many-body problems. Sparse modeling refers to methodologies for finding a small number of relevant parameters that well explain a given dataset. This concept reminds us physics, where the goal is to find a small number of physical laws that are hidden behind complicated phenomena. Sparse modeling extends the target of physics from natural phenomena to data, and may be interpreted as "physics for data". The first half of this review introduces sparse modeling for physicists. It is assumed that readers have physics background but no expertise in data science. The second half reviews applications. Matsubara Green's function, which plays a central role in descriptions of correlated systems, has been found to be sparse, meaning that it contains little information. This leads to (i) a new method for solving the ill-conditioned inverse problem for analytical continuation, and (ii) a highly compact representation of Matsubara Green's function, which enables efficient calculations for quantum many-body systems.

preprint2016arXiv

Optimization of Markov process violates detailed balance condition

We consider the optimization of Markovian dynamics to pursue the fastest convergence to the stationary state. The brachistochrone method is applied to the continuous-time master equation for finite-size systems. The principle of least action leads to a brachistochrone equation for the transition-rate matrix. Three-state systems are explicitly analyzed, and we find that the solution violates the detailed balance condition. The properties of the solution are studied in detail to observe the optimality of the solution. We also discuss the counterdiabatic driving for the Markovian dynamics. The transition-rate matrix is then divided into two parts, and the state is given by an eigenstate of the first part. The second part violates the detailed balance condition and plays the role of a counterdiabatic term.

preprint2016arXiv

Quantum Monte Carlo simulation of a particular class of non-stoquastic Hamiltonians in quantum annealing

Quantum annealing is a generic solver of the optimization problem that uses fictitious quantum fluctuation. Its simulation in classical computing is often performed using the quantum Monte Carlo simulation via the Suzuki--Trotter decomposition. However, the negative sign problem sometimes emerges in the simulation of quantum annealing with an elaborate driver Hamiltonian, since it belongs to a class of non-stoquastic Hamiltonians. In the present study, we propose an alternative way to avoid the negative sign problem involved in a particular class of the non-stoquastic Hamiltonians. To check the validity of the method, we demonstrate our method by applying it to a simple problem that includes the anti-ferromagnetic XX interaction, which is a typical instance of the non-stoquastic Hamiltonians.

preprint2016arXiv

Statistical mechanical models of integer factorization problem

We formulate the integer factorization problem via a formulation of the searching problem for the ground state of a statistical mechanical Hamiltonian. The first passage time required to find a correct divisor of a composite number signifies the exponential computational hard- ness. Analysis of the density of states of two macroscopic quantities, i.e. the energy and the Hamming distance from the correct solutions, leads to the conclusion that the ground state (the correct solution) is completely isolated from the other low energy states, with the distance being proportional to the system size. In addition, the profile of the microcanonical entropy of the model has two peculiar features which are each related to two dramatic changes in the energy region sampled via Monte Carlo simulation or simulated annealing. Hence, we find a peculiar first-order phase transition in our model.

preprint2015arXiv

Full-order fluctuation-dissipation relation for a class of non-equilibrium steady states

Acceleration of relaxation toward a fixed stationary distribution via violation of detailed balance was reported in the context of a Markov chain Monte Carlo method recently. Inspired by this result, systematic methods to violate detailed balance in Langevin dynamics were formulated by using exponential and rotational nonconservative forces. In the present paper, we accentuate that such specific nonconservative forces relate to the large deviation of total heat in an equilibrium state. The response to these nonconservative forces can be described by the intrinsic fluctuation of the total heat in the equilibrium state. Consequently, the fluctuation-dissipation relation for nonequilibrium steady states is derived without recourse to a linear response approximation.

preprint2015arXiv

Hysteretic transition of coarse-grained power grid model on small-world network

We study synchronous phenomena of a coarse-grained power grid model, the swing equation, on small-world networks. We show that its steady state, which stands for the normal operation of the power systems, can be realized even if the phases are disordered. In addition, we prove the linear stability of steady state with small-different phases between the adjacent oscillators. On the other hand, a trigger of instantaneous power failure, which is described by the hysteretic transition, might disappear on an appropriate small-world network. This result suggests that the small-world connection would potentially prevent the massive blackouts.

preprint2015arXiv

Inverse Renormalization Group Transformation in Bayesian Image Segmentations

A new Bayesian image segmentation algorithm is proposed by combining a loopy belief propagation with an inverse real space renormalization group transformation to reduce the computational time. In results of our experiment, we observe that the proposed method can reduce the computational time to less than one-tenth of that taken by conventional Bayesian approaches.

preprint2015arXiv

L_1-regularized Boltzmann machine learning using majorizer minimization

We propose an inference method to estimate sparse interactions and biases according to Boltzmann machine learning. The basis of this method is $L_1$ regularization, which is often used in compressed sensing, a technique for reconstructing sparse input signals from undersampled outputs. $L_1$ regularization impedes the simple application of the gradient method, which optimizes the cost function that leads to accurate estimations, owing to the cost function's lack of smoothness. In this study, we utilize the majorizer minimization method, which is a well-known technique implemented in optimization problems, to avoid the non-smoothness of the cost function. By using the majorizer minimization method, we elucidate essentially relevant biases and interactions from given data with seemingly strongly-correlated components.

preprint2015arXiv

Langevin dynamics neglecting detailed balance condition

An improved method for driving a system into a desired distribution, for example, the Gibbs-Boltzmann distribution, is proposed, which makes use of an artificial relaxation process. The standard techniques for achieving the Gibbs-Boltzmann distribution involve numerical simulations under the detailed balance condition. In contrast, in the present study we formulate the Langevin dynamics, for which the corresponding Fokker-Planck operator includes an asymmetric component violating the detailed balance condition. This leads to shifts in the eigenvalues and results in the acceleration of the relaxation toward the steady state. The numerical implementation demonstrates faster convergence and shorter correlation time, and the technique of biased event sampling, Nemoto-Sasa theory, further highlights the efficacy of our method.

preprint2015arXiv

Statistical-mechanical analysis of pre-training and fine tuning in deep learning

In this paper, we present a statistical-mechanical analysis of deep learning. We elucidate some of the essential components of deep learning---pre-training by unsupervised learning and fine tuning by supervised learning. We formulate the extraction of features from the training data as a margin criterion in a high-dimensional feature-vector space. The self-organized classifier is then supplied with small amounts of labelled data, as in deep learning. Although we employ a simple single-layer perceptron model, rather than directly analyzing a multi-layer neural network, we find a nontrivial phase transition that is dependent on the number of unlabelled data in the generalization error of the resultant classifier. In this sense, we evaluate the efficacy of the unsupervised learning component of deep learning. The analysis is performed by the replica method, which is a sophisticated tool in statistical mechanics. We validate our result in the manner of deep learning, using a simple iterative algorithm to learn the weight vector on the basis of belief propagation.

preprint2015arXiv

Stochastic gradient method with accelerated stochastic dynamics

In this paper, we propose a novel technique to implement stochastic gradient methods, which are beneficial for learning from large datasets, through accelerated stochastic dynamics. A stochastic gradient method is based on mini-batch learning for reducing the computational cost when the amount of data is large. The stochasticity of the gradient can be mitigated by the injection of Gaussian noise, which yields the stochastic Langevin gradient method; this method can be used for Bayesian posterior sampling. However, the performance of the stochastic Langevin gradient method depends on the mixing rate of the stochastic dynamics. In this study, we propose violating the detailed balance condition to enhance the mixing rate. Recent studies have revealed that violating the detailed balance condition accelerates the convergence to a stationary state and reduces the correlation time between the samplings. We implement this violation of the detailed balance condition in the stochastic gradient Langevin method and test our method for a simple model to demonstrate its performance.

preprint2014arXiv

Detection of cheating by decimation algorithm

We expand the item response theory to study the case of "cheating students" for a set of exams, trying to detect them by applying a greedy algorithm of inference. This extended model is closely related to the Boltzmann machine learning. In this paper we aim to infer the correct biases and interactions of our model by considering a relatively small number of sets of training data. Nevertheless, the greedy algorithm that we employed in the present study exhibits good performance with a few number of training data. The key point is the sparseness of the interactions in our problem in the context of the Boltzmann machine learning: the existence of cheating students is expected to be very rare (possibly even in real world). We compare a standard approach to infer the sparse interactions in the Boltzmann machine learning to our greedy algorithm and we find the latter to be superior in several aspects.

preprint2014arXiv

High-precision phase diagram of spin glasses from duality analysis with real-space renormalization and graph polynomials

We propose a duality analysis for obtaining the critical manifold of two-dimensional spin glasses. Our method is based on the computation of quenched free energies with periodic and twisted periodic boundary conditions on a finite basis. The precision can be systematically improved by increasing the size of the basis, leading to very fast convergence towards the thermodynamic limit. We apply the method to obtain the phase diagrams of the random-bond Ising model and $q$-state Potts gauge glasses. In the Ising case, the Nishimori point is found at $p_N = 0.10929 \pm 0.00002$, in agreement with and improving on the precision of existing numerical estimations. Similar precision is found throughout the high-temperature part of the phase diagram. Finite-size effects are larger in the low-temperature region, but our results are in qualitative agreement with the known features of the phase diagram. In particular we show analytically that the critical point in the ground state is located at finite $p_0$.

preprint2013arXiv

A challenge for critical point of spin glass in ground state

We show several calculations to identify the critical point in the ground state in random spin systems including spin glasses on the basis of the duality analysis. The duality analysis is a profound method to obtain the precise location of the critical point in finite temperature even for spin glasses. We propose a single equality for identifying the critical point in the ground state from several speculations. The equality can indeed give the exact location of the critical points for the bond-dilution Ising model on several lattices and provides insight on further analysis on the ground state in spin glasses.

preprint2013arXiv

Belief propagation with multipoint correlations and its application to Inverse problem

We give explicit formulas of the Bethe approximation with multipoint correlations for systems with magnetic field. The obtained formulas include the closed form of the magnetization and the correlations between adjacent degrees of freedom. On the basis of our results, we propose a new iterative algorithm of the improved Bethe approximation. We confirm that the proposed technique is available for the random spin systems and indeed gives more accurate locations of the critical points. We discuss the possibility of the application of our method to the Inverse Ising model by use of these equations.

preprint2013arXiv

Duality with real-space renormalization and its application to bond percolation

We obtain the exact solution of the bond-percolation thresholds with inhomogenous probabilities on the square lattice. Our method is based on the duality analysis with real-space renormalization, which is a profound technique invented in the spin-glass theory. Our formulation is a more straightforward way compared to the very recent study on the same problem [R. M. Ziff, et. al., J. Phys. A: Math. Theor. 45 (2012) 494005]. The resultant generic formulas from our derivation can give several estimations for the bond-percolation thresholds on other lattices rather than the square lattice.

preprint2013arXiv

Nonequilibrium Work Relation in Macroscopic System

We reconsider a well-known relationship between the fluctuation theorem and the second law of thermodynamics by evaluating a probability measure-valued process. In order to establish a bridge between microscopic and macroscopic behaviors, we consider the thermodynamic limit of a stochastic dynamical system following the fundamental procedure often used in statistical mechanics. The thermodynamic path characterizing a macroscopic dynamical behavior can be formulated as an infimum of the action functional for the probability measure-valued process. In our formulation, the second law of thermodynamics can be derived by symmetry of the action functional, which is generated from the fluctuation theorem. We find that our formulation not only confirms that the ordinary Jarzynski equality in the thermodynamic limit can be rederived, but also enables us to establish a nontrivial nonequilibrium work relation for metastable states.

preprint2013arXiv

Variational Principle in Langevin Processes

The recent work, Nemoto and Sasa [Phys. Rev. E, 83: 030105(R) (2011)], has shown that large deviations of the current characterizing a nonequilibrium system are obtained by observing the typical current for a modified system specified by a variational principle. In the present study, we will give a generalized version of the Nemoto-Sasa study by extracting a hidden mathematical structure from the fluctuation-response relation which is well-known in statistical mechanics. Here, the minimization of the Kullback-Leibler divergence plays an essential role.

preprint2013arXiv

Violation of detailed balance accelerates relaxation

Recent studies have experienced the acceleration of convergence in Markov chain Monte Carlo methods implemented by the systems without detailed balance condition (DBC). However, such advantage of the violation of DBC has not been confirmed in general. We investigate the effect of the absence of DBC on the convergence toward equilibrium. Surprisingly, it is shown that the DBC violation always makes the relaxation faster. Our result implies the existence of a kind of thermodynamic inequality that connects the nonequilibrium process relaxing toward steady state with the relaxation process which has the same probability distribution as its equilibrium state.

preprint2012arXiv

Duality analysis on random planar lattice

The conventional duality analysis is employed to identify a location of a critical point on a uniform lattice without any disorder in its structure. In the present study, we deal with the random planar lattice, which consists of the randomized structure based on the square lattice. We introduce the uniformly random modification by the bond dilution and contraction on a part of the unit square. The random planar lattice includes the triangular and hexagonal lattices in extreme cases of a parameter to control the structure. The duality analysis in a modern fashion with real-space renormalization is found to be available for estimating the location of the critical points with wide range of the randomness parameter. As a simple testbed, we demonstrate that our method indeed gives several critical points for the cases of the Ising and Potts models, and the bond-percolation thresholds on the random planar lattice. Our method leads to not only such an extension of the duality analyses on the classical statistical mechanics but also a fascinating result associated with optimal error thresholds for a class of quantum error correction code, the surface code on the random planar lattice, which known as a skillful technique to protect the quantum state.

preprint2012arXiv

Error threshold estimates for surface code with loss of qubits

We estimate optimal thresholds for surface code in the presence of loss via an analytical method developed in statistical physics. The optimal threshold for the surface code is closely related to a special critical point in a finite-dimensional spin glass, which is disordered magnetic material. We compare our estimations to the heuristic numerical results reported in earlier studies. Further application of our method to the depolarizing channel, a natural generalization of the noise model, unveils its wider robustness even with loss of qubits.

preprint2012arXiv

Fluctuation Theorems on Nishimori Line

The distribution of the performed work for spin glasses with gauge symmetry is considered. With the aid of the gauge symmetry, which leads to the exact/rigorous results in spin glasses, we find a fascinating relation of the performed work as the fluctuation theorem. The integral form of the resultant relation reproduces the Jarzynski-type equation for spin glasses we have obtained. We show that similar relations can be established not only for the distribution of the performed work but also that of the free energy of spin glasses with gauge symmetry, which provides another interpretation of the phase transition in spin glasses.

preprint2012arXiv

Measurement-Based Quantum Computation on Symmetry Breaking Thermal States

We consider measurement-based quantum computation (MBQC) on thermal states of the interacting cluster Hamiltonian containing interactions between the cluster stabilizers that undergoes thermal phase transitions. We show that the long-range order of the symmetry breaking thermal states below a critical temperature drastically enhance the robustness of MBQC against thermal excitations. Specifically, we show the enhancement in two-dimensional cases and prove that MBQC is topologically protected below the critical temperature in three-dimensional cases. The interacting cluster Hamiltonian allows us to perform MBQC even at a temperature an order of magnitude higher than that of the free cluster Hamiltonian.

preprint2012arXiv

Strong Resilience of Topological Codes to Depolarization

The inevitable presence of decoherence effects in systems suitable for quantum computation necessitates effective error-correction schemes to protect information from noise. We compute the stability of the toric code to depolarization by mapping the quantum problem onto a classical disordered eight-vertex Ising model. By studying the stability of the related ferromagnetic phase both via large-scale Monte Carlo simulations and via the duality method, we are able to demonstrate an increased error threshold of 18.9(3)% when noise correlations are taken into account. Remarkably, this agrees within error bars with the result for a different class of codes-topological color codes-where the mapping yields interesting new types of interacting eight-vertex models.

preprint2011arXiv

Accuracy thresholds of topological color codes on the hexagonal and square-octagonal lattices

Accuracy thresholds of quantum error correcting codes, which exploit topological properties of systems, defined on two different arrangements of qubits are predicted. We study the topological color codes on the hexagonal lattice and on the square-octagonal lattice by the use of mapping into the spin glass systems. The analysis for the corresponding spin glass systems consists of the duality, and the gauge symmetry, which has succeeded in deriving locations of special points, which are deeply related with the accuracy thresholds of topological error correcting codes. We predict that the accuracy thresholds for the topological color codes would be $1-p_c = 0.1096-8 $ for the hexagonal lattice and $1-p_c = 0.1092-3$ for the square-octagonal lattice, where $1-p$ denotes the error probability on each qubit. Hence both of them are expected to be slightly lower than the probability $1-p_c = 0.110028$ for the quantum Gilbert-Varshamov bound with a zero encoding rate.

preprint2011arXiv

Jarzynski Equality for an Energy-Controlled System

The Jarzynski equality (JE) is known as an exact identity for nonequillibrium systems. The JE was originally formulated for isolated and isothermal systems, while Adib reported an JE extended to an isoenergetic process. In this paper, we extend the JE to an energy-controlled system. We make it possible to control the instantaneous value of the energy arbitrarily in a nonequilibrium process. Under our extension, the new JE is more practical and useful to calculate the number of states and the entropy than the isoenergetic one. We also show application of our JE to a kind of optimization problems.

preprint2011arXiv

Jarzynski Equality for an Energy-Controlled System (Proceedings of nanoPHYS'11)

We extend the Jarzynski equality, which is an exact identity between the equilibrium and nonequilibrium averages, to be useful to compute the value of the entropy difference by changing the Hamiltonian. To derive our result, we introduce artificial dynamics where the instantaneous value of the energy can be arbitrarily controlled during a nonequilibrium process. We establish an exact identity on such a process corresponding to the so-called Jarzynski equality. It is suggested that our formulation is valuable in a practical application as in optimization problems.

preprint2011arXiv

Nonequilibrium work on spin glasses in longitudinal and transverse fields

We derive a number of exact relations between equilibrium and nonequilibrium quantities for spin glasses in external fields using the Jarzynski equality and gauge symmetry. For randomly-distributed longitudinal fields, a lower bound is established for the work done on the system in nonequilibrium processes, and identities are proven to relate equilibrium and nonequilibrium quantities. In the case of uniform transverse fields, identities are proven between physical quantities and exponentiated work done to the system at different parts of the phase diagram with the context of quantum annealing in mind. Additional relations are given, which relate the exponentiated work in quantum and simulated (classical) annealing. It is also suggested that the Jarzynski equality may serve as a guide to develop a method to perform quantum annealing under non-adiabatic conditions.

preprint2011arXiv

Nonequilibrium work performed in quantum annealing

Quantum annealing is a generic solver of classical optimization problems that makes full use of quantum fluctuations. We consider work statistics given by a repetition of quantum annealing processes by employing the Jarzynski equality proposed in nonequilibrium statistical physics. In particular, we analyze the distribution of the work performed by a transverse field. A special symmetry, gauge symmetry, leads to a non-trivial relationship between quantum annealing toward different targets in the theory of spin glasses. We believe that our results will be a step toward an alternative realization of efficient quantum computation as well as our better understanding of nonequilibrium behavior of systems under quantum control.

preprint2011arXiv

Universality in phase boundary slopes for spin glasses on self dual lattices

We study the effects of disorder on the slope of the disorder--temperature phase boundary near the Onsager point (Tc = 2.269...) in spin-glass models. So far, studies have focused on marginal or irrelevant cases of disorder. Using duality arguments, as well as exact Pfaffian techniques we reproduce these analytical estimates. In addition, we obtain different estimates for spin-glass models on hierarchical lattices where the effects of disorder are relevant. We show that the phase-boundary slope near the Onsager point can be used to probe for the relevance of disorder effects.

preprint2010arXiv

Non-equilibrium Relations for Spin Glasses with Gauge Symmetry

We study the applications of non-equilibrium relations such as the Jarzynski equality and fluctuation theorem to spin glasses with gauge symmetry. It is shown that the exponentiated free-energy difference appearing in the Jarzynski equality reduces to a simple analytic function written explicitly in terms of the initial and final temperatures if the temperature satisfies a certain condition related to gauge symmetry. This result is used to derive a lower bound on the work done during the non-equilibrium process of temperature change. We also prove identities relating equilibrium and non-equilibrium quantities. These identities suggest a method to evaluate equilibrium quantities from non-equilibrium computations, which may be useful to avoid the problem of slow relaxation in spin glasses.

preprint2010arXiv

Nonequilibrium relations in spin glasses

The applications of nonequilbrium relations such as the Jarzynski equality and the fluctuation theorem to spin glasses are considered. The spin glass is a basic platform where we consider an application of an approximate solver of combinatorial optimization problems, simulated annealing. We find a novel relationship between an average through a nonequilibrium process where the temperature changes as in simulated annealing and a thermal average in equilibrium with different amounts of quenched randomness. The results shown in the present study may serve as an alternative way to overcome critical slowing down in spin glasses. It means that this way may mitigate difficulties in several hard optimization problems.

preprint2010arXiv

Quantum annealing with Jarzynski equality

We show a practical application of the Jarzynski equality in quantum computation. Its implementation may open a way to solve combinatorial optimization problems, minimization of a real single-valued function, cost function, with many arguments. We consider to incorpolate the Jarzynski equality into quantum annealing, which is one of the generic algorithms to solve the combinatorial optimization problem. The ordinary quantum annealing suffers from non-adiabatic transitions whose rate is characterized by the minimum energy gap $Δ_{\rm min.}$ of the quantum system under consideration. The quantum sweep speed is therefore restricted to be extremely slow for the achievement to obtain a solution without relevant errors. However, in our strategy shown in the present study, we find that such a difficulty would not matter.

preprint2010arXiv

Quantum Annealing with Jarzynski Equality (CCP2009)

We show a practical application of an well-known nonequilibrium relation, the Jarzynski equality, in quantum computation. Its implementation may open a way to solve combinatorial optimization problems, minimization of a real single-valued function, cost function, with many arguments. It has been disclosed that the ordinary quantum computational algorithm to solve a kind of hard optimization problems, has a bottleneck that its computational time is restricted to be extremely slow without relevant errors. However, by our novel strategy shown in the present study, we might overcome such a difficulty.

preprint2010arXiv

Quantum annealing: An introduction and new developments

Quantum annealing is a generic algorithm using quantum-mechanical fluctuations to search for the solution of an optimization problem. The present paper first reviews the fundamentals of quantum annealing and then reports on preliminary results for an alternative method. The review part includes the relationship of quantum annealing with classical simulated annealing. We next propose a novel quantum algorithm which might be available for hard optimization problems by using a classical-quantum mapping as well as the Jarzynski equality introduced in nonequilibrium statistical physics.

preprint2010arXiv

Zero-Temperature Complex Replica Zeros of the $\pm J$ Ising Spin Glass on Mean-Field Systems and Beyond

Zeros of the moment of the partition function $[Z^n]_{\bm{J}}$ with respect to complex $n$ are investigated in the zero temperature limit $β\to \infty$, $n\to 0$ keeping $y=βn \approx O(1)$. We numerically investigate the zeros of the $\pm J$ Ising spin glass models on several Cayley trees and hierarchical lattices and compare those results. In both lattices, the calculations are carried out with feasible computational costs by using recursion relations originated from the structures of those lattices. The results for Cayley trees show that a sequence of the zeros approaches the real axis of $y$ implying that a certain type of analyticity breaking actually occurs, although it is irrelevant for any known replica symmetry breaking. The result of hierarchical lattices also shows the presence of analyticity breaking, even in the two dimensional case in which there is no finite-temperature spin-glass transition, which implies the existence of the zero-temperature phase transition in the system. A notable tendency of hierarchical lattices is that the zeros spread in a wide region of the complex $y$ plane in comparison with the case of Cayley trees, which may reflect the difference between the mean-field and finite-dimensional systems.

preprint2009arXiv

Analytical evidence for the absence of spin glass transition on self-dual lattices

We show strong evidence for the absence of a finite-temperature spin glass transition for the random-bond Ising model on self-dual lattices. The analysis is performed by an application of duality relations, which enables us to derive a precise but approximate location of the multicritical point on the Nishimori line. This method can be systematically improved to presumably give the exact result asymptotically. The duality analysis, in conjunction with the relationship between the multicritical point and the spin glass transition point for the symmetric distribution function of randomness, leads to the conclusion of the absence of a finite-temperature spin glass transition for the case of symmetric distribution. The result is applicable to the random bond Ising model with $\pm J$ or Gaussian distribution and the Potts gauge glass on the square, triangular and hexagonal lattices as well as the random three-body Ising model on the triangular and the Union-Jack lattices and the four dimensional random plaquette gauge model. This conclusion is exact provided that the replica method is valid and the asymptotic limit of the duality analysis yields the exact location of the multicritical point

preprint2009arXiv

Multicritical point of spin glasses

We present a theoretical framework to accurately calculate the location of the multicritical point in the phase diagram of spin glasses. The result shows excellent agreement with numerical estimates. The basic idea is a combination of the duality relation, the replica method, and the gauge symmetry. An additional element of the renormalization group, in particular in the context of hierarchical lattices, leads to impressive improvements of the predictions.