Quantum implications of non-extensive statistics

Exploring the analogy between quantum mechanics and statistical mechanics we formulate an integrated version of the Quantropy functional [1]. With this prescription we compute the propagator associated to Boltzmann-Gibbs statistics in the semiclassical approximation as $K=F(T) \exp\left(i S_{cl}/\hbar\right)$. We determine also propagators associated to different non-additive statistics; those are the entropies depending only on the probability $S_{\pm}$ [2] and Tsallis entropy $S_q$ [3]. For $S_{\pm}$ we obtain a power series solution for the probability vs. the energy, which can be analytically continued to the complex plane, and employed to obtain the propagators. Our work is motivated by [4] where a modified q-Schr\"odinger equation is obtained; that provides the wave function for the free particle as a q-exponential. The modified q-propagator obtained with our method, leads to the same q-wave function for that case. The procedure presented in this work allows to calculate q-wave functions in problems with interactions; determining non-linear quantum implications of non-additive statistics. In a similar manner the corresponding generalized wave functions associated to $S_{\pm}$ can also be constructed. The corrections to the original propagator are explicitly determined in the case of a free particle and the harmonic oscillator for which the semi-classical approximation is exact.


INTRODUCTION
Nonextensive entropies depending only on the probabilities have been obtained in [1]. They belong to a family of nonextensive statistical mechanics, relevant for nonequilibrium systems. Renowned examples are Tsallis (S q ) [2,3] and Sharma-Mital [4]; all of them can be obtained within the framework of Superstatistics [5].
For the entropies depending only on the probability, there are two entropy functionals [1]: where the index l runs over the states of the system and p l denotes the probability of the state l. These expressions can be considered as building blocks for nonextensive entropies without parameters. For example, one can consider S 1 (S + + S − )/2. These entropies are noticeably distinct to Boltzmann-Gibbs (BG) entropy for systems with few degrees of freedom; however, when the number of degrees of freedom goes to the thermodynamic limit, they match perfectly with BG statistics Cabo [6]. This is when the probability p l is small, the leading term in S + and S − expansions is the BG entropy. Therefore, in this limit, BG statistics is recovered. They belong to the class of Superstatistics of [5] with an intensive parameter χ 2 distribution found in [1]. These entropies have also been studied in [6][7][8][9][10].
There is a universality of the Superstatistics family [5]. As it has been shown, for several distributions of the temperature, the Boltzmann factor essentially coincides up to the first expansion terms. This has as a consequence that also the entropies associated to these Boltzmann factors have all of them basically the same first corrections to the usual entropy. Furthermore, the three entropies listed here that depend only on the probability are expanded only on the parameter y plnp; this is always smaller than 1 giving correction terms to the entropy which at any order are smaller than the previous ones. So, that any function of y proposed as another generalized entropy, depending only on this parameter, when expanded in y will basically coincide with one of the three ones studied here; clearly demanding that the first term in the expansion is −y giving BG entropy. Thus, the entropies S + , S − , and their linear combinations can be considered as building blocks to compute any possible modified entropies depending only on the probability.
We are motivated by the concept of Quantropy developed by Baez and Pollard [11] and by nonlinear quantum systems with modified wave functions based on Tsallis statistics in [3,12,13]. For example, the work [13] developed a nonlinear quantum mechanics with q-mathematics motivated by Tsallis entropy. In recent years, there have also been other interesting developments in the connections of nonextensive entropies and quantum mechanics [14][15][16][17][18][19][20][21][22][23][24][25]. Also, the work [26] showed extensions of nonlinear quantum equations arising from an effective one particle treatment of many-body physics, such that the nonlinearity represents the interactions, obtaining wave function solutions that are q-distributions and including the harmonic oscillator potential. There exists as well a connection between nonlinear quantum equations and nonlinear diffusion and Fokker-Planck equations [18,20] that also is noticed in [26]. Moreover, interesting applications of nonextensive entropies to compute statistical and thermodynamical properties of graphene and 2-dimensional quantum structures [27][28][29][30] have been developed. We develop a version of Quantropy in terms of the propagator of a quantum mechanical theory. Our generalized propagators could be connected to the appropriate quantum equations. Baez and Pollard's Quantropy is a functional of the amplitude on the path integral a, with the same functional form as the entropy in terms of the probability Q − X a(x)lna(x)dx. Giving the functional where α and λ are Lagrange multipliers and x is a path in the space of all possible paths X. From the search of extrema of this functional, restricted to values of a normalized and an average of the action, Baez and Pollard obtained the relation a exp(iS/h − 1 − α) with λ 1 ih . Then, a ∼ exp(iS/h ) with the normalization fixed by the Lagrange multiplier α. This scheme deepens on the relation between Quantum Mechanics and Statistical Physics, which has also been studied in different approaches [31,32]. For example, in [33], the Fisher Information measure is employed to explore this relation.
In Baez and Pollard's approach, the energy is mapped to the action S and the temperature to ih . We consider the same identification but instead we identify E with the classical action S cl . Thus, we consider as the analog of the entropy a functional in terms of the propagator, instead of the amplitude a. This is an extrapolation of the Quantropy [11] to an integration over all classical paths. It is worth to mention that the analog to the microstate in statistical mechanics is a particle path in quantum mechanics. Such that as the partition function in statistical mechanics is the sum over all the microstates, the quantum mechanical analog is the sum over all the paths of the particle (Feynman path integral). The standard expression for the propagator is given semiclassically by K ∼ exp(iS cl /h ). We use this fact to define a kind of integrated Quantropy functional now in terms of the propagator for BG statistics, which we extend to the modified statistics S + , S − , and S q .
This article is organized as follows. In Section 2, we obtain a series expansion for the probabilities versus βE for the generalized entropies depending only on the probabilities S + and S − . In Section 3, we continue these expansions to the complex plane. In Section 4, we present a version of Quantropy for BG statistics, S + and S − and S q . In Section 5, we study in particular the case of the free particle propagator, obtained from the extrema of the Quantropy in the cases of S + and S − and S q for q < 1 and q > 1. We show that the K q propagator results exactly in the q-exponential that defines the q-wave function for the free particle [3]. In a similar manner, we argue that the corresponding generalized propagators K + and K − provide us with a procedure to construct Ψ + and Ψ − for the free particle. Our method however gives the possibility to construct K q , K + , and K − also for problems with interactions and by this mean to identify the corresponding wave functions. We also provide a way to perform the normalization inspired in Feynman and Hibbs work [34]. Section 6 is devoted to the analysis of the harmonic oscillator, and we exemplify with the case corresponding to K + . Finally, in Section 7, we study the K + propagator for the particle in an infinite potential barrier. In Section 8, we summarize our results and present the conclusion. In a Supplementary Appendix we present a numerical study of the propagators. the contrary to BG statistics, for a system subjected to S + extremization, probability normalization, and energy conservation, there is not a simple inverse function of the probabilities p vs. the values of the energy state E. We overcome this difficulty by finding a series solution to the extremum equation. There are other possible series solutions, but we discuss here one that has a good convergence. At the end of the section, we give also the probability expansion for the entropy S − , which is obtained by an equivalent Ansatz.
The functional to maximize the S + entropy subjected to probability normalization and averaged energy is given by [7,8]: (2) β and γ are Lagrange multipliers and E l is the energy of the state l, with probability p l . The average values of energy and the normalization value have been omitted for simplicity. The extrema of (2) given by δΦ + δp l 0 gives a relation between the energy E and the probabilities p (we have omitted the index l): Notice that we omit the subindex l from the quantities. Setting the Lagrange multiplier γ to −1, we first expand the previous equation around p 0 that accounts to consider the expansion around y plnp 0. That is, for the exponential of minus equation (3), one gets e −βE p − p 2 lnp 2 + 1 2p 3 ln p 2 + 2lnp 3 + ln p 4 +1 6p 4 − 3lnp 2 − 8lnp 3 − 9lnp 4 − 6lnp 5 − ln p 6 +1 24p 5 12 lnp 2 + 44 lnp 3 + 70 ln p 4 + 68ln p 5 + 42ln p 6 + 12ln p 7 + ln p 8 + ... .
where c n can be functions of βE. Plugging (5) in (4), to have the LHS equal to the RHS, the coefficients multiplying the powers of e −nβE with n > 1 have to vanish. This gives a recurrent expression for the coefficients c n which for the first four coefficients is solved as We have denoted βE as x. The coefficients (6) give the following approximate solution for the probabilities versus βE: In Figure 1, we compare the exact value of p vs. βE with the power series solution (7) till 3rd order and with the Boltzmann distribution e −βE .
Probability Expansion for the Entropy S − -Consider the other generalized entropy dependent only on the probabilities S − . For this entropy, the functional to extremize Φ − reads β and γ are Lagrange multipliers, and E l is the energy of the state l and p l its probability. Finding the extrema of (8) as δΦ − δp l 0 and proposing the same Ansatz (5), we obtain a set of equations that can be solved to give the recursive probability solution:

MODIFIED AMPLITUDE EXPANSIONS
In this section, we use the series solutions for the probabilities in terms of the energy obtained in the previous section, to perform an analytic continuation to the complex plane. Considering a as the amplitude of a path, this is a new complex variable substituting the probability p l , and A as the action replacing βE l . This identification will allow to study modified Quantropy functionals, for the definition of Baez and Pollard (equation (1)), as well as our definition (18). The usual Quantropy solution will give an exponential a ∼ e iA h . In our approach, this would be the propagator K ∼ e iS cl h . We want to analyze the new statistics S + and S − . We will find a functional dependence of a vs. A (K vs. S cl ) that deviates from the exponential dependence.
The main idea is to complexify first the power expansion solution (7) since the amplitude is a complex number, such that we have a solution to the extrema of the modified Quantropy. The functional to extremize reads Finding the extrema of (10) w.r.t. a, i.e., solving δΦBP,+ δa(x)  standard propagator (yellow line), for the free particle for the modified statistics S + . We set the mass and the Planck constant to unity. The quantum regime is given by S cl ≈ h . Imposing S cl (h which translates for fixed x 1 in ta1/2, for fixed t 1 translates in x 2 (2. Notice that in the quantum regime S cl ≈ h , there are differences between the standard and the modified propagator. In the classical regime, there are many oscillations caused by the series expansion that should sum up when computing more terms. When S cl ≪ h , then both results coincide. The difference between the results of the modified propagator and the standard free particle is that the modified propagator result could be interpreted as the particle with an effective potential; this would give a spatially bounded wave function.

0, one gets
Frontiers in Physics | www.frontiersin.org May 2021 | Volume 9 | Article 634547 4 The range of validity of the propagators expressions depends on the convergence of the imaginary series solution to this equation. The series is obtained by doing the replacement βE by A ih and p by a in (3). Also, the Lagrange multipliers have to be mapped: β to λ and α to γ. The minus sign gives the right sign after a rotation on the argument of the exponential (similar to a Wick rotation). As the series solution continuation of (7), we obtain the following expression: Since A has units of action, the argument of the exponentials and the terms on the expansion are adimensional. Substituting this expression on the constraint equation (11), we obtain the real and imaginary parts of F a A ih . In Figure 2, the relevant difference of the propagators K + obtained employing (12) with respect to the standard one can be observed in the region of S cl ≈ h . Using the parameter λ 1 ih in (12), the expression for the amplitudes becomes It is not difficult to observe that all terms of this expansion can be written as derivatives with respect to the parameter λ. If we derive with respect to λ, the usual amplitude we obtain is z zλ e −λA −Ae −λA . Higher derivates can be written as where m and n are positive integers. Thus, we rewrite (3) as One can compute the corrections to any order. Those corrections to the usual amplitude a 0 e iS/h can be interpreted as higher order interactions of the action at different frequencies 1 . Now, one can apply the same method to determine the distribution arising from the Quantropy with statistics S − . We also have to perform the extension to the complex plane. The amplitude distribution for the modified Quantropy coming from S − is given by the following:

QUANTROPY IN TERMS OF THE PROPAGATOR
In this section, we present as an alternative proposal a kind of integrated version of the Quantropy of [1]. First, we do it for the BG entropy and then for S + , S − , and S q . The change in distribution probabilities which arise from modified entropies in statistics is now reflected in the quantum arena as modifications to the propagators. The propagators between points in space-time (x a , t a ) and (x b , t b ) in quantum mechanics determine the probability amplitude of particles to travel from certain position to another position in a given time. As modified entropies in statistical physics lead to modified probability distributions, distinct probabilities of propagation over all paths from (x a , t a ) and (x b , t b ) will arise from a modified Quantropy. In the work [11], the Quantropy functional associated with BG statistics was formulated, and its maximization leads to the weight on the path integral a ∼ exp(−λS) with λ 1/ih . We propose another functional, which in a sense constitutes an integrated version of Quantropy. Its maximization leads to the propagator K(x) ∼ exp(−λS cl (x)). The Wentzel-Kramers-Brillouin (WKB) method [35] allows to compute the wave function in a semiclassical approximation. In a sense, this is linked to our approach, in which we maximize a functional which determines the propagator with a semiclassical approximation in terms of the classical action S cl . This is exact for the free particle, the harmonic oscillator as well as for other cases [36,37]. The procedure is applied to generalized entropy functionals, giving a modified propagator. For the Tsallis statistics, we obtain K q (x) ∼ exp q (−λS cl (x)). This structure is the same as the one of the wave function for the free particle of the Tsallis statistic Ψ q (x) exp q (i(kx − wt)), which has been proposed as solution of the nonlinear quantum equations of [13]. According to Feynman arguments, one can start with the free particle propagator and determine the corresponding wave function, as discussed in [34]. Thus, our procedure allows to find a propagator which can be identified with the wave function of interest. We should note that the propagator resulting from our procedure will describe not only the free particle but to a good approximation any other problem with its corresponding classical action. Our method should give the wave function solution for the problem of interest.
With the same method, we write functionals for S + and S − and obtain probability distributions; we can write the corresponding Quantropies and obtain the propagators K + and K − and correspondingly extrapolate them to the wave functions Ψ + and Ψ − . This would give us the quantum behavior for a given action. The obtained propagators can be related to nonlinear quantum systems studied in the literature [26].
To define our functionals, we use the semiclassical limit to compute the propagator; this is K(x) F(a, b)e iS cl (x) h , denoting the classical action as S cl (x) and being F(a, b) a constant depending on the time difference t b − t a . For the free particle and the harmonic oscillator as well as other physical problems [36,37], this is an exact result.
The integrated Quantropy functional for the new S + statistic is studied in [1,[6][7][8][9][10] is given by the following: The extrema condition δΦ+ δK(y) 0 with functional derivatives leads to the equation: and by changing the variable notation y to x, this can be written as Using our knowledge to solve this type of equation from the statistical physics case, presented in Section 2, this gives for the modified propagator the following series solution: where N + is the normalization constant. This is obtained by taking the normalization α −1. A different normalization would change the coefficients in the expansion (21).
The maximization constraint for the new S − statistic is given by The extrema condition δΦ − δK(x) 0 gives the equation: Using our knowledge of this type of equation from the statistical physics case, we obtain for the modified propagator the series solution: N − is the normalization constant. In the case of Tsallis statistics, the functional is given by the following: and the solution to We have still to discuss the normalization of the different Kernels. This q-propagator is related to the q-wave function for the free particle nonlinear quantum mechanics of [3]. We explore this case which has been studied by other means in the literature [3,13,24].

FREE PARTICLE PROPAGATORS
In this section, we write a modified propagator up to third order for the free particle in the case of the statistics S + , S − , and S q for q 1 − δ and q 1 + δ with δ > 0. The values of q less or equal than one are considered in order to compare the different propagators. We determine when the corrections to the usual propagator play an important role which turns to be in the quantum regime characterized by S cl ≈ h . First, we describe the procedure; then, we describe the normalization; and in the last subsection, we summarize our results.

Superposition of Kernels
Now, we proceed to describe a generalized Kernel. The generalized complex probability distribution given by expansion (21) can be regarded as a superposition of Kernels. Furthermore, the superposition will carry to the wave functions. In order to normalize the superposition, we consider that the total Kernel expansion integration is the same as the usual (1 for the free particle), as is explicit in the Quantropy functional (4). We show that this coincides with the result for the normalization obtained from propagating the wave function [34]. For the free particle, the unnormalized Kernel is as follows: where n is the number of divisions of the time interval and ϵ is an infinitesimal time parameter that satisfies t ϵn. This expression arises from computing the path integral to get the following: with A m 2ϵh . This calculation is performed dividing the integral in multiple Gaussian integrals (see [34]). The normalization constant is given by N 2πih ϵ m − n 2 ; hence, the normalized propagator reads We define the unnormalized Kernel for the free particle as and the first two corrections in K + are given by Thus, the generalized Kernel associated with S + entropy is given by K + (x, t) N + (k(x, t; 1) + k(x, t; 2) + k(x, t; 3) + . . . ).
The normalization constant is determined by the requirement ∞ −∞ K + (x, t)dx 1, and up to the first corrections is given by 0.883.... The reason for this normalization is also understood by an argument presented in the following, motivated by Feynmann and Hibbs procedure [34].
Let us also discuss the normalization of the modified propagator with respect to the usual one, as shown in the standard case [34]. We start considering the original unnormalized Kernel for the free particle computed from the path integral: To determine the normalization constant in the Feynman and Hibbs method, we can apply formulas  and (4-3) on their book [34] to write the new infinitesimal Kernel between position x i and x i+1 , with Δx i x i+1 − x i , in a time ϵ as follows: The method consists in writing the wave function at a position x at a time t + ϵ in terms of the wave function at position y x + η at a time t, explicitly In the quantum standard theory, the normalization constant can be determined by expanding the LHS of (31) ψ(x, t + ϵ) ψ(x, t) + ϵz t ψ and the RHS ψ(x + η) ψ(x, t) + ηz x ψ + η 2 2 z 2 x ψ and exp(−iϵV/h ) 1 − iϵV h + . . .; then, we compare the leading term ϵ 0 . This implies that In a similarly fashion, one gets for the first correction to K + written in (30): It is worth to mention that the more important contribution to (31) is given for small η′s, as well as in our generalized case. It is necessary to check this argument; in order to verify, let us consider the following integrals: The first correction gives the following relation: The previous normalization factor is a general feature to apply to any potential V(x, t), in particular is valid for the cases discussed here. A similar expression holds for the normalization of K − normalization, and this will be calculated in the next section.

Analysis of the Propagators
Here, we summarize the propagators obtained with the normalization methods described in previous subsections. The results for K + can be extrapolated to K − and K q because the method applied to obtain all of the propagators is the same.
Recall the standard propagator of the free particle from the space-time point (0, 0) to (x, t) is given by The constant w.r.t. x is as follows: N 0 m 2πih t . For the case of S + and S − statistics, the first two contributions to the modified propagators (21) and (24) read: with N + ∼ N 0 . For the Tsallis statics, the associated propagator (26) is given by the expression: We calculate the normalization constants for K ± up to the first correction and exactly K q to get the following: In the quantum regime S cl ≈ h , the differences between the propagators K + , K − , K q , and K 0 are shown in Figures 2-4. In Figure 2 the K − propagator is compared to the usual one. Also, Figure 3 shows a comparision of the K − propagator with the standard one. The last Figure 4 shows a comparison between the propagators for the S q statistics K q for q < 1 and q > 1 with the usual one. The region of interest is the quantum regime with S cl ≈ h . Furthermore, in the classical regime, the oscillations of the standard propagator grow averaging to zero as discussed in [34]. Thus, we are interested in comparing the corrections arising from different statistics in the quantum region of interest and, in the plots, the differences between the standard and the modified propagator can be observed. In the classical regime, there are oscillations that should sum up when computing more terms. When S cl ≪ h , then both results coincide as shown in Figure 5. The modified propagators could be interpreted as describing a particle with an effective potential; this would give a spatially bounded wave function.

THE HARMONIC OSCILLATOR
In this section, we apply the formulation of our modified Quantropy of Section 4 for the case of the harmonic oscillator. We compute the modified propagator constructed by a superposition as it was done previously. The extension of quantum systems employing the modified q-statistics has been made only for the case of the free particle [13] with different arguments. Our proposal allows to search the manifestation of nonextensive statistics in nonlinear quantum systems for generic potentials. We illustrate the procedure calculating only K + , and the other propagators K − and K q could be similarly calculated.
For the harmonic oscillator with Lagrangian L m 2 _ x 2 − mω 2 2 x 2 , the path integral Kernel reads Next, following a similar procedure as for the free particle, we compute the generalized Kernel and normalize it. The unnormalized Kernel is given by the following: Now, we compute the next terms of the Kernel up to third order, which are as follows: Thus, the total normalized propagator up to third order reads where N + 1 cos(ωT) . The relative normalization of the modified propagator with respect to the normalization of the usual propagator is a result obtained in Section 5 [see formulas (31)(32)(33)]. This result is universal, i.e., independent of the action. In Figure 6, we compare the propagator for the harmonic oscillator for the standard Quantropy and for the one based on S + and S − statistics. We are interested in the quantum regime given by S cl ≈ h . There are noticeable effects in that regime. Outside the quantum region, oscillations grow as in the usual case [34]. This behavior occurs in the classical region in which the modified Kernels will also not contribute. In this section, we apply the formulation of Quantropy developed previously to compute the propagators associated to a particle in an infinite potential barrier given as follows: The standard unnormalized propagator for this problem is given as follows [38]: where x 0 is the initial position. It is important to specify the initial position since the particle cannot be located at x ≤ 0; thus, we set x 0 ϵ, where ϵ is a small positive non zero parameter. This allows to make a comparison with the free particle case. Notice that if we set x 0 0, the propagator vanishes since the x dependence is quadratic. Now, we compute the nonlinear propagator associated to the statistics S + : K + (x, t; ϵ, 0) up to third order, which is given by expression (24). Since our theory is nonlinear, the superposition principle is not valid, i.e., we cannot consider the difference between the propagators of two free particles, like is done for the usual propagator K(x, t; x 0 , 0), rather in the semiclassical regime, we consider the logarithm: This allows to substitute S cl into the series expansion (21) K + (x): Note that the propagator is ill-defined if ϵ is zero. This propagator oscillates faster than the usual propagator in x, and its amplitude is greater. However, the global behavior of both propagators is quite similar, and the oscillations and the amplitude grow as x increases. For the time dependence in both cases, the propagator tends to zero as t grows (see the set of graphics in Figure 7).
In the usual case, the propagator K(x, t; x 0 , 0) has the following interpretation, the first part corresponds to the classical path of the free particle from (x 0 , 0) to (x, t), while the second part corresponds to the classical path of a free particle from (x 0 , 0) bouncing off the wall and going to (x, t). The modified propagator expression suggests a similar interpretation, i.e., the whole propagator can be considered as the sum of both classical paths with the leading terms given by the standard free particle and nonlinear corrections which can be interpreted as an effective potential.

FINAL REMARKS
In this work, we explore the novel concept of Quantropy in Quantum Mechanics (Q.M.), which constitutes the analog of the entropy in Statistical Mechanics (S.M.). Mathematically, Quantropy can be regarded as an analytical continuation of We establish a new definition of Quantropy, with the energy mapped to the classical action E → S cl , i.e., we consider that the main entity is the propagator K(x) instead of the amplitude a(x) of the path. Thus, we construct the propagator, a kind of integrated version of the Quantropy Q 0 − X K(x)lnK(x)dx; in this way, the functional corresponding to the BG under the maximization procedure leads to K 0 ∼ exp(iS cl /h ).
We applied this concept to find the quantum mechanical implications of modified entropies S + , S − , and S q , and the associated quantropies lead to generalized propagators, which imply a modified wave functional quantum mechanics (modified Schroedinger equation). Let us point out that the small probabilities limit of modified entropies S + and S − leads to Boltzmann-Gibbs entropy. Similarly in the limit of small path amplitude, one recovers the original Quantropy from the modified functionals. This formalism could be a novel framework to study nonlinear quantum mechanics as those consider in [4,12,13]. We also provide an understating of modified propagators associated with Tsallis statistics, leading to wave functions corresponding to q-distributions. Also, the result for the Tsallis statistics implies a propagator K q ∼ exp q (iS cl /h ), where exp q (x) is the q-exponential; this is relevant since it makes contact with the Tsallis result of modified wave function for the free particle, whose solution is ψ q ∼ exp q (i(kx − ωt)). The connection is due to the arguments of [34] in the discussion of the propagator for the free particle K 0 ∼ exp(iS cl /h ). They show that the propagator K 0 corresponds to the free particle wave function ψ 0 ∼ exp(i(kx − ωt)). Thus, analogously K q will lead to ψ q . As a further work, we need to explore the relations in the case of the modified propagators K + and K − . They will give rise to wave functions whose dependencies will be given by Ψ ± exp ± (i(kx − ωt)), also for the free particle. In this case, we have a recurrent series solution but we do not have exact expressions for these generalized exponentials. As discussed, our proposal provides also generalized propagators K + , K − , and K q for problems with interactions; we illustrated this by considering the K + associated with the harmonic oscillator and the infinite potential barrier.
There are hints from previous studies that the modified entropies considered here can be interpreted as linked with modified effective potentials. Therefore, these modifications to the free particle could be related to a usual quantum mechanics with an effective potential [9]. However, these effects could also lead to nonlinear quantum equations explored in the literature with modified wave functions [3,12,13,26]. Furthermore, what we found here based on the concept of Quantropy could be linked to results for quantum systems in terms of usual entropy vs. the density matrix [10]. A system governed by a modified statistics (S + , S − or S q ) will lead to modified density-matrix distributions. Moreover, the modified "propagators" K + , K − , and K q actually are strictly no longer standard propagators because they lack the usual propagation property. This means that is not equivalent to propagate the particle from (0, 0) to (t 2 , x 2 ), than to first propagate it from (0, 0) to (t 1 , x 1 ) and then from (t 1 , x 1 ) to (t 2 , x 2 ). This occurs because, for example, K + (0, 0; x 1 , t 1 )K + (x 1 , t 1 ; x 2 , t 2 )dx 1 ≠ K + (0, 0; x 2 , t 2 ).
This relates to the fact that in a quantum open systems where these generalized entropies are motivated, the nature of the processes is non-Markovian. Those systems in consideration are modeled with Master Equations (Stochastic) [39]. We consider that this formalism could be a natural framework to study nonlinear quantum mechanics.
We would like to explore further processes where the modified statistics in Quantropy play a central role. This could be done via modified wave functions, which could be interpreted as the usual quantum mechanics with an effective interaction [40] or from nonlinear quantum equations. The modified wave functions will correspond to the modified propagators obtained in this work. In this work, we obtained the modified propagators for the free particle, harmonic oscillator, and the infinite potential barrier associated to the different statistics K + , K − , and K q , discussing the associated quantum behavior. We would also like to explore in future work the quantum mechanical evolution of other physical systems.

DATA AVAILABILITY STATEMENT
The original contributions presented in the study are included in the article/Supplementary Material; further inquiries can be directed to the corresponding author.