Hamiltonian quantum simulation with boundedstrength controls
Abstract
We propose dynamical control schemes for Hamiltonian simulation in manybody quantum systems that avoid instantaneous control operations and rely solely on realistic boundedstrength control Hamiltonians. Each simulation protocol consists of periodic repetitions of a basic control block, constructed as a suitable modification of an “Eulerian decoupling cycle,” that would otherwise implement a trivial (zero) target Hamiltonian. For an open quantum system coupled to an uncontrollable environment, our approach may be employed to engineer an effective evolution that simulates a target Hamiltonian on the system, while suppressing unwanted decoherence to the leading order. We present illustrative applications to both closed and opensystem simulation settings, with emphasis on simulation of nonlocal (twobody) Hamiltonians using only local (onebody) controls. In particular, we provide simulation schemes applicable to Heisenbergcoupled spin chains exposed to general linear decoherence, and show how to simulate Kitaev’s honeycomb lattice Hamiltonian starting from Isingcoupled qubits, as potentially relevant to the dynamical generation of a topologically protected quantum memory. Additional implications for quantum information processing are discussed.
pacs:
03.67.Lx, 03.65.Fd, 03.67.amitctp 4504
1 Introduction
The ability to accurately engineer the Hamiltonian of complex quantum systems is both a fundamental control task and a prerequisite for quantum simulation, as originally envisioned by Feynman [1, 2, 3]. The basic idea underlying Hamiltonian simulation is to use an available quantum system, together with available (classical or quantum) control resources, to emulate the dynamical evolution that would have occurred under a different, desired Hamiltonian not directly accessible to implementation [4]. From a controltheory standpoint, the simplest setting is provided by openloop Hamiltonian engineering in the time domain [5, 6], whereby coherent control over the system of interest is achieved solely based on suitably designed timedependent modulation (most commonly sequences of control pulses), without access to ancillary quantum resources and/or measurement and feedback. While openloop Hamiltonian engineering techniques have their origin and a long tradition in nuclear magnetic resonance (NMR) [8, 7], the underlying physical principles of “coherent averaging” have recently found widespread use in the context of quantum information processing (QIP), leading in particular to dynamical symmetrization and dynamical decoupling (DD) schemes for control and decoherence suppression in open quantum systems [9, 10, 11, 12, 13, 14].
As applications for quantum simulators continue to emerge across a vast array of problems in physics and chemistry, and implementations become closer to experimental reality [3, 15, 16], it becomes imperative to expand the repertoire of available Hamiltonian simulation procedures, while scrutinizing the validity of the relevant control assumptions. With a few exceptions (notably, the use of socalled “perturbation theory gadgets” [17]), openloop Hamiltonian simulation schemes have largely relied thus far on the ability to implement sequences of effectively instantaneous, “bangbang” (BB) control pulses [18, 19, 20, 21, 22, 23, 24, 25]. While this is a convenient and often reasonable first approximation, instantaneous pulses necessarily involve unbounded control amplitude and/or power, something which is out of reach for many control devices of interest and is fundamentally unphysical. In the context of DD, a general approach for achieving (to at least the leading order) the same dynamical symmetrization as in the BB limit was proposed in [26], based on the idea of continuously applying boundedstrength control Hamiltonians according to an Eulerian cycle, socalled Eulerian DD (EDD). From a Hamiltonian engineering perspective, EDD protocols translate directly into boundedstrength simulation schemes for specific effective Hamiltonians – most commonly, the trivial (zero) Hamiltonian in the case of “nonselective averaging” for quantum memory (or “timesuspension” in NMR terminology). More recently, EDD has also served as the starting point for boundedstrength gate simulation schemes in the presence of decoherence, socalled dynamically corrected gates (DCGs) for universal quantum computation [27, 28, 29, 30].
In this work, we show that the approach of Eulerian control can be further systematically exploited to construct boundedstrength Hamiltonian simulation schemes for a broad class of target evolutions on both closed and open (finitedimensional) quantum systems. Our techniques are deviceindependent and broadly applicable, thus substantially expanding the control toolbox for programming complex Hamiltonians into existing or nearterm quantum simulators subject to realistic control assumptions.
The content is organized as follows. We begin in Sect. II by introducing the appropriate controltheoretic framework and by reviewing the basic principles underlying openloop simulation via average Hamiltonian theory, along with its application to Hamiltonian simulation in the BB setting. Sect. III is devoted to constructing and analyzing simulation schemes that employ boundedstrength controls: while Sec. III.A reviews required background material on EDD, Sec. III.B introduces our new Eulerian simulation protocols for a generic closed quantum system. In Sec. III.C we separately address the important problem of Hamiltonian simulation in the presence of slowlycorrelated (nonMarkovian) decoherence, identifying conditions under which a desired Hamiltonian may be enacted on the target system while simultaneously decoupling the latter from its environment, and making further contact with DCG protocols. Sect. IV presents a number of illustrative applications of our general simulation schemes in interacting multiqubit networks. In particular, we provide explicit protocols to simulate a large family of twobody Hamiltonians in Heisenbergcoupled spin systems additionally exposed to depolarization or dephasing, as well as to achieve Kitaev’s honeycomb lattice Hamiltonian starting from Isingcoupled qubits. In all cases, only local (singlequbit, possibly collective) control Hamiltonians with bounded strength are employed. A brief summary and outlook conclude in Sec. V.
2 Principles of Hamiltonian simulation
2.1 Controltheoretic framework
We consider a quantum system , with associated Hilbert space , whose evolution is described by a timeindependent Hamiltonian . As mentioned, Hamiltonian simulation is the task of making evolve under some other timeindependent target Hamiltonian, say, . Without loss of generality, both the input and the target Hamiltonians may be taken to be traceless. Two related scenarios are worth distinguishing for QIP purposes:
Closedsystem simulation, in which case coincides with the quantum system of interest, (also referred to as the “target” henceforth), which undergoes purely unitary (coherent) dynamics;
Opensystem simulation, in which case is a bipartite system on , where represents an uncontrollable environment (also referred to as bath henceforth), and the reduced dynamics of the target system is nonunitary in general.
In both cases, we shall assume the target system to be a network of interacting qudits, hence , for finite and . In the general opensystem scenario, the joint Hamiltonian on may be expressed in the following form,
(1) 
where the operators ( and () act on () respectively, and all the bath operators are assumed to be normbounded, but otherwise unspecified (potentially unknown). The closedsystem setting is recovered from Eq. (1) in the limit . Likewise, we may express the target Hamiltonian in a similar form, with two simulation tasks being of special relevance: , in which case the objective is to realize a desired system Hamiltonian while dynamically decoupling from its bath , thereby suppressing unwanted decoherence [11]; or, more generally, and , where the simulated, dynamically symmetrized error generators may allow for decoherencefree subspaces or subsystems to exist [13, 31].
Formally, the dynamics is modified by an openloop controller acting on the target system according to
(2) 
where the operators and the (real) functions represent the available control Hamiltonians and the corresponding, generally timedependent, control inputs respectively. Clearly, if the Hamiltonian is contained in the admissible control set, the corresponding control problem is trivial and the desired timeevolution,
can be exactly simulated continuously in time. However, this level of control need not be available in settings of interest, including open quantum systems where control actions are necessarily restricted to the target system alone, in Eq. (2). Following the general idea of “analog” quantum simulation [3], we shall assume in what follows a restricted set of control Hamiltonians (in a sense to be made more precise later) and focus on the task of approximately simulating the desired time evolution at a final time , or more generally, stroboscopically in time, that is, at instants , where
and is a fixed minimum time interval. Choosing sufficiently small allows in principle any desired accuracy in the approximation to be met, with the limit formally recovering the continuous limit.
Specifically, let and denote the unitary propagators associated to the total and the control Hamiltonians in Eq. (2), respectively:
(3)  
(4) 
where we have set and indicates timeordering, as usual. Then, for a given pair , we shall provide sufficient conditions for to be “reachable” from and, if so, devise a cyclic control procedure such that the resulting controlled propagator
(5) 
where is the cycle time of the controller, that is, . In general, we shall allow for to differ from , corresponding to an overall scale factor in the simulated time, as it will become apparent later. If, for a fixed input Hamiltonian , arbitrary target Hamiltonians are reachable for given control resources, the simulation scheme is referred to as universal. In this case, complete controllability must be ensured by the tunable Hamiltonians in conjunction with the “drift” [6]. In contrast, we shall be especially interested in situations where control over is more limited.
Similar to DD protocols, Hamiltonian simulation protocols are most easily constructed and analyzed by effecting a transformation to the “toggling” frame associated to in Eq. (4) [7, 11, 14]. That is, evolution in the toggling frame is generated by the timedependent, controlmodulated Hamiltonian
(6) 
with the corresponding togglingframe propagator being related to the physical propagator in Eq. (3) by . Since the control propagator is cyclic and is timeindependent, it follows that and, furthermore, acquires the periodicity of the controller, . Thus, the stroboscopic controlled dynamics of the system is determined by
(7) 
Average Hamiltonian theory [7, 35] may then be invoked to associate an effective timeindependent Hamiltonian to the evolution in the togglingframe:
(8) 
where is determined by the Magnus expansion [32], Explicitly, the leadingorder term, determining evolution over a cycle up to the first order in time, is given by
(9) 
with (absolute) convergence being ensured as long as [34]. Subject to convergence condition, higherorder corrections for evolution over time can also be upperbounded by (see Lemma 4 in [33])
(10) 
Ideally, one would like to achieve , so that equality would hold in Eq. (5) for all . In what follow, we shall primarily focus on achieving firstorder simulation instead, by engineering the control propagator in such a way that
(11) 
whereby, using Eq. (10) with ,
(12) 
In general, the accuracy of the approximation in Eq. (11) improves as the “fast control limit”, , is approached. Physically, this corresponds to requiring that the shortest control time scale (pulse separation) involved in the control sequence be sufficiently small relative to the shortest correlation time of the dynamics induced by [35, 36]. While the problem of constructing general highorder Hamiltonian simulation schemes is of separate interest, we stress that secondorder simulation can be readily achieved, in principle, by ensuring that is timesymmetric, namely, for . Since all oddorder Magnus corrections vanish in this case [35], it follows (by using again Eq. (10), with ), that , correspondingly boosting the accuracy of the simulation.
2.2 Hamiltonian simulation with bangbang controls
BB Hamiltonian simulation provides the simplest control setting for achieving the intended objective, given in Eq. (5). Two main assumptions are involved: (i) First, we must be able to express the target Hamiltonian as
(13) 
where are unitary operators on and the nonnegative real numbers (not all zero). (ii) Second, the available control resources include a discrete set of instantaneous pulses on , whose application results in a piecewiseconstant control propagator over , with corresponding togglingframe propagators , , [9, 14]. Assumptions (i)(ii) together allow for the timeaverage in Eq. (9) to be mapped to a convex (positiveweighted) sum. Eq. (13) may be interpreted as a sufficient condition for the target Hamiltonian to be reachable from given openloop unitary control on alone. Reachable Hamiltonians must thus be at least as “disordered” as the input one in the sense of majorization [21, 14].
Specifically, Eq. (13) leads naturally to the following BB simulation scheme. Given simulation weights , define the following simulation intervals and timing pattern:
(14) 
A piecewiseconstant control propagator for the basic simulation block to be repeated may then be constructed as follows:
(15) 
By using Eq. (9), it is immediate to verify that
(16) 
resulting in the desired controlled evolution, Eqs. (11)(12), provided that the convergence conditions for firstorder simulation under are obeyed. Since, in practice, technological limitations always constrain the cycle duration to a finite minimum value , such conditions ultimately determine the maximum simulated time up to which evolution under may be reliably simulated using the physical Hamiltonian .
In analogy with BB DD schemes, realizing the prescription of Eq. (15) requires to discontinuously change the control propagator from to , via an instantaneous BB pulse at the th endpoint . As a result, despite its conceptual simplicity, BB simulation is unrealistic whenever large control amplitudes are not an option, and the evolution induced by during the application of a control pulse must be considered from the outset. This demands redesigning the basic control block in such a way that the actions of and are simultaneously accounted for.
3 Hamiltonian simulation with bounded controls
3.1 Eulerian simulation of the trivial Hamiltonian
The key to overcome the disadvantages of BB Hamiltonian simulation is to ensure that the control propagator varies smoothly (continuously) in time during each control cycle. We achieve this goal by relying on Eulerian control design [26]. To introduce the necessary grouptheoretical background, we begin by revisiting how, for the special case of a target identity evolution (that is, , also corresponding to a “noop” gate, in terms of the endtime simulated propagator), EDD can be naturally interpreted as a boundedstrength simulation scheme.
In the Eulerian approach, the available control resources include a discrete set of unitary operations on , say, , , which are realized over a finite time interval through application of boundedstrength control Hamiltonians , with . That is,
(17) 
Note that the choice of the control Hamiltonians is not unique, which allows for implementation flexibility. The unitaries are identified with the image of a generating set of a finite group under a faithful, unitary, projective representation [26]. That is, let be a finite group of order , such that each element may be written as an ordered product of elements in a generating set of order , be the representation map [37], and . The Cayley graph of relative to can be thought of as pictorially representing all elements of as strings of generators in . Each vertex represents a group element and a vertex is connected to another vertex by a directed edge “colored” (labeled) with generator if and only if . The number of edges in is thus equal to . Because a Cayley graph is regular, it always has an Eulerian cycle that visits each edge exactly once and starts (and ends) on the same vertex [38, 39]. Let us denote with the ordered list of generators defining an Eulerian cycle on which, without loss of generality, starts (and ends) at the identity element of .
Once a control Hamiltonian for implementing each generator as in Eq. (17) is chosen, an EDD protocol is constructed by assigning a cycle time as and by applying the control Hamiltonians sequentially in time, following the order determined by the Eulerian cycle . Thus,
(18) 
where is the image of the generator labeling the th edge in . As established in [26], the lowestorder average Hamiltonian associated to the above EDD cycle has the form , where for any operator acting on , the map
(19) 
projects onto the centralizer of (i.e., commutes with all ), and
(20) 
implements an average of over both the control interval and the group generators. Accordingly, boundedstrength simulation of is achieved provided that the following DD condition is obeyed:
(21) 
By Schur’s lemma, this is automatically ensured if the group representation acts irreducibly on . Formally, the BB limit may be recovered by letting for all [26], reflecting the ability to directly implement all the group elements (with no overhead, as if ) and corresponding to uniform simulation weights .
3.2 Eulerian simulation protocols beyond noop: Construction
We show how the Eulerian cycle method can be extended to boundedstrength simulation of a nontrivial class of target Hamiltonians. We assume that may be expressed as a convex unitary mixture of the group representatives ,
(22) 
We construct the desired control protocol starting from an Eulerian cycle on . Specifically, the idea is to append to each of the control slots that define an EDD scheme a freeevolution (or “coasting”) period of suitable duration, in such a way that the net simulated Hamiltonian is modified from to as given in Eq. (22). A pictorial representation of the basic control block is given in Fig. 1. As in Eq. (17), let denote the minimum time duration required to implement each generator, hence, to smoothly change the control propagator from a value to along the cycle. While such “ramping up” control intervals have all the same length, each “coasting” interval is designed to keep the control propagator constant at for a duration determined by the corresponding weight . Since the control is switched off during coasting, continuity of the overall control Hamiltonian may be ensured, if desired, by requiring that
(23) 
in addition to the boundedstrength constraint.
An Eulerian simulation protocol may be formally specified as follows. As before, let the th time interval be denoted as , , with and defining the cycle time . For each , let as in the BB case. The duration of the th coasting period is then assigned as
(24) 
resulting in the following timing pattern [compare to Eq. (14)]:
(25) 
As the expression for the cycle times makes it clear, the resulting protocol may be equivalently interpreted in two ways: starting from an EDD cycle, corresponding to and , we introduce the coasting periods to allow for nontrivial simulated dynamics to emerge; or, starting from a BB simulation scheme for , corresponding to , we introduce the rampingup periods to allow for control Hamiltonians to be smoothly switched over . Either way, boundedstrength protocols imply a timeoverhead relative to the BB case, recovering the BB limit as as expected. Explicitly, the control propagator for Eulerian simulation has the form:
(26)  
(27) 
The resulting firstorder Hamiltonian under Eulerian simulation is derived by evaluating the timeaverage in Eq. (9) with the control propagator given by Eqs. (26)(27). Direct calculation along the lines of [26] yields:
where the last equality follows from two basic properties of Eulerian cycles: firstly, the list (and also ) of the vertices that are being visited contains each element precisely times; secondly, in traversing the Cayley graph, each group element is left exactly once by a labeled edge for each generator . Thus, by recalling the definitions given in Eqs. (19) and (20), we finally obtain
(28) 
which indeed achieves the intended firstorder simulation goal, Eqs. (11)(12), as long as convergence holds and the DD condition of Eq. (21) is obeyed.
The simulation accuracy may be improved by symmetrizing in time. In analogy to symmetrized EDD protocols [9], this can be easily accomplished by running the protocol and then suitably running it again in reverse. Specifically, let the duration of the coasting interval be changed as . Run the protocol as described above until time . Then, from time until time , modify Eqs. (26)(27) as follows:
for . Provided that one is able to implement , we again obtain
while satisfying for , hence ensuring that .
3.3 Eulerian simulation while decoupling from an environment
The ability to implement a desired Hamiltonian on the target system , while switching off (at least to the leading order) the coupling to an uncontrollable environment , is highly relevant to realistic applications. That is, with reference to Eq. (1), the objective is now to simultaneously achieve , by unitary control operations acting on alone. Because the firstorder Magnus term is additive [recall Eq. (9)], it is appropriate to treat each summand of individually, leading to a relevant average Hamiltonian of the form
where for a generic operator on we let
We can then apply the analysis of Sec. 3.2 to the internal system Hamiltonian () and each error generator () separately, to obtain in both cases a simulated operator of the form given in Eq. (28):
Since the task is to decouple from while maintaining the nontrivial evolution due to , the reachability condition of Eq. (22) must now ensure that
(29)  
(30) 
Accordingly, it is necessary to extend the DD assumption of Eq. (21) to become
(31)  
(32) 
such that holds for each of the summands in . Altogether we recover
It is interesting in this context to highlight some similarities and differences with DCGs [27], which also use Eulerian control as their starting point and are specifically designed to achieve a desired unitary evolution on the target system while simultaneously removing decoherence to the leading [27, 28, 30] or, in principle, arbitrarily high order [29]. By construction, the opensystem simulation procedure just described does provide a firstorder DCG implementation for the target gate : in particular, the requirement that Eqs. (29)(30) be obeyed together (for the same weights ) is effectively equivalent to evading the “nogo theorem” for blackbox DCG constructions established in [28], with the coasting intervals and the resulting “augmented” Cayley graph playing a role similar in spirit to a (firstorder) “balancepair” implementation. Despite these formal similarities, a number of differences exist between the two approaches: first, an obvious yet important difference is that DCG constructions focus directly on synthesizing a desired unitary propagator, as opposed to a desired Hamiltonian generator. Second, while the internal system Hamiltonian, , is a crucial input in a Hamiltonian simulation problem, it is effectively treated as an unwanted error contribution in analytical DCG constructions, in which case complete controllability over the target system must be supplied by the controls alone. Although in more general (optimalcontrol inspired) DCG constructions [30], limited external control is assumed and may become essential for universality to be maintained, emphasis remains, as noted above, on endtime synthesis of a target propagator. Finally, a main intended application of DCGs is realizing lowerror single and twoqubit gates for use within faulttolerant quantum computing architectures, as opposed to robust Hamiltonian engineering for manybody quantum simulators which is our focus here.
3.4 Eulerian simulation protocols: Requirements
Before presenting explicit illustrative applications, we summarize and critically assess the various requirements that should be obeyed for Eulerian simulation to achieve the intended control objective of Eq. (5) in a closed or, respectively, opensystem setting:

Time independence. Both the internal Hamiltonian and the target Hamiltonian are taken to be timeindependent (and, without loss of generality, traceless).

Reachability. The target Hamiltonian must be reachable from , that is, there must be a control group , with a faithful, unitary projective representation mapping , such that Eq. (22) holds. For dynamicallycorrected Eulerian simulation in the presence of an environment, this requires, as noted, that for the same weights , the desired system Hamiltonian is reachable from while the trivial (zero) Hamiltonian is reachable from each error generator separately, such that both Eqs. (29)(30) hold together.

Bounded control. For each generator of the chosen control group , we need access to bounded control Hamiltonians , such that application of over a time interval of duration realizes the group representative , additionally subject (if desired) to the continuity condition of Eq. (23).

Timeefficiency. If the choice of is not unique for given , the smallest group should be chosen, in order to keep the number of intervals per cycle, , to a minimum. In particular, efficient Hamiltonian simulation requires that (hence also ) scales (at most) polynomially with the number of subsystems .
The key simplification that the timeindependence Assumption (1) introduces into the problem is that the periodicity of the control action is directly transferred to the togglingframe Hamiltonian of Eq. (6), allowing one to simply focus on singlecycle evolution. Although this assumption is strictly not fundamental, general timedependent Hamiltonians may need to be dealt with on a casebycase basis (see also [40, 41, 42]). A situation of special practical relevance arises in this context for open systems exposed to classical noise, in which case and the systembath interaction in Eq. (1) is effectively replaced by a classical, timedependent stochastic field. Similar to DD and DCG schemes, Eulerian simulation protocols remain applicable as long as the noise process is stationary and exhibiting correlations over sufficiently long time scales [9, 43].
The reachability Assumption (2) is a prerequisite for Eulerian Hamiltonian simulation schemes. Although BB Hamiltonian simulation need not be groupbased, most BB schemes follow this design principle alike. Assumption (3), restricting the admissible control resources to physical Hamiltonians with bounded amplitude (thus finite control durations, as opposed to instantaneous implementation of arbitrary group unitaries as in the BB case) is a basic assumption of the Eulerian control approach. As remarked, our premise is that the available Hamiltonian control is limited, restricted to only the target system if the latter is coupled to an environment, and typically nonuniversal on ; in particular, we cannot directly express and apply , or else the problem would be trivial. In addition to errorcorrected Hamiltonian simulation in open quantum systems, scenarios of great practical interest may arise when the control Hamiltonians are subject to more restrictive locality constraints than the system and target Hamiltonians are (e.g., twobody simulation with only local controls, see also Sec. 4.1).
The required decoupling conditions in Assumption (4) are automatically obeyed if the representation acts irreducibly on . This follows directly from Schur’s lemma, together with the fact that the map defined in Eq. (20) is tracepreserving, and both and can be taken to be traceless. While convenient, irreducibility is not, however, a requirement. When the representation is reducible, care must be taken in order to ensure that Assumption (4) is nevertheless obeyed. It should be stressed that this is possible independently of the target Hamiltonian . Therefore, if the choice () works for one Eulerian simulation scheme (whether is irreducible or not), then it can be used for Eulerian simulation with any target that belongs to the reachable set from , that is, that can satisfy Eq. (22).
We close this discussion by recalling that it is always possible for a finitedimensional target system to find a control group for which both Assumptions (2) and (4) are satisfied, by resorting to the concept of a transformer [22, 14]. A transformer is a pair , where is a finite group and is a faithful, unitary, projective representation such that, for any traceless Hermitian operators and on with , one may express
We illustrate this general idea in the simplest case of a single qubit, . Let denote the Pauli matrices and the unitary matrix defined by
(33) 
which corresponds to a rotation by an angle about an axis . Direct calculation shows that and that conjugation by cyclically shifts the Paulimatrices, i.e., , and . Consider now the group given by the presentation
Using the defining relations of this group, its elements can always be written as , where and . Clearly, the assignment given by yields a faithful, unitary, irreducible representation since the Pauli matrices commute up to phase. It is shown in [22] that the pair defines a transformer in the sense given above, namely, any traceless matrix may be reached from any fixed traceless, nonzero matrix , for suitable nonnegative weights . The irreducibility property for any transformer pair can be easily established by contradiction [44].
A drawback of the transformer formalism is that general transformer groups tend to be large, making purely transformerbased simulation schemes inefficient. In practice, given the native system Hamiltonian , the challenge is to find a group that grants a reasonably efficient scheme while satisfying Assumptions (2) and (4), and subject to the ability to implement the required control operations. As we shall see next, transformerinspired ideas may still prove useful in devising simulation schemes in the presence, for instance, of additional symmetry conditions.
4 Illustrative applications
In this section, we explicitly analyze simple yet paradigmatic Hamiltonian simulation tasks motivated by QIP applications. While a number of other interesting examples and generalizations may be envisioned (as also further discussed in the Conclusions), our goal here is to give a concrete sense of the usefulness and versatility of our Eulerian simulation approach in physically realistic control settings. In particular, we focus on achieving nonlocal Hamiltonian simulation using only boundedstrength local (singlequbit) control, in both closed and open multiqubit systems.
4.1 Eulerian simulation in closed Heisenbergcoupled qubit networks
Let us start from the simplest case of a system consisting of qubits, interacting via an isotropic Heisenberg Hamiltonian of the form
where has units of energy and the second equality defines an equivalent compact notation. We are interested in a class of target XYZ Hamiltonians of the form
(34) 
For instance, , corresponds to an isotropic XX model, whereas if with , an XXZ interaction is obtained, the special value corresponding to the important case of a dipolar Hamiltonian. The construction of a simulation protocol starts from observing that Hamiltonians as in Eq. (34) are reachable from , in the sense of Eq. (22), based on singlequbit control only.
Specifically, let , and let the representation map to . That is, is mapped to the following set of unitaries:
(35) 
Choosing the generators of to be and , we assume that we have access to the control Hamiltonians
where the control inputs and satisfy and , for . Recalling Eq. (17), this yields the control propagators
with and (up to phase), as desired.
Note that for any singlequbit Hamiltonians and , averaging over the unitary group in Eq. (35) results in the following projection superoperator:
(36) 
In general, the map is tracepreserving and, in this case, it acts nontrivially only on the first qubit. Thus, is tracepreserving on the first qubit. Since each term in is traceless in the first qubit, the decoupling condition follows directly from Eq. (36), even though the relevant representation is, manifestly, reducible.
Having satisfied our main requirements for Eulerian simulation, reachability of XYZ Hamiltonians as in Eq. (34) is equivalent to the existence of a solution to the following set of conditions:
(37)  
for nonnegative weights . While infinitely many choices exist in general, minimizing the total weight keeps the simulation time overhead to a minimum. For instance, it is easy to verify that a dipolar Hamiltonian of the form
may be simulated with minimum time overhead by choosing weights
The Cayley graph associated with the resulting Eulerian simulation protocol is depicted in Fig. 2, with the explicit timing structure of the control block as in Fig. 1 and control segments per block. It is worth observing that although the weights and are zero in the particular case at hand, all group members of are nonetheless required, and the unitaries and still show up in the simulation scheme (during the rampingup subintervals, as evident from Eq. (26)). This is crucial to guarantee that the unwanted term is projected out.
The above analysis and simulation protocols can be easily generalized to a chain of qubits (or spins), subject to nearestneighbor (NN) homogeneous Heisenberg couplings, that is, a Hamiltonian of the form
where for later reference we have introduced the standard compact notation and we assume for concreteness that is even. In this case, we need only change the unitary representation of to be defined by the two generators and , resulting in the set of unitaries [42]
Physically, the required generators and correspond to control Hamiltonians that are still just sums of 1local terms, and that act nontrivially on odd qubits only:
We expect that the design of Eulerian simulation schemes for more general scenarios where both the input and the target are arbitrary twobody Hamiltonians (including, for instance, longrange couplings) will greatly benefit from the existence of combinatorial approaches for constructing efficient DD groups [45, 41]. A more indepth analysis of this topic is, however, beyond our current scope.
4.2 Errorcorrected Eulerian simulation in open Heisenbergcoupled qubit networks
Imagine now that the Heisenbergcoupled system considered in the previous section is coupled to an environment , and the task is to achieve the desired XYZ Hamiltonian simulation while also removing arbitrary linear decoherence to the leading order. The total input Hamiltonian has the form
(38) 
where and , for each and , are operators acting on , whose norm is sufficiently small to ensure convergence of the relevant Magnus series, similar to firstorder DCG constructions [27, 28]. The target Hamiltonian then reads
in terms of suitable couplingstrength parameters as in Eq. (34). As before, we start by analyzing the case of qubits in full detail. Our strategy to synthesize a dynamically corrected simulation scheme involves two stages: (i) We will first decouple from