Many-body state engineering using measurements and fixed unitary dynamics

We develop a scheme to prepare a desired state or subspace in high-dimensional Hilbert-spaces using repeated applications of a single static projection operator onto the desired target and fixed unitary dynamics. Benchmarks against other control schemes, performed on generic Hamiltonians and on Bose-Hubbard systems, establish the competitiveness of the method. As a concrete application of the control of mesoscopic atomic samples in optical lattices we demonstrate the near deterministic preparation of Schr\"{o}dinger cat states of all atoms residing on either the odd or the even sites.

Introduction -Quantum state engineering aims at preparing a desired target quantum state [1]. When the system Hamiltonian is sufficiently controllable, tailored control pulses can produce every state via unitary control [2]. Paradigmatic proofs of controllability rely on the decomposition of the desired unitary transformation into two-level unitaries [3] [ Fig. 1(a)], or on two sufficiently non-commuting Hamiltonians that are switched in a "bang-bang" fashion [4] [ Fig. 1(b)]. Although tailored applications of optimal control to many-body systems such as ultacold atoms in optical lattices have lately been applied to specific applications such as the superfluid to Mott-insulator phase-transition via control of the lattice depth [5,6], the implementation of generic approaches to full controllability remains impractical [7].
Besides unitary control via the system Hamiltonian, quantum state engineering can also be exerted for a fixed Hamiltonian via the back-action induced by measurements [8][9][10][11][12][13][14]. On the one hand, measurement operators that are adapted on previous measurement outcomes can steer a quantum state into a target state [8,9]. For example, spin correlations [15] can be induced by a sequence of standing wave probes [16]. The implementation of measurements for such adaptive schemes is, however, impractical for high-dimensional systems with natural experimental restrictions.
On the other hand, given a projector onto the desired stateP t = |ψ t ψ t |, an initial state |ψ 1 can be steered into the target: Control-free control [13] uses multiple evenly distributed observables (MEDO), i.e. the system is projected step-wise onto states that are evenly distributed between the initial and the target state [ Fig. 1(c)]. Alternatively, if -besides the projector onto the target stateP t = |ψ t ψ t | -only one auxiliary pro-jectorP U can be used, it is advisable to chooseP U to be mutually unbiased with respect toP t to maximise the probability to find the target state, i.e. to perform mutually unbiased measurements (MUM) [ Fig. 1(d)] [17]. These schemes assume no unitary evolution between subsequent measurements.
Here, we propose a hybrid of unitary and controlfree control by investigating quantum state engineering via optimised fixed unitary evolution and measurements (FUMES), which imposes fairly weak requirements on the experimental infrastructure: The system Hamilto-nianĤ of finite dimension d is fixed and uncontrollable. The only other means to steer the system is provided by a unique fixed measurement operator. The measurement operator signals whether the target subspace is reached (eigenstate ofP t , "success"), or into which other subspace the state was projected ("failure"). A crucial feature of this protocol is the ability to model precisely the unitary many-body evolution and to determine the precise times that maximise the success probability of the subsequent measurement. Even provided with these minimal tools, desired states can be produced under mild assumptions on the measurement operator and the system Hamiltonian. To benchmark the procedure, FUMES applied with a randomly chosen Hamiltonian is first compared to MUM and MEDO. Further below, we apply it to Bose-Hubbard multi-mode systems, which also permit benchmarking against unitary bang-bang control. Finally, as a practical example of FUMES, we demonstrate the nearly deterministic preparation of macroscopic superposition states. Quantum state engineering by FUMES -We denote the d eigenstates of the constant Hamiltonian H by |ψ 1 , . . . , |ψ d . The goal is to prepare an eigenstate of a given K-dimensional the aim is to prepare a state orthogonal to some |φ ⊥ . The system is initially in the ground state, |ψ 1 .
The application of the projector at a random moment in time has low success probability [18]. We therefore optimise numerically the moment in time at which a measurement is performed, maximising the probability to populate a desired eigenstate ofP t . A measurement that does not yield the desired outcome described byP t has failed and projects the state into one of the M different possible outcomes that indicate failure (1 ≤ M ≤ d−K). The procedure is repeated until successful, i.e. for each failed measurement, a new optimal waiting time is chosen before the next measurement is applied.
For a binary measurement (M = 1), failure leads to state collapse onto an eigenstate of the d−K-dimensional projectorQ = 1 −P t , and only marginal information is gained. A fully granular measurement with M = d − K reveals the precise state |η m that the system collapsed onto.
In general, a target can be reached by FUMES if and only if no conserved quantity can be constructed out of the failed projectorsQ j , i.e. every sum of m projectors (1 ≤ m ≤ M ) does not commute withĤ, for any set of {j l ...j k }. In this case, every state that the system can be projected on evolves into the target subspace. Alternatively, if a combination of theQ j com-mutes withĤ, states in the corresponding eigenspace will never evolve into the target subspace.
Benchmarking -To prove the general applicability of FUMES, we benchmark it against MEDO and MUM with a random target state |φ t that is prepared using a random Hamiltonian sampled from the Gaussian unitary ensemble, which ensures uniform distribution in the space of Hamiltonians with a particular dynamical time-scale [19]. Our figure of merit is the probability to prepare the target state |φ t after r measurements, p(r) [28].
For MEDO [13] [ Fig. 1(c)], we assume a set of r operators, each of which projects onto a state where |φ ⊥ t is the component of |ψ 1 orthogonal to |φ t , θ is the angle between the initial and the target state with |ψ 1 = cos (θ) |φ t +sin (θ) |φ ⊥ t and j ∈ {1, ..., r}. Applying the projection operators |φ j φ j | subsequently from j = 1 to j = r, the probability that all measurements be successful is independently of the system dimension d.
The two maximally unbiased projectors used in MUM [17] [ Fig. 1 That is, each failed measurement ofP t is followed by a measurement ofP U t , which makes a subsequent successful outcome ofP t occur with probability 1/d. The success probability after r measurements becomes i.e. the probability of r consecutive failures 1−p(r) decays exponentially on a scale defined by d. The probability of r consecutive failures 1 − p(r), i.e. the probability of not yielding the desired target state after r steps, is shown for MEDO, MUM and FUMES in Fig. 2(a). The two latter clearly show exponential scaling, and FUMES consistently outperforms MUM: To understand why, consider a random time-evolution after which the projectorP t is applied [18]. Such a procedure yields, on average, a target-state occupation probability 1/d, just as for MUM. For FUMES, however, each state |η j characterised by a failed measurement is accompanied by an optimal waiting time that maximises the probability p j to populate the target state, and, on average, p j > 1/d. Since the target state can be populated after each measurement with finite probability, 1 − p FUMES (r) features an exponential decay, which allows us to define the effective dimensionality d eff via which quantifies the performance of FUMES with respect to MUM: For a system of dimension d, FUMES requires as many measurements as MUM does for the effective (and smaller) dimension d eff . The effective dimension d eff /d decreases with the system dimension d [see Fig. 2(b)], because optimising in a higher-dimensional space is more likely to yield a target-state occupation probability p j 1/d. Thanks to the r projectors onto states of the form (2), the performance of MEDO is independent of the system dimension. As illustrated in Fig. 2(a), the exponential character of FUMES guarantees that it always outperforms MEDO if the target failure probability is below a certain threshold. For large system dimensions d > 100, the target failure probability needs to be below 10 −2 before FUMES becomes favourable, for which, however, the actual feasibility of MEDO becomes improbable.
Fock-state generation -We aim at engineering quantum states of multi-mode systems of ultra-cold bosons described by the Bose-Hubbard Hamiltonian [20] where J is the inter-well tunnelling and U is the collisional interaction strength [ Fig. 3(a)]. The standard strategies of unitary control [3,4] or control-free engineering [8][9][10][11][12][13] require modifications in this setting: Changing the Hamiltonian parameters abruptly, as for bangbang control or the SU(2)-decomposition, compromises the lowest-band approximation due to vibrational excitations [7]. Continuous changes impose a speed-limit on achievable operations [5]. On the other hand, implementing non-destructive measurements of several precisely chosen non-commuting observables [8-11, 13, 17] in one experimental realisation remains impractical.
As an initial demonstration of control, we assume to be equipped with quantum non-demolition measurements of the local atom-numbers [16,[21][22][23][24], which leads to statecollapse onto a Fock-state [ Fig. 3(c)] or onto a superposition of Fock-states with certain properties [ Fig. 3(b,d)].
which counts the number of tunnelling events required to obtain | n starting from | m . The distance D motivates an intermediate level of granularity M = D between binary (M = 1) and granular (M = d − K) measurements, which we refer to as subspace measurements.
In anticipation of the experimentally relevant infrastructure discussed below, we illustrate state engineering by FUMES for the target |φ t = |0, 2, 0, 2, 0, 2 . The probability 1 − p(r) to remain unsuccessful after r measurements is shown in Fig. 3(e) for J/U = 1.5, for the three levels of granularity (binary: solid, subspace: dot-dashed, and granular: dashed). For all three granularities, numerical optimisation yields an average time between measurements of 2/J, which corresponds to 0.076 ms for 87 Rb in a conventional 512 nm optical lattice [25]. A finer granularity facilitates faster state-engineering, since it gives more detailed information about the current state of the system, which helps to choose the optimal time to applyP t . Cumulating a success rate of 99% with granular measurements thus takes approximately 30 ms.
To set the results into context, we also consider a random Hamiltonian with d = 462 and fully granular measurements (dotted). The structure of the Bose-Hubbard Hamiltonian implies that states with large distance D need many tunnelling events to be connected, while random Hamiltonians typically couple every pair of states. Therefore, FUMES performs better for random Hamiltonians than for the Bose-Hubbard Hamiltonian.
We also implement bang-bang control [4], alternating the HamiltoniansĤ 1 =Ĥ(J/U = 0) andĤ 2 = H(J/U = 100) (squares). We count the number r of cycles e −iĤ2t2 e −iĤ1t1 and interpret the fidelity of the target state preparation as the success probability p(r). While the non-commutativity of the two Hamiltonians guarantees controllability, bang-bang control [4] requires more control cycles than measurement cycles needed for FUMES: A measurement in FUMES strongly perturbs the system, which facilitates the fast population of the target state. In practice, a bang-bang control sequence needs to be chosen in advance, i.e. before the actual start of the procedure, while FUMES allows just-in-time optimisation: Every time a measurement fails to produce the target state, the timing for the next attempt is optimised.
Schrödinger cat generation -As a prominent application of FUMES, we demonstrate the near-deterministic generation of Schrödinger-cat states in the Bose-Hubbard model by the experimentially feasible measurements of the atomic number difference [ Fig. 3(d)] between the even and odd sites [22,23],Ẑ = L j=1 (−1) jn j /2, wheren j counts the atoms in site j. For vanishing interaction U → 0, macroscopic superpositions can be generated probabilistically [23], conditioned on a nonvanishing measurement result ofẐ. However, this approach typically yields states of very low macroscopicity [26], as quantified by where is the measurement operatorŜ( w) = L j=1 w jnj = w · n, where we restrict w to {±1} L , accounting the impossibility to directly measure superpositions of different particle numbers without auxiliary reservoirs [27]. The macroscopicity N eff (|ψ ) is the minimal number of particles for which the validity of quantum mechanics is required for a description of the observed macroscopic fluctuations [26]. A Schrödinger-cat state fulfils N eff = N , the groundstate of the Bose-Hubbard Hamiltonian for J/U = 1.5 and d = N = 6 gives N eff = 0.87 ≈ N/7, while single Fock-states do not carry any macroscopic entanglement and yield N eff = 0. The macroscopicity within subspaces of constant Z has a small spread and is given in Table. I. Thus, finding Z = N/2 is necessary to achieve Schrödinger-cat-like macroscopicity characterised by N eff ≈ N . However, measuring Z = N/2 in the ground state of the Bose-Hubbard model is highly unlikely: For a superfluid state of 100 bosons in 100 sites, the probability is less than 10 −30 , and even for N = 6 particles and L = 6 sites with J/U = 1.5, the probability remains less than 2%.
Using FUMES, we can exploit the tunnelling dynamics induced by the Hamiltonian to steer the initial state into the otherwise improbable subspace Z = N/2 [ Fig. 4]. The target space Z = N/2 = 3 has dimension 56, which facilitates state engineering in comparison to Fock-state generation [ Fig.3 (e)]. Larger values of J/U come with super-fluid-like eigenstates that mix more Fock-states than low values leading to Mott-insulator-like eigenstates. This makes larger values of J/U more favourable, consistent with Fig 4. Eight measurements then suffice to accumulate a success rate of more than 80% for J/U = 1.5, while the typical time between two measurements for J/U = 1.5 remains around 2/J, as for the Fock-state target.
Conclusion -The minimal requirements for unitary control are the availability of two non-commuting Hamil-toniansĤ 1 andĤ 2 [4]. We translated this setup into the realm of control-free engineering: A fixed Hamiltonian H and a projector onto the desired target subspaceP t then permit FUMES as long as the dynamics induced by the Hamiltonian effectively steers quantum states into the target subspace. In situations in which the Hamiltonian is uncontrollable or no tailored measurements can be applied, FUMES emerges as a natural way to perform state engineering. Its computational costs are manageable even for large system dimensions, since the waiting time until the next measurement can be optimised in a just-in-time fashion. Thanks to the optimisation of the waiting time, it outperforms methods for which the success probability of the measurement scales as 1/d [17,18].
Our greedy strategy that always aims at the desired target in the next measurement may be refined further: For granular measurements, the trajectory of a state steered by FUMES corresponds to a classical random walk on d − K nodes: d − K − 1 nodes correspond to nontarget states, one represents the target subspace. The probability to jump from one node to the other is then optimised to achieve the largest probability to populate the target node. For a state with very low -albeit optimised -success probability, it can be advisable to induce a detour and rather aim at another non-target state, provided it comes with a large probability to reach the target in a subsequent measurement. Besides the assessment of such more sophisticated strategies, it remains to be studied which ensembles of states and Hamiltonians can be engineered, and how the strategy performs on systems beyond the Bose-Hubbard-Hamiltonian.