Böttcher-Wenzel inequality for weighted Frobenius norms and its application to quantum physics

Abstract

By employing a weighted Frobenius norm with a positive matrix $\omega$, we introduce natural generalizations of the famous Böttcher-Wenzel (BW) inequality. Based on the combination of the weighted Frobenius norm $\|A\|_{\omega}:=\sqrt{\tr(A^{\ast}A\omega)}$ and the standard Frobenius norm $\|A\|:=\sqrt{\tr(A^{\ast}A)}$, there are exactly five possible generalizations, labeled (i) through (v), for the bounds on the norms of the commutator $[A,B]:=AB-BA$. In this paper, we establish the tight bounds for cases (iii) and (v), and propose conjectures regarding the tight bounds for cases (i) and (ii). Additionally, the tight bound for case (iv) is derived as a corollary of case (i). All these bounds (i)-(v) serve as generalizations of the BW inequality. The conjectured bounds for cases (i) and (ii) (and thus also (iv)) are numerically supported for matrices up to size $n=15$. Proofs are provided for $n=2$ and certain special cases. Interestingly, we find applications of these bounds in quantum physics, particularly in the contexts of the uncertainty relation and open quantum dynamics.

1 Introduction

The seminal Böttcher-Wenzel (BW) inequality [1, 2, 3] provides the bound of the norm of the commutator $[A,B]:=AB-BA$ of $n\times n$ complex matrices $A,B\in M_{n}(\mathop{\mathbb{C}}\nolimits)$ in the form

$$\|[A,B]\|\leq\sqrt{2}\|A\|\|B\|,$$ (1)

where $\|A\|:=\sqrt{\tr(A^{\ast}A)}$ is the Frobenius norm. Here $\tr A$ and $A^{\ast}$ are the trace and the Hermitian conjugate of the matrix $A$, respectively. The bound (1) is tight in the sense that there exist non-zero matrices $A,B$ that attain the equality. This inequality was then generalized in several directions, e.g., with Schatten $p$-norm, Ky Fan $(p,k)$ norm [4, 5, 6], or with the $q$-deformed commutator [7, 8].

In this paper we consider generalizations of BW inequality by replacing the Frobenius norm defined in terms of the Hilbert Schmidt inner product $\langle A|B\rangle:=\tr(A^{\ast}B)$, i.e. $\|A\|^{2}=\langle A|A\rangle$, by $\|A\|_{\omega}^{2}:=\langle A|A\rangle_{\omega}$, where the new inner product is defined as follows $\langle A|B\rangle_{\omega}:=\tr(A^{*}B\omega)$ and $\omega\in M_{n}(\mathop{\mathbb{C}}\nolimits)$ is a positive (definite) matrix. In what follows we call $\omega$-weighted Frobenius norm

$$\|A\|_{\omega}:=\sqrt{\tr(A^{\ast}A\omega)}$$ (2)

the $\omega$-norm. The $\omega$-norm satisfies the axioms of the norm and provides a generalization of the Frobenius norm when $\omega$ is the identity matrix $\mathop{\mathbb{I}}\nolimits$. In the following discussion, $\omega$ is always assumed to be a positive matrix.

Depending on the combinations of the $\omega$-norm and the Frobenius norm, there are exactly six types of bounds to be considered, including the BW inequality itself: For positive constants $c_{i}(\omega)\ (i=1,\ldots,6)$, dependent on $\omega$, it holds that, for any $A,B\in M_{n}(\mathop{\mathbb{C}}\nolimits)$,

$$\rm{(i)}\ \|[A,B]\|_{\omega} \leq c_{1}(\omega)\|A\|_{\omega}\|B\|_{\omega},$$ (i)

$$\rm{(ii)}\ \|[A,B]\|_{\omega} \leq c_{2}(\omega)\|A\|_{\omega}\|B\|,$$ (ii)

$$\rm{(iii)}\ \|[A,B]\|_{\omega} \leq c_{3}(\omega)\|A\|\|B\|,$$ (iii)

$$\rm{(iv)}\ \|[A,B]\| \leq c_{4}(\omega)\|A\|_{\omega}\|B\|_{\omega},$$ (iv)

$$\rm{(v)}\ \|[A,B]\| \leq c_{5}(\omega)\|A\|_{\omega}\|B\|,$$ (v)

$$\rm{(vi)}\ \|[A,B]\| \leq c_{6}(\omega)\|A\|\|B\|.$$ (vi)

In particular, we are interested in the tightest bounds $\tilde{c_{i}}(\omega)\ (i=1,\ldots,6)$, which are the minimum values of $c_{i}(\omega)$ ($i=1,\ldots,6$) for which each inequality in Eqs. (3) is satisfied. Alternatively, $\tilde{c_{i}}(\omega)$ can be characterized by the following optimization problems:

$$\tilde{c_{1}}(\omega) =\max_{A,B\neq 0\in M_{n}(\mathop{\mathbb{C}}\nolimits)}\frac{\|[A,B]\|_{\omega}}{\|A\|_{\omega}\|B\|_{\omega}},$$ (i)

$$\tilde{c_{2}}(\omega) =\max_{A,B\neq 0\in M_{n}(\mathop{\mathbb{C}}\nolimits)}\frac{\|[A,B]\|_{\omega}}{\|A\|_{\omega}\|B\|},$$ (ii)

$$\tilde{c_{3}}(\omega) =\max_{A,B\neq 0\in M_{n}(\mathop{\mathbb{C}}\nolimits)}\frac{\|[A,B]\|_{\omega}}{\|A\|\|B\|},$$ (iii)

$$\tilde{c_{4}}(\omega) =\max_{A,B\neq 0\in M_{n}(\mathop{\mathbb{C}}\nolimits)}\frac{\|[A,B]\|}{\|A\|_{\omega}\|B\|_{\omega}},$$ (iv)

$$\tilde{c_{5}}(\omega) =\max_{A,B\neq 0\in M_{n}(\mathop{\mathbb{C}}\nolimits)}\frac{\|[A,B]\|}{\|A\|_{\omega}\|B\|},$$ (v)

$$\tilde{c_{6}}(\omega) =\max_{A,B\neq 0\in M_{n}(\mathop{\mathbb{C}}\nolimits)}\frac{\|[A,B]\|}{\|A\|\|B\|}.$$ (vi)

Note that, given inequalities (3), the tightness can be also shown by providing non-zero matrices $A,B$ such that the equalities are attained.

We remark that type (vi) corresponds to BW bound (1), so that we already know

$$\tilde{c_{6}}(\omega)=\sqrt{2}.$$ (5)

In this paper, we first show that the tight bounds for cases (iii) and (v) are give by

$$\tilde{c_{3}}(\omega)=\sqrt{2\lambda_{\rm M}},\ \tilde{c_{5}}(\omega)=\sqrt{\frac{2}{\lambda_{\rm m}}},$$ (6)

where $\lambda_{\rm M}$ and $\lambda_{\rm m}$ are the largest and the smallest eigenvalues of $\omega$. Second, we give conjectures for cases (i) and (ii) that the tight bounds are given respectively by

$$\tilde{c_{1}}(\omega)=\sqrt{\frac{\lambda_{\rm m}+\lambda_{\rm sm}}{\lambda_{\rm m}\lambda_{\rm sm}}}$$ (7)

and

$$\tilde{c_{2}}(\omega)=\sqrt{\frac{\lambda_{\rm m}+\lambda_{\rm M}}{\lambda_{\rm m}}},$$ (8)

where $\lambda_{\rm sm}$ is the second smallest eigenvalue of $\omega$. Third, the tight bound for case (iv) is then given, as a corollary of conjecture (7), by

$$\tilde{c_{4}}(\omega)=\sqrt{\frac{\lambda_{\rm m}+\lambda_{\rm sm}}{\lambda_{\rm m}^{2}\lambda_{\rm sm}}}.$$ (9)

Note that all these tight bounds for cases (i)-(v) generalize the BW inequality by setting $\omega=\mathop{\mathbb{I}}\nolimits$, since the eigenvalues of $\mathop{\mathbb{I}}\nolimits$ are all 1, allowing (6), (7), (8), and (9) to recover the BW bound (5).

To summarize, the logical relations between the tight bounds for cases (i)-(vi) are as follows:

In this regard, it becomes evidence that the verification of conjectures (i) and (ii) holds considerable importance. Furthermore, besides the intrinsic mathematical interest in the norm bounds of the commutator, we find that there are direct applications of both (i) and (ii) in the filed of quantum physics. By applying a quantum state, i.e., a density operator as $\omega$, we demonstrate that the conjectured bounds for type (i) and type (ii) respectively introduce a novel uncertainty relation between observables $A$ and $B$, and impose a non-trivial constraint on relaxation rates in the general quantum Markovian dynamics.

In the following discussion, we will frequently use the Dirac notation, as commonly used in the field of quantum physics: a vector in $\mathop{\mathbb{C}}\nolimits^{n}$ is denoted by a ket vector, e.g., $\ket{\psi}$; For vectors $\ket{\psi}=(x_{i})_{i=1}^{n},\ket{\phi}=(y_{i})_{i=1}^{n}\in\mathop{\mathbb{C}}\nolimits^{n}$, the symbol $\langle\psi|\phi\rangle:=\sum_{i}\overline{x_{i}}y_{i}$ denote the (complex Euclidean) inner product, while the symbol $\outerproduct{\psi}{\phi}:=[x_{i}\overline{y}_{j}]_{i,j}$ denotes the matrix in $M_{n}(\mathop{\mathbb{C}}\nolimits)$ such that its action to a vector is given by: $\outerproduct{\psi}{\phi}\ket{\xi}=\langle\phi|\xi\rangle\ket{\psi}$.

The structure of this paper is as follows: In Sec. 2 (and also in Appendix), we explore several generalizations of the BW inequality, each characterized by a combination of the $\omega$-norm and the Frobenius norm. Sec. 3 demonstrates the practical applications of our conjectures within the realm of quantum physics. The paper concludes with Sec. 4.

2 Six types of bounds of the commutators with respect to $\omega$-norm

To derive bounds for types (i) through (v) in (3), we start with a simple observation about the relationship between the $\omega$-norm and the Frobenius norm. In what follows, we denote by $\lambda_{i}>0$ $(i=1,\ldots,n)$ the eigenvalues of a positive matrix $\omega\in M_{n}(\mathbb{C})$, arranged in ascending order: $\lambda_{1}\leq\lambda_{2}\leq\cdots\leq\lambda_{n}$, so that

$$\lambda_{\rm m}=\lambda_{1},\lambda_{\rm sm}=\lambda_{2}\ {\rm and}\ \lambda_{\rm M}=\lambda_{n}.$$

The corresponding unit eigenvectors are denoted by $\ket{\lambda_{i}}$ $(i=1,\ldots,n)$, i.e., $\omega\ket{\lambda_{i}}=\lambda_{i}\ket{\lambda_{i}}$ and $\langle\lambda_{i}|\lambda_{i}\rangle=1$. Considering the expressions $\|A\|^{2}_{\omega}=\tr(A^{\ast}A\omega)=\sum_{i=1}^{n}\lambda_{i}\langle\lambda_{i}|A^{\ast}A\lambda_{i}\rangle$ and $\|A\|^{2}=\tr(A^{\ast}A)=\sum_{i=1}^{n}\langle\lambda_{i}|A^{\ast}A\lambda_{i}\rangle$, it follows that

$$\lambda_{\rm m}\|A\|^{2}\leq\|A\|^{2}_{\omega}\leq\lambda_{\rm M}\|A\|^{2}.$$ (10)

By utilizing these inequalities, we derive that

$$\|[A,B]\|_{\omega}^{2}\leq\lambda_{\rm M}\|[A,B]\|^{2}\leq 2\lambda_{\rm M}\|A\|^{2}\|B\|^{2}\leq\frac{2\lambda_{\rm M}}{\lambda_{\rm m}}\|A\|_{\omega}^{2}\|B\|^{2}\leq\frac{2\lambda_{\rm M}}{\lambda_{\rm m}^{2}}\|A\|_{\omega}^{2}\|B\|_{\omega}^{2},$$

where we have used the BW inequality (1) to estimate the second inequality. Now the inequalities clearly imply bounds (i)-(v) in (3) for

$$c_{1}(\omega)=\frac{\sqrt{2\lambda_{\rm M}}}{\lambda_{\rm m}},c_{2}(\omega)=\sqrt{\frac{2\lambda_{\rm M}}{\lambda_{\rm m}}},c_{3}(\omega)=\sqrt{2\lambda_{\rm M}},c_{4}(\omega)=\frac{\sqrt{2}}{\lambda_{\rm m}},c_{5}(\omega)=\sqrt{\frac{2}{\lambda_{\rm m}}}.$$ (11)

It turns out that this simple observation yields the tight bounds for cases (iii) and (v), since there are non-zero matrices $A,B$ that attain the equalities of (iii) and (v) with $c_{3}(\omega)$ and $c_{5}(\omega)$ in (11). Examples of such matrices include, for instance, $A=\outerproduct{\lambda_{1}}{\lambda_{n}}$ and $B=\frac{1}{\sqrt{2}}(\outerproduct{\lambda_{n}}{\lambda_{n}}-\outerproduct{\lambda_{1}}{\lambda_{1}})$ for (iii) and $A=\outerproduct{\lambda_{n}}{\lambda_{1}}$ and $B=\outerproduct{\lambda_{1}}{\lambda_{n}}$ for (v).

Put differently, we have shown the following:

Proposition 1

Let $\omega\in M_{n}(\mathop{\mathbb{C}}\nolimits)$ be a positive matrix. For any matrices $A,B\in M_{n}(\mathop{\mathbb{C}}\nolimits)$,

$$\|[A,B]\|_{\omega}\leq\sqrt{2\lambda_{\rm M}}\|A\|\|B\|$$ (12)

and

$$\|[A,B]\|\leq\sqrt{\frac{2}{\lambda_{\rm m}}}\|A\|_{\omega}\|B\|,$$ (13)

where $\lambda_{\rm M},\lambda_{\rm m}$ are the largest and the smallest eigenvalues of $\omega$. Both bounds are tight, i.e., there are non-zero matrices $A$ and $B$ that attain the equalities.

Meanwhile, we have conducted numerical optimizations of (4) for cases (i), (ii), and (iv), and the results suggest that the bounds $c_{1}(\omega),c_{2}(\omega)$ and $c_{4}(\omega)$ in (11) are not tight (See dashed lines in Fig. 1). Instead, we conjecture for the tight bounds as specified in (7), (8), and (9). In other words, our conjectures for cases (i), (ii) and (iv) read:

Conjecture 1

For any matrices $A,B\in M_{n}(\mathop{\mathbb{C}}\nolimits)$,

$$\|[A,B]\|_{\omega}\leq\sqrt{\frac{\lambda_{\rm m}+\lambda_{\rm sm}}{\lambda_{\rm m}\lambda_{\rm sm}}}\|A\|_{\omega}\|B\|_{\omega},$$ (14)

$$\|[A,B]\|_{\omega}\leq\sqrt{\frac{\lambda_{\rm m}+\lambda_{\rm M}}{\lambda_{\rm m}}}\|A\|_{\omega}\|B\|,$$ (15)

and (as a corollary of (14))

$$\|[A,B]\|\leq\sqrt{\frac{\lambda_{\rm m}+\lambda_{\rm sm}}{\lambda_{\rm m}^{2}\lambda_{\rm sm}}}\|A\|_{\omega}\|B\|_{\omega}.$$ (16)

All bounds are tight, i.e., there are non-zero matrices $A$ and $B$ that attain the equalities.

Some remarks are in order: First, these conjectures are strongly supported by numerical evidence: For randomly generated positive matrices $\omega$, we have confirmed that numerical optimizations of (4) for cases (i), (ii), and (iv) perfectly match the conjectured bounds (7), (8), and (9) up to size $n=15$. For illustration, consider a one parameterized matrix

$$\omega(p)=\text{diag}[\sin(2p),\sin(2p^{2}),\sin(2p^{3}),\sin(2p^{4}),\sin(2p^{5})]\in M_{5}(\mathbb{C})$$

with parameter $p\in(0,1]$ as a positive matrix $\omega$. (This matrix is just one example of $\omega$, but note that due to the positivity of $\omega$, it can always be diagonalized, so considering $\omega$ as a diagonal matrix does not result in any loss of generality.) Figure 1 illustrates the comparisons between numerical optimizations (dotted points, red in the online version) and the conjectured bounds (solid lines), as well as the loose bounds (dashed lines) given in (11).

Figure 1: As an illustration, using for a $5\times 5$ positive matrix $\omega(p)={\rm diag}[\sin(2p),\sin(2p^{2}),\sin(2p^{3}),\sin(2p^{4}),\sin(2p^{5})]$ for $p\in(0,1]$, we compare numerical optimizations (dotted points/red online) with our conjectures. In graphs (i), (ii), (iv), conjectured bounds (7), (8), and (9) are plotted with solid lines, while bounds given in (11) are plotted with dashed lines. For references, in graphs (iii) and (v), the tight bounds given in (11) are plotted with dashed lines.

Second, one can find non-zero matrices $A,B$ that attain equalities in (14), (15), and (16): For type (i), let, for instance, $A=\outerproduct{\lambda_{\rm m}}{\lambda_{\rm sm}}$ and $B=A^{\ast}$. Then, it holds that $\|[A,B]\|^{2}_{\omega}=\lambda_{\rm m}+\lambda_{\rm sm},\|A\|^{2}_{\omega}=\lambda_{\rm sm},\|B\|^{2}_{\omega}=\lambda_{\rm m}$ and the equality in (14) is attained. Similarly, attainability can be shown by setting, for type (ii), $A=\outerproduct{\lambda_{\rm M}}{\lambda_{\rm m}}$, $B=A^{\ast}$, and for type (iv), $A=\lambda_{\rm m}\outerproduct{\lambda_{\rm sm}}{\lambda_{\rm sm}}-\lambda_{\rm sm}\outerproduct{\lambda_{\rm m}}{\lambda_{\rm m}},B=\outerproduct{\lambda_{\rm sm}}{\lambda_{\rm m}}$, respectively. These facts imply that the proofs of the part concerning the tightness of the bounds are complete. Third, (16) can be obtained as a corollary of (14) by using the left inequality in (10) for the commutator part. Forth, we provide proofs of the conjectures for some special cases (see A); At the end of this section, we give proofs for case $n=2$. Last but not least, we emphasize that the forms of our conjectures are far from trivial. For instance, the fact that the tight bounds depend solely on the eigenvalues of $\omega$, and more specifically, on its largest, smallest, or second smallest eigenvalues, as presented in forms (7), (8), (9), was completely beyond our initial expectations. We arrived at these conjectured forms only after conducting numerous numerical experiments and trials.

2.1 Proofs of Conjecture 1 for $n=2$

Here we prove inequalities (14) and (15), hence Conjecture 1, for the case $n=2$. Before addressing the specific case of $n=2$, we begin with several general observations applicable to matrices of any size $n$. First, in order to prove the conjecture for types (i) and (ii) (indeed also for all other cases), we note that, without loss of generality, $\omega$ can be assumed to be a diagonal matrix with unit trace. This follows from the facts that, for any unitary matrix $U$, $\|UAU^{\ast}\|_{U\omega U^{\ast}}=\|A\|_{\omega}$ (thus $\|UAU^{\ast}\|=\|A\|$) and $[UAU^{\ast},UBU^{\ast}]=U[A,B]U^{\ast}$. Moreover, for any positive constant $p$, $\|A\|_{p\omega}=\sqrt{p}\|A\|_{\omega}$, and therefore $\sqrt{p}\tilde{c_{1}}(p\omega)=\tilde{c_{1}}(\omega)$ and $\tilde{c_{2}}(p\omega)=\tilde{c_{2}}(\omega)$ for (7) and (8). Note that these are also true for (i) and (ii).

Second, note that $M_{n}(\mathbb{C})$ constitutes an inner product space with respect to the $\omega$-inner product, defined as $\langle A|B\rangle_{\omega}:=\tr(A^{*}B\omega)$, where the $\omega$-norm is induced from it. As the identity matrix is a unit vector if $\omega$ has a unit trace, we have

$$\|A\|^{2}_{\omega}\geq\|A-\langle I|A\rangle_{\omega}I\|^{2}_{\omega}=\|A-\tr(A\omega)I\|^{2}_{\omega}.$$ (17)

Applying $\omega=\mathop{\mathbb{I}}\nolimits/n$, we also have

$$\|A\|^{2}\geq\|A-\frac{\tr(A)}{n}I\|^{2}.$$ (18)

Considering these observations, we now proceed to prove (14) and (15) in the case where $n=2$. In what follows, let $\omega={\rm diag}[\lambda_{1},\lambda_{2}]$ with $\lambda_{1}+\lambda_{2}=1$ and $\lambda_{1}\leq\lambda_{2}$, and $A=[a_{ij}]_{i,j},B=[b_{ij}]_{i,j}\in M_{2}({\mathbb{C}})$.

[Proof of (14)] Note that inequality (17) along with the fact $[A+\alpha\mathop{\mathbb{I}}\nolimits,B+\beta\mathop{\mathbb{I}}\nolimits]=[A,B]$ for any $\alpha,\beta\in{\mathbb{C}}$ allows us to assume that $A$ and $B$ take the forms $A-\tr(A\omega)I\ {\rm and}\ B-\tr(B\omega)I$, respectively. This is equivalent to assuming $a_{22}=-\frac{\lambda_{1}}{\lambda_{2}}a_{11}$ and $b_{22}=-\frac{\lambda_{1}}{\lambda_{2}}b_{11}$. The direct computation then yields the identity

$$\frac{\lambda_{\rm m}+\lambda_{\rm sm}}{\lambda_{\rm m}\lambda_{\rm sm}}\|A\|^{2}_{\omega}\|B\|^{2}_{\omega}-\|[A,B]\|^{2}_{\omega}=\left|\frac{1}{\lambda_{2}}\sqrt{\frac{\lambda_{1}}{\lambda_{2}}}a_{11}\bar{b}_{11}+\sqrt{\frac{\lambda_{1}}{\lambda_{2}}}a_{21}\bar{b}_{21}+\sqrt{\frac{\lambda_{2}}{\lambda_{1}}}a_{12}\bar{b}_{12}\right|^{2},$$ (19)

where we have used $\lambda_{\rm m}=\lambda_{1},\lambda_{\rm sm}=\lambda_{2}$ and $\lambda_{1}+\lambda_{2}=1$. This clearly implies inequality (14).

[Proof of (15)] We may assume $B$ to take the form $B-\frac{\tr(B)}{n}I$ by using (18), while $A$ retains the form $A-\tr(A\omega)I$ as above. Thus, $a_{22}=-\frac{\lambda_{1}}{\lambda_{2}}a_{11}$ and $b_{22}=-b_{11}$. Noting that $\lambda_{\rm m}=\lambda_{1},\lambda_{\rm M}=\lambda_{2}$, $\lambda_{1}+\lambda_{2}=1$ and $\lambda_{2}\geq\frac{1}{2}\geq\lambda_{1}$, we have

$$\frac{\lambda_{\rm M}+\lambda_{\rm m}}{\lambda_{\rm m}}\|A\|^{2}_{\omega}\|B\|^{2}-\|[A,B]\|^{2}_{\omega}=\frac{2}{\lambda_{2}}|a_{11}|^{2}|b_{11}|^{2}+|a_{21}|^{2}|b_{21}|^{2}+\frac{\lambda_{2}}{\lambda_{1}}|a_{12}|^{2}|b_{12}|^{2}$$

$$\quad+\frac{\lambda_{2}-\lambda_{1}}{\lambda_{1}}|a_{12}|^{2}|b_{21}|^{2}+2(1-2\lambda_{1})|a_{21}|^{2}|b_{11}|^{2}+\frac{\lambda_{2}-\lambda_{1}}{\lambda^{2}_{2}}|a_{11}|^{2}|b_{21}|^{2}+\frac{2\lambda_{2}(1-2\lambda_{1})}{\lambda_{1}}|a_{12}|^{2}|b_{11}|^{2}$$

$$\quad+2{\rm Re}a_{12}\bar{a}_{21}\bar{b}_{12}b_{21}+\frac{4\lambda_{1}}{\lambda_{2}}{\rm Re}a_{11}\bar{a}_{21}\bar{b}_{11}b_{21}+4{\rm Re}a_{11}\bar{a}_{12}\bar{b}_{11}b_{12}$$

$$\geq\frac{2}{\lambda_{2}}|a_{11}|^{2}|b_{11}|^{2}+|a_{21}|^{2}|b_{21}|^{2}+\frac{\lambda_{2}}{\lambda_{1}}|a_{12}|^{2}|b_{12}|^{2}$$

$$\quad+2{\rm Re}a_{12}\bar{a}_{21}\bar{b}_{12}b_{21}+\frac{4\lambda_{1}}{\lambda_{2}}{\rm Re}a_{11}\bar{a}_{21}\bar{b}_{11}b_{21}+4{\rm Re}a_{11}\bar{a}_{12}\bar{b}_{11}b_{12}$$

$$\geq\frac{4\lambda_{1}}{\lambda_{2}}|a_{11}|^{2}|b_{11}|^{2}+\frac{\lambda_{1}}{\lambda_{2}}|a_{21}|^{2}|b_{21}|^{2}+\frac{\lambda_{2}}{\lambda_{1}}|a_{12}|^{2}|b_{12}|^{2}$$

$$\quad+2{\rm Re}a_{12}\bar{a}_{21}\bar{b}_{12}b_{21}+\frac{4\lambda_{1}}{\lambda_{2}}{\rm Re}a_{11}\bar{a}_{21}\bar{b}_{11}b_{21}+4{\rm Re}a_{11}\bar{a}_{12}\bar{b}_{11}b_{12}$$

$$=\left|2\sqrt{\frac{\lambda_{1}}{\lambda_{2}}}a_{11}\bar{b}_{11}+\sqrt{\frac{\lambda_{1}}{\lambda_{2}}}a_{21}\bar{b}_{21}+\sqrt{\frac{\lambda_{2}}{\lambda_{1}}}a_{12}\bar{b}_{12}\right|^{2}$$

which implies inequality (15).

3 Applications to quantum physics

This section is dedicated to exploring how both conjectures (14) and (15) can be effectively applied in the realm of quantum physics, demonstrating their utility and relevance. In Sec. 3.1, we will see that conjecture (14) unveils a new type of uncertainty relation for quantum observables, and in Sec. 3.2, we will discuss how conjecture (15) provides a constraint on relaxation rates in quantum Markovian dynamics.

3.1 Application of Conjecture (14)

Our primary motivation for introducing the $\omega$-norm was to apply BW inequality to the uncertainty relations in quantum physics. In the context of a $n$-level quantum system $\mathop{\mathbb{C}}\nolimits^{n}$, physical quantities (observables) are represented by Hermitian matrices $A,B$, where a quantum state is represented by a density matrix $\rho$, a positive semidefinite matrix with a unit trace. Under a state $\rho$, the expectation value of $A$ is given by $E(A)_{\rho}=\tr\rho A$, and the variance of $A$ is given by

$$V(A)_{\rho}=\tr\rho A^{2}-(\tr\rho A)^{2}=\tr\rho(A-(\tr A\rho)\mathop{\mathbb{I}}\nolimits)^{2}.$$ (20)

The famous Robertson uncertainty relation is

$$V(A)_{\rho}V(B)_{\rho}\geq\frac{1}{4}|\tr\rho[A,B]|^{2}.$$ (21)

Therefore, in the field of quantum physics, it is understood that the non-commutativity of physical quantities results in a trade-off between their uncertainties.

On the other hand, it is straightforward to observe that the BW inequality offers a comparable uncertainty relation for observables when a quantum system is in the maximally mixed state, denoted as $\rho_{\max}=\frac{\mathop{\mathbb{I}}\nolimits}{n}$. Specifically, for this state, the variances can be represented through the Frobenius norm by

$$V(A)_{\rho_{\max}}=\frac{1}{n}\|A-\frac{\tr A}{n}\mathop{\mathbb{I}}\nolimits\|^{2},V(B)_{\rho_{\max}}=\frac{1}{n}\|B-\frac{\tr B}{n}\mathop{\mathbb{I}}\nolimits\|^{2}.$$

Thus, noting that $[A-\frac{\tr A}{n}\mathop{\mathbb{I}}\nolimits,B-\frac{\tr B}{n}\mathop{\mathbb{I}}\nolimits]=[A,B]$, the BW inequality provides the following uncertainty relation:

$$V(A)_{\rho_{\max}}V(B)_{\rho_{\max}}\geq\frac{1}{2n^{2}}\|[A,B]\|^{2}.$$ (22)

Applying the maximally mixed state to the Robertson relation (21), the bound on the right-hand side is given by $\frac{1}{4n^{2}}|\tr[A,B]|^{2}$. Consequently, Eq. (22) yields a similar, yet distinct, uncertainty relation. As an example, consider $A=\sigma_{x}$ and $B=\sigma_{y}$, where $\sigma_{x}=\begin{pmatrix}0&1\\ 1&0\end{pmatrix}$, $\sigma_{y}=\begin{pmatrix}0&-i\\ i&0\end{pmatrix}$, and $\sigma_{z}=\begin{pmatrix}1&0\\ 0&-1\end{pmatrix}$ are the Pauli matrices for a two-level quantum system. Since $[\sigma_{x},\sigma_{y}]=2i\sigma_{z}$ and $\tr\sigma_{z}=0$ the Robertson bound yields a trivial bound $0$, while the bound in (22) is $1$.

Unfortunately, the relation (22) is only valid for the maximally mixed state. To extend this to an arbitrary state, the $\omega$-norm becomes essential: By defining $\tilde{A}:=A-(\tr\rho A)\mathop{\mathbb{I}}\nolimits$ and $\tilde{B}:=B-(\tr\rho B)\mathop{\mathbb{I}}\nolimits$, the variance of $A$ and $B$ under a state $\rho$ can be represented by

$$V(A)_{\rho}=\|\tilde{A}\|^{2}_{\rho},V(B)_{\rho}=\|\tilde{B}\|^{2}_{\rho}.$$ (23)

Thus, type (i) bound (Eq. (i)) — as well as type (iv) — introduces another uncertainty principle that is valid for any quantum state $\rho$. Notably, conjecture (14) leads to a novel form of the uncertainty relation: For any faithful quantum state $\rho$ (i.e., a positive definite matrix with unit trace) and for any observables $A$ and $B$,

$$V(A)_{\rho}V(B)_{\rho}\geq\frac{\lambda_{\rm m}\lambda_{\rm sm}}{\lambda_{\rm m}+\lambda_{\rm sm}}\|[A,B]\|^{2}_{\rho},$$ (24)

where $\lambda_{\rm m}$ and $\lambda_{\rm sm}$ denote the smallest and the second smallest eigenvalues of $\rho$, respectively. For comparison, if we use a bound $c_{1}$ (which is already proven) given in (11), we have the looser uncertainty relation:

$$V(A)_{\rho}V(B)_{\rho}\geq\frac{\lambda_{\rm m}^{2}}{2\lambda_{\rm M}}\|[A,B]\|^{2}_{\rho}.$$ (25)

As illustrations, let’s consider $A=\sigma_{x}$ and $B=\sigma_{y}$ again and a quantum state which is a probabilistic mixture of the maximally mixed state and the eigenstate $\ket{0}$ of $\sigma_{z}$ corresponding to the eigenvalue $1$, i.e., $\rho(p)=p\frac{\mathop{\mathbb{I}}\nolimits}{2}+(1-p)\outerproduct{0}{0}\ (p\in(0,1])$. Simple computations shows that the bounds in (21), (24) and (25) are given respectively by $(1-p)^{2},p(2-p)$, and $\frac{p^{2}}{2-p}$. Therefore, the bound (24) outperforms the Robertson bound (21) if $p>(2-\sqrt{2})/2\simeq 0.293$ (See Fig. 2). Interestingly, even the looser bound (25) surpasses the Robertson bound if $p\gtrsim 0.547$. In [10], we conduct a comprehensive comparison of these new uncertainty relations with the standard ones, including not only the Robertson relation but also the Schrödinger relation. It is demonstrated that the bounds in (24) and (25) unveil entirely new types of trade-offs in quantum uncertainty that were previously undetected by the standard relations, particularly in cases where the state exhibits a higher degree of mixedness.

Figure 2: The Robertson bound (21) and our bound (24) between observables $\sigma_{x}$ and $\sigma_{y}$ under a qubit state $\rho(p)=p\frac{\mathop{\mathbb{I}}\nolimits}{2}+(1-p)\outerproduct{0}{0}\ (p\in(0,1])$ are plotted with dashed and solid lines, respectively. For comparison, the looser bound, using $c_{1}$ in (11) is also plotted with dot-dashed line.

3.2 Application of Conjecture (15)

Although conjecture (15) employs an asymmetric combination of the $\omega$-norm and the Frobenius norm, making it appear artificial and pedantic, it finds direct application in the field of open quantum dynamics.

In the field of open quantum system, it is widely recognized that general quantum Markovian dynamics is described by a completely positive (CP) dynamical semigroup:

$$\frac{d\rho_{t}}{dt}=\mathop{\mathcal{L}}\nolimits\rho_{t},$$

where $\rho_{t}$ is a quantum state at time $t$ and $\mathop{\mathcal{L}}\nolimits$ is the generator of time evolution. One of the seminal works in this area is the representation theorem of the generator: $\mathop{\mathcal{L}}\nolimits$ qualifies as a generator for a CP dynamical semigroup if and only if it has the form

$$\mathop{\mathcal{L}}\nolimits(\rho)=-i[H,\rho]+\frac{1}{2}\sum_{k}\gamma_{k}(2L_{k}\rho L_{k}^{\dagger}-L_{k}^{\dagger}L_{k}\rho-\rho L_{k}^{\dagger}L_{k}),$$ (26)

where $\gamma_{k}\geq 0$ and jump operators $L_{k}$ are normalized as $\|L_{k}\|=1$. Since this was independently discovered by Gorini, Kossakowski, Sudarshan [11] and by Lindblad [12], the generator is now referred to as the GKLS (or GKSL) generator (See also [13]).

Letting $\mathop{\mathcal{L}}\nolimits^{\ddagger}$ be the dual of $\mathop{\mathcal{L}}\nolimits$ via $\tr(X\mathop{\mathcal{L}}\nolimits(Y))=\tr(\mathop{\mathcal{L}}\nolimits^{\ddagger}(X)Y)$ and $\ell_{\alpha}$ and $Y_{\alpha}$ be an eigenvalue and an eigenvector of $\mathop{\mathcal{L}}\nolimits^{\ddagger}$: $\mathop{\mathcal{L}}\nolimits^{\ddagger}(Y_{\alpha})=\ell_{\alpha}Y_{\alpha}$. Note that a relaxation rate $\Gamma_{\alpha}$ is given by the real part of the eigenvalue $\ell_{\alpha}$. Let $\omega$ be a faithful stationary state, i.e., $\mathop{\mathcal{L}}\nolimits(\omega)=0$. In [14], a relaxation rate $\Gamma_{\alpha}$ for any GKLS generator has been characterized by

$$\Gamma_{\alpha}=\frac{1}{2\|Y_{\alpha}\|^{2}_{\omega}}\sum_{k}\gamma_{k}\|[L_{k},Y_{\alpha}]\|^{2}_{\omega}.$$ (27)

Now, by applying conjecture (15) and recalling that $\|L_{k}\|=1$, we have

$$\Gamma_{\alpha}\leq\frac{1}{2}\Bigl{(}1+\frac{\lambda_{\rm M}}{\lambda_{\rm m}}\Bigr{)}\sum_{k}\gamma_{k},$$

where $\lambda_{\rm M}$ and $\lambda_{\rm m}$ are the maximal and the minimal eigenvalues of the stationary state $\omega$. On the other hand, as we have $\sum_{k}\gamma_{k}=\frac{1}{n}\sum_{\beta=1}^{n^{2}-1}\Gamma_{\beta}$ [15, 16], we obtain the following constraints on relaxation rates for GKLS master equation:

$$\Gamma_{\alpha}\leq\frac{1}{2n}\Bigl{(}1+\frac{\lambda_{\rm M}}{\lambda_{\rm m}}\Bigr{)}\sum_{\beta=1}^{n^{2}-1}\Gamma_{\beta}.$$ (28)

Note that relaxation rates are observable in experiments, the above constraint can be directly tested in experimental setups.

4 Concluding remarks

In this paper, we investigated generalizations of the BW inequality, utilizing combinations of the $\omega$-norm and the Frobenius norm. There are six types of generalization, types (i)-(vi), with type (vi) being the BW inequality itself. We have established the tight bounds for types (iii) and (v), and offered conjectures for the tight bounds of types (i) and (ii), with type (iv) emerging as a corollary of type (i). Both conjectures for types (i) and (ii) are backed by numerical evidence and have been validated for matrices of size $n=2$ and certain special cases for general $n$. Additionally, we showed that both types (i) and (ii) have direct applications to quantum physics. The type (i) bound introduces a novel uncertainty relation that reveals a previously undetected trade-off in the uncertainties between non-commuting observables. The type (ii) bound imposes a significant constraint on relaxation rates within general quantum Markovian dynamics, offering a directly testable prediction in experimental settings.

We hope that our results open a new direction of generalizations of BW inequality, enriching both the mathematical and physical perspectives on the subject.

Appendix A Proofs of conjectures for special cases

In this appendix, we give proofs of (14) and (15) in some restricted cases.

[The case where $A$ is normal and commutes with $\omega$]

Here we give proofs of (14) and (15) in the case where $A$ is normal and commutes with $\omega$. In this case, we can assume that both $A$ and $\omega$ are diagonal: Let $A={\rm diag}[a_{1},\ldots,a_{n}]\ (a_{i}\in{\mathbb{C}})$ and $\omega={\rm diag}[\lambda_{1},\ldots,\lambda_{n}]$ $(\lambda_{1}\leq\lambda_{2}\leq\cdots\leq\lambda_{n})$ and $B=[b_{ij}]\in M_{n}({\mathbb{C}})$.

Now, a direct computation yields:

$$\frac{\lambda_{\rm m}+\lambda_{\rm sm}}{\lambda_{\rm m}\lambda_{\rm sm}}\|A\|^{2}_{\omega}\|B\|^{2}_{\omega}-\|[A,B]\|^{2}_{\omega}$$

$$\quad=\sum_{i,j}|b_{ij}|^{2}\lambda_{j}\left(\frac{\lambda_{1}+\lambda_{2}}{\lambda_{1}\lambda_{2}}\sum_{k}|a_{k}|^{2}\lambda_{k}-|a_{i}-a_{j}|^{2}\right)$$

$$\quad\geq\sum_{i,j}|b_{ij}|^{2}\lambda_{j}\left(\left(\frac{\lambda_{i}}{\lambda_{1}}+\frac{\lambda_{i}}{\lambda_{2}}-1\right)|a_{i}|^{2}+2{\rm Re}(a_{i}\bar{a}_{j})+\left(\frac{\lambda_{j}}{\lambda_{1}}+\frac{\lambda_{j}}{\lambda_{2}}-1\right)|a_{j}|^{2}\right).$$ (29)

Note that $T(i):=\frac{\lambda_{i}}{\lambda_{1}}+\frac{\lambda_{i}}{\lambda_{2}}-1$ is larger than or equal to $1$ when $i\neq 1$, and satisfies $T(i)T(1)\geq 1$ for any $2\leq i\leq n$. This clearly shows that (29) is non-negative, confirming (14).

Similarly, (15) can also be demonstrated as

$$\frac{\lambda_{\rm M}+\lambda_{\rm m}}{\lambda_{\rm m}}\|A\|^{2}_{\omega}\|B\|^{2}-\|[A,B]\|^{2}_{\omega}$$

$$\quad=\sum_{i,j}|b_{ij}|^{2}\left(\frac{\lambda_{n}+\lambda_{\rm 1}}{\lambda_{\rm 1}}\sum_{k}|a_{k}|^{2}\lambda_{k}-|a_{i}-a_{j}|^{2}\lambda_{j}\right)$$

$$\quad\geq\sum_{i,j}|b_{ij}|^{2}\left(\left(\lambda_{i}+\frac{\lambda_{n}}{\lambda_{\rm 1}}\lambda_{i}-\lambda_{j}\right)|a_{i}|^{2}+2{\rm Re}(a_{i}\bar{a}_{j}\lambda_{j})+\frac{\lambda_{n}}{\lambda_{\rm 1}}\lambda_{j}|a_{j}|^{2}\right)$$

$$\quad\geq\sum_{i,j}|b_{ij}|^{2}\left(\lambda_{i}|a_{i}|^{2}+2{\rm Re}(a_{i}\bar{a}_{j}\lambda_{j})+\frac{\lambda_{j}}{\lambda_{i}}\lambda_{j}|a_{j}|^{2}\right)$$

$$\quad=\sum_{i,j}|b_{ij}|^{2}\left|\sqrt{\lambda_{i}}a_{i}+\frac{\lambda_{j}}{\sqrt{\lambda_{i}}}a_{j}\right|^{2}\geq 0.$$ (30)

[The case where $B$ commutes with $\omega$]

When matrix $B$ commutes with $\omega$, conjecture (15) can be easy to show. Moreover, the tight bound in this case coincides with BW bound, i.e., $\|[A,B]\|_{\omega}\leq\sqrt{2}\|A\|_{\omega}\|B\|$. To show this, we use the relation:

$$\|A\|_{\omega}=\|A\sqrt{\omega}\|,$$ (31)

which follows from the cyclic property of trace as $\|A\|^{2}_{\omega}=\tr A^{\ast}A\omega=\tr\sqrt{\omega}A^{\ast}A\sqrt{\omega}=\tr(A\sqrt{\omega})^{\ast}A\sqrt{\omega}=\|A\sqrt{\omega}\|^{2}$. Given that $B$ commutes with $\omega$, and therefore with $\sqrt{\omega}$, the relation $[B,A]\sqrt{\omega}=[B,A\sqrt{\omega}]$ is satisfied. By using (31) and BW inequality, we have

$$\|[A,B]\|^{2}_{\omega}=\|[B,A]\|^{2}_{\omega}=\|[B,A]\sqrt{\omega}\|^{2}=\|[B,A\sqrt{\omega}]\|^{2}\leq\sqrt{2}\|B\|^{2}\|A\sqrt{\omega}\|^{2}=\sqrt{2}\|A\|^{2}_{\omega}\|B\|^{2}.$$

[The case where $B=\outerproduct{\lambda_{1}}{\lambda_{n}}$]

Our numerical experiments indicate that the optimal matrix $B$ for achieving the bound $\tilde{c_{2}}(\omega)$, as specified in (4), invariably takes the form $B=\outerproduct{\lambda_{1}}{\lambda_{n}}$ (after the normalization). Once we assume this fact, we can prove conjecture (15) as follows. Denoting $a_{ij}:=\langle\lambda_{i}|A\lambda_{j}\rangle$, a direct computation yields $\norm{[A,B]}_{\omega}^{2}=\lambda_{n}\sum_{j}|a_{j1}|^{2}-2\lambda_{n}\real(a_{nn}\overline{a_{11}})+\sum_{j}\lambda_{j}|a_{nj}|^{2}$ , $\norm{A}_{\omega}^{2}=\sum_{j}\lambda_{j}\sum_{k}|a_{kj}|^{2}$, $\|B\|^{2}=1$. Consequently, it follows that

$$(\lambda_{\rm m}+\lambda_{\rm M})\norm{A}^{2}_{\omega}\norm{B}^{2}-\lambda_{\rm m}\norm{[A,B]}_{\omega}^{2}$$

$$= \lambda_{1}\sum_{\begin{subarray}{c}j\\ k\neq n\end{subarray}}\lambda_{j}|a_{kj}|^{2}+\lambda_{n}\sum_{\begin{subarray}{c}j\neq 1\\ k\end{subarray}}\lambda_{j}|a_{kj}|^{2}+2\lambda_{1}\lambda_{n}\real(a_{nn}\overline{a_{11}})$$ (32)

$$= \Bigl{|}\lambda_{1}a_{11}+\lambda_{n}a_{nn}\Bigr{|}^{2}+\lambda_{1}^{2}\sum_{\begin{subarray}{c}k\neq 1,n\end{subarray}}|a_{k1}|^{2}+\lambda_{1}\sum_{\begin{subarray}{c}j\neq 1\\ k\neq n\end{subarray}}\lambda_{j}|a_{kj}|^{2}+\lambda^{2}_{n}\sum_{\begin{subarray}{c}k\neq n\end{subarray}}|a_{kn}|^{2}+\lambda_{n}\sum_{\begin{subarray}{c}j\neq 1,n\\ k\end{subarray}}\lambda_{j}|a_{kj}|^{2}\geq 0$$ (33)