Superior Resilience of Non-Gaussian Entanglement against Local Gaussian Noises

Entanglement distribution task encounters a problem of how the initial entangled state should be prepared in order to remain entangled the longest possible time when subjected to local noises. In the realm of continuous-variable states and local Gaussian channels it is tempting to assume that the optimal initial state with the most robust entanglement is Gaussian too; however, this is not the case. Here we prove that specific non-Gaussian two-mode states remain entangled under the effect of deterministic local attenuation or amplification (Gaussian channels with the attenuation factor/power gain κi and the noise parameter μi for modes i=1,2) whenever κ1μ22+κ2μ12<14(κ1+κ2)(1+κ1κ2), which is a strictly larger area of parameters as compared to where Gaussian entanglement is able to tolerate noise. These results shift the “Gaussian world” paradigm in quantum information science (within which solutions to optimization problems involving Gaussian channels are supposed to be attained at Gaussian states).


Introduction
Quantum entanglement is a feature of genuinely quantum correlations that underly many interesting physical phenomena and are a primary resource in quantum information theory enabling advantageous protocols of information processing and transmission [1]. A typical scenario to share this resource among distant parties is to prepare an entangled state in a laboratory and send entangled parts of the state toward the recipients via quantum communication lines [2]. Upon receiving the respective components of the entangled state the recipients can utilize them in entanglement-enabled protocols such as deviceindependent quantum key distribution [3][4][5] or further convert these entangled degrees of freedom into other degrees of freedom at their disposal, thus creating an entanglement among particles or modes that have never interacted before (entanglement swapping) [6]. The crucial practical challenge in this scenario is that the quantum communication lines are far from being ideal and introduce noise that degrades entanglement. The noise is local because the signals propagate to recipients via different communication lines whose environments do not interact with each other. There are some known techniques on how to locally compensate the negative effect of noise and distill some number of highly entangled states from many poorly entangled ones [7]. However, if the noise destroys entanglement completely so that the recipients actually receive a fully separable quantum state, no local strategy is able to cure this state and make it entangled. Therefore, there exists a fundamental limitation on the noise intensity exceeding which no entanglement can survive-a manifistation of the entanglement annihilation phenomenon [8].
The maximally permissible noise level for entanglement preservation apparently depends on the physical nature of the quantum states sent, the degrees of freedom used, and the nature of the noise that affects those degrees of freedom. In the case of two recipients A and B a local noise is mathematically described by a completely positive trace preserving map Φ A 1 ⊗ Φ B 2 . The quantum channels Φ 1 and Φ 2 can differ in general as the noises in the communication lines do not have to be identical due to, e.g., a different length or a different physical environment, say, atmosphere [9], optical fiber [10,11], photonic chip [12] or even a quantum memory cell, which can also be considered as a quantum channel. Historically, the first fully studied case of entanglement dynamics is the case of one-sided noise for which Φ 1 = Id, the identity transformation. In this case the local map Id A ⊗ Φ B 2 degrades entanglement of any state if and only if Φ 2 has a measure-andprepare (entanglement-breaking) structure [13][14][15]. For the one-sided noise, the optimal initial state for entanglement preservation is maximally entangled. The case of a two-sided noise Φ A 1 ⊗ Φ B 2 is more complicated to study even for simple systems (such as qubits) and relatively simple noise models (such as T 1 -T 2 qubit decoherence) because the optimal initial state for entanglement preservation is not known a priori. The entanglement annihilation parameters for a local two-qubit depolarizing noise were found in Ref. [16]. Local two-qubit unital noises were studied in Ref. [17] and then generalized with the help of the quantum Sinkhorn theorem to the case of an arbitrary two-qubit local noise Φ A 1 ⊗ Φ B 2 [18]. The case of trace decreasing qubit maps (corresponding to the loss of particles) was analyzed in Ref. [19]. The higher dimensional systems (qudits) were considered in Refs. [20][21][22], with a full characterization of entanglement annihilation parameters being obtained in the case of local qudit-depolarizing noises [22].
The continuous-variable quantum states of electromagnetic radiation (in contrast to discrete-variable polarization states) naturally incorporate the loss of photons within their description and encode quantum information in the field amplitude. The vacuum, coherent, squeezed, and thermal states are typical examples of Gaussian continuous-variable quantum states whose characteristic function is Gaussian [23][24][25]. These states are conventionally sensed by means of the homodyne and heterodyne measurements, which in turn are Gaussian quantum measurements [26]. The losses in quantum communication lines (potentially accompanied by an admixture of extra noise) are described by a quantum attenuator channel Φ(κ, µ), where κ ∈ [0, 1) is the intensity attenuation factor and µ ≥ 1 2 (1 − κ) is the noise parameter precisely defined in terms of characteristic functions in Section 2. Similarly, a deterministic phase-insensitive amplification of optical signal is described by a quantum linear amplifier Φ(κ, µ), where κ ∈ (1, ∞) is the power gain and µ ≥ 1 2 (κ − 1) is the noise parameter precisely defined in terms of characteristic functions in Section 2. For both channels the minimal noise (µ QL ≡ 1 2 |κ − 1|) corresponds to the quantum limited operation, the extra noise is a = µ − µ QL . The fact that the noise parameter µ > 0 is a consequence of the canonical commutation relation for the photon creation and annihilation operators which is to be respected by the deterministic attenuation and amplification. The channel Φ(1, 0) is the ideal channel (identity transformation), whereas Φ(1, µ) describes an addition of classical noise. All the described channels are Gaussian and transform any Gaussian state into another Gaussian state. So does the tensor product 2 ) destroys any entanglement between modes A and B (i.e., the channel Φ(κ 2 , µ 2 ) is entanglement breaking) if µ 2 ≥ 1 2 (κ 2 + 1). The region of parameters κ 1 , µ 1 , κ 2 , µ 2 , where the two-sided noise Φ A (κ 1 , µ 1 ) ⊗ Φ B (κ 2 , µ 2 ) annihilates entanglement is not fully characterized yet, and this is a goal of this paper to advance toward better understanding of this general case.
Given the fact that Gaussian states, Gaussian channels, and Gaussian measurements cover many practically relevant scenarios, it is not surprising the subfield of Gaussian quantum information emerged [27]. The subfield is often thought of to be closed in the sense that an optimization problem for a Gaussian channel or a Gaussian measurement should have a solution within the class of Gaussian states or the Gaussian ensembles of Gaussian states. This is true, e.g., with regard to the maximal reliable communication rate through a quantum channel Φ(κ, µ) (for all κ) [28] as well as the maximal communication rate by using a Gaussian homodyne or heterodyne measurement [29][30][31]. Therefore, it is very tempting to conjecture that the most robust entanglement with respect to Gaussian noises is exhibited by a Gaussian state [32] (with the experiments being focused on the Gaussian states too [33,34]). The claim of Ref. [35] is that if one adopts a bona fide measure of entanglement, which is continuous and strongly superadditive, such as the distillable entanglement; then at time t in the Markovian Gaussian evolution e Lt (with generator L), the bipartite Gaussian state e Lt [ G 12 (0)] with energyn(t) is the most entangled among all possible continuous-variable states e Lt [ 12 (0)] with the same energy, as a simple consequence of the maximum entropy property of Gaussian states [36,37]. We believe this claim of Ref. [35] is incorrect for the following reason. Ref. [36] correctly mentions that if only the covariance matrix of some density operator 1 is known, then there exists a Gaussian density operator G 1 with the same covariance matrix (hence, the same energy) such that the entropy of G 1 is greater than the entropy of 1 . If 1 and G 1 are the reduced density operators of pure bipartite states 12 and G 12 , then the state 12 is indeed less entangled than the state G 12 [37]. However, this line of reasoning cannot be generalized to the evolution states G 12 (t) and 12 (t) without breaking the equalities G 12 (t) = e Lt [ G 12 (0)] and 12 (t) = e Lt [ 12 (0)]. In other words, the maximum entropy consideration does not guarantee that the states G 12 (t) and 12 (t) are obtained from some legitimate states G 12 (0) and 12 (0) in the same dynamical evolution e Lt . In fact, Proposition 1 in Ref. [36] is in accordance with our findings (that we present in this paper) as it states that a continuous and strongly superadditive entanglement measure E satisfies the inequality E( G 12 ) ≤ E( 12 ) for every density operator 12 with finite covariance matrix and any Gaussian state G 12 with the same covariance matrix, i.e., Gaussian states give a lower bound for the entanglement E( 12 ). Finally, the conjecture on the superior robustness of Gaussian entanglement was shown to be false in 2011 by Sabapathy et al. [38]. They showed that none of two-mode Gaussian states can preserve entanglement under the Gaussian transformation Φ(κ, µ) ⊗ Φ(κ, µ) if µ ≥ 1 2 , whereas there exits a non-Gaussian state 12 (0) that remains entangled if κ > κ * ≈ 1 2 and µ = 1 2 . Moreover, for such parameters the preserved entanglement of the output state 12 (t) is distillable [i.e., E 12 (t) > 0 whereas E G 12 (t) = 0 for any initial Gaussian state G 12 (0) with an arbitrary energy] because for the entanglement detection Sabapathy et al. consider an effective two-qubit subspace (spanned by two-mode Fock states |00 , |0n , |n0 , |nn ) and any non-zero two-qubit entanglement is known to be distillable [39].
The result of Sabapathy et al. [38] was further strengthened in Ref. [40], where it was shown that the one-photon non-Gausian state |ψ AB * = 1 √ 2 (|01 − |10 ) outperforms all Gaussian states (with an arbitrarily high energy) and remains entangled when affected upon by a symmetric noise Φ(κ, µ) ⊗ Φ(κ, µ) with µ < 1 2 √ 1 + κ 2 for all κ. In the case of asymmetric local Gaussian noise Φ A (κ 1 , µ 1 ) ⊗ Φ B (κ 2 , µ 2 ), the question has remained open so far because the state |ψ AB * cannot generally outperform any Gaussian state in terms of entanglement robustness. It is not hard to see with the help of Simon's criterion [41] that the Gaussian channel Φ A (κ 1 , µ 1 ) ⊗ Φ B (κ 2 , µ 2 ) may preserve a two-mode Gaussian entanglement (i.e., entanglement of some two-mode Gaussian state) if and only if [40] The results of Ref. [40] indicate that the state can hardly remain entangled for all the parameters satisfying the inequality (1), thereby leaving an open problem on whether Gaussian states may potentially be optimal for entanglement distribution through asymmetric local Gaussian channels. Here we finally resolve this question in negative by demonstrating a channel-dependent non-Gaussian state | ψ AB c κ 1 ,µ 1 ,κ 2 ,µ 2 that outperforms any Gaussian state and remains entangled under the transformation which is a strictly larger region than (1). This result explicitly shows that the non-Gaussian entanglement is indeed more resilient to local Gaussian noises as compared to Gaussian states. Non-Gaussian entanglement can be useful not only in the entanglement distribution but also in other tasks, for example, the multi-component cat states were shown to outperform two-mode squeezed vacuum states in phase estimation [42]. The paper is organised as follows. In Section 2, we briefly review the formalism of characteristic functions and describe noisy attenuator and amplifier channels. In Section 3, the entanglement witness used is presented. Section 4 justifies the main result of the paper, the inequality (2), and presents a family of channel-dependent non-Gaussian states | ψ AB c κ 1 ,µ 1 ,κ 2 ,µ 2 . Summary of the results is given in Section 5.

Gaussian Attenuation and Amplification Channels
A continuous-variable quantum state can be alternatively described by a density operator (positive semidefinite operator with unit trace) or any of numerous phasespace functions, for instance, the characteristic function ϕ(z) = tr[ W(z)], where the Weyl operator W(z) for N modes is expressed through the canonical operators q i and p i (satisfying the canonical commutation relation [q i , p j ] = iδ ij ) through where (x 1 , y 1 , . . . , . This implies that the Gaussian state is fully determined by the first and second moments of the canonical operators. The Gaussian channel is a completely positive and trace preserving map that maps any Gaussian state into a Gaussian one and acts as follows in terms of input and output characteristic functions ϕ in (z) and ϕ out (z) [26]: where the real scaling matrix K and the real noise-quantifying symmetric matrix M must satisfy the inequality M ≥ i 2 (∆ − K ∆K) for the map to be completely positive. The full characterization of one-mode Gaussian channels is given in Ref. [43]. We focus on the most important (from the viewpoint of physical applications) deterministic Gaussian channels: phase-insensitive attenuators and amplifiers. These two types of Gaussian channels can be both described within unified formulas: where 0 ≤ κ < 1 for the attenuator, κ > 1 for the amplifier, and the noise parameter µ ≥ 1 2 |1 − κ| ≡ µ QL in both cases. We will denote such one-mode Gaussian channels by Φ(κ, µ). The minimal noise µ QL corresponds to the quantum-limited operation in which the admixed noise originates from vacuum fluctuations of the environment [44]; however, in some physically relevant systems the noise µ can be much higher that µ QL due to a hightemperature environment [45], for instance, in microwave quantum experiments [46][47][48].
The channel Φ(κ, µ) is known to be entanglement-breaking (so that the one-sided noise Φ(1, 0) ⊗ Φ(κ, µ) destroys any two-mode entanglement) if and only if the noise parameter µ exceed the threshold value, µ ≥ 1 2 (κ + 1) [49]. This value is depicted in Figure 1. Clearly, the quantum channel Φ(κ, µ) with fixed parameters κ and µ represents a snapshot of the dynamical map at a particular time moment t, which may correspond to a finite propagation time through a communication line. In a true time evolution the parameters κ and µ become functions of time t, κ(t) and µ(t). For instance, in the semigroup attenuation or amplification dynamics Φ = e Lt with the generator L ( [50], Section 3.4.6) we obtain the following dependencies: where the sign + (−) describes amplification (attenuation), Γ ≥ 0 is the process rate, and n is the average number of thermal photons in the environment. For such a semigroup dynamics the one-parameter family of maps Φ κ(t), µ(t) is associated with a straight line in the parameter space (κ, µ), see Figure 1.
The diagonal sum representation for the channel Φ(κ, µ) is found in Ref. [51] in terms of the Fock states {|n } ∞ n=0 and in Ref. [40] in terms of the coherent states {|α } α∈C . We use the latter one as it sheds more light on the phase space picture [52,53]. The integral form of the diagonal sum representation reads where the Kraus operators A αβ (κ, µ) rather nontrivially depend on parameters κ and µ through auxiliary expressions The peculiar form of the presented Kraus operators originates from the concatenation formula Φ(κ, µ) = Φ QLτ • Φ QLη that shows the channel Φ(κ, µ) can always be thought of as a sequential application of the quantum limited attenuator Φ QLη (with the attenuation factor η) and the quantum limited amplifier Φ QLτ (with the power gain τ). Using the diagonal sum representation above it is not hard to see the effect of the channel Φ(κ, µ) on the outer product of coherent states, namely, Since any density operator can be decomposed as = π −2 d 2 γd 2 δ γ| |δ |γ δ|, Equations (9) and (10) define the effect of the Gaussian channel Φ(κ, µ) on any state (including non-Gaussian ones).

Entanglement Witness
A bipartite state AB is called entangled if it cannot be represented by a convex sum of density operators A k and B k for individual subsystems A and B. Otherwise the state is separable. In our case each subsystem is a mode of electromagnetic radiation, so we analyze a two-mode entanglement. Although simple and physically clear entanglement criterions are known for Gaussian two-mode systems [41,54], they are still missing for non-Gaussian states. For this reason we resort to an entanglement witness formalism that enables us to detect entanglement within non-Gaussian states too [40].
Consider a Hermitian operator which reduces to a swapping operator if λ = 0. The important feature of this operator with λ ∈ R is that its average value is nonnegative in any separable state because for all vectors |φ A and |χ B . Therefore, W λ is an entanglement witness in the sense that the negativity in the average value tr[ AB W λ ] can only arise from the entangled state AB , so tr[ AB W λ ] < 0 unambiguously indicates the entanglement between subsystems A and B. If λ > 0, then W λ is an unbounded operator; however, the average value tr[ AB W λ ] may still be finite. Suppose the local Gaussian channel Φ A (κ 1 , µ 1 ) ⊗ Φ B (κ 2 , µ 2 ) acts on some two-mode pure state |ψ AB ψ AB |, then the output state is verified to be entangled if tr[ AB W λ ] < 0 for some λ ∈ R. One can actually adjust the input state |ψ AB ψ AB | to be dependent on the channel parameters κ 1 , µ 1 , κ 2 , µ 2 in order to find a robust entangled state for a given channel.

Robust Non-Gaussian Entanglement
Here we present a non-Gaussian state |ψ AB c ψ AB c |, which exhibits high resilience to the local Gaussian noise Φ A (κ 1 , µ 1 ) ⊗ Φ B (κ 2 , µ 2 ). A real scaling coefficient c ∈ R describes asymmetry in the skewed two-mode "cat state" where |γ and |cγ are coherent states and |0 is a vacuum state. At the end of this section we will make the coefficient c dependent on parameters κ 1 , µ 1 , κ 2 , µ 2 so as to present the noise-specific resilient entangled state; however, for now we treat it as an independent parameter. If c = 1, then we get a symmetric state whose entanglement is robuster against symmetric Gaussian noise Φ A (κ, µ) ⊗ Φ B (κ, µ) than that of any Gaussian state [40]. When we deal with asymmetric noise Φ A (κ 1 , µ 1 ) ⊗ Φ B (κ 2 , µ 2 ), the robustest input state is expected to be asymmetric too and this is the physical reason to consider the skewed state (13).
Each of the operators (14)- (17) is readily calculated via Equation (9) and then used to calculate the Hilbert-Schmidt scalar product with the entanglement witness W λ : The straightforward integral calculations show that all four expressions are finite if Suppose γ → 0, then it is not hard to see via the Taylor series that the condition is fulfilled if the following inequality is met for the leading term (proportional to |γ| 2 ): where τ i , η i are expressed through κ i , µ i via Equation (7), i = 1, 2. Some elementary algebra yields that the inequality (24) is fulfilled for some λ satisfying (22) if Since the limit γ → 0 has been used in the derivation above, we actually get the one-photon continuous-variable state where |1 is a single-photon Fock state. Inequalities (28) give sufficient conditions under which the output state is entangled. These inequal-ities can be rewritten in terms of the additional noise on top of the quantum-limited operation, a i = µ i − 1 2 |κ i − 1| ≥ 0, i = 1, 2, as the following result. µ 1 ) and Φ B (κ 2 , µ 2 ) are attenuators (κ 1,2 ≤ 1) and • Φ A (κ 1 , µ 1 ) is an attenuator (κ 1 ≤ 1) and Φ B (κ 2 , µ 2 ) is an amplifier (κ 2 ≥ 1) and • both Φ A (κ 1 , µ 1 ) and Φ B (κ 2 , µ 2 ) are amplifiers (κ 1,2 ≥ 1) and If c = 1 and we deal with a symmetric state | ψ 1 = 1 √ 2 (|1 |0 − |0 |1 ), then Proposition 1 reproduces the known results from Ref. [40]; however, if c = 1, then Proposition 1 gives a valuable generalization. In Figure 2, we depict typical noise levels tolerated by the state | ψ c without losing entanglement. If we choose specific values of κ 1 , κ 2 , and c, then the shaded region in Figure 2 shows the extra noises a 1 , a 2 , which do not destroy entanglement of | ψ c . One can clearly see that the non-Gaussian state | ψ c outperforms all Gaussian ones (dashed line) if we deal with a two-side (asymmetric or symmetric) amplification ( Figure 2c); however, a single state | ψ c with some fixed value of c cannot outperform all Gaussian states (dashed line) in the case of asymmetric attenuation (κ 1,2 < 1, κ 1 > κ 2 , Figure 2d) or in the case of amplification-attenuation channel (κ 1 > 1 > κ 2 , Figure 2b). If a 2 tends to zero, then Φ B (κ 2 , µ 2 ) approaches the quantum limited attenuation and the high-energy Gaussian entanglement tolerates this noise. Dashed line in Figure 2 corresponds to the infinite-energy two-mode squeezed state [40]. Nonetheless, varying the parameter c, we cover the area in the extra noise space (a 1 , a 2 ) which is strictly larger than the area for the Gaussian states (dashed line, Equation (1)) and fully comprises the latter. Some elementary algebra shows that the envelope curve for the former area (the arc of an ellipsoid in Figure 2) is exactly given by the equation κ 1 µ 2 2 + κ 2 µ 2 1 = 1 4 (κ 1 + κ 2 )(1 + κ 1 κ 2 ). Therefore, if the inequality (2) is fulfilled, then there exists some c = c κ 1 µ 1 κ 2 µ 2 such that the Example 1. Consider the Gaussian local attenuation dynamics e L 1 t ⊗ e L 2 t with the time-dependent parameters (5), where the rates Γ 1 = Γ 2 = Γ andn 1 = 0.5,n 2 = 1.5. Then In such dynamics, the fundamental threshold (1) implies that no Gaussian entanglement can survive if Γt ≥ − ln 2 3 ≈ 0.405. On the other hand, the inequality (2) implies that the non-Gaussian entanglement of states (13) may be preserved up to the threshold Γt = − ln 10− √ 19 9 ≈ 0.467. In Figure 3, we show both thresholds as well as the dynamics of entanglement for Gaussian and non-Gaussian states. Non-Gaussian entanglement is characterized by the introduced entanglement witness, namely, the quantity −tr[W λ (t)] with λ = 0.74. As initial states we use the non-Gaussian state (13) with γ = 1 and c = 2, which has the energy 2.7 photons, and the single-photon non-Gaussian state (29) with c = 2. Figure 3 suggests that these non-Gaussian states remain entangled after passing the Gaussian threshold, with the single-photon state exhibiting stronger resilience for the entanglement measure used. The latter fact is in agreement with the earlier observations on the low-energy entanglement robustness [37]. (Non-monotonic behavior of our entanglement measure for the higher-energy state should not be surprising as the witness operator W λ assigns higher values to coherent states with higher amplitude. As the noise increases, the overlap with high-amplitude states increases too. In this sense, our entanglement measure is not an entanglement monotone but it is rather a very sensitive entanglement detector for slightly entangled states.) As examples of Gaussian states we consider two-mode squeezed vacuum states with the same initial energy (2.7 photons and 1 photon). Unfortunately, our entanglement measure diverges for these states so we resort to the Simon criterion (necessary and sufficient condition for two-mode Gaussian entanglement [41]), namely, the entanglement measure is the minimal eigenvalue of the matrix V(p 2 → −p 2 ) − i 2 ∆, where V is the covariance matrix and V(p 2 → −p 2 ) is its modified version obtained via inversion of the momentum operator for the second mode. Figure 3 illustrates that the entanglement of such Gaussian states vanishes before the Gaussian threshold. The entanglement death time approaches the Gaussian threshold if the energy of the two-mode squeezed state tends to infinity. Whatever low the non-Gaussian entanglement can be after passing the Gaussian threshold, Ref. [38] suggests that this entanglement is distillable in general, so it can be converted into a useful form. In principle, there may exist non-Gaussian states that remain entangled even after the established non-Gaussian threshold, because this threshold was derived for a specific family of non-Gaussian states {|ψ AB c } and a specific entanglement measure −tr[W λ (t)].  To make the analysis complete and more applicable to practice we specify how the advantageous parameter c κ 1 µ 1 κ 2 µ 2 should be chosen for given channel parameters κ 1 , µ 1 , κ 2 , µ 2 . Suppose the channel parameters κ 1 , µ 1 , κ 2 , µ 2 satisfy the inequality (2), then we fix the amplification/attenuation parameters κ 1,2 and find the maximally permissible noise levels µ 1,2 by projecting the point (µ 1 , µ 2 ) onto the boundary of the ellipsoid (see Figure 4), namely, Then the advantageous parameter c κ 1 µ 1 κ 2 µ 2 is chosen in such a way that a corner point for the noise-toleration figure coincides with the point ( a 1 , a 2 ) in Figure 4, a i = µ i − 1 2 |κ i − 1|. The algebraic expression reads The obtained advantageous state | ψ AB c κ 1 µ 1 κ 2 µ 2 has an interesting feature regarding the distribution of energy between the modes ( 1 1+c 2 quanta in the first mode, c 2 1+c 2 quanta in the second mode). In Figure 2b we observe that in the case when the first mode is amplified and the second one is attenuated, sometimes it is beneficial to put more energy in the first (amplified) mode to protect entanglement. Another especially counterintuitive result takes place if both modes are asymmetrically amplified (κ 1,2 ≥ 1) and κ 2 > κ 1 + 2: then c κ 1 µ 1 κ 2 µ 2 > 1 irrespective to the noise admixed, i.e., it is always beneficial to put more energy in the second mode with a higher power gain to preserve entanglement. , whose entanglement is resilient against a given extra noise (a 1 , a 2 ) (green dot) in the two-mode Gaussian attenuator/amplifier with the attenuation factor/power gain κ 1 , κ 2 in the corresponding mode. Straight solid line passes through the ellipsoid center and the green dot, thus intersecting the ellipsoid (2) in the point ( a 1 , a 2 ) that determines c κ 1 µ 1 κ 2 µ 2 via Formula (34).

Conclusions
The current study shows that the entanglement of specific low-energy non-Gaussian states has stronger resilience against local Gaussian noises as compared to that of arbitraryenergy Gaussian ones. A similar observation was made earlier in the case of symmetric Gaussian noise Φ A (κ, µ) ⊗ Φ B (κ, µ); however, there was a loophole in the case of asymmetric Gaussian noises Φ A (κ 1 , µ 1 ) ⊗ Φ B (κ 2 , µ 2 ). In this study, we have completely closed this loophole by a rigorous proof on the existence of an entangled non-Gaussian state whose entanglement can tolerate as high noises µ 1 and µ 2 as meet the requitement (2). The Guassian states do lose entanglement beyond the region (1), which is a subset of (2). Moreover, we have explicitly characterized the advantageous non-Gaussian single-photon skewed state | ψ c κ 1 µ 1 κ 2 µ 2 that exhibits resilient entanglement and discussed peculiarities of the energy distribution among the modes in this state. Our results imply that the optimization problem of maximizing time t during which the state Φ A (κ 1 (t), µ 1 (t)) ⊗ Φ B (κ 2 (t), µ 2 (t))[|ψ AB ψ AB |] remains entangled definitely has a solution |ψ AB beyond the set of Gaussian states. Conceptually, these results stimulate further research beyond the "Gaussian world" paradigm in quantum information science (within which solutions to optimization problems involving Gaussian channels and Gaussian measurements are supposed to be attained at Gaussian states and Gaussian ensembles).