Welcome to visit Communications in Theoretical Physics,
Quantum Physics and Quantum Information

Quantum carpets: efficiently probing fractional revivals in position-dependent mass systems

  • Tooba Bibi ,
  • Sunia Javed ,
  • Shahid Iqbal , *
Expand
  • Department of Physics, School of Natural Sciences, National University of Sciences and Technology, H-12 Islamabad, Pakistan

*Author to whom all correspondence should be addressed.

Received date: 2022-09-02

  Revised date: 2022-10-22

  Accepted date: 2022-11-01

  Online published: 2022-12-22

Copyright

© 2022 Institute of Theoretical Physics CAS, Chinese Physical Society and IOP Publishing

Abstract

Position-dependent-mass systems are of great importance in many physical situations, such as the transport of charge carriers in semiconductors with non-uniform composition and in the theory of many-body interactions in condensed matter. Here we investigate, numerically and analytically, the phenomenon of fractional revivals in such systems, which is a generic characteristic manifested by the wave-packet evolution in bounded Hamiltonian systems. Identifying the fractional revivals using specific probes is an important task in the theory of quantum measurement and sensing. We numerically simulate the temporal evolution of probability density and information entropy density, which manifest self-similarly recurring interference patterns, namely, quantum carpets. Our numerical results show that the quantum carpets not only serve as an effective probe for recognizing the fractional revivals of various order but they efficiently describe the effect of spatially-varying mass on the structure of fractional revivals, which is manifested as a symmetry breaking in their designs.

Cite this article

Tooba Bibi , Sunia Javed , Shahid Iqbal . Quantum carpets: efficiently probing fractional revivals in position-dependent mass systems[J]. Communications in Theoretical Physics, 2023 , 75(1) : 015102 . DOI: 10.1088/1572-9494/ac9f0d

1. Introduction

The idea of position-dependent mass (PDM) can emerge in many physical situations, such as the effective theory of many-body interactions in the field of condensed-matter physics [14]. Typically, a charge carrier in a crystal may carry effective mass, meff, significantly different from its corresponding free-state mass m0. Moreover, in many practical situations, meff may become anisotropic, which leads to spatially varying mass. Because of advancements in nanotechnology to fabricate ultra-smart semiconductor devices, the theory of PDM has become significantly important to study various physical properties of semiconductor structures [516]. However, the quantum theory of PDM systems is challenging due to the non-uniqueness of the quantum Hamiltonian within the same many-body approximation and because of the mathematical complexity in finding the analytic solutions. Conventionally, the solutions of PDM systems are obtained by directly solving the underlying Schrödinger equation [17, 18], however, various other procedures are also available in the literature, for instance, potential algebras using the notion of shape invariance [1925], the path integral formalism [26] and by employing the point canonical transformations [2730]. Nonetheless, the possibility of finding analytical solutions of PDM systems has triggered numerous other important explorations, such as the construction of generalized coherent states [3135] and the calculation of information entropy [3642], which may be beneficial in quantum technologies.
In bounded systems, wave-packet dynamics manifest the phenomena of quantum revivals and fractional revivals as their inherent features, which have been well studied in various physical contexts [4353] and [5457], respectively. In particular, the phenomenon of fractional revivals has enormous significance in the theory of quantum measurement and quantum sensing, since they can retrieve information about the particle even if the corresponding wave packet collapses well before its first full quantum revival time [54]. Moreover, the notion of fractional revivals can be useful in mapping the quantum phase of molecular wave-packet in the context of spectroscopy [55], in developing the scanning microscope [56] and even to factorize the prime numbers [57]. However, identifying the wave-packet fractional revivals with better resolution using various physically observable quantities as measurement probes [46, 51] is still an important task. A conventional approach to studying fractional revivals is to analyze the autocorrelation function [46], which measures the overlap of the time-evolved wave packet onto the initial wave packet. Nonetheless, some other measures, such as the time-evolution of probability density [58] and information entropy [51, 59] have also been introduced in this context. In this letter, we are aimed to explore the phenomenon of wave-packet fractional revivals in PDM systems, which has not been much discussed, except in some recent studies [30, 53, 59, 60]. However, it has been shown recently [59, 60], that the effect of PDM on the structure of fractional revivals cannot be fully captured by usual probes, for instance, by the autocorrelation function. However, in this context, it has been shown [59] that the information entropy does have a better sensitivity to elucidate the effect of spatially-varying mass.
In this work, we investigate the wave-packet fractional revivals in PDM systems by numerically simulating the temporal evolution of probability and Shannon information entropy densities, both in the position as well as momentum space, which manifests the recurrence of self-similar interference patterns, termed as the quantum carpets. The fractional revivals of the wave packet are then identified by the recurrence order of self-similar structures in the design of such quantum carpets. The concept of quantum carpets has been used in numerous contexts, such as, to explore wave-packet quantum recurrences [58], for studying decoherence effects [61], to explain the Talbot effect [62] and to investigate various probability density profiles [6365]. Here, in addition to conventional position-time (spatio-temporal) quantum carpets, we extend this notion to momentum-time (momento-temporal) and entropy-time (information) quantum carpets. It has been shown by our numerical results that the quantum carpets not only serve as an effective probe for identifying the fractional revivals of various orders but they efficiently elucidate the effect of spatially-varying mass on the underlying phenomenon, which is manifested as symmetry breaking in their designs. Moreover, it has been observed that, under the influence of the PDM, the fractional revivals are shifted in position-space towards the region where the particle experiences a higher mass.
The rest of the paper is organized as follows. In section 2 we present the quantization procedure for PDM systems and then present the analytic solutions of the Schrödinger equation corresponding to our model PDM system. Then in section 3 we discuss the analytic formalism regarding probability and information entropy densities required to numerically simulate the respective quantum carpets. We present, in section 4 the numerical results and discussions regarding fractional revivals of the PDM wave packet. Finally, in section 5 we present the conclusion of our work.

2. Schrödinger equation with PDM

While quantizing the PDM Hamiltonian, H = p2/2m(x) + V(x), the operators concerning p and m(x) in the kinetic energy term encounter an ambiguity with respect to their mutual ordering. In this context, a general procedure was proposed by von Roos [10, 11], and the corresponding kinetic energy operator is given as
$\begin{eqnarray*}\hat{T}=\displaystyle \frac{-1}{4}\left[{m}^{\alpha }(x)\displaystyle \frac{{\rm{d}}{m}^{\beta }(x)}{{\rm{d}}x}\displaystyle \frac{{\rm{d}}{m}^{\gamma }(x)}{{\rm{d}}x}+{m}^{\gamma }(x)\displaystyle \frac{{\rm{d}}{m}^{\beta }(x)}{{\rm{d}}x}\displaystyle \frac{{\rm{d}}{m}^{\alpha }(x)}{{\rm{d}}x}\right],\end{eqnarray*}$
where α, β, γ are so-called ambiguity parameters satisfying the constraint α + β + β = −1. It is to be noted that by choosing various sets of values for α, β and γ, one finds different non-equivalent quantum Hamiltonians [12]. Nonetheless, the special choice, α = γ = 0 and β = −1, results in a symmetric ordering of the underlying operators [12]. In this case, the corresponding Schrödinger equation becomes [17, 18] as
$\begin{eqnarray}\displaystyle \frac{-1}{2}\displaystyle \frac{{\rm{d}}}{{\rm{d}}x}\left[\displaystyle \frac{1}{m(x)}\displaystyle \frac{{\rm{d}}\psi (x)}{{\rm{d}}x}\right]+V(x)\psi (x)=E\psi (x).\end{eqnarray}$
The particular solutions of this equation depend on the choice of the profile of PDM m(x) and on the type of confining potential V(x). Various authors have attempted to obtain the exact solutions of this equation for various choices of m(x) and V(x) [17, 18]. Here we consider a particle with spatially-varying mass
$\begin{eqnarray}m(x)=\displaystyle \frac{{m}_{0}}{{\left(1+\tau x\right)}^{2}},\end{eqnarray}$
trapped inside a one-dimensional infinitely deep square well of length L, defined as $V(\hat{x})=0$ for 0 < x < L and $V(\hat{x})=\infty $ otherwise. In equation (2) above, the dependence of the particle's mass on position is governed by a strength parameter, τ > 0, such that, m(x) → m0 when τ → 0, where m0 is the corresponding constant free-state mass of the particle, which will be taken as m0 = 1 in our later calculations. Upon following parametric transformation
$\begin{eqnarray}\begin{array}{rcl}\psi (x) & = & m{\left(x\right)}^{\displaystyle \frac{1}{4}}\varphi (z),\\ \displaystyle \frac{{\rm{d}}z}{{\rm{d}}x} & = & m{\left(x\right)}^{\displaystyle \frac{1}{2}}=\displaystyle \frac{1}{(1+\tau x)},\end{array}\end{eqnarray}$
the Schrödinger equation, given in (1), reduces into a simpler differential equation as
$\begin{eqnarray}-\displaystyle \frac{1}{2}\displaystyle \frac{{{\rm{d}}}^{2}\varphi (z)}{{\rm{d}}{z}^{2}}+\displaystyle \frac{{\tau }^{2}\varphi (z)}{8}=E\varphi (z),\end{eqnarray}$
where, $z=\mathrm{ln}(1+\tau x)$, as is evident from equation (3). It is straightforward to find the solutions to equation (4), by applying the boundary conditions, ψ(0) = ψ(L) = 0. As a result, we find the normalized eigenfunctions as
$\begin{eqnarray}{\psi }_{n}^{(\tau )}(x)=\sqrt{\displaystyle \frac{2\tau }{(1+\tau x)\mathrm{ln}(1+\tau L)}}\sin \left(\displaystyle \frac{n\pi \mathrm{ln}(1+\tau x)}{\mathrm{ln}(1+\tau L)}\right),\end{eqnarray}$
and the corresponding eigenenergies as
$\begin{eqnarray}{E}_{n}^{(\tau )}=\displaystyle \frac{{\tau }^{2}}{8}+\displaystyle \frac{{n}^{2}{\pi }^{2}{\tau }^{2}}{2{\mathrm{ln}}^{2}(1+\tau L)},\,\,\,n=0,1,2,...,\end{eqnarray}$
where the ground-state energy takes the value as E0 = τ2/8. However, it is to be noted from equation (5) that the eigenfunctions exist only for ${E}_{n}^{(\tau )}\gt {\tau }^{2}/8$ (i.e. ground-state eigenfunction, ${\psi }_{0}^{(\tau )}(x)$, vanishes and constitutes the so-called zero-curvature solution [66]).
The eigenfunctions and eigenenergies, given in equations (5) and (6) respectively, enable us to construct the wave-packet for the PDM system and to study the wave-packet time-evolution through probability and information entropy densities both analytically and numerically.

3. Quantum carpets woven by PDM: identifying fractional revivals

As mentioned earlier, the quantum carpets are woven by self-interference of the probability density of a wave packet during its temporal evolution in a given bounded system. Consequently, the self-similarly recurring interference pattern can serve as a probe to identify the structure of fractional revivals. Mostly, position-time (spatio-temporal) quantum carpets have been discussed in the literature [58, 61, 6365], which are woven by the temporal evolution of position-space probability density. In addition, analogously, we also present the momentum-time and the entropy-time quantum carpets in the following subsections. However, the solutions of PDM Schrödinger equation (1), obtained in the previous section 2 provide us the starting point to obtain the analytical expressions for probability and entropy densities in position and momentum spaces for our chosen PDM system.

3.1. Position-time quantum carpets

The time-dependent wave packet for the given PDM system can be expressed as a linear superposition of the corresponding eigenfunctions
$\begin{eqnarray}{{\rm{\Psi }}}^{(\tau )}(x,t)=\sum _{n=0}^{\infty }{c}_{n}{\psi }_{n}^{(\tau )}(x){{\rm{e}}}^{-{\rm{i}}{E}_{n}^{(\tau )}t/{\hslash }},\end{eqnarray}$
where ${E}_{n}^{(\tau )}$ and ${\psi }_{n}^{(\tau )}(x)$ are obtained by the solutions of PDM Schrödinger equation (1), given in equations (6) and (5), respectively. Here, the expansion coefficient, cn, can be obtained as ${c}_{n}={\int }_{0}^{L}{{\rm{\Psi }}}^{(\tau )}(x,t=0){\psi }_{n}^{* (\tau )}(x){\rm{d}}x$, where $\Psi$(τ)(x, t = 0) is the initial profile of the wave packet in position space. For the sake of convenience, we shall consider an initial wave packet with a Gaussian profile, centered at x0 and with initial momentum p0, which is defined as
$\begin{eqnarray}{{\rm{\Psi }}}^{(\tau )}(x,0)=\displaystyle \frac{1}{{\left({\sigma }^{2}\pi \right)}^{\tfrac{1}{4}}}\exp \left(\displaystyle \frac{-{\left(x-{x}_{0}\right)}^{2}}{2{\sigma }^{2}}\right)\exp \left(\displaystyle \frac{{\rm{i}}{p}_{0}x}{{\hslash }}\right),\end{eqnarray}$
where σ is the spread of the wave packet at FWHM. Using equation (7), the time-dependent probability density in position space is calculated as
$\begin{eqnarray}\begin{array}{l}{P}^{(\tau )}(x,t)=| {{\rm{\Psi }}}^{(\tau )}(x,t){| }^{2}\\ \quad =\displaystyle \sum _{n,m=1}^{\infty }{c}_{n}{c}_{m}^{* }{\psi }_{n}^{(\tau )}(x){\psi }_{m}^{* (\tau )}(x){{\rm{e}}}^{-{\rm{i}}({E}_{n}^{(\tau )}-{E}_{m}^{(\tau )})t/{\hslash }}.\end{array}\end{eqnarray}$
During its temporal dynamics, the wave packet undergoes a series of constructive and destructive interference, which leads to manifesting revivals at various time scales. The natural time scales of these recurrences are defined [46] as
$\begin{eqnarray}\begin{array}{rcl}{T}_{j}^{(\tau )} & = & \frac{2\pi }{{\omega }_{j}^{(\tau )}},\ \ \ \ \ \ \mathrm{where}\\ {\omega }_{j}^{(\tau )} & = & {\left({\left.\frac{1}{j!}\left|\frac{{{\rm{d}}}^{j}{E}_{n}^{(\tau )}}{{\rm{d}}{n}^{j}}\right|\right|}_{n={n}_{0}}\right)}^{-1},\end{array}\end{eqnarray}$
such that, ${T}_{1}^{(\tau )}={T}_{\mathrm{cl}}^{(\tau )}$, is the classical period and ${T}_{2}^{(\tau )}={T}_{\mathrm{rev}}^{(\tau )}$ is the quantum revival time, which follow the hierarchy as ${T}_{\mathrm{cl}}^{(\tau )}\lt {T}_{\mathrm{rev}}^{(\tau )}\lt ...$. In the present case, the quantum revival time, ${T}_{\mathrm{rev}}^{(\tau )}$, appears to be the highest time scale of the PDM system.
Using equation (6), the classical period and quantum revival time, for PDM particles trapped in an infinite square well, take the values as
$\begin{eqnarray}{T}_{\mathrm{cl}}^{(\tau )}=\displaystyle \frac{2{\mathrm{ln}}^{2}(1+\tau L)}{{n}_{0}\pi {\tau }^{2}},\,\,\,\,\,\,\,{T}_{\mathrm{rev}}^{(\tau )}=\displaystyle \frac{4{\mathrm{ln}}^{2}(1+\tau L)}{\pi {\tau }^{2}},\end{eqnarray}$
respectively. To recognize the effect of PDM on the structure of quantum recurrence, let us consider the corresponding constant mass system having energy ${E}_{n}^{(0)}=\tfrac{{n}^{2}{\pi }^{2}{{\hslash }}^{2}}{2{m}_{0}{L}^{2}}$. In this case, the ${T}_{\mathrm{cl}}^{(0)}$ and the ${T}_{\mathrm{rev}}^{(0)}$ take the values as
$\begin{eqnarray}{T}_{\mathrm{cl}}^{(0)}=\displaystyle \frac{2{L}^{2}}{{n}_{0}\pi },\qquad \qquad {T}_{\mathrm{rev}}^{(0)}=\displaystyle \frac{4{L}^{2}}{\pi }.\end{eqnarray}$
It can be seen from equations (11) and (12), that the ratio, ${T}_{\mathrm{rev}}^{(0)}/{T}_{\mathrm{cl}}^{(0)}=2{n}_{0}={T}_{\mathrm{rev}}^{(\tau )}/{T}_{\mathrm{cl}}^{(\tau )}$, for both constant mass (CM) and PDM systems are the same. This fact implies that the effect of PDM on classical periodicity and quantum revival time is nothing more than a constant scaling of ${T}_{\mathrm{cl}}^{(\tau )}$ and ${T}_{\mathrm{rev}}^{(\tau )}$ with same magnitude.
However, in addition to full revivals, the wave-packet revivals occur partially at times $t=(r/s){T}_{\mathrm{rev}}^{(\tau )}$, where r and s are relative prime numbers. The numerical simulation of probability density, P(τ)(x, t), given in equation (9), as a function of time t and position x leads to generating the position-time quantum carpets, which exhibit the structure of fractional revivals as displayed in figures 14.
Figure 1. The position-space probability density, P(τ)(x, t), as defined in equation (9), is plotted for different values of PDM strength parameter τ, after evolution time: (a) $t={T}_{\mathrm{rev}}^{(\tau )}/4;$ (b) $t={T}_{\mathrm{rev}}^{(\tau )}/2$. These plots show the effect of spatial dependence (defined by strength parameter τ) on a given fractional revival. Initially, the Gaussian wave packet has width σ = 0.1 which is centered at the middle of the well x0 = L/2.

3.2. Momentum-time quantum carpets

In order to obtain an analogous description of quantum carpets in momentum space, we first need to transform the eigenstates of the system into momentum space. This can be done using the Fourier transform of ${\psi }_{n}^{(\tau )}(x)$ as
$\begin{eqnarray}{\phi }_{n}^{(\tau )}(p)=\displaystyle \frac{1}{\sqrt{2\pi {\hslash }}}{\int }_{-\infty }^{\infty }{\psi }_{n}^{(\tau )}(x){{\rm{e}}}^{-\tfrac{{\rm{i}}{px}}{{\hslash }}}{\rm{d}}x.\end{eqnarray}$
It is worth noting that the Fourier Integral (13) vanishes outside the interval $\left[0,L\right]$ because of the boundary conditions chosen to solve the PDM Schrödinger equation in equation (4).
Using equation (13), the time-dependent wave packet in momentum space is given as
$\begin{eqnarray}{{\rm{\Phi }}}^{(\tau )}(p,t)=\sum _{n=0}^{\infty }{c}_{n}{\phi }_{n}^{(\tau )}(p){{\rm{e}}}^{-{\rm{i}}{E}_{n}^{(\tau )}t/{\hslash }},\end{eqnarray}$
and the corresponding time-dependent probability density in momentum space is now given as
$\begin{eqnarray}\begin{array}{l}{Q}^{(\tau )}(p,t)=| {{\rm{\Phi }}}^{(\tau )}(p,t){| }^{2}\\ \quad =\displaystyle \sum _{n,m=0}^{\infty }{c}_{n}{c}_{m}^{* }{\phi }_{n}^{(\tau )}(p){\phi }_{m}^{* (\tau )}(p){{\rm{e}}}^{-{\rm{i}}({E}_{n}^{(\tau )}-{E}_{m}^{(\tau )})t/{\hslash }}.\end{array}\end{eqnarray}$
For the momentum-time (momento-temporal) quantum carpets, we numerically compute, Q(τ)(p, t), as a function of p and t and the corresponding results are displayed in figures 5 and 6.

3.3. Entropy-time quantum carpets: quantum carpets of information

Having defined the probability density of the wave packet in position and momentum space, as given in equations (9) and (15), it is instructive to calculate the corresponding information entropies. For instance, the position-space Shannon information entropy is defined [59, 67] as
$\begin{eqnarray}{S}_{x}^{(\tau )}(t)=-\int {P}^{(\tau )}(x,t)\mathrm{ln}{P}^{(\tau )}(x,t){\rm{d}}x,\end{eqnarray}$
where P(τ)(x, t) is given above in equation (9). Analogously, using equation (15), the corresponding momentum space information entropy is defined as
$\begin{eqnarray}{S}_{p}^{(\tau )}(t)=-\int {Q}^{(\tau )}(p,t)\mathrm{ln}{Q}^{(\tau )}(p,t){\rm{d}}p.\end{eqnarray}$
It is worth mentioning that the information entropy ${S}_{\zeta }^{(\tau )}$ quantifies the amount of uncertainty associated with the localization of the wave packet in a given space ζ (where ζ denotes either x or p space). For instance, the smaller is the information entropy, and sharply localized is the wave packet. This further implies that the lower the uncertainty, the greater the accuracy) in predicting the particle's localization. Hence, the temporal evolution of the information entropy, ${S}_{\zeta }^{(\tau )}(t)$, may have the potential to encode the information of the wave-packet fractional revivals. Moreover, it is found that the sum, ${{\rm{\Sigma }}}^{(\tau )}(t)={S}_{x}^{(\tau )}(t)+{S}_{p}^{(\tau )}(t)$, satisfies the relation [59, 67]
$\begin{eqnarray}{{\rm{\Sigma }}}^{(\tau )}(t)\geqslant 1+\mathrm{ln}\pi ,\end{eqnarray}$
which is an entropy-based version of the Heisenberg uncertainty relation, analogous to the usual variance-based one ΔxΔp/2. However, it is found that, in addition to ${S}_{x}^{(\tau )}(t)$ and ${S}_{p}^{(\tau )}(t)$, the sum Σ(τ)(t) can also be used to identify the structure of fractional revivals. In our recent work [59], we have studied the fractional revivals for PDM tapped in an infinite square well [59] using information theoretic measures defined above in equations (16), (17) and (18).
In the present context of quantum carpets, let us define the entropy density as a function of position and time,
$\begin{eqnarray}{{\rm{\Gamma }}}^{(\tau )}(x,t)=-{P}^{(\tau )}(x,t)\mathrm{ln}{P}^{(\tau )}(x,t),\end{eqnarray}$
which is essentially the integrand of the integral (16) defining the position-space information entropy. Analogously, using equation (17), the entropy density in momentum space can be defined as
$\begin{eqnarray}{{\rm{\Lambda }}}^{(\tau )}(p,t)=-{Q}^{(\tau )}(p,t)\mathrm{ln}{Q}^{(\tau )}(p,t).\end{eqnarray}$
It is important to note that these entropy densities may have the sensitivity to detect the local effect of PDM on the structure of fractional revivals. By numerically simulating information entropy densities Γ(τ)(x, t) and Λ(τ)(p, t), given above in equations (19) and (20), we can generate the information entropy-time quantum carpets which we referred to as information quantum carpets or simply, information carpets. The snapshots of these entropy densities, exhibiting the fractional revivals, are displayed in figure 7.

4. Numerical results and discussions

In order to characterize the phenomenon of fractional revivals using quantum carpets as probes for their identification, we numerically compute the probability and entropy densities defined above in equations (9), (15), (19) and (20). As mentioned earlier in previous section 3, for our numerical simulations, we consider an initial Gaussian wave packet, defined in equation (8), having width σ = 0.1 and centered at the middle of the well x0 = L/2. Moreover, for convenience of our numerical computation, we choose the arbitrary units for which = m0 = L = 1. To better understand the structure of quantum carpets in a given space, we first take the snapshots of corresponding probability and information entropy density after different evolution times, specifically, chosen at fractional revival time of various orders, defined as $t=(r/s){T}_{\mathrm{rev}}^{(\tau )}$ with r and s are prime numbers.
For instance, we numerically compute the position-space probability density, P(τ)(x, t), defined by equation (9), as a function of position x after different fractional revivals and for various values of PDM strength parameter τ. The corresponding results are displayed in figures 1 and 2. In particular, the plots in figure 1 show that at any chosen fractional revival time, the probability density peaks shift toward the left half of the well where PDM experiences a higher mass. Whereas, figure 2 shows a comparison of the main fractional revivals of the PDM particle with the ones of the corresponding constant mass (CM) particle. It can be seen from these plots that the symmetric probability density peaks, in the case of fractional revivals of CM particle (τ = 0), get shifted towards the left region when the spatial dependence of the mass is turned on (τ > 0).
Figure 2. The snapshots of position-space probability density, P(τ)(x, t), as defined in equation (9), are plotted after different evolution time: (a) $t={T}_{\mathrm{rev}}^{(\tau )}/8;$ (b) $t={T}_{\mathrm{rev}}^{(\tau )}/4;$ (c) $t={T}_{\mathrm{rev}}^{(\tau )}/2;$ (d) $t=3{T}_{\mathrm{rev}}^{(\tau )}/4;$ (e) $t=7{T}_{\mathrm{rev}}^{(\tau )}/8;$ (f) $t={T}_{\mathrm{rev}}^{(\tau )}$. These plots show a comparison of fractional revivals of the PDM system (τ > 0) with those of the corresponding constant mass system (τ = 0). The values of characteristic parameters of initial Gaussian are chosen the same as in figure 1.
Figure 3 shows the generation of position-time quantum carpet, where the probability density function, P(τ)(x, t), has been plotted as a function of x and t for short-term temporal evolution (${\rm{\Delta }}t=4{T}_{\mathrm{cl}}^{(\tau )}$). The strength of probability density has been represented by means of colors, such that, brightness (darkness) indicates the high (low) probability density, as shown by the color scaling in figure 3. It can be seen from these plots that an initially well-localized (Gaussian) wave packet, centered at the middle of the box, starts its time evolution tending to follow a classical trajectory between two boundaries. However, after a few classical periods dephasing among various constituents eigenstates of the wave-packet dominates and as a result wave packet experiences a collapse. Later on, fractional or full revivals of the wave packet occur when the phase-matching condition is regained, respectively, either partially or completely. Moreover, figure 4 shows a comparison of the position-time quantum carpet of the CM particle with that of the PDM particle. It is clearly seen from these plots that the symmetry of quantum carpets gets broken, such that, the probability of finding the particle increases in the region where it experiences a higher mass.
Figure 3. The short-term temporal evolution (${\rm{\Delta }}t=4{T}_{\mathrm{cl}}^{(\tau )}$) of position-space probability density P(τ)(x, t), defined in equation (1), indicates the formation of position-time quantum carpet: (a) constant mass (CM) particle; (b) PDM particle. The initial profile of the wave packet is the same as in figure 1.
Figure 4. The position-time quantum carpets woven by position-space probability density, P(τ)(x, t), of the wave-packet during its long-time evolution (${\rm{\Delta }}t={T}_{\mathrm{rev}}^{(\tau )}/2$): (a) CM particle; (b) PDM particle. The initial profile of the wave packet is the same as in figure 1.
Similarly, to visualize the analogous description of fractional revivals in momentum-space, we numerically compute the momentum-space probability density, Q(τ)(p, t), as function p and t, for different values of strength parameter τ. The results regarding the snapshots of Q(τ)(p, t) at different fractional revivals have been displayed in figure 5 whereas, the momentum-time quantum carpet for an evolution time interval ${\rm{\Delta }}t={T}_{{\rm{rev}}}^{(\tau )}/2$ is shown in figure 6. From these plots, it is evident that probability peaks become asymmetric as soon as the mass of the particle becomes position-dependent (τ > 0). Likewise, we numerically simulate information entropy densities in position and momentum spaces, given in equations (19) and (20), and the corresponding results are displayed in figure 7. A detailed analysis of fractional revivals using information entropy densities has recently been presented [59].
Figure 5. The momentum-space probability density, Q(τ)(p, t), as defined in equation (15), is plotted for different values of PDM strength parameter τ, after evolution time: (a) t = Trev/8; (b) t = Trev/4. These plots are aimed to show the effect of PDM on fractional revivals in momentum space. The initial profile of the wave packet is kept the same as in figure 1.
Figure 6. The momentum-time quantum carpets woven by momentum-space probability density, Q(τ)(p, t), during its long-time evolution ${\rm{\Delta }}t={T}_{\mathrm{rev}}^{(\tau )}/2$: (a) CM particle; (b) PDM particle. The initial profile of wave-packet is the same as in figure 1.
Figure 7. The information entropy densities, defined in equations (19) and (20), are plotted for different values of PDM strength parameter τ: (a) position-space entropy density after evolution time t = 3Trev/4; (b) momentum-space entropy density after evolution time t = Trev/4. The initial profile of the wave packet is kept the same as in figure 1.
From the above analysis, it can be seen from the position-time quantum carpets that under the influence of PDM, the fractional revivals get shifted in position-space towards the region where the particle experiences a higher mass. However, in any given space (whether position or momentum), the symmetry of the quantum carpets' designs is broken when the position-dependence of the mass is turned on (τ > 0). It is therefore concluded that quantum carpets, in any given space, not only serve as an effective probe for recognizing the fractional revivals but they efficiently elucidate the effect of spatially-varying mass on the underlying phenomenon, as indicated by the symmetry breaking in their designs.

5. Summary and conclusions

In the present work, we have investigated the phenomenon of fractional revivals in the varying-mass systems with respect to position. The PDM systems are of great significance in many physical situations of practical interest, such as the transport of charge carriers in semiconductors of non-uniform chemical composition and in the theory of many-body interactions in condensed matter physics. However, quantization of such systems is a challenging task because of the non-uniqueness of quantum Hamiltonian within the same many-body approximation and due to the mathematical complexity in finding the analytic solutions. In our work, we first presented a general procedure for the quantization of a PDM Hamiltonian by symmetrically ordering the operators concerning momentum and PDM. For a specific choice of PDM, we then solved the PDM Schrödinger equation for a special case of a PDM particle trapped in a one-dimensional infinitely deep square well. Such toy models have been very useful in many practical situations, for instance, in the fabrication of semiconductor structures, and are important to understand numerous fundamental notions of quantum mechanics [68].
Then we investigated the characteristic features of wave-packet dynamics in our chosen PDM system. Generally, detecting the fractional revivals in a given system is an important task in the theory of quantum measurement and sensing, which requires the use of some physically observable quantity such as a probe to analyze the underlying phenomena. In our work, we have computed, analytically and numerically, the temporal evolution of probability and information entropy densities in position and momentum spaces which lead to the formation of self-similarly recurring interference patterns known as quantum carpets. The fractional revivals of the wave packet have been identified by the recurrence order of self-similar patterns in the design of these quantum carpets. For our analysis, in addition to conventional position-time (spatio-temporal), we have simulated the momentum-time and entropy-time quantum carpets, which we would like to be named information carpets. It has been shown by our numerical results that the quantum carpets not only serve as an effective probe to identify the fractional revivals but they efficiently describe the effect of spatially-varying mass on the structure of fractional revivals, which is manifested as symmetry breaking in their designs. Moreover, it has been observed that, under the influence of spatially-varying mass, the fractional revivals are shifted in position-space towards the region where the particle experiences a higher mass.

Financial support from Higher Education Commission (HEC) of Pakistan, under Grant No. 20-14808/NRPU/R&D/HEC/2021 2021, is gratefully acknowledged.

1
Bastard G 1988 Wave Mechanics Applied to Semiconductor Heterostructures Les Editions de Physique Paris Les Ulis Cedex

2
Bastard G 1981 Superlattice band structure in the envelope-function approximation Phys. Rev. B 24 5693

DOI

3
Smith D L Mailhiot C 1990 Theory of semiconductor superlattice electronic structure Rev. Mod. Phys. 62 173

DOI

4
Trzeciakowski W Gurioli M 1991 Electric-field effects in semiconductor quantum wells Phys. Rev. B 44 3880

DOI

5
de Saavedra F A Boronat J Polls A Fabrocini A 1994 Effective mass of one 4He atom in liquid 3He Phys. Rev. B 50 4248

DOI

6
Geller M R Kohn W 1993 Quantum mechanics of electrons in crystals with graded composition Phys. Rev. Lett. 70 3103

DOI

7
Serra L Lipparini E 1997 Spin response of unpolarized quantum dots Europhys. Lett. 40 667

DOI

8
Peter A J Navaneethakrishnan K 2008 Effects of position-dependent effective mass and dielectric function of a hydrogenic donor in a quantum dot Physica E 40 2747

DOI

9
Keshavarz A Zamani N 2013 Optical properties of spherical quantum dot with position-dependent effective mass Superlattices Microstruct. 58 191

DOI

10
von Roos O 1983 Position-dependent effective masses in semiconductor theory Phys. Rev. B 27 7547

DOI

11
von Roos O Mavromatis H 1985 Position-dependent effective masses in semiconductor theory II Phys. Rev. B 31 2294

DOI

12
Lévy-Leblond J M 1995 Position-dependent effective mass and Galilean invariance Phys. Rev. A 52 1845

DOI

13
Morrow R A 1987 Establishment of an effective-mass Hamiltonian for abrupt heterojunctions Phys. Rev. B 35 8074

DOI

14
Yung K C Yee J H 1994 Derivation of the modified Schrödinger equation for a particle with a spatially varying mass through path integrals Phys. Rev. A 50 104

DOI

15
Einevoll G T 1990 Operator ordering in effective-mass theory for heterostructures: II. Strained systems Phys. Rev. B 42 3497

DOI

16
Koç R Koca M Şahinoğlu G 2005 Scattering in abrupt heterostructures using a position dependent mass Hamiltonian Eur. Phys. J. B 48 583

DOI

17
Amir N Iqbal S 2014 Exact solutions of Schrödinger equation for the position-dependent effective mass harmonic oscillator Commun. Theor. Phys. 62 790

DOI

18
Dong S H Lozada-Cassou M 2005 Exact solutions of the Schrödinger equation with the position-dependent mass for a hard-core potential Phys. Lett. A 337 313

DOI

19
Amir N Iqbal S 2015 Ladder operators and associated algebra for position-dependent effective mass systems Europhys. Lett. 111 20005

DOI

20
Amir N Iqbal S 2016 Algebraic solutions of shape-invariant position-dependent effective mass systems J. Math. Phys. 57 062105

DOI

21
Chaturvedi S Dutt R Gangopadhyaya A Panigrahi P Rasinariu C Sukhatme U 1998 Algebraic shape invariant models Phys. Lett. A 248 109

DOI

22
Balantekin A B 1998 Algebraic approach to shape invariance Phys. Rev. A 57 4188

DOI

23
Hounkonnou M N Sodoga K Azatassou E S 2005 Factorization of Sturm-Liouville operators: solvable potentials and underlying algebraic structure J. Phys. A 38 371

DOI

24
Dong S H Pena J J Pacheco-Garcia C Garcia-Ravelo J 2007 Algebraic approach to the position-dependent mass Schrödinger equation for a singular oscillator Mod. Phys. Lett. A 22 1039

DOI

25
Marques G D A Bezerra V B Dong S H 2013 Solution of the dirac equation with position-dependent mass in a coulomb and scalar fields in a conical spacetime Mod. Phys. Lett. A 28 1350137

DOI

26
Pak N K Sökmen I 1984 Exact path integral solution of a class of potentials related to the rigid rotator Phys. Lett. A 103 298

DOI

27
Rasinariu C Mallow J Gangopadhyaya A 2007 Exactly solvable problems of quantum mechanics and their spectrum generating algebras: A review Cent. Eur. J. Phys. 5 111

DOI

28
De R Dutt R Sukhatme U 1992 Mapping of shape invariant potentials under point canonical transformations J. Phys. A 25 L843

DOI

29
Yu J Dong S H 2004 Exactly solvable potentials for the Schrödinger equation with spatially dependent mass Phys. Lett. A 325 194

DOI

30
Schmidt A G M 2006 Wave-packet revival for the Schrödinger equation with position-dependent mass Phys. Lett. A 353 459

DOI

31
Amir N Iqbal S 2020 Coherent states of position-dependent mass trapped in an infinite square well J. Math. Phys. 61 082102

DOI

32
Amir N Iqbal S 2015 Coherent states for nonlinear harmonic oscillator and some of its properties J. Math. Phys. 56 062108

DOI

33
Amir N Iqbal S 2014 Comment on Coherent states for the nonlinear harmonic oscillator J. Math. Phys. 55 0114101

DOI

Amir N Iqbal S 2012 Comment on Coherent states for the nonlinear harmonic oscillator J. Math. Phys. 53 062104

DOI

34
Amir N Iqbal S 2016 Barut-girardello coherent states for nonlinear oscillator with position-dependent mass Commun. Theor. Phys. 66 41

DOI

35
Amir N Iqbal S 2016 Generalized coherent states for position-dependent effective mass systems Commun. Theor. Phys. 66 615

DOI

36
Falaye B J Serrano F A Dong S H 2016 Fisher information for the position-dependent mass Schrödinger system Phys. Lett. A 380 267

DOI

37
Yañez-Navarro G Sun G H Dytrych T Launey K D Dong S H Draayer J P 2014 Quantum information entropies for position-dependent mass Schrödinger problem Ann. Phys. 348 153

DOI

38
Sun G H Popov D Camacho-Nieto O Dong S H 2015 Shannon information entropies for position-dependent mass Schrödinger problem with a hyperbolic well Chin. Phys. B 24 100303

DOI

39
Sun G H Aoki M A Dong S H 2013 Quantum information entropies of the eigenstates for the Pöschl-Teller-like potential Chin. Phys. B 22 050302

DOI

40
Valencia-Torres R Sun G H Dong S H 2015 Quantum information entropy for a hyperbolical potential function Phys. Scr. 90 035205

DOI

41
Valencia-Torres R Sun G H Dong S H 2015 Quantum information entropy for a hyperbolical potential function Phys. Scr. 90 035205

DOI

42
Song X D Sun G H Dong S H 2015 Shannon information entropy for an infinite circular well Phys. Lett. A 379 1402

DOI

43
Alber G Ritsch H Zoller P 1986 Generation and detection of Rydberg wave packets by short laser pulses Phys. Rev. A 34 1058

DOI

44
Parker J Stroud C R Jr. 1986 Coherence and decay of Rydberg wave packets Phys. Rev. Lett. 56 716

DOI

45
Yeazell J A Mallalieu M Stroud C R Jr. 1990 Observation of the collapse and revival of a Rydberg electronic wave packet Phys. Rev. Lett. 64 2007

DOI

46
Robinett R W 2004 Quantum wave packet revivals Phys. Rep. 392 1

DOI

47
Buchleitner A Delande D Zakrzewski J 2002 Non-dispersive wave packets in periodically driven quantum systems Phys. Rep. 368 409

DOI

48
Saif F 2005 Classical and quantum chaos in atom optics Phys. Rep. 419 207

DOI

49
Iqbal S Qurat-ul-Ann Saif F 2006 Quantum recurrences in driven power-law potentials Phys. Lett. A 356 231

DOI

50
Riviére P Iqbal S Rost J M 2014 Electronic excitation by short x-ray pulses: from quantum beats to wave packet revivals J. Phys. B 47 124039

DOI

51
Romera E de los Santos F 2007 Identifying wave-packet fractional revivals by means of information entropy Phys. Rev. Lett. 99 263601

DOI

52
Iqbal S 2020 Dynamical interference of shannon information entropy: identification of wave-packet fractional revivals J. Russ. Laser Res. 41 451

DOI

53
Amir N Iqbal S 2017 Coherent states of nonlinear oscillators with position-dependent mass: temporal stability and fractional revivals Commun. Theor. Phys. 68 181

DOI

54
Ghosh S Banerji J 2007 A time-frequency analysis of wave packet fractional revivals J. Phys. B 40 3545

DOI

55
Schubert A Renziehausen K Engel V 2010 Mapping of quantum phases by two-dimensional vibronic spectroscopy of wave-packet revivals Phys. Rev. A 82 013419

DOI

56
Saif F 2006 Recurrence tracking microscope Phys. Rev. A 73 033618

DOI

57
Merkel W Averbukh S I Girard B Paulus G G Schleich W P 2006 Factorization of numbers with physical systems Fortschr. Phys. 54 856

DOI

58
Yousaf I Iqbal S 2016 Quantum carpets: a probe to identify wave-packet fractional revivals J. Russ. Laser Res. 37 328

DOI

59
Iqbal S 2022 Information entropy, fractional revivals and Schrödinger equation with position-dependent mass Phys. Scr. 97 015001

DOI

60
Iqbal S 2022 Quantum dynamical recurrences in position-dependent mass systems J. Russ. Laser Res. 43 96

DOI

61
Kazemi P Chaturvedi S Marzoli I O'Connell R F Schleich W P 2013 Quantum carpets: a tool to observe decoherence New J. Phys. 15 013052

DOI

62
Saif F 2012 Talbot effect with matter waves Las. Phys. 22 1874

DOI

63
Grossmann F Rost J M Schleich W P 1997 Spacetime structures in simple quantum systems J. Phys. A 30 L277

DOI

64
Marzoli I Saif F Bialynicki-Birula I Friesch O M Kaplan A E Schleich W P 1998 Quantum carpets made simple Acta Phys. Slovaca 48 323

DOI

65
Friesch O M Marzoli I Schleich W P 2000 Quantum carpets woven by Wigner functions New J. Phys. 2 4

DOI

66
Belloni M Doncheski M A Robinett R W 2005 Zero-curvature solutions of the one-dimensional Schrödinger equation Phys. Scr. 72 122

DOI

67
Bialynicki-Birula I Mycielski J 1975 Uncertainty relations for information entropy in wave mechanics Commun. Math. Phys. 44 129

DOI

68
Iqbal S Saif F 2008 Quantum computing with particle in a driven square well J. Russ. Laser Res. 29 466

DOI

Outlines

/