Uncertainty Relations for Coarse–Grained Measurements: An Overview

Uncertainty relations involving incompatible observables are one of the cornerstones of quantum mechanics. Aside from their fundamental significance, they play an important role in practical applications, such as detection of quantum correlations and security requirements in quantum cryptography. In continuous variable systems, the spectra of the relevant observables form a continuum and this necessitates the coarse graining of measurements. However, these coarse-grained observables do not necessarily obey the same uncertainty relations as the original ones, a fact that can lead to false results when considering applications. That is, one cannot naively replace the original observables in the uncertainty relation for the coarse-grained observables and expect consistent results. As such, several uncertainty relations that are specifically designed for coarse-grained observables have been developed. In recognition of the 90th anniversary of the seminal Heisenberg uncertainty relation, celebrated last year, and all the subsequent work since then, here we give a review of the state of the art of coarse-grained uncertainty relations in continuous variable quantum systems, as well as their applications to fundamental quantum physics and quantum information tasks. Our review is meant to be balanced in its content, since both theoretical considerations and experimental perspectives are put on an equal footing.


Introduction
The physics of classical waves distinguishes itself from that of a classical point particle in several ways. Waves are spread-out packets of energy moving through a medium, while a particle is localized and follows a well-defined trajectory. It was thus most surprising when it was discovered in the early 20th century that quantum objects, such as electrons and atoms, could exhibit behavior that at times was best described according to wave mechanics. Moreover, it was shown that either wave or particle behavior could be observed depending almost entirely upon how an observer chooses to measure the system. This complementarity of wave and particle behavior played a key role in the early debates concerning the validity of quantum theory [1], and has been linked to several interesting and fundamental phenomena of quantum physics [2][3][4][5]. Though several complementarity relations have been cast in quantitative forms [6,7], perhaps complementarity is most frequently observed in terms of quantum uncertainty relations. In words, uncertainty relations establish the fact that the intrinsic uncertainties associated to measurement outcomes of two incompatible observations of a quantum system can never both be arbitrarily small. We note that this type of behavior appears in classical wave mechanics, for example in the form of time-bandwidth uncertainty relations, which are quite important in communications and signal processing [8]. In contrast, there is no aspect of a classical physics that prohibits us from measuring all of the relevant properties of a classical point particle, at least in principle.
In addition to quantum fundamentals, quantum uncertainty relations play an important role in several interesting tasks associated to quantum information protocols, such as the detection of quantum correlations and the security of quantum cryptography [9]. In this paper, we focus on continuous variable (CV) quantum systems [10,11]. Though many interesting results have been found for discrete systems, they are outside the scope of this manuscript. We refer the interested reader to Reference [9], being a comprehensive unification and extension of two older reviews on entropic uncertainty relations, more focused on the physical [12] and information-theoretic [13] side respectively. However, since the coarse-grained scenario situates itself somehow in-between the discrete and continuous description, we make a short introduction to discrete entropic uncertainty relations before discussing their coarse-grained relatives.
In CV systems, one encounters a fundamental problem when performing measurements. That is, the eigenspectra of the corresponding observables are infinite dimensional, and can be continuous or discrete. Since any measurement device registers measurement outcomes with a finite precision and within a finite range of values, the experimental assessment of CV observables can be quite different from theory. Of course, one can consider a truncation of the relevant Hilbert space [14], as well as some type of binning or coarse graining of the measurement outcomes. This is similar to the idea of coarse graining that was discussed by Gibbs [15] and used by Paul and Tanya Ehrenfest [16,17] in the early 20th century to account for imprecise knowledge of dynamical variables in statistical mechanics [18]. Coarse graining has also appeared in the quantum mechanical context as an attempt to describe the quantum-to-classical transition, where the idea is that measurement imprecision could be responsible for the disappearance of quantum properties [19][20][21][22][23]. Though this is quite an intuitive notion, it was recently shown that one can always find an uncertainty relation that is satisfied non-trivially for any amount of coarse graining [24]. That is, quantum mechanical uncertainty is always present in this type of "classical" limit. This motivates the formulation of coarse-grained uncertainty relations.
In addition to the necessity of coarse graining, there could be practical advantages: for tasks such as entanglement detection, it might be interesting to perform as few measurements as possible, advocating the use of coarse-grained measurements. However, improper handling of coarse graining can result in false detections of entanglement [25,26], pseudo-violation of Bell's inequalities or the Tsirelson bound [27,28], and sacrifice security in quantum key distribution [29], for example. Thus, the proper formulation and application of uncertainty relations for coarse-grained observables is both interesting and necessary.
In the present contribution we review the current state of the art of uncertainty relations (URs) for coarse-grained observables in continuous-variable quantum systems. In Section 2 we review the concept of uncertainty of continuous variable (CV) quantum systems in more depth and introduce several prominent URs. In Section 3 we discuss the utility of CV URs in quantum physics and quantum information, in particular for identifying non-classical states and quantum correlations. Section 4 presents the problem of coarse graining of CVs in detail, and two coarse-graining models are provided. The current status of URs for these coarse-graining models is reviewed in Section 5, where we present a series of coarse-grained URs previously reported in the literature [12,24,[30][31][32]. In addition, we extend the validity of some of these URs to general linear combinations of canonical observables. Section 6 is devoted to the experimental investigation and application of coarse-grained URs in quantum physics and quantum information. Concluding remarks are provided in Section 7.

Uncertainty Relations
The history of uncertainty relations traces back to the early days of the formalization of quantum theory and begins with the celebrated work by Heisenberg in 1927 [33] (see [1] for an English version).
The work discussed what later became known as Heisenberg's uncertainty principle. The first mathematical formulation for this principle, in [33], essentially reads: where ∆x and ∆p are the uncertainties of the position and linear momentum of a particle, respectively, and h is the Planck constant. Although the existence of such a principle is ultimately due to the non-commutativity of the position and momentum observables, it took almost 80 years for all the physical meanings, scope and validity of this principle to be elucidated [34]. Distinct physical meanings emerge from different definitions for "uncertainty" of position or momentum, and in each case a proper multiplicative constant makes the lower bound sharp. All of these inequalities are known by the generic name of Uncertainty Relations, from the beginning of this review referred to as URs.
Even though the inception of the URs was made in the context of position and momentum of a particle, their existence can be extended to the "uncertainties" associated with any pair of non-commuting observables in discrete or continuous variable quantum systems. Thus, generically we can define the URs as inequalities that stem from the fact that the measured quantities involved are associated to non-commuting observables. Nowadays, we can say that it is clear that there are three conceptually distinct types of URs [34]: (i) URs associated with the statistics of the measurement results of non-commuting observables after preparing the system repeatedly in the same quantum state, or statistical URs for short, (ii) the error-disturbance URs, also known as noise-disturbance URs, for the relation of the imprecision in the measurement of one observable and the corresponding disturbance in the other, and, (iii) the joint measurement URs associated with the precision of the joint measurements of non-commuting observables. The error-disturbance URs has two main contributions: one in References [35][36][37] that present state-dependent error-disturbance URs and the other in References [38][39][40] that argue for a state-independent characterisation of the overall performance of measuring devices as a measure of uncertainty that satisfies an UR of the form given in Equation (1). There was a certain controversy involving these two contributions and we recommend the work [41] that discusses the limitations of the state-dependent error-disturbance URs. The development of joint measurement URs has an early contribution in Reference [42] and further developments were given in References [38,[43][44][45][46][47][48].
The statistical URs are also referred to in the literature as preparation URs. This is because it is impossible to prepare a quantum system in a state for which two non-commuting observables have sharply defined values. However, here we prefer to call them statistical URs, as they express the limits to the amount of information that can be obtained about incompatible observables of a quantum system when it is repeatedly measured after being prepared in the same initial state in each round of the measurement process. We emphasize that there is not any attempt to measure the two non-commuting observables simultaneously. In each round of the measurement process only one observable is measured, the choice of which could be made randomly. In this sense the "uncertainties" contained in the statistical URs are of the statistical type: the more certain the sequence of outcomes of one observable is in a given state, then the more uncertain is the sequence of outcomes of the other non-commuting observable(s) considered.
This review focuses on statistical URs that are valid for coarse-grained measurements in continuous variable quantum systems, although a similar approach can be made for the other two types of URs mentioned above. There are two types of quantum mechanical degrees of freedom: the ones that can be described by a Hilbert space of quantum states with finite dimension and the others in which it has infinite dimension. In particular, we are interested only in continuous variable (CV) systems where the Hilbert space, H, of pure states, |ψ , has an infinite dimension. The CV systems that we consider consist of a finite set of n bosonic modes, sometimes called "qumodes" [10], so that H := H 1 ⊗ . . . ⊗ H n . Each mode is described by a pair of canonically conjugate operators,x j andp j , such that [x j ,p k ] = ih1δ jk . (2) jk . Therefore, the separable Hilbert space of each mode, H j , has a enumerable basis {|n j } n j =1,...,∞ consisting of eigenstates of the number operator, viz. n j |n j = n j |n j , evidencing the infinite dimensionality of the Hilbert space of the quantum states. In the case of mixed states we use density operators represented by greek letters with a hat, i.e.,ρ,σ etc. Important examples of CV systems are the motional degrees of freedom of atoms, ions and molecules, wherex j andp j are the components of the position and linear momentum of the particles (in this caseh in Equation (2) is the usual reduced Planck constant, i.e.,h = h/2π); the quadrature modes of the quantized electromagnetic field wherex j andp j are canonically conjugate quadratures (in this caseh in Equation (2) is justh = 1 ) [10]; and the transverse spatial degrees of freedom of single photons propagating in the paraxial approximation (in this caseh in Equation (2) ish = λ/2π where λ is the photon's wave length [49]).
In what follows, we summarize the principal statistical URs in CV systems that have been generalised to coarse-grained measurements. The corresponding coarse-grained URs will be presented in Section 5.

Heisenberg (or Variance) Uncertainty Relation
Let us consider two operators: where T means transposition and we define the 2n-dimensional vector of operators, as well as the arbitrary real vectors, d = (a, a ) T = (a 1 , . . . , a n , a 1 , . . . , a n ) T and The commutation relation ofû andv is where J is the 2n × 2n-dimensional matrix of the symplectic norm [50]: and the n × n matrices in the blocks are the identity matrix I and the null matrix O. In this review, matrices of an arbitrary shape not treated as quantum-mechanical operators are denoted in bold and without a hat. The parameter γ in definition Equation (6) is a scalar that in some sense quantifies the non-commutativity ofû andv. Commutation relations such as Equation (6) are called Canonical Commutation Relations (CCR) (sometimes the name CCR is used in the case when γ = 1, however, ashγ can be interpreted as an effective Planck constant, so the name CCR here is well justified). However, a CCR between two operatorsû andv does not guarantee that they are necessarily Canonically Conjugate Operators (CCOs). For this to be true we additionally need that the eigenvectors ofû andv must be connected by a Fourier Transform. In such a case we callû andv CCOs (also note that when two operators like the ones defined in Equation (3) have their eigenstates connected by a Fourier Transform, they necessary satisfy a commutation relation like in Equation (6), as can be easily shown. However the converse is not true. Take for example the single mode operatorsû =x and v =x +p, which satisfy [û,v] = [x,p] = ih but are not a Fourier pair).
Every pair of operators,û andv, that obey a CCR also satisfies the statistical UR: where σ 2 P u := û 2 − û 2 , and σ 2 are the variances of the marginal probability distribution functions (pdf): where we have defined . . . := Tr(. . .ρ), withρ being an arbitrary n−mode quantum state. We call the UR in Equation (8) the Heisenberg UR, or variance-product UR. For one mode CCOs, such asû =x andv =p (therefore γ = 1), the Heisenberg UR in Equation (8) was first proved by Kennard in 1927 [51], inspired by the inequality in Equation (1) of Heisenberg's seminal paper of the same year [33]. Later, it was also proved by Weyl in 1928 [52]. In 1929 Robertson [53] extended the Heisenberg UR for any pair of Hermitian operatorsÂ andB: This result extends the Heisenberg UR in Equation (8) toû andv that are not CCOs. For every variance-product UR in Equation (12) there is an associated linear UR: In fact, this UR is a consequence of Equation (12) and the trivial inequality (σ P A − σ P B ) 2 ≥ 0, so that where it also follows that the linear UR is weaker than the variance product UR. In 1930 Schrödinger [54] improved the lower bound in Equation (12), so the new stronger UR reads: where {· · · , · · · } is the anti-commutator. One interesting property of the Heisenberg UR in Equation (8) is that the lower bound is independent of the quantum stateρ under consideration. Another property is that it can be seen as a bona fide condition on the covariance matrix of an n−mode quantum stateρ, viz the matrix of second moments of the CCOs, contained in the vectorx, of the stateρ: Indeed, in [55,56] it was shown that the bona fide condition on the covariance matrix V of a quantum stateρ is, where the inequality means that the matrix on the left hand side is positive semi-definite, viz. all of its eigenvalues are greater or equal to zero. Applying the inequality in Equation (15) to the canonical conjugate operatorsx andp, we have, For one mode systems, this inequality is equivalent to the bona fide condition in Equation (17). However, for multimode systems it is not enough. For multimode systems, a way to verify the bona fide of the covariance matrix was given in [57,58]. It was shown that testing the condition in Equation (17) is equivalent to verify the linear UR in Equation (13) for all the operators,û andv, defined in Equation (3). Therefore, using Equation (14) we can write the series of implications: Thus, it is enough to verify the violation of the Heisenberg UR for some pair of operatorsû and v to confirm that the bona fide condition on the covariance matrix of some n−mode operatorρ is not satisfied.

Entropic URs
The use of entropy functions to quantify uncertainty of a probabilistic variable dates back to the early work of Shannon [59]. Since then, several different entropy functions have been defined, with distinct relations to meaningful characteristics of the probability distributions considered. A number of these entropy functions have found use in quantum mechanics and, in particular, in QIT [9]. Here we outline the application of these functions to uncertainty relations between non-commuting observables.

Shannon-entropy UR
The UR based on the differential Shannon entropy for operators defined in Equation (3) is: where P u and P v are the marginal pdf defined in Equation (10) and the differential Shannon entropy of a pdf, P, is defined as [60]: Forû andv as CCOs, this uncertainty relation was first proved in 1975 by Bialynicki-Birula and Mycielski [61]. In their derivation the authors used the L p -L q norm inequality for the Fourier transform operator obtained by Beckner [62]. Please note that in the literature this inequality is sometimes referred to as the Babenko-Beckner inequality (Equation 1.104 from [12] provides an extension of this inequality to the case of arbitrary mixed states, using two variants of the Minkowski inequality), because Babenko [63] had proved it before Beckner, but only for certain combinations of (p, q) parameters. For the sake of completeness, we should also mention that Hirschman [64] had derived a weaker version of (20) with the constant eπ inside the logarithm replaced by 2π. The extension of the validity for operatorsû andv that are not CCOs was provided very recently in References [58,65].
The Shannon-entropy UR is in general stronger than the Heisenberg UR as the former implies the latter. This can be seen by using the inequality for a pdf P [60]: where σ 2 is the variance of P. Therefore, we can write the chain of inequalities: that compress the URs in Equations (8) and (20). It is clear from Equation (23) that the verification of the Shannon-entropy UR for any pair of the operators in Equation (3) is enough to guarantee the bona fide condition in Equation (17) [58]. When the quantum stateρ is Gaussian, viz when the Wigner function ofρ is a multivariate Gaussian probability distribution [11], the marginal pdfs, P u and P v , are also Gaussians. Remembering that the differential Shannon entropy of a Gaussian pdf P, with variance σ 2 P , is h[P] = (1/2) ln 2πeσ 2 P [60], we can see that Gaussian states saturate the first inequality in Equation (23). Therefore, for Gaussian states the Heisenberg UR and the Shannon-entropy UR are completely equivalent. As we will see in Section 5 this is not the case for the coarse-grained versions of these URs.

Rényi-Entropy URs
The UR based on the differential Rényi entropy for the operators defined in Equation (3) that are CCOs is given by the inequality: where 1/α + 1/β = 2 with 1/2 ≤ α ≤ 1 and γ = 1 since we deal with CCO operators. As before, P u and P v are the marginal pdfs defined in Equation (10) and the differential Rényi entropy of order α relevant for an arbitrary pdf, P, is defined as [60]: The Rényi-entropy UR was proved recently (in 2006) by Bialynicki-Birula [31] (see also [12]) again with the help of the powerful mathematical tools developed in [62]. Please note that in the limit α → 1 we also have β → 1, and consequently α 1 (2−2α) β 1 (2−2β) → 1/e. Therefore, in the limit α → 1 we have , so the expression in Equation (24) reduces to the Shannon-entropy UR in Equation (20) for γ = 1. As far as we know, in contrast to the Shannon-entropy UR, the extension of the Rényi-entropy UR to the general case of operators that are not necessarily CCOs is still a challenge for the future. A first attempt in this direction was provided in Reference [65], where the authors show that the Rényi UR in Equation (24) is still valid when the eigenvectors ofû andv are connected by a Fractional Fourier Transform [8], which corresponds to rotation in phase space.
All of the URs mentioned in this section (this is a general pattern though) can be cast in a general form where F is an uncertainty functional [left hand side of inequalities Equations (8), (20) and (24) for example] and f represents its respective lower bounds. In particular, we do not pay much attention here to the Tsallis entropy and URs associated with it. Again such URs can be cast in the general form stated above and their derivation is usually very similar in spirit to the case of the Rényi entropy. In Section 3 we will summarise the relevance of the statistical UR in general and in particular the URs in Equations (8), (20) and (24). In Section 5 we will present versions of the Heisenberg, Shannon-entropy and Rényi-entropy URs for coarse-grained measurements.

Utility of Uncertainty Relations in Quantum Physics
Uncertainty relations can be applied in several useful and interesting ways. First, they provide a way to test if experimental results are consistent with quantum mechanics, since data from the measurement of incompatible observables must verify any valid quantum UR. This is particularly helpful in identifying systematic errors in the measurement process, in testing the experimental reconstruction of density matrices, phase-space distributions (quantum state tomography), as well as the covariance matrix [66], or any other set of moments of the CCOs of the modes.
URs can also be used to characterize non-classical states of light, such as squeezed states [67]. In this case observation of the variance σ 2 P u ≤h/4 whereû is a phase-space quadrature in Equation (3), indicates noise fluctuations in this quadrature that are smaller than the vacuum state. As a consequence of the Heisenberg UR, the noise fluctuations in the conjugate quadrature must be larger or equal tō h/4σ 2 P u . In a similar fashion, in Reference [68] it was shown that violation of one out of an infinite hierarchy of inequalities involving normally ordered quadrature moments is sufficient to demonstrate non-classicality. We note that σ 2 P u ≤h/4 corresponds to the lowest-order inequality of this set. Related techniques have been developed based on the quantum version of Bochner's theorem for the existence of a positive semi-definite characteristic function [69,70]. Both of these methods have been used experimentally in Reference [71]. More recently, these two techniques were unified into a single criteria involving derivatives of the characteristic function [72], and put to test on a squeezed vacuum state.
To our knowledge, the first application of URs to identify quantum correlations was described in Reference [73], in which the authors proposed a Heisenberg-like UR, similar to that in Equation (8), to identify non-classical correlations between both the phases and intensities of the fields produced by a non-degenerate parametric oscillator. It was shown by M. Reid [74] that these measurements provide a method to demonstrate correlations for which the seminal Einstein-Podolsky-Rosen (EPR) argument [75] is valid. An experiment using this UR-based method to demonstrate EPR-correlations between light fields was realized shortly therafter [76]. It was later shown by Wiseman et al. [77,78] that the Reid EPR-criterion was indeed a method to identify quantum states that violate a "local hidden state" model of correlations. This type of correlation has been called "EPR-steering", or just "steering" [79], as this was the terminology used by Schrödinger when he discussed EPR correlations in 1935 [80]. Since 2007, EPR-steering has been understood to make up part of a hierarchy of quantum correlations, situated between entanglement [81,82] and Bell non-locality [83]. In addition to methods utilizing variance-based URs [84], entropic URs, such as those in Section 2.2, can be used to identify EPR-steering [85,86] and to quantify high-dimensional entanglement [87,88]. Some of these URs can be used to test security in continuous variable quantum cryptography [89,90], and it has been shown that violation of entropic EPR-steering criteria are directly related to the secret key rate in one-sided device independent cryptography [91]. We also highlight techniques based on a matrix-of-moments approach [92]. Continuous-variable EPR-steering has been observed in intense fields [76,93,94] as well as photon pairs [85,[95][96][97].
Perhaps one of the most important tasks in quantum information is identifying quantum entanglement. In this respect, URs have also found widespread use in simple and experimentally friendly entanglement detection methods, as we will now describe. Several early entanglement criteria for bipartite CV systems were developed using URs [98][99][100][101]. A particularly convenient method to construct entanglement criteria is to use the Peres-Horedecki positive partial transposition argument [102,103] (PPT), and apply it to uncertainty relations [82,[104][105][106][107]. The PPT argument is as follows. A bipartite separable stateσ 12 can be written as [108] whereρ 1i andρ 2i are bona fide density operators of subsystems 1 and 2, respectively. The transpose of the stateρ 2i , here denotedρ T 2i , is still a positive operator, since full transposition preserves the eigenspectrum. Thus, partial transposition (with respect to second subsystem) ofσ 12 gives the valid quantum state:σ On the other hand, partial transposition of an entangled stateˆ 12 , which cannot be written in the form (27), can lead to a non-physical density matrix since partial transposition may not preserve the positivity of the eigenspectrum. Thus, one can identify entanglement in a bipartite density operator by calculating the partial transposition and searching for negative eigenvalues, and even quantify the amount of entanglement via the negativity [109]. However, applications of this method in experiments requires quantum state tomography and reconstruction of the density operator, which involves a large number of measurements. A more experimentally friendly method to identify entanglement is to evaluate an UR applied to the partial transposition ofˆ 12 , which we describe in the next paragraph. The PPT-argument is only a sufficient entanglement criteria in a general bipartition of m × (n − m) modes, but is necessary and sufficient in the particular case of bipartitions of 1 × (n − 1) modes in CV Gaussian states [10,110]. Thus, there are no Gaussian states which are PPT entangled states in bipartitions of the form 1 × (n − 1). However, there do exist entangled CV Gaussian states that are PPT in general bipartitions of the type m × (n − m). These are called bound entangled states [111]. In Gaussian states, this set of bound entangled states coincides with the set of all states whose entanglement in a bipartition m × (n − m) cannot be distilled using local operations and classical communication [112][113][114]. However, to our knowledge, for non-Gaussian states it is conjectured that the set of bound entangled states in a given bipartition is only a sub-set of the set of undistillable states in that bipartition.
For continuous variables, Simon showed that transposition is equivalent to a momentum reflection, taking the single mode Wigner phase-space distribution W (x, p) −→ W T 2 (x, p) = W (x, Tp) [57], where T is a diagonal matrix whose elements are +1 for non-transposed modes, and −1 for the transposed ones. Thus, evaluating the "transposed" Wigner function is the same as evaluating the original Wigner function with a sign change in the reflected p variables.
For simplicity, we consider now the particular example of global operators of a bipartite state: andv ± =v 1 ±v 2 .
We note that operators with the same sign satisfy the commutation relations [û ± ,v ± ] = 2ihγ , so that these non-commuting operators after being an input to the uncertainty functionals fulfill the UR of the aforementioned form [note the factor of 2 in the argument of f (·)] Using the transformation of the Wigner function under partial transposition described above, one can evaluate the uncertainty functional of the partially transposed stateˆ T 2 12 via measurements on the actual state 12 using the relation which can be lower than f (2h|γ|) since the operators with different signs do commute. This possibility, when experimentally confirmed, indicates that T 2 12 is not a bona fide density operator, and thus the bipartite quantum state 12 is entangled.
Building on this general reasoning (PPT argument applied to an UR) several entanglement criteria have been developed. A comprehensive list of the criteria contains those based on the variances [115,116] and higher-order moments [117,118], Shannon entropy [105], Rényi entropy [106], characteristic function [119] as well as the triple product variance relation [120]. Particularly powerful is the formalism developed by Shchukin and Vogel, which provides an infinite set of inequalities involving moments of the bipartite state [121], such that violation of a single inequality indicates entanglement. We note that some of these criteria can be applicable to any non-commuting global operators. Uncertainty-based approaches (using the PPT method directly or not) have been developed for multipartite systems [122,123], and a general framework to construct entanglement criteria for multipartite systems based on the "PPT+UR" interrelation was presented in Reference [107]. The Shchukin-Vogel hierarchy of moment inequalities has also been applied to the multipartite case [124].
The PPT+UR approach has been used to identify continous variable entanglement experimentally in several systems, including entangled fields from parametric oscillators and amplifiers [94,125,126] as well as spatially entangled photon pairs produced from parametric down conversion [96,120,127], and time/frequency entangled photon pairs [128,129]. A higher-order inequality in the Shchukin-Vogel criteria [121] has been used to observe genuine non-Gaussian entanglement [130].

Realistic Coarse-Grained Measurements of Continuous Distributions
Coarse graining of observables with continuous spectra is a consequence of any realistic measurement process. In the laboratory, an experimentalist is given the task of designing projective measurements in order to recover information about probability densities of a continuous variable quantum system. Naturally, only partial information about the underlying continuous structure of the infinite-dimensional physical system is retrieved in a laboratory experiment. Whichever measurement design is chosen, the experimentalist is faced with two main difficulties, namely the finite detector range and finite measurement resolution, related to the size of the total region of possible detection events and the precision in which events are registered, respectively. The detector range problem [25,29] results from the finite amount of resource available to the experimentalist. For instance, consider a position discriminator based on a multi-element detector array. The array has a spatial reach (in a single spatial dimension) that increases linearly with the number of detectors. In a similar fashion, the sampling time of a single element detector used in raster scanning mode increases linearly with the chosen detection range. Continuous variables such as the position are also inevitably affected by the inherent finite resolution of the measurement apparatus [32], such as the size of each individual detector in the array, or the pixel size of a camera. Altogether, the finite detector range and measurement resolution restrict the capability to probe the detection position, limiting the experimentalist to a coarse-grained sample of the underlying CV degree of freedom.
The constraints imposed by the finite spatial reach and resolution of the measurement apparatus are then important features that must be considered in the experiment design. Ideally, the experimentalist would chose measurement settings producing the finest coarse-grained sample possible. As a trade-off, the increased resolution entails the sampling of a greater number of pixels (if the range of detection is preserved), increasing the amount of resources used in data acquisition and analysis. The compromise between the used resource and chosen resolution depends on the specific design and measurement technique. A single raster scanning detector is inherently inefficient and leads to acquisition times that grow with the number of scanned outcomes. On the other hand, the acquisition time is dramatically reduced by the use of multi-element detector arrays [131][132][133][134]. Other techniques such as position-to-time multiplexing [135,136] allow the sampling of multiple position outcomes with single element detectors, but at the expense of an increased dead-time between consecutive detections. We have exemplified the finite detector range and finite measurement resolution problems in terms of a detector that registers the position of a particle. However, similar considerations are valid for any detection system that registers a digitalized value of a continuous physical parameter.
Under constraints of resource utilisation-such as the number of detectors and/or sampling time-the experimentalist needs to set the number of possible detection outcomes for their coarse-grained measurements. Therefore, a natural question that arises regards the coarse-graining design allowing the extraction of the desired information. Naively, one might think that usual quantum mechanical features learnt from physics textbooks would be directly observable from the coarse-grained distributions obtained in the laboratory. The most prominent counter-example is the experimental observation of the Heisenberg UR in Equation (8). As shown in Reference [32], coarse-grained distributions of conjugate continuous variables do not necessarily satisfy the well known UR valid for continuous distributions. In order to accurately inspect the uncertainty product of the measured distributions in accordance with the Heisenberg UR, the latter must be modified to account for the detection resolution of the measurement apparatus. Another important quantum mechanical feature that one usually fails to observe from standard coarse-grained distributions is the mutual unbiasedness [137] relation between measurement outcomes of the incompatible observables. That is, eigenstates of-say-the coarse-grained position operator do not necessarily present a uniform distribution of outcomes for coarse-grained momentum measurements. In addition, some authors [138][139][140][141] have demonstrated that one can define functions of incompatible observables that indeed commute. Interestingly, it was shown in Reference [142] that one can indeed enjoy full quantum mechanical unbiasedness using a specific periodic coarse-graining design rather than the standard one. Other practical issues regarding false positives in entanglement detection [26,29] and cryptographic security [25,29] must also be reconsidered when one deals with realistic coarse-grained distributions.
In this section, we will introduce the projective measurement operators both for the standard and the periodic models of coarse graining. Practical features such as measurement resolution, detector range and positioning degrees of freedom in the measurement design will be discussed. We will also briefly discuss relations of mutual unbiasedness between coarse-grained measurement outcomes in domains of incompatibles observables. A detailed discussion of uncertainty relations for coarse-grained distributions will be presented in the next section.

Coarse-Graining Models
A laboratory experiment necessarily yields a discrete, finite set of measurement outcomes of any observable in any physical system. This is also the case for an experiment probing a continuous degree of freedom,û, for which measurement outcomes {u k } labeled by the discrete integer index k ∈ Z relate to the underlying continuous real variable u ∈ R corresponding to the eigenspectra ofû. In the most general scenario, a coarse-graining model is obtained from an arbitrary partition of the set of real numbers R, in intervals R k with u k ∈ R k . The orthogonality of the measurement outcomes requires the subsets to be mutually disjoint: R k ∩ R k = ∅, ∀ k = k . Even though the continuous variable can be formally discretised into an infinite number of outcomes (with k an unbounded integer), the experiment can only probe a finite range of the continuous variable. Thus, the detection range, R range , can be formally defined by the union of the disjoint subsets associated with the probed outcomes: This relation limits the set of possible values of k to a finite subset of integers Z k ⊂ Z. Due to the finite range, R range , of the measurement process it is important to secure under reasonable experimental conditions that the underlying probability density is supported within the chosen range of detection [25,29]. Mathematically, a faithful coarse-grained measurement design should ensure that where P u is the marginal pdf defined in Equation (10).
The probability p (u) k that the outcome u k is produced writes as an integral of the marginal probability density, P u , for the continuous variable: where the integration is performed in the interval R k . Due to the faithful coarse-grained condition in Equation (34) we have We can define projective operators associated with the coarse-grained measurements: so that the probabilities (35) can be written as with P u (u) = u|ρ|u . In order to study mutual unbiasedness and uncertainty relations, we shall later in this and the following sections define coarse-grained operators like those in Equation (37) for conjugate variables of the quantum state, such as the position and the linear momentum of a quantum particle.

Standard Coarse Graining
The standard model of coarse graining describes, for example, the typical projective measurements performed with an array of adjacent, rectangular detectors. A conventional example of such an apparatus is the image sensor of a digital camera, for which the pixel size stands for the detection resolution whereas the length of the full sensor embodies the range of detection. In the current analysis, we shall consider a linear detector array along a single spatial dimension rather than the two-dimensional area of a typical image sensor, as illustrated in Figure 1. The coarse-graining interval representing the detection window of the k-th pixel of the linear array is then: where ∆ is the detector or pixel size-also commonly referred to as the coarse-graining width or the bin width. Using the definition Equation (39), the discretised outcomes u k represent the u value of the center of the corresponding bin: Multi-element detector array Detection range The parameter u cen sets the position of the central bin of the array, whose outcome label is k = 0, yielding u 0 = u cen . To illustrate the effect of the coarse-graining design on measured distributions, we plot in Figure 2 coarse-grained distributions (blue bars) obtained using 3 different resolutions: ∆ = 2 (left colum), ∆ = 1 (central column) and ∆ = 1/2 (right column). For each resolution, we plot two distinct distributions obtained using u cen = 0 (top row) and u cen = ∆/2 (bottom row). In other words, the coarse-graining bins of the distributions plotted at the bottom part of the figure are displaced by half a "pixel" in relation to the distributions at the top. Clearly, the distribution obtained using a fixed resolution is not unique, but the effect of small displacements (smaller than the bin width) gets less important as the resolution is increased. For comparison, the generating continuous distribution is plotted in red. We shall now use this model for standard coarse graining to explicitly define the discretised counterparts of the position and momentum operators given in Equation (3).
where the projectorĈ k is defined in Equation (37) having an equivalent definition forv measurements), and we used ∆ (δ) as the detection resolution forû (v) measurements. According to the definition in Equation (35), as a result of the the coarse-grained measurement ofû andv we obtain the discrete probabilities, p (u) ∆,k and p (v) δ,l .The discrete variances associate with these discrete probabilities are: where we define the set of discrete probabilities: One can see from the definitions (42) that if the bin widths ∆ and δ are such that p δ,l are sufficiently close to unity for some value of k and l, we have σ 2 Thus, naive application of any of the variance-based URs given in Section 2.1 would indicate a false violation of a UR. It has been shown in Reference [32] that the same argument applies to discretized versions of entropic URs, such as those of Section 2.2. Thus, proper treatment of standard coarse-grained measurements is essential in order to take advantage of the practical application of URs in QIT and quantum physics in general. In Section 5 we show how this can be done.

Periodic Coarse Graining
A distinct model of coarse graining discussed in the literature [142,143] is refereed to as periodic coarse graining (PCG). In this model, the partition of the whole set of real numbers R is performed in a periodic manner, leading to a finite number d of subsets R k , with k = 0, · · · , d − 1. The resulting discretization utilizes the index k as a direct label for the detection outcomes, in a similar fashion to what is usually defined for finite-dimensional quantum systems. The subsets R k are defined as [142]: where s u plays the role of a bin width similar to the resolution ∆ used for the standard coarse graining. In the definition Equation (44), bins of size s u are arranged periodically with the parameter T u representing the period, as illustrated in Figure 3 for the particular design using d = T u /s u = 5 detection outcomes. It is important to notice that this coarse graining design do not distinguish detections in distinct bins associated with the same detection outcome k (ranging from 0 to 4 in Figure 3). For example, a detection within any bin colored in red in (44) would lead to the same detection outcome k = 1. An interesting feature of the PCG model is that the number of detection outcomes is utterly adjustable by the choice of the parameters T u and s u , regardless of the chosen detection range. For instance, doubling the range of detection allows one to design PCG measurement using twice as much periods in its design, while maintaining the same number d = T u /s u of detection outcomes. As with the standard model, the reference coordinate u cen sets the center of the detection range also for the PCG design. Using the subset definition given in Equation (44), we can explicitly write the projector operators, Equation (37), for the PCG model aŝ where we extend the sum in n over Z without loss of generality, assuming that Equation (34) is satisfied. Analogously, we also define the PCG projective operators over the conjugate variable v: where we define s v and T v as the bin width and periodicity used in the PCG measurements of v.

Mutual Unbiasedness in Coarse-Grained Measurements
If a quantum system , with finite dimension, is described as an eigenstate of a given observable, the measurement outcomes of complementary observables are completely unbiased: each one of them occurring with equal probability, 1/d, where d is the dimension of the quantum system's Hilbert space. This unbiasedness relation is an important feature of quantum mechanics with no classical counterpart, and is usually cast in terms of the basis vectors constituting the eigenstates of two (or more) complementary observables. To be more precise, two orthonormal bases {|a k } and {|b l } are said to be mutually unbiased if and only if | a k |b l | 2 = 1/d for all k, l = 0, · · · , d − 1 [137]. The observation of unbiased measurement outcomes is customary in experiments with finite dimensional quantum systems. Not only routine, measurements in mutually unbiased bases (MUB) constitute a key procedure in several quantum information processing tasks, such as verification of cryptographic security [9], certification of quantum randomness [144], detection of quantum correlations [145][146][147] and tomographic reconstruction of quantum states [148,149].
Mutual unbiasedness is also extendable to continuous variables quantum systems [150], for which bases {|u } and {|v } such [û,v] = ihγ, satisfy | u|v | 2 = 1/(2πh|γ|), i.e., the overlap of the basis vectors |u and |v is independent (no bias) of their eigenvalues, u and v (note, however, that even thoughû andv are mutually unbiased observables, this does not imply that they are complementary, as would be the case for operators in a discrete quantum system [151]. In continuous variable quantum mechanics, mutual unbiasedness does not imply thatû andv are maximally incompatible [152]. In this case, complementary observables are typically defined as CCOs, that is, forming a Fourier transform pair). For CV systems, nevertheless, this relation is rather a theoretical definition than an experimentally observable fact, since the experimentalist has neither the capability to prepare nor to measure the (infinitely squeezed) eigenstates of theû andv. Instead, both the preparation and measurement procedures are limited to the finite resolution of the experimental apparatus. As discussed previously in this section, measurements of a CV degree of freedom render discretized, coarse-grained outcomes whose probabilities, Equation (35), are provided by a coarse-graining model described by the projective operators given in Equation (37). These coarse-grained probabilities obtained experimentally do not in general preserve the mutual unbiasedness complied by the underlying continuous variables.
To elaborate the issue, let us consider sets of projectors {Ĉ l } defining coarse-graining measurements in the complementary domains u and v of a continuous variable quantum systemρ. We assume measurement designs providing a number d of outcomes in each domain. In this scenario, the requirement for mutual unbiasedness is thus that the coarse-grained probabilities for measurements of one variable are evenly spread between all discretized outcomes whenever the quantum state is localized with respect to the coarse graining applied to its conjugate variable (and vice-versa). The subtlety in this requirement is the (infinite) degeneracy of normalizable quantum states that can be localized with respect to the chosen coarse graining. To emphasize this degeneracy, we refer to the outcome probabilities, Equation (35), with explicit dependency on the quantum state in order to mathematically phrase the condition for mutual unbiasedness in coarse-grained CV: the outcomes l } are mutually unbiased if for all quantum statesρ and k 0 , l 0 = 0, · · · , d − 1 we have [142]: p where, again, we stress that p l ), as in Equation (35). Having formulated the conditions for mutual unbiasedness, Equation (47), it is easy to perceive that the adjacent, rectangular subsets defining the standard coarse graining [Equation (39)] will not lead to unbiased measurement outcomes. Any CV distribution localized in a single coarse-graining bin (for example in the u variable) generates a probability density that decays in the Fourier domain (the v variable) along the adjacent bins within the detection range. This decay generates a non constante coarse-grained distribution that, by definition, is biased. Furthermore, the number d of detection outcomes in the standard design depends directly on the selected detection range, as well as on the chosen resolution. As a consequence, even though a particular localized distribution could lead to approximately unbiased coarse-grained outcomes in the Fourier domain, an extended detection range would increase the number of outcomes, thus spoiling the unbiasedness.
It is thus evident that in order to retrieve unbiased outcomes from coarse-grained measurement, a more contrived coarse-graining design is needed. As it turns out, it was shown in Reference [142] that the PCG design exactly fulfils the requirements for unbiased measurements of finite cardinality stated in Equation (47). A relation between the periodicities T u and T v used in the PCG of the conjugate variables u and v was analytically derived as a single condition for unbiased coarse-grained measurements: The unbiasedness condition stated in Equation (48) establishes infinite possibilities for the pair of periodicities T u and T v that can be used to design the mutually unbiased pair of PCG measurements defined in Equations (45) and (46), respectively. For instance, the simplest and most important case is the condition with m = 1, since it is valid for all d and provides the best trade-off between experimentally accessible periodicities: T u T v = (2πh)d. Conditions with m > 1 are also possible but are not general since they depend on the chosen number of outcomes d [142]. For example, for d = 4, valid conditions are found using m(mod d) = 1, 3 whereas for d = 5, valid conditions are found using m(mod d) = 1, 2, 3, 4. Importantly, the case with m(mod d) = 0 is always excluded, since in this case the PCG projectors describe commuting sets, [138][139][140]. In other words a joint eigenstate of the productΠ l existis for all k and l whenever T u T v = 2πh/c with c ∈ N [153]. It is also interesting to note that using the periodicity definition from the PCG design (T = ds), it is possible to write the unbiasedness condition given in Equation (48) in alternative, equivalent ways: Finally, in Reference [143] these results were generalized for PCG measurements applied to an arbitrary pair of phase space variables other than the conjugate pair formed by position and momentum. What is more, a triple of unbiased PCG measurements was also shown to exist for rotated phase space variables, along the same lines as the demonstration of a MUB triple in the continuous regime done in Reference [150]. Experimental demonstrations of unbiased PCG measurements were also carried out in References [142,143], both of them utilizing the transverse spatial variables of a paraxial light field.

UR for Coarse-Grained Observables
A kind of a paradigm shift in the theory of uncertainty relations was brought by the observation that everything can be efficiently characterized solely by means of probability distributions. As a result, tools known from information theory, such as information entropy, Fisher information and other measures, came into play. Additionally, the notion of uncertainty for discrete systems could better be captured that way. Since products of variances calculated for observables such as the spin are bounded in a state-dependent manner (so that the ultimate lower bound typically assumes the trivial value of 0), information entropies provide an attractive alternative [154]. Written already in the Rényi form, the above equation is a discrete counterpart of Equation (25), which corresponds to the discrete counterpart of Equation (21) when α = 1.
In the finite-dimensional case given by an arbitrary stateρ acting on a d-dimensional Hilbert space H, and a pair of non-degenerate, non-commuting observables,Â andB, one usually defines the probabilities associated to projective measurements: where |a i and b j , i, j = 1, . . . , d are the eigenstates of the operators associated with both observables. Disctrete entropic URs for the above probability distributions are of the general form with U ∈ U (d) being a unitary matrix with matrix elements U ij = a i b j . We denote P A substantially more renowned Maassen-Uffink (MU) bound [155] derived in 1988, is B MU αβ = − ln c 1 . This bound is however valid only for the conjugate parameters 1/α + 1/β = 2. Very recently, a plethora of new results [41,[156][157][158][159][160][161][162][163] improving the celebrated MU bound has been obtained. In particular, an approach based on the notion of majorization (suitable from the perspective of resource theories and quantum thermodynamics [164]) provides a significant qualitative novelty [156,157,159,163], which will also be touched upon in this section.
In this review we are concerned with the case in which continuous probability distributions P u (u) and P v (v) are replaced (viz. they were measured this way) by their discrete counterparts (k, l ∈ Z). According to the discussion in Section 4 we can use the definitions in Equations (35) and (39), and the condition in Equation (33), to write the discrete probabilities: with k ∈ Z k ⊂ Z. In the following we describe a series of URs for these discrete probabilities that are known as coarse-grained URs, derived in [24,[30][31][32]. These are the coarse-grained counterpart of the Heisenberg, Shannon entropy and Rényi entropy URs in Equations (8), (20) and (24) respectively. Here, we will closely follow the treatment in [24,32]; however, before we start we give a short historical overview and discuss a path towards extensions going beyond CCOs. The idea that generic quantum uncertainty could be quantified by the sum of Shannon entropies evaluated for discretized position and momentum probability distributions for the first time appeared in the contribution by Partovi [165]. He also derived the first coarse-grained UR which in the form is reminiscent to the Deutsch bound for finite-dimensional systems [154] (please note that both papers [154,165] have been published in 1983; however, Partovi in his first sentence refers to a "recent letter" by Deutsch). Both bounds [154,165] were obtained by means of a direct optimization, independently applied to every logarithmic contribution. Symmetry in developments of the URs for finite-dimensional and coarse-grained systems happened to be much deeper as the second coarse-grained result, by Bialynicki-Birula [30], is a counterpart of the MU bound [155]. The former result is an application of the continuous variant of the Shannon entropy UR (so the L p -L q norm inequality by Beckner [62]) supported by the Jensen inequality for convex functions, while the MU bound is a direct consequence of the Riesz theorem for the l p -l q norms. Please note that relatively often, integration limits in (53) were chosen as "from k∆ to (k + 1)∆" and "from lδ to (l + 1)δ"; however this choice causes a formal pathology in the limit of infinite coarse graining [166]. Thus, sticking to terminology of Equation (39), in theory it is better to avoid borderline settings for the position of the central bin, i.e., u cen = ∆/2.
To briefly report later developments, one shall mention that Partovi reconsidered the problem he had posed several years ago, pioneering applications of majorizaiton techniques [167]. Also Schürmann and Hoffmann [168] discussed the Shannon entropy UR from the perspective of the integral equation associated to it, while the first author conjectured an improvement (later mentioned in detail) which agrees with his numerical tests [169]. Finally, we mention (without details) an erroneous improvement of [31] by Wilk and Wlodarczyk [170,171], mainly devoted to the case of the Tsallis entropy.
Although originally the URs were derived for CCOs,û andv, here we show which of the URs in [24,32] can be valid also for operatorsû andv that are arbitrary linear combinations of all positions and momenta of the n−bosonic modes like the ones defined in Equation (3), viz. operators that are not necessarily CCOs. In the general case, we stress that there is always a unitary metaplectic transformation (soÛ S belongs to the metaplectic group Mp(2n, R) and it is always associated with a matrix S that belongs the symplectic group Sp(2n, R) [50]),Û S , that connectsû andv, viz.v =Û † SûÛ S . However, this metaplectic transformation is not necessarily a π/2 rotation, which would be the case ifû andv were CCOs. In order to see this, we first define two sets of operators (û,û ) T = (û = u 1 , . . . ,û n ,û 1 , . . . ,û n ) T = √ γSx and (v,v ) T = (v =v 1 , . . . ,v n ,v 1 . . . ,v n ) T = √ γ S x, whereS and S are some matrices belonging to the symplectic group Sp(2n, R), with the only restriction that the first rows ofS and S correspond to the real coefficients d and d in Equation (5), respectively, which define the operatorsû andv in Equation (3). Due to the properties of symplectic matrices, all the pairsû i andû j , and alsov i andv j , satisfy CCRs, viz. [û i ,û j ] = ihγδ ij and [v i ,v j ] = ihγδ ij with i, j = 1, . . . , n. However, it is immediate to see that (v,v ) T = S(û,û ) T where the matrix S := S S −1 is a generic symplectic matrix. Then the Stone-von-Neumann theorem guarantees that the change (û,û ) T → (v,v ) T is unitarily implementable by a metaplectic transformationÛ S [50]. In particular we haveÛ † SûÛ S = (Sx) 1 =:v.

URs Proved Only for CCOs
The key concept behind the treatment of coarse-grained URs in [24,32] is the introduction of the piece-wise continuous probability density functions: where D ∆ (u, u k ) and D δ (v, v l ) are called the histogram functions (HF) with u k (and v l in an analogous way) defined in Equation (40). Generically, these functions are defined such that they are normalized in each bin: and approach the Dirac delta distribution for infinitesimal bin size: Therefore, in the limit Z k , Z l → Z and ∆, δ → 0 we have Q ∆,u (u) → P u (u) and Q δ,v (v) → P v (v). We shall stress here that the HF can, in general, have any functional form as long as it is non-negative, normalized and fulfills Equation (56). However, the most common histogram function is the rectangular HF: with an equivalent definition for D R δ (v, v l ). In Figure 2 we show an example of coarse-grained probability distributions functions Q ∆,u (u) (the area beneath these functions are displayed in full) using rectangular histogram functions and for different size bins ∆.

z n 4 U t z Z 2 + A F e A l e g x 5 4 D w 7 A V 3 A E R g C D n + A X + A 3 + t J 6 3 P r U + t w 6 X 6 O Z G E f M M l F a r / w 9 9 Q 2 K d < / l a t e x i t > < l a t e x i t s h a _ b a s e 6 4 = " l X s d j + F v / o s A e y l r Y N 4 I T u T N o G Y = " > A A A G k n i c f Z R N b 9 M w G M e 9 w c o o b x u I E x e L a o g D m l q E B I j L Y B c O b T U m u k a S u U 4 b m r q 2 J b t b I m 8 f A e u 8 M 3 4 N j h t y p a X z l K S R 8 / / 9 7 z Z k X 3 J q D b d 7 t + N z T t 3 t r 3 t u + 3 H z x 8 9 P j J z u 7 T E y i h c k I C y b U m Y 8 0 Y Z S T k a G G k T O p C I p 8 R k 7 9 + W G u n 4 Q p a n g 3 0 0 q y T h C I a d T i p F x r h M P M T l D k 5 O d 7 + 7 W L B u 9 A q j A 4 p N N n d u v Q C g e O I c I M Z 0 v q 8 5 V m b J E y F D O S t b Y E 4 n w H I X E L p r M 4 J 5 z B X A q l H u 4 g Q t v i U O R m n k O z J C Z q a r W u 5 s L
The discrete Rényi entropy is always positive, and we have with the last line being valid because lim x→∞ (2x/π)R 2 00 (x, 1) = 1 (Equation (28) in [174] reads: This result is based on the appropriate asymptotic expansion [175] valid for z → ∞). This results show that the coarse-grained UR in Equation (63) is non-trivially satisfied for an arbitrary (even very large) values of the coarse-graining widths. However, this desired property is not enjoyed by the UR first derived in [31]. This UR corresponds to Equation (63) in the coarse-grained regime Γ/4 1.79 in which ε 1 (Γ/4) = 1/e. Obviously, this is not a mere coincidence, as Equation (63) subsumes (65). This is clearly visible inside the definition of ε which involves the minimum of two different bounds.
When Γ/4 > 1.79 the lower bound in Equation (65) is negative so this UR is trivially satisfied, since the discrete entropy is always non-negative. From the above considerations we can obtain an UR for the variances, σ 2 Q ∆,u and σ 2 , if we set α = 1 in Equation (58) and use the inequality (22): where h[·] stands for the Shannon entropy. Now, we can use the decompositions: where the variances of the discrete probability distributions were defined in Equation (42), while σ 2 D ∆ and σ 2 D δ , are the variances of the generic HFs. Therefore, applying the above splitting to Equation (66) we arrive at the lower bound [24]: When the HF are rectangular, and in the coarse-grained regime Γ/(4|γ|) 1.79 where ε 1 (Γ/4|γ|) = 1/e, we recover the UR [32]: where we have used the fact that in this case Both Equations (68) and (69) are the coarse-grained versions of the Heisenberg UR in Equation (8). It is important to emphasize that Equation (69) cannot be obtained by the simple substitution σ 2 Although both σ 2 D ∆ and σ 2 D δ are the variances of a generic HF, viz. D ∆ (u, u k ) and D δ (v, v k ) for any value of k, it is interesting to associate them to the respective central bins, namely those that contain the mean value of the probability distributions P u and P v . By doing this, together choosing the origins of the coordinates in the middle of the central bin, we can see that the variances σ 2 P (u) ∆ and σ 2 P (v) δ are free from contributions associated with the statistics relevant for the central bins. Thus, if the widths of the coarse graining increase in the measurement ofû andv, the respective central bin-widths grow, so that the variances σ 2 only involve contributions from the tails of the probability distributions Q ∆,u and Q δ,v . Therefore, for large coarse grainings, the variances σ 2 D ∆ and σ 2 D δ become more important in the inequalities Equations (68) and (69). Thus, in the regime when: both Equations (68) and (69) are satisfied trivially. Note, that in Equation (71) we have used the relation 4π 2 ≥ e 2(h[D ∆ ]+h[D δ ]) /e 2 σ 2 D ∆ σ 2 D δ > 0 which can be obtained from the inequality in Equation (22).
However, Equation (68) is only the starting point for the second construction, proposed in [24], that is free from the above limitation, and cannot be trivially satisfied. This improved UR reads: where K(t) is implicitly defined as x 0 e −y 2 dy being the error function and M −1 (t) denoting the inverse of the invertible function The idea behind derivation of the coarse-grained UR in Equation (72) is the following. Let us rewrite Equation (68) in the form:  [24] is a Gaussian with support restricted to the central bin and whose variance is an appropriate function of σ 2 D ∆ (σ 2 D δ ) (for details see [24].) Therefore, for this optimal HF its Shannon entropy , which results in the left hand side product in Equation (72).
According to the discussion above Equation (71) the coarse-grianed UR in Equation (72) has no contributions from the statistics corresponding to the central bin. In the limit when ∆, δ → 0 we recover the Heisenberg UR in Equation (8) thanks to the identities [24] lim In the opposite limit of infinite coarse graining, viz ∆, δ → ∞, we have σ 2 It is important to note that since π 2 Γ 2 ε 2 1 (Γ/4) whenever both ∆ and δ are finite, it is forbidden to set σ 2 as simultaneously equal to zero, as it would contradict the coarse-grained UR (72). This means that any quantum state (pure or mixed) cannot be localised in both observablesû andv that are CCOs. In other words, the associated probability distributions cannot simultaneously have compact support. This remarkable conclusion somehow threatens the scientific program to recover classical mechanics solely from coarse-grained averaging, physically originating from the finite-precision of the observations [19,176,177]. Indeed, quantum features can be observed in the measurement ofû andv irrespective of the precision of the detectors. However, for very large coarse-graining widths the variances σ 2 Thus, as these probabilities are likely very small, they would be particularly susceptible to statistical fluctuations and it would in general require very long acquisition times to collect the sufficient amount of data necessary to verify the UR (72) in the regime of extremely large coarse graining. (3) If we let α = 1 in Equation (58), use rectangular HFs such that Equation (60) is valid and restrict the size of the involved bins such that ε 1 (Γ/4|γ|) = 1/e-this is the regime of the coarse graining when Γ/4 1.79-we obtain the simplified coarse-grained UR of the form:

URs Valid for General Observables,û andv, Defined in Equation
Because the coarse-grained UR in Equation (58) was derived only for CCOs,û andv, a priori it is not clear why the above UR could remain valid also for generalized observables defined in Equation (3). This fact, however, can be proved with the help of the Shannon-entropy UR (20), that has properly been extended to the desired observables, and the inequalities: whose detailed derivation based on the Jensen inequality is relegated to Appendix B. Passing to the discrete entropies we find the coarse-grained UR: which looks the same as the one derived in [30] for CCOs. Here, the validity of this UR has been extended for any observablesû andv as defined in Equation (3). Also, following the same arguments that lead from Equation (66) to the UR in Equation (69) we can see that the UR for the discrete variances is also valid for generalû andv as defined in Equation (3).
To briefly summarize, entropic uncertainty relations for coarse-grained probability distributions were almost only considered for position and momentum variables. As far as we know, the only exceptions are given in References [58,65]. However, as we have shown here, the generalization of entropic URs for differential probabilities associated with general observablesû andv, which are linear combinations of position and momentum, can be done in many cases. However, in each case a careful analysis should be carried out to verify that the related coarse-grained URs are also valid for these generalised operators. Here, we have done this only in the simple cases.

Coarse-Grained URs Merged with the Majorization Approach
In [174] the coarse-grained scenario has been discussed with the help of the results obtained in [156,157,159], namely the majorization-based approach to quantification of uncertainty. To say it briefly, a majorization relation x ≺ y between two arbitrary d-dimensional probability distributions means that for every n ≤ d the inequality ∑ n k=1 x ↓ k ≤ ∑ n k=1 y ↓ k holds, with an equality (normalization) for n = d. Traditionally, by "↓" we denote the decreasing order, so that x ↓ k ≥ x ↓ l , for all k ≤ l. The Rényi entropy (and also others, such as the Tsallis entropy) is Schur-concave, which implies In the context of coarse-grained probability distributions it was conceptually simpler to consider the so-called direct-sum majorization introduced in [159]. An advantage of the majorization approach is that it covers a regime of (α, β) parameters, β = α to be precise, which in some way is perpendicular to the conjugate choice 1/α + 1/β = 2. In [174] an infinite hierarchy of majorization vectors, depending on a single parameter Γ = ∆δ/h, has been derived. The discussion is conducted for CCOs, thus one can easily recognize the dimensionless Γ parameter as those which appears in all previous URs with γ = 1.
The main result, namely a family of lower bounds denoted as R (n) α (∆δ/h) for n = 2, . . . , ∞, has been presented in Equation (27) from [174], however, we refrain from providing its detailed construction here. It seems enough to say that the bound in question is a function of R 2 00 (j 0 Γ/4, 1) with j 0 being certain positive integers. In other words, in spirit, the majorization bound is close to that derived in [24] and extensively discussed above. A comparison of the new bound and (63) for α = 1 = β-the only value of both parameters for which the involved bounds describe the same situation--showed that R 1 outperforms (63) in the regime when the R 00 -term does contribute to ε 1 . Asymptotic behavior of the new and previous coarse-grained bounds shows that for α = 1 = β and large Γ, all R (n) 1 bounds improve (63) by a divergent factor Γ/4. Moreover, the typical behavior of discrete majorization bounds has been confirmed in the coarse-grained setting. In the discrete case, the majorization relations almost surely dominate the MU bound, with an exception being a small neighborhood of the point for which the unitary matrix U is the Fourier matrix. The analog of the Fourier matrix in the coarse-grained scenario is the continuous limit Γ → 0. This probably intuitive fact has been rigorously shown by means of the asymptotics of R (∞) 1 for small Γ, which is equal to − 1 2 ln Γ.

Other Coarse-Grained URs
At the very end of this long section we would like to touch on a few coarse-grained URs which go beyond the standard position-momentum conjugate pair. First of all, Bialynicki-Birula also provided his major Shannon entropy UR in the case of angle and angular momentum [30], as well as (together with Madajczyk) to the variables on the sphere [178]. Coarse graining in these physical settings is only relevant for the periodic CVs (angle on a circle and two angles on a sphere), as the conjugate variables are discrete (though infinite dimensional).
Also, the coarse-grained scenario has been developed [179] in relation to the memory-assisted UR [180] relevant for quantum key distribution. The result, even though non-trivial, differs from Equation (63) in a similar fashion as the MU bound differs from the UR in the presence of quantum memory by Berta et al [180].
Going in a completely different direction, Rastegin [181] in his recent contribution proposed an extension of (65) to the case of a modified CCR, which assumes the form [x,p] = ih(1 + βp 2 ). The parameter β is related to the so-called minimal length predicted by certain variants of string theory and similar approaches (not to be confused with β playing the role of a conjugate parameter in the MU bound and similar URs for the Rényi entropies).
Last but not least, some of us have very recently derived an inequality (see Equations 9-12 from [182]), which could be understood as an UR (valid for CCOs) in the setting relevant for periodic coarse graining discussed in Section 4.1.2. As this UR involves additional averaging of p (x) k (ρ) and p (p) l (ρ) defined below Equation (47) with respect to the positioning degrees of freedom, we do not provide further details of this construction encouraging the interested reader to consult [182].

Applications of Coarse-Grained Measurements and Coarse-Grained Uncertainty Relations
As discussed above, when detecting the position and momentum of particles such as photons or individual atoms, coarse-grained measurements are not just necessary but can be much more practical. In this regard, URs that deal with coarse-grained measurements can be useful for several applications, such as those discussed in Section 3.
Section 3 discussed the use of URs along with the PPT arguement for the convenient detection of quantum entanglement in continous variable quantum systems. However, sufficient care must be taken in regards to coarse-grained measurements. The pitfalls of applying the usual entanglement criteria for continuous variables to coarse-grained measurements was discussed in Reference [26], where it was argued that this can lead to false-positive identifications of entanglement, such that the entanglement criteria based on uncertainty relations discussed in Section 3 can be (falsely) violated even for separable states. For a simple illustration of this, consider the most trivial example of a separable continuous variable state, the two-mode vacuum state [10]. Even though the state is separable, improper application of entanglement criteria without correctly taking account of coarse graining can lead to erroneous results. A demonstration of this is shown in Figure 5. We consider the results from coarse grained measurements, and apply entanglement criteria based on an ideal continuous variable UR and its coarse-grained version. The red circles show a variance-based entanglement criteria based on the variance product UR Equation (8), using the global operators defined in Equations (29) and (30), as developed in Reference [99]. Here we have subtracted the lower bound from the product of variances, so that a negative value indicates entanglement. As can be seen, when the coarse graining is large, one would erroneously conclude that the quantum state is entangled. On the other hand, the coarse-grained variance product UR (69) applied to the global operators (29) and (30) never indicates that the state is entangled, as indicated by the blue squares in Figure 5. Similar results hold for other UR-based entanglement criteria.
To show how coarse-grained data should be properly handled to identify entanglement, an experimental study was performed in a system of spatially-entangled photons [26]. In particular, the same variance criteria based on (69) was tested for the global operators defined in Equations (29) and (30), in which case entanglement was identified for a wide range of coarse-graining widths. It was also shown that coarse-grained entropic entanglement criteria, for example based on inequality Equation (58) (α = β = 1) applied to operators (29) and (30), can be superior to coarse-grained variance-based criteria, identifying entanglement when variance criteria do not, even for the case of Gaussian states. This is due to the fact that the coarse-grained probability distributions functions such as those shown in Figure 2 are not Gaussian functions, even when the quantum state under investigation is Gaussian.
An advantage of coarse-graining is that the measurement time can be drastically reduced. In Reference [86] EPR-steering was tested for discrete distributions of measurements made from standardized binning on the two-photon state produced from spontaneous parametric down-conversion, using a coarse-grained version of the EPR-steering criteria of Reference [85]. Bi-dimensional steering was observed for sample sizes ranging from 8 × 8 to 24 × 24, representing a considerable reduction in measurement overhead when compared with the quasi-continuous measurements reported in Reference [85], which sampled about 100 data points per cartesian direction (about 10 4 total measurements) to evaluate entropic EPR-steering criteria of continuous variables.
Standard coarse graining has been studied in the context of quantum state reconstruction of single and two-mode Gaussian states, and the quantum to classical transition [183]. Two scenarios were considered: direct reconstruction of the covariance matrix alone, and full reconstruction of the state using maximum likelihood estimation. The reconstructed coarse-grained functions were compared to those of Gaussian states subject to thermal squeezed reservoirs, indicating that in this context coarse graining does not produce a thermalized (decohered) Gaussian state. Numerical results testing entanglement criteria for the two-mode vacuum state, a separable pure state. The entanglement criteria are based on URs following the PPT argument outlined in Section 3. The criteria are evaluated as a function of the bin widths ∆ = δ, which are given in units of the standard deviations σ P u and σ P v . We note that σ P u = σ P v for the two-mode vacuum state. The red circles show the variance product UR Equation (8), where we apply the naive approach in which the variances of the continous variables are calculated from the discretized data using Equation (42). One can see that in this case we obtain a false-positive for entanglement when the coarse-graining widths are large. The blue squares show the coarse-grained variance product UR Equation (69), both applied to the global operators Equations (29) and (30). Here the lower bounds for both inequalities have been subtracted, so that a negative value indicates entanglement. The lines are merely guides for the eye.
The work mentioned above considered standard coarse graining, as described in Section 4.1. In some cases it is interesting to consider different models, such as that of periodic coarse graining described in Section 4.1.2. The mutual unbiasedness of periodic coarse graining described in Section 4.2 has been tested experimentally for two [142] and even three [143] phase-space directions. It was shown that mutual unbiasedness appears when the appropriate bin widths of the two or three conjugate variables are chosen. Periodic coarse graining has also been used in the detection of spatial correlations of photon pairs from SPDC [182]. Using a novel entanglement criteria based on the UR for characteristic functions [153], it was possible to identify entanglement with as few as 2 × 2 measurements in position and momentum (8 total), representing a considerable reduction in measurement overhead.
Simple binary binning of homodyne measurements has been proposed as a means to test dichotomic Bell's inequalities in CV systems, while allowing for high detection efficency [184][185][186][187]. Other types of non-standard coarse graining have been proposed as a means to violate Bell's inequality using homodyne measurements on non-Gaussian states [188]. Though it was shown that one could achieve maximal violation in principle, exotic non-Gaussian states are required. In Reference [27] it was shown that imperfect binning could result in false violations of Bell's inequalities, and even in violations of Cirelson's bound for quantum Bell correlations.
A closely related subject to periodic coarse graining of CVs is that of the so called modular variables [189][190][191], for which phase-space variables u are rewritten as u = n u +ū, where n u is the integer component andū the modular component, such that 0 ≤ū < . Here is a scaling parameter of appropriate dimension. For two CCOs, such asx andp for example, the integer operator of one observable-say-n x and the modular operator of the other observable-p satisfy URs that closely resemble those of the angle and angular momentum variables [30]. The modular variable construction was first introduced by Aharanov et al. [138,189] as a method to identify non-locality in quantum mechanics. Since then, several interesting applications have been developed. Variance-based URs for the modular variable construction were proposed as a method to identify a novel type of squeezing, as well as entanglement in pairs of atoms [192]. This entanglement criteria was used in References [193], along with one based on entropic uncertainty relations, to identify spatial entanglement of photon pairs that have passed through multiple slit apertures. Application to multiple-photon states was studied in Reference [194]. It is worth noting that in this case the usual CV entanglement criteria as discussed in Section 3 are incapable of detecting entanglement. Modular variables have been proposed as a way to test for the Greenberger-Horne-Zeilinger paradox in CV systems [195], as well as quantum contextuality [196][197][198] and as a method to construct algebras resembling that of discrete systems [190,191,199].
Finally, we briefly mention that URs play an important role in the attempt to unify quantum theory with general relativity. In this case, the Heisenberg uncertainty principle is modified to become a generalized uncertainty principle, taking into account Planck scale effects, which impose coarse-graining that is a fundamental part of nature, leading to minimum and maximum length quantum mechanics. An extensive amount of literature exists on the subject, for two recent reviews, see References [200,201].

Conclusions
Uncertainty relations play an important role in quantum physics, which is two-fold: on the one hand they have historically represented the difference between classical and quantum physics, while on the other hand they are a tool that can be used to identify and even quantify interesting quantum properties. Beginning with the seminal work of Heisenberg in 1927, several uncertainty relations have been developed for continuous variable quantum systems. However, in a realistic experimental setting, one never has access to the infinite dimensional spectrum associated to these observables. Thus, coarse graining is imposed by the detection apparatus to account for the measurement precision and range.
Here we have provided a review of several quantum mechanical uncertainty relations tailored specifically to coarse-grained measurement of continuous quantum observables. Our aim was to survey the state-of-the-art of the subject, from both the theoretical advances to experimental application of coarse-grained uncertainty relations. We also extend the validity of some of the coarse-grained URs, already in the literature, to general linear combinations of canonical observables in n-mode bosonic systems.
Several interesting open questions remain. First, it would be interesting to see the generalization of all the coarse-grained URs presented here for pairs of observables that are connected by general unitary metaplectic transformations. Second, one can consider applying coarse graining to URs not mentioned explicitely here, such as the triple variance product criteria [120,150], UR for characteristic functions [153], among others, as well the plethora of moment inequalities arising from tests for non-classicality [68,72] and entanglement [117,118,121]. Third, and more important, a deep discussion of the role of coarse-grained URs within the scientific program to recover classical mechanics solely from coarse-grained averaging should be developed. We hope that this review encourage this discussion.
Author Contributions: All authors contributed equally to this work.

Acknowledgments:
The authors acknowledge financial support from the Brazilian Funding Agencies CNPq, CAPES (PROCAD2013 project) and FAPERJ, and the National Institute of Science and Technology-Quantum Information. Ł.R. acknowledges financial support by Grant number 2015/18/A/ST2/00274 of the National Science Center, Poland.

Conflicts of Interest:
The authors declare no conflict of interest. The founding sponsors had no role in the design of the study, in the collection, analyses, or interpretation of data, in the writing of the manuscript, and in the decision to publish the results.
After using the condition in Equation (36) which is the desired result.