Non-adiabatic control of quantum energy transfer in ordered and disordered arrays

An elementary excitation in an aggregate of coupled particles generates a collective excited state. We show that the dynamics of these excitations can be controlled by applying a transient external potential which modifies the phase of the quantum states of the individual particles. The method is based on an interplay of adiabatic and sudden time scales in the quantum evolution of the many-body states. We show that specific phase transformations can be used to accelerate or decelerate quantum energy transfer and spatially focus delocalized excitations onto different parts of arrays of quantum particles. We consider possible experimental implementations of the proposed technique and study the effect of disorder due to the presence of impurities on its fidelity. We further show that the proposed technique can allow control of energy transfer in completely disordered systems.


Introduction
The experiments with ultracold atoms and molecules trapped in optical lattices have opened a new frontier of condensed-matter physics research. The unique properties of these systems -in particular, large (> 400 nm) separation of lattice sites, the possibility of tuning the tunnelling amplitude of particles between lattice sites by varying the trapping field and the possibility of controlling interparticle interactions with external electric or magnetic fields -offer many exciting applications ranging from quantum simulation of complex lattice models [1,2,3,4,5,6,7,8,9,10] to the study of novel quasi-particles [13] that cannot be realized in solid-state crystals. In the limit of strong trapping field, ultracold atoms or molecules on an optical lattice form a Mott insulator phase, in which each lattice site is populated by a fixed number of particles. With one particle per lattice site, this phase represents a periodic, crystal-like structure. Such a system can be thought of as a prototype of an ordered system, in which a single lattice site (or a small number of lattice sites) can be individually addressed by an external field of a focused laser beam. This can be exploited for engineering the properties of quantum many-body systems by changing the energy of particles in individual lattice sites [14].
In the present work, we consider the generic problem of energy transfer -i.e. the time evolution of an elementary quantum excitation -in such a system. In particular, we explore the possibility of controlling energy transfer through an array of coupled quantum monomers by applying monomer-specific external perturbations. This is necessary for several applications. First, collective excitations in molecular arrays in optical lattices have been proposed as high-fidelity candidates for quantum memory [15]. The ability to manipulate collective excitations is necessary for building scalable quantum computing networks [16]. Second, ultracold atoms and molecules in optical lattices can be perturbed by a disorder potential with tunable strength [17]. Engineering localized and delocalized excitations in such systems can be used to investigate the role of disorder-induced perturbations on quantum energy transfer, a question of central importance for building efficient light-harvesting devices [18]. Third, the possibility of controlling energy transfer in an optical lattice with ultracold atoms or molecules can be used to realize inelastic scattering processes with both spatial and temporal control. Finally, control over energy transfer in quantum systems can be used for studying condensed-matter excitations and energy transport without statistical averaging.
An excitation of a coupled many-body system generates a wave packet representing a coherent superposition of single-particle excitations. The method proposed here is based on shaping such many-body wave packets by a series of sudden perturbations, in analogy with the techniques developed for strong-field alignment and orientation of molecules in the gas phase [19]. Alignment is used in molecular imaging experiments and molecular optics [19,20,21,22], and is predicted to provide control over mechanical properties of molecular scattering [23,24]. Here, we consider the use of similar techniques for controlling quantum energy transfer in a many-body system. When applied to a completely ordered system, the proposed method is reminiscent of the techniques used to move atoms in optical lattices, where a uniform force is applied for a short period of time [25]. The conceptual difference comes from the fact that in the present case the momentum is acquired by a quasi-particle -a collective excitation distributed over many monomers. During the subsequent evolution, the particles do not move -rather, the excitation is transferred from one monomer to another. In order to control such excitations, we exploit an interplay of the adiabatic and sudden time scales, which correspond to single-monomer and multi-monomer evolution. We also exploit the wavelike nature of the excitation wave function to draw on the analogy with wave optics. This analogy, too, is not complete due to the discrete nature of the lattice.
In order to emphasize the generality of the proposed method, we formulate the problem and present the results in Sections 2 and 3 in terms of the general Hamiltonian parameters. Section 4 then describes how the external perturbations corresponding to the results presented can be realized in experiments with ultracold atoms and molecules. Section 5 discusses controlled energy transfer in systems with, specifically, dipole -dipole interactions. Section 6 considers the effects of lattice vacancies on the possibility of controlling energy transfer and Section 7 extends the proposed technique to control of excitation dynamics in strongly disordered arrays with a large concentration of impurities. Section 8 presents the conclusions.

Sudden phase transformation
Consider, first, an ensemble of N coupled identical monomers possessing two internal states arranged in a one-dimensional array with translational symmetry.
The Hamiltonian for such a system is given by where |g n and |e n denote the ground and excited states in site n, ∆E e−g is the monomer excitation energy and α(n − m) represents the coupling between two monomers at sites n and m. The singly excited state of the system is In general, the expansion coefficients C n are complicated functions of n determined by the properties of the system, in particular, the translational invariance or lack thereof as well as the strength of disorder potential. If an ideal, periodic system with lattice constant a is excited by a single-photon transition, the expansion coefficients are C n = e iakn / √ N and |ψ exc ⇒ |ψ exc (k) represents a quasi-particle called Frenkel exciton, characterized by the wave vector k [26]. The magnitude of the wave vector k is determined by the conservation of the total (exciton plus photon) momentum. The energy of the exciton is given by E(k) = ∆E e−g + α(k) with α(k) = n α(n)e −iakn . In the nearest neighbor approximation, where α = α(1). With atoms or molecules on an optical lattice, it is also possible to generate a localized excitation placed on a single site (or a small number of sites) by applying a gradient of an external electric or magnetic field and inducing transitions in selected atoms by a pulse of resonant electromagnetic field [27]. The presence of a disorder potential, whether coming from jitter in external fields or from incomplete population of lattice sites, also results in spatial localization. Similar to how Eq. (2) defines the collective excited states in the basis of lattice sites, any localized excitation |ψ can be generally written as a coherent superposition of the exciton states |ψ exc (k) with different k: Control over energy transfer in an ordered array can be achieved by (i) shifting the exciton wave packets in the momentum representation (which modifies the group velocity and the shape evolution of the wave packets) and (ii) focusing the wave packets in the coordinate representation to produce localized excitations in an arbitrary part of the lattice. To achieve this, we propose to apply a series of site-dependent perturbations that modify the phases of the quantum states of spatially separated monomers. These phase transformations change the dynamics of the time evolution of the collective excitations. Here we consider the transformations leading to acceleration or deceleration of collective excitations, while the focusing phase transformations are described in Section 3.
For modifying the group velocity of a collective excitation, the essential idea is to add a factor e iδan to each term in the expansion (2), so that each |ψ exc (k) component in a wave packet is transformed into |ψ exc (k + δ) . This transformation shifts the wave packets by δ in k-space while preserving their shape. As a result, one can engineer wave packets probing any part of the dispersion E(k) leading to different group velocity and shape evolution. The feasibility of such transformation in an ensemble of atoms or molecules on an optical lattice is discussed below and in Section 4.
Adding a site-dependent phase to the excitonic wavefunction exploits an interplay of the adiabatic and sudden time scales. Consider the n-th monomer subjected to an external field E n (t) which varies from 0 to some value and then back to 0 in time T . If the variation is adiabatic with respect to the evolution of the free monomer states, T /∆E e−g , each eigenstate |f of the monomer acquires a state-dependent phase shift [28] |f n (T ) = e −iφ f n |f n (0) , is the instantaneous eigenenergy and f can be e or g. Now consider the action of such phase change on the collective excitation state (2).
If T /α, the change is sudden with respect to the excitation transfer between monomers and the state (2) acquires a site-dependent phase Φ n = φ e n − φ g n . If Φ n = Φ 0 + naδ, then the momentum δ is imparted onto the excitonic wavefunction. By analogy with "pulsed alignment of molecules" [19], we call this transformation a "phase kick" or "momentum kick". Its action is also similar to that of a thin prism on a wavefront of a monochromatic laser beam.
In order to illustrate the shifting of exciton wave packets in the momentum space, we solve numerically the time-dependent Schrödinger equation with the unperturbed Hamiltonian (1), subjected to a transient site-dependent external perturbation that temporarily modulates ∆E e−g . We choose the parameters ∆E e−g , α and the lattice constant a that correspond to an array of polar molecules trapped in an optical lattice, as described in Section 4. The time-dependent perturbation, chosen to vary almost linearly along the lattice, has the form of a short pulse with the duration T = 3 µs. The phase acquired by the particles during this time is given by Φ n Φ 0 − 1.29n, which can be achieved with a focused laser beam, as described in Section 4.
The excitation at t = 0 is described by a Gaussian wave packet of the exciton states |ψ exc (k) , with the central wavevector k = 0. Fig. 1 shows that the entire wave packet acquires momentum during the external perturbation pulse (left panels). This is manifested as a phase variation in the coordinate representation, and as a shift of the central momentum in the k-representation. After the external perturbation is gone, the wave packet does not evolve in the k-representation and moves with the acquired uniform velocity in the coordinate representation.

Focusing of a delocalized excitation
In order to achieve full control over excitation transfer, it is desirable to find a particular phase transformation that focuses a delocalized many-body excitation onto a small part of the lattice, ideally a single lattice site. In optics, a thin lens focuses a collimated light beam by shifting the phase of the wavefront, thus converting a plane wave to a converging spherical wave. Similarly, a phase kick can serve as a time domain "lens" for collective excitations: an excitation initially delocalized over a large number of monomers can be focused onto a small region of the array after some time. By analogy with optics, a concave symmetric site-dependent phase Φ(n) applied simultaneously to all monomers may turn a broad initial distribution C n (t = 0) into a narrow one.
The dynamics of the excitation state in the lattice is determined by the time dependence of the coefficients C n (t) in Eq. (2). In order to find the expression for C n (t), we expand the amplitudes at t = 0 in a Fourier series and apply the propagator e −iE(q)t/ to each q-component with E(q) representing the exciton energy given by Eq. 14 GHz, which corresponds to LiCs molecules trapped on an optical lattice with lattice constant a = 400 nm and subjected to a homogeneous DC field of 1 kV/cm directed perpendicular to the intermolecular axis (for details, see Section 4). The kicking potential leading to a phase transformation Φ n Φ 0 − 1.29n is provided by a λ = 1064 nm Gaussian laser beam, with the propagation direction along the array axis, focused to 5 µm, with the intensity at the focus equal to 10 7 W/cm 2 . The laser pulse is on between 0 and 3 µs. The molecules are placed on the beam axis with the first molecule 5 µm away from the focus.
representation then yields where Φ(n) is a site-dependent phase applied at t = 0, as described in the previous section. Note that the phase Φ(n) does not have to be applied instantaneously. The phase φ(n, τ ) can be applied continuously over an extended time interval as long as the accumulated phase gives the desired outcome T 0 φ(n, τ )dτ = Φ(n). As Eq. (7) shows, the focusing efficiency is determined by the phase transformation and the shape of the dispersion curve E(k). Given the cosine dispersion of excitons (3), is it possible to focus a delocalized excitation onto a single lattice site? To answer this question, we assume that C n=n 0 (t = 0) = 1 and apply Eq. (7) to calculate the coefficients C m (t) at t = −τ . Using the expansion of an exponent in Bessel functions e −ia cos x = n e i(x−π/2)n J n (a) (8) and the orthonormality of the Bessel functions we find that the wave packet (2) with the expansion coefficients focuses, upon coherent evolution, in time τ on a single site n 0 . This shows that a phase transformation alone is, generally, not sufficient to create a collective excitation state that focuses onto a single lattice site. The best focusing must involve both the phase and amplitude modulations, which may be difficult to realize in experiments. A simpler procedure can be implemented if the phase transformations are restricted to a particular part of the exciton dispersion. From wave optics, waves with quadratic dispersion can be focused, while those with linear dispersion propagate without changing the wave packet shape [29,30]. It is this interplay of the quadratic (at low k) and linear (at k ≈ ±π/2a) parts of the cosine-like exciton dispersion (3) that precludes perfect focusing of a general collective excitation. In order to avoid the undesirable amplitude modulations, it may be possible to focus delocalized excitations by a phase transformation that constrains the wave packet (4) to the quadratic part of the dispersion E(k). For such wave packets, adding a quadratic phase Φ(n) = Φ 0 (n − n 0 ) 2 must lead to focusing around site n 0 . Below we illustrate the effect of the quadratic phase transformation for two types of initial states.
First, consider a broad Gaussian wave packet (2) a is the initial width. The corresponding width in the wave vector space is given by σ k = 1/σ x . The application of an inhomogeneous phase Φ(n) = Φ 0 (n − n 0 ) 2 at t = 0 results in additional broadening of the initial state, and the total width of the wave packet in the wave vector space with the account of the phase-induced contribution becomes [29,30] By analogy with optics, one should expect better focusing with larger Φ 0 (the width of the wave packet in real space is σ x (Φ 0 ) = 1/σ k (σ x , Φ 0 )). However, large values of Φ 0 may take the wave packet outside the quadratic part of the dispersion, impeding the focusing. To find the optimal phase Φ * 0 that keeps the wave packet within the quadratic dispersion while focusing it, we use the condition ∆ k = aσ k 1, which yields Φ * 0 = ±a/2σ x for the optimal focusing. At time the wave packet is most focused and has a width For the time t * in Eq. (12) to be positive, α and Φ * 0 must have the same sign. Therefore, a convex quadratic phase profile Φ(n) with Φ 0 > 0 must focus collective excitations in a system with repulsive couplings between particles in different lattice sites (α > 0), and a concave quadratic phase profile Φ(n) with Φ 0 < 0 must focus excitations in a system with attractive couplings (α < 0).
Second, consider a completely delocalized excitation (2) where Θ(z) = 1 if z is true and zero otherwise. In order to derive Eq. (14), we used the approximate equality obtained by approximating the error function of a complex argument Erf( √ ix) by the sign function, which is accurate for large argument x.
At time t * = 1/4αΦ 0 , the terms quadratic in q in Eq. (14) are canceled, and the sum over q reduces to a delta-function, if the summation limits are from −π/a to π/a. Therefore, the where ν k is the index of the initial wave vector k = 2πν k /N a, quantized due to the discreteness of the lattice. According to Eq. (14), the dimensionless width of the wave packet in the wave vector space is ∆ k (Φ 0 ) ≡ aσ k (Φ 0 ) ≈ 2N Φ 0 . When Φ 0 = π/N , the wave packet spreads over the entire Brillouin zone, including the linear parts of the exciton dispersion. Using Eq. (15) we find that for an arbitrary value of ∆ k (Φ 0 ), the site amplitudes at the time of focusing t * = 1/4αΦ 0 are In order to keep the linear part of the dispersion spectrum unpopulated, we choose the optimal focusing phase Eqs. (13) and (16) are valid for a many-body system with nearest neighbor interactions only. In most physical systems, the energy dispersion is modified by longrange couplings. In order to confirm that the above predictions are also valid for systems with long-range interactions and illustrate the focusing of delocalized excitations, we compute the time evolution of the wave packets by solving the wave equation numerically for a system with long-range dipole-dipole interactions. Figure 2 illustrates the focusing dynamics of a completely delocalized excitation (panels a and b) and a broad Gaussian wave packet (panels c and d) in a system with all (first neighbour, second neighbour, etc.) couplings explicitly included in the calculation. The results show that the collective excitations can be focused to a few lattice sites. The role of the long-range coupling will be explicitly discussed in Section 6.
The focusing scheme demonstrated above can be generalized to systems of higher dimensionality. To illustrate this, we repeated the calculations presented in Figures  2c and 2d for a delocalized excitation placed in a square 2D lattice with an external potential that modulates the phase as a function of both x and y. Figure 3 shows the focusing of an initially broad wave packet onto different parts of a 2D lattice induced by the quadratic phase transformation Φ(x, y) where n x and n y are the lattice site indices along the x and y directions. The calculations include all long-range couplings as in Figure 2. The comparison of Figures 2(c,d) and 3 illustrates that the focussing efficiency in 2D is greater. The results also demonstrate that the delocalized excitations can be effectively focused on different parts of the lattice simply by varying the reference site (n x 0 , n y 0 ) in the phase transformation.

Controlled excitations of ultracold atoms and molecules
The techniques proposed in Sections 2 and 3 can be realized with ultracold atoms or molecules trapped in an optical lattice in a Mott insulator phase [31]. There are three general requirements that must be satisfied: • (i) The time required for a simple phase transformation must be shorter than the spontaneous decay time of the excited state.
• (ii) The overall coherence of the system must be preserved on the time scale of the excitonic evolution in the entire array, set by K /α, where K is the number of monomers participating in the dynamics of the collective excitation.
• (iii) The lattice constant must be large enough to allow considerable variation of the external perturbation from site to site.
Optical lattices offer long coherence times (> 1 sec) and large lattice constants (> 400 nm) [32]. The lifetime of the collective excitations depends on the internal states of the particles used in the experiment and the momentum distribution of the excitonic states in the wave packet (4).
For ultracold alkali metal atoms in an optical lattice, an optical excitation may generate collective states (2), as discussed in Refs. [33,34,35,36,37,38,39,40,41]. The lifetime of these excited states is limited by the spontaneous emission of the electronically excited atoms and is in the range of 10 -30 ns. However, the collective excited states can be protected from spontaneous emission if the wave vector range populated by excitons in the wave packet (4) is outside of the light cone, so that k > ∆E e−g / c [26]. These states do not readily radiate as energy and momentum conservation cannot be simultaneously satisfied for them. The emission of photons may occur on a much longer time-scale at the array boundaries or due to perturbations breaking the translational symmetry. Due to the same conservation laws, single-photon excitation of atomic ensembles always generates excitons with k ≈ 0. Once these excited states are created, the phase-kicking technique introduced in Section 2 can be used to shift the excited states in the wave vector space away from k = 0 (cf. Figure 1) and thus protect the excited states from fast spontaneous decay. This phase transformation can be induced by a pulse of an off-resonant laser field E AC , detuned from the e ↔ g resonance by the value δω, leading to the AC Stark shift (see e.g. Ref. [29]) where V eg is the matrix element of the dipole-induced transition. By choosing V eg = 1 a.u., δω = 3V eg , and the laser intensity I = 5 × 10 10 W/cm 2 , we obtain that the shift φ = ∆E AC × T pulse = π can be achieved in less than 1 ns.
This shift brings a wave packet initially centered at k = 0 to the "dark" edge of the Brillouin zone, where the dispersion of excitons is still quadratic and all the focusing schemes discussed in Section 3 can be applied. Another phase transformation can bring the excited state back to the k ≈ 0 region, where it can be observed via fast spontaneous emission. The experiments with ultracold atoms have demonstrated the creation of a Mott insulator phase with the lattice filling factor reaching 99 % [42,43,44]. The phase transformations proposed here can be used to stabilize excitonic states in ultracold atomic ensembles against spontaneous emission for multiple interesting applications [33,34,36,38,40].
The spontaneous decay problem can be completely avoided by using rotational excitations in an ensemble of ultracold polar molecules trapped in an optical lattice. The rotational states are labeled by the quantum number of the rotational angular momentum J and the projection M J of J on the space-fixed quantization axis Z. We choose the rotational ground state |J = 0, M J = 0 as |g and the rotational excited state |J = 1, M J = 0 as |e . The state |J = 1, M J = 0 is degenerate with the states |J = 1, M J = ±1 . This degeneracy can be lifted by applying a homogeneous DC electric field, making the |g and |e states an isolated two-level system. The molecules in different lattice sites are coupled by the dipole-dipole interaction V dd (n − m). The magnitude of the coupling constant α(n − m) = e n , g m |V dd (n − m)|g n , e m between molecules with the dipole moment 1 Debye separated by 500 nm is on the order of 1 kHz [11]. Due to the low value of ∆E e−g , the spontaneous emission time of rotationally excited molecules exceeds 1 second.
For molecules on an optical lattice, one can implement the phase kicks by modifying the molecular energy levels with pulsed AC or DC electric fields. The rotational energy levels for 1 Σ molecules in a combination of weak AC and DC electric fields are given by [45] where B is the rotational constant, G(0, 0) = −1/3, G(1, 0) = 1/5, F (0, 0) = −1/3, F (1, 0) = −3/5, E AC is the envelope of the quickly oscillating AC field, α and α ⊥ are the parallel and perpendicular polarizabilities and µ is the permanent dipole moment of the molecule. The momentum shift of the exciton wave packets can be achieved by applying a time-varying DC electric field E(t) = E * + E(n) sin 2 (πt/T ), where E(n) is linear with respect to n. Assuming that E(n) = (n − n 0 )A and E(n) E * , and using Eqs. (5) and (18), gives δ = 4AE * µ 2 T /15 Ba. We have confirmed this result by a numerical computation showing that for LiCs molecules in an electric field of E * = 1 kV/cm, an electric field pulse with A = 7.434 × 10 −4 kV/cm and T = 1 µs results in a kick of δ = π/2a, bringing an excitonic wave packet from the k = 0 region to the middle of the dispersion zone.
Non-adiabatic control of quantum energy transfer in ordered and disordered arrays 13 An alternative strategy is to use a pulse of an off-resonance laser field, as for atoms. The phase transformations can be induced by a Gaussian laser beam with the intensity profile where I 0 is the light intensity at the beam center, r is the radial distance from the center axis of the beam, z is the axial distance from the beam center, z R = πw 0 /λ is Rayleigh range, w 0 is the beam waist and λ is the wavelength. With the 1D molecular array arranged along the z-axis, the laser field intensity can be made to vary linearly along the array, where I c is the intensity at the center of the wave packet. This can be achieved if x a 0.5z R , where z 0 is the distance between the center of the wave packet and the beam center, and σ (2d) x is the width (in the coordinate representation) of the two-dimensional wave packet. Using Eqs. (5), (18) and (19), we estimate the momentum kick by such a pulse as δ = − √ 3T I 0 (α − α ⊥ )/80z R . The results presented in Figure 1 were obtained for a 1D array of LiCs molecules on an optical lattice with a = 400 nm and the external perturbation given by the laser field pulse (20) with parameters I c and I 1 derived from Eq. (19) with z 0 = 45 µm and z R = 73.8 µm. The numerical results deviate from the analytical prediction for δ by less than 7 %.
The Gaussian intensity profile (19) can be used also to implement the quadratic phase transformations needed for focusing of collective excitations. To achieve this, a 2D molecular array must be arranged in the z = 0 plane, with the x-axis defined to be along the polarization direction of a linearly polarized field. If the dimension of the molecular array is smaller than one third of the beam waist, the Gaussian intensity profile in Eq. (19) can be approximated by This is a concave quadratic intensity profile which can be used to focus a wave packet in a system with negative couplings α (see Section 3).

Control of energy transfer in dipolar systems
Dipolar interactions play a central role in the study of long-range interaction effects using ultracold systems [46]. While, in general, the coupling constant α in Eq. (1) can be determined by a variety of interactions, the dominant contribution to α for atoms and molecules on an optical lattice is determined by the matrix elements of the dipoledipole interaction. It is therefore particularly relevant to discuss the specifics of energy transfer in systems with dipolar interactions.
The dipolar interactions are long-range and anisotropic. The long-range character of the dipolar interactions manifests itself in the modification of the exciton dispersion (3). While Eq. (3) is valid for a system with nearest neighbour couplings only, higherorder couplings in the case of α(n − m) ∝ 1/(n − m) 3 modify the exciton dispersion leading to a cosine-like, but non-analytic dispersion relation, both in 1D and 2D. To investigate the effect of this nonanalyticity in dispersion curve, we have performed a series of calculations with the long-range couplings neglected after a certain lattice site separation n − m for the 1D system. The results become converged (to within 0.2 %) when each molecules is directly coupled with 20 nearest molecules. While the calculations with only the nearest neighbor couplings are in good agreement with the analytical predictions given by Eqs. (12) and (13), the full calculations reveal that longrange couplings somewhat decrease the focusing efficiency. The long-range couplings also decrease the focusing time, by up to a factor of 2. The dynamics of collective excitations leads to interference oscillation patterns clearly visible in panels b and d of Figure 2. These oscillations are much less pronounced when all but nearest neighbor couplings are omitted. The numerical results of Figures 1 -3 are particularly important because they demonstrate that the phase transformations introduced in the present work are effective for systems with dipolar interactions.
The anisotropy of the dipolar interactions can be exploited for controlling energy transfer in dipolar systems by varying the orientation of a dressing external DC electric field. For example, for polar molecules on an optical lattice, the matrix elements α(n − m) = e n , g m |V dd (n − m)|g n , e m depend not only on the choice of the states |g and |e , but also on the magnitude and orientation of an external dc electric field [13,11]. Since the value of α determines the exciton dispersion (3), the exciton properties can be controlled by varying the angle θ between the intermolecular axis and the applied DC field. This is illustrated in Figure 4.
The calculations presented in Figure 4 are for a 1D array of LiCs molecules in a lattice with a = 400 nm. As before, |g is the absolute ground state of the molecule and |e is the rotationally excited state that adiabatically correlates with the rotational state |J = 1, M J = 0 in the limit of vanishing electric field. The upper panel of Figure  4 shows that the angle θ between the electric field vector and the molecular array axis determines the sign and magnitude of α, and therefore the shape of the dispersion curve. This enables control over the sign and magnitude of the group velocity of an excitonic wave packet containing contributions with k = 0. Dynamically tuning θ, one can propagate a localized excitation to different parts of the lattice, as shown in Figure  4b.
In a 2D lattice, the intermolecular interactions depend on an additional azimuthal angle φ that describes the rotation of the electric field axis around the axis perpendicular to the lattice. The numerical calculations presented in Figure 5 show that the energy flow in two dimensions can be controlled by varying both θ and φ. In addition to the phase transformation discussed earlier, this allows for a dynamical energy transfer in quantum many-body systems with anisotropic interparticle interactions.

Energy transfer in the presence of vacancies
While experiments with ultracold atoms have produced a Mott insulator phase with 99% of lattice sites filled [42,43,44], the latest experiments with molecules yield lattice-site populations about 10% [31]. Multiple experiments are currently underway to produce a Mott insulator phase of polar molecules with close to the full population of the lattice. However, lattice vacancies may be unavoidable in the best experiments. In this section, we examine the effect of vacancies on the possibility of focusing collective excitations to a desired region of the lattice by the phase transformations discussed in Section 3. For concreteness, we perform calculations for the system described in Section 4, namely a 2D array of LiCs molecules on a square optical lattice with a = 400 nm.
To explore the effect of vacancy-induced interactions, we performed simulations respectively. The initial wavepacket is a 2D Gaussian distribution centered around ak x = ak y = π/2 and has a width of ∼60 lattice sites in coordinate space. The magnitude of the DC field is fixed to 6 kV/cm while its direction is changing. The calculations are done for a 2D array of LiCs molecules in a lattice with a = 400 nm.
for different vacancy numbers using the same parameters for molecule-field and intermolecular interactions as in the calculations presented in Figure 3b. For each vacancy concentration, we carried out 48 calculations with random distributions of empty lattice sites. The quadratic phase transformations are applied, as described in Section 3, in order to focus the collective excitation at time t * to the molecule in the middle of the 2D array. Vacancies disturb the translational symmetry of the system and produce an effective disordered potential that tends to localize collective excitations [47]. Because the natural time evolution of the wave packet in a disordered potential may lead to enhancement of the probability in certain regions of the lattice, it is necessary to distinguish the effect of the vacancy-induced localization and the effect of the focusing phase transformation. To quantify these two effects, we define two factors: the enhancement of the probability at the target molecule with respect to the initial value, and the ratio of the probability to find the excitation on the target molecule with (p ) and without (p) the focusing phase transformation, The time t * is the focusing time predicted in Section 3 for an ideal, vacancy-free system. The quantity η illustrates the actual enhancement of the probability to focus a collective excitation, while the quantity χ illustrates the effect of the focusing phase transformation. Figure 6 presents the values of η and χ as functions of the vacancy concentration. It illustrates two important observations. First, the disorder potential with vacancy concentrations > 20 % renders the phase transformation uneffective. In the presence of strong disorder, the dynamics of the system is entirely determined by the disorder potential and the energy transfer becomes highly inefficient (however, see Section 7). On the other hand, vacancy concentrations of less than 10 % appear to have little effect on the efficacy of the focusing phase transformation.
Our calculations indicate that the focusing time may be somewhat modified by the disorder potential, even if the concentration of vacancies is less than 10 %. Figure 7 depicts the excitation wave functions at the time of the maximal enhancement on the target molecule, chosen as molecule (71,71). Figure 7 shows that despite the presence of multiple vacancies, the focusing transformation enhances the probability to find the excitation on the target molecule by 16 times.

Focusing in the presence of strong disorder
Although the focusing method demonstrated in Sections 3 and 6 appears to be robust in the presence of a disorder potential induced by a small concentration of vacancies, it is important for practical applications to also consider controlled energy transfer in quantum arrays under a strong disorder potential. To consider focusing in a strongly disordered system, we employ an analogy with the "transfer matrix" methods for focusing of a collimated light beam in opaque medium [48,49,50,51,52,53,54,55,56,57].
In optics, a collimated laser beam passing through an opaque medium results in a random pattern of speckles arising from random scattering of light inside the medium [58]. Likewise, the random distribution of empty sites in an optical lattice with molecules scatters the exciton wavepackets, resulting in a completely random excited state. However, in optics, the randomness of the scattering centers inside the opaque medium can be compensated for by shaping the incident wavefront with a spatial light modulator such that the contributions from various parts of the medium can add constructively upon exit from the medium, producing a focus. We suggest that the same can be achieved with a many-body system on a lattice by separating the entire lattice into multiple blocks and applying proper phase transformations to those individual blocks.
The initial state for an ensemble of molecules on a lattice with multiple vacancies can be written as where |i = |e i j =i |g j (25) and the indexes i and j run over all occupied sites. After a long evolution time T , the probability amplitude for the excitation to reside on a particular target molecule is given by where U o,i (t) = o| exp[−iH exc t]|i is a matrix element of the time evolution operator. In a disordered system, the transfer coefficients U o,i are not a-priori known and depend on the disorder potential. The phasors c oi (T ) have quasi-random amplitudes and phases. While the amplitude of each phasor cannot be controlled experimentally, their phases are controllable via the phases of the coefficients c i at t = 0, which can be tuned using the phase-kicking transformations introduced above. To achieve the highest probability at the target molecule, it is necessary to ensure that the contribution c oi = U o,i (T )c i (t = 0) from every site i has the same phase so that they add up constructively. In a practical implementation, it may be difficult to control the phase of each molecule in each individual site. It may be more desirable to work with blocks of several lattice sites. Assuming that the entire array of molecules can be divided into M blocks, each containing many molecules, and that the blocks can be perturbed individually, the excitation probability amplitude at the target molecule at time T is where This equation implies that the contributions from different blocks can be made to interfere constructively by adding a phase exp(−iφ γ ) to each occupied site in block γ. For M blocks in the array and quasi-random evolution matrix, simply setting all the phases equal must lead to ∼ M -fold increase of the excitation probability at the target molecule, as compared to a sum of M quasi-random phasors in Eq. (27) [48]. Similarly to optics, the phases −φ γ which must be added in each block, can be found experimentally provided that the same (or similar) realization of disorder persists in a series of trials. A straightforward optimization would scan through the strengths of phase kicks applied to different blocks. In each experiment one would measure the excitation probability at the target molecule |c o (T )| 2 , e.g. via resonance fluorescence from the target molecule at the end of the experiment. More sophisticated optimization techniques, aimed at fast focusing multi-frequency light in optical systems, are currently under rapid development [49,50,51,52,53,54,55,56,57].
For a proof-of-principle calculation, we consider a 2D lattice of size 101×101 with 60% of sites vacant and each non-vacant site occupied by a single LiCs molecule. Due to time reversibility of the time evolution operator U (T ), The matrix element U γ j,o (−T ) can be calculated by performing a backward time propagation starting from a local excitation at site "o" and calculating the coefficient c j (t) at time −T . Alternatively, one can propagate the evolution equations forward in time, finding c j (T ): Since the Hamiltonian (1) is real, its eigenfunctions are real, and the evolution matrix U is symmetric, U o,j = U j,o . Thus we find since c o (0) = 1 and all other coefficients are zero. For a completely delocalized initial state, we assume that all coefficients in Eq. (24) are equal, so that the phases φ γ required for block γ are where the index j runs over all occupied sites in block γ. Figure 8 shows that this choice of phases leads to effective focusing of the collective excitation in a strongly disordered system. To illustrate the efficiency of the focussing method described above, we have carried out a series of calculations with different vacancy concentrations. For each vacancy concentration, we performed 48 calculations with random distributions of empty lattice sites. The phase transformations are calculated individually for each random distribution of vacancy sites as described above. The results are shown in Figure 9. As can be seen, the transformations proposed above are effective for vacancy concentration < 70%. At higher concentrations of vacancies, the excited states become strongly localized and immobile. The focusing efficiency at vacancy concentrations 10% and 20% appears to be higher than that in the absence of vacancies, which we attribute to the effect of the boundaries. (Color online) Enhancement factors η (red symbols) as a function of vacancy percentage in a 2D lattices. η has the same definition as in Eq. (22) except the time t * is arbitarily chosen to be 4 ms here. The error bars are for 95% of confidence interval.

Conclusion
We have proposed a general method for controlling the time evolution of quantum energy transfer in ordered 1D and 2D arrays of coupled monomers. Any elementary excitation in an aggregate of coupled monomers can be represented as a coherent superposition of Frenkel exciton states. We propose shaping the exciton wave packets using nonadiabatic perturbations that temporarily modulate the energy levels of the monomers leading to monomer-dependent linear phase transformation and a displacement of the wave packets in the wave vector representation. This, combined with the possibility of focusing a collective excitation on a particular part of the lattice by a quadratic phase transformation and with the directed propagation of collective excitations, allows for control of energy transfer in the lattice.
We have presented numerical calculations for an ensemble of polar molecules trapped on an optical lattice that demonstrate the feasibility of both momentum-shifting and focusing of collective excitations by applying external laser fields, with parameters that can be easily achieved in the laboratory. We also investigated the effect of disorder potential arising from incomplete population of the lattice. Our results show that the phase transformations leading to focusing of collective excitations on different regions of a 2D lattice remain effective in the presence of vacancies with concentrations not exceeding 10 %. For systems with larger concentrations of vacancies and affected by strong disorder potentials, we propose an alternative procedure based on engineering constructive interference of the wave function contributions arising from difference parts of the lattice.
The momentum-shifting technique proposed here can be used to protect collective excitations of ultracold atoms from spontaneous emission. The spontaneous decay processes, which in the case of an ordered many-body system must satisfy both the energy and wave vector conservation rules, can be restricted by shifting the exciton wave packets to a region of the dispersion curve, where the wave vector conservation cannot be satisfied. If performed faster than the spontaneous emission time, such phase transformations should create collective excitations with much longer lifetimes, which opens a variety of new applications for ultracold atoms on an optical lattice.
Control over excitation transfer is needed for creating networks of quantum processors where information is transmitted over large distances with photons and stored in arrays of quantum monomers via one of the quantum memory protocols [59]. Momentum kicking can be used for information transport within a single array. Focusing excitonic wave packets enables local storage of information, while directed propagation combined with controlled interactions of multiple excitons [13] or excitons with lattice impurities [60] may be used to implement logic gates. Controlled energy transfer in molecular arrays may also be used for the study of controlled chemical interactions for a class of reactions stimulated by energy excitation of the reactants. Directing quantum energy to a particular lattice site containing two or more reagents can be used to induce a chemical interaction [61], an inelastic collision or predissociation [62] with the complete temporal and spatial control over the reaction process.
Finally, the present work may prove to be important for simulations of open quantum systems. We have recently shown [11,12] that the rotational excitations of ultracold molecules in an optical lattice can, by a suitable choice of the trapping laser fields, be effectively coupled to lattice phonons. The exciton -phonon couplings can be tuned from zero to the regime of strong interactions [11,12]. The possibility of shaping (accelerating, decelerating and focusing) collective excitations as described in the present work combined with the possibility of coupling these excitations to the phonon bath opens an exciting prospect of detailed study of controlled energy transfer in the presence of a controllable environment. Of particular interest would be to study the effect of the transition from a weakly coupled Markovian bath to a strongly coupled non-Markovian environment on energy transfer with specific initial parameters.
We note that the effect of site-dependent phase transformations on quantum transport was independently considered in Ref. [63] from the point of view of timereversal symmetry breaking. The authors of Ref. [63] propose an experimental realization based on ions in a linear Paul trap. Their method relies on the possibility of tuning time-dependent phases, leading to new effects. The present work and Ref. [63] should be considered complementary.