Minimal Products of Coordinate and Momentum Uncertainties of High Orders: Significant and Weak High-Order Squeezing

We consider the problem of minimization of products of mean values of the high powers of operators x and p. From this point of view, we study several two-term superpositions of the Fock states, as well as three popular families of infinite superpositions: squeezed states, even/odd coherent states, and orthogonal even coherent states (or compass states). The new element is the analysis of products of the corresponding (co)variances and the related generalized (Robertson–Schrödinger) intelligent states (RSIS). In particular, we show that both Fock and pure Gaussian homogeneous states are RSIS for the fourth powers (but not for the sixth ones). We show that lower bounds of the high-order uncertainty products can be significantly below the vacuum values. In this connection, the concept of significant and weak high-order squeezing is introduced.


Introduction
The exact inequality for the product of the second order statistical moments of the dimensionless coordinate and momentum operators (in the system of units withh = 1), is known since the paper by Kennard in 1927 [1]. Moreover, it is well known that the equality in (1) is achieved for the vacuum state |0 of the quantum harmonic oscillator with x 2 = p 2 = 1/2. Using the inequality Â 2 ≥ Â 2 (which holds for any Hermitian operatorÂ) one immediately arrives at the inequality Π (4) ≡ x 4 p 4 ≥ 1/16.
Soon after paper [2], Lynch and Mavromatis [4] looked for the minimal value of the product Π (N) ≡ x N p N in the restricted class of states with equal mean values, x N = p N . They found that the minimal value of Π (N) for such states coincides with the minimal eigenvalue squared λ 2 of the equation Using approximate variational schemes, they obtained the minimal value (Π (4) ) LM = (0.6984) 2 = 0.4878 = 0.8672 (Π (4) ) vac (5) in the specific superposition of the Fock states of the form |ψ K = ∑ K n=0 c n |4n , with K = 6 and certain values of coefficients c n (not given in that papers). Recently, Equation (4) was studied in [5], where the same value λ (4) min = 0.6984 was confirmed. The value (5) is the best lower limit for the product Π (4) known until now.
An interesting infinite-dimensional analog of the state |ψ K is the superposition of four coherent states with equal amplitudes and phases shifted by π/2, Such kinds of states were studied for a long time from different points of view, due to their interesting properties. They were named "four-photon states" [6], "orthogonal-even coherent states" [7], "pair cat states" [8], "compass states" [9], and "four-headed cat states" [10]. The state (6) is a special case of large families of "circular states" [11][12][13] and "crystallized cat states" [14]. For extensive lists of references one can consult [15,16]. Following [7], we use the abbreviation OECS for the state (6). It was shown in [7] that the value (Π (4) ) min ≈ 0.49 can be obtained in the state (6) with |α| ≈ 0.67. Later, this result was confirmed in [17,18]. It is worth noting that the relative probability of the state |4 in the superposition (6), p 4 /p 0 = (0.67) 8 /24 = 0.00169, is very close to the relative probability 0.00166 of the state |4 in the superposition |ψ 4 (3).
One of our goals is to analyze some known families of quantum superposition states, in order to see whether they can give the values of Π (4) smaller than 0.49, without the restriction x N = p N . The simplest natural candidates are superpositions of two coherent states or superpositions of states |ψ k (3), containing even Fock states only. This is explained by the fact that operatorsx 2 andp 2 mix the Fock states with the same parity only. Because we need the vacuum state |0 as the basic element, only even Fock states can lead to useful interference effects. In principle, one can consider various families of states, which are characterized by arbitrary numbers of independent parameters. However, explicit expressions for the quantity Π (4) are so complicated in the majority of cases that we confine ourselves to the cases of no more than four real parameters.
However, let us look at the ratio F (2n) ≡ (Π (2n) ) min /(Π (2n) ) vac , where (Π (2n) ) min is the minimal value of the product Π (2n) . Using the values of (Π (2n) ) min calculated in [4] (the best values known up to now), we obtain the following numbers: We see that this ratio diminishes with increase of n, but the decrease rate slows down for big values n.
We study the asymptotic behavior of F (2n) for the two-Fock superpositions (3) in Section 8, comparing our results with those that were obtained in [4]. The subject that was not addressed in the available literature is the minimal value of the product of variances of operatorsx N andp N . Remember that variances and the covariance of two Hermitian operators,Â andB, are defined, as follows, Let us define the variance uncertainty product (VUP) of the order 2n as Obviously, Π (2n) σ > 0 for any normalized quantum state. It is also obvious that Π σ ? The answer is well known for n = 1: in view of the Heisenberg-Kennard inequality [1] σ x σ p ≥ 1/4, one has the minimal value (Π (2) σ ) min = 1/4, and this minimum is achieved for all coherent states of the harmonic oscillator. However, the minimum of Π (2n) σ with n > 1 was never discussed. This issue is considered in Sections 2, 5 and 6 for n = 2. The generalization of inequality σ x σ p ≥ 1/4 for arbitrary Hermitian operatorsÂ andB was found by Robertson and Schrödinger [19,20]. It can be written as follows, Two equivalent formulas for R, which can be useful in some cases, are The problem with inequality (8) consists in the fact that the quantity Y AB is state dependent for many operators, different fromx andp. Moreover, Y AB can be equal to zero in many quantum states. Therefore, the following non-negative quantities can be studied in connection with inequality (8): The name "intelligent states" (IS) was coined in [21] for the spin (angular momentum) states satisfying the equality Z = 0. The difference between "intelligent" and "minimum uncertainty" states for the angular momentum operators was discussed in [22]. A little later, the same name was accepted in connection with the non-compact group SU(1, 1) [23]. Further extensions of the concept of "intelligent" states were made, e.g., in [24][25][26][27]. Of course, the word "intelligent" is nothing but a label, but this name has been adopted nowadays by the majority of researchers. Later on, the states satisfying a more general equality R = 0 were named "generalized intelligent states" (GIS) in [28]. The names "Schrödinger intelligent states" and "Robertson intelligent states" were also used in [29,30]. Many references to the papers that are devoted to IS and GIS for various sets of operators can be found in [15]. The physical meaning of IS and methods of their generation were discussed, e.g., in papers [31][32][33][34]. "Constrained" intelligent states were introduced in study [35]. We shall use the name "ordinary intelligent state" (OIS) in the case of Z = 0, and the name "Robertson-Schrödinger intelligent state" (RSIS) when R = 0.
In Section 2, we consider the operatorsÂ =x N andB =p N and Gaussian pure (squeezed) states, trying to answer the following questions: (1) When the equalities R = 0 or Z = 0 can be satisfied? (2) What is the minimal value of the quantity D? In Section 3, we show that the Fock states are RSIS (but not OIS) for N = 2. Sections 4-8 are devoted to the detailed analysis of the products Π (2n) and Π (2n) σ in various superposition states, described by a small number (from 2 to 4) of real parameters. The main goal is to minimize these products with respect to the chosen parameters. The results of that sections naturally lead us to some thoughts about the best definition of the concept of the high-order squeezing. This subject is discussed in Section 9. The last Section 10 contains a brief discussion of main results that were obtained in the paper, together with outlines of possible directions for further studies.

Robertson-Schrödinger Relations with High-Order Moments for the Gaussian States
It is known since [1] that R ≡ 0 for A = x and B = p in any pure Gaussian state, described by means of the normalized wave function where a, b, c, d are real numbers with a > 0. If b = 0 and a = 1, then (10) is the well known Schrödinger's non-dispersive wave packet, called nowadays the Glauber-Sudarshan coherent state. The states with arbitrary a and b = 0 were named "correlated coherent states" in [36]. The states with c = d = 0 are nothing but the "squeezed vacuum states". For a detailed history of these states, one can consult review [15]. Coefficients of the wave function (10) can be written in terms of mean values and variances: The high-order moments of operatorsx andp in the state (10) are given by the formulas (see Appendix A for the proof) where H n (z) is the Hermite polynomial. The right-hand sides of expressions in (12) are always real, due to the even/odd parity of the Hermite polynomials for even/odd values of power n. If x = p = 0, then A simple formula can also be obtained for the mean value x npn , if x = p = 0: where P n (z) is the Legendre polynomial. The consequence of (14) for n = 1 is the known formula for the covariance in the Gaussian state (10): σ xp = −b/(2a).

N = 2
ForÂ =x 2 andB =p 2 , we have In view of the x − p symmetry, we can write Therefore, the minimum of Π (4) σ is achieved for the states with b = c = d = 0, i.e., for all uncorrelated vacuum squeezed states (here, we restore the Planck constant): In addition, where α = a + ib and β = c + id. Using the real part of (16), together with (A4), we obtain a very simple formula for the covariance: Because [x 2 ,p 2 ] = 2i (xp +px), we immediately obtain the value Y AB = 2 σ xp + x p . Therefore, This quantity equals zero for c = d = 0 and arbitrary values of a and b. Consequently, any squeezed vacuum state is the Robertson-Schrödinger intelligent state not only for the pair of operators (x,p), but also for the pair (x 2 ,p 2 ). Actually, this fact can be derived from the results of study [28], since linear combinations of operatorsx 2 ,p 2 andxp +px are generators of the su(1, 1) algebra. The new result obtained here is that the presence of linear terms in the argument of the exponential function (10) breaks the fourth-order generalized "intelligence", in contradistinction to the second-order "intelligence". Another interesting formula is The second and third terms in the right-hand side are always non-negative, going to zero simultaneously for x = p = 0 only. On the other hand, σ 2 . Consequently, the correlated state (10) is the OIS for the pair (x 2 ,p 2 ), if c = d = 0 and b = ±a, when σ xp = ∓1/2 and the correlation coefficient This result is new. Note that σ AB = 0, in this case, so that the state is uncorrelated from the point of view of the (x 2 ,p 2 ) pair.

Arbitrary N, Vacuum Squeezed States
For the vacuum squeezed (uncorrelated) states (b = c = d = 0), we obtain, from (13), the uncertainty product of the 2nth order, as given by Equation (7). The same value corresponds to Π (2n) σ , if n is any odd number. If n = 2m, then However, using the known Stirling formula for n 1, n! ≈ √ 2πn(n/e) n , one can see that the relative contribution of the second term inside the figure brackets rapidly decreases with m as 4 −m √ 2.

Fock States as Intelligent States
Simple formulas can be obtained for various combinations of products of the second-order moments of operatorsÂ =x 2 andB =p 2 in an arbitrary Fock state |n : They follow from the standard expressions for the canonical dimensionless operatorsx andp in terms of the bosonic annihilation and creation operators, Consequently, In addition, Y AB = 0, as soon as xp +px = 0. Therefore, D = R = 0 for any Fock state. This means that all of the Fock states are Robertson-Schrödinger intelligent states for the pair of operatorsÂ =x 2 andB =p 2 . However, they are not ordinary intelligent states for this pair. The product σ A σ B equals 1/4 for n = 0, and it grows rapidly with increase of n. However, we can expect that some superpositions of the vacuum and a few first excited states can give smaller values of this product. This conjecture is analyzed in the next sections.

Nth Order Products in the Two-Term Even Fock Superpositions: General Formulas
In this section, we consider the mean values of even powers of the coordinate and momentum in the two-term superpositions of the even Fock states |ψ 2j (3). In this case, we can write All of the coefficients (19) are real non-negative numbers. This is the consequence of formulas (17).
Using the form of wave functions in the coordinate and momentum representations, with the same coefficients as in (18). According to Equation (18), the minimum of ψ 2j |x 2n |ψ 2j with a fixed value of |β| is achieved for Re(α * β) = −|αβ|. In order to find the absolute minimum, we have to minimize the function Here, we dropped all indexes. Calculating the derivative d f /dz, we arrive at the quadratic equation g 2 z 2 − g 2 z + ρ 2 = 0, where g 2 = 4ρ 2 +ν 2 . It has two real solutions. We need the solution z * , which is close to zero. Afterwards, These formulas are exact. For small values of |β| (which are most relevant for our purposes, as will be shown in the following sections), one can use the parabolic approximation f (|β|) = 1 +ν|β| 2 − 2ρ|β|, which results in simple formulas |β| min =ρ/ν and f min = 1 −ρ 2 /ν. If j is odd, then the minimum of ψ 2j |x 2n |ψ 2j corresponds to the maximum of ψ 2j |p 2n |ψ 2j , according to Equations (18) and (20). In this case, (Π (2n) ) min = χ min (Π (2n) ) vac , where χ min is the minimum of function χ(z) = (1 +νz) 2 − 4ρ 2 z(1 − z). The exact position of minimum z * and exact minimal value χ min are as follows: Explicit expressions for the coefficients (19) in the most general case are as follows (see Appendix B): ν (2n) 2j where F(a, b; c; z) is the Gauss hypergeometric function.

Fourth-Order Variance Products in States |ψ 2 and |ψ 4
General formulas for the variances σ A and σ B are rather involved for arbitrary numbers n and j. In this section, we consider in detail the case of n = 2, when all calculations can be performed in a relatively simple way.

Superpositions of the Vacuum and Fourth Fock States
For n = j = 2, the general scheme of the preceding section results in the following mean values: Therefore, For the fixed values of |α| and |β|, the minimum of (26) is achieved for Re(α * β) = −|αβ|. Thus, we arrive at the function Its minimum in the interval 0 < |β| < 1 can be found numerically:

Superpositions of Three Even Fock States
Now, let us consider the three-parameter states of the form where all three coefficients C j can be assumed to be real and non-negative. General formulas in this case are rather involved. For example, If C 1 = 0, then both products, Π (4) or Π (4) σ , attain their minimal values for χ = π. Using this value, we can write Π (4) as Π (4) The partial derivative of this function with respect to C 2 1 at C 2 1 = 0 equals Using the value C 2 = 0.041 (when function Π (4) (C 2 , C 1 = 0) attains the minimum), we obtain ≈ 12.7 − 10.5 cos 2 (φ). This derivative is positive for any phase φ. Consequently, the addition of the second Fock state |2 can only spoil the result obtained for the two-term superposition |ψ 4 . The same conclusion holds for the product Π σ (although the calculations are more cumbersome).

Two Coherent States
The first simplest superpositions of two coherent states were introduced in [37] under the names even/odd coherent states: Here, |α is the well known Glauber-Sudarshan coherent state, satisfying the equationâ|α = α|α withâ = (x + ip) / √ 2. Remember that wave functions of the coherent states in the coordinate and momentum representations have the form Therefore, Note that |0 + is the ground Fock state |0 , whereas |0 − is the first excited Fock state |1 . The states (28) and their generalizations were discussed in detail, e.g., in reviews [15,38,39]. For our purposes, the following properties of the states (28) are useful: Using these properties, the following mean values can be found: where a = |α| 2 , α = |α|e iφ , t ± = [tanh(a)] ±1 , c n = cos(nφ), s n = sin(nφ), so that c 4 = c 2 2 − s 2 2 . Hence, we have to minimize the following functions: Note that functions Z ± do not depend on the phase φ. Their minimal values (found numerically) are: Z min + = 0.21259 for |α| ≈ 0.52963, Z min − = 1.76386 for |α| ≈ 1.03430.
The same minimal values are obtained for the products (Π (4) σ ) ± , if cos 2 (2φ) = 1. The following mean values can be found for |α| = 0.52963 and c 2 = −1: The product Π ± , as given by Equation (31), is a quadratic polynomial with respect to the variable ξ = c 2 2 . It is easy to verify that the derivative of this polynomial with respect to ξ is negative in the interval 0 ≤ ξ ≤ 1. Consequently, the minimum is achieved at ξ = 1. Thus, we have to minimize the function f (a) = 9 16 + 9 2 at ± − 6a 2 + 9a 2 t 2 ± + 4a 4  − ) min is smaller than the value (Π (4) ) 1 = 225/16 = 14.0625 in the first Fock state |1 , the difference is only about 7%. This example shows that, removing the vacuum state, one cannot achieve a small value of Π (4) . Figure 3 shows the functions ψ + (x) and ψ + (p) (30) for two extremal values of parameter α. . Wave functions of the even coherent states ψ + (x) and ψ + (p) (30) for two extremal values of parameter α, as compared with the ground-state wave functions ψ 0 (x) and ψ 0 (p).

Asymptotic Values of High-Order Products in the Two-Term Fock Superpositions
It is interesting to see the behavior of (Π (2n) ) min for n 1. This can be easily done for the two-Fock superpositions, while using Equations (22) and (23). Coefficients (24) and (25) have the following asymptotic forms: where we neglect corrections of the relative order n −2 inside square brackets. Afterwards, Equations (22) and (23) assume the form If j is an even number, then ψ 2j |x 2n |ψ 2j = ψ 2j |p 2n |ψ 2j , so that F (2n) ≡ (Π (2n) ) min /(Π (2n) ) vac = f 2 min . Therefore, this ratio is minimized for j = 2, i.e., in the state |ψ 4 , which has preferences over other-two term superpositions for n = 2 and n = 3, as we saw in the preceding sections: The authors of [4] minimized the mean value of operatorx N +p N for the state |ψ 4 with an arbitrary real coefficient α. We re-write formula (33) from that paper in a more convenient form: If N 1, then (35) assumes the form √ F (N) ≈ 8/N, which coincides with (34) for N = 2n. However, it is worth noticing that such a simple asymptotic formula gives accurate enough values of the ratio F (N) for very big values of N only. Note that the authors of [4] interpreted Formula (35) incorrectly, saying that the uncertainty product goes to zero as N → ∞. Of course, this was an obvious oversight: although the ratio F (N) tends to zero as N → ∞, the total uncertainty product Π (N) = (Π (N) ) vac F (N) remains very big in this limit. Figure 6 shows the ratio (Π (2n) ) 2j /(Π (2n) ) vac ≡ f 2 , calculated using exact Formula (21), as function of |β| for φ = π, in the state |ψ 4 (i.e., j = 2), for different values of power 2n. We see that the asymptotic Formula (34) gives approximate minimal values higher than actual ones. Exact and asymptotic values of f 2 min are compared in Figure 7. If j is an odd number, then F (2n) = χ min ∼ 4j 2 /n. Therefore, the state |ψ 4 is the best one among all two-term Fock superpositions from the point of view of the minimization of Π (2n) . However, the state |ψ 2 is preferable from the point of view of the high-order squeezing: see Section 9. We do not study the variance product Π

Significant and Weak High-Order Squeezing
The concept of high-order squeezing was introduced by Hong and Mandel [40]. According to their definition, the state |ψ is squeezed to the 2nth order with respect tox, if the mean value ψ|(∆x) 2n |ψ is less than the mean value of (∆x) 2n in the coherent state, i.e., (∆x) 2n < Λ vac 2n ≡ 2 −n (2n − 1)!!. In particular, for n = 2 we have the requirement (∆x) 4 < 3/4. Hong and Mandel showed that the usual squeezed states are squeezed to any even order 2n. There are other definitions of the higher-order squeezing [41][42][43][44], but they are not related directly to the powers of coordinate and momentum operators.
The choice of the coherent (or vacuum) state as the reference one is natural for the usual (second-order) squeezing, because the minimal value of the uncertainty product Π (2) σ (as well as Π (2) ) is achieved in this state. However, minimal values of high-order products are achieved in non-Gaussian states. Therefore, we suggest to say that the x variable is significantly squeezed to the 2nth order, if (∆x) 2n < Λ 2n , where Λ 2n = (Π (2n) ) min . If Λ vac 2n > (∆x) 2n > Λ 2n , then the 2nth order squeezing can be considered as weak. In particular, using the results from [4], we can use the following values: In addition, we can consider the σ-squeezing of the fourth order, defining its weak form by the inequality σ x 2 < 1/2, and the significant form by the inequality σ x 2 < (Λ σ ) 4 = (Π (4) σ ) min ≈ 0.44. Note that, according to our definition, the simultaneous significant high-order squeezing is impossible.
For example, all of the Gaussian states with b = c = d = 0 (see Section 2) possess some amount of weak high-order squeezing with respect to x for any value a > 1 and any value 2n. However, the significant fourth-order squeezing can be achieved for a > a 4 = 1.035, and the significant σ-squeezing of the fourth order exists for a > 1.066. For higher values of 2n, the minimal value of parameter a resulting in the significant squeezing is given by the formula a 2n = [F (2n) ] −1/(2n) , according to Equation (13). This minimal value a 2n slowly grows initially with the power 2n, but for 2n ≥ 16 it slowly decays to unity as n → ∞. In particular, using the best available values of Λ 2n from Equation (36), we obtain a 4 ≈ 1.035, a 6 ≈ 1.083(1.0715), a 8 ≈ 1.12(1.093), a 10 ≈ 1.105, a 16 ≈ 1.112.
The values in parentheses for 2n = 6 and 2n = 8 are obtained using Equation (22), derived for the two-Fock superpositions. For other values of 2n, the results of the two-Fock approximation practically coincide with that of [4]. The asymptotic formula is a 2n ≈ exp n −1 ln(n/4) ≈ 1 + n −1 ln(n/4) for n 4. Therefore, homogeneous Gaussian states with parameter a > 1.12 (or σ x < 0.45) are significantly squeezed with respect to the x-coordinate to any order 2n. Figure 8 shows minimal values a 2n , calculated with the aid of Equation (22), as compared with two asymptotic expressions.  Figure 10 shows x 6 as function of parameter |β| for the two-Fock superpositions |ψ 2 and |ψ 6 (3) with φ = π. In the case of |ψ 2 , the minimal value x 6 min = 0.510146 is achieved at |β| = 0.169102. For |ψ 6 , the minimal value x 6 min = 1.85904 is achieved at |β| = 0.00475743. In this case, we only observe a very weak sixth-order squeezing.
The Hong-Mandel squeezing of the 4th and 6th order in the even coherent states was studied, e.g., in [45][46][47]. Arbitrary values of 2n were considered in [48]. The most general superpositions of two coherent states were considered in [49,50] in connection with the 4th order squeezing, and in [51,52] in connection with the 2nth order Hong-Mandel squeezing. The enhancement of the 2nth order Hong-Mandel squeezing in superpositions of even coherent states and the vacuum state was demonstrated in [53]. Two-mode generalizations of even/odd coherent states were studied in [54,55]. The high-order Hong-Mandel squeezing in the Gaussian mixed states (squeezed thermal states) was studied, e.g., in paper [56]. Squeezed number states were studied from this point of view in [57]. Multiphoton non-Gaussian states were discussed in [58]. Paper [59] showed the existence of high-order squeezing in specific four-photon states. Combinations of squeezed and number states were considered in [60]. Two-mode squeezed states were considered in [61]. Higher-order squeezing in the photon-added coherent states and their generalizations was studied in [62,63]. "Semi-coherent" states and their generalizations were investigated in this connection in [64,65]. The binomial states and their generalizations were considered in [66,67]. The generation of high-order squeezing in multiphoton processes was studied in [68].
The concept of significant squeezing can be especially important for n 1, when the significant squeezing boundary Λ 2n is approximately n/4 times lower than the vacuum (Hong-Mandel) boundary Λ vac 2n . Figure 11 shows the ratio ψ 2 |x (2n) |ψ 2 / 0|x (2n) |0 ≡ f as function of parameter |β| for the two-Fock superposition |ψ 2 (3) with φ = π. Function f is given by Equation (21) with j = 1. The normalized significant squeezing levelΛ 2n = f min is calculated while using the exact Formula (22) with j = 2. The high-order σ-squeezing practically coincides with the simple high-order squeezing for n 1, 2 ∼ (2j)!2 n (n/2) −j 1 in this case. However, the boundaries of two kinds of squeezing are quite different if n is not very big.

Conclusions
Let us emphasize the main results of our study. We have analyzed the high-order moments ofx andp operators in the most general pure inhomogeneous Gaussian (squeezed) states: see, e.g., Formulas (12) and (14). We have found that the minimal product of the fourth order variances Π (4) σ for this specific family of quantum states is attained for all of the uncorrelated vacuum squeezed states: see Equation (15). However, this product can be diminished by more than 20% in the two-Fock superpositions |ψ 4 (see Equation (27)) and in the compass (four-photon) states [see Section 6.2].
We have shown that all of the Fock states are Robertson-Schrödinger intelligent states (RSIS) for the pair of operatorsx 2 andp 2 , although they are not ordinary intelligent states for this pair.
Additionally, we have discovered that the only Gaussian ordinary intelligent state for the pairx 2 ,p 2 is the correlated coherent state with the correlation coefficient r = ±1/ √ 2. We have analyzed in detail the uncertainty products of the 4th and 6th orders in the two-Fock and three-Fock superpositions. In this connection, we have found new explicit expressions for the mean values and matrix elements between arbitrary Fock states for arbitrary powers of operatorx (as well as forp): see Equations (A7) and (A11) in Appendix B. Using these expressions, we have obtained the asymptotic Formula (34) for the uncertainty product Π (2n) with n 1 in the two-Fock superpositions, clarifying the earlier results of paper [4].
We have introduced the new concept of weak and significant high-order squeezing. It was illustrated in the examples of Gaussian states and two-Fock superpositions. The distinguished role of the superposition of the vacuum and the fourth Fock state was demonstrated.
One of the most impressive results is that very small deformations of the vacuum state can result in big changes of such quantities as the high-order uncertainty product or high-order squeezing. This is an indication to the insufficiency of the concept of quantum "fidelity" (the scalar product between two states in the case of pure quantum states) as a measure of distinguishability between quantum states in certain physical situations. See, in this connection, e.g., earlier papers [69][70][71][72][73][74][75].
Several interesting questions and problems arise in connection with our study. For example, is it possible to obtain more precise limits for the high-order uncertainty products for several degrees of freedom, generalizing known results for the products of variances [2,29,[76][77][78][79][80][81][82]? Probably, applications to the problem of intensity-difference squeezing [83,84] can be found. Or, how to generalize our results to the case of mixed quantum states with the fixed value of quantum purity (see, e.g., studies [85,86] for the second-order uncertainty products)? It would be interesting to verify the new relations in experiments, for example, using the scheme of measuring high-order moments of canonical variables proposed in [87] (some experiments related to uncertainty relations for the second-order moments were performed in [88]; experimental verifications of the uncertainty relations for the angular momentum operator and other quantum systems in finite-dimensional Hilbert spaces were reported in [89][90][91][92][93][94]). Other interesting areas, where high-order uncertainty relations can be useful, include entanglement, steerability, and Bell nonlocality (see, e.g., studies [95][96][97]). However, all of these problems need separate investigations.

Conflicts of Interest:
The authors declare no conflict of interest.

Appendix A. Details of Calculations for the Gaussian States
To derive the first formula in Equation (12), we start from the expressions The second integral is well known: see, e.g., Equation 3.462.4 in [98]. It results in the formula where H n (z) is the Hermite polynomial, defined according to the formula H n (y) = (−1) n e y 2 d n dy n e −y 2 . (A2) Using the relations (11), we obtain from (A1) the first equality in Equation (12). Then, the second equality in (12) is the consequence of the fact that the Fourier transform of any Gaussian function is a Gaussian function again. A formula similar to (A1) for the powers of linear momentum has the form p n = i −n (a 2 + b 2 )/(4a) n/2 H n i(ad − bc) a (a 2 + b 2 ) .
In particular, the special cases of (A1) and (A3) for n = 1 and n = 2 yield the formulas x = c/a, σ x = (2a) −1 , p = (ad − bc)/a, σ p = (a 2 + b 2 )/(2a), confirming the known fact that all moments in the Gaussian states are determined by the first-order mean values and variances. We can write (A1) and (A3) also in terms of the complex numbers α = a + ib and β = c + id: x n = 2i Re(α) A simple formula for the mean value x npn can be obtained, if x = p = 0. Indeed, using the relation −α * = α − 2a and the definition (A2), one can see that x npn = i n √ a/π(α/2) n/2 ∞ −∞ x n e −ax 2 H n x √ α/2 . Then formula 7.385.2 from [98] yields Equation (14).