I outline the construction of Dyson Brownian motion which governs the evolution of the eigen-values of an -dimensional stochastic process of Hermitian matrices. For instance, if is such a process, then:
where in an -dimensional random Hermitian matrix drawn from the Gaussian Unitary Ensemble ().
Why study the eigen-values of a stream of Hermitian matrices? At first glance, this seems like a rather obscure mathematical object. Before I can answer this question, in Section 2 I first define what a Hermitian matrix is and discuss how I would select this random matrix . Then, in Section 3 I can give some practical examples in which Dyson Brownian motion would be a useful construction. I also give an alternative interpretation of Dyson Brownian motion related to non-intersecting Brownian processes and explain the use of complex matrices in an economic context. Finally, in Section 4, I construct Dyson Brownian motion.
2. Mathematical Foundation
First things first: “What is a Hermitian matrix?”
Definition (Hermitian Matrix): A square matrix is called Hermitian if it is self-adjoint:
Each element of a Hermitian matrix is the complex conjugate of element in . Thus, the diagonal elements have to be real. Consider an example Hermitian matrix :
Hermitian matrices are just the complex extension of real symmetric matrices. For instance, the matrix below is a real instance of a Hermitian matrix:
Next, when I defined the stochastic process above, I characterized each lurch forward by the addition of a random matrix scaled by the square root of the time interval. I pull this random matrix from the Gaussian Unitary Ensemble.
Definition (Gaussian Unitary Ensemble): The Gaussian Unitary Ensemble is a probability space over the vector space of dimensional Hermitian matrices governed by the measure defined as:
So, each upper triangular element will be drawn from while each element along the diagonal will be drawn from . To get a feel for what this definition really means, consider a concrete example. Suppose that we want to know the probability for the example above. Using these paramters, I find that:
The analysis below will follow through if you consider only adding random matrices drawn from an ensemble of real symmetric matrices. This ensemble is known as the Gaussian Orthogonal Ensemble.
Note that each of the elements of the matrices do not follow their own independent Brownian motion processes. In a loose sense, this would be “too little” structure. The main thrust of the Dyson Brownian motion construction is that the eigen-values of this process follow a Brownian motion plus a twist term. The eigen-values are an attractive summary statistic for reasons. First, we know that they have a simple spectrum due to the fact that at each point in time is a Hermitian matrix. What’s more, the Ky Fan inequality tells us that the eigen-values should have a smooth transition function over time.
3. Motivating Examples
With these terms in hand, I can now ask: “Why worry about matrix valued stochastic processes?”
First, consider some financial applications. Financial theory is built around the -pricing models which measure the correlation between the returns of assets and various risk factors. Models which assume that these measures remain constant for long periods of time do poorly in empirical tests.1 It would be nice to characterize the evolution of these variance-covariance matrices.2 Alternatively, suppose that you were in the business of identifying the principle components of stock returns–either explicitly or by hunting for additional factors.3 Here, you might want to check to work out how likely it is that the largest principle component has moved by in order to test your model.
In an entirely different context, Dyson Brownian motion can also be thought of as characterizing the evolution of Brownian motions that have been restricted to never intersect. The problem of modelling the eigen-values of Hermitian matrices and non-intersecting Brownian processes are not link a priori. However, constructing these non-intersecting processes is hard and an elegant solution emerges via solving the eigen-value process problem for Hermitian matrices which have a simple spectrum. For an example of the economic usefulness of such a trick, conside modelling the real option of a worker to switch jobs.4 At each point in time, he has a next best option but the exact nature of that next best option will change over time. Rather than keeping track of all possibilities, you could just model the evolution of the best option via Dyson Brownian motion.
Finally, I want to make a quick note about the use of complex valued rather than real matrices. Physicists declare that real symmetric matrices preserve time reversal symmetry while complex Hermitian matrices do not. For instance, in the financial application above where each of the entries in the variance-covariance matrix process have to be real, you can always undo the last step of the stochastic process by hitting by a well chosen inverse. However, when using complex valued matrices, this inverse is no longer possible as complex numbers are periodic; i.e.,
To see the implications of this fact in a macroeconomic setting, consider a complex valued extension of a Leontif production model as follows. Suppose that prices are local5 and form the real part of each entry while the magnitude of the transaction forms the complex part. So, for instance, a transaction of tons of steel between a builder and a steel maker at a price of dollars per ton would manifest itself as and . Thus, by introducing an additional dimension to the Leontif matrix, the physical properties of the process it represents changes dramatically.
Finally, I actually get around to defining Brownian motion. Below, I state the result:
Theorem (Dyson Brownian Motion): Let , and be the spectrum of eigen-values of the Hermitian matrix valued process . Then, we have:
for all , where and are independent Brownian motion processes.
In words, this theorem says that the eigen-values of a stochastic process of Hermitian matrices behave like independent Brownian motions plus a repulsion force which is inversely proportional to the distances between any eigen-values. What’s more, this repulsive force is not-localized. Each eigen-value is pushed an pulled by and but also and as well as each and every eigen-level in between.
To formulate the construction, I use a Lemma from Hadamard given below:
Lemma (Hadamard Operator): The eigen-values of have the following first and second derivatives with respect to time :
Proof (Hadamard Operator):
Finally, I use this lemma together with the properties of Hermitian matrices to finish the construction of Dyson Brownian motion.
Proof (Dyson Brownian Motion):