## 1 Introduction

The Schrödinger equation with varying mass has gained great attention in solid state physics, and been successfully used to to study electronic properties of semiconductor hetero-structures [27, 24, 11, 29]

. For example, it can describe localized defects in crystalline media, which may yield bound states localized to the defect. It is also related to describe non-compact line defect (edge) perturbations, e.g., by interpolating two-dimensional honeycomb structures via domain walls

[13, 12]. Edge modes may be produced by such perturbations, which propagate in the direction parallel to the edge and localize transverse to the edge.We are interested in deriving the asymptotic limit of the following varying-mass Schrödinger equation

(1.1) |

where , with and is the rescale Planck constant. The varying mass , which can be time-dependent, is assume to be random and highly oscillatory, with a given covariance matrix in time and space. We shall assume that decays fast enough at infinity to validate all the derivations. One goal of the paper is to show that in the regime, the Wigner transform of the solution converges to a special radiative transfer equation.

The problem is motivated by a fact that simulating (1.1) is extremely challenging in the semi-classical regime (). The challenges are two-folded. In deterministic regime (meaning is a deterministic highly oscillatory function in ), standard numerical solvers require to the small wavelength of to be resolved, for example, a mesh size and time step of order is required when finite difference methods are used [22, 23]. The time-splitting spectral method [4, 6] can improve the mesh size to be of order , however, has limitations to compute the Schrödinger equation with varying mass. A bigger problem comes from the randomness in . Since only the covariance of

is given, numerically one has to find many realizations and compute the deterministic Schrödinger equation before finding the ensemble mean/variance of the solution. The number of realizations, however, increases as

, as details in the random fluctuation become more and more important.Disregarding the challenges from the stochasticity, merely for the deterministic system, alternative approaches have been developed. These included the WKB-type methods, e.g., Gaussian beam methods [18, 17] and frozen Gaussian approximation [15, 19, 20]. The idea is to apply the WKB-type ansatz

and derive the eikonal equation for and transport-like equation for , where both and are functions of large scale, i.e., -independent. To our best knowledge, no such types of methods have been applied to efficiently solve (2.1) in the literature yet. And even for standard Schrödinger equation with random potential term, the application of the methods have not been fully understood.

In the paper, we shall systematically derive asymptotic equations for (2.1) by the Wigner transform [14], which is a main tool in semiclassical theory parallel to the WKB-type methods mentioned above. The literature on deriving the asymptotic equations for wave propagation in random media [2, 3, 1, 26, 21, 10, 9, 8, 7] is rich, most of which started with the Schrödinger equation with constant mass, and the randomness and high oscillations are introduced through the potential term. When it is the effective mass term that is random and highly oscillatory, the process of the derivation is rather similar but much more delicate, as will be detailed later in our paper. As a proof of concept, we numerically verify the derived radiative transfer equation by carefully computing and comparing its solution to the one of VMSE (1.1), and show that the two solutions agree.

The rest of the paper is organized as follows. To better illustrate the derivation, we start with a simpler case with being deterministic and only spatially dependent, and derive the limiting radiative transfer equation by the Wigner transform in Section 2. In Section 3, we systematically introduce the derivation of the limiting equation for the varying-mass Schrödinger equation (1.1) with random heterogeneities. We present our numerical validation in Section 4 and make conclusive remarks in Section 5.

## 2 Wigner transform of VMSE in the deterministic setting

As a preparation, we first investigate the semi-classical limit for (1.1) with deterministic mass in this section, which is to consider:

(2.1) |

The varying mass is a real function of . It is assumed to be deterministic and time-independent. We study the Cauchy problem, i.e., the is assumed to decay at infinite. is a complex function and we typically care about its physical observables such as the energy density and the energy flux , defined respectively as

Wigner transform is a technique explored in [25] for the Schrödinger equation with random potential, and has been demonstrated as a very powerful tool for investigating the semi-classical limit. It defines a function on the phase space:

(2.2) |

Here is the complex conjugate of

. This definition is essentially the Fourier transform of

on the variable.

The Wigner transform loses phase information: if is perturbed to

, the Wigner transform is kept the same. However, the physical observables can be recovered, namely, the first and second moments of

provide the energy density and the energy flux:(2.3) |

It is not guaranteed that is positive, and thus it does not serve directly as the energy density on the phase space. By plugging in the Schrödinger equation, we derives the equation satisfied by in the following Lemma.

Let satisfy the VMSE (2.1), then its Wigner transform (2.2) satisfies:

(2.4) |

where the three operators are defined to be:

(2.5) |

(2.6) |

(2.7) | ||||

where the Fourier transform of in space is defined by

(2.8) |

Since the two terms and are conjugate with for the second term, we only study the first one. With integration by parts:

(2.11) | ||||

We treat the and respectively in the following. Perform integration by parts again to

(2.12) | ||||

Note that and the last term can be combined so that a complete gradient of is available, namely one arrives at a formula for

(2.13) | ||||

For in (2.11), integration by parts against produces

(2.14) | ||||

On the other hand, integration by parts against gives

(2.15) | ||||

Note that and can be combined after another integration by parts

(2.16) | ||||

and can be combined similarly

(2.17) | ||||

Hence using (2.14)-(2.17) and the trick in (2.13), one derives the formula for in (2.11)

(2.18) | ||||

Finally from (2.13) and (2.18), one gets

(2.19) | ||||

All the terms can be explicitly expressed by the Wigner transform (2.2). In particular:

(2.20) | ||||

We use as an example to show this. Recalling:

we have

(2.21) | ||||

where we used the fact that

(2.22) |

Using (2.20), we get

(2.23) | ||||

By the conjugate argument, one gets, setting :

(2.24) | ||||

Finally, substitute (2.23) and (2.24) into (2.10), and we arrive at the Wigner equation in (2.4).

[Formal] With certain regularity, the formal limit of (2.4) is a transport equation:

(2.25) |

The trajectory of particles follows:

(2.26) |

and

(2.27) |

To prove the limiting equation (2.25), one only needs to derive the limiting behavior of . Indeed, formally, assuming :

and that

One then arrives at (2.25) by plugging them in the original equation (2.4). To show (2.27) with , one simply takes the time derivative along the trajectory:

according to the trajectory equation (2.26), and thus the quantity is conserved along the trajectory.

This proposition essentially guarantees the positivity of the solution to the limiting equation (2.25).

## 3 Semi-classical limit for VMSE with random perturbation

We consider the VMSE where the effective mass involves random perturbation, namely:

(3.1) |

where the effective mass is

(3.2) |

While the leading order is assumed to be deterministic and smooth, we allow the random perturbation to present small scales at . Furthermore we assume it is mean-zero and stationary in both and with the correlation function :

(3.3) |

Taking the Fourier transform of the function in both time and space, one has:

(3.4) |

then it is straightforward to show:

(3.5) |

and

We dedicate this section to the derivation of the semi-classical limit of the equation above. We will show that In the zero limit of , the Wigner transform of , the solution to the VMSE (2.1) with varying random mass (3.2), solves the radiative transfer equation:

(3.6) |

In view of (2.4) in Lemma 2, noting that , the Wigner equation (3.1) is transformed to

(3.7) |

where the operators are defined in (2.5)-(2.7), and are their counterparts defined by :

(3.8) |

(3.9) |

and

Comments

There are no comments yet.