Hierarchy of continuous-variable quantum resource theories

Connections between the resource theories of coherence and purity (or non-uniformity) are well known for discrete-variable, finite-dimensional, quantum systems. We establish analogous results for continuous-variable (CV) systems, in particular Gaussian systems. To this end, we define the concept of maximal coherence at fixed energy, which is achievable with energy-preserving unitaries. We show that the maximal Gaussian coherence (where states and operations are required to be Gaussian) can be quantified analytically by the relative entropy. We then propose a resource theory of non-uniformity, by considering the purity of a quantum state at fixed energy as resource, and by defining non-uniformity monotones. In the Gaussian case, we prove the equality of Gaussian non-uniformity and maximal Gaussian coherence. Finally, we show a hierarchy for non-uniformity, coherence, discord and entanglement in CV systems.


Introduction
Quantum resource theories [1] describe the resources of quantum states in a quantitative way.The set of states is divided into free states (having no resource) and resource states.Quantum operations are called free when they transform any free state into a free state, i.e. free operations cannot increase the resources.
Different resource theories have different sets of free states and operations.For instance, the resource theory of coherence [2][3][4][5][6][7] identifies states that are diagonal in a certain basis as free, while the resource theory of purity [8][9][10] considers the maximally mixed state as the only free state.The resourcefulness of a quantum state can be quantified by a resource monotone, which is a function that is non-increasing under free operations.In particular, relevant monotones for several resource theories are based on the relative entropy [2,9,11,12].
For discrete-variable (DV) systems, an important connection between coherence, purity, entanglement and discord was found in [10].Here, it was shown that the purity of a quantum state is -for an appropriate resource monotone -equal to the maximal coherence that can be obtained by applying unitary operations to the state.This quantity then upper-bounds the maximal quantum discord and entanglement of the state.This result cannot be straightforwardly extended to continuous-variable (CV) systems.Infinite-dimensional Hilbert spaces are structurally different from their finite-dimensional counterpart [13][14][15][16], and this difference influences the mathematical definition of the physical quantities themselves and their resource theories [17].In particular, infinite-dimensional Hilbert spaces allow the generation of infinite resources via unitary operations [18,19].
A common strategy to replicate DV results is to introduce valid physically and experimentally motivated constraints.The first restriction is to consider finite energy [18,20,21] and to explore the use of energy-preserving unitaries, as those operations are easily available in laboratories.The second restriction is to focus on Gaussian states and operations, as most of the relevant phenomena in quantum information and quantum optics can be described by at most quadratic Hamiltonians.
Equipped with these assumptions, we investigate the resource theories of CV coherence [6,7].We discuss both general and Gaussian coherence, the latter by restricting the set of quantum states and operations to be Gaussian.We define the concept of maximal coherence of a quantum state at fixed energy, as the coherence that can be obtained via applying energypreserving unitaries.In the case of Gaussian coherence, we find the structure of the states with maximal coherence and discuss the form of the maximizing unitary in some specific cases.We derive an analytical expression for the relative entropy of the maximal Gaussian coherence.We then propose a resource theory of non-uniformity, to describe purity at fixed energy as resource, considering states that maximise the entropy at fixed energy as free.However, our resource theory studies the interactions of a quantum system with a noisy thermal environment, therefore it is connected to the resource theories of a-thermality [22,23], where states out of thermal equilibrium are identified as resources.Our theory emphasises the entropic exchanges between the system and the environment, rather than the energetic ones.
Finally, we establish a connection between non-uniformity, coherence, quantum discord and entanglement, by identifying a hierarchy between them.In particular, the maximal Gaussian coherence is bounded by the Gaussian non-uniformity and both upper-bound the maximal discord and the entanglement.Our results represent an extension to infinite dimensions of the hierarchy found in [10].
We begin our work by introducing our setting together with the basic notions of CV quantum information in Sec. 2. We review coherence and Gaussian coherence in Sec. 3. In Sec. 4, we define the maximally coherent mixed state at fixed energy and derive its properties for the Gaussian case.We assemble the resource theory of non-uniformity in Sec. 5 and illustrate the connections between non-uniformity, coherence, discord and entanglement in Sec. 6.Finally, we summarise our results in Sec. 7.

Notation and preliminaries
In the following, we will indicate vectors and matrices as bold lowercase and uppercase letters, respectively.We shall consider systems with a finite number of discrete spectral and spatial modes, which refer to the frequency and location of the mode, respectively.We label the mode operators of a mode with two indices: an index ω for the spectral degrees of freedom and an index j for the spatial degrees of freedom.
We sort our mode operators by gathering all mode operators with the same frequency, i.e.
where M f is the number of different frequencies, and for each frequency ω âω = (â ω;1 , âω;2 , . . ., âω;Ms ) T , with âω;j being the annihilation operator for a mode with spectral label ω and spatial label j, and M s being the total number of spatial labels.Without loss of generality, we assume that M s is the same for all frequencies.A graphical depiction of our mode labeling is drawn in Fig. 1.The operators satify the usual bosonic commutation relations âω;j , âω ;j = â † ω;j , â † Using the notation of Eqs. ( 1) and ( 2), the free Hamiltonian of the system reads In the following, we will always assume that the system has finite mean energy, i.e. that the Hamiltonian satisfies Ĥ < ∞.This is a natural and physically reasonable assumption, and a mathematically necessary precondition for the trace-norm continuity of many functionals [18,20,21].The total number of modes is M = M f M s .We shall label the modes with the index m ≡ ω; j whenever there is no need to distinguish between spectral and spatial modes.We use the order m = (ω 1 ; 1), (ω 1 ; 2), . . ., (ω 1 ; M s ), (ω 2 ; 1), . . ., (ω M f ; M s ).
We now recall basic concepts of Gaussian quantum information, inspired from [13][14][15][16] and written according to our notation.
For each mode m, we define the canonical conjugate operators qm : . We group them in a vector r := (q 1 , p1 , q2 , p2 , . . ., qM , pM ) T .The first and second moment of a state ρ is the displacement vector d and the covariance matrix V , respectively.Their components d m and V mm read in terms of components rm of r as [15]: Both moments are real, and V is symmetric and positive definite.Gaussian states are represented by a Gaussian quasi-probability distribution in the phase space and are fully characterised by the first and second moments.The 2M -dimensional displacement vector and the 2M × 2M covariance matrix of any Gaussian state can be written in the following block form: where d m are 2-dimensional vectors, and V m and ∆ mm are 2 × 2 real matrices.In particular, d m (V m ) corresponds to the displacement vector (the covariance matrix) of the reduced state ρ m = Tr m\m [ ρ ] after partial trace of all modes but the m-th, while ∆ mm is related to the correlations between the modes m and m [13].
The total average occupation number can be derived as where nm is the average occupation number of the m-th mode.This expression can be obtained by writing Tr[V m ] in terms of Eq. ( 7) and the mode operators âm and â † m Eq. ( 9) follows by setting N =

M m=1
â † m âm .Notice that this formula holds regardless of the presence of correlations between different modes, for both Gaussian and non-Gaussian states.By Williamson's theorem [24], any covariance matrix V can be brought into a diagonal form: where D is a diagonal matrix and S is a symplectic matrix, i.e. a real matrix that satisfies The variables ν m ≥ 1 are called symplectic eigenvalues and obey the Bose-Einstein statistics: where κ is the Boltzmann constant and T ω;j is the temperature of the m-th mode, with m = (ω; j).The symplectic eigenvalues are used to express several properties of Gaussian states.
For instance, the von-Neumann entropy of a Gaussian state reads [25]: Any unitary that preserves the Gaussianity of quantum states is called Gaussian unitary.In terms of the moments, a Gaussian unitary acts as [16] where v is a 2M -dimensional vector and S is a 2M × 2M symplectic matrix.
Theorem 2.1 (Bloch-Messiah decomposition [26]).Any 2M × 2M symplectic matrix can be decomposed as where the 2M × 2M matrices O 1 and O 2 are symplectic and orthogonal (we generally denote symplectic orthogonal matrices as O), and the 2 × 2 matrix Z(r m ) is a single-mode squeezer with squeezing parameter r m , that is For r m = 0 (absence of squeezing), Z(r m ) becomes the identity I. Therefore, in Eq. ( 17), S = O 1 O 2 is orthogonal, being the product of two orthogonal matrices.Definition 2.2.A passive unitary [13,14] is a Gaussian unitary ÛO that is represented in the phase space (Eq.( 16)) by v = 0 and a symplectic orthogonal matrix O in the form of where O ω are 2M s × 2M s symplectic orthogonal matrices acting on the subset of modes with frequency ω.Conversely, any Gaussian unitary that is not passive, is called active.
Active unitaries are associated with linear displacements and squeezing.Passive unitaries are realised by linear-optics circuits, that is any multiport interferometer made of beam splitters and phase shifters.They preserve the total average occupation number (see Eq. ( 9)): Since N commutes with the Hamiltonian Ĥ of the system, passive Gaussian unitaries preserve the average energy, too.Passive Gaussian unitaries are the only energy-preserving Gaussian unitaries [15,16].They create correlations between spatial modes with the same frequency ω.
They do not allow interactions between modes with different frequencies (see Eq. ( 19), and [27]).Gaussian unitaries are not the only operations that preserve the Gaussianity of a quantum state.A Gaussian channel is a completely positive trace-preserving (CPTP) operation that maps Gaussian states into Gaussian states.In terms of the moments, a Gaussian channel acts as [28] where T , N are 2M × 2M real matrices, N ≥ 0, and v is a 2M -dimensional vector.
3 Resource theory of (Gaussian) coherence for an arbitrary set of non-negative probabilities { p n 1 ...n M }.
We denote the set of all incoherent states by I.The resource theory of coherence admits different sets of free operations.The maximal set of free operations for I are called maximally incoherent operations (MIO) [2].These are all maps that cannot create coherence, i.e., An extensive study of MIO in infinite-dimensional Hilbert spaces has not been carried out so far.We shall not investigate this here, referring to [5] for a general review of coherence and to [6] for the CV case.Definition 3.1.For continuous-variable systems, a function C(ρ) is a suitable measure of coherence with respect to a chosen set of free operations, e.g.MIO (see Eq.( 23)), if it satisfies the following properties [3,6,7]: (C1) Positivity: C(ρ) ≥ 0 for any density operator ρ and C(ρ) = 0 iff ρ ∈ I; (C2) Monotonicity under the chosen set of free operations, e.g. for MIO: (C3) Convexity: (C4) Finite coherence for systems with finite energy: The condition (C4) is specific for CV systems.Denoting with S(ρ τ ) the quantum relative entropy between ρ and τ , S(ρ τ the relative entropy of coherence is defined as: This measure satisfies all conditions of Def.3.1 [6].
The generic CV coherence has not been deeply studied, mainly because of theoretical and experimental difficulties associated with general bosonic Hilbert spaces.We shall therefore address now the relevant Gaussian subclass, in which all states and operations are Gaussian1 .
In the realm of Gaussian states, a state is diagonal in the Fock basis if and only if it is a thermal state [7], defined as where n = (n 1 , . . ., nM ) and nm is the average occupation number of the m-th mode (see Eq. ( 9)).The subscript M in τ M denotes the number of modes of τ M , and is omitted for singlemode thermal states.Thermal states have a zero displacement vector and a diagonal covariance matrix We denote the subset of all incoherent Gaussian states by I G .Xu [7] introduced incoherent Gaussian operations (IG).Definition 3.2.Incoherent Gaussian operations (IG) are defined as all Gaussian channels Λ IG that map thermal states (Eq.( 29)) into thermal states, i.e. the maximal set of free operations in this scenario.A generic IG can be written in the form of Eq. ( 21), where with ω m ≥ 0; • T IG is composed of M f submatrices T ω that act on single frequency sectors, namely: and each T ω can be generated as follows: 1. Take M s real coefficients t ω;j ∈ R; 2. Take M s 2 × 2 orthogonal matrices O ω;j , which do not need to be symplectic; 3. T ω is given by a permutation of the columns of  29).This function quantifies the Gaussian coherence, which is, in general, an upper bound for the general coherence (since In particular, the relative entropy of Gaussian coherence C G rel (ρ) can be defined as the relative entropy of coherence (see Eq. ( 28)) by performing the minimization over I G .For M -mode Gaussian systems, it reads [7]: where nm is the reduced average occupation number of ρ (see Eq. ( 9)), τ M ( nρ ) represents the thermal state with the same nm of ρ and the von-Neumann entropy S(ρ) is given by Eq. ( 15).

Maximally coherent mixed states at fixed energy
Coherence is a basis-dependent quantity, therefore it is affected by unitary operations.For a given DV state ρ, the maximally coherent mixed state (MCMS) [29,30] is defined as ρ max = V ρ V † , where V is the unitary that maximises the coherence of ρ.
For CV systems, this definition is not applicable, since the coherence depends on the energy of the system.This can be seen, for instance, in Eq. (33) for the relative entropy of Gaussian coherence.Therefore, energy non-preserving unitaries can in principle increase the coherence indefinitely.
From an experimental point of view, energy-preserving unitaries are easier to realise and do not require interaction with an external source of energy.This is the motivation to define a family of maximally coherent mixed states at fixed energy: where ÛEP are energy-preserving unitaries.If we consider Gaussian states and operations in Eq. ( 34), then ρ max is the maximally coherent mixed Gaussian state (MCMGS) at fixed energy with respect to C.
Let us now focus on the Gaussian case.Passive Gaussian unitaries ÛO (see Def. 2.2) are the only energy-preserving Gaussian unitaries, i.e. they preserve N ω = j nω;j , and thus N = ω N ω .However, the interaction between modes with the same frequency ω results in a redistribution of nω;j .
Let us consider a generic Gaussian state ρ and call ρ ω = Tr ω\ω (ρ) the state obtained by tracing out all modes with any frequency but ω.From its definition in Eq. ( 28), the relative entropy of Gaussian coherence of ρ can be written as the sum of the relative entropies for all ρ ω : Therefore, the maximal Gaussian coherence can be obtained by maximizing the Gaussian coherence for each ρ ω .
Theorem 4.2.Among all Gaussian states ρ ω with a given symplectic spectrum { ν 1 , . . ., ν Ms } and a given average occupation number N ω = j nω;j , the states with equidistributed reduced average occupation numbers, i.e. nω;j = N ω /M s , ∀ j, are the maximally coherent mixed Gaussian states with respect to the relative entropy of Gaussian coherence C G rel (ρ ω ) (see Eq. ( 33)).The proof is given in Appendix A, using Lagrange multipliers.Combining this result with Eq. ( 33), it follows that the maximal relative entropy of Gaussian coherence of a Gaussian state ρ reads where τ M ( nρmax ) is the M -mode thermal state with occupation numbers nρmax (see Eq. ( 29)).
In Appendix B, we provide an alternative analytical expression for C G rel; max (see Theorem 4.3), that will play a pivotal role in the next sections.
Finally, in Appendix C, we provide examples of passive unitaries that maximise the Gaussian coherence for two generic classes of Gaussian states.

Resource theory of (Gaussian) non-uniformity
The resource theory of purity (or non-uniformity) [9] belongs to a family of resource theories of quantum thermodynamics in which states out of some form of equilibrium are considered as resources [31][32][33][34].Usually, this equilibrium is given by assuming the environment at a certain background temperature T : the free states are thermal states at the same temperature and the free operations are those which are generated by an energy-preserving unitary acting on the system and the environment.
For DV systems, the resource theory of purity arises when the Hamiltonian is fully degenerate at any temperature.Then all unitaries become energy preserving (hence free operations) and the exchanges between the system and the environment are purely entropic [9].The state representing informational equilibrium becomes the maximally mixed state (MMS) I/d, where d is the dimension of the system.The MMS is the only free state, as every other state possesses some non-uniformity.
The DV theory cannot be straightforwardly extended to Gaussian systems, because a proper MMS is nonphysical, as it is associated with infinite energy in infinite-dimensional Hilbert spaces [16].
While several resource theories of Gaussian states out of thermal equilibrium exist [22,23,35], we choose a different approach, that emphasises the informational aspects of the interactions between the system and the environment over the energetic ones.We consider purity at given energy as a resource and we refer to this resource theory as non-uniformity: with a similar argument as given by Gour et al [9], we use this term in place of "purity" because we shall consider pure states at different energy as states with different resource content.
Consider an M -mode state ρ, with M f frequencies and M s spatial labels (see Sec. 2).Consider also a Hamiltonian Ĥ in the form of Eq. ( 5).In terms of the frequencies its mean energy can be written as where E ω = ωN ω (see Eq.( 5)) and any other contribution is set to zero by a suitable choice for the zero-point energy.In the DV resource theory of purity, a dimension d for the set of all states was fixed.In our resource theory, we fix a set of frequencies ω = ω 1 , ω 2 , . . ., ω M f and the energy in each frequency mode E ω , thus fixing Ĥ as in Eq. (38).The states and modes can have any temperature that is compatible with E ω , i.e. the thermal component of the energy cannot be higher than E ω for any frequency sector.For a single frequency ω, the state that maximises the entropy is the natural CV counterpart of the maximally mixed state in the DV case.We call it the uniform state at frequency ω.In our setting, with M f different frequencies, we consider as free state the tensor product of all uniform states at frequency ω over all frequencies.We call it the uniform state.
From Eq. ( 9), we see that the average energy E ω = ω N ω is a function of the first and second moments of ρ, for both Gaussian and non-Gaussian states.It is well known that Gaussian states attain the maximum von-Neumann entropy among all states having the same displacement vector and covariance matrix [25].Therefore, even if we consider the set of all CV states, we can search for the uniform state in the subset of Gaussian states.
For single-mode Gaussian systems at a fixed energy E ω = ωN ω , the von-Neumann entropy is maximised by Gaussian thermal states with average occupation number N ω [36].For M s spatial modes, we prove in Appendix D the following result: Theorem 5.1.For a quantum system of M s spatial modes, the uniform state at frequency ω, i.e. the state that maximises the entropy, is the Gaussian thermal state (see Eq. ( 29)) with equal single-mode occupation numbers, i.e.
where δ ω = (δ ω , δ ω , . . ., δ ω ), δ ω := N ω /M s , and N ω is the total occupation number for all spatial modes with frequency ω.Considering both spatial and frequency modes, the uniform state becomes where M f is the total number of frequencies.The covariance matrix of τ M (δ) reads This result has an intuitive explanation.The von-Neumann entropy of a Gaussian state (see Eq. ( 15)) depends solely on the symplectic eigenvalues, and the m-th eigenvalue is a function of the m-th mode's temperature (see Eq. ( 14)).To maximise the entropy, we need to consider thermal states.Among thermal states, the uniform state is defined as the state with the most homogeneous distribution of single-mode energies.
Both the general and the Gaussian version of this resource theory have the same set of free states.We could, in principle, distinguish them via the set of free operations.
We call the Gaussian channels in UP the uniformity-preserving Gaussian operations (UPG).
Clearly U P G ⊆ U P , but we do not know whether this inclusion is strict.Concerning Gaussian operations, a more practical set of free operations is that of Gaussian noisy operations (GN), i.e.Gaussian channels Λ GN that admit the following decomposition: where τ M E (δ) is the uniform state (see Theorem 5.1) for M E environmental modes and with the same δ of the system (see Eq.( 40)), and Gaussian noisy operations preserve the equilibrium state.This can be seen in phase space representation, since, for every frequency sector, the covariance matrix of τ M (δ) ⊗ τ M E (δ) is proportional to the identity (see Eq. ( 41)), and the symplectic matrix of O is orthogonal.Clearly, GN ⊆ U P G, but also here we do not know whether this inclusion is strict.
We introduce a quantifier for the resource of non-uniformity as follows: Definition 5.3.A function P, mapping density operators to real numbers, is a non-uniformity monotone if (P1) P is non-negative and vanishes for the uniform state (see Theorem 5.1).
(P2) P does not increase under the chosen set of free operations, for instance UP (see Def. 5.2), i.e.
In analogy with coherence, we introduce the relative entropy of non-uniformity: This function clearly satisfies (P1).The property (P2) follows from the contractivity of the relative entropy, Restricting ourselves to Gaussian states and operations, we find results for the relative entropy of Gaussian non-uniformity P G rel (ρ) (that is the relative entropy of non-uniformity for Gaussian states).
Theorem 5.4.The relative entropy of Gaussian non-uniformity (see Eq. ( 45)) of a Gaussian state ρ is equal to its maximal coherence (see Eq. ( 36)): This result follows from Theorem 4.3: and establishes a strong connection between coherence and non-uniformity for Gaussian systems, in analogy to DV systems [10].
We conclude this section by noticing two additional properties of the relative entropy of Gaussian non-uniformity, which can be found by employing Theorem 5.4: • Pure Gaussian states |ψ G ψ G | maximise the relative entropy of Gaussian non-uniformity among the states with given average occupation number N ω and total number of spatial modes M s : which follows from Eq. ( 33).
• The relative entropy of Gaussian non-uniformity is invariant under passive Gaussian unitaries, i.e.
This property follows by noticing that the maximal Gaussian coherence cannot be increased via passive Gaussian unitaries (see Eq. ( 34)).

Hierarchy of quantum resources in CV systems
The relative entropy also quantifies, for CV systems, multipartite entanglement [18] and symmetric quantum discord [37,38]: where Z and S denote the sets of zero-discord and separable states, respectively.The former contains all mixtures of pure, locally orthonormal projectors [37,39], while the latter contains all convex combinations of arbitrary product states [36], i.e. where Note that Z is non-convex, since the convex combination of two sets of orthonormal projectors is not, in general, orthonormal.
Using the relative entropy, we can extend the ordering of resources for discrete-variable to continuous-variable systems: This relation directly follows by noting that τ M (δ) ∈ I ⊂ Z ⊂ S, and holds for all quantum states ρ (see Fig. 2).Let us now consider the Gaussian case.Let D G rel and E G rel be the relative entropies of Gaussian discord and entanglement, respectively.They are obtained with Eq. ( 51) and (52) by performing the minimization over the Gaussian subsets Z G and S G of Z and S, respectively.While S G is defined analogously to Eq. (54), by taking Gaussian states, Z G is formed by product Gaussian states [38,40], i.e.
We have discussed in Sec. 4 how passive unitaries can generate coherence.It is well established that they can also generate entanglement [41] and discord [42].
Let us introduce We prove in Appendix E the following hierarchy between the mentioned CV resources (see Fig. 3): Theorem 6.1.The relative entropy of Gaussian non-uniformity P G rel (Eq.(45)) of any Gaussian state ρ is equal to the maximal relative entropy of Gaussian coherence C G rel; max (Eq.( 36)), and this quantity upperbounds the maximal relative entropies of Gaussian symmetric discord D G rel; max (Eq.(58)) and Gaussian entanglement E G rel; max (Eq.(59)): Here we discussed the action of energy-preserving unitaries, in particular passive Gaussian unitaries.An energy non-preserving unitary can, in principle, increase the energy indefinetely and create infinite resources.However, for a fixed finite energy, the ordering of Eq. ( 55) is preserved, because the ordering of the sets remains.In the Gaussian scenario, we conjecture that active unitaries exist that keep the hierarchic ordering in Eq. ( 60).The verification of this claim is an interesting open question.

Conclusions
In this manuscript, we extended a hierarchy of dicrete-variable quantum resources to continuousvariable systems, under the condition of fixed energy.Considering Gaussian states and operations and using quantifiers based on the relative entropy, we found that the Gaussian nonuniformity is equal to the maximal Gaussian coherence, and we provided an analytical expression for this quantity.This means that, if we quantify the resources with the relative entropy, any amount of Gaussian non-uniformity can be converted into Gaussian coherence by means of a suitable energy-preserving Gaussian unitary.To quantify the non-uniformity, we designed a resource theory by identifying purity at fixed energy as resource.We also considered generic (non-Gaussian) states and found that the non-uniformity always upper-bounds the coherence.Finally, we showed that, for Gaussian states the non-uniformity and the maximal coherence provide upper bounds on the maximal symmetric quantum discord and the maximal entanglement.Our results advance the field of continuous-variable resource theories and establish a further connection between quantum thermodynamics and quantum information theory.
Our work leaves some interesting questions open.A possible next step could be to study the hierarchy of resources in the presence of energy-nonpreserving Gaussian unitaries, up to a finite maximum energy.In addition, one should investigate whether the equality of maximal coherence and non-uniformity also holds in the general non-Gaussian case.In order to achieve this a deeper understanding of non-Gaussian resource theories is required.
The theorem can be proven with a constrained optimization.Let n = { n1 , . . ., nMs } and L(n, λ) be the Lagrangian function where Ms is the number of modes for which nj = 0. Since S(ρ) depends only on the symplectic spectrum (see Eq. 15), it holds The condition ∂L/∂ nj = 0 for nj = 0 is equivalent to The above relations are satisfied by any state ρ * with Ms reduced occupation numbers nj = N/ Ms ∀ j = 1, . . .Ms and the others equal to zero.Clearly 1 ≤ Ms ≤ M s .
The Gaussian coherence of any ρ * reads: Taking the derivative of this expression with respect to Ms , i.e.
one can verify that this function is monotonically increasing with Ms .The minimum is therefore obtained for Ms = 1, i.e. when N is contained in a single reduced spatial mode, and the maximum is obtained for Ms = M s , i.e. when N is equally distributed.

B Proof of Theorem 4.3
Let ÛC = ω M f ω=ω 1 ÛCω be the passive unitary that maximises the relative entropy of Gaussian coherence for ρ, i.e.
Let τ M ( nρmax ) be the thermal state with the same nω;j as ρ max .Using Theorem 4.2, τ M ( nρmax ) reads Using Eq. ( 33) we get Here, we have used the invariance of entropy under unitary operation and the diagonality of τ ( nω ρmax ) in the Fock basis.We then notice that This can be proven by using the phase space representation, since the covariance matrix of τ Ms ( nω ρmax ) is proportional to the identity (consider equal nm in Eq. ( 30)), and passive Gaussian unitaries are associated to symplectic orthogonal matrices, by Def.2.2.It follows that

C Maximal Gaussian coherence for specific states
In this section, we will consider modes with the same frequency and drop the subscript ω.By Theorem 4.2, we can search for a passive unitary that equally distributes the average occupation number N of a Gaussian state ρ among its modes: this unitary maximises the relative entropy of coherence of ρ.As a first case, let us consider a generic two-mode Gaussian state ρ, with mode operators â1 and â2 .We now apply a 50 : 50 beam splitter of phase φ (to be specified later): Then we have b This proof is similar to that of Theorem 4.2, and also here we drop the subscript ω in nω;j and N ω .
The covariance matrix of a Gaussian thermal state (see Eq. ( 30)) is diagonal and coincides with the diagonal matrix D in Williamson's theorem (see Eq. ( 11)).From Eq. ( 15 one can verify that this function is monotonically increasing with Ms .The minimum is therefore obtained for Ms = 1, i.e. when N is contained in a single reduced spatial mode, and the maximum is obtained for Ms = M s , i.e. when N is equally distributed.

E Proof of Eq. (60)
Let ÛE be the passive Gaussian unitary that achieves E G rel; max (ρ) in Eq. ( 59).Then we have Similarly, let ÛD be the Gaussian unitary that achieves D G rel; max (ρ) in Eq. ( 58).Then Finally, using C G rel (ρ max ) = P G rel (ρ) from Theorem 5.4, we obtain the desired result.

Figure 1 :
Figure 1: Graphical depiction of the labeling of the spectral and spatial modes.The modes are cataloged first in terms of their frequency ω = ω 1 , ω 2 , . . .ω M f (represented by a distinct colour) and then in terms of their spatial label j = 1, 2, . . .M s .

Definition 4 . 1 .
A state ρ max is a maximally coherent mixed state (MCMS) at fixed energy with respect to a coherence monotone C (see Def. 3.1) if

Figure 2 :
Figure 2: Graphical depiction of the relative entropy of non-uniformity P rel (black line), coherence C rel (red line), symmetric quantum discord D rel (green line) and entanglement E rel (blue line) for a quantum state ρ.The uniform state τ M (δ) is an element of the incoherent set I, which is a convex subset of the zero-discord set Z, which in turn is a non-convex subset of the separable set S.

Figure 3 :
Figure 3: Graphical depiction of Eq. (60).The dotted circle represents all the states that can be obtained from ρ via passive unitaries ÛO .The red line, connecting the uniform state τ M (δ) to the MCMGS, is the maximal Gaussian coherence C G rel; max (ρ) = P G rel (ρ).The green and blue lines are the maximal Gaussian discord D G rel; max and entanglement E G rel; max , respectively.The uniform state τ M (δ) is an element of the Gaussian incoherent set I G , which is a convex subset of the Gaussian zero-discord set Z G , which in turn is a non-convex subset of the Gaussian separable set S G .
j + 1) log(n j + 1) − nj log nj ] + λ Ms j=1 nj − λN, in this case M s = 2), thus maximising the coherence.Let us now consider a generic M s -mode product state ρ = 1 ⊗ • • • ⊗ Ms with d = 0. We prove that the quantum Fourier transform (QFT)âj → bj := Ms (j−1)(k−1) âj(73)is the passive Gaussian unitary that maximises the coherence.After the action of the QFT, the occupation number for the mode j readsMs (j−1)(k −k) â † k âk .(74)Weseparate the sum into two parts, with k = k and k = k :Ms (j−1)(k −k) â † k âk .(75)Since ρ is a product state and d = 0, we conclude the proof by noticing transformation is an extension of the DV unitary to CV.The unitary that maximises the coherence of an arbitrary DV state ρ for any MIO monotone reads[10] Ms (n−1)(k−1) |k ρ n | , (77)where d is the dimension of ρ, |ρ n are the eigenstates of ρ, and |k are the elements of the incoherent basis.Notice, however, that the CV result only holds for product states with d = 0.D Proof of Theorem 5.1 ), it follows that the von-Neumann entropy of τ Ms ( n) readsS (τ Ms ( n)) = Ms j=1 [(n j + 1) log(n j + 1) − nj log nj ] ,(78)where n = (n 1 , . . ., nMs ) and Ms is the number of modes for which nj = 0. Let L(n, λ) be the Lagrangian functionL(n, λ) :=S (τ Ms ( n)) − λ j + 1) log(n j + 1) − nj log nj ] + λ condition ∂L/∂ nj = 0 for nj = 0 is equivalent to λ = − log nj + 1 nj , ∀ j.(81)The above relations are satisfied by any state τ Ms ( n * ) with Ms reduced occupation numbers nj = N/ Ms ∀ j = 1, . . .Ms and the others equal to zero.Clearly 1 ≤ Ms ≤ M s .The entropy of any τ Ms ( n * ) reads:S (τ Ms ( n * )) = (N + Ms ) log N + Ms Ms − N log N Ms .(82)Taking the derivative of this expression with respect to Ms dS (τ Ms ( n * )) Ms j=1 t ω;j O ω;j .Definition 3.3.A function C G (ρ) is a suitable measure of Gaussian coherence with respect to IG (I G ) as free operations (free states), if it satisfies the properties (C1)-(C4) of Def.3.1.Here, IG are defined in Def.3.2 and states in I G are defined by Eq. (