Modelling assisted tunneling on the Bloch sphere using the Quantum Composer

The Bloch sphere representation is a geometric model for all possible quantum states of a two-level system that can be used to describe the time dynamics of a qubit. As explicit application, we consider the time dynamics of a particle in a double-well potential. In particular, we adopt a recent method for off-resonant excitations, the so-called SUPER principle (Swing-UP of the quantum emitter population) driven by periodic electromagnetic fields, to the context of quantum tunnelling. We show that the tunnelling probability can be enhanced significantly when an appropriate oscillation of the potential height is introduced. Driven by a collaborative approach we call educator-developer dialogue, an updated version of the software Quantum Composer is presented. For educational purposes, we map the two lowest energy states of the 1D-Schr\"odinger equation to the Bloch sphere representation, leading to a rather clear and intuitive physical picture for the pertinent time dynamics.


Introduction
With the advent of quantum technologies, the importance of two-level systems, also denoted as qubits, is increasing not only from a technological, but also from a didactical perspective.In many standard textbooks [1], the topic is addressed in a rather formal manner.Generally speaking, multiple representations are an important tool in science education [2].This holds true in particular in quantum physics, as the basic concepts are rather abstract and complex.Using the software tool Quantum Composer [3], the mathematical formalism comes to life by interactive graphical visualisations.
In this contribution, as an illustration, we discuss tunnelling in a double well as an explicit example.In particular, using an update of the Quantum Composer, we apply the so-called Bloch-sphere representation to the double-well system.The (single-particle) quantum double-well system is of great importance to many areas of quantum physics as it enables the implementation of a two-level quantum system, i.e. a qubit.Tunnelling in the double-well potential is extensively explored for the symmetric case [4] and the asymmetric (tilted) case [5].Periodically driven potentials as examined in this paper are usually described using the Floquet formalism [6][7][8][9].
The Bloch-sphere representation is of particular importance, as it can be applied to all two-level systems, independent of the particular physical interpretation.A geometric model for the time dynamics of qubits on the Bloch sphere is -for timeindependent Hamiltonians -well known: Just choose a rotation axis ⃗ n, which amounts to representing the 2 × 2 Hamiltonian as H = ϑ 2 ⃗ σ • ⃗ n.Then, the time dynamics of any initial state |ψ⟩ = α|0⟩ + β|1⟩ on the Bloch sphere translates to a rotation of this state around the axis ⃗ n on the Bloch sphere.There are two eigenstates, namely those parallel/anti-parallel to ⃗ n.This geometric picture is universal in the sense that it is independent of the explicit physical interpretation of the degrees of freedom described by the Hamiltonian.
A relevant concept to the time dynamics of two-level systems are so-called Rabi oscillations [10].These oscillations describe the behaviour of a two-level quantum systems with an external periodic driving force.This time-dependent driving force leads to an effective coupling between the two states.If the driving force is in resonance with the two-level system, the Rabi oscillations are full, i.e. the particle will oscillate between the two levels completely.Otherwise, the particle will only partially reach the other level.We emphasise that throughout this paper, the term "Rabi oscillations" is used metaphorically as we don't explore systems with time-dependent external driving.Rather, in the double well potential, a time-independent coupling between the states localised in the left and the right well is present due to the tunnelling, which can only be changed by changing the shape of the potential.Due to the formal equivalence of the time-indepedent and the time-dependent situation, we adopt the notion of "Rabi oscillations" also to the case of tunnelling.
Multiple representations introduced for Rabi oscillations can be useful for undergraduate courses on quantum physics.We provide explicit animations and quantum composer files for full Rabi oscillation (Level 1) and partial Rabi oscillation (Level 2).Next, for graduate level and beyond, we consider a time dependent doublewell potential (Level 3).
This paper is organised as follows: In Sect.2, the Quantum Composer by Quatomic is introduced.In Sect.3, we review the well-known general theory for Rabi oscillations of two-level systems.We discuss the usual co-moving frame approach and contrast it with the situation of time-independent Hamiltonians, where this transformation is not needed.In Sect.4, the tilted double-well potential is introduced as an explicit model system for a time independent Hamiltonian.We compare analytic results for the degenerate case with numerics from the Quantum Composer (Level 1).Full Rabi oscillations, that is, full oscillation between the left and the right minimum of the double-well potential emerge in this case of a symmetrical double well.In the tilted potential, we have partial oscillation between the left and the right minimum of the double-well.In Sect.5, we introduce the so-called SUPER principle on the Bloch-sphere, which generalises the simple model for the time dynamics of a qubit to a situation where the rotation axis ⃗ n(t) on the Bloch sphere itself is a certain function of time, allowing for the transition from a given initial state to a given final state with almost 100% probability by adjusting ⃗ n(t) accordingly.In Sect.5, we apply the SUPER principle, which has been introduced in the context of off-resonant excitations [11], [12] to the situation of quantum tunnelling in the tilted double well as model system, which is the main result of the present contribution (Level 3).Note that in contrast to [11], we do not use the co-moving frame ansatz.In Sect.6, we discuss further extentions and applications of this work, and its possible didactical merits.

The Quantum Composer
In this work, we will use the digital tool Quantum Composer for numerical simulation of the mechanics of a tilted double well potential.Composer was developed in the group of Prof. Sherson at Aarhus University [13].Unlike many other fixedinterface applet approaches [14][15][16], which are limited to showcasing particular sets of phenomena, Composer provides an intuitive and user-friendly re-programmable, node-based graphical interface, as shown in Fig. 1.These nodes can be connected to construct quantum systems which are simulated by the underlying QEngine, a C++ library for the simulation and control of one-dimensional quantum systems [17].These systems under investigation can be generated completely from scratch, or take as a basis pre-defined scenarios loaded into the program using .flowfiles.This flexibility makes the tool suitable for both education (use by teaching staff and students), and for research purposes [3,18,19], with a low barrier to entry due to the drag-and-drop nature of the programming.The Quantum Composer can be downloaded for free for most desktop platforms [13].Due to its hands-on re-programmability, Quantum Composer is uniquely suited to inquiry-based learning (IBL) [20], in which students interact with didactic material through active discovery, questioning, evaluation, and sensemaking.When used for teaching quantum mechanics, IBL has shown potential for improving understanding and reducing overconfidence [21,22].Despite this, there are relatively few examples of publicly available resources for inquiry-based quantum mechanics teaching.
There are many simulation and visualisation tools [23] already available that can be invaluable in promoting understanding, interest, and intuition for quantum concepts [24].However, those which are most widely used for teaching, such as QuVis [14] from and QuILTs [25] offer a fixed set of scenarios covering different phenomena and allow students to vary parameters.While these are effective to enable educators to set up guided and partially inquiry-based teaching situations, they do not allow for truly open and exploratory inquiry, in which the system is entirely set up from scratch by the participants.Thus, such tools may miss out on some of the benefits associated with more open-ended IBL in understanding and promoting scientific thinking [26].Composer, however, being entirely reprogrammable, allows for teaching material to be developed which can be anywhere on the spectrum of fully guided to completely open-ended and exploratory.An example IBL approach was previously explored using Quantum Composer to devise student exercises based on the Predict, Explore and Revise approach [27].
For researchers, Composer may be useful as a relatively quick and simple means to set up, visualise, and solve analytically intractable quantum systems.The drag-anddrop interface and parameters which can be adjusted on the fly makes it a valuable alternative to more conventional approaches taken by researchers, such as the use of MATLAB or Python code which can be less intuitive to use than the node-based UI of Composer [3,18].The tool has previously been used to investigate the dynamics of atoms in atomic lattices [28] and to solve quantum optimal control problems [17,19].
In the context of this work, however, we wish to highlight another aspect of the Composer re-programmability: the ability of teachers and instructors to seamlessly create custom demonstrators to suit the particular student needs.Here, instructors in advance prepare the relevant .flowfiles and either use them as lecture demonstrators or distribute them to the students with suitably highlighted experimentation parameters and visualisations preconfigured.In addition, the Composer design team is continually available to support the addition of new interaction and visualisation features.The current work illustrates exactly such interactive collaboration between two instructors (J.B. and S.H.) with the Composer design team resulting in the addition of new features to Composer.
In addition, the Composer design team is continually available to support the addition of new interaction and visualisation features.The current work illustrates a form of collaboration between educators (J.B, S.H), and the Composer developers (V.M, S.G, and J.S) which we call educator-developer dialogue.In this interactive approach to co-design of educational material, the digital tool is able to act as an instrument to respond to community needs.
There are two possible avenues for this dialogue.The first consists of educators designing material with close input from the developers of the tool.This can be invaluable in identifying and preparing validation criteria, as in reference [3].The second mechanism, which is demonstrated in this article, is where developers implement new features responding to the needs of the educators, in this case the Bloch-sphere representation.Here we demonstrate the value of Composer by extending the body of existing educational material, using the example of Rabi oscillations in the double well potential.Quantum Composer files are provided for three different levels: full Rabi oscillations (Sect.4.1), partial Rabi oscillations (Sect.4.2) and Swing Up with the SUPER principle (Sect.5).The corresponding files can be downloaded from [29].By combining the Bloch sphere representation with the representation of wave functions in the double-well, Composer helps us to develop an intuitive physical picture of the time dynamics in the time-dependent double well using multiple representations.This is the first application of an update of the Quantum Composer incorporating the Bloch sphere representation, driven by an educator-developer dialogue.

A brief review of the general theory of Rabi oscillations
In this section, we briefly review the general theory of Rabi oscillations, which holds for many different applications.Crucially, we distinguish between two-level systems driven by periodic fields, and those that are not driven by periodic fields.

Two-level systems driven by periodic fields
We start our discussion with the usual textbook example of a two-level systems driven by a periodic field.The unperturbed states are defined as |g⟩ (ground state) and |e⟩ (excited state) at energies E g , E e .Let ∆e iΩt be the interaction between the driving field and the two-level system that introduces a coupling between the two eigenstates, where Ω is the driving frequency of the periodic field and A the (complex) amplitude.This system is described by the following Hamiltonian In matrix notation, we obtain (with mean energy Ē = (E e + E g )/2 and energy gap η ≡ E e − E g between excited and ground state) where I is the unit matrix.The Schrödinger equation describing the dynamics of the system is given by iℏ∂ t |ψ⟩ = H 1 (t)|ψ⟩.An elegant way to compute the time dynamics of the state |ψ(t)⟩ = z g (t)|g⟩+z e (t)|e⟩ is the transformation to the so-called co-moving frame: Using the unitary operation , where σ 3 is the third Pauli matrix, we define by |ϕ⟩ ≡ C(t)|ψ⟩ the state in the co-moving frame.Within the co-moving frame, the time-dependence of the state |ϕ(t)⟩ is then given by with In the co-moving frame, the Hamiltonian H CM becomes time-independent, while the diagonal terms ±η are shifted to ±(η − ℏΩ).That is, the time dependence ∆e iΩt of the interaction is transformed to a change in the energy shift in the co-moving frame.
Resonance emerges for η ≡ E e − E g = ℏΩ which can then be interpreted as absorption or emission of a photon in the transition between the states |e⟩ ↔ |g⟩.

Two-level systems without external driving by periodic fields
In contrast to the textbook example for Rabi oscillations induced by external periodic fields, in this contribution, we consider two-level systems without external driving.In other words, the Hamiltonian is time-independent from the very beginning, and no transformation to the co-moving frame is necessary in order to achieve this.Due to the formal equivalence between (2) and ( 4), all mathematical results can be adopted, while the interpretation changes significantly.Note that in our approach to tunnelling, we therefore do not work in the co-moving frame.Consider a time-independent two-level system with states |L⟩ and |R⟩ at energies E L , E R with mean energy Ē and energy difference δ = E L − E R .Let ∆ be the (real) interaction between the two states (which will later on be interpreted as transition probability between the left and right potential minimum, see Eq. ( 19)).This system is described by the following time-independent Hamiltonian First, we want to exploit the general results following from this Hamiltonian.Later on, we use this Hamiltonian to approximate the two lowest energy states of the doublewell potential.We introduced the so-called Rabi frequency ℏω R = √ δ 2 + ∆ 2 and the rotation angle θ, which describes the rotation axis of the corresponding time development in the unitary operator ⃗ σ are the three Pauli-matrices.Indeed, on the Bloch sphere, this operator generates a rotation of any given initial state around the rotation axis On the Bloch sphere, the eigenstates point in the ±⃗ n directions.In Hilbert space, these eigenstates are given by The eigenvalues are given by For a given initial state |ψ(0)⟩, the fidelity under time development 2 .For a general superposition state with |α| 2 + |b| 2 = 1, the time dynamics is given by (up to a global phase) and the fidelity, i.e. the tunnelling probability, reads For |ψ(0)⟩ = |L⟩, we obtain For 2θ = π/2, the Rabi oscillation is maximal.The fidelity is then given by The state |ψ 0 ⟩ = α| + ⃗ n⟩ + β| − ⃗ n⟩ can be mapped on the Bloch sphere via Hopf mapping (the qubit vector in C 2 is mapped to the sphere S 2 , see e.g.[30]).Explicitly, On the Bloch sphere, there is a simple geometric interpretation for the behaviour of the fidelity.In Fig. 2, we show the rotation axis ⃗ n of the unitary time development, and the orbits of states at angle θ = 0, π (the eigenstates), and θ k = kπ/6, k = 1, . . . 5.
In analogy to the time dependent system described, for 2θ = π/2, |L⟩ and |R⟩ are part of the same orbit.Therefore, the transition between the two states is certain, leading to the fidelity function (16).We therefore call this full Rabi oscillation, in analogy to the time dependent system described.

Tunnelling in the double well
Quantum tunnelling plays a major role in various areas of modern physics.Early applications have been e.g. the calculation of lifetimes for alpha-decay of nuclei, or for nuclear fusion.In recent times, from a technological perspective, tunnelling in solid state physics, and in particular in semiconductor physics and in spintronics [31] are of interest both for the first and in the second quantum revolution.Furthermore, model systems such as Bose-Einstein condensates in a double-well potential are part of actual research [32,33].In the past decades, periodically driven tunnelling systems also These two values of λ will be used to realise assisted tunnelling with the SUPER principle which is shown in Sect. 5.The theoretical prediction for the fidelity 1 − F (t) = 4|α| 2 |β| 2 sin 2 (ω R t/2) coincides with the numerical result and cannot be distinguished within the resolution of these plots.
have been widely studied [6].Since tunnelling is an important topic in physics since the advent of quantum physics, it is not astonishing that also in physics education, quantum tunnelling is part of the standard curriculum not only at university, but also in high school [34].Thus, it is desirable to explore model systems for quantum tunnelling at all levels of education, and also for research.

Full Rabi oscillations in the double well (Level 1)
As a model for the tunnelling system, we consider the (tilted) double-well potential If we consider the system to consist of two subsystems L and R, we find respective ground states |L⟩ and |R⟩.First, we want to find a relation between the parameters (δ, ∆) of the general theory of Sect.3, and the parameters (λ, a) of V (x).For this purpose, we investigate the degenerate case µ = δ = 0 where the ground state energies E L and E R are equal.We will calculate ∆ = ∆(λ, a) within this model potential.In the degenerate case, we obtain full Rabi oscillation between the state on the left |L⟩ = 1/ √ 2(|E 2 ⟩ − |E 1 ⟩) and on the right side of the potential |R⟩ = 1/ √ 2(|E 2 ⟩ + |E 1 ⟩) (see Eq. 9 and 10 for θ = π/4).The eigenenergies are given as In this case, the coupling ∆ can be computed as tunnelling amplitude either using WKB or the instanton ansatz.In both approaches, we obtain [35,36] Here, x inst is the so-called instanton solution, a classical path in imaginary time τ connecting the left and the right potential minimum.The instanton action is given by S with ω H.O. = (8a 2 λ)/m.In such a way, we find explicit and analytic expressions for all parameters in the general Hamiltonian (1) in terms of the model potential (18) for the degenerate case µ = 0.

Partial Rabi oscillations in the tilted double well
As an approximation, we can consider the two wells to be completely separate, i.e. the system to consist of two independent subsystems L and R. Let E L be the energy of the ground state |L⟩ of the left minimum, and E R be the energy of the ground state |R⟩ of the right minimum.Concerning these two ground state energies, we may approximate the respective potentials as harmonic oscillators.We find in the vicinity of x ≃ a with y ≡ (x − a) V (y) ≃ µ(y + a) + (2a) 2 λy 2 + . . .
In first order approximation, the µ-dependence is irrelevant for the ground state apart from a vertical shift by −µa or µa, respectively, as the term µ(x ± a) is odd around x = ∓a, and the harmonic oscillator ground state is even around its origin.Thus, we can compare the ground state energies with the usual ground state energy of the harmonic oscillator.We find E L = 1 2 ℏω H.O. − µa and E R = 1 2 ℏω H.O. + µa with ℏω H.O. = (8a 2 λ/m).Thus, we find in leading order δ = (E L − E R ) ≃ −2µa.While it is possible to refine the (rather bad) approximations for E L/R using asymptotic series [35], we refrain from doing so here as our main concern is the introduction of assisted tunnelling using the Quantum Composer.
We will discuss to what extent we can verify the theory of Rabi oscillations using the Quantum Composer.This is the second level of studying Rabi oscillations in the double-well potential with the Quantum Composer.Natural units (ℏ = 1 and m = 1) are used throughout the rest of this paper.Since the eigenvalues for the energy are and the mean energy by E 1 + E 2 = 2 Ē, obtaining E 1 and E 2 numerically using the Quantum Composer.In Fig. 3, we compare the theoretical expression for the fidelity 1 − F (t) = 4|α| 2 |β| 2 sin 2 (ω R t/2) (16) with the results from the quantum composer for two different potentials (λ 1 = 35, λ 2 = 17).We find excellent agreement.
In the tilted potential, the instanton calculations that were shown in Sect.4.1 for full Rabi oscillations are still possible, but much more involved [35].We will not embark into these troubled waters of theory not to loose track to our main purpose, the introduction of assisted tunnelling.However, it would be interesting for future work to link also these analytic results with numerical results obtained with the Quantum Composer.
To conclude, Rabi oscillations can be studied in an elegant manner with the Quantum Composer for the tilted double well potential and can give valuable insights into some important aspects of the underlying theory.By swapping between them, we can go from initial state |A⟩ to state |B⟩ following the red path, or we can go from |B⟩ to |A⟩ following the pink path.When falling out of the swinging rhythm, there is the danger to 'loose track'.Once the state is at |B⟩, if the Swing Up process is not stopped it will follow the pink path and start decaying back to |A⟩.Note that in a full Swing Up, a state viewed in 2D would evolve along the red path and then "back" along the pink path to end in the desired state, orthogonal to the initial state.

Application of the SUPER-principle to assisted tunnelling (Level 3)
The SUPER principle was originally introduced in [11], while first experimental realisations have been reported in [12] for two level-systems driven by periodic fields.Note that the rotating frame transformation must be applied in this case.
In the rotating frame, the time-independent Hamiltonian (5) arises.As the tunnelling system is of the type of equation ( 5) from the very beginning, no transformation to the rotating frame is necessary.This necessity is the only difference between the original SUPER principle and the ideas applied here.This allows for the application to assisted tunnelling, as we will show in what follows.
The key idea of the SUPER principle is the following: For given rotation axis ⃗ n, the orbit of all initial states under the action of (6) form a "ring system" as shown in Fig. 4. If the rotation axis ⃗ n(t) oscillates between two different, but near positions, the time development of an initial state switches between these two "ring systems".
In [29], we show an animation of the time evolution of the state |L⟩ to the state |R⟩ based on the super principle, applied to the situation shown in Fig. 5. Indeed, as we show in this section, if the oscillation period is chosen nearby the average Rabi frequency of the two different unitary operators (6) with rotation axis ⃗ n 1 and ⃗ n 2 , a "Swing Up" leads to an enhancement of the Rabi amplitude even in the case of a tilted well.From our point of view, this result can be seen as analogous to [11] in the case of Rabi oscillations with an external, time-dependent driving field, where it was shown that the modulation frequency has to be chosen close to the mean Rabi frequency of the two detunings.In order to apply the SUPER principle to the particular situation of tunnelling in a double-well, we consider a situation where the potential height V (0) = λa 4 for constant distance a is periodically shifted from a 'high' value to a 'low' value by changing λ.The potentials used, and the initial and desired states are shown in Fig. 6.We consider the situation where we start in state For large tilting parameter µ, the maximum fidelities rather low (see Fig. 3 with P 1 = 0.304 and P 2 = 0.362, respectively).
shown in blue, with eigenenergies in red.The corresponding Bloch-sphere representation is shown schematically in Fig. 5 Next, we modify the system and allow for an oscillation between V 1 or V 2 , where the potential oscillation frequency is about the mean Rabi frequency.We introduce a smoothstep function to model the transition between the potentials.After some fine tuning that is described in 5.1, the maximum Fidelity was increased drastically to P = 0.994.The result is shown in Fig. 7, confirming the Swing up mechanism, leading to assisted tunnelling.
Besides the SUPER principle, other approaches to enhance the tunnelling rate are known.Here, the notion of parametric resonance is used in a variety of different situations.Our approach differs from parametric resonance as introduced in [37], as we consider a tilted potential, and change the coupling between the two wells in an appropriate way by changing the barrier height.We do not use the rotating wave approximation.In [38], the amplitude of an external driving field is changed to enhance the tunnelling rate.In both cases, parametric resonance works with a weakly nonlinear, parametrically modulated oscillator described in the rotating wave approximation with tilting µ = 0.
Other approaches to enhance tunnelling -also denoted parametric resonancehave been proposed in [39] and [40].Here, similarly to our approach, the potential height is varied periodically to enhance tunnelling.In this case, a Bose-Einstein condensate is discussed, where the time dynamics is governed by the Gross-Pitaevskii equation.
In the future, it will be interesting to study the relation of the SUPER principle to these approaches for parametric resonance in more detail.

Practical hints towards realising assisted tunnelling with the Quantum Composer
Assisted tunnelling by swapping between two double-well potentials with different barrier heights λ 1 and λ 2 is not an easy task.The parameters have to be chosen carefully to achieve convincing results.We start by choosing two sets of parameters for two different potential functions where a is fixed and would typically be given by the physical system, λ 1 ≈ 2λ 2 and µ 1 = µ 2 .The closer λ 1 and λ 2 are, the closer the two corresponding rotation axes on the Bloch-sphere.Then, the orbits (see Fig. 4) are closer together, resulting in less distance travelled per oscillation period.This translates to more oscillations needed for a full Swing Up.The fidelity graph becomes less steep, and in turn the Swing Up process takes longer.Thus, we choose substantially different λ.
The Rabi frequencies ω R and ω R can be found via E (1) 2 and E 2 , respectively, that are calculated numerically in Quantum Composer.Using Eq. ( 7) and double-angle formulas, we find δ/ω R = cos 2 (θ) − sin 2 ((θ)) = 2 cos 2 (θ) − 1 = 1 − 2 sin 2 (θ).With (9), we calculate the initial state |L 1 ⟩ localised in the left well as = ω (1) R − 2µa 2ω using the approximation δ ≈ 2µa.The desired state can be either chosen as or depending on which of the two potential functions should be final.In our example Swing Up, we chose |R 2 ⟩.Next, the potential swap time T = 2π ω ( where ω is the potential oscillation frequency) should be chosen.T /2 is the time of one potential being active before swapping to the other.For an optimal fidelity, we observed that T should be about 5% larger than the average Rabi period T R of the two potentials, T R = 4π/(ω Next, the time offset t 0 at the start of the swing-up can be optimised (see Fig. 8).T /4 + t 0 is the time of the first swap of the potential.
Using these parameters, we achieve a high fidelity for assisted tunnelling in the following manner: First, one should optimise the time offset t 0 for the starting potential swap time T .Next, we optimise T for this offset and repeat.If the resulting fidelity is not satisfying, different values for λ 1 , λ 2 or µ should be chosen.In this manner, a maximum fidelity of over F=0.994 was achieved as is shown in more detail in Fig. 9.
The process is challenging in the sense that it is very sensitive to T , t 0 and the potential parameters λ and µ.This could be because the Rabi frequency is not constant during the Swing Up process.Rather, it increases the further |ψ(t)⟩ travels towards |R⟩ as shown in Fig. 10.We have no analytical results that predict this behaviour and can only refer to the numerical results.This can lead to falling out of the swinging rhythm, but can be solved by careful choice of parameters, as shown in this paper.

Limits of the Bloch-sphere representation for the double-well system
The 2 × 2 Hamiltonian (5) includes only the two lowest-lying states of an infinite number of eigenstates in the double-well.In this section, we check the validity of the two-level system approximation ( 22), (24) for the two lowest-lying states localised in the left and the right well, respectively.
In general, these states are a superposition not only of the first two eigenstates |E 2 ⟩,|E 1 ⟩, but include more states |χ⟩ R,L with higher energy.Thus, the exact wave functions in position space reads These states are fully localised either in the left or the right valley, leading to In Fig. 11, we show the expression 1 − |χ R (x)| 2 depending on the values for λ, µ.In the region with ω R > 2µa, the deviation from the two-level system is below 5%.Note that for ω R = 2µa, the amplitude for |E 1 ⟩ in (22) (and |E 2 ⟩ in (24)) vanishes.After the transition to ω R < 2µa, the amplitude becomes imaginary, and the error increases in an asymmetric manner.While the error stays below 15% for the left valley, it increases to up to 35% in the right valley.While the visualisation of the SUPER principle on the Bloch sphere is rather intuitive and helps understanding the mechanism, the SUPER principle itself is not restricted to the 2 × 2 level system and indeed includes higher order states during the transition shown in Fig. 7.

Outlook and discussion
Our conclusions can be drawn from a didactical perspective, and from the perspective of recent research in applied physics.While periodic oscillations of potentials have already been studied [6], to our knowledge, a systematic enhancement of tunnelling using the SUPER-mechanism has not yet been demonstrated.
In the context of quantum tunnelling, it seems possible to find experimental realisations for a "high" value and a "low" value of the potential barrier.Explicitly, we may think of a Bose-Einstein condensate trapped in a double-well potential.An experimental realisation, where the atoms are magnetically coupled to a single-mode of the microwave field inside a superconducting resonator, has been reported in [41].24).In the region ω R > 2µa, the deviation from the two-level system is rather small < 10%.
By modulating the magnetic field strength, the trap potential could be manipulated to embark for assisted tunnelling.
As another model system, we may think of a double-well potential realised for an individual fluxon trapped in a Josephson tunnel junction.Again, parameters of the potential can be manipulated by the strength of a magnetic field [42].Since the double-well potential is the smallest subpart of the general class of Hubbard-model systems in condensed matter physics, many other applications of assisted tunnelling might be of interest.
Our main concern however, is physics education: Since the Quantum Composer is available for free, it is possible to discuss Rabi oscillations in an intuitive and direct manner using this powerful modelling tool on almost any computer system.The result of the educator-developer dialogue is the updated version of the Composer.With that, the explicit visualisation of the Bloch sphere dynamics is possible.Using the Bloch-sphere representation, the SUPER principle has a simple and general geometric interpretation as shown in Figs. 2 and 4 and even a knot theoretic interpretation [43].Translated back to the double-well, we may apply this mechanism to realise assisted tunnelling by introducing an appropriate oscillation of the potential height.
Explicitly, at high school, the dependence of the Rabi frequency on the potential parameters (λ, a, µ) can qualitatively be discussed using just the graphical interface provided by the Quantum Composer.On the undergraduate level, the quantitative relations explored in Sect. 4 can be accessed, in particular, the dependence of S inst.∝ √ λa 3 on the level spacing of the two lowest eigenvalues, (E 1 − E 2 ) ∝ exp (−S inst.).On the graduate level, assisted tunnelling and the SUPER principle can be explored numerically as shown in Sect. 5. Here, numeric results can be used to provide a qualitative understanding of the SUPER-mechanism as shown in Fig. 4. The multiple representation that the Quantum Composer offers, like the state evolution plot and the Bloch sphere, can enhance understanding of this concept.All Quantum Composer files we used in our work can be downloaded from [29].
Concerning research in physics, still many questions remain open.It would be desirable to understand the relation between the numerical results from the Quantum Composer and analytic results, such as 19, even better.Moreover, analytic results for assisted tunnelling are not yet available, and it would be desirable to analytically explore this approach further in accordance with [44].
An empirical study on solving tasks with multiple representations in quantum composer has been conducted recently in an eye-tracking study [45].Similarly, a possible direction of further didactical research is the exploration of the benefits and challenges of teaching approaches using visualisation tools such as the Quantum Composer, with various target groups and in different contexts like the one described in this paper.

Figure 1 .
Figure 1.An illustration of the drag-and-drop, node-based, graphical user interface of Quantum Composer.

Figure 2 .
Figure 2. Coloured Orbits on the Bloch sphere for different rotation axis angles 2θ.The rotation axis ⃗ n, defined by the first two energy eigenstates of the system E 1 and E 2 , is shown in red.Left: Partial Rabi oscillations, time development of the initial state |L⟩ for 2θ = π/6, π/4, π/3 and 5π/12.Right: Full Rabi oscillation, 2θ = π/2.In this case, the transition probability to state |R⟩ is 100% after a π−flip.

Figure 4 .
Figure 4.The Swing Up principle illustrated in 2D.Two different potentials V 1 and V 2 have two different rotation axis ⃗ n 1 and ⃗ n 2 pointing towards the reader.By swapping between them, we can go from initial state |A⟩ to state |B⟩ following the red path, or we can go from |B⟩ to |A⟩ following the pink path.When falling out of the swinging rhythm, there is the danger to 'loose track'.Once the state is at |B⟩, if the Swing Up process is not stopped it will follow the pink path and start decaying back to |A⟩.Note that in a full Swing Up, a state viewed in 2D would evolve along the red path and then "back" along the pink path to end in the desired state, orthogonal to the initial state.

Figure 5 .
Figure5.Relation between potential barrier height and the corresponding time dynamics on the Bloch sphere.By switching between these two potentials, on the Bloch-sphere, the SWING up mechanism can be realised as shown in Fig.4

Figure 9 .
Figure 9. Various maximum fidelities and the times to reach these fidelities for different potential swap times T /2 and different time offsets t 0 holding the potential swap time constant (T /4 + t 0 is the time of the first potential swap, see Fig. 8).Swapping between two potentials with µ = 2.8571, a = 0.7, V 1 with λ 1 = 35 and V 2 with λ 2 = 17.Left: The maximum fidelities reached for various potential swap times T /2 at time offsets t 0 = 0.1T /2.Shown are the Rabi oscillation periods of the two potentials T (1) R /2 = π/ω (1) R and T (2) R /2 = π/ω (2) R and the mean Rabi oscillation period T R /2.The maximum fidelity F = 0.987 is achieved at Tmax/2 = 1.052(TR /2) = 1.01.This can be further improved to F = 0.994 by holding the potential constant at λ 2 = 17 after t = 15 as shown in Fig. 7. Right: Maximum fidelities for various different time offsets 2t 0 /Tmax at potential swap time Tmax/2.

Figure 10 .
Figure 10.Fidelity F(t) as in Fig. 7, top middle and modulation function f (t) as in Fig. 8 in red.f (t) = 0 corresponds to λ = λ 1 = 17 and f (t) = 1 corresponds to λ = λ 2 = 35.The frequency of the Swing Up increases during the process: the first 4 peaks are 6.35 time units apart while the next 4 peaks are only 5.44 time units apart. 0