The Generalized OTOC from Supersymmetric Quantum Mechanics: Study of Random Fluctuations from Eigenstate Representation of Correlation Functions

The concept of out-of-time-ordered correlation (OTOC) function is treated as a very strong theoretical probe of quantum randomness, using which one can study both chaotic and non-chaotic phenomena in the context of quantum statistical mechanics. In this paper, we define a general class of OTOC, which can perfectly capture quantum randomness phenomena in a better way. Further we demonstrate an equivalent formalism of computation using a general time independent Hamiltonian having well defined eigenstate representation for integrable supersymmetric quantum systems. We found that one needs to consider two new correlators apart from the usual one to have a complete quantum description. To visualize the impact of the given formalism we consider the two well known models viz. Harmonic Oscillator and one dimensional potential well within the framework of supersymmetry. For the Harmonic Oscillator case we obtain similar periodic time dependence but dissimilar parameter dependences compared to the results obtained from both micro-canonical and canonical ensembles in quantum mechanics without supersymmetry. On the other hand, for one dimensional potential well problem we found significantly different time scale and the other parameter dependence compared to the results obtained from non-supersymmetric quantum mechanics. Finally, to establish the consistency of the prescribed formalism in the classical limit, we demonstrate the phase space averaged version of the classical version of OTOCs from a model independent Hamiltonian along with the previously mentioned these well cited models.


Abstract
The concept of out-of-time-ordered correlation (OTOC) function is treated as a very strong theoretical probe of quantum randomness, using which one can study both chaotic and non-chaotic phenomena in the context of quantum statistical mechanics. In this paper, we define a general class of OTOC, which can perfectly capture quantum randomness phenomena in a better way. Further we demonstrate an equivalent formalism of computation using a general time independent Hamiltonian having well defined eigenstate representation for integrable supersymmetric quantum systems. We found that one needs to consider two new correlators apart from the usual one to have a complete quantum description. To visualize the impact of the given formalism we consider the two well known models viz. Harmonic Oscillator and one dimensional potential well within the framework of supersymmetry. For the Harmonic Oscillator case we obtain similar periodic time dependence but dissimilar parameter dependences compared to the results obtained from both microcanonical and canonical ensembles in quantum mechanics without supersymmetry. On the other hand, for one dimensional potential well problem we found significantly different time scale and the other parameter dependence compared to the results obtained from non-supersymmetric quantum mechanics. Finally, to establish the consistency of the prescribed formalism in the classical limit, we demonstrate the phase space averaged version of the classical version of OTOCs from a model independent Hamiltonian along with the previously mentioned these well cited models.
Keywords: OTOC, Supersymmetry, Out-of-equilibrium quantum statistical mechanics. 1 Introduction The concept of out-of-time-ordered-correlators (OTOC) first introduced by the author duo, Larkin and Ovchinnikov to describe the semi-classical correlation in the context of superconductivity [1], which was mostly used in various condensed matter systems to study various out-of-equilibrium phenomena in the quantum regime [2]. However, recently it has attracted the attention of theoretical physicists from other branches in very different context finding applications in the finite temperature extension of quantum field theories, bulk gravitational theories, quantum black holes, and many more sensational topics in the list [3][4][5][6][7][8]. It is considered to be one of the strongest theoretical probe for quantifying quantum chaos in terms of quantum Lyapunov exponent [9], quantum theories of stochasticity and randomness among the theoretical physics community. Besides playing a key role in investigating the holographic duality [10][11][12][13] between a strongly correlated quantum system and a gravitational dual system, it also characterizes the chaotic behaviour and information scrambling [14][15][16][17][18][19] in the context of many-body quantum systems [21][22][23]. The detailed study of OTOCs reveal an intimate relationship between three entirely different physical concepts, namely holographic duality, quantum chaos and information scrambling.
The key idea of OTOCs can be best understood as the growth of the non-commutativity of quantum mechanical operators 5 which are defined at different time scales, and hence can be described using the Poisson Brackets for its classical counterpart. Not only that, but also the quantum mechanical thermal ensemble average or equivalently the quantum mechanical trace operation can be described by using the phase space average in the classical limit. It is considered as the quantum mechanical analogue of the classical sensitiveness to the initial conditions in the time dynamics of a quantum system. The exponential growth of these correlators indicates the presence of chaos in the quantum system, which has led to discussions of the "butterfly effect" in black holes [24][25][26] with a saturation bound on the quantum Lyapunov exponent and for various spin models [27,28].
There has been a growing interest to understand the behaviour of OTOC even for systems where quantum chaos is expected to be absent, the most relevant example being the study of OTOC in the quantum Ising spin chain model, where power law growth of OTOC's is observed as opposed to the exponential growth in non-integrable models in support of non-chaoticity [29][30][31]. Another interesting revelation came from a recent study of OTOC's for a quantum system with discrete energy levels, weakly coupled to a non-adiabatic dissipative thermal environment. This type of system is commonly known as open quantum system (OQS), where the OTOC was found to saturate exponentially in contrast to the exponential growth for a quantum chaotic system [32]. OTOC's has been of prime theoretical interest for diagnosing also the rate of growth of chaos with respect to the different time scales involved in the quantum system through the operators and hence for studying the scrambling of quantum information in black holes and strongly correlated quantum mechanical systems. It serves as a strong theoretical probe for investigating various bulk gravitational dual theories in the framework of AdS/CFT correspondence. Among others, the existence of shock waves inside black holes [33][34][35] and the maximum saturation bound of the quantum version of the Lyapunov exponent are the most famous examples where out of time ordered correlation functions have proven to be useful in the Ads/CFT correspondence. This maximum saturation bound is famously known as the Maldacena-Shenker-Stanford (MSS) bound [36]. The SYK model  is a well known model which saturates this well known bound and shows the signature of maximal chaos both in its 1D CFT and gravitational dual 2D balck hole counterpart. In the ref. [68] this bound on quantum Lyapunov exponent was further generalized for many body systems using the well known Eigenstate Thermalization Hypothesis (ETH). Very recently in the ref. [69], the author used the tools and techniques of computing OTOC in the context of Cosmology by following the underlying slogan Cosmology meets Condensed Matter Physics to study the quantum mechanical correlation functions from random primordial fluctuations appearing in the context of cosmological perturbation theory of background spatially flat FLRW metric. Specifically, these fluctuations are appearing in the context of stochastic particle production during inflation, during the epoch of reheating and also for the primordial phenomena which is governed by the quantum generalization of Brownian motion i.e. for the cosmological epochs in the time line of the universe where the physics of out-of-equilibrium phenomena play significant role.
In ref. [70], the earlier discussed fact regarding the exponential growth of the OTOCs in the associated time scales to describe the chaotic fluctuations for non-integrable systems has been established for various well known quantum mechanical systems. A study of the same for integrable models suggest non-chaotic quantum mechanical fluctuations in the quantum regime. In the present context, the phrase quantum randomness describes a physical phenomena which describes chaotic or non-chaotic i.e. in principle any random behaviour of a system with respect to the associated time scales of the system. For the physical systems such quantum randomness can be described by the following two-fold formalisms: 1. Formalism I: The first approach is based on the construction and the mathematical from of the solution of the Fokker Planck equation, using which various stochastic phenomena in the quantum out-of-equilibrium regime can be studied. One of the famous ex-ample is the stochastic cosmological particle production phenomena which can be directly mapped to a problem of solving Schröinger equation with an impurity potential within the framework of quantum mechanics, which is actually describing propagation of electrons inside an electrical wire in presence of an impurity or defect. Within the framework of quantum statistical mechanics instead of solving the Schröinger problem directly or may be solving the dynamical equation for the quantum mechanical fluctuations during the stochastic particle production one can think about a cumulative probability distribution function of of this stochastic process, P(n, τ ) which depend on two crucial quantities, which are -the number density of the produced particle and the associated time scale of the stochastic process. Using a detailed physical arguments and computation one can explicitly show that this probability distribution function, P(n, τ ), satisfy Fokker Planck equation, which is given by: At infinite temperature (β → 0) : 1 µ ∂P(n, τ ) ∂τ = n(n + 1) ∂ 2 P(n, τ ) where µ represents the mean stochastic particle production rate and V (n) associated potential which is only significant at finite temperature. By solving these set of equations in presence of appropriate initial conditions one can get to know about the related profile of the stochastic process semi classically in the present context. Just not only that, but also one can treat these versions of the Fokker Planck equation as the semi classical statistical moment generating equations because by replacing the profile function P(n, τ ) with the appropriate moment generating function F(n) one can compute all the moments. To serve this purpose one need to use the following fundamental equation: Statistical average of moment generator : F(n) := dn F(n) P(n, τ ), (1.3) which physically represents the expectation value or the statistical average value of the number density dependent moment generating any arbitrary function F(n).
Further substituting the appropriate form of this function in the moment dependent Fokker Planck equations one can explicitly compute the expression for all the physically relevant statistical moments i.e. n , n 2 , · · · explicitly without explicitly knowing about the particular mathematical structure of the profile of the related stochastic dynamical process at infinite or finite temperature for a given structure of number density potential function. These moments are extremely important in the present context of study as all of them semi classically compute the expressions for all the equal time quantum correlation functions required to study the out-of-equilibrium aspects, such as stochastic effects, disorder, random fluctuations etc. both at infinite and finite temperatures. See the references [2,[71][72][73][74] where all authors have studied the physical impact of this formalism to describe out-of-equilibrium aspects in various different contexts.
Now let us speak about some applicability issues related to this particular formalism. Since this formalism only allows us to know about the effect of the semi-classical correlations at equal time that might be not very interesting when we are actually thinking of doing the computation of the correlations and its rate of change at different time scales associated with the quantum mechanical system of study. For an example, if we are interested to compute the any general N -point semi-classical correlation function as defined by the following expression: n(τ i ) = n(τ 1 )n(τ 2 ) · · · n(τ N ) ∀ τ i (i = 1, 2, · · · , N ) contain disorder, (1.4) then this particular formalism will not work, as using this formalism one cannot capture the effect of disorder effect in the associated time scale of the system. On the other by following the usual tools and techniques one can only compute the above mentioned N -point correlators either in time ordered sense, where τ 1 > τ 2 > · · · > τ N or in the anti-time ordered sense, where τ 1 < τ 2 < · · · < τ N . So from the technical ground defining this N -point correlators including the effect of disorder in the time scale at any arbitrary temperature is also a very important topic of research and here comes the crucial role of the next formalism where we can allowed to define and explicitly compute such quantum effects at the level of correlation functions. Another important aspect we want to point out here that, the present formalism don't bother about any Lagrangian or the Hamiltonian formulation of the associated quantum mechanical system under consideration. So if we really interested to know about the effect of time disordering in the expressions for the quantum mechanical correlation functions which are defined in terms of the fundamental operators appearing in the quantum version of the Lagrangian or the Hamiltonian of the system under consideration and also want to explicitly know about time variation with respect to different time scales associated with the system, which are actually the source of time disordering, then instead of the present formalism it is obviously technically correct and easier to think about the implementation of the second formalism, which gives us the better understanding of time scale disordering. In the next point and in the rest of the paper we will follow the second formalism to compute the quantum correlation functions from the fundamental operators from the quantum mechanical systems under study which can explicitly capture the effect of disordering in the time scale. Not only us, but also the present trend in the research suggesting to make use of the next formalism to get better understanding of time disordering phenomena in quantum mechanical system.

Formalism II:
The second approach is based on finding quantum correlation functions including the time disordering effect and through out the paper we have followed this formalism to study effects of out-of-equilibrium physics in physical systems [69]. The present computational methodology helps us to know more about the underlying unexplored physical facts regarding the quantum mechanical aspects of various stochastic random process where time ordering or anti time ordering is not at all important and instead of that disorder in the time scale can be captured in the quantum correlations at very early time scale. This method not only helps us to know about the early time behaviour of quantum correlations in the out-of-equilibrium regime of the quantum statistical mechanics, but also give crucial information regarding the late time equilibrium behaviour of the quantum correlations of a specific quantum system. However for all the systems in nature the above interpretation of the quantum mechanical aspects of the randomness phenomena are not same. Based on all these types of quantum systems one can categorize the random time disordering phenomena as, A. Chaotic system which shows exponential growth in the quantum correlators, B. Non-chaotic system which shows periodic or aperiodic or irregular random fluctuations in the quantum correlators. The best possible theoretical measure of all such time disorder averaging phenomena for various statistical ensembles, micro-canonical and canonical ensembles are described by out-of-time-ordered correlation (OTOC) function within the framework of quantum statistical mechanics. Let us define, a set of operators, O i (t j ) ∀ i, j = 1, 2 with i = j or i = j possibilities. The time disorder thermal average over statistical ensemble is described by the following expression: where the thermal density matrix ρ β is defined as: In the present context, C (ij) N (t 1 , t 2 ), represent three possible types of OTOC out of which only i = j possibility which will describe only one OTOC have been explored in earlier works in this area. The other two OTOCs which are appearing from the i = j possibility will be explicitly studied in this paper. The prime objective is to incorporate two more type of OTOCs along with the well known other OTOC is to study the all possible signatures of time disordering average from a quantum mechanical system. Our expectation is all these three types of OTOCs can able to describe the more general structure of stochastic randomness or any simple type of random process in the quantum regime. Though, this idea was revived by Kitaev, then followed by Maldacena, Shenker and Stanford (MSS) and many more to study the quantum mechanical signature of chaos, which is actually the i = j case in the above definition, but the mathematical structure of the other two OTOCs represented by the i = j case suggests that any non-chaotic behaviour, such as, periodic or aperiodic time dependent behaviour, any time dependent growth in the correlators which are different from any type of exponential growth, any type of decaying behaviour in the correlators can be explained in a better way compared to just studying the time dependent behaviour from the well known OTOC which are commonly used in the literature. So in short, to give a complete picture of any kind of time disordering phenomena it is better to study all these possible three types of OTOCs to finally comment on the properties of any physical systems in quantum mechanical regime. Few other important things we want to point out here for better understanding the structure of all these OTOCs capturing the underlying physics of disordering averaging phenomena. First of all here we have to mention that in the definition that we have provided in this paper using which we can able to compute three sets of 2N point OTOCs. Though in the further computation we have restricted our study in the paper by considering N = 1 and N = 2 cases, but to study the general time disorder averaging process one may study the any even multipoint (i.e. 2N point) correlation functions from the present definition. Now here the N = 1 case is basically representing a non-zero UTCR and can be treated as the building block of the full computation as this particular case is mimicking the computation of the Green's function in presence of time disordering. More technically, one can interpret that this contribution is made up of two disconnected time disorder averaged thermal correlator. These disconnected contributions are extremely significant if we wait for a large time scale, in literature we usually identified this time scale as a dissipation time scale on which one can explicitly factorize any higher 2N point correlators in terms of the non-vanishing disconnected contributions. For this specific reason one can treat N = 1 case result as the building block of the any higher 2N point thermal correlators. However, for most of the quantum systems the N = 1 case shows random but decaying behaviour with respect to the associated time scales which are explicitly appearing in the quantum operators of the theory. For this reason study of any N > 1 play significant role to give a better understanding of the time disordering phenomena. For this purpose next we have studied the N = 2 case, which represents the four point thermal correlator in the present context of study and one of the most significant quantity in the present day research of this area which can capture better information regarding the time disorder averaging compared to the N = 1 case. In a future version of this work we have a plan to extend the present computation to study the physical implications of N > 2 quantum correlators to better understand the time disorder averaging phenomena. Now we will comment on the technical side of the present formalism using which one can explicitly compute these OTOCs in the present context. First of all we talk about the time independent Hamiltonian of a quantum system which have their own eigenstates with specific energy eigenvalue spectrum. In this case, construction of the OTOCs describing the time disorder thermal averaging over a canonical ensemble is described by two crucial components, the Boltzmann factor on which the general eigenstate dependent spectrum appears and also the temperature independent micro-canonical part of the OTOCs. At the end we need to take the sum over all possible eigenstates, which will finally give a simplified closed expressions for OTOCs in the present context. Due to the appearance of the eigenstates from the time independent Hamiltonian this particular procedure will reduce the job extremely to study the time dependent behaviour of all the previously mentioned OTOCs that we have defined earlier in this paper. In the rest of the paper we have followed this prescription which is only valid for time independent Hamiltonian which have their own well defined eigenstates. For more details see the rest of the computations and related discussion that we have studied in this paper. Most importantly, using this formalism we can compute all of these OTOCs in a very simple model -independent way. The other technique is more complicated than the previously discussed one. In this case one starts with a time dependent Hamiltonian of the theory and using the well known, Schwinger Keldysh formalism, which is a general path integral framework at finite temperature for the study of the time evolution of a quantum mechanical system which is in the out-of-equilibrium state. At the early time scale once a small perturbation or a response is provided to a quantum system, then it is described by a out-of-equilibrium process within the framework of quantum statistical mechanics and the present formalism provides us the sufficient tools and techniques using which one can compute the expressions for the OTOCs. Not only that, the late time behaviour of such OTOC is described by a saturation behaviour for chaotic systems from which one can compute the various characteristic features of large time equilibrium behaviour from these OTOCs.

Formalism III:
The third approach is based on the circuit quantum complexity, [75][76][77][78][79][80][81] which is relatively a very new concept and physically defined as the minimum number of unitary operators, commonly known as quantum gates, that are specifically required to construct the desired target quantum state from a suitable reference quantum mechanical state. In amore generalized physical prescription, quantum mechanical complexity can serve as one of several strong diagnostics for probing the time disorder averaging phenomena of a quantum mechanical system or quantum randomness. The underlying physical concept of circuit complexity can essentially provide essential information regarding various aspects of quantum mechanical randomness, such as the concept of scrambling time, Lyapunov exponent 6 , etc, which are particularly the key features of the study of quantum mechanical chaos. One can further compare between the physical outcomes of the two strongest measures of quantum randomness, which are appearing from out-of-time-order correlators (OTOCs) and the quantum mechanical circuit complexity and comment on further that for a specific quantum system which one is capturing the more information regarding the description of quantum mechanical randomness.
In this paper we generalize the study of OTOCs for investigating the phenomena of quantum randomness in various supersymmetric integrable quantum mechanical models. The main motivation behind introducing the concept of supersymmetry lies in the well established fact that for any Supersymmetric quantum mechanical model, the original Hamiltonian is always associated with a partner Hamiltonian, which in general is widely different from its original counterpart in terms of eigenstates. Though not always but it is possible that the quantum mechanical model under consideration, attains vastly different properties in the context of quantum randomness, due to the introduction of supersymmetry within the framework of quantum mechanics. It is our expectation that this generalization of the study of all the classes of OTOCs would provide an understanding about the significant role that supersymmetry plays in modifying the randomness behaviour of the quantum mechanical models under investigation. Most importantly, from the present study it will be clear that the additional inclusion of symmetries in the form of supersymmetry will at all effect, if so then how much it will effect the time disorder averaging phenomena for a canonical and micro-canonical statistical ensemble studied within the framework of out-of-equilibrium quantum statistical mechanics. In the eigenstate representation of the OTOCs we have actually studied three types of OTOC in this paper, out of that one of them is commonly studied in this literature mostly used to explain the phenomena of 6 The associated time scale when the quantum circuit complexity starts to grow is usually identified to be the scrambling time scale and in the representative plot with respect to the time scale particularly the magnitude of the slope of the linear portion of the curve is physically interpreted as the Lyapunov exponent for the specific systems where the general quantum randomness or the time disorder averaging phenomena is described by quantum chaos. quantum mechanical randomness (not only chaotic behaviour, but also a general feature including any types of non-chaotic behaviour), and the other two OTOCs that we have included in this paper might not be completely independent physical information of each other, but extremely important to capture the complete quantum mechanical effect in the quantum mechanical correlation functions to give a more general physical interpretation and a complete and detailed description of time disorder averaging phenomena for various quantum statistical ensembles within the framework of supersymmetric quantum mechanics. Now if we can't able to find out the eigenstate representation of a given time independent Hamiltonians within the framework of supersymmetry then the prescribed methodology for computing the general class of OTOCs in terms of the simplest eigenstate representation will not work. This can happen for the physical systems which is actually described by the time dependent Hamiltonians. In that case to compute all of these previously mentioned general class of OTOCs one need to use the quantum mechanical path integral generalization of the present framework, which is commonly described by the, Schwinger Keldysh formalism, at finite temperature within the framework of supersymmetric quantum mechanics. The good news is physical outcomes of such a generalization also have not studied yet, which we have future plan to look into in detail and we are also hopeful to get non-trivial and better understanding of various supersymmetric quantum mechanical models out of these computations.
Organization of the Paper is as follows- • In Sec. 3 we provide a brief review of the concept of Supersymmetric Quantum Mechanics.
• In Sec. 4, we explain how the phenomenon of Quantum Randomness can be diagnosed through the out of time ordered correlators.
• In Sec. 5, we provide a model independent eigenstate representation of the 2 and the 4-pt correlators of all the three kinds defined equally well for any QM model with well defined eigenstates.
• In Sec. 6 we explicitly calculate the correlators for the supersymmetric Harmonic Oscillator.
• In Sec. 7 we provide the numerical calculations of the correlators for the Supersymmetric 1D infinite potnetial well.
• In Sec. 8 we discuss the semiclassical analogue results for the two Supersymmetric QM model.
• Finally we conclude with the most important observations from our analysis of the considered Supersymmetric QM models.
Energy difference between the nth and mth energy eigenstate.
Un-normalized 2-pt correlator of third kind.
Un-normalized 4-pt canonical correlator of first kind.
Un-normalized 4-pt canonical correlator of second kind.
Un-normalized 4-pt canonical correlator of third kind.
Normalized 2-pt correlator of third kind.
Normalized 4-pt correlator of first kind.
Normalized 4-pt correlator of second kind.
Normalized 4-pt correlator of third kind.

A short review of Supersymmetric Quantum Mechanics
The theory of Supersymmetric Quantum Mechanics [83][84][85][86][87][88] relates quantum eigenfunctions and the corresponding eigenvalues between two partner Hamiltonians through an intertwining relationship using the so called charge operators. It generally uses the technique of factorizing the Hamiltonian in terms of the intertwining operators and hence to determine the superpotential using the well known Riccati equation. The idea of factorizing generally allows one to express the superpotential in terms of the ground state wave-function of the original Hamiltonian of the quantum mechanical model under consideration. The process of factorization is mathematically described by the following equation: where A and A † are the intertwining operators which are defined from the superpotential as: Within the framework of supersymmetry, the ground state energy is usually taken to be zero, which is well justified because it is only the relative energy difference that matters. For a zero energy ground state, the Schrödinger equation can be written as: Substituting the expressions for A and A † in the above equation it is not very hard to derive the Riccati equation, which further gives a way of writing the potential in terms of the superpotential given by the following expression: where corresponds to d/dx in the above equation. It is often an overrated fact that one needs to know the form of the potential guiding the Hamiltonian to have an idea about the wave functions of the quantum mechanical system and the fact that the knowledge of the ground state wave function allows one to exactly know the potential associated with the system is overlooked. However, in Supersymmetric quantum mechanics one generally utilizes this unappreciated fact to construct the potential from the known ground state wave function with zero modes: The knowledge of superpotential allows one to determine the supersymmetric partner potential via the following equation: It is easy to show that the knowledge of the eigen functions of H 1 can be used to derive the eigen functions of H 2 using the A operator and the eigen functions of H 1 from that of H 2 using the A † operator. The role of the operators A and A † , apart from the conversion of an eigenfunction of the original Hamiltonian to that of its partner Hamiltonian with the same energy, it also destroys or creates one node in the eigenfunction. This justifies the absence of the zero energy or the ground energy state of the partner Hamiltonian. One can put this argument simply by stating that the operator A converts an energy state of the original Hamiltonian into a lower energy state of the partner Hamiltonian keeping the energy value of the state constant. A † on the other hand does the opposite conversion i.e takes an energy state of the partner Hamiltonian and converts it into an higher energy state of the original Hamiltonian keeping the energy eigenvalue fixed.
A supersymmetric quantum mechanical model is generally described by a Hamiltonian having the following form: In general, for such a Hamiltonian a quantum state is represented by: where, |ψ (1) and |ψ (2) are the wave functions of the original and the partner Hamiltonian respectively. The prime objective of this paper is to provide an eigenstate representation of desired OTOCs, that we have already defined in the introduction, using which study the various well studied quantum mechanical models in the context of supersymmetry to study the general aspects of time disorder averaging phenomena. With this aim one would generally look for an eigenstate of the Hamiltonian under inspection. Remembering the relation between the energy eigenvalues and eigenstates of the original and the partner Hamiltonian it can be easily verified that the wave function given by Eq. (3.9) is not be an eigenstate of the Hamiltonian having structure as given in Eq. (3.8). We therefore take the wave function of the Hamiltonian to be of the form which indeed represents a normalized eigenfunction of the Hamiltonian of the supersymmetric quantum mechanical systems considered in this work.
mechanical system under study. Thermal average over the canonical statistical ensemble of a quantum operator is a very powerful technique using which one can explicitly study the time dependent exponential growth (chaotic) or some other time dependent non-chaotic random behaviour of an operator for a quantum system that is in out of equilibrium after giving an external response.For a very longer time it was not very clear that how actually one can quantify these quantum correlation functions within the framework of out-of-equilibrium quantum statistical physics. Following the previous set of works the present work helps us to quantify as well as to physically understand the impact of them in the present context. In this paper we are actually interested in three specific kinds of OTOCs which are described by six set of correlators which are given by: where x is the quantum position operator and p is the associated canonically conjugate momentum and most importantly both the quantum operators are defined at different time scales, which is one of the prime requirement to study the effect of time disorder averaging phenomena through the above set of OTOCs. Also, it is important to note that, the symbol · · · actually represents the thermal average of a time dependent quantum operator over a canonical ensemble within the framework of quantum mechanics, which is technically defined as: where the partition function Z and thermal density matrix operator ρ β in terms of a quantum system Hamiltonian, H, are already defined in the introduction of this paper.
Since we are dealing with quantum mechanical operators using which we are trying to understand the impact of random features in the quantum regime, it is quite expected to start with the fact that thermal one point function of the position operator x and momentum operator p defined at a specific time scale are zero, which can be technically demonstrated as, where t is the associated time scale on which both of the quantum operators are evolving. For this specific reason, the explicit study and the computation of these one point functions are not at all important in the present context of discussion. On the other hand, due to the time translational symmetry in these thermal correlators, which is actually described by the well known, Kubo Martin Schwinger condition, all the odd point OTOCs appearing in the present context will be trivially zero and for that reason not the object of interest in the present context of study. This can be further technically demonstrated as, Odd x correlator : x(t 1 )x(t 2 )x(t 3 ) · · · odd number of terms = 0, (4.3) Odd p correlator : · · · odd number of terms = 0, (4.4) One p & even x correlator : · · · odd number of terms = 0, (4.5) One x & even p correlator : x(t 1 )p(t 2 )p(t 3 ) · · · odd number of terms = 0. (4.6) This implies that we are only left with all even order OTOCs out of which in this paper we are explicitly computing the physical outcomes from the sex set of time dependent correlators, describe by the previously mentioned two point and the four point all possible OTOCs. Among these correlators two are made of different operators which will show the perturbation of one operator measured at one time scale to the other measured at a different time and vice versa. The other four operators are the new sets of two and four point correlators which are defined in terms of the same quantum operators, which are basically capturing the quantum effect of the self-correlation of one operator on itself having any arbitrary time dependent profile in general to describe the phenomena of time disordering within the framework of quantum mechanics. In specific cases it may happen that these newly defined operators show exponential growth or some other kind of growth which is periodic or aperiodic in the corresponding associated time scales of the quantum mechanical operator on which we are interested in the present context. In a more general context one can see by studying different kinds of physical systems available in the literature where the two-point self-correlation will decay exponentially with an associated time scale of t d , widely known as "dissipation time scale". This particular time scale is playing the role of transition scale in the present context after reaching that the four-point correlators can be factorised into the product of two two-point correlation functions representing disconnected diagrams within the framework of quantum field theory. Also, it is important to note that, here at the "dissipation time scale" all other terms are exponentially suppressed by the factor e −t/t d , which will completely disappear from the factorized version of the four point correlators on which we are interested in this paper in the large time limit given by, t t d with t → ∞. In usual prescription the "dissipation time scale", t d is identified with the inverse temperature β i.e. t d ∼ β = 1/T , where we use Boltzmann constant k B = 1 and T physically represents the equilibrium temperature of the quantum statistical ensemble on which we are interested in. So one can translate the large time limit in terms of the associated equilibrium temperature as, tT 1, which is not obviously true for zero temperature case, but can be justifiable in any (small or large) temperature of the system under consideration. For example, we now look into a specific four point thermal correlation function at the vicinity of the previously mentioned dissipation time scale, around which one can factorize it in the following specific form: Equal time 2−point correlator Equal time 2−point correlator where in this above mentioned factorization O(e −t/t d ) are the higher order correction terms which are actually sub-dominant at the vicinity of the dissipation time scale. One more thing we can observe from the above mentioned factorization is that, the individual contributions of the two point contributions don't mix up the time scales and for this reason they can be written in terms of the product of two equal time two point correlators. Also for this particular example when we are thinking of doing the computation with two different operators, after doing the factorization one can easily observe that the two different operators don't mix with each other at the level of two-point correlators. So, we are mainly interested in the terms that can not be written in time ordered form and those terms provide more insight to the randomness present in a system. Now we would like to go one step further and normalize the OTOC, which is basically related to this factorization process of the four-point correlators in terms disconnected equal time two point contributions. The process of normalization actually helps us to reduce the unwanted fluctuations from the computed OTOCs which further allows us to give a clearer picture of the time dependent behaviour of the desired OTOCs in which we are interested in this paper. We do not normalize the two-point correlation functions as they are the main building blocks of our computation of OTOCs and hence we only normalize the four-point OTOC as given by the following expression. We can write − [x(t 1 ), p(t 2 )] 2 in the following simplified mathematical form: After imposiing the previously mentioned constrained obtained at the vicinity of the dissipation time scale one can further simplify the expression for the mentioned correlator as given by the following expression: Now we normalize this above mentioned quantity using two different equal time two-point correlators which we have obtained from the factorization in terms of the disconnected pieces. Consequently, we get the following simplified form of the normalized OTOC: The first expression in Eq. (4.10) is a universal contribution, which will always appear for the quantum systems where the previously mentioned factorization process works at the , p(t 2 )] 2 , and by following the same logical argument at the vicinity of the dissipation time scale we can normalize them as well and write them in the following simplified mathematical forms: and Now, since these OTOCs acts as a theoretical probe to know about the generic chaotic or non-chaotic time dependent behaviour of quantum system, it has to satisfy the following constraints in terms of the four point correlations, which survived in the vicinity of the previously introduced dissipation time scale, as given by: The above mentioned expression captures all the possibilities which one can observe in different quantum mechanical systems available in nature. Here A 12 , A 12 and A 12 are the quantum mechanical model dependent pre-factors which shows exponential growth (chaotic behaviour) in the four point correlator with respect to the time scales associated with the system under study. On the other hand, B 12 and B (3) 12 are the quantum mechanical model dependent pre-factors which shows any type of time dependent fluctuations (non-chaotic behaviour). Also, for the general prescription the quantum Lyapunov exponents, λ 1 , λ 2 and λ 3 are not same for which the MSS bound on quantum chaos from these three cases are also not the same. Consequently, the equilibrium saturation temperatures at the late time scales from these three OTOCs also differ from each other i.e. β 1 = β 2 = β 3 . Also it is important to point out that, the mathematical structure of the time dependent functions f 1 (t 1 , t 2 ), f 2 (t 1 , t 2 ) and f 3 (t 1 , t 2 ) are also different for general quantum mechanical set ups. When we are considering the quantum mechanical models are described by time dependent Hamiltonian in that case these expectations and all sorts of predictions work very well. However, when we are thinking about particularly quantum mechanical models which are described by the time independent Hamiltonian and the eigenstate representation of OTOCs, in that situations one might have further simplifications. There might be possibility to have underlying connection between the two functions f 2 (t 1 , t 2 ) and f 3 (t 1 , t 2 ) in the eigenstate representation of the OTOCs in the non-chaotic case and for this reason they might not be capturing completely independent information of the time disorder averaging. On the other hand, in the chaotic case if three of the OTOCs independently show exponential growth in the time scale the quantum Lyapunov exponents and the related equilibrium saturation temperatures are not at all same even in the eigenstate representation. But if the first OTOC is showing the chaotic behaviour and other two are not, in that case the previous connection between the two functions f 2 (t 1 , t 2 ) and f 3 (t 1 , t 2 ) holds good in the eigenstate representation.
In figure, 3 and 4, we have presented the diagrammatic representations of all types of OTOCs in which we are interested in this paper. Particularly in figure 3, we have explicitly depicted the possible two point OTOCs. Here we have three possibilities, which are given by, ] . Since each OTOC are made up of commutator bracket in the quantum mechanical description, for each case we have two different contributions having overall opposite signatures. Further to draw the representative diagrams we need to consider the flow of time scale from t 1 to t 2 or t 2 to t 1 . In all the representative diagrams, the two vertical solid thick line corresponds to the specific time slice having time, t = t 1 and t = t 2 respectively. It is understandable from the mathematical structure of the mentioned two point OTOCs that since the correlator involve only two time scales , that is why drawing two vertical parallel lines are physically justifiable in the present context. Because of the previously mentioned flow of time scale from t 1 to t 2 or t 2 to t 1 for each two point OTOCs we have two possible diagrams. So as a whole for the two-point OTOCs one can draw six possibilities. Also by studying each of the diagrammatic representations we can also observe that each of the contributions of the two point functions are represented by separate lines with representative arrows on that which is completely depend on the structure of the individual two point correlators. Two differentiate between these two contributions we have used red dotted line and blue solid line in the representative diagrams. Further, in figure 4, we have explicitly shown the possible four point OTOCs in the present context of discussion. Here we can draw three possible representative diagrams, which are coming from the three possible OTOCs, as given by, Here as we can see that the each of the OTOCs are made up of commutator bracket squared contributions in the quantum description, for each case we have four contributions if we expand them. Out of these four four point thermal correlators two of them having positive signature and other two of them overall negative signature in the front of each contributions. Just like the two point OTOCs here also we need to consider the flow of time scale from t 1 to t 2 or from t 2 to t 1 . It is understandable from the mathematical structure of the mentioned four point OTOCs that since the correlator involve only two time scales instead of four different time scales , that is why drawing two vertical parallel lines to represent the time slice at t = t 1 and t = t 2 are physically justifiable in the present context. Now, because of the time scale flow each of the four point OTOCs have two contributions in the diagrammatic representation. Also for a given time flow we have four possible diagrams which we shown in a single diagram for the sake of simplicity. So as whole consider both the possibilities of the time flow we have cumulatively 24 diagrams from the four point OTOCs. Like the previous case, here also to differentiate between each of the individual terms for a given OTOC with a specified time flow we have used red dotted line, blue dotted line and red thick line respectively in the representative diagrammatic representations. Last but not the least, here it is important to point that these set of diagrams are the simplest version of the well known Feynman diagrams as appearing in the context of quantum field theory literature. However, within the present framework we don't have exactly similar Feynman diagrams, but to understand the structure of the previously mentioned OTOCs the present version of the diagrammatic representations play a very crucial role. All of the arrows appearing in all the diagrams representing an underlying time disordering upon which we have to take finally average in the present context.

Eigenstate Representation of thermal OTOCs in Supersymmetric Quantum Mechanics
In this section our prime objective is to study various thermal correlators, OTOCs and see how they can be expressed in a model independent manner in the framework of quantum mechanics. We also demarcate clearly the effect of Supersymmetry and how it modifies the functional form of the correlators at the end.
To perform this explicit computation we will follow the prescription outlined in the ref. [70] to compute all the thermal OTOCs we have mentioned in the previous section of this paper 7 . According to this prescription for any time-independent Hamiltonian we define the expectation value of quantum mechanical operators as thermal expectation value in the present context, which one can easily apply for a canonical quantum statistical ensemble. Let's say, a quantum mechanical time dependent operator A(t) is defined on a specific Hilbert Space H with an associated Hamiltonian H (which is time independent obviously) and the eigenvalues (eigen energy spectrum) of H corresponding to a infinite tower of eigenkets |Ψ n ∀ n = 0, 1, · · · , ∞, where the corresponding energy levels are characterized by the index n. This energy eigen spectrum is represented by E n ∀ n = 0, 1, · · · , ∞. Then for a canonical quantum statistical ensemble the thermal expectation value of the quantum mechanical time dependent operator A(t) at inverse temperature β, considering the Boltzmann constant k B = 1, is defined as: where the thermal partition function is represented by Z such that Z = Tr e −βH . Here our job is to represent this mathematical trace operation in terms of the sum over all possible eigenstates starting from the ground state (n = 0). Once we can able to express this operation clearly then the rest of the story is very usual and following this one can easily compute all of the mentioned OTOCs in the eigen state representation. It pays to use the thermal representation for expectation values of quantum mechanical operators because there are many physical models and physically relevant toy models which have wellstudied structure of time independent Hamiltonian and hence we can utilize this property to give an eigenstate representation to the correlators as presented elaborately in this work. Consequently, the analysis presented in this paper is physically justifiable, applicable and believable, as such, to quantum mechanical models with well defined Hamiltonian. In the eigenstate representation the thermal expectation value of a quantum mechanical time dependent operator A(t) can be written as : Here we work in the Heisenberg Picture where the operators evolve with time as given by, where O H (t) represents operator O H at time scale t in Heisenberg Picture and O S (0) represents its Schrödinger Picture representation at all time scales since operators do not evolve with time in Schrödinger Picture. The later one is denoted as O S to simplify the notation and make it easier to read.
What we learn here?
To demonstrate the computation further let us consider a general time independent Hamiltonian of the following form: where we have used the fact that the mass of all N number of particles are same and given by, m i = m = 1/2 to make the further computation simpler. Now here it is very easy to prove the following relation (see [70]), which relates the quantum mechanical momentum operator matrix elements with that of the position and the energy operator matrix elements by the following expression: Furthermore, the role Supersymmetry plays in modification of non-supersymmetric quantum mechanical observables can be seen vis-a-vis the matrix elements of observables. Let's consider a simple matrix element for the position operator x. In the following chart we show exactly show how the matrix elements in non-supersymmetric quantum mechanical theories and supersymmetric quantum mechanical theories are connected: The modifications due to supersymmetry can clearly be traced back to the fact that we could formulate a wave function of total Hamiltonian as given in Eq. (3.10) from the two partner Hamiltonians H 1 and H 2 , as appearing from the bosonic and the fermion sectors in supersymmetry. In fact, its the factorization property of H 1 and H 2 which generalizes to a more powerful setting in terms of supersymmetric generalized description of the theory under consideration. One might consider the power of the supersymmetric description in two fold ways: 1. First, one can merely consider it as a tool to solve non-trivial potentials by means of solutions of their partner ones, provided that the partner ones are easily solvable.
2. Secondly, one can consider it as a more unifying description of a more beautiful theory based on the principles of symmetries of nature. It is this second philosophy which has been considered in this work.
We will consider the following six correlators in this work as stated below : Correlation Functions We will also consider normalized 4-pt Correlation Functions and the discussion pertaining to those is provided in the previous section:

Partition Function from Supersymmetric Quantum Mechanics
In the context of Supersymmetric quantum mechanics the partition function, Z, can be can be expressed in terms of the eigenstate by the following expression: Here the underlying relation between supersymmetric total system and the component (or partner) systems, as discussed in Sec. 3, is established through the following relations : So, the supersymmetric quantum partition function in the eigenstate representation can be explicitly written as : Ground State Contribution As shown in Sec. 3 in supersymmetric quantum mechanical theories we have the following constraint: with the requirement that the ground state is eigenstate of H 1 only and all other states are doubly degenerate. It is also important to note that the ground state always has zero energy eigenvalue. It is to be noted that E m refers to the energy eigenvalue of the total Hamiltonian H SUSY . Using this sets of requirements, the above mentioned expression for the quantum partition function within the framework of supersymmetry can be further simplified as: which is different compared to usual quantum mechanics result without supersymmetry as in the present context ground state energy eigen value E 0 = 0 which is not in general zero for the other case. So this implies that if we separately write down the contribution from the ground state and from all other excited states, instead of writing sum over all eigenstates together, then one can clearly visualize the difference between the results obtained for quantum partition function in the eigenstate representation in both the cases. In the next section, we will provide the summary of all the obtained general model independent results in which we will implement the above mentioned fact to explicitly show that our obtained results for all types of OTOCs in the framework of supersymmetric quantum mechanics are different compared to the results obtained from usual quantum mechanics without supersymmetry.
The first two point OTOC is given by the thermal average of the operator −[x(t 1 ), p(t 2 )], which is described in the eigenstate representation as: Using Eq. (5.2) for Heisenberg representation for x(t 1 ) and p(t 2 ) and inserting the identities between the operators this two-point correlator can be written in terms of the microcanonical correlator, which shows the temperature independent behaviour of the system: where the micro-canonical two point OTOC is defined as: Expanding the commutator in the definition of the micro-canonical correlator, and using the above relation with appropriate insertion of identities between the operators, it can be shown that the eigenstate representation of the micro-canonical correlator y m (t 1 , t 2 ) takes the following form.
where we define, E mk/m,k = E m − E k and x mk/m,k = Ψ m |x|Ψ k .
Substituting the above mentioned expression for the micro-canonical OTOC in the definition of the canonical correlator, Y (1) (t 1 , t 2 ) we get the following expression : (5.20) The second two point OTOC is given by the thermal average of the operator −[x(t 1 ), x(t 2 )], which is described in the eigenstate representation as: Using Eq. (5.2) for Heisenberg representation for x(t 1 ) and x(t 2 ) the two point correlator in terms of the temperature independent micro-canonical correlator can be written as: where the micro-canonical two point OTOC is defined as: Expanding the commutator in the definition of the micro-canonical correlator and inserting identity in between the operators, it can be shown that the eigenstate representation for micro-canonical correlator y m (t 1 , t 2 ) takes the following form: Substituting Eq. (5.23) in Eq. (5.21), the eigenstate representation for the canonical correlator Y (2) (t 1 , t 2 ) takes the form : The third two point OTOC is given by the thermal average of the operator −[p(t 1 ), p(t 2 )], which is described in the eigenstate representation as: Proceeding as previously mentioned we obtain the following expression for the OTOC which can be expressed in terms of the temperature independent micro-canonical correlator as: where the micro-canonical two point OTOC is defined as: Expanding the commutator in the definition of the micro-canonical correlator and inserting identity in between the operators, it can be shown that the eigenstate representation for micro-canonical correlator y m (t 1 , t 2 ) is given by the following expression: Substituting the expression for the micro-canonical correlator obtained from Eq. (5.27) in Eq. (5.25), the eigenstate representation for the canonical correlator Y (3) (t 1 , t 2 ) can be written as : (5.28)

5.5
Representation of four point OTOC : In this section, we provide the eigenstate representation of the four-point OTOCs that are mainly used for studying the quantum mechanical analogue of phenomena of time disorder averaging , related randomness in quantum regime. Similar to the two-point correlators we define three kinds of four-point correlators. Generally in the literature people study the four-point correlator C (1) (t 1 , t 2 ) (as per our notation) which is the thermal expectation value of the square of the commutator made up of the dynamical operators of different kinds at different time scales. The four-point correlator, C (1) (t 1 , t 2 ) is proposed as a quantifier of quantum chaos but chaos is a specific kind of randomness. Hence, it is important to consider the correlators constructed from similar operators at different times for having a complete understanding of the underlying phenomenon of randomness. To understand this time disordering phenomena explicitly we have also introduced two more new correlators, C (2) (t 1 , t 2 ) and C (3) (t 1 , t 2 ), described in detail in the later half of this paper.
The first four point OTOC is given by the thermal average of the operator −[x(t 1 ), p(t 2 )] 2 , which is described in the eigenstate representation as: is the micro-canonical correlator and is responsible for the temperature independent behaviour of the correlator. The temperature dependence in the canonical correlator is actually appearing due to the exponential thermal Botzmann factor in the eigenstate representation of the OTOC. Once we will take the sum over all possible eigen states (finite or infinite in number) we get a cumulative dependence on time, temperature and energy eigenstates. In the present context the micro-canonical four-point correlator for Eq. (5.29) is given by the following simplified expression: where we define a time dependent matrix element, b (1) Using the Heisenberg picture evolution equation for an operator and inserting identity between the operators after expanding the commutator, it is not hard to verify that b (1) mk can be written as: Substituting the above expression of b (1) * mk (t 1 , t 2 ) in Eq. (5.30), it can be shown that after simplification the eigenstate representation of the micro-canonical correlator for Eq. (5.29) is given by the following expression: So, the eigenstate representation of canonical correlator from Eq. (5.29) using Eq. (5.31) is given as : The normalized first four-point OTOC is from the un-normalized one expressed in Eq. (5.12) can be expressed as: We have already calculated the numerator C (1) (t 1 , t 2 ) as given in Eq. (5.32). So here our only job is compute the two sets of equal time thermal correlators in the eigenstate representation. Since we know the formalism very well that we have developed in this paper, computing these disconnected pieces of equal time correlators which are extremely significant in the large time dissipation limit are not at all complicated. Now we are going to show that how one can compute these contributions.
To serve this purpose we need to compute the following expressions in the eigenstate representation: Consequently, the desired canonical equal time thermal correlators can be computed as: Finally, the normalized first four-point OTOC can be expressed by the following simplified expression in the eigenstate representation, as given by: (5.37)

5.6
Representation of four point OTOC : The second four point OTOC is given by the thermal average of the operator −[x(t 1 ), x(t 2 )] 2 , which is described in the eigenstate representation as: is the micro-canonical correlator and is responsible for the temperature independent behaviour of the correlator. The temperature dependence in the canonical correlator is actually appearing due to the exponential thermal Botzmann factor in the eigenstate representation of the OTOC. Once we will take the sum over all possible eigen states (finite or infinite in number) we get a cumulative dependence on time, temperature and energy eigenstates. In the present context the micro-canonical four-point correlator for Eq. (5.38) is given by the following simplified expression: where we define a time dependent matrix element, b mk (t 1 , t 2 ), which is given by: Using the Heisenberg picture for the evolution of operators and inserting identity between the operators, it can be shown that b (2) mk (t 1 , t 2 ) can be written as: Substituting Eq. (5.40) in Eq. (5.39) and simplifying the eigenstate representation for the temperature independent micro-canonical correlator can be expressed by the following simplified expression: The canonical or the temperature dependent correlator can be calculated by substituting the expression of Eq. (5.38) in Eq. (5.41) and it after applying some algebraic manipulation it can be shown that the canonical correlator can be expressed by the following two expressions : The normalized first four-point OTOC is from the un-normalized one expressed in Eq. (5.13) can be expressed as: We have already calculated the numerator C (2) (t 1 , t 2 ) as given in Eq. (5.42). So here our only job is compute the two sets of equal time thermal correlators in the eigenstate representation. Since we know the formalism very well that we have developed in this paper, computing these disconnected pieces of equal time correlators which are extremely significant in the large time dissipation limit are not at all complicated. Now we are going to show that how one can compute these contributions.
To serve this purpose we need to compute the following expressions in the eigenstate representation: Consequently, the desired canonical equal time thermal correlators can be computed as: Finally, the normalized first four-point OTOC can be expressed by the following simplified expression in the eigenstate representation, as given by: (5.45)

5.7
Representation of four point OTOC : The third four point OTOC is given by the thermal average of the operator −[p(t 1 ), p(t 2 )] 2 , which is described in the eigenstate representation as: is the micro-canonical correlator and is responsible for the temperature independent behaviour of the correlator. The temperature dependence in the canonical correlator is actually appearing due to the exponential thermal Botzmann factor in the eigenstate representation of the OTOC. Once we will take the sum over all possible eigen states (finite or infinite in number) we get a cumulative dependence on time, temperature and energy eigenstates. In the present context the micro-canonical four-point correlator for Eq. (5.46) is given by the following simplified expression: where we define a time dependent matrix element, b mk (t 1 , t 2 ), which is given by: Using the Heisenberg picture for the evolution of operators and inserting identity between the operators, it can be shown that b mk (t 1 , t 2 ) can be simplified into the following form: Substituting the expression of b The eigenstate representation of the canonical correlator from Eq. (5.46) using Eq. (5.49) after simplification can be expressed by the following two equations : The normalized first four-point OTOC is from the un-normalized one expressed in Eq. (5.14) can be expressed as: We have already calculated the numerator C (3) (t 1 , t 2 ) as given in Eq. (5.50). So here our only job is compute the two sets of equal time thermal correlators in the eigenstate representation. Since we know the formalism very well that we have developed in this paper, computing these disconnected pieces of equal time correlators which are extremely significant in the large time dissipation limit are not at all complicated. Now we are going to show that how one can compute these contributions.
To serve this purpose we need to compute the following expressions in the eigenstate representation: Consequently, the desired canonical equal time thermal correlators can be computed as: Finally, the normalized first four-point OTOC can be expressed by the following simplified expression in the eigenstate representation, as given by: (5.53)

Eigenstate Representation for Micro-Canonical Correlators
Eigenstate Representation for Canonical Correlators without normalization Eigenstate Representation for Canonical Correlators with normalization We consider a system described by Hamiltonian H 1 with the harmonic oscillator potential as given by: where ω is the natural frequency of the oscillator and we have assumed that the mass, m = 1, for the sake of algebraic simplification. Now, we represent Ψ (1) n (x) and E (1) n as the eigenstates and eigenvalues corresponding to Hamiltonian H 1 which can be obtained by solving the Schrödinger equation as: where H m (y) appearing in the eigenfunctions are the well-known Hermite polynomials of order m. As discussed in Sec. 3 the superpotential can be calculated once the ground state wavefunction is known. Hence, using the ground state wavefunction obtained by substituting n = 0 in ??, the superpotential can be calculated using Eq. (3.5) and for the case of supersymmetric one dimensional harmonic oscillator it is obtained to be the following: Thus, we can construct the partner Hamiltonian which is of the following form: where V 2 is the partner potential and can be computed using Eq. (3.6). For the supersymmetric one dimensional harmonic oscillator we get the following expression for the partner potential: More precisely, in the context of supersymmetric harmonic oscillator, the original Hamiltonian is usually identified as the Bosonic Hamiltonian and its associated super-partner is identified as the Fermionic Hamiltonian, which are constructed by subtracting off the ground state energy as: Bosonic Hamiltonian : Fermionic Hamiltonian : Now since this is just a constant shift in the energy, the corresponding eigenstates of the partner Hamiltonians remain the same, and can be represented as: Bosonic eigenstate :

9)
Fermionic eigenstate : However, in this case the energy eigenvalues get shifted by an amount E (1) 0 and can be described as:

11)
Fermionic eigenspectrum : where we have used the property of supersymmetric quantum mechanical theories:

Partition Function
The partition function in the eigenstate representation is defined by the following expression: In the framework of supersymmetric quantum mechanics as already discussed the ground state is only bosonic. The structure of the normalized eigenfunctions are also different for the ground state and the higher energy state. Hence, the partition function can be written in terms of the ground state contribution and in terms of the contribution coming from the excited state as: Excited State Contribution (6.15) Further using the fact that the bosonic and the fermionic Hamiltonians n th energy eigenfunction are exactly identical, the partition function for the supersymmetric one dimensional harmonic oscillator can be computed as: 6.3 Computation of two point OTOCs We compute the 2-pt correlators for the Supersymmetric Harmonic oscillator. We begin with the correlator constructed from two different dynamical operators which is defined as: where y 0 (t 1 , t 2 ) and y (1) n (t 1 , t 2 ) represent the micro-canonical correlators for the ground state and the excited states and the overall temperature dependent normalization function is appearing from the expression for the inverse of the thermal partition function for canonical quantum statistical ensemble. The purpose of taking the ground state and excited states separately is due to the fact that in supersymmetric quantum mechanics the ground state is always bosonic having zero energy eigenvalue, whereas for the excited states there are contributions from the bosonic as well as the fermionic states. It is easy to see that the ground state contribution can be calculated by expanding the commutator and inserting the completeness relation, k |Ψ k Ψ k | = 1. However, here it is important to note that there is no contribution appearing from the k = 0 term as for the ground state we have, Ψ 0 |x(t 1 )|Ψ 0 = 0. Hence, from the ground state we finally get the following time dependent non-trivial contribution, which is given by: Now further to determine the contribution from the excited states we need to expand the commutator and insert the completeness relation between the desired time dependent quantum operators which gives the following expression: Therefore, the micro-canonical correlators for the ground state and the excited states are given by the following temperature independent simplified expressions: cos ω(t 1 − t 2 ) 1 + δ n1 + n(n + 1) − n(n − 1) . (6.21) Further, using Eq. (6.20) and Eq. (6.21) the canonical two point thermal correlator can be obtained by the following expression: It can be seen that the obtained OTOC that it is just a periodic function indicating the absence of chaos. However, it is quite interesting that the obtained result has an explicit temperature dependence which mainly comes from the contribution of the ground state.
Though the two point OTOC is not of much significance. It stills give us a pretty good idea that at the level two point correlation function is that in the context of supersymmetric one dimensional harmonic oscillator is non-chaotic in nature.
Micro-Canonical 2-pt Correlator for Ground State : Micro-Canonical 2-pt Correlator for Excited States : Canonical 2-pt Correlator : As discussed earlier the correlator Y (1) (t 1 , t 2 ) alone does not provide the complete knowledge about the quantum randomness of the system. Hence, we need to calculate two other correlators constructed from similar operators at different time scales. The other two correlators will be : − [x(t 1 ), x(t 2 )] β and − [p(t 1 ), p(t 2 )] β which we will compute in the next part of the paper.
The second correlator is constructed from the commutator of the two position operators defined at different time scales and can be represented as: The reason for separately considering the ground state from the other states has already been discussed in the calculation of Y (1) (t 1 , t 2 ). Following a similar procedure we get: sin ω(t 1 − t 2 ) 1 + δ n1 + n(n + 1) − n(n − 1) (6.28) Using Eq. (6.27) and Eq. (6.28) the canonical correlator can be obtained as: e −βω(n− 1 2 ) 1 + n(n + 1) − n(n − 1) (6.29) Micro-Canonical 2-pt Correlator for Ground State :

Canonical 2-pt Correlator
: The third correlator is constructed from the commutator of the two momentum operators defined at different time scales and can be represented as: The reason for separately considering the ground state from the other states has already been discussed in the calculation of Y (1) (t 1 , t 2 ). Following a similar procedure we get the following expressions for the micro-canonical correlators for the ground state and excited states as given by: sin ω(t 1 − t 2 ) 1 + δ n1 + n(n + 1) − n(n − 1) (6.35) Using Eq. (6.34) and Eq. (6.35) the canonical correlator can be obtained as: e −βω(n− 1 2 ) 1 + n(n + 1) − n(n − 1) .

Computation of un-normalized four point OTOCs
Though the building block two point OTOCs gives us some basic insight about the nature of randomness present in the quantum system, it is actually the four point OTOCs that are considered to be the prime observables for determining the degree and the nature of randomness present in the quantum mechanical systems. Similar to the two point case we consider here different types of OTOCs for the four point case from all possible three combinations of the dynamical operators characterising the quantum mechanical system under consideration, here it is one dimensional supersymmetric harmonic oscillator.
The corresponding four point desired OTOC of first kind is defined by the following expression: where the ground and the excited contribution of the micro-canonical part of the four point temperature independent OTOCs are defined by the following expressions: and m0 (t 1 , t 2 ) and b (1) 00 (t 1 , t 2 ) used in the above expression are given by: Now an extremely important fact to keep in mind while performing this calculation is that the ground state is always bosonic, this enforces considering the ground state separately each time an identity is inserted. On expanding the commutator and inserting the identity operator (I := k |Ψ k Ψ k |) and considering the k = 0 term separately from the k > 0 term above sets of equations can be written as: and a similar expression will be obtained for b (1) m0 (t 1 , t 2 ) and b (1) 00 (t 1 , t 2 ). Keeping all the above discussed facts it is not difficult to show that for supersymmetric one dimensional harmonic oscillator we get: cos ω(t 1 − t 2 ) δ n1 + 1 + n(n + 1) − n(n − 1) δ nm , Here one can explicitly show that for the supersymmetric one dimensional harmonic oscillator case the following contributions are trivially vanish: The above equations can be used to calculate the micro-canonical correlator which is the temperature independent part of the total OTOC. The temperature dependent part comes from the canonical part of the correlator which can be calculated by substituting Eq. (6.48) and Eq. (6.49) in Eq. (6.78) and for the supersymmetric one dimensional harmonic oscil-lator we get: Hence the temperature dependent canonical part of the four point OTOC of the first kind for a supersymmetric harmonic oscillator is given by substituting Eq. (6.52) and Eq. (6.53) in Eq. (6.40) to obtain the following simplified result: Canonical 4-pt Correlator : (6.57)

Computation of
The corresponding four point desired OTOC of second kind is defined by the following expression: where the ground and the excited contribution of the micro-canonical part of the four point temperature independent OTOCs are defined by the following expressions: and m0 (t 1 , t 2 ) and b (2) 00 (t 1 , t 2 ) used in the above expression are given by: Now an extremely important fact to keep in mind while performing this calculation is that the ground state is always bosonic, this enforces considering the ground state separately each time an identity is inserted. On expanding the commutator and inserting the identity operator (I := k |Ψ k Ψ k |) and considering the k = 0 term separately from the k > 0 term above sets of equations can be written as: and a similar expression will be obtained for b (2) m0 (t 1 , t 2 ) and b (2) 00 (t 1 , t 2 ). Keeping all the above discussed facts it is not difficult to show that for supersymmetric one dimensional harmonic oscillator we get: Here one can explicitly show that for the supersymmetric one dimensional harmonic oscillator case the following contributions are trivially vanish: (2) * n0 (t 1 , t 2 ) = 0 ∀ n > 0. (6.69) The above equations can be used to calculate the micro-canonical correlator which is the temperature independent part of the total OTOC and the corresponding ground and excited state contributions are explicitly given by the following expressions: The canonical part of the four point thermal OTOC of the second kind for a supersymmetric Harmonic Oscillator is given by substituting Eq. (6.70) and Eq. (6.71) in Eq. (6.58) to obtain the following simplified result: 2ω 2 e − βω 2 sinh βω 2 sin 2 ω(t 1 − t 2 ) × 1 + e −βω − n>0 e −βωn 1 + n(n + 1) − n(n − 1) e −βωn 1 + n(n + 1) − n(n − 1) 2 .

Micro-Canonical 4-pt Correlator for Excited States
: Canonical 4-pt Correlator : The corresponding four point desired OTOC of third kind is defined by the following expression: where the ground and the excited contribution of the micro-canonical part of the four point temperature independent OTOCs are defined by the following expressions: and m0 (t 1 , t 2 ) and b Now an extremely important fact to keep in mind while performing this calculation is that the ground state is always bosonic, this enforces considering the ground state separately each time an identity is inserted. On expanding the commutator and inserting the identity operator (I := k |Ψ k Ψ k |) and considering the k = 0 term separately from the k > 0 term above sets of equations can be written as: and a similar expression will be obtained for b m0 (t 1 , t 2 ) and b (3) 00 (t 1 , t 2 ). Keeping all the above discussed facts it is not difficult to show that for supersymmetric one dimensional harmonic oscillator we get: Using the above equations the ground state and the excited state contributions to the micro-canonical OTOC can be calculated as follows: Then the canonical part of the thermal four point OTOC of the third kind for a supersymmetric harmonic oscillator is given by substituting Eq. (6.87) and Eq. (6.88) in Eq. (6.76) to obtain the following simplified result: Micro-Canonical 4-pt Correlator for Excited States : Canonical 4-pt Correlator : (6.92)

Computation of normalized four point OTOCs
In this section our prime objective is to normalize all of the derived un-normalized three types of OTOCs with appropriate normalization factors, which we have already introduced in the earlier half of the paper for model independent eigenstate representation. In this section we will explicitly derive the expressions for the normalization factors for supersymmetric ine dimensional harmonic oscillator model in its eigenstate representation and derive all of the possible three types of OTOCs after normalization.
To normalize the obtained OTOC C (1) (t 1 , t 2 ) we need to compute the appropriate factors which we are going compute in this subsection.
First of all, we need to evaluate the following two point equal time thermal correlator, which is given by: e −βEn Ψ n |x(t 1 )x(t 1 )|Ψ n Excited State contribution (6.93) Inserting the completeness relation between the operators and using the Heisenberg picture equation for the evolution of an operator, the normalization factor involving the position operator can be written as: A similar calculation can carried out for computing the two point equal time thermal correlator involving the momentum operators, which is given by: Ground state contribution + n>0 ψ n |p(t 2 )p(t 2 )|ψ n Excited state contribution (6.95) One generally needs to consider the ground state separately from the other higher energy states due to the fact that in Supersymmetric QM, the ground state has contributions only from the original Hamiltonian. There is no contribution of the associated partner Hamiltonian in the ground state. Keeping this interesting fact in mind and adopting a similar approach as taken in the previous case, the normalization factor involving the thermal expectation value of the momentum operators are given by: e −βωn [2n + n(n + 1) + n(n − 1)] (6.96) Therefore, the normalization factor N 1 for the 4-point correlator C (1) (t 1 , t 2 ) is given by Consequently, the normalized OTOC of the first kind can be expressed as: : e −βωn 1 + n(n + 1) − n(n − 1) (6.98) To normalize the obtained OTOC C (2) (t 1 , t 2 ) we need to compute the appropriate factors which we are going compute in this subsection.
Similarly the normalization factor N 2 for the four point OTOC C (2) (t 1 , t 2 ) can be computed explicitly, which is given by the following simplified expression: Consequently, the normalized OTOC of the second kind can be expressed as: Normalized C (2) (t 1 , t 2 ) : To normalize the obtained OTOC C (3) (t 1 , t 2 ) we need to compute the appropriate factors which we are going compute in this subsection.

Model II: Supersymmetric One Dimensional Potential Well
The one dimensional infinite potential well is characterized by the following potential: The eigenfunctions and the corresponding energy eigenvalues associated with the Hamiltonian H 1 for this potential is a well known result [88] and is given by the following expressions: where we have replaced the energy level tagging from n to n + 1 compared to the regular expressions one encounters in typical textbooks. This does not change anything physically and with this convention n can take values from 0 instead of 1. Also, for the sake of simplicity we consider = L = 2m = 1. Furthermore, as we have previously mentioned we need the ground state energy to be zero, so that after subtracting off the ground state energy we get the following simplified results: ψ (1) n = √ 2 sin[(n + 1)πx] and E (1) n = n(n + 2) π 2 for n ∈ {0, 1, 2, ...}.

(7.2)
To obtain the partner potential associated with the original potential, the superpotential needs to be calculated which has been done in [88] and is given by the following expression: Once the superpotential of a supersymmetric quantum mechanical model is known it can be used to obtain the eigenspectrum and the associated partner potential as discussed in Sec. 3. For the supersymmetric one dimensional potential well it can be very easily verify that the associated partner potential is given by the following equation: A look at Eq. (7.4) immediately suggests that the partner potential is remarkably different from the original one unlike the harmonic oscillator case whose partner potential is identical to that of the original one. The energy eigenfunctions and eigenvalues associated with the partner potential can easily be calculated as: n (x) = 2 (n + 2) 2 − 1 (n + 2) cos((n + 2)πx) − cot(πx) sin((n + 2)πx) (7.4) E (2) n = (n 2 + 4n + 3) π 2 for n ∈ {0, 1, 2, . . . } (7.5) To compute the correlators we need the partition function and the matrix elements of the position operator between any two arbitrary energy states. However while computing the matrix elements if one of the state is the ground state then the expression can be written in a closed form as in the ground state there is no contribution of the partner Hamiltonian.
(1 + m) cos(1 + m)πx − cot πx sin(1 + m)πx (7.8) Substituting the above expressions in the eigenstate representation of the correlators obtained in Sec. 5, the correlators for 1D SUSY potential well can be calculated. Now it is important to mention here that the OTOCs computed from this particular model only can be presented in terms of integrals which at the end we need to evaluate using numerical computation. So to avoid writing complicated mathematical expressions in terms of huge size integrals and instead of presenting a detailed calculation here we provide and discuss the results in the later Sec. 11 of this paper.

General remarks on the classical limiting interpretation of OTOCs
In this section our prime objective is to study the classical limit of the thermal OTOCs derived explicitly in the previous sections of this paper in the context of supersymmetric quantum mechanical systems. This computation is essential to understand the time and the temperature dependent behaviour of the two and four point thermal OTOCs in the classical limit. By styudying the behaviour in this limit one can check the consistency of the result obtained of the quantum randomness from the computed thermal OTOCs in the previous sections.
The strategy adopted in calculating the classical limit of OTOCs is usually replacing the quantum mechanical commutator bracket with the usual classical Poisson bracket defined by the following equation: In the above equation, the q i and the p i are the generalized coordinates and momenta and f and g are functions of these coordinates and momenta i.e.
It can be readily checked that the commutator bracket of two quantum mechanical operators satisfy the same properties as that of the classical Poisson bracket and it can be viewed as the outcome of the following limit on the commutator bracket: The thermal average of the correlators is carried out by the trace operation in the quantum case which can be further simplified in the eigenstate representation. For it's classical counterpart the trace operation is replaced by the phase space integral in the classical limit. For a quantum mechanical model in the context of supersymmetry usually one would expect that there are two generalised coordinates and momenta, one coming from the original Hamiltonian of the system and the other one from its supersymmetric part. Below we write down the generic expressions for the Poisson Brackets involving the position and the momentum operators and provide a general expression for the classical limit of the two and four point classical versions of the correlators.

Classical limit of four point Canonical OTOCs
As discussed earlier any supersymmetric quantum mechanical Hamiltonian is associated with a partner Hamiltonian. For a Supersymmetric Harmonic oscillator the potential associated with the original and the partner Hamiltonian are exactly equal apart from a constant factor as shown in Sec. 6.1. Hence, the classical solutions of the dynamical operators takes identical forms and is obtanied by trivially solving the following differential equation: In solving the above differential equation we take the initial position and momentum to be x(0) and p(0). The classical solutions of the operators obtained by solving the above differential equation and the mentioned initial conditions are given by: We want to calculate the classical limit of OTOC using the position operators at different times. Using the distributive property of Poisson Bracket {x(t 1 ), p(t 2 )} can be expanded for our supersymmetric case in the following way : Each of the terms of the above equation can be evaluated using the definition of classical Poisson Bracket and has been done below: (9.6) which implies that the contributions from both the bosonic and the fermionic part of the Poisson brackets are exactly identical for supersymmetric one dimensional harmonic oscillator. Therefore, the Poisson Bracket of the position and the momentum at different times is given by: We want to calculate the classical limit of OTOC using the position variables at different times. Using the distributive property of Poisson Bracket {x(t 1 ), x(t 2 )} can be expanded for our supersymmetric case in the following way: Each of the terms of the above equation can be evaluated using the definition of the classical Poisson Bracket and has been done below: . (9.10) which implies that the contributions from both the bosonic and the fermionic part of the Poisson brackets are exactly identical for supersymmetric one dimensional harmonic oscillator. Therefore, the classical Poisson Bracket of the position variables at different times is given by: We want to calculate the classical limit of OTOC using the momentum variables at different times. Using the distributive property of Poisson Bracket {p(t 1 ), p(t 2 )} can be expanded for our supersymmetric case and the non zero contribution are given by the following equation: Each of the terms of the above equation can be evaluated using the definition of Poisson Bracket and has been done below: which implies that the contributions from both the bosonic and the fermionic part of the Poisson brackets are exactly identical for supersymmetric one dimensional harmonic oscillator. Therefore, the classical Poisson Bracket of the momentum variables at different times is given by: Finally, the classical limit of OTOC of two point thermal classical version of the OTOCs are given by the following expressions: Classical limit of two point OTOCs for SUSY 1D HO Classical limit of four point OTOCs for SUSY 1D HO C (1) (t 1 , t 2 ) = 4 cos 2 (ω(t 1 − t 2 )), From the classical limit result of the Supersymmetric Harmonic oscillator it is clear that the classical statistics do not produce the quantum result. Though we get a similar time dependence which is periodic both in the classical and the quantum case, the results are not identical to each other. The prime difference in the result lies in the fact that the quantum result depends on the energy eigenstate. This explicit dependence of the OTOC on the energy eigenstates is what prevents the quantum results to give the classical result in the high temperature limit. The important factor to note here is the appearance of the factor 2 in the two point classical correlators. It was already discussed in Sec. 3 that for every potential in supersymmetry there is an associated partner potential. For the case of Supersymmetric Harmonic oscillator it can be seen that the structure of the partner potential is exactly similar to the original potential differing in only an overall constant factor. So it is very easy to understand that the classical solutions for both the potentials will be exactly identical. Now while calculating the Poisson Bracket in the context of Supersymmetric, the bosonic and the fermionic part can be considered as two degrees of freedom. Due to the similar solutions of the dynamical variables the contribution of the bosonic and the fermionic degrees of freedom for the Supersymmetric harmonic oscillator are exactly identical, which adds up to give twice the result obtained from one degree of freedom. Another important point to note that the classical limit of the correlators do not depend on the initial conditions of the dynamical variables.

Classical Limit of OTOC for Supersymmetric 1D box
The case of supersymmetric infinite potential well is not as trivial as that of the Harmonic oiscillator. The associated partner potential is not identical to that of original potential. As has been derived in the paper [88] the partner potential for a 1D box of unit length is given by whereas the well known original potential is given by The classical solutions of the dynamical operators for the original potential are pretty trivial to calculate and is obtained by trivially solving the following differential equation The classical solutions of the operators in this case can be explicitly written as: where x 1 (0) and p 1 (0) are the initial position and momentum of the particle moving in the potential V 1 (x). However, the classical solution of the particle moving in the potential V 2 (x) can be obtained by solving the following differential equation: The above differential equation can be solved explicitly to give the following solutions of the dynamical variables: The constants c and c 1 can be fixed from the initial conditions by taking the initial position and momentum to be x 2 (0) and p 2 (0) respectively and can be written as: It is to be noted that the classical solution written above is valid before the particle bounces at a boundary. After it experienxes a bounce at the boundary the momentum changes its direction i.e p(t) → −p(t). We take this fact into account by considering the infinitesimal deviation of the initial position and fixing the momentum as: The bouncing of the particle at the boundary is given by the time evolution i.e.
δx(t) = (−1) n δx(0) after the n th bounce. We are interested in calculating the Poisson Bracket of the position operator at a certain time with the momentum operator at another time which gives us the classical limit of the correlator of first kind. Similarly the classical limit of the correlator of the second kind is obtained from the Poisson Bracket of the position operators at two different times. The classical limit of the correlator of the third kind is given by the Poisson Bracket of the momentum operators at two different times. The Poisson Bracket for the classical limit of first kind of correlator is given by: Using the classical solutions obtained in Eq. (10.4) the Poisson Bracket involving the position and momentum of the particle moving in the potential V 1 can be written as: In a similar manner, the Poisson Bracket involving the classical position variables of the particle associated with the original potential can be calculated as: The Poisson Bracket involving the classical momentum variables of the particle associated with the original potential can be calculated as In the same manner the Poisson Bracket relations of the particle associated with the partner potential can be calculated in the following way: where the expressions of X 2 and Q 2 are given in Appendix B. The Poisson Bracket involving the classical position variables of the particle associated with the partner potential can be calculated as: where the symbols K 1 and K 2 represent the following terms Similarly, the Poisson Bracket involving the classical momentum variables of the particle associated with the partner potential can be calculated as The explicit expressions for the symbols used in the above equations has been written in the Appendix B. Therefore, the classical limit of OTOC of 2-pt correlators are given by Classical limit of two point Canonical Correlators for SUSY 1D infinite potential well Classical limit of four point Canonical Correlators for SUSY 1D infinite potential well The solutions of the classical dynamical variables obtained for the partner potential associated with the 1D infinite well potential is not trivial as obtained in Eqs. (10.7) and (10.8). The Poisson Bracket involving the partner potential degrees of freedom have complicated terms whose explicit expressions are given in the Appendix B. An important observation here is that the contribution of the two degrees of freedom,( one from the original potential and the other from the partner potential) are not identical, which is expected as the structure and hence the classical solutions of the dynamical variables are not same for both the cases. On careful observation of the Poisson Bracket relations of the partner potential degrees of freedom, it can be seen that there is an explicit dependence on the initial values of the dynamical variables.

Numerical Results
In this section we do the following studies to ascertain the randomness properties of microcanonical and canonical correlators of the two integrable models we have considered : • Study A : The time dependence of 2-pt and 4-pt micro-canonical correlators for four different states : m = 1, m = 2, m = 5, m = 10. This demonstrates the comparative behavior of micro-canonical correlators for different states under time evolution.
• Study B : The time dependence of 2-pt and 4-pt canonical correlators for three different temperatures : T = 10, T = 50, T = 100 at fixed time t = 0.5. We have chosen units such that the Boltzmann Constant, k B = 1, so that the inverse temperature, β = 1/T . This demonstrates the comparative behaviour of canonical correlators for different temperature under time evolution.
• Study C: The temperature dependence of the 2-pt and 4-pt canonical correlators in the temperature range : 10 ≤ T ≤ 100. This demonstrates the explicit temperature dependence of the canonical correlators.
Technically, one has to take the full infinite dimensional Hilbert Space associated with the Supersymmetric 1D Potential Well for computing the correlators but it is not possible to do that in a numerical evaluation and hence one must choose some finite number of states in the Hilbert Space. This choice of a finite value of total number of states will result in an error and this kind of error is known as Truncation Error. The terminology refers to the fact that the error is arising because we have truncated the number of states to a finite value. Here, we have chosen the truncation to be such that all states ≤ N trunc = 10 for numerical evaluation.
From the eigenstate representations for the correlators : y (t 1 , t 2 ) and C (1,2,3) (t 1 , t 2 ) given in Sec. 5.8 we know that they depend only on t 1 − t 2 and hence we have defined t = t 1 − t 2 . So, the negative values of t refer to the case of t 2 > t 1 and so on. We proceed to the discussion of obtained numerical results.
To make any bold comments about whether we are observing randomness or not we need to take the commutator brackets as operators O i & then calculate : and check whether ∆O i depends on time t or not. If it doesn't then what we have is just statistical fluctuation which is arising from the inherent quantum nature of the systems and not from randomness. For example : Consider a QHO which as ∆x = 0 but it is also independent of time, this means that at each instant of time if we take a large number of copies then there will be some statistical variation in the values of x across different copies but that variation will be found exactly the same at each instant of time. If ∆O i = f (t) then it is indeed a true signature of quantum randomness perhaps depending on functional form of f (t). Since, we have calculated both O 2 i β and O i β it is really easy to check for any signature of randomness. Furthermore, we can use the same method to calculate higher moments than ∆O i and get more and more sensitive measurements for randomness which is helpful if the f (t) above turns out to be some simple periodic function like f (t) ∼ cos (at) and we are not satisfied with it because we have really sensitive and amazing technology which can probe for even finer randomness signatures.
To elaborate further consider the following. We don't really know what is quantum uncertainty and hence to work with it we replace it by statistical uncertainty using ensembles and so on. Then uncertainty of an operator O i represents its fluctuation across ensembles and it is not really a true measure of randomness. The fluctuation is merely a sampling fluctuation. Say, we take an ensemble consisting of 10000 copies of a system then the uncertainty tells us the variation in the measured value of O i across the copies / ensemble at a particular time. Now, if we obtain the ∆O i as a function of time then we can easily probe for randomness depending on the sensitivity of our instrument. If our instrument is really sensitive then we can even use higher moments to check for stronger and stronger conditions of randomness. Obtaining ∆O i as a function of time means that the fluctuation itself is changing and hence this can be a true signature of randomness this ofcourse depends on what kind of function we get.

Supersymmetric 1D Infinite Potential Well
For numerical evaluation we have chosen : L = 1 & 2m = 1, where L is the length of the box in which a particle of mass m is confined. We also consider = 1.
• In Figs. 5, 6 and 7 we perform the Study A on the 2-pt micro-canonical correlators y for Supersymmetric 1D Infinite Potential Well.
-We observe that the correlators y  m (t 1 , t 2 ) and we plan to do the same in a work which is to appear very shortly.
-The amplitude of the correlators increases with increasing m which primarily comes about because we have : in the correlator's eigenstate representation which increases with increasing state number.
-It is also observed that with higher and higher excited state the number of nodes increases and takes on the shape of a wave-packet. This is so because as we go to higher and higher excited states the high frequency modes become more and more prominent. So the comparison of the scaling of the number of nodes here with non-Supersymmetric case is also important job to perform in near future.            • We present the results of performing Study B and Study C on 2-pt canonical correlators Y (i) (t 1 , t 2 ) as follows -In Figs. 8, 9 and 10 we perform Study B on the 2-pt canonical correlators Y (1,2,3) (t 1 , t 2 ). We observe that the correlators shows periodic behaviour for the different chosen values of temperature. We observe that for Y 1 (t 1 , t 2 ), the mid temperature value 50 shows the minimum amplitude, whereas lower value of temperature (10)   in the amplitude of the correlator for temperatures in the higher value range as can be seen from Fig. 8. In correlators Y 2 (t 1 , t 2 ) and Y 3 (t 1 , t 2 ) however we follow a gradual pattern of decreasing amplitude with increasing in temperature as observed from Fig. 9 and Fig. 10. To have a better understanding of the temperature dependence of the 2-pt correlators we plot them with varying temperature keeping the time constant.
-In Figs. 11, 12 and 13 we perform Study C on the 2-pt canonical correlators Y (1,2,3) (t 1 , t 2 ). We plot, respectively, Y (1,2,3) (t 1 , t 2 ) which are the thermal or canonical correlators corresponding to y (1,2,3) m (t 1 , t 2 ) respectively with respect to temperature. It is clearly visible that for very low temperatures the canonical correlators are constant. However after a certain value of the temperature the correlators decays rapidly and falls off to zero within a small temperature range   m (t 1 , t 2 ) the periodicity is ∆t 0.35 which is roughly half of the corresponding 2-pt micro-canonical correlator. We note that this is approximately the same periodicity, within numerical error, observed in the case of non-Supersymmetric 1D Infinite Potential Well as obtained by Hashimoto et al. [70]. Hence, we conclude that introducing Supersymmetry in integrable QMcal models does not affect the periodicity of four-pt micro-canonical correlators. At present, there are no studies for the non-Supersymmetric case for c (2,3) m (t 1 , t 2 ) and we plan to do the same in a work which is to appear very shortly.
are symmetric about t = 0 which means that to these 4-pt microcanonical correlators it does not matter whether t 1 > t 2 or t 1 < t 2 .  (t 1 , t 2 ) respectively with respect to temperature. It is observed that for low temperatures the 4-pt correlators show negligible value. However after a certain threshold temperature the 4-pt correlators increases and then saturates to a certain finite value. plots suggests that 2-pt correlators are a better probe for understanding Quantum randomness at low temperatures whereas at high temperatures it is actually the 4-pt correlators which is more significant. However in the mid-temperature range the 2-pt and the 4-pt correlators have exactly opposite behavior. Hence, to understand the significance of temperature in this range on any Supersymmetric Quantum mechanical model having knowledge of both 2 and 4-pt correlators is of utmost importance.

Supersymmetric 1D Harmonic oscillator
For numerical evaluation we have chosen : ω = 1 & 2M = 1, where ω is the frequency of the oscillator in which a particle of mass M is confined. We also consider = 1.
-We observe that the correlators y   for m=1. However with further increase of m the amplitude of the correlator shows negligible change and the amplitudes of the higher states almost overlaps. This suggests that for the lower energy states the micro-canonical correlators depend on the energy state in which it is calculated. However this state dependency goes away when calculated for the higher energy states. This can also be understood from the analytical expression obtained for the micro-canonical correlators obtained in Sec. 6 (calculated for Harmonic Oscillator of unit mass i.e M = 1). The micro-canonical correlators have a non trivial state dependence in the form of the factor 1 + m(m + 1) − m(m − 1) which reduces simply to 1 for higher energy states.
-In Figs. 23, 24 and 25 we have also plotted y               thermal or canonical correlators corresponding to y (1,2,3) m (t 1 , t 2 ) respectively with respect to temperature. It is clearly visible that for very low temperatures the canonical correlators are constant. After a certain value of the temperature the amplitude of the correlator shows a gradual increase. It reaches a maximum for a particular value of the temperature and then decays exponentially to zero. -We observe that the correlators c . This is because for the 4-pt correlators there is a mass square dependence unlike the 2-pt correlators which has just mass dependence. The amplitudes of the respective 4-pt correlators can also be found to be exactly half of the amplitudes of its 2-pt counterpart. This is obvious from the time dependent functions appearing in the case of Supersymmetric Harmonic oscillator.
-All c (i) m (t 1 , t 2 ) are symmetric about t = 0 which means that to these 4-pt microcanonical correlators it does not matter whether t 1 > t 2 or t 1 < t 2 .
-In Figs  C (1,2,3) (t 1 , t 2 ). We observe that the correlators shows periodic behaviour for the different chosen values of temperature. We observe that each of the 4pt correlators behave identically with respect to temperature which is exactly opposite in character from the 2-pt canonical correlators. The amplitude of each of them increases with increasing temperature. To have a better understanding of the temperature dependence of the 4-pt correlators we plot them with varying temperature keeping the time constant.
-In Figs. 38, 39 and 40 we present the results of performing Study C on 4pt canonical correlators C (i) (t 1 , t 2 ). Here we plot C (1,2,3) (t 1 , t 2 ) which are the thermal or canonical correlators corresponding to y                   the temperature the amplitude of the correlator starts increasing and finally saturates to a finite value.
• The temperature dependent plots shows the significance of computing the 4-pt and the 2-pt correlators to study the phenomenon of quantum randomness for a Supersymmetric quantum mechanical model. From the plots it can be seen that for very low temperatures the 2-pt correlators shows a certain finite value whereas the 4-pt correlators are almost negligible. On the other hand at very high temperatures the 2pt correlators are almost zero whereas the 4-pt correlators shows certain finite value. This suggests that to understand quantum randomness for a Supersymmetric model at very low temperature, the results from the 4-pt correlators can be misleading and similarly at high temperatures the 2-pt correlators might not be an appropriate quantity to study randomness. However to understand quantum randomness at the mid temperature range both 2 and 4-pt correlators play significant role. We feel this shows the necessity for computing the 2-pt correlators.

Conclusions
To summarize, in this work, we have addressed the following issues to study the OTOC in the context of supersymmetric quantum mechanics: • We apply the computational techniques of recently developing methodology of out of time ordered correlators (OTOCs) to study the phenomenon of time disorder averaging for a given quantum statistical ensemble or quantum randomness for two very well known integrable one dimensional quantum mechanical models viz. Harmonic Oscillator and 1D potential well in the context of one dimensional supersymmetric quantum mechanics. We show that to develop a complete understanding of the underlying randomness in the quantum system, not only the correlators constructed from different operators at different times scales are important, but also the correlators constructed from similar quantum mechanical operators at different time scales, play a pivotal role.
• We have constructed all the temperature independent micro-canonical and temperature dependent canonical un-normalized and normalized version of these OTOCs in the eigenstate representation of the supersymmetric time independent Hamiltonian of the quantum system and represent all of them in a general model independent way. From the previous study it is expected that the OTOCs in the eigenstate representation one should not expect any chaotic behaviour i.e. the exponential growth with respect to the time scale in the correlators. However, one can get a broader knowledge of some other aspects of quantum randomness, which might capture the random behaviour in the correlators in terms of non-chaoticity. From our analysis it is expected that a large class of quantum mechanical models will be covered which provide the signature of quantum mechanical randomness in general.
• The explicit calculation of the 4-pt correlator C (1) (t 1 , t 2 ) for the Supersymmetric Harmonic Oscillator, shows the significance of the introduction of supersymmetry within the context of quantum mechanics compared to the framework of quantum mechanics without having any supersymmetry. Supersymmetry introduces an energy state dependence on the temperature independent micro-canonical correlators which is usually not appear in the framework without having supersymmetry in the quantum mechanics description. Apart from the dependence on the energy states the canonical correlators have an additional dependence on temperature, which is also a different and notable feature compared to the results obtained from quantum mechanical set up without having supersymmetry. This energy state and temperature dependence of the correlators differentiates a Supersymmetric and a Non-Supersymmetric Harmonic oscillator [70] though the time dependence in the OTOCs remains same. Also, particularly for the Supersymmetric Harmonic Oscillator we have found that the normalized four point OTOCs that are made up of same operators at different time scales show exactly same behaviour, which implies they are not independent of each other. However, this statement might not be true for other integrable models. On the other hand, in the un-normalized version of these two correlators we have found exact same time dependence, but the overall frequency dependent normalization factors will be different, which implies that they are proportional to each other in this case.
• The classical limit however matches with the non supersymmetric case apart from a factor of 2 which can be inferred from the fact that in a supersymmetric quantum mechanical model there are two degrees of freedom, one from the original potential and the other associated with the partner potential. The time dependence of the classical and the quantum version of the correlators are exactly identical. However the temperature dependence observed in the quantum case vanishes for its classical counterpart, which is obviously a important finding from our computation.
• We observe a similar temperature and state dependence on the correlators for the Supersymmetric 1D Infinite potential well. However, it is interesting to note that [70] also observed this state and temperature dependence for Non-Supersymmetric Infinite potential well. The behaviour of the only correlator studied in [70] is exactly identical to what we observe for the Supersymmetric case. The correlator showed increase in the amplitude with increasing state number and higher temperatures which is exactly what we observe here. Hence we conclude that introduction of supersymmetry does not introduce new features in the case of 1D Infinite potential well.
• The significance of supersymmetry in 1D potential well however can be realised from its classical counterpart which is markedly different from what is obtained in the nonsupersymmetric case. The classical limit of OTOC for the 1D non supersymmetric infinite potential well is independent of time and is merely a constant, whereas for the supersymmetric case there is a non-trivial time dependence, which is obviously a new finding from our computations.
The future prospects of this work are as follows: • In this paper we have restricted ourselves in considering only 2 and 4-pt correlation functions, to study quantum randomness in various Supersymmetric QM models. However a more generalised approach would be to consider the multipoint correlation functions to have a better understanding about the randomness underlying the system. We have an immediate plan to carry forward the work along this direction very soon.
• The study of OTOCs can be used to understand quantum randomness in various quantum mechanical models which are of prime significance in various condensed matter, nuclear and atomic physics models. Particularly the time dependent Hamiltonians we have not studied in this paper where this eigenstate formalism and related simplification to determine OTOC will not work. In that case one needs to use the general definition and representation of OTOC in terms of the well known Schwinger Keldysh formalism. We have some future plan on that direction as well.
• Use of the other types of correlators defined in this paper can be used to study some of the well understood QM models, to have an insight on the lost information of quantum randomness. We are very hopeful that incorporating the study of these additional correlators which we have defined and evaluated in this paper can able to capture more broader perspective on time disorder averaging phenomena through quantum randomness, rather only give insight about quantum mechanical chaos from the temporal growth of the correlators.
• Very recently in [110][111][112][113][114] the authors have studied various relevant measures for an entangled open quantum system. The study of OTOCs for such type of entangled systems [115][116][117][118] will be quite relevant for understanding the rtime disorder averaging phenomena and chaos for such an entangled OQS.
• Last but not the least, very recent we have proposed a detailed mechanism and framework using which one can able to compute the OTOC within the framework of primordial cosmological perturbation theory by making use of the gauge invariant scalar perturbations and its associated canonically conjugate momenta [69] and finally found out the chaotic like behaviour in the representative cosmological version of OTOC. However, in that paper we have not reported anything about the other possible two operators which we have introduced in this paper. At present we are working on that direction and very hopeful to get the certain non-trivial features of time disorder averaging which might have application to explain various cosmologically relevant phenomena within the evolution history of our universe, which was not explored earlier at all.

A Derivation of the normalization factors for the Supersymmetric HO
To normalize the obtained OTOC we divide it by the thermal average of the dynamical variables considered in calculating the OTOC.
x(t 1 )x(t 1 ) β = 1 Z Tr(e −βH x(t 1 )x(t 1 )) = 1 Z n e −βnω Ψ n |p(t 1 )p(t 1 )|Ψ n (A.1) Now keeping in mind that the ground state is only bosonic whereas all the other eigenstates has a fermionic part associated with it. Hence, separating the ground state from the other states, the above equation can be written as e −βnω ψ n |x(t 1 ) |ψ k ψ k | x(t 1 )|ψ n

Term B
Now we are going to explicitly show the calculation of term A and term B. Expanding the position operators using the Heisenberg picture for the evolution of operators Term A =   x nk (t 1 )x kn (t 1 ) = e −βω 4ω + 1 4ω n>0 e −βnω 2n + n(n + 1) + n(n − 1) (A.4) x(t 1 )x(t 1 ) β = 1 4Zω 1 + e −βω + n>0 e −βEn [2n + n(n + 1) + n(n − 1)] (A.5) A similar calculation is carried out for the thermal average of the product of the momentum operators.