Error-Mitigated Quantum Simulation of Interacting Fermions with Trapped Ions

Quantum error mitigation has been extensively explored to increase the accuracy of the quantum circuits in noisy-intermediate-scale-quantum (NISQ) computation, where quantum error correction requiring additional quantum resources is not adopted. Among various error-mitigation schemes, probabilistic error cancellation (PEC) has been proposed as a general and systematic protocol that can be applied to numerous hardware platforms and quantum algorithms. However, PEC has only been tested in two-qubit systems and a superconducting multi-qubit system by learning a sparse error model. Here, we benchmark PEC using up to four trapped-ion qubits. For the benchmark, we simulate the dynamics of interacting fermions with or without spins by applying multiple Trotter steps. By tomographically reconstructing the error model and incorporating other mitigation methods such as positive probability and symmetry constraints, we are able to increase the fidelity of simulation and faithfully observe the dynamics of the Fermi-Hubbard model, including the different behavior of charge and spin of fermions. Our demonstrations can be an essential step for further extending systematic error-mitigation schemes toward practical quantum advantages.


INTRODUCTION
Recently, quantum computation in the noisy intermediate scale quantum (NISQ) regime has been actively developed for the possibility of reaching practical quantum advantages without quantum error correction [1,2].In this direction, quantum advantages with random gate sampling algorithms in superconducting systems have been demonstrated [3,4].The possibility of surpassing conventional classical computation in solving useful and practical problems has been seriously explored [5].In particular, a quantum-classical hybrid method that combines the advantages of classical and quantum computation, such as the variational quantum algorithm (VQA), has been theoretically and experimentally explored to reach higher performance than with only classical computers [6][7][8][9].One of the critical issues in NISQ computation is to obtain results with high precision from quantum devices without quantum error correction [10,11].
Various error-mitigation methods without using additional quantum resources for quantum error corrections have been proposed and applied to improve the computation results of quantum devices [12][13][14][15].Representative error-mitigation methods include zero-noise extrapolation [12][13][14][15][16][17], probabilistic error cancellation (PEC) [12][13][14][15][18][19][20][21], virtual distillation [22], etc.In particular, PEC is considered the general and systematic method for mitigating errors, which are not dependent on details of the physical implementation as long as errors are properly characterized.PEC consists of characterizing errors of each quantum gate, decomposing the inverse of the errors into combinations of basis operations, and sampling the operations with the characterized weights.In this way, we can estimate the expectation values of the desired observables without errors.PEC has been theoretically analyzed in detail and shown to be optimal for dephasing dominant noise situations [12,13,20,21].However, it has been experimentally tested with only two-qubit systems of superconducting circuits and trapped ions [18,19] and a multi-qubit system of the superconducting circuit via fitting a sparse error model motivated by the processor topology [23].
Systematic error-mitigation schemes like PEC would enhance the performance of quantum computation in the NISQ regime drastically.Many important quantum algorithms have been developed based on the fact that quantum dynamics can be efficiently simulated in quantum computation, which also inspires VQA designs [6][7][8][9].In digital quantum simulation, the time evolution of a certain Hamiltonian is typically implemented by the Trotter-Suzuki expansion [24].Besides the errors in the expansion, imperfect primitive quantum gates also degrade the fidelity of the quantum simulation and eventually lead to faulty results.Due to this reason, it is arXiv:2302.10436v1[quant-ph] 21 Feb 2023 an essential and urgent task to design and benchmark error-mitigation techniques that are both systematic and efficient in the NISQ era.As an instance for benchmark, the Fermi-Hubbard model [25] is particularly interesting.This model is one of the fundamental models for interacting electrons, and it has been extensively studied in classical manners, either analytically [26][27][28] or numerically [29].As experimental quantum technologies develop, understanding and insights for strongly-correlated quantum systems and high-temperature superconductivity [30] can be explored through quantum simulation of the Fermi-Hubbard model eventually beyond the classical computational capability.This model has been utilized in experimental demonstration of various quantum platforms [31][32][33][34].Recently, a combination of non-PEC methods has been applied to enhance simulation performance and observe the dynamics of the Fermi-Hubbard model in a superconducting system [35].
In this paper, we experimentally study the possibility and limitations of PEC by applying the method to a fourqubit system with trapped ions [36,37].In particular, we focus on simulating the dynamics of the 1D Fermi-Hubbard model by using multiple Trotter steps.We find the PEC method with a tomographically reconstructed error model significantly improves the fidelity of the simulation with up to three ion qubits.However, with four qubits, the improvement of fidelity by the PEC method alone is much worse than those for two-and three-qubit systems.This is because the tomographic characterization of the noisy gates is not precise enough to capture the time correlation and the cross-talk errors.With the application of total-spin and particle-number conservation [15,38,39] as well as positivity constraint, however, the enhanced PEC method provides recovery for the ideal dynamics of the 1D spin-Fermi-Hubbard model, where the spin and the charge behave differently.

FERMI-HUBBARD MODEL AND EXPERIMENTAL REALIZATION WITH TRAPPED IONS
Here we consider the extended Fermi-Hubbard (eFH) model, which describes a system of fermions moving in a one-dimensional lattice, as illustrated in Fig. 1a and Fig. 1b for one-component and two-component fermions.The one-component fermions can be considered spinless ones.The Hamiltonian reads with nl,λ ≡ ĉ † l,λ ĉl,λ , where ĉl,λ (ĉ † l,λ ) annihilates (creates) a fermion with spin λ (λ =↑, ↓) on the l-th site.Here J is the nearest-neighbor tunneling strength, and U and V respectively quantify the strength of the on-site and the nearest-neighbor interaction.
We can map the fermion system to a qubit system by using the Jordan-Wigner transformation, ĉl,↑ = ĉl,↓ = with σ± n = (σ x n ± iσ y n ) /2 and σx,y,z n being the Pauli operators on the n-th qubit.Thus generally speaking, we need N = 2L qubits to simulate the dynamics of an L-site fermionic chain.
As shown in Fig. 1c, we use 171 Yb + ions for qubit systems.The 171 Yb + ions are confined in a linear Paul trap.The hyperfine levels in the 2 S 1/2 ground-state manifold, i.e. |F = 0, m F = 0 and |F = 1, m F = 0 , are encoded as the qubit {|0 , |1 }, which can be initialized to |0 by optical pumping and disgustingly detected by fluorescence measurement via multi-channel photo-multiplier tube [40].
The qubit-system Hamiltonian becomes Note that according to the orientation of the Pauli operators, we divide the qubit-system Hamiltonian Ĥqb into three parts, each of which contains mutually commutative terms.
To extract dynamical properties, the essential part is to simulate the evolution operator Û (t) = exp −i Ĥqb t with a universal gate set available in a certain quantum computational platform.One of the well-known solutions is to use the Trotter-Suzuki decomposition.The firstorder Trotter-Suzuki decomposition of Û (t) reads with the Trotter step being δt = t/M .Figs. 1d and e show examples of the circuits for simulating the dynamics of 3-site spinless and 2-site spin-1/2 fermionic systems, with the native gate set including the Mølmer-Sorensen gates, i.e.XX ϕ , Y Y ϕ and ZZ ϕ , and single-qubit rotations.By applying the amplitude-shaped [41] global and individual laser beams on the desired ions and driving the transverse motional modes [42,43], we can realize the native Mølmer-Sørensen Y Y -gate (Y Y ϕ ).The gate XX ϕ and ZZ ϕ are composed with Y Y ϕ and corresponding single qubit rotations, as shown in Fig. 1f.
In this experiment, we consider two scenarios, one for spinless fermions, and the other for two-component fermions.In the former case, the on-site interaction strength naturally vanishes, i.e.U = 0, due to the fermionic anti-commutation relation.In this case, we use N = L ion-qubits to simulate the dynamics of a spinless fermionic chain with L sites.While for the latter case, we only consider the on-site interaction, since the interaction strength decreases as the distance increases and thus U V in most circumstances.

PROBABILISTIC ERROR CANCELLATION IN TRAPPED ION SYSTEM
Here, we introduce the PEC error-mitigation method, which systematically recovers ideal expectation values of target quantum circuits with the experimental noisy gates [14].The essential idea is to first characterize experimental noisy gates, then identify the noise parts by comparing the ideal and the noisy gates, and then decompose the inverse of the noise parts into combinations of basis operations.Then, each of the ideal gates in target quantum circuits is replaced by a corresponding noisy one followed by basis operations sampled with probability distributions obtained from the inverse-noise decomposition.At last, the ideal expectation value can be obtained by adding up those of the sampled circuits with appropriate weights.
In circuit quantum computation, a computational task can be compiled into a quantum circuit consisting of single-and two-qubit gates.In the trapped ion system, the error of the two-qubit gate is larger than the single-qubit one [19], and we can mitigate only the errors in two-qubit gates.To characterize the experimental noisy gates, we perform quantum process tomography (QPT), under the Pauli error assumption, to obtain the Pauli transfer matrices (PTMs) of the Mølmer-Sørensen Y Y gates on different qubit pairs, i.e.R Y Y ns ϕ (m,n) for the qubit pair (Q m , Q n ).Together with the ideal PTM R Y Yϕ , we can formally define an error operator for each of the noisy gates, with the PTM being Note that under the Pauli-error assumption, both R E and its inverse R −1 E can be decomposed into linear superpositions of PTMs of two-qubit Pauli operators P i ∈ {II, IX, IY, ..., ZZ}.As a result, ideal PTMs can be rewritten in terms of the PTMs of the noisy gate and the corresponding error operators, with the inverse-error decomposition Here the quasi-probabilities q i (m, n) are real and satisfy i q i (m, n) = 1.The fact that q i (m, n) can be negative indicates that the inverse of the error is not physical and thus can not be implemented by a deterministic quantum operation.To cancel the effect of the error operator, we rewrite the decomposition as follows where and sgn(•) is the sign function.
Then, we can use the Monte-Carlo sampling to compute the PEC error-mitigated results.Here we consider a target quantum circuit consisting of N g number of twoqubit gates Y Y ϕ (m, n), with the g-th gate labeled as G g , each of which is fully characterized and the inverse-error decomposition is written as R Gg = i q i,g R Pi R G ns g with q i,g being the quasi-probability of the two-qubit Pauli operator P i .The full set of the random circuits to be implemented in the experimental device denoted as S i : i = i 1 , . . ., i Ng , contains 16 Ng different circuits, each of which is obtained by adding P ig right after G g .To obtain the ideal expectation values, each random circuit S i is assigned with a probability ∝ g q ig,g and a sign sgn i ≡ g sgn q ig,g .As it is infeasible to implement all these circuits, we use the Monte Carlo sampling to obtain average values over the random circuits.Specifically, we generate N s random circuits by sampling P ig from the probability distribution ∝ q ig,g and obtain the expectation value µ i of an observable µ for the circuit S i by averaging 300 repetitions.The ideal expectation value µ is then calculated by where the cost is given by C ≡ Ng g=1 C g with the cost of the g-th gate being C g ≡ 15 ig=0 q ig,g .

ERROR MITIGATION AND SIMULATION OF TWO SPINLESS FERMIONIC MODES
We now discuss the simulation of two spinless fermionic modes, where the qubit-system Hamiltonian Eq. ( 3) reduces to With the Trotter-Suzuki expansion of Eq. ( 4), the evolution of Hamiltonian for two spinless fermionic modes can be experimentally realized with the specific quantum circuit shown in Fig. 2a.For the singly occupied fermion in one of the sites, the fermion oscillates between two sites caused by the nearest-neighbor tunneling.For the doubly occupied fermions, there will be no dynamics due to the Pauli exclusion principle of fermions.We apply M = 8 Trotter steps and each step consists of three Y Y ϕ gates and ten single-qubit rotations.The whole quantum circuit consists of twenty-four two-qubit gates and eighty-one single-qubit gates, where the additional one is for the initial state preparation.We set the nearestneighbor tunneling and interaction strength as V = 2J.Errors caused by gate imperfection accumulate as the number of Trotter steps increases and the fidelity of evolved dynamics would be degraded.We mitigate the errors and improve the results of dynamics.Since it is a two-qubit system, only the one native Y Y π/4 gate needs to be mitigated.The ideal Y Y π/4 gate can be decom-  2c and d, where the average gate fidelity of R Y Y ns π/4 is 0.9811 ± 0.0047 [44,45].Here, we assume the initial state preparation is perfect and apply the detection error correction method [46] for the measurement results, which makes state preparation and measurement (SPAM) errors negligible.The quasiprobabilities, i.e. the coefficients q i of the two-qubit Pauli The experimental results of fermionic dynamics are shown in Figs.2g and h for the cases without and with error mitigation, respectively.We use the initial two-qubits state of (|11 + |10 )/ √ 2, which contains the superposition of two occupied fermions and a single occupied fermion.The initial states simultaneously reveal the time evolutions of both single and double fermions, since no interference occurs between both dynamics.As shown in both Figs.2g and h, the state |11 and |00 do not change and there is an oscillation between the state |01 and |10 , which is the hallmark of fermion dynamics [32].We can see that without error mitigation, the contrast between the experimental populations (solid points) clearly decays compared with the numerical simulation under the Trotter step (hollow points) as shown in Fig. 2g, and the error-mitigated results are closer to the ideal numerical simulations as shown in Fig. 2h.The population fidelities depending on Trotter steps without and with the three error-mitigation schemes are shown in Fig. 2i.The fidelity of the gate without and with error mitigation are 0.9895±0.0021and 0.9975±0.0009,which are obtained by exponential fitting of data.
We note that the main improvement in fidelities comes from the PEC method.It is because our characterization of the gate is a good representation of the actual gate.We confirm it by simulating the dynamics of two spinless modes with the experimentally constructed PTM of R Y Y ns π/4 .The comparisons between the simulation with the noisy PTM and experimental data clearly validate the Pauli error assumption and accuracy of measured R Y Y ns π/4 (see Appendix B).We also note that MLE reduces the error bar of the results after PEC, which may compensate for the drawback of PEC in the aspect of variance.All the error bars are estimated by bootstrapping method (see Appendix A).

SIMULATION OF THREE SPINLESS FERMIONIC MODES
We implement the time evolution of three spinless fermionic modes as shown in Fig. 1a.The quantum circuit in Fig. 1d, obtained by the Trotter-Suzuki expansion, simulates the dynamics of the system.We initialize the three-qubit state to (|101 + |110 )/ √ 2, which includes two fermions in three sites.The fermions propagate in the system caused by the nearest-neighbor tunneling to other sites.The dynamical evolution would be restricted in the subspace spanned by {|110 , |101 , |011 } because of the conservation law of the number of fermions.
The dynamics can be different depending on the existence of nearest-neighbor interaction different from the previous case with two sites.In the experiment, we choose V = 0 and V = 2J for the scenarios without and with nearest-neighbor interactions, respectively.In Figs.3a, b, d, and e, dashed-dot lines show the results of ideal numerical simulation.When V = 0, the initial state |110 mostly evolves to the state |101 because of the Pauli exclusion principle of fermions.On the other hand, the initial state |101 can evolve to both |110 and |011 states.Therefore, in the beginning, the overall population of the state |110 reduces, and the population of the state |011 increases as shown in Figs.3a and b.When V = 2J, we can see the population of |110 is not rapidly decreasing and the state |101 is not populated compared to the case of V = 0.It is because the interaction is effectively attractive in our experiment.
In experiment, we set ϕ =  Appendix B).For the case of V = 0, the experimental results without and with error mitigation are shown in Figs.3a and b, respectively.We can see the error-mitigated results are more consistent with ideal numerical simulations.The fitted population fidelities per gate without and with error mitigation are 0.9792±0.0161and 0.9942±0.0067,as shown in Fig. 3c.We note that the main improvement comes from the PEC method.This is consistent with the comparison results between the simulation with the noisy PTMs of both Y Y gates and experimental data, which indicates the measured PTMs properly characterize the two-qubit gate errors (see Appendix C).Similar conclusions are obtained for the scenario with the nearest-neighbor interaction of V = 2J, as shown in Figs.3d, e and f.The fitted population fidelities per gate without and with error mitigation are 0.9827±0.0142and 0.9889 ± 0.0080.The two-qubit gate errors in the threequbit system are larger than that in the two-qubit system, which leads to the larger cost as C Y Y π/8 (1,2) = 1.157 and C Y Y π/8 (2,3) = 1.171.Therefore, we sample more random circuits as N s = 1500 to reduce the error bars of the final error-mitigated results.
We note that the fitted population fidelities after PEC could be unphysical because of the possible negative quasi-probabilities and the cost being over 1.As shown in Figs.3c and f, fidelities greater than 1 after applying PEC exist at the beginning of the three-fermion experiments.With the help of the MLE method based on the general properties of probabilities, we can make the population fidelities physical and reasonable, which is, below 1 as shown in green dots of Figs.3c and f.In addition, MLE reduces the error bars of the results after PEC, which may compensate for the drawback of PEC in the aspect of variance.After applying PEC, the standard deviation of the fourth Trotter step become unexpectedly large.However, the MLE method suppresses the standard deviation by a few factors in this case.

SIMULATION OF FOUR FERMIONIC MODES
We implement the simulation of four fermionic modes with spins by encoding qubits Q 1 Q 3 and Q 2 Q 4 as two fermionic sites, where each site can contain two fermionic modes with different spins, as shown in Fig. 1b.We initialize the four qubits state to (|1001 + |1010 )/ √ 2, which is corresponding to the superposition state between |↑ 1 |↓ 2 and |↑ 1 |↓ 1 as shown in Fig. 1b.The two fermions interact and exchange on the four modes caused by the nearest-neighbor tunneling and on-site interaction on subspace states {|0101 , |0110 , |1001 , |1010 } with conserved fermions and spins.The corresponding quantum circuit is illustrated in Fig. 1e.Here, we expect the different dynamic behaviors of spin and charge depending on the on-site interaction strength.Spin and charge of the fermions are defined as spin = n j,↑ − n j,↓ and charge = n j,↑ + n j,↓ with n j,↑(↓) being the fermion number of spin-up (down) on the j-th site.
To demonstrate the different behavior of spin and charge, we consider two scenarios U = 0 and U = 2J, which are related without and with on-site interactions, respectively.For the case of no on-site interaction, the spin and the charge of fermions should show the same dynamics.But if there is on-site interaction, the dynamics for spin and charge should show a different behavior, which is related to spin-charge separation for low energy excitation [35,[47][48][49].For the case of no onsite interaction, U = 0, the dynamics of fermions with spins are similar to those without spins, which are discussed in Figs.2g and h ing back and forth with the other state |↑ 2 |↓ 1 (|0110 ).Likewise, the state |↑ 1 |↓ 1 (|1010 ) is oscillating with the other state |↑ 2 |↓ 2 (|0101 ).This leads to equal amounts of charge and spin change at each site.However, when there exists on-site interaction, the dynamics of those fermionic states are no longer simple oscillations, and the state |↑ 1 |↓ 2 , which is not affected by the on-site interac-tion, and the state |↑ 1 |↓ 1 , which is strongly influenced by the on-site interaction, show completely different dynamics.
In experiment, we set ϕ = Jδt 2 = π 8 and apply M = 4 Trotter steps for both scenarios.With on-site interactions, each step consists of six Y Y ϕ gates and twenty single-qubit rotations.The whole quantum circuit con-sists of twenty-five two-qubit gates, where the additional one is for the initial state preparation, and eighty singlequbit gates.In order to apply error mitigations, in particular, the PEC scheme, we characterize four types of Y Y gates on different pairs of qubits, i.e., Y Y π/8 (1,2), Y Y π/8 (3,4), Y Y π/8 (1,3), and Y Y π/8 (2,4) For the case of U = 0, the experimental results without and with error mitigation are shown in Figs.4a and b, respectively.We can see the error-mitigated results are more consistent with the numerical simulation.The fitted population fidelities per gate without and with error mitigation are 0.9658±0.0426and 0.9993±0.0053,respectively, as shown in Fig. 4c.In Fig. 4d, the same behaviors of spin and charge, where the amounts of charge and spin changes are the same in the dynamics, are clearly shown with error mitigation, where the net charge and the net spin at site 1 are displayed.We can observe similar improvements in the dynamics with on-site (U = 2J) interaction, as shown in Figs.4e, and f.In Fig. 4g, the fitted population fidelities per gate without and with error mitigation are 0.9551 ± 0.0285 and 0.9925 ± 0.0081, respectively.With on-site interaction, the different behaviors of spin and charge are more clearly shown after error mitigations as shown in Fig. 4h.
We note that the PEC scheme alone does not recover the dynamics for the four fermionic modes as shown in Figs.4c and g different from the simulations with two and three fermionic modes.The main improvements in fidelities come from the post-selection methods based on the fermion-number-conservation assumption.One reason is that the experimental PTMs of two-qubit gates in the four-qubit system do not properly characterize the imperfections of those gates.It can be seen from the inconsistency of the numerically simulated data with the PTMs of two-qubit gates and experimental results (see Appendix C).The PTMs in our experiments cannot accommodate infidelities from cross-talks between qubits and time-correlated errors, that is, the performance of the gates during characterization can be different from those of actual gates in the quantum simulation.We note that coherent errors beyond the Pauli-error model assumption also are not captured by our experimental PTMs, which partially characterize errors related to Pauli errors.However, errors from cross-talks and time-correlations would be more important factors because coherent errors have no obvious reason to be seriously worse with system size.We also note that the number of sampling is greatly increased for applying the PEC method in the four-qubit system.With the measured fidelities, the costs increase to (2,4) = 1.223.The experimental Monte-Carlo sampling as N s = 2000 is not sufficient, which is limited by the time-consuming running of massive random circuits in physical hardware.

DISCUSSION AND OUTLOOK
We have applied the PEC method to improve the simulation of interacting fermions on up to four trapped ion qubits.We have observed the improvement of the simulation close to an order of magnitude in fidelities.With four qubits, we are able to observe the different dynamics of spin and charges for the interacting fermions.
Although PEC is the fundamental method in the experiment, several other error migration methods have been developed and applied.One of them is using the positive probability constraint.One intrinsic problem in the PEC method is that the probabilities of quantum states can be negative and not necessarily normalized.We impose constraints such that the probabilities of fermion states in the dynamics are positive and normalized, and determine the probabilities through the maximum-likelihood method.This method guarantees that observables are physically reasonable and can suppress statistical errors.
We also apply the symmetry constraint such that the total number of fermions and spins are conserved.In the simulation with four qubits, the majority of infidelity comes from the populations outside the numberconserving states.The constraints of total fermions conservation resolve the problem of population leakages.However, it is questionable if the number-conservation constraints would be scalable since the leakage can increase with the system and circuit depth.
From our experimental demonstration, although PEC incorporating constraints recoveries the dynamics of four trapped-ion qubits, it is clear that we need further improvements on the method for large-scale implementation.In particular, the drift of laser parameters and spectral crowding causes correlations resulting in remaining errors after PEC.As demonstrated in this work, combining with constraint methods can ease this issue.Recently, several proposals and experiments have been reported to tackle correlations such as using sparse Pauli-Lindblad models [23], learning-based error-mitigation [50], matrix product operator representation [51], etc.With all these techniques combined together, NISQ computers can potentially achieve the accuracy required for gaining the quantum advantage in practical applications, such as model inference for nuclear magnetic resonance (NMR) spectroscopy [52] and large-depth QAOA [53].To characterize our noisy two-qubit gates, a partial QPT with Pauli error assumption is applied with 15 measurements out of 144.To verify the assumption that the errors of our gates are consistent with the increasing of Trotter steps, we perform the numerical simulation of the extended Fermi-Hubbard model with the PTMs of the noisy gates.The simulation results are shown in Figs.9-11 for two spinless fermions, three spinless fermions, and two fermions with spins, respectively.For the simulation of two spinless fermions, with no additional leakage to other states out of the selected subspace, we can get a consistent simulation result with the experimental data, which further convinced the significant improvement of fidelity with the PEC method as shown in the introduction of Trotter errors Fig. 2(i).For the simulation results for three spinless fermions and two fermions with spins, the population in other states out of the ideal subspace is introduced because of the Trotter errors.The leakage will be further increased by the noisy gates in the experiment, which needs to be mitigated by MLE and PS methods.As a result, the simulation results shown in Fig. 10 and 11 have obvious deviations from the experimental data, which can be considered as a result of the unexpected errors of gates such as crosstalk error.

FIG. 1 .
FIG. 1.(a) Spinless and (b) spin-1 2 fermions moving in a one-dimensional lattice.The system dynamics can be unified by the extended Fermi-Hubbard model in Eq. (1), with the tunneling strength denoted by J and the on-site and the nearestneighbor interaction strengths by U and V , respectively.(c) Trapped-ion system with four ion-qubits.The ions form a 1-D crystal and each ion represents a qubit.The ion qubits can be manipulated by Raman laser beams, consisting of a global beam and individual beams, with the former covering the whole ion crystal and each of the latter focusing on a single ion.Mølmer-Sørensen gate Y Yϕ(m, n) ≡ exp (iϕσ y m ⊗ σy n ) can be implemented between any qubit pair (Qm, Qn) by modulating the phase of the individual beams.(d) Quantum circuit for simulating the dynamics of a three-site single-component extended Fermi-Hubbard model.With the Trotter-Suzuki expansion, the evolution is approximately implemented by M identical steps.(e) Quantum circuit for a two-site double-component Fermi-Hubbard model.(f) Quantum circuit compilation.While the two-qubit entangling gate provided by the native gate set of the trapped-ion system is Y Yϕ, XXϕ and ZZϕ can be constructed by Y Yϕ and appropriate single-qubit rotations, such as √ Z ≡ exp −i π 4 σz and √ X ≡ exp −i π 4 σx .The single-qubit gates (yellow squares) in (d, e) are the σz-rotations, i.e.Z2ϕ ≡ exp (−iϕσz).

FIG. 2 .
FIG. 2. (a) Quantum circuit of Trotter expansion for the simulation of dynamics of two spinless fermions.Two qubits are initialized to (|11 + |10 )/ √ 2 after the state preparation, then we apply M times of Trotter steps consisting of three two-qubit gates Y Y ns in each step, which is not perfect.Here, J and V represent the strength of nearest-neighbor tunneling and nearest-neighbor interaction.(b) Probabilistic error cancellation (PEC) scheme.The imperfections of the gate Y Y ns ϕ can be described by error channels acting on the ideal gate.The ideal gate Y Y ϕ can be realized by including the inverse of the error channels, which can be decomposed into combinations of basis operations with quasi-probabilities.Under the Pauli error assumption, there are 16 basis operations, where the list is shown in (c),(d),(e).(c) Pauli transfer matrix and (d) Deviation of the experimental Y Y ns ϕ gate, which is characterized through QPT.(e) Quasi-probabilities in the decomposition of the inverse error operations of experimental Y Y ns ϕ gate.(f ) Error-mitigated quantum circuits by the PEC.We implement the inverse error operations by sampling with the renormalized probabilities of quasi-probabilities as pi = |qi|/ 15 i=0 |qi|, where i indicates what number it corresponds to among basis operations.(g) Experimental data and (h) Error-mitigated data for V = 2J.The solid points represent the state populations measured from the experiment and the error-mitigated data after applying PEC, MLE, and PS methods.The dashed-dot lines represent ideal numerical simulation data with only Trotter errors, where the same colors correspond to the same experimental data points.The shadows around the data points represent error bars.(i) Population fidelity of experiment data and error-mitigated data.The fidelity is calculated by | k P k,ideal P k | 2 , which k represents the state basis.The population fidelities after each error-mitigation method are presented.All shadows in (g)-(i) represent the standard deviation of 1000 samples generated from the raw experimental data, using the bootstrapping method.

FIG. 3 .
FIG. 3. Dynamics of three spinless fermionic modes with the initial state of (|101 + |110 )/ √ 2 (a)-(c) for no interaction V = 0 and (d)-(f ) for V = 2J.(a),(d) Experimentally measured populations without error mitigation and (b),(e) with error mitigation.The solid points represent the experimentally-measured state populations and the dashed-dot lines represent the numerical simulation with the given Trotter steps.(c),(f ) The population fidelities after each error-mitigation method as PEC, MLE, and PS.All shadows in (a)-(f ) represent the standard deviation of 1000 samples generated from the raw experimental data, using the bootstrapping method.
FIG. 4. Dynamics of four fermionic modes mapped by two fermion sites and spins with the initial state of (|1001 + |1010 )/ √ 2 (a)-(d) for no interaction V = 0 and (e)-(h) for V = 2J.(a),(e) Experimentally measured population without error mitigations and (b),(f ) with error mitigation.The solid points represent the experimentally-measured state populations and the dashed-dot lines represent (c),(g) The population fidelities after each error-mitigation method as PEC, MLE, and PS.(d),(h) Fermionic spin and charge population.The hollow and filled symbols represent the experimental data without and with error mitigations, respectively.The dashed-dot lines represent the ideal numerical simulation with Trotter errors.All shadows represent the standard deviation of 1000 samples generated from the raw experimental data, using the bootstrapping method.

FIG. 5 .
FIG. 5. Fidelity distribution at the fourth Trotter step in the simulation of three spinless fermions of all the samples generated from the bootstrap method after (a) PEC method, (b) MLE method, and (c) PS method.The black lines in all the figures represent the original calculated results.Here each figure contains 1000 samples.

FIG. 9 .
FIG. 9. Experimental raw data compared with error-included simulation for two spinless fermions.We simulated the dynamics with noisy gates measured from QPT.(a) Populations of experimental raw data and simulation.The solid points represent the state populations measured from the experiment, and the dashed-dot lines represent numerical simulation data with noisy gates.(b) Population fidelity of experimental raw data and simulation data shown in (a).All the shadows in (a)-(b) represent the standard deviation of 1000 samples generated from the raw experimental data, using the bootstrapping method.

FIG. 10 .
FIG. 10.Experimental raw data compared with error-included simulation for three spinless fermions.We simulated the dynamics of V = 0 case with noisy gates measured from QPT.(a) Populations of experimental raw data and simulation.The solid points represent the state populations measured from the experiment, and the dashed-dot lines represent numerical simulation data with noisy gates.(b) Population fidelity of experimental raw data and simultion data shown in (a).(c)-(d) population and fidelity for V = 2J case, the representation is the same as (a)-(b).All the shadows in (a)-(d) represent the standard deviation of 1000 samples generated from the raw experimental data, using the bootstrapping method.

FIG. 11 .
FIG. 11.Experimental raw data compared with error-included simulation for two fermions with spins.We simulated the dynamics of U = 2J case with noisy gates measured from QPT.(a) Populations of experimental raw data and simulation.The solid points represent the state populations measured from the experiment, and the dashed-dot lines represent numerical simulation data with noisy gates.(b) Population fidelity of experimental raw data and simultion data shown in (a).(c)-(d) population and fidelity for U = 2J case, the representation is the same as (a)-(b).All the shadows in (a)-(d) represent the standard deviation of 1000 samples generated from the raw experimental data, using the bootstrapping method.
4 under Pauli error assumption, which is visualized as Fig. 2b.The experimental Y Y ns π/4 gate is characterized by QPT.The experimental results of the PTMs of the Y Y ns π/4 gate and the deviation from the ideal one obtained by R Y Y ns π/4 −R Y Y π/4 are shown in Figs.
Trotter steps for both scenarios.With nearest-neighbor interactions, each Trotter step consists of six Y Y ϕ gates and twenty single-qubit rotations.The whole quantum circuit consists of twentyfive two-qubit gates, where the additional one is for the initial state preparation, and eighty single-qubit gates.Then, we apply the error-mitigation methods to the dynamic evolution of both scenarios.There are two types of Y Y gates on different pairs of qubits, i.e., Y Y π/8 (1, 2) and Y Y π/8 (2, 3), both of which require characterization and inverse-error decomposition.In the PEC scheme, both Y Y π/8 (1, 2) and Y Y π/8 (2, 3) gates are individually characterized by QPT and the inverse errors are decomposed independently.The average gate fidelities are 0.9779 ± 0.0088 for R Y Y ns π/8 (1,2) and 0.9748 ± 0.0085 for R Y Y ns π/8 (2,3) (see