Skip to article frontmatterSkip to article content
Site not loading correctly?

This may be due to an incorrect BASE_URL configuration. See the MyST Documentation for reference.

The-Quantum-Foundations-Project

There is no general consensus as to what its fundamental principles are, how it should be taught, or what it really “means.” Every competent physicist can “do” quantum mechanics, but the stories we tell ourselves about what we are doing are as various as the tales of Scheherazade, and almost as implausible (Griffiths (2017)).

Analogy with Classical Mechanics

Newton’s Second Law

A law of physics can be understood as a rule that describes—usually by means of a mathematical formulation—how certain physical systems behave(Siddiqui (2014)). Newton’s second law, for example, can be expressed by the equation

F(x,t)=md2dt2x(t)F\left(x,t\right)=m\frac{d^{2}}{dt^{2}}x\left(t\right)

where x(t)x\left(t\right) represents the position of a particle of mass mm relative to the origin at the instant tt. When we solve this differential equation under certain initial conditions, we obtain as a solution a function x(t)x\left(t\right). Thus, once the forces acting on the particle and its initial conditions are specified, Newton’s second law determines its temporal evolution: solving the equation provides us with the equation of motion of the particle.

For example, for a free particle(F=0F=0), we can obtain its equation of motion from Newton’s second law:

0=md2xdt20=m\frac{d^{2}x}{dt^{2}}

Defining the velocity as v(t)≡dxdtv\left(t\right)\equiv\frac{dx}{dt}, we have:

0=d2xdt20=dvdt0⋅∫t0tdt′=∫t0tdvdt′dt′0=v(t)−v(t0)\begin{align} 0 & =\frac{d^{2}x}{dt^{2}}\\ 0 & =\frac{dv}{dt}\\ 0\cdot\int^{t}_{t_{0}}dt' & =\int^{t}_{t_{0}}\frac{dv}{dt'}dt'\\ 0 & =v\left(t\right)-v\left(t_{0}\right) \end{align}

Therefore, v(t)=v(t0)v\left(t\right)=v\left(t_{0}\right), that is, the particle’s velocity is constant. Denoting this constant velocity by vv, we have:

v(t)=v(t0)dxdt=v∫t0tdxdt′dt′=v∫t0tdt′x(t)−x(t0)=v(t−t0)\begin{align} v\left(t\right) & =v\left(t_{0}\right)\\ \frac{dx}{dt} & =v\\ \int^{t}_{t_{0}}\frac{dx}{dt'}dt' & =v\int^{t}_{t_{0}}dt'\\ x\left(t\right)-x\left(t_{0}\right) & =v\left(t-t_{0}\right) \end{align}

If x(t0)=x0x\left(t_{0}\right)=x_{0} and t0=0t_{0}=0, we obtain:

x(t)=x0+vtx\left(t\right)=x_{0}+vt

In this case, x(t)x\left(t\right) is easily interpretable as the position of the particle at the instant tt.

Schrödinger Equation

The quantum mechanical analogue of Newton’s second law is the Schrödinger equation. In one dimension and in the non-relativistic case, the Schrödinger equation is:

iℏ∂∂tΨ(x,t)=−ℏ22m∂2∂x2Ψ(x,t)+V(x,t)Ψ(x,t)i\hbar\frac{\partial}{\partial t}\Psi\left(x,t\right)=-\frac{\hbar^{2}}{2m}\frac{\partial^{2}}{\partial x^{2}}\Psi\left(x,t\right)+V\left(x,t\right)\Psi\left(x,t\right)

The corresponding equation for a free particle (V(x,t)=0V\left(x,t\right)=0) is:

iℏ∂Ψ∂t=−ℏ22m∂2Ψ∂x2\begin{align} i\hbar\frac{\partial\Psi}{\partial t} & =-\frac{\hbar^{2}}{2m}\frac{\partial^{2}\Psi}{\partial x^{2}} \end{align}

Using the method of separation of variables, that is, assuming that the solution can be separated into a spatial part Φ(x)\Phi\left(x\right) and a temporal part Π(t)\Pi\left(t\right), that is, Ψ(x,t)=Φ(x)Π(t)\Psi\left(x,t\right)=\Phi\left(x\right)\Pi\left(t\right). Then:

iℏ∂∂t(Φ(x)Π(t))=−ℏ22m∂2∂x2(Φ(x)Π(t))iℏ1Π∂Π∂t=−ℏ22m1Φ∂2Φ∂x2\begin{align} i\hbar\frac{\partial}{\partial t}\left(\Phi\left(x\right)\Pi\left(t\right)\right) & =-\frac{\hbar^{2}}{2m}\frac{\partial^{2}}{\partial x^{2}}\left(\Phi\left(x\right)\Pi\left(t\right)\right)\\ i\hbar\frac{1}{\Pi}\frac{\partial\Pi}{\partial t} & =-\frac{\hbar^{2}}{2m}\frac{1}{\Phi}\frac{\partial^{2}\Phi}{\partial x^{2}} \end{align}

Since one side depends only on the spatial dimension (xx) and the other only on the temporal dimension (tt), the only way this can be true is if both sides are equal to a constant: EE. Solving then for the temporal part:

iℏΠ∂Π∂t=EdΠdt=−iEℏΠ\begin{align} \frac{i\hbar}{\Pi}\frac{\partial\Pi}{\partial t} & =E\\ \frac{d\Pi}{dt} & =-\frac{iE}{\hbar}\Pi \end{align}

Since the derivative of the function Π(t)\Pi\left(t\right) is equal to the function itself multiplied by a constant, we have an exponential solution Π(t)=Aexp⁡(−iEℏt)\Pi\left(t\right)=A\exp\left(-\frac{iE}{\hbar}t\right). For the spatial part:

1Φ∂2Φ∂x2=−2mEℏ2∂2Φ∂x2=−2mEℏ2Φ\begin{align} \frac{1}{\Phi}\frac{\partial^{2}\Phi}{\partial x^{2}} & =-\frac{2mE}{\hbar^{2}}\\ \frac{\partial^{2}\Phi}{\partial x^{2}} & =-\frac{2mE}{\hbar^{2}}\Phi \end{align}

Since this is a second-order linear homogeneous differential equation with constant coefficients, we can look for solutions of the form Φ(x)=e±ikx\Phi\left(x\right)=e^{\pm ikx}. In this case, eikxe^{ikx} and e−ikxe^{-ikx} are two linearly independent solutions to d2Φdx2=−k2Φ\frac{d^{2}\Phi}{dx^{2}}=-k^{2}\Phi. By the linearity of the equation, its general solution is a linear combination of these solutions:

Ψ(x,t)=Π(t)Φ(x)=Aexp⁡(−iEℏt)(Bexp⁡(−i2mEℏ2x)+Cexp⁡(i2mEℏ2x))\Psi\left(x,t\right)=\Pi\left(t\right)\Phi\left(x\right)=A\exp\left(-\frac{iE}{\hbar}t\right)\left(B\exp\left(-i\sqrt{\frac{2mE}{\hbar^{2}}}x\right)+C\exp\left(i\sqrt{\frac{2mE}{\hbar^{2}}}x\right)\right)

Or, writing k2=2mE/ℏ2k^{2}=2mE/\hbar^{2} and ω=Eℏ\omega=\frac{E}{\hbar}:

Ψ(x,t)=K1e−i(kx+ωt)+K2ei(kx−ωt)\begin{align} \Psi\left(x,t\right) & =K_{1}e^{-i\left(kx+\omega t\right)}+K_{2}e^{i\left(kx-\omega t\right)} \end{align}

where K1K_{1} and K2K_{2} are constants determined by the conditions imposed on the solution. Here, some questions already arise. When we solved Newton’s second law, we obtained x(t)x\left(t\right), which described the position of the particle. When we solve the Schrödinger equation, we obtain Ψ(x,t)\Psi\left(x,t\right), but what is Ψ(x,t)\Psi\left(x,t\right)? What does Ψ(x,t)\Psi\left(x,t\right) describe?

A somewhat more in-depth discussion of this can be found here.

Postulates of quantum mechanics

Notation

To move forward, we need to know the postulates[1] of quantum mechanics. The first point to note is that they can be formulated in different ways from one book to another. We will use Professor Jonas Maziero’s material (Maziero (2026)). We need to begin, then, with some comments on Dirac notation.

In this notation, we write a vector as ∣v⟩\left|v\right\rangle, called a ket. Taking the conjugate transpose of this vector (or adjoint), we obtain ⟨v∣=∣v⟩†=(∣v⟩T)∗\left\langle v\right|=\left|v\right\rangle ^{\dagger}=\left(\left|v\right\rangle ^{T}\right)^{*}, which is called a bra. If we represent the vector ∣v⟩\left|v\right\rangle in an orthonormal basis {∣bj⟩}\left\{ \left|b_{j}\right\rangle \right\}, we can write it as a linear combination of the basis vectors: ∣v⟩=∑i⟨bi∣v⟩∣bi⟩\left|v\right\rangle =\sum_{i}\left\langle b_{i}|v\right\rangle \left|b_{i}\right\rangle. In matrix representation, we have

∣v⟩=[⟨b1∣v⟩⋮⟨bn∣v⟩]\left|v\right\rangle =\left[\begin{array}{c} \left\langle b_{1}|v\right\rangle \\ \vdots\\ \left\langle b_{n}|v\right\rangle \end{array}\right]

The product of a bra and a ket is called the inner product and results in a scalar:

⟨w∣v⟩=∣w⟩†∣v⟩=[⟨b1∣w⟩∗…⟨bn∣w⟩∗][⟨b1∣v⟩⋮⟨bn∣v⟩]=∑j=1n⟨bj∣w⟩∗⟨bj∣v⟩\left\langle w|v\right\rangle =\left|w\right\rangle ^{\dagger}\left|v\right\rangle =\left[\begin{array}{ccc} \left\langle b_{1}|w\right\rangle ^{*} & \dots & \left\langle b_{n}|w\right\rangle ^{*}\end{array}\right]\left[\begin{array}{c} \left\langle b_{1}|v\right\rangle \\ \vdots\\ \left\langle b_{n}|v\right\rangle \end{array}\right]=\sum^{n}_{j=1}\left\langle b_{j}|w\right\rangle ^{*}\left\langle b_{j}|v\right\rangle

On the other hand, the product of a ket and a bra is called the outer product and results in an operator, whose matrix representation is:

∣v⟩⟨w∣=∣v⟩∣w⟩†=[⟨b1∣v⟩⋮⟨bn∣v⟩][⟨b1∣w⟩∗…⟨bn∣w⟩∗]=[⟨b1∣v⟩⟨b1∣w⟩∗…⟨b1∣v⟩⟨bn∣w⟩∗⋮⋮⋮⟨bn∣v⟩⟨b1∣w⟩∗…⟨bn∣v⟩⟨bn∣w⟩∗]\left|v\right\rangle \left\langle w\right|=\left|v\right\rangle \left|w\right\rangle ^{\dagger}=\left[\begin{array}{c} \left\langle b_{1}|v\right\rangle \\ \vdots\\ \left\langle b_{n}|v\right\rangle \end{array}\right]\left[\begin{array}{ccc} \left\langle b_{1}|w\right\rangle ^{*} & \dots & \left\langle b_{n}|w\right\rangle ^{*}\end{array}\right]=\left[\begin{array}{ccc} \left\langle b_{1}|v\right\rangle \left\langle b_{1}|w\right\rangle ^{*} & \dots & \left\langle b_{1}|v\right\rangle \left\langle b_{n}|w\right\rangle ^{*}\\ \vdots & \vdots & \vdots\\ \left\langle b_{n}|v\right\rangle \left\langle b_{1}|w\right\rangle ^{*} & \dots & \left\langle b_{n}|v\right\rangle \left\langle b_{n}|w\right\rangle ^{*} \end{array}\right]

With this in mind, we can now move on to the postulates of quantum mechanics.

State Postulate

In every physical theory, we seek to describe the possible states of the system of interest, the changes in these states (temporal evolution), and the observable quantities and their measurements. For example, in Newtonian Mechanics (NM), the state of a system at a certain instant of time is specified by the values of position and velocity relative to a certain inertial reference frame, and its temporal evolution is given by Newton’s second law. The observables are essentially the positions, from which we obtain the velocities using the time parameter. It is implicit in NM that the values of the observables are determined even before they are measured, and that all these observables can have well-defined values at a given instant of time. Next, we will see the postulates of QM for describing states, changes in these states, observables, and their measurements. We will see that the measurement part differs greatly from that of NM, mainly because of the existence of observables that cannot have well-defined values at the same time (Maziero (2026)).

The state of a quantum system is described by a normalized vector in a Hilbert space H\mathcal{H}. For simplicity, we will consider H\mathcal{H} as a vector space equipped with an inner product function ⟨⋅∣⋅⟩:H→C\left\langle \cdot|\cdot\right\rangle :\mathcal{H}\rightarrow\mathbb{C}.

Furthermore, a Hilbert space is a complete complex vector space. To understand what it means to be complete, we first need to understand the concept of a Cauchy sequence. A Cauchy sequence is a sequence whose elements become arbitrarily close to one another as we advance through the sequence. A complete metric space[2] is one in which every Cauchy sequence converges to a limit that also belongs to the space itself.

By system, we refer to a degree of freedom[3]. For example, the spin of an electron is a system, another system would be its position, and so on. When we consider two or more degrees of freedom as a single entity, we call it a composite system, and if we divide the composite system into subsystems, its state is described by a normalized vector in a Hilbert space obtained by taking the tensor product of the individual spaces:

H=Ha⊗Hb⊗…\mathcal{H}=\mathcal{H}_{a}\otimes\mathcal{H}_{b}\otimes\dots

Measurement Postulate

Observable (measurable) quantities are represented by linear Hermitian operators, O=O†O=O^{\dagger}, defined on H\mathcal{H}. When we measure an observable, we will obtain one of its possible values, which are given by the eigenvalues of the Hermitian operator.

Unlike classical mechanics (CM), quantum mechanics (QM), in general, does not allow us to predict with certainty the outcome of a measurement. What QM provides us with are the probabilities of obtaining the different outcomes, calculated by means of the Born rule:

P(o∣ψ)=∣⟨ψ∣o⟩∣2=∣⟨o∣ψ⟩∣2P\left(o|\psi\right)=\left|\left\langle \psi|o\right\rangle \right|^{2}=\left|\left\langle o|\psi\right\rangle \right|^{2}

This is the probability of obtaining the value oo when measuring the observable OO of a system that is in the state ∣ψ⟩\left|\psi\right\rangle immediately before the measurement, and ∣o⟩\left|o\right\rangle is the eigenvector of OO corresponding to the eigenvalue oo. O=O†O=O^{\dagger} are defined on H\mathcal{H}. When we measure an observable, we obtain one of its possible values, which are represented by the eigenvalues of the Hermitian operator.

Before moving forward, we can stop to reflect for a moment. When we were working with classical mechanics and Newton’s law, one of the observables of a particle was its position, and Newton’s law directly described how this quantity varied in time.

However, when we approach quantum mechanics, the Schrödinger equation describes how the wave function varies in time (and space). With Dirac notation, as we will see when studying the dynamics postulate, we have a mathematical description of how the state vector of the system varies in time. However, the state vector is not itself an observable.

Even temporarily ignoring the probabilistic aspect, we use the state vector to calculate quantities associated with observables. This already constitutes a fundamental difference between the classical and quantum descriptions of systems. But if we calculate scalar quantities associated with observables from the state vector (or the wave function) of the system, what is the physical meaning of the state vector itself? Does it even have a direct physical meaning?

This is one of the points on which different interpretations of quantum mechanics diverge. There is still no consensus on which interpretation provides the correct description of physical reality, although some interpretations may be ruled out or considered problematic in light of certain theoretical arguments and experimental results.

Likewise, different interpretations diverge on the question of whether all relevant physical information about a system is necessarily contained in its state vector. We will explore different interpretations later on. For now, it is enough to emphasize that the physical meaning of the state vector and the status of the information it contains are profound interpretational questions on which there is no consensus.

Repeatability of measurements

However, there is one more detail that, although it is not explicitly contained in the postulates we have presented so far, is significantly different when we are in the quantum domain: performing a measurement affects the system.

When we observe the position of a soccer ball that moves according to classical mechanics, we can, in principle, make this observation without significantly altering the state of the ball. We can therefore consider that the observation was performed approximately “externally” to the system itself, leaving it practically unchanged.

In the quantum domain, however, the situation is different. When we perform a measurement, the system interacts with the measurement apparatus and, in general, its state is significantly altered. One consequence of this is that, if we measure the same observable twice immediately in succession, the result of the second measurement will necessarily be equal to that obtained in the first.

Indeed, suppose that the first measurement of the observable (OO) produced the result (oo). After this measurement, the state of the system becomes a state associated with the result (oo). In the case where the eigenvalue (oo) is non-degenerate[4], we can write

P(o∣ψ)=1,→∣ψ⟩=eiϕ∣o⟩P\left(o|\psi\right)=1,\quad\rightarrow\left|\psi\right\rangle =e^{i\phi}\left|o\right\rangle

A second immediately subsequent measurement of the same observable will therefore once again produce the result (oo).

Note that, theoretically, before the measurement is performed, quantum mechanics provides us only with the probabilities associated with the possible outcomes. In the laboratory, however, an individual realization of the measurement necessarily produces a definite outcome. This leads us to another fundamental question: the probabilistic nature of quantum mechanics. Its probabilistic predictions concern the frequencies and distributions of outcomes obtained from ensembles of systems prepared in equivalent ways and subjected to the same experimental procedure, and they do not, in general, determine which specific outcome will be obtained in an individual realization.

Before the measurement is performed, we have only the probabilities associated with the different values we may obtain. In this situation, in what state is the system? When dealing with classical mechanics, deterministic dynamics allows us, in principle, to determine the state of a system at an earlier instant t′<tt'<t by knowing its state at an instant tt. Can we extend this line of reasoning to quantum mechanics? Is the probabilistic character of the theory a fundamental feature of nature, or merely an expression of our ignorance about the state of the system? Again, we have no conclusive answers to this question, and different interpretations of quantum mechanics offer different answers.

And finally, still related to the probabilistic nature of measurements, misguided interpretations often support the idea that quantum mechanics is a scientific formulation that naturally incorporates the possibility of free will. However, the probabilistic character of the theory, by itself, provides no support for this conclusion and does not make quantum mechanics compatible with free will. We can discuss this misconception in greater detail when we address the dynamics postulate.

Quantum incompatibility

Let AA and BB be two Hermitian operators defined on a finite-dimensional Hilbert space of dimension dd. By the spectral theorem, we can write them in their respective spectral decompositions as

A=∑j=1daj∣aj⟩⟨aj∣B=∑j=1dbj∣bj⟩⟨bj∣\begin{align} A= & \sum^{d}_{j=1}a_{j}\left|a_{j}\right\rangle \left\langle a_{j}\right| & B= & \sum^{d}_{j=1}b_{j}\left|b_{j}\right\rangle \left\langle b_{j}\right| \end{align}

The numbers aja_{j} and bjb_{j} are the respective eigenvalues, and the sets{∣aj⟩}j=1d\left\{ \left|a_{j}\right\rangle \right\} ^{d}_{j=1} and {∣bj⟩}j=1d\left\{ \left|b_{j}\right\rangle \right\} ^{d}_{j=1} form orthonormal bases of eigenvectors of AA and BB, respectively, that is, they diagonalize these operators, so that:

A∣aj⟩=aj∣aj⟩B∣bj⟩=bj∣bj⟩\begin{align} A\left|a_{j}\right\rangle & =a_{j}\left|a_{j}\right\rangle & B\left|b_{j}\right\rangle = & b_{j}\left|b_{j}\right\rangle \end{align}

Assuming that the two operators share the same basis ∣aj⟩=∣bj⟩≔∣j⟩\left|a_{j}\right\rangle =\left|b_{j}\right\rangle \coloneqq\left|j\right\rangle, then, with the commutator defined as [A,B]=AB−BA\left[A,B\right]=AB-BA, we have:

[A,B]=(∑j=1daj∣j⟩⟨j∣)(∑k=1dbk∣k⟩⟨k∣)−(∑k=1dbk∣k⟩⟨k∣)(∑j=1daj∣j⟩⟨j∣)=∑j,k=1dajbk∣j⟩⟨j∣k⟩⟨k∣−∑j,k=1dbkaj∣k⟩⟨k∣j⟩⟨j∣=∑j,k=1dajbk∣j⟩⟨k∣δjk−∑j,k=1dbkaj∣k⟩⟨j∣δjk=∑j=1d(ajbj∣j⟩⟨j∣−bjaj∣j⟩⟨j∣)=0\begin{align} \left[A,B\right] & =\left(\sum^{d}_{j=1}a_{j}\left|j\right\rangle \left\langle j\right|\right)\left(\sum^{d}_{k=1}b_{k}\left|k\right\rangle \left\langle k\right|\right)-\left(\sum^{d}_{k=1}b_{k}\left|k\right\rangle \left\langle k\right|\right)\left(\sum^{d}_{j=1}a_{j}\left|j\right\rangle \left\langle j\right|\right)\\ & =\sum^{d}_{j,k=1}a_{j}b_{k}\left|j\right\rangle \left\langle j|k\right\rangle \left\langle k\right|-\sum^{d}_{j,k=1}b_{k}a_{j}\left|k\right\rangle \left\langle k|j\right\rangle \left\langle j\right|\\ & =\sum^{d}_{j,k=1}a_{j}b_{k}\left|j\right\rangle \left\langle k\right|\delta_{jk}-\sum^{d}_{j,k=1}b_{k}a_{j}\left|k\right\rangle \left\langle j\right|\delta_{jk}\\ & =\sum^{d}_{j=1}\left(a_{j}b_{j}\left|j\right\rangle \left\langle j\right|-b_{j}a_{j}\left|j\right\rangle \left\langle j\right|\right)\\ & =0 \end{align}

Then:

0=[A,B]=AB−BA⟹AB=BA0=\left[A,B\right]=AB-BA\Longrightarrow AB=BA

Thus, we have that if the operators share the same set of eigenvectors, they commute. And to prove that if they commute, they share the same basis of eigenvectors:

A(B∣aj⟩)=(AB)∣aj⟩=B(A∣aj⟩)=B(aj∣aj⟩)=ajB∣aj⟩\begin{align} A\left(B\left|a_{j}\right\rangle \right) & =\left(AB\right)\left|a_{j}\right\rangle \\ & =B\left(A\left|a_{j}\right\rangle \right)\\ & =B\left(a_{j}\left|a_{j}\right\rangle \right)\\ & =a_{j}B\left|a_{j}\right\rangle \end{align}

That is, in this case, we must have that B∣aj⟩B\left|a_{j}\right\rangle and ∣aj⟩\left|a_{j}\right\rangle belong to the eigenspace of AA corresponding to the eigenvalue aja_{j}. If aja_{j} is non-degenerate, its eigenspace has dimension 1. Therefore, any eigenvector of AA associated with aja_{j} is linearly dependent on ∣aj⟩\left|a_{j}\right\rangle. Since B∣aj⟩B\left|a_{j}\right\rangle also belongs to this eigenspace, we necessarily have B∣aj⟩=bj∣aj⟩B\left|a_{j}\right\rangle =b_{j}\left|a_{j}\right\rangle.

Now, if aja_{j} is degenerate, that is, if there are two or more linearly independent eigenvectors corresponding to the same eigenvalue aja_{j}, then there exists a subspace of H\mathcal{H} generated by the eigenvectors of AA corresponding to aja_{j}. It is a little more complex, but it is the same idea applied to a subspace of dimension greater than 1. To understand this, we will base our discussion on another source (Massachusetts Institute of Technology (2013)).

For the degenerate case, we then have more vectors associated with aja_{j}. Let us call this subspace Hj\mathcal{H}_{j} and, with {∣ajk⟩}k=1dj\left\{ \left|a^{k}_{j}\right\rangle \right\} ^{d_{j}}_{k=1} being any basis of Hj\mathcal{H}_{j}, then A∣ajk⟩=aj∣ajk⟩A\left|a^{k}_{j}\right\rangle =a_{j}\left|a^{k}_{j}\right\rangle. Now, since AA and BB commute, we have B∣aj⟩∈HjB\left|a_{j}\right\rangle \in\mathcal{H}_{j}, that is, BB transforms vectors from Hj\mathcal{H}_{j} into vectors that also belong to Hj\mathcal{H}_{j}. In particular, BB transforms an eigenvector of AA associated with aja_{j} into a vector that remains in the same eigenspace of AA. If B∣ajk⟩≠0B\left|a^{k}_{j}\right\rangle \neq0, this vector is also an eigenvector of AA associated with aja_{j}.

Therefore, each B∣ajk⟩B\left|a^{k}_{j}\right\rangle can be written as a linear combination of the basis vectors of Hj\mathcal{H}_{j}. Thus, we can consider the action of BB “internally” within the subspace Hj\mathcal{H}_{j}. Since BB is Hermitian, its restriction to Hj\mathcal{H}_{j} is also Hermitian. And since Hj\mathcal{H}_{j} has finite dimension (and is Hermitian), we can diagonalize this restriction of BB in an orthonormal basis of Hj\mathcal{H}_{j}. Therefore, there exists a basis of Hj\mathcal{H}_{j}, which we can denote by {∣ajk⟩}\left\{ \left|a^{k}_{j}\right\rangle \right\}, such that B∣ajk⟩=bj∣ajk⟩B\left|a^{k}_{j}\right\rangle =b_{j}\left|a^{k}_{j}\right\rangle. Since this new basis remains a basis of Hj\mathcal{H}_{j}, its vectors continue to be eigenvectors of AA associated with aja_{j}.

Thus, the eigenvectors are simultaneously eigenvectors of AA and BB. Repeating this procedure for each eigenspace of AA, we obtain an orthonormal basis of H\mathcal{H} consisting of common eigenvectors of the two operators.

In other words, if the Hermitian operators AA and BB commute, then they share the same orthonormal basis of eigenvectors. To better understand this proof, we can represent the operator BB as a matrix in the basis of eigenvectors of AA. Using the completeness relation, I=∑j∣aj⟩⟨aj∣\mathbb{I}=\sum_{j}\left|a_{j}\right\rangle \left\langle a_{j}\right|, we have B=IBI=∑j,k(⟨aj∣B∣ak⟩)(∣aj⟩⟨ak∣)B=\mathbb{I}B\mathbb{I}=\sum_{j,k}\left(\left\langle a_{j}\left|B\right|a_{k}\right\rangle \right)\left(\left|a_{j}\right\rangle \left\langle a_{k}\right|\right).

Thus, the elements of the matrix representing BB in this basis are Baj,ak=⟨aj∣B∣ak⟩B_{a_{j},a_{k}}=\left\langle a_{j}\left|B\right|a_{k}\right\rangle, or in matrix form:

B=[⟨a1∣B∣a1⟩⟨a1∣B∣a2⟩…⟨a2∣B∣a1⟩⟨a2∣B∣a2⟩…⋮⋮⋱]B=\left[\begin{array}{ccc} \left\langle a_{1}\left|B\right|a_{1}\right\rangle & \left\langle a_{1}\left|B\right|a_{2}\right\rangle & \dots\\ \left\langle a_{2}\left|B\right|a_{1}\right\rangle & \left\langle a_{2}\left|B\right|a_{2}\right\rangle & \dots\\ \vdots & \vdots & \ddots \end{array}\right]

Then we can think of BB as a block-diagonal matrix:

B=[B1000B2000⋱]B=\left[\begin{array}{ccc} B_{1} & 0 & 0\\ 0 & B_{2} & 0\\ 0 & 0 & \ddots \end{array}\right]

where each BjB_{j} is a submatrix representing the internal action of BB acting on the subspace Hj\mathcal{H}_{j}. If the eigenvalue aja_{j} is non-degenerate, then the block BjB_{j} is just a 1×11\times1 matrix. If aja_{j} is degenerate, then the block BjB_{j} has dimension dj×djd_{j}\times d_{j}.

Since BB is Hermitian, each block BjB_{j} is also a Hermitian matrix and, therefore, can be diagonalized. Thus, we can diagonalize BB within each subspace Hj\mathcal{H}_{j}, obtaining a basis of eigenvectors of BB that are also eigenvectors of AA. Thus, we prove that if and only if two operators AA and BB commute, then they share a common basis of eigenvectors.

Example 1

In summary, when we have a non-degenerate eigenvalue aja_{j}, the situation is very simple, since in this case, the only eigenvector of AA, up to a multiplicative factor, is automatically also an eigenvector of BB.

When aja_{j} is degenerate, its eigenvectors form a subspace Hj\mathcal{H}_{j}, corresponding to the eigenvalue aja_{j}, in which every nonzero vector is also an eigenvector of AA associated with aja_{j}. Since AA and BB commute, applying BB to any vector in Hj\mathcal{H}_{j} results in another vector belonging to Hj\mathcal{H}_{j}. Thus, it is sufficient to find, within this subspace, the vectors that are also eigenvectors of BB.

This possibility is guaranteed by the fact that BB is Hermitian: its action restricted to Hj\mathcal{H}_{j} can be represented by a Hermitian submatrix, which can be diagonalized by an orthonormal basis of eigenvectors[5]. Therefore, we can choose an orthonormal basis of Hj\mathcal{H}_{j} formed by eigenvectors of BB. Since all vectors in Hj\mathcal{H}_{j} are also eigenvectors of AA associated with aja_{j}, these vectors are common eigenvectors of AA and BB.

Repeating this procedure for each eigenspace of AA, we obtain an orthonormal basis of H\mathcal{H} formed by common eigenvectors of AA and BB.

With the proof of this theorem, we know that if two observables AA and BB of a system do not commute, they cannot share all eigenvectors. If, by measuring the observable AA, we prepare the system in a state that is not an eigenvector of the observable BB, this vector must necessarily be a linear combination of the eigenvectors of with two or more nonzero probability amplitudes.

Example 2

Uncertainty Relation

Here we will use the quantum kinematic formalism to obtain general mathematical relations that provide limitations on the product or sum of the “uncertainties” (standard deviations) of two observables.

Heisenberg-Robertson-Schrödinger Uncertainty Relation (HRSUR)

The HRSUR :

Var(A)Var(B)≥(CovQ(A,B))2+2−2∣⟨[A,B]⟩∣2,\text{Var}\left(A\right)\text{Var}\left(B\right)\ge(\text{CovQ}\left(A,B\right))^{2}+2^{-2}|\langle[A,B]\rangle|^{2},

with the variance defined by Var(A)=⟨(A−⟨A⟩)2⟩.\mathrm{Var}(A)=\langle(A-\langle A\rangle)^{2}\rangle. The quantities on the right-hand side of the inequality will be defined below.

To prove the HRSUR, we will use the Cauchy-Schwarz inequality (CSI) as one of the main tools:

⟨ψ∣ψ⟩⟨ϕ∣ϕ⟩≥∣⟨ψ∣ϕ⟩∣2=⟨ψ∣ϕ⟩⟨ϕ∣ψ⟩,\langle\psi|\psi\rangle\langle\phi|\phi\rangle\ge|\langle\psi|\phi\rangle|^{2}=\langle\psi|\phi\rangle\langle\phi|\psi\rangle,

with equality obtained if and only if ∣ψ⟩|\psi\rangle and ∣ϕ⟩|\phi\rangle are collinear (∣ψ⟩∝∣ϕ⟩|\psi\rangle\propto|\phi\rangle). This result is valid for any pair of vectors ∣ψ⟩,∣ϕ⟩|\psi\rangle,|\phi\rangle in a Hilbert space H.\mathcal{H}. Consider AA and BB to be two observables of a physical system prepared in the state ∣ξ⟩|\xi\rangle. We will use

⟨X⟩:=⟨X⟩ξ=⟨ξ∣X∣ξ⟩\langle X\rangle:=\langle X\rangle_{\xi}=\langle\xi|X|\xi\rangle

to denote the expectation value of any operator XX and II to denote the identity operator in H\mathcal{H}. Before proceeding, I would like to explore this result further. Let us consider a system in the state ∣ψ⟩=∑ncn∣ψ⟩\left|\psi\right\rangle =\sum_{n}c_{n}\left|\psi\right\rangle, and an operator XX with eigenvalues xnx_{n} associated with the eigenvectors ∣xn⟩\left|x_{n}\right\rangle. By the definition of a weighted average, the mean value of the results is ⟨O⟩=∑nonP(on)\left\langle O\right\rangle =\sum_{n}o_{n}P\left(o_{n}\right), and since P(on)=∣cn∣2P\left(o_{n}\right)=\left|c_{n}\right|^{2}, we can write ⟨O⟩=∑non∣cn∣2\left\langle O\right\rangle =\sum_{n}o_{n}\left|c_{n}\right|^{2}. Since ⟨on∣ψ⟩=cn\left\langle o_{n}|\psi\right\rangle =c_{n}, then ⟨O⟩=∑non⟨ψ∣on⟩⟨on∣ψ⟩\left\langle O\right\rangle =\sum_{n}o_{n}\left\langle \psi|o_{n}\right\rangle \left\langle o_{n}|\psi\right\rangle, since ∣⟨ψ∣on⟩∣=⟨on∣ψ⟩⟨ψ∣on⟩\left|\left\langle \psi|o_{n}\right\rangle \right|=\sqrt{\left\langle o_{n}|\psi\right\rangle \left\langle \psi|o_{n}\right\rangle }. Then, by the spectral decomposition O=∑non∣on⟩⟨on∣O=\sum_{n}o_{n}\left|o_{n}\right\rangle \left\langle o_{n}\right|, we have:

⟨O⟩=⟨ψ∣(∑non∣on⟩⟨on∣)∣ψ⟩=⟨ψ∣o∣ψ⟩\left\langle O\right\rangle =\left\langle \psi\right|\left(\sum_{n}o_{n}\left|o_{n}\right\rangle \left\langle o_{n}\right|\right)\left|\psi\right\rangle =\left\langle \psi\left|o\right|\psi\right\rangle

We will also define the following vectors:

∣ψ⟩=(A−⟨A⟩I)∣ξ⟩,∣ϕ⟩=(B−⟨B⟩I)∣ξ⟩,\begin{align} |\psi\rangle & =(A-\langle A\rangle I)|\xi\rangle,\\ |\phi\rangle & =(B-\langle B\rangle I)|\xi\rangle, \end{align}

which we will replace at CSI. Let us note that

⟨ψ∣ψ⟩=((A−⟨A⟩I)∣ξ⟩,(A−⟨A⟩I)∣ξ⟩)=(∣ξ⟩,(A−⟨A⟩I)†(A−⟨A⟩I)∣ξ⟩)=⟨ξ∣(A−⟨A⟩I)2∣ξ⟩=Var(A),⟨ϕ∣ϕ⟩=Var(B),\begin{align} \langle\psi|\psi\rangle & =\left((A-\langle A\rangle I)|\xi\rangle,(A-\langle A\rangle I)|\xi\rangle\right)\\ & =\left(|\xi\rangle,(A-\langle A\rangle I)^{\dagger}(A-\langle A\rangle I)|\xi\rangle\right)\\ & =\langle\xi|(A-\langle A\rangle I)^{2}|\xi\rangle\\ & =\mathrm{Var}(A),\\ \langle\phi|\phi\rangle & =\mathrm{Var}(B), \end{align}

with Var(X)=Var(X)ξ=⟨ξ∣(X−⟨X⟩I)2∣ξ⟩=⟨(X−⟨X⟩I)2⟩\mathrm{Var}(X)=\mathrm{Var}(X)_{\xi}=\langle\xi|(X-\langle X\rangle I)^{2}|\xi\rangle=\langle(X-\langle X\rangle I)^{2}\rangle being the variance of XX. We also use (A−⟨A⟩I)†(A−⟨A⟩I)=(A−⟨A⟩I)2(A-\langle A\rangle I)^{\dagger}(A-\langle A\rangle I)=(A-\langle A\rangle I)^{2} since A=A†A=A^{\dagger} and ⟨A⟩=⟨A⟩∗\left\langle A\right\rangle =\left\langle A\right\rangle ^{*}. We can see that this follows from the fact that AA (and II) are Hermitian, and:

⟨A⟩∗=(⟨ψ∣A∣ψ⟩)∗=⟨ψ∣A†∣ψ⟩=⟨ψ∣A∣ψ⟩=⟨A⟩\begin{align*} \left\langle A\right\rangle ^{*} & =\left(\left\langle \psi\left|A\right|\psi\right\rangle \right)^{*}=\left\langle \psi\left|A^{\dagger}\right|\psi\right\rangle =\left\langle \psi\left|A\right|\psi\right\rangle =\left\langle A\right\rangle \end{align*}

Furthermore, we have to

⟨ψ∣ϕ⟩=((A−⟨A⟩I)∣ξ⟩,(B−⟨B⟩I)∣ξ⟩)=(∣ξ⟩,(A−⟨A⟩I)(B−⟨B⟩I)∣ξ⟩)=(∣ξ⟩,(AB−A⟨B⟩−⟨A⟩B+⟨A⟩⟨B⟩I)∣ξ⟩)=⟨ξ∣(AB−A⟨B⟩−⟨A⟩B+⟨A⟩⟨B⟩I)∣ξ⟩=⟨ξ∣AB∣ξ⟩−⟨ξ∣A∣ξ⟩⟨ξ∣B∣ξ⟩−⟨ξ∣A∣ξ⟩⟨ξ∣B∣ξ⟩+⟨ξ∣A∣ξ⟩⟨ξ∣BI∣ξ⟩=⟨AB⟩−⟨A⟩⟨B⟩−⟨A⟩⟨B⟩+⟨A⟩⟨B⟩=⟨AB⟩−⟨A⟩⟨B⟩\begin{align} \langle\psi|\phi\rangle & =((A-\langle A\rangle I)|\xi\rangle,(B-\langle B\rangle I)|\xi\rangle)\\ & =(|\xi\rangle,(A-\langle A\rangle I)(B-\langle B\rangle I)|\xi\rangle)\\ & =(|\xi\rangle,(AB-A\langle B\rangle-\langle A\rangle B+\langle A\rangle\langle B\rangle I)|\xi\rangle)\\ & =\left\langle \xi\left|\left(AB-A\langle B\rangle-\langle A\rangle B+\langle A\rangle\langle B\rangle I\right)\right|\xi\right\rangle \\ & =\left\langle \xi\left|AB\right|\xi\right\rangle -\left\langle \xi\left|A\right|\xi\right\rangle \left\langle \xi\left|B\right|\xi\right\rangle -\left\langle \xi\left|A\right|\xi\right\rangle \left\langle \xi\left|B\right|\xi\right\rangle +\left\langle \xi\left|A\right|\xi\right\rangle \left\langle \xi\left|BI\right|\xi\right\rangle \\ & =\left\langle AB\right\rangle -\left\langle A\right\rangle \left\langle B\right\rangle -\left\langle A\right\rangle \left\langle B\right\rangle +\left\langle A\right\rangle \left\langle B\right\rangle \\ & =\left\langle AB\right\rangle -\left\langle A\right\rangle \left\langle B\right\rangle \end{align}

And since:

[A,B]=AB−BA,{A,B}=AB+BA\begin{align} \left[A,B\right] & =AB-BA,\\ \left\{ A,B\right\} & =AB+BA \end{align}

are, respectively, the commutator and anticommutator of AA and BB. Summing them:

[A,B]+{A,B}=AB−BA+AB+BA=2AB2−1([A,B]+{A,B})=AB\begin{align*} \left[A,B\right]+\left\{ A,B\right\} & =AB-BA+AB+BA\\ & =2AB\\ 2^{-1}\left(\left[A,B\right]+\left\{ A,B\right\} \right) & =AB \end{align*}
AB=2−1({A,B}+[A,B]).AB=2^{-1}(\{A,B\}+[A,B]).

Then:

⟨ψ∣ϕ⟩=⟨AB⟩−⟨A⟩⟨B⟩=⟨2−1({A,B}+[A,B])⟩−⟨A⟩⟨B⟩=2−1(⟨{A,B}⟩+⟨[A,B]⟩)−⟨A⟩⟨B⟩\begin{align} \left\langle \psi|\phi\right\rangle & =\left\langle AB\right\rangle -\left\langle A\right\rangle \left\langle B\right\rangle \\ & =\left\langle 2^{-1}(\{A,B\}+[A,B])\right\rangle -\left\langle A\right\rangle \left\langle B\right\rangle \\ & =2^{-1}\left(\left\langle \{A,B\}\right\rangle +\left\langle [A,B]\right\rangle \right)-\left\langle A\right\rangle \left\langle B\right\rangle \end{align}

So far, then, we have:

⟨ψ∣ψ⟩⟨ϕ∣ϕ⟩≥∣⟨ψ∣ϕ⟩∣2Var(A)Var(B)≥∣2−1(⟨{A,B}⟩+⟨[A,B]⟩)−⟨A⟩⟨B⟩∣2\begin{align} \left\langle \psi|\psi\right\rangle \left\langle \phi|\phi\right\rangle & \ge\left|\left\langle \psi|\phi\right\rangle \right|{}^{2}\\ \text{Var}\left(A\right)\text{Var}\left(B\right) & \geq\left|2^{-1}\left(\left\langle \{A,B\}\right\rangle +\left\langle [A,B]\right\rangle \right)-\left\langle A\right\rangle \left\langle B\right\rangle \right|{}^{2} \end{align}

Since {A,B}\{A,B\} and [A,B][A,B] are, respectively, Hermitian ({A,B}†={A,B}\{A,B\}^{\dagger}=\{A,B\}) and anti-Hermitian ([A,B]†=−[A,B][A,B]^{\dagger}=-[A,B]) operators, their expectation values are, respectively, purely real and purely imaginary numbers. We have already verified that the expectation value of a Hermitian operator is a real number; now, for an anti-Hermitian operator:

⟨A⟩∗=(⟨ψ∣A∣ψ⟩)∗=⟨ψ∣A†∣ψ⟩=−⟨ψ∣A∣ψ⟩=−⟨A⟩\begin{align*} \left\langle A\right\rangle ^{*} & =\left(\left\langle \psi\left|A\right|\psi\right\rangle \right)^{*}=\left\langle \psi\left|A^{\dagger}\right|\psi\right\rangle =-\left\langle \psi\left|A\right|\psi\right\rangle =-\left\langle A\right\rangle \end{align*}

Now, if ⟨A⟩=a+ib\left\langle A\right\rangle =a+ib, its conjugate ⟨A⟩∗=a−ib\left\langle A\right\rangle ^{*}=a-ib, and if ⟨A⟩∗=−⟨A⟩\left\langle A\right\rangle ^{*}=-\left\langle A\right\rangle, then a+ib=−(a−ib)a+ib=-\left(a-ib\right), that is, a+ib=−a+iba+ib=-a+ib. For this to be equal, we must have a=0a=0. Thus, we are left with ⟨A⟩=ib\left\langle A\right\rangle =ib. Since ⟨A⟩\left\langle A\right\rangle and ⟨B⟩\left\langle B\right\rangle are also real numbers, we can verify this for [A,B]\left[A,B\right] and {A,B}\left\{ A,B\right\}.

(A,B)=AB∓BA(A,B)†=(AB∓BA)†=(AB)†∓(BA)†=B†A†∓A†B†=∓A†B†+B†A†\begin{align*} \left(A,B\right) & =AB\mp BA\\ \left(A,B\right)^{\dagger} & =\left(AB\mp BA\right)^{\dagger}\\ & =\left(AB\right)^{\dagger}\mp\left(BA\right)^{\dagger}\\ & =B^{\dagger}A^{\dagger}\mp A^{\dagger}B^{\dagger}\\ & =\mp A^{\dagger}B^{\dagger}+B^{\dagger}A^{\dagger} \end{align*}

That is, if [A,B]=AB−BA\left[A,B\right]=AB-BA then [A,B]†=−[A,B]\left[A,B\right]^{\dagger}=-\left[A,B\right]. And if {A,B}=AB+BA\left\{ A,B\right\} =AB+BA then {A,B}†={A,B}\left\{ A,B\right\} ^{\dagger}=\left\{ A,B\right\}. Thus, we can write:

⟨ψ∣ϕ⟩=ℜ(⟨ψ∣ϕ⟩)+iℑ(⟨ψ∣ϕ⟩)=[2−1(⟨{A,B}⟩)−⟨A⟩⟨B⟩]+[2−1⟨[A,B]⟩]\begin{align*} \left\langle \psi|\phi\right\rangle & =\Re\left(\left\langle \psi|\phi\right\rangle \right)+i\Im\left(\left\langle \psi|\phi\right\rangle \right)\\ & =\left[2^{-1}\left(\left\langle \left\{ A,B\right\} \right\rangle \right)-\left\langle A\right\rangle \left\langle B\right\rangle \right]+\left[2^{-1}\left\langle \left[A,B\right]\right\rangle \right] \end{align*}

We then have:

ℜ(⟨ψ∣ϕ⟩)=2−1⟨{A,B}−2⟨A⟩⟨B⟩I⟩,ℑ(⟨ψ∣ϕ⟩)=−i2−1⟨[A,B]⟩.\begin{align} \Re(\langle\psi|\phi\rangle) & =2^{-1}\langle\{A,B\}-2\langle A\rangle\langle B\rangle\mathbb{I}\rangle,\\ \Im(\langle\psi|\phi\rangle) & =-i2^{-1}\langle[A,B]\rangle. \end{align}

Note: Since 2−1⟨[A,B]⟩=ib2^{-1}\langle[A,B]\rangle=ib is a purely imaginary number, multiplying by −i-i: −i⋅(ib)=b-i\cdot\left(ib\right)=b, gives us its real coefficient. Thus, considering that

∣⟨ψ∣ϕ⟩∣2=⟨ψ∣ϕ⟩∗⟨ψ∣ϕ⟩=(ℜ(⟨ψ∣ϕ⟩)−iℑ(⟨ψ∣ϕ⟩))(ℜ(⟨ψ∣ϕ⟩)+iℑ(⟨ψ∣ϕ⟩))=ℜ(⟨ψ∣ϕ⟩)2+ℑ(⟨ψ∣ϕ⟩)2\begin{align*} \left|\left\langle \psi|\phi\right\rangle \right|^{2} & =\left\langle \psi|\phi\right\rangle ^{*}\left\langle \psi|\phi\right\rangle \\ & =\left(\Re\left(\left\langle \psi|\phi\right\rangle \right)-i\Im\left(\left\langle \psi|\phi\right\rangle \right)\right)\left(\Re\left(\left\langle \psi|\phi\right\rangle \right)+i\Im\left(\left\langle \psi|\phi\right\rangle \right)\right)\\ & =\Re\left(\left\langle \psi|\phi\right\rangle \right)^{2}+\Im\left(\left\langle \psi|\phi\right\rangle \right)^{2} \end{align*}

We then have:

Var(A)Var(B)≥(ℜ(⟨ψ∣ϕ⟩))2+(ℑ(⟨ψ∣ϕ⟩))2\begin{align*} \mathrm{Var}(A)\mathrm{Var}(B) & \geq(\Re(\langle\psi|\phi\rangle))^{2}+(\Im(\langle\psi|\phi\rangle))^{2} \end{align*}

Now let’s open it:

ℑ(⟨ψ∣ϕ⟩)2=(−i2−1⟨[A,B]⟩)2=−2−2⟨[A,B]⟩2\Im\left(\left\langle \psi|\phi\right\rangle \right)^{2}=\left(-i2^{-1}\langle[A,B]\rangle\right)^{2}=-2^{-2}\langle[A,B]\rangle^{2}

Furthermore, since it is purely imaginary, then:

⟨[A,B]⟩2=(−ib)2=−b2=−∣b∣2=−∣⟨[A,B]⟩∣2\langle[A,B]\rangle^{2}=\left(-ib\right)^{2}=-b^{2}=-\left|b\right|^{2}=-\left|\left\langle \left[A,B\right]\right\rangle \right|^{2}

We are then left with ℑ(⟨ψ∣ϕ⟩)2=2−2∣⟨[A,B]⟩∣2\Im\left(\left\langle \psi|\phi\right\rangle \right)^{2}=2^{-2}\left|\left\langle \left[A,B\right]\right\rangle \right|^{2}. And our inequality:

⟨ψ∣ψ⟩⟨ϕ∣ϕ⟩≥∣⟨ψ∣ϕ⟩∣2Var(A)Var(B)≥ℜ(⟨ψ∣ϕ⟩)2+ℑ(⟨ψ∣ϕ⟩)2Var(A)Var(B)≥ℜ(⟨ψ∣ϕ⟩)2+2−2∣⟨[A,B]⟩∣2\begin{align*} \left\langle \psi|\psi\right\rangle \left\langle \phi|\phi\right\rangle & \ge\left|\left\langle \psi|\phi\right\rangle \right|{}^{2}\\ \text{Var}\left(A\right)\text{Var}\left(B\right) & \geq\Re\left(\left\langle \psi|\phi\right\rangle \right)^{2}+\Im\left(\left\langle \psi|\phi\right\rangle \right)^{2}\\ \text{Var}\left(A\right)\text{Var}\left(B\right) & \geq\Re\left(\left\langle \psi|\phi\right\rangle \right)^{2}+2^{-2}\left|\left\langle \left[A,B\right]\right\rangle \right|^{2} \end{align*}

It now remains for us to demonstrate that:

ℜ(⟨ψ∣ϕ⟩)2=CovQ(A,B)2=(Cov(A,B)+Cov(BA)2)2\Re\left(\left\langle \psi|\phi\right\rangle \right)^{2}=CovQ\left(A,B\right)^{2}=\left(\frac{Cov\left(A,B\right)+Cov\left(BA\right)}{2}\right)^{2}

Where the covariance is Cov(X,Y)=⟨XY⟩−⟨X⟩⟨Y⟩\mathrm{Cov}(X,Y)=\langle XY\rangle-\langle X\rangle\langle Y\rangle between the observables XX and YY, therefore:

ℜ(⟨ψ∣ϕ⟩)2=(⟨AB⟩−⟨A⟩⟨B⟩+⟨BA⟩−⟨B⟩⟨A⟩2)2ℜ(⟨ψ∣ϕ⟩)2=(⟨AB⟩+⟨BA⟩−2⟨B⟩⟨A⟩2)2=(⟨AB+BA⟩2−⟨B⟩⟨A⟩)2=(2−1⟨{AB}⟩−⟨B⟩⟨A⟩)2\begin{align*} \Re\left(\left\langle \psi|\phi\right\rangle \right)^{2} & =\left(\frac{\langle AB\rangle-\langle A\rangle\langle B\rangle+\langle BA\rangle-\langle B\rangle\langle A\rangle}{2}\right)^{2}\\ \Re\left(\left\langle \psi|\phi\right\rangle \right)^{2} & =\left(\frac{\langle AB\rangle+\langle BA\rangle-2\langle B\rangle\langle A\rangle}{2}\right)^{2}\\ & =\left(\frac{\left\langle AB+BA\right\rangle }{2}-\langle B\rangle\langle A\rangle\right)^{2}\\ & =\left(2^{-1}\left\langle \left\{ AB\right\} \right\rangle -\langle B\rangle\langle A\rangle\right)^{2} \end{align*}

And this ℜ(⟨ψ∣ϕ⟩)=2−1⟨{A,B}⟩−⟨A⟩⟨B⟩\Re\left(\left\langle \psi|\phi\right\rangle \right)=2^{-1}\left\langle \left\{ A,B\right\} \right\rangle -\left\langle A\right\rangle \left\langle B\right\rangle is the result we had obtained. We can then write ℜ(⟨ψ∣ϕ⟩)2=CovQ(A,B)2\Re\left(\left\langle \psi|\phi\right\rangle \right)^{2}=\text{CovQ}\left(A,B\right)^{2} and:

Var(A)Var(B)≥ℜ(⟨ψ∣ϕ⟩)2+2−2∣⟨[A,B]⟩∣2Var(A)Var(B)≥CovQ(A,B)2+2−2∣⟨[A,B]⟩∣2\begin{align*} \text{Var}\left(A\right)\text{Var}\left(B\right) & \geq\Re\left(\left\langle \psi|\phi\right\rangle \right)^{2}+2^{-2}\left|\left\langle \left[A,B\right]\right\rangle \right|^{2}\\ \text{Var}\left(A\right)\text{Var}\left(B\right) & \geq\text{CovQ}\left(A,B\right)^{2}+2^{-2}\left|\left\langle \left[A,B\right]\right\rangle \right|^{2} \end{align*}

It is worth noting that if AA and BB are compatible observables, i.e., if [A,B]=0[A,B]=0, then we will have CovQ(A,B)=Cov(A,B)\mathrm{CovQ}(A,B)=\mathrm{Cov}(A,B).

Let us now analyze the problem of the possible triviality of the HRSUR, whose general objective is to determine how small the product of the variances of a pair of observables can be. For this, let us assume that the system is prepared in a state that is an eigenvector of the observable AA, i.e.,

∣ξ⟩=∣a⟩|\xi\rangle=|a\rangle

with A∣a⟩=a∣a⟩A|a\rangle=a|a\rangle and a∈Ra\in\mathbb{R}. In this case, it can be verified that

Var(A)∣a⟩=CovQ(A,B)∣a⟩=⟨[A,B]⟩∣a⟩=0.\mathrm{Var}(A)_{|a\rangle}=\mathrm{CovQ}(A,B)_{|a\rangle}=\langle[A,B]\rangle_{|a\rangle}=0.

And with that, HRSUR gets back to us:

0Var(B)≥0.0\mathrm{Var}(B)\ge0.

Therefore, in this case, the HRSUR provides no information whatsoever about the possible incompatibility of AA and BB. The same type of limitation applies if ∣ξ⟩=∣b⟩|\xi\rangle=|b\rangle for B∣b⟩=b∣b⟩.B|b\rangle=b|b\rangle. In any case, let us verify. For the commutator:

⟨ak∣[A,B]∣ak⟩=⟨ak∣AB∣ak⟩−⟨ak∣BA∣ak⟩=⟨ak∣∑nan∣an⟩⟨an∣B∣ak⟩−⟨ak∣B∑nan∣an⟩⟨an∣∣ak⟩=∑nan[⟨ak∣an⟩⟨an∣B∣ak⟩−⟨ak∣B∣an⟩⟨an∣ak⟩]=ak[⟨ak∣B∣ak⟩−⟨ak∣B∣ak⟩]=0\begin{align*} \left\langle a_{k}\left|\left[A,B\right]\right|a_{k}\right\rangle & =\left\langle a_{k}\left|AB\right|a_{k}\right\rangle -\left\langle a_{k}\left|BA\right|a_{k}\right\rangle \\ & =\left\langle a_{k}\left|\sum_{n}a_{n}\left|a_{n}\right\rangle \left\langle a_{n}\right|B\right|a_{k}\right\rangle -\left\langle a_{k}\left|B\sum_{n}a_{n}\left|a_{n}\right\rangle \left\langle a_{n}\right|\right|a_{k}\right\rangle \\ & =\sum_{n}a_{n}\left[\left\langle a_{k}|a_{n}\right\rangle \left\langle a_{n}\left|B\right|a_{k}\right\rangle -\left\langle a_{k}\left|B\right|a_{n}\right\rangle \left\langle a_{n}|a_{k}\right\rangle \right]\\ & =a_{k}\left[\left\langle a_{k}\left|B\right|a_{k}\right\rangle -\left\langle a_{k}\left|B\right|a_{k}\right\rangle \right]\\ & =0 \end{align*}

For the variance, we will adopt a different strategy: instead of writing out the spectral decomposition, we will use the fact that A∣a⟩=a∣a⟩A|a\rangle=a|a\rangle:

Var(A)∣a⟩=⟨an∣(A−⟨an∣A∣an⟩I)2∣an⟩=⟨an∣(A−an⟨an∣an⟩I)2∣an⟩=⟨an∣(A−anI)2∣an⟩=⟨an∣A2∣an⟩+an2⟨an∣an⟩−2an⟨an∣A∣an⟩=an2+an2−2an2=0\begin{align*} \mathrm{Var}(A)_{|a\rangle} & =\langle a_{n}|(A-\langle a_{n}\left|A\right|a_{n}\rangle I)^{2}|a_{n}\rangle\\ & =\langle a_{n}|(A-a_{n}\langle a_{n}|a_{n}\rangle I)^{2}|a_{n}\rangle\\ & =\langle a_{n}|(A-a_{n}I)^{2}|a_{n}\rangle\\ & =\langle a_{n}|A^{2}|a_{n}\rangle+a^{2}_{n}\langle a_{n}|a_{n}\rangle-2a_{n}\langle a_{n}|A|a_{n}\rangle\\ & =a^{2}_{n}+a^{2}_{n}-2a^{2}_{n}\\ & =0 \end{align*}

And for Quantum Covariance, we’ve already seen that ⟨an∣A∣an⟩=an\langle a_{n}\left|A\right|a_{n}\rangle=a_{n}, ⟨ak∣AB∣ak⟩=⟨ak∣BA∣ak⟩=ak⟨ak∣B∣ak⟩\left\langle a_{k}\left|AB\right|a_{k}\right\rangle =\left\langle a_{k}\left|BA\right|a_{k}\right\rangle =a_{k}\left\langle a_{k}\left|B\right|a_{k}\right\rangle, then:

CovQ(A,B)∣a⟩=⟨ak∣AB∣ak⟩+⟨ak∣BA∣ak⟩−2⟨ak∣B∣ak⟩⟨ak∣A∣ak⟩2=ak⟨ak∣B∣ak⟩+ak⟨ak∣B∣ak⟩−2ak⟨ak∣B∣ak⟩2=0\begin{align*} \mathrm{CovQ}(A,B)_{|a\rangle} & =\frac{\langle a_{k}\left|AB\right|a_{k}\rangle+\langle a_{k}\left|BA\right|a_{k}\rangle-2\langle a_{k}\left|B\right|a_{k}\rangle\langle a_{k}\left|A\right|a_{k}\rangle}{2}\\ & =\frac{a_{k}\left\langle a_{k}\left|B\right|a_{k}\right\rangle +a_{k}\left\langle a_{k}\left|B\right|a_{k}\right\rangle -2a_{k}\langle a_{k}\left|B\right|a_{k}\rangle}{2}\\ & =0 \end{align*}

Robertson and Heisenberg Uncertainty Relations (HUR)

Ignoring[6] the term related to quantum covariance in the HRSUR and taking the square root of the remaining expression, we obtain the Robertson uncertainty relation:

Var(A)Var(B)≥CovQ(A,B)2+2−2∣⟨[A,B]⟩∣2Var(A)Var(B)≥2−2∣⟨[A,B]⟩∣2σAσB≥2−1∣⟨[A,B]⟩∣\begin{align*} \text{Var}\left(A\right)\text{Var}\left(B\right) & \geq\text{CovQ}\left(A,B\right)^{2}+2^{-2}\left|\left\langle \left[A,B\right]\right\rangle \right|^{2}\\ \sqrt{\text{Var}\left(A\right)}\sqrt{\text{Var}\left(B\right)} & \geq\sqrt{2^{-2}\left|\left\langle \left[A,B\right]\right\rangle \right|^{2}}\\ \sigma_{A}\sigma_{B} & \ge2^{-1}\left|\left\langle \left[A,B\right]\right\rangle \right| \end{align*}

with, e.g., σA=Var(A).\sigma_{A}=\sqrt{Var(A)}. When applied to position and linear momentum, it is one of Heisenberg’s uncertainty relations:

σXσP≥2−1∣⟨[X,P]⟩∣=ℏ/2,\sigma_{X}\sigma_{P}\ge2^{-1}|\langle[X,P]\rangle|=\hbar/2,

since [X,P]=iℏI[X,P]=i\hbar I.

Here we are in a position to engage in some further discussion. Heisenberg’s uncertainty relation is constantly the subject of misunderstandings. The first concerns the status of the relation itself. Although it is traditionally called an uncertainty principle, it does not need to be postulated as a principle. As we have seen, Heisenberg’s relation can be obtained as a particular case of the Heisenberg-Robertson-Schrödinger Uncertainty Relation, which, in turn, can be derived from the Cauchy–Schwarz inequality, together with the mathematical properties of observables represented by operators in a Hilbert space.

The second point is that it is incorrect to interpret the uncertainty relation as stating that a particle “does not have a definite position and momentum” or that “we cannot measure position and momentum in the same preparation”. The first statement is a possible interpretation, but it is not the mathematical content of the relation, whereas the second is simply incorrect. The uncertainty relation does not say that it is impossible to obtain position and momentum results in the same preparation; it establishes a relationship between the variances of the results of these measurements.

In other words, we can repeatedly prepare the system in the same state and perform position and momentum measurements. From the results, we can determine the respective distributions and, consequently, their variances. The uncertainty relation establishes that the product of these variances — or, more precisely, of the standard deviations — must satisfy a given relation.

This contrasts with classical mechanics. We can imagine, for example, a particle at rest at a perfectly determined position. Ignoring experimental limitations, it would be possible to conceive of a classical state in which position and velocity were simultaneously well defined, for example, a particle at rest at the origin (x=0,v=0x=0,v=0). In this case, if we repeated the measurements, we would always obtain the same values and, therefore, we would have σx=σv=0\sigma_{x}=\sigma_{v}=0.

Dynamics Postulate

Changes of states in QM are described by a unitary transformation:

∣ψt⟩=Ut∣ψ0⟩,|\psi_{t}\rangle=U_{t}|\psi_{0}\rangle,

i.e.,

UtUt†=UtUt†=I.U_{t}U^{\dagger}_{t}=U_{t}U^{\dagger}_{t}=\mathbb{I}.

The time-evolution operator being obtained from the Schrödinger equation:

iℏ∂tUt=HtUt,i\hbar\partial_{t}U_{t}=H_{t}U_{t},

with the initial condition:

U0=I,U_{0}=\mathbb{I},

where HtH_{t} is the Hamiltonian operator of the system at the instant of time tt,ℏ\hbar is Planck’s constant divided by 2π2\pi and we use the notation ∂t=∂∂t\partial_{t}=\frac{\partial}{\partial_{t}}. The two equations above are equivalent to the Schrödinger equation for the state

iℏ∂t∣ψt⟩=Ht∣ψt⟩.i\hbar\partial_{t}|\psi_{t}\rangle=H_{t}|\psi_{t}\rangle.

We verify this result as follows:

iℏ∂t∣ψt⟩=iℏ∂t(Ut∣ψ0⟩)=(iℏ∂tUt)∣ψ0⟩=HtUt∣ψ0⟩=Ht∣ψt⟩.\begin{align} i\hbar\partial_{t}|\psi_{t}\rangle & =i\hbar\partial_{t}(U_{t}|\psi_{0}\rangle)\\ & =(i\hbar\partial_{t}U_{t})|\psi_{0}\rangle\\ & =H_{t}U_{t}|\psi_{0}\rangle\\ & =H_{t}|\psi_{t}\rangle. \end{align}

Example: When Ht=HH_{t}=H

In this type of situation—i.e., when the Hamiltonian is time-independent—we can verify that the following general form holds:

Ut=e−iHt/ℏ.U_{t}=e^{-iHt/\hbar}.

We will do this first using the spectral decomposition of the Hamiltonian,

H=∑jEj∣Ej⟩⟨Ej∣,H=\sum_{j}E_{j}|E_{j}\rangle\langle E_{j}|,

and applying the operator function:

iℏ∂tUt=iℏ∂texp⁡(−iHtℏ)=iℏ∂texp⁡(∑j(−iEjtℏ)∣Ej⟩⟨Ej∣)\begin{align} i\hbar\partial_{t}U_{t} & =i\hbar\partial_{t}\exp\left(-\frac{iHt}{\hbar}\right)\\ & =i\hbar\partial_{t}\exp\left(\sum_{j}\left(-\frac{iE_{j}t}{\hbar}\right)|E_{j}\rangle\langle E_{j}|\right) \end{align}

Here we should pause to understand the mathematics. We have seen that for any Hermitian operator TT we have exp⁡(T)=∑n=0∞1n!Tn\exp(T)=\sum^{\infty}_{n=0}\frac{1}{n!}T^{n}. First, we can see that H2H^{2}:

H2=(∑jEj∣Ej⟩⟨Ej∣)(∑nEn∣En⟩⟨En∣)=∑j,nEjEn∣Ej⟩⟨Ej∣En⟩⟨En∣=∑jEj2∣Ej⟩⟨Ej∣\begin{align*} H^{2} & =\left(\sum_{j}E_{j}|E_{j}\rangle\langle E_{j}|\right)\left(\sum_{n}E_{n}|E_{n}\rangle\langle E_{n}|\right)=\sum_{j,n}E_{j}E_{n}|E_{j}\rangle\langle E_{j}|E_{n}\rangle\langle E_{n}|=\sum_{j}E^{2}_{j}|E_{j}\rangle\langle E_{j}| \end{align*}

where we used the fact that ⟨Ej∣En⟩=δjn\left\langle E_{j}|E_{n}\right\rangle =\delta_{jn}. In general, we have Hn=∑jEjn∣Ej⟩⟨Ej∣H^{n}=\sum_{j}E^{n}_{j}|E_{j}\rangle\langle E_{j}|. Writing then:

exp⁡(H)=∑n=0∞1n!Hn=∑n=0∞1n!(∑jEjn∣Ej⟩⟨Ej∣)=∑j(∑n=0∞Ejnn!)∣Ej⟩⟨Ej∣\begin{align*} \exp(H) & =\sum^{\infty}_{n=0}\frac{1}{n!}H^{n}\\ & =\sum^{\infty}_{n=0}\frac{1}{n!}\left(\sum_{j}E^{n}_{j}|E_{j}\rangle\langle E_{j}|\right)\\ & =\sum_{j}\left(\sum^{\infty}_{n=0}\frac{E^{n}_{j}}{n!}\right)|E_{j}\rangle\langle E_{j}| \end{align*}

And since exp⁡(x)=∑n=0∞xnn!\exp(x)=\sum^{\infty}_{n=0}\frac{x^{n}}{n!}, then: exp⁡(H)=∑jeEj∣Ej⟩⟨Ej∣\exp(H)=\sum_{j}e^{E_{j}}|E_{j}\rangle\langle E_{j}|. In our case, we then have:

iℏ∂tUt=iℏ∑j∂texp⁡(−iEjtℏ)∣Ej⟩⟨Ej∣=iℏ∑j(−iEjℏ)exp⁡(−iEjtℏ)∣Ej⟩⟨Ej∣=iℏ(−iℏ)∑jexp⁡(−iEjtℏ)Ej∣Ej⟩⟨Ej∣=∑je−iEjtℏH∣Ej⟩⟨Ej=H(∑je−iEjtℏ∣Ej⟩⟨Ej∣)\begin{align} i\hbar\partial_{t}U_{t} & =i\hbar\sum_{j}\partial_{t}\exp\left(-\frac{iE_{j}t}{\hbar}\right)|E_{j}\rangle\langle E_{j}|\\ & =i\hbar\sum_{j}\left(-\frac{iE_{j}}{\hbar}\right)\exp\left(-\frac{iE_{j}t}{\hbar}\right)|E_{j}\rangle\langle E_{j}|\\ & =i\hbar\left(-\frac{i}{\hbar}\right)\sum_{j}\exp\left(-\frac{iE_{j}t}{\hbar}\right)E_{j}|E_{j}\rangle\langle E_{j}|\\ & =\sum_{j}e^{-\frac{iE_{j}t}{\hbar}}H|E_{j}\rangle\langle E_{j}\\ & =H\left(\sum_{j}e^{-\frac{iE_{j}t}{\hbar}}|E_{j}\rangle\langle E_{j}|\right) \end{align}

where we used H∣Ej⟩=Ej∣Ej⟩H|E_{j}\rangle=E_{j}|E_{j}\rangle throughout the entire proof. But now, since

∑je−iEjtℏ∣Ej⟩⟨Ej∣=e−itℏ∑jEj∣Ej⟩⟨Ej∣=e−iHtℏ∣\sum_{j}e^{-\frac{iE_{j}t}{\hbar}}|E_{j}\rangle\langle E_{j}|=e^{-\frac{it}{\hbar}\sum_{j}E_{j}|E_{j}\rangle\langle E_{j}|}=e^{-\frac{iHt}{\hbar}|}

Thus iℏ∂tUt=He−iHtℏ∣i\hbar\partial_{t}U_{t}=He^{-\frac{iHt}{\hbar}|} and since iℏ∂tUt=HtUt,i\hbar\partial_{t}U_{t}=H_{t}U_{t}, therefore Ut=e−iHt/ℏU_{t}=e^{-iHt/\hbar}.

OBS. Note that we could obtain UtU_{t} by solving the Schrödinger equation for each Hamiltonian HH. However, having this general formula, we can begin the solution by calculating the exponential directly, without having to solve the differential equation.

Now we can address another question related to “determinism” as opposed to free will. There are at least two senses of the word determinism. One is more closely associated with classical physics and consists of the idea that, if we know perfectly the current state of a system and the laws governing its evolution, its future and past states are, in principle, determined. The other is a broader sense, of a philosophical nature, which concerns the denial of free will.

Regarding the first sense, the usual formulation of quantum mechanics introduces a probabilistic element that contrasts with the determinism of classical mechanics. The results of measurements are, in general, not uniquely determined by the previous quantum state, but are described by probability distributions. In this sense, the theory contains an indeterministic element. It is worth noting, however, that the Schrödinger equation governing the evolution of the quantum state is deterministic. Given the initial state and the Hamiltonian of the system, the equation uniquely determines the quantum state at any subsequent time. The probabilities of measurement outcomes are then obtained from this state through the Born rule.

Whether nature is ultimately deterministic and the probabilistic nature of quantum mechanics results from some incomplete description, or whether quantum indeterminacy is fundamental, is a question for which we do not have a consensual answer. Different positions regarding this question give rise, in part, to different interpretations of quantum mechanics.

However, it is a mistake to think that this probabilistic mathematical element would, by itself, serve as support for a defense of free will. There is no formulation in the known laws of physics that incorporates free will. We can approximate the chance of obtaining heads or tails when tossing an approximately symmetric coin to 50% for each outcome. This is a probabilistic event, but we would not say, for that reason, that the coin possesses free will. The introduction of probabilities into the description of nature, therefore, does not confer any more free will upon nature than the probabilistic tossing of a coin.

Furthermore, not only is the evolution of the quantum state deterministic, but the probability distribution itself at any given instant imposes constraints on the possible outcomes and their probabilities. In my view, these constraints are incompatible with a conception of “free” will in which the agent should possess unrestricted freedom to determine any of the possible outcomes.

Footnotes
  1. A postulate is a proposition assumed to be true within a theory, serving as a starting point for the construction of its concepts and results.

  2. A metric space is, essentially, a set of objects in which we define a way of measuring the distance between two objects. In the case of a vector space with an inner product, the distance can be defined from the norm induced by the inner product.

  3. A degree of freedom corresponds to an independent physical variable that can vary and that can be used to characterize the state of a system. For example, the position and spin of a particle correspond to different degrees of freedom.

  4. A non-degenerate eigenvalue is one to which a single eigenvector corresponds, up to a global phase, that is, a complex phase factor eiϕe^{i\phi} that multiplies the state vector and does not affect the measurement outcomes, since eiϕ(eiϕ)∗=1e^{i\phi}\left(e^{i\phi}\right)^{*}=1.

  5. The fact that a matrix is Hermitian guarantees that we can choose, in the space on which it acts, an orthonormal basis formed by its eigenvectors (Delft University of Technology (2026), on this page.)

  6. That is, if a≥b+c≥ca\geq b+c\geq c. With b≥0b\geq0, then necessarily a≥ca\geq c.

References
  1. Griffiths, D. J. (2017). Introduction to Quantum Mechanics. Cambridge University Press. https://books.google.com.br/books?id=0h-nDAAAQBAJ
  2. Siddiqui, S. (2014). Why are differential equations used for expressing the laws of physics? 10.48550/arXiv.1406.1112
  3. Maziero, J. (2026). Notas de aula de Mecânica Quântica. https://github.com/jonasmaziero/mecanica_quantica/tree/main/notas_de_aula
  4. Massachusetts Institute of Technology. (2013). Quantum Physics I. MIT OpenCourseWare, 8.04 Quantum Physics I. https://www.ocw.mit.edu/courses/8-04-quantum-physics-i-spring-2013/9ccc9fd132124dd102e1a7863e6e12ce_MIT8_04S13_OnCommEigenbas.pdf
  5. Delft University of Technology. (n.d.). Mathematics for Quantum Physics. Retrieved August 31, 2026, from https://interactivetextbooks.tudelft.nl/mqp-v2/