Multiple scattering induced negative refraction of matter waves.

Starting from fundamental multiple scattering theory it is shown that negative refraction indices are feasible for matter waves passing a well-defined ensemble of scatterers. A simple approach to this topic is presented and explicit examples for systems of scatterers in 1D and 3D are stated that imply negative refraction for a generic incoming quantum wave packet. Essential features of the effective scattering field, densities and frequency spectrum of scatterers are considered. Additionally it is shown that negative refraction indices allow perfect transmission of the wave passing the ensemble of scatterers. Finally the concept of the superlens is discussed, since it is based on negative refraction and can be extended to matter waves utilizing the observations presented in this paper which thus paves the way to 'untouchable' quantum systems in analogy to cloaking devices for electromagnetic waves.

Negative refraction of electromagnetic (EM) waves was first discovered by V. G. Veselago in1 and in particular it has been noted that at the negative refractive material’s surface waves entering the material are bent towards negative angles. Decades later meta-materials with this property were developed and they were found to give rise to the phenomenon of sub-wavelength localisation (superlensing). In addition to perfect localisation of the EM waves those materials imply complete transmission of the incoming wave passing through those so called superlenses234. For a long time negative refraction was considered a mathematical curiosity rather than a feasible physical phenomenon shown in experiments, it is only due to the invention of artificial meta-materials that the necessary conditions could be provided experimentally345 and by the means of those materials sub-wavelength localisation was shown6. Due to the huge impact of sub-wavelength localisation of EM waves a natural question to ask is: Are negative refractive materials feasible for matter waves as well and what is necessary to refract matter waves in a similar way as EM waves? To give answers to those questions in this work I will utilize multiple scattering theory, which has been widely applied to derive refraction indices for matter waves78910. It is worth noting that multiple quantum scattering theory is an application of the full quantum many-body problem for fixed particle numbers that allows a more practical viewpoint when single particle scattering on a cloud of potentially randomly positioned scatterers is of interest10. The refraction indices n derived so far by this approach are however approximations and their validity is particularly confined to the case . So one may ask: Are matter waves restricted to refract positively, perhaps due to a property of the governing many-body Schrödinger equation? In fact, this work discusses the feasibility of negative refraction in a way that is consistent with the quantum many-body problem and conditions for negative refraction are stated, i.e. scenarios of ensembles of scatterers implying negative refraction. Finally based on the concept of negative refraction the concept of superlenses34 for matter waves is discussed. Although the presented theory is new I point out that the topic of negative refraction for matter waves was considered before in11 and a specific proposal in terms of the single particle Schrödinger equation including external magnetic fields was regarded as a candidate for negative refraction, while the results presented in11 were found to be inaccurate12. In addition only recently Veselago-type lensing was successfully implemented utilizing ultra cold atoms confined to optical lattices by switching between the positive and negative energy branches of the many-body states13. In contrast to the proposal11 and the realisation13 this work follows a different path showing that negative refraction is a property of certain closed quantum many-body systems (with fixed particle numbers) without the necessity to apply external fields, feasible only due to a proper arrangement and choice of scatterers.

Multiple scattering theory

The fundamental multiple scattering equations for quantum matter waves were first introduced by L. L. Foldy and extend the classical treatment based on the Boltzmann integro-differential equations when interference effects of the scattered particle become important14. The set of equations usually used to describe the scenario for point scatterers where the strength of the scattered wave from a scatterer is proportional to the external field acting on it read as follows10141516,

Here the l.h.s. is the wave function generated by multiple scattering processes that are induced by a system of N randomly arranged particles/obstacles and

describes the local field which generates the scattered wave at positions . is the incident wave, i.e. the wave without the effects of scatterers, is the local field at where the ith scatterer e.g. a noble gas atom is located at a fixed position - importantly can be approximated by when all multiple scattering processes at scatterer i are neglected14, i.e. only the first order scattering from each scatterer i of the ensemble of multiple scatterers is considered. is the scattering amplitude/coefficient that relates the scattered wave to the scatterer and which in general is a function of the wave vector and position in space. In 3D the behavior of the particle close to the ith scatterer is assumed to be101415 , which is the Green’s function describing the propagation of the wave scattered from the ith particle between the scattering particles. This form assumes a spherically symmetric scattering process (s-wave scattering)14. The corresponding Green’s function for a 2D system is the Hankel function and for the 1D system it is a plane wave, which will be introduced later. It should be noted that from (1) and (2) a plethora of derivations have been made particularly for refractive indices of matter waves and comparison with experiments show excellent agreement with data e.g. in experiments with monatomic sodium gases78910141516.

Statement of the problem

Let us assume and to approximately resemble plane waves, which can be regarded a valid assumption for waves far away from the scattering processes and they form a basis for general localised wave packets. Given that and utilizing the mathematical notions from above we can state the scenario of interest as follows: For

do amplitudes , local fields and configurations of scatterers exist that extinguish the initial field and generate a new field which moves in opposite direction? Thus formally we are interested in states where the incoming wave with wave vector k implies an outgoing wave with wave vector K so that

This equation defines the refraction index for the matter wave n10, which by definition of the problem is negative. Note that although at first sight negative refraction with could be confused with the properties of a perfect mirror with reflection coefficient , but in contrast waves refracted negatively occupy the space ‘behind’ the scatterers as well. As shown later there even is no loss of amplitude of the scattered wave after passing the scatterers for certain stated parameters. The remainder of this work is dedicated to show which local fields and amplitudes are necessary to imply the refractive indices to be negative and of absolute value 1 and further properties of those systems. The motivation to investigate the case of negative refraction stems from2, where meta-materials allow sub-wavelength localisation for electromagnetic matter waves and the argumentation is based on the extinction theorem for matter waves14. Mathematically the l.h.s. of (3) consists of a set of homogeneous linear algebraic equations and to make statements about etc. approximations can be made. Furthermore as plane waves are simplifications we will also discuss wave packets, i.e. superpositions of plane waves, moving in opposite directions, such that for each component of the wave packet’s spectrum (4) is satisfied.

Simplified state equation

Alternatively to solving eq. 1 directly one can proceed as follows. As discussed in detail by14 eq. 1 can be simplified when considering the statistical average over all possible random positions of scatterers, which is our first assumption, while the case of fixed positions is trivially included. Then by further supposing the scatterers not to have internal degrees of freedom and to occupy a finite volume V one effectively considers

Here is the scatterer induced field acting on the ith scatterer, which is averaged over all possible configurations of all the other scatterers14. This field can be approximately replaced by an averaged field, which is a fair approximation when the number of scatterers N is large14 and supposing that we substitute . So (5) reduces to the simplified state equation

where is the coherent wave, i.e. the average matter wave function for all possible configurations of scatterers. denotes the average local field induced by the scatterers and the average amplitude - both functions are given by1415

and is the number density of scatterers. As for scatterers close to each other compared to other length scales of the system, those sums can be approximated as well as integrals.

Positive refraction

While this paper follows the aim to derive wave functions for negative refraction indices we give a short overview how positive refraction can be derived easily from (6). The macroscopic one-body Schrödinger equation equation for the average wave function is15

where E is the incident particle energy and m its mass. In terms of eq. 9 (6) is a solution for the optical potential

when assuming i.e. to be proportional by a number c and for the incident wave energy given by 15. From (10) one derives the simple relation and thus a refraction index , with denoting the wave vector of the generated wave. This derivation of the index of refraction, however cannot address the case of negative refraction as only a formula for the squares of the wave vectors is given. In order to extend the discussion to negative refraction indices we may proceed as follows.

Wave packets

To develop the formalism our starting point is a general scalar matter wave function describing a quantum particle in D spatial dimensions

with being the distribution function of different frequency components of the wave. Correspondingly the de Broglie and Einstein relations which associate the wave properties with that of a massive ‘projectile’ are17

Here p denotes the momentum of the projectile, v its velocity, v the frequency and E the energy of the particle to be scattered. We consider the pair of initial and outgoing wave formally given by

when incoming − and outgoing waves + are stationary.

Results

Negative refraction in 1D

We are interested in the scenario where the initial wave packet is moving in opposite direction as the outgoing wave packet and we first confine our considerations to . To model two counter-propagating waves we consider by using eq. 6 the expression

The superscript of the Green’s function indicates its spacial dimension. Here we have introduced the abbreviation for the unknown , which we determine as follows. Let us rewrite the r.h.s. of (15) as

and assume to be compactly supported on V, which is a fairly reasonable statement in physical terms particularly as the density of scatterers is only supported in V. Now for a 1D system the Green’s function of the inhomogeneous Helmholtz equation

corresponding to isotropic point scatterers at is

which provides the behaviour of the particle close to the scatterer with incident wave number . Clearly linear combinations of Green’s functions are Green’s functions as well and so we obtain by the above solutions e.g.

i.e. a stationary wave solution and similarly one gets .

By combining (15) and (16) we obtain the condition for negative refraction in 1D when writing the two complex exponentials as sine,

We assume the amplitude f to be independent of , i.e. the strength of the scattered wave is only proportional to the Green’s function and the number density of scatterers. In principle scatterers that obey a position dependency of the scattered wave beyond that of the Green’s function of each scatterer, i.e. in f can be simplified by effectively choosing the number distribution of scatterers accordingly. As is a product of the density of scatterers and the local effective field there is in principle freedom to generate a variety of experimentally by arranging scatterers appropriately, which fixes . It is thus the choice of the experimenter first to use scatterers that indeed obey a position independency in f and secondly to effectively compensate this dependency by the arrangement of scatterers. Note that the nuclear scattering amplitude has an expansion of the form15

with the bound scattering length .

Examples for 1D scattering

Example 1

In the following we will give explicit examples of scenarios yielding negative refraction indices: For the simplest example we set

resembling scatterers at the position with effective field . Further we choose the stationary Green’s function for the scatterers to be

which can be constructed from the elementary solutions of the Helmholtz equation with point source by taking the superpositions

So the material consists of scatterers inducing an anisotropic wave (23). The frequency dependent amplitude is set

Thus we have given an example for and a condition on the frequency spectrum of the scatterers, which implies a negative refractive material for the generic incoming wave function defined by (20). Wave packets can be obtained by integrating (20) over k with appropriate weighting factors. On the other hand for given f we can model wave functions via g that satisfy (25).

Example 2: Anisotropic scatterers

The above Green’s functions can be generalised to anisotropic scatterers, so that scattering events depend on the direction of the incoming wave. In 1D the corresponding anisotropic Helmholtz equation for each scatterer is generalised to1819

We assume

with solutions

where are chosen to satisfy continuity. Setting we get a decaying solution for . For two different species of scatterers positioned at the same points in space (each Green’s function approximately being solutions of a Helmholtz equation) we get the superposed effective Green’s function

Utilizing two species of scatterers is applicable for high as this decreases the relevance of the decaying part of the wave functions and we formally arrive at the previous case. Thus we have generated a similar Green’s function as in (23), but by very different physics involving two different scattering profiles of two species of scatterers occupying the same space.

Extended scattering domain

While the δ distributed scatterers only resemble an approximation to realistic formations of scatterers, we extend the 1D analysis to fields,

where the plane wave part is due to the spatial dependence of the effective field and the remainder models an extended distribution of scatterers. Hence the l.h.s. (20) for a stationary Greens function solution (23) is given by,

Further

Again choosing the frequency spectrum f accordingly we can match the r.h.s. in (20) to obtain negative refraction for or for or for given f we confine the consideration to appropriately related g.

Negative refraction in 3D

Next we turn to the case of negative refraction in 3D materials. Green’s functions of the 3D Helmholtz equation for different frequencies are

with , called respectively the outgoing and ingoing spherical waves solutions. Linear combinations of the incoming and outgoing wave solutions for the 3D Helmholtz equation are e.g. the stationary waves

and

Similar as in the 1D case we set all scatterers to rest at the origin, but in contrast with a certain spatial distribution given by

while using the which neglects all additional multiple scattering processes at this particular scatterer14. W.l.og. we set and in addition we assume the scatterers to induce a Green’s function (35)

Starting from (6) we require as condition for negative refraction in 3D:

We assume . On the other hand we get using (36) and (37)

Next we apply the paraboloidal wave approximation, that is valid for . We are interested in the far field where , so that the approximation is satisfied. Thus we obtain

Choosing by comparison of (40) with (38) concludes the proof for the far field, i.e. in the limit.

Implicit negative refraction

Scattering in 1D

While the previous examples provide a general guide to induce negative refraction, we now turn to a specific example feasible in systems satisfying the stated assumptions. For the sake of simplicity we consider the 1D case. For each k there is a solution with no incoming wave from the right - a standard solution familiar in calculations for reflection and transmission from a generic scattering potential (with as modelling a single finite scatterer. While here we make specific choices note that the above general framework applies for an ensemble of many (generally randomly distributed) scatterers. Now, we consider an incoming wave

and consequently we have an outgoing wave

where the reflection probability is and the transmission probability - for unitary scattering the conservation of probability applies, i.e.

We assume scatterers to act on finite range, . Using this generic framework we explicitly construct an example where due to multiple scattering events and in contrast to perfect mirrors the solutions are not confined to the area in front of the perfect mirror.

Refracted waves in 1D

To further elucidate the feasibility of negative refraction in an ensemble of scatterers we proceed as follows. We consider the incoming wave moving from the l.h.s. towards the scatterers which then scatters on the meta-material spanning from a to c,

The induced wave function between two scatterers, say A and B at positions a and b, given there is another scatterer C at c with and subsequent scatterers and scattering events are neglected, is

where and with and are the corresponding transmission or reflection coefficients from the left or right of the three scatterer system. At the same time each pair of transmission and reflection coefficients has to conserve probability (43) and if the transmitted wave from the r.h.s. has to be taken into account by consider the effective incoming wave

To get the scattering behavior that induces a negatively refracted wave between A and B, i.e. , as implied in (3) we find using (45) and (46) the necessary conditions

and

Note that the setup trivially extends to N scatterers, denoting the jth scatterer. By letting the first at position and the last scatterer at we have generated a scenario consistent with (20) when choosing accordingly. Formula (45) is an approximation considering only the first and second order reflections and transmissions - iteratively one would obtain the exact form including all reflection and transmission processes.

Numerical example

Let us consider an electron of mass with velocity corresponding to a wave vector (12) . Furthermore we consider the three scatterers and specify their scattering properties by choosing , , , , , , , to satisfy (47) and (48). Then using the above characteristics we obtain the wave functions

where and note that the wave is independent of the actual positions a, b and c.

Complete tunnelling in negative refraction media

Next we analyse the transmittance of a matter wave passing a negative refraction index material. We assume the material with refraction index n to span from to . Therefore the wave function before, within and after the material is

correspondingly, where is in accordance to789. Here some attenuation σ of the wave within the material and a change in phase has been introduced and the transmission amplitude is again t. At the points and L continuity of the wave function requires

and continuity of its first derivative

Consequently the transmission amplitude is

with . Further using and writing we obtain the transmittance

In Fig. 1 the transmittance is shown as a function of the real part of the refraction index, i.e. . For the given parameters we observe that a negative refraction medium for matter waves allows their amplification within the material in analogy to the concept of electromagnetic waves3. However the conservation of probability provides a natural upper bound to the transmittance in a single-body Schrödinger equation picture. Recall the relation , where A corresponds to the absorption due to material properties. In case that necessarily and we consider the scenario of perfect transmission of a negatively refracted wave through the medium. In contrast in positive refraction indices materials usually waves decay - this is particularly shown by the graph of Fig. 1 for . Hence, in this setting complete transmission of a matter wave is feasible only within a negative refraction material. Furthermore it is known that negative refraction of waves allows localisation to a single point given that as it has been shown in346 for electromagnetic waves, which is a consequence of the modified wave propagation within those meta-materials and applies in one-to-one analogy. The reason can be understood by the argument that the optical path through a negative refraction material is zero and so the initial point source is mapped on itself. I refer to34 and references therein for a complete discussion of this topic.

Figure 1

Transmittance T as a function of Re(n) for decay σ = 0.95, c = 1.94 showing amplification of the matter wave for certain negative refraction indices n.

Discussion

While negative refraction for electromagnetic waves was regarded a mathematical curiosity until the emergence of meta-materials, it now is an established and well-tested physical phenomenon. The results presented in this paper show the feasibility of negative refraction materials for matter waves due to the variety of free parameters of the many-body problem under consideration namely, the generally complex-valued scattering amplitude, the real-valued number density of local scatterers and the effective local field, which itself is a function of the former two and can be approximated by a plane wave when neglecting multiple scattering processes at each scatterer but considering only the leading order. The general theoretical approach presented here extends to statistical ensembles, but trivially includes the case of a definite choice of an ensemble configuration of scatterers. One striking feature of negative refraction is the possibility of perfect transmission through the medium with - in contrast to a perfect mirror the outgoing wave occupies the space ‘behind’ the negative refraction material, but with opposite momentum. These features imply the feasibility of static and statistical super-lenses for matter particles with promising applications such as cloaking of objects regarding influx of specific matter - the object becomes ‘untouchable’ - and perfect focusing of atom beams to a single point.

In addition to the general framework we have stated an explicit scenario in terms of specific scatterers with realistic numerical values in a formally simpler framework, which provides a first guide for developing meta-materials for negative matter scattering. The stated examples are chosen to be as simple as possible. In the first explicit example the consideration is confined to 1D and scattering of a single frequency is assumed while wave packets are merely integrals thereof. Furthermore we suppose that the scatterers are located around a single point. The effective field there corresponds to a complex function. Then for the given choice of the incoming wave the associated wave packet is negatively refracted. Extensions thereof show the flexibility of the presented scheme. In addition feasibility has been shown for a 3D scenario when considering a far field approximation. Limitations of the scheme are clearly given by that fact that only specific particle distributions ρ with a certain fixed frequency dependence of the scattering amplitude and scattering behaviour of their generating Green’s function G imply negative scattering. Further when considering statistical ensembles, instead of a definite choice, the considerations made only apply in the large N limit.

Methods

I utilized multiple scattering theory to analyse the feasibility of negative refraction for matter waves systems due to multiple scattering events. After introducing a general widely accepted and experimentally well-tested approach analytical expressions for the effective local field, Green’s functions, density and the frequency behaviour of the scatterers were stated, that consequently yield an outgoing counter-propagating/negatively refracted wave in analogy to the scenario observed for EM waves. The deductive reasoning/mathematical analysis is based on standard analytical tools, e.g. approximations in the large N limit, utilizing continuity properties of the wave function at the surfaces of the negative refractive material - an essential property in quantum mechanics. The latter implied the complete transmission of the quantum particle for particular . Finally the proposed negative refractive materials for sub-wavelength localisation of matter waves are an extension from light waves and are justified by analog behaviour of waves in materials with refraction index n independent of their physical nature.

Additional Information

How to cite this article: Pinsker, F. Multiple scattering induced negative refraction of matter waves. Sci. Rep. 6, 20751; doi: 10.1038/srep20751 (2016).