Posted: April 24th, 2025
The Exploration and Exploration Discussion provide students with practical motivation, and insight into common application areas.
Within the Module 3 Exploration, after reading through the Exploration tabs, and checking your knowledge in the Check Understanding tab, click on the Discuss! tab. Choose one of the applications, and one of the numbered Discussion Topics for that application to explore with your classmates.
Module 03/M3_Exploration.html
Exploration of Vector Spaces and Subspaces
Introduction
Crystal Lattices
Markov Chains
Check Understanding
Discuss!
References
Vector spaces and subspaces
Outcomes
Recognize the general vector space framework with basic definitions and problems. (D)
Demonstrate how closely other vector spaces resemble Rn. (D)
Use vector space terminology to tie together important facts about rectangular matrices. (A, B, D)
Apply theory to Markov Chains models in biology, business, and engineering. (A, B, C, D, I)
In these Explorations, we will explore a few introductory applications of vector spaces and subspaces. Click on the tabs to explore each application, and read and watch all of the examples. Then, answer the questions on the Check Understanding tab.
The Discuss! tab contains a few discussion questions about the applications in this module, for use in the Exploration Discussion. The Exploration Discussion will provide you with an opportunity to expand further upon the ideas presented here.
Crystal Lattices
Note: If equations don’t render correctly when initially viewed, refresh your browser page and wait for it to load completely before visiting any of the tabs.
Some applications have a physically natural coordinate system that may be different from the standard coordinate system with standard unit basis vectors in Rn{“version”:”1.1″,”math”:”Rn”}.
In crystallography, a crystal lattice in R3{“version”:”1.1″,”math”:”R3″} is composed of unit cells, which are minimal parallelepipeds, that when translated by each of the lattice vectors along three adjacent edges, make up the crystal lattice, as shown.
Fig. 3.1.1 (generated from): Averill, B., & Eldredge, P. (2007). Chemistry: Principles, patterns, and applications. San Francisco, CA: Pearson Benjamin Cummings.
Figure 3.1.1: A three-dimensional unit cell with its three lattice vectors, and the resulting crystal lattice.
The arrangement of atoms within a unit cell is repeated throughout the crystal lattice. Therefore, knowledge about a crystal lattice can be reduced to knowledge about the unit cell and the arrangement of atoms within it. There are 14 basic types of unit cells comprising all known crystal lattice structures.
Scientists use the lattice vectors that make up the adjacent edges of a unit cell as their working basis of R3{“version”:”1.1″,”math”:”R3″}, forming a natural coordinate system. This gives the atomic arrangement within each unit cell its natural coordinates, whereas using the standard unit basis and resulting unit cubes might be useful for some calculations, but unintuitive for others, since they do not fit the crystal’s natural structure.
By using this natural coordinate system, scientists can combine atom configurations with formulas derived from the relevant type of unit cell to determine such things about a given material as: in what direction it is likely to fracture, where it is easiest to make a clean cut, in what direction it has the greatest tensile strength, or the greatest compressive resistance, and so on. These directions are naturally aligned with the directions of the lattice vectors, so using the lattice vectors as a basis is not only natural, but also a very useful choice.
Example 3.1.1: Finding a Coordinate Vector of an Atom Relative to a Lattice Basis for Graphite
A lattice basis for graphite is given by LB=a3/2-1/20, a3/21/20, c001{“version”:”1.1″,”math”:”LB=a3/2-1/20, a3/21/20, c001″}, where a=2.46 Ao{“version”:”1.1″,”math”:”a=2.46 Ao”} (angstroms), and c=6.71 Ao{“version”:”1.1″,”math”:”c=6.71 Ao”}.
An atom within a unit cell is located at v=-a/23-a/2c/2{“version”:”1.1″,”math”:”v=-a/23-a/2c/2″}. The following video shows how to find the coordinates of this atom relative to the lattice basis.
Video 3.1.1: Finding the Coordinates of an Atom Relative to a Lattice Basis
Your browser does not support the video tag.
Video by Dr. Lisa Korf: Finding the Coordinates of an Atom Relative to a Lattice Basis [© CCCOnline]
Video 3.1.1 Transcript: Finding the Coordinates of an Atom Relative to a Lattice Basis
It is also important to be able to convert from a lattice basis back to the standard unit basis or to another basis that may be useful for other computations. One such basis is called the reciprocal lattice basis.
Example 3.1.2: Converting from the Lattice Basis for Graphite to the Reciprocal Lattice basis
We begin again with the lattice basis for graphite given by LB=a3/2-1/20, a3/21/20, c001{“version”:”1.1″,”math”:”LB=a3/2-1/20, a3/21/20, c001″}.
This time, we are given a coordinate vector with respect to the lattice basis, vLB=-201{“version”:”1.1″,”math”:”vLB=-201″}, and we want to know the coordinates of v{“version”:”1.1″,”math”:”v”}
with respect to the reciprocal lattice basis given by RLB=2πa1/3-10, 2πa1/310, 2πc001.
Video 3.1.2: Converting Coordinates Between Two Lattice Bases
Your browser does not support the video tag.
Video by Dr. Lisa Korf: Converting the Coordinates Between Two Lattice Basis [© CCCOnline]
Video 3.1.2 Transcript: Converting Coordinates Between Two Lattice Bases
Crystallography has many areas of application, such as to metallurgy, solid state electronics, and geology. In all of these fields, it is important to be able to perform change of basis operations from one coordinate system to another.
Markov Chains
Consider a system that can be in exactly one of n possible states. Suppose that the system evolves in time such that the probability of being in a particular state at a given time period k depends only on the system’s state at the previous time period, k-1{“version”:”1.1″,”math”:”k-1″}.
For example, suppose you present your pet cat, Cosmo, with a Box Game, in which there are three boxes that she can open. Inside box 1 is a cat-chasing Chihuahua, while box 2 is empty, and box 3 contains an enticing mouse. If Cosmo opens box 1 and finds the Chihuahua, there is only a 2% probability that she will open that box the next time. If she opens box 2, finding it empty, there is a 20% chance she will open that box the next time, and if she opens box 3, there is a 90% chance she will open that box for a repeat game of cat and mouse. We’ll assume that Cosmo only remembers what happened the previous time she opened a box, but not any of the times prior to that.
Such a scenario can be modeled with a Markov Chain.
Definition 3.2.1. A probability vector is a vector of states whose entries sum to 1. A stochastic matrix is a square matrix with probability vectors as its columns, i.e. the entries of every column sum to one. A Markov Chain is a sequence of probability vectors x0, x1, …, xk, xk+1, …{“version”:”1.1″,”math”:”x0, x1, …, xk, xk+1, …”}, together with a stochastic matrix P, such thatx1=Px0, x2=Px1, …, xk+1=Pxk, ….The probability vectors represent a sequence of states, and the stochastic matrix P has entries Pij that give the transition probabilities of going from state j to state i.
Example 3.2.1: The Stochastic Matrix and Probability Vectors for the Cat and Box Game
If Cosmo the cat picks box 1 on her first try, then the probability that she will pick box 1 on her second try is only 2%, so the stochastic matrix P for this example would have entry P11=0.02. Assuming that the remaining probability of 0.98 is split evenly between the remaining boxes in that column, and similarly for the remaining columns, we can write down the stochastic matrix for Cosmo: From Box: P= 1 2 3 0.020.40.050.490.20.050.490.40.9To Box:123If Cosmo has a 1/3 probability of opening each of the three boxes at time period 0, her initial state vector would be x0=1/31/31/3{“version”:”1.1″,”math”:”x0=1/31/31/3″}. The following video computes the first 6 state vectors that result after Cosmo chooses which box to open at 6 consecutive time periods.
Video 3.2.1: Computing the First 6 State Vectors of a Markov Chain
Your browser does not support the video tag.
Video by Dr. Lisa Korf: Computing the First 6 State Vectors of a Markov Chain [© CCCOnline]
Video 3.2.1 Transcript: Computing the First 6 State Vectors of a Markov Chain
From the video, it appears that as Cosmo continues to choose boxes, her state vector starts to hone in on a single vector. This is not a coincidence, due to the following theorem.
Definition 3.2.2. A stochastic matrix P is regular if there is any power of P whose entries are all strictly positive.
Theorem 3.2.1. If P is a regular stochastic matrix, then there is a unique probability vector, q, such that Pq = q. The vector q is called the steady state vector for P. In addition, no matter what the initial probability vector x0 is, the sequence of probability vectors {xk} will converge to q as k tends to infinity.
Example 3.2.2: The Steady State Vector for the Cat and Box Game
To find a steady state vector q, we know from the theorem above that it must be the unique probability vector that solves the equation Px=x. So let’s use linear algebra to find Cosmo’s steady state vector!
Video 3.2.2: Finding the Steady State Vector
Your browser does not support the video tag.
Video by Dr. Lisa Korf: Finding the Steady State Vector [© CCCOnline]
Video 3.2.2 Transcript: Finding the Steady State Vector
We can now see that in the long run, Cosmo will pick box 1 with probability 40/481, or about 8.316% of the time, box 2 with probability 49/481, or about 10.187% of the time, and box 3 with probability 392/481, or about 81.497% of the time. Based on these results, do you think that in reality, Cosmo’s memory would really only go back one time period, or would she be more likely to avoid box 1, perhaps indicating a longer-term memory?
Image Attributions: Dog Silhouette: https://pixabay.com/en/dog-puppy-animal-black-silhouette-159481/ , Cat Silhouette: http://www.clipartbest.com/clipart-ndT88AKie , Mouse Silhouette: (mirrored image) http://www.rgbstock.com/bigphoto/mSdQGMK/silhouette+gerbil
check your understanding
Crystal Lattices:
Suppose that a crystal lattice has a lattice basis given by LB = { u, v, w }, and that an atom is located at x = u + 3v + w. What are the coordinates of the atom with respect to the lattice basis LB?
See Answer
Answer: xLB=131{“version”:”1.1″,”math”:”xLB=131″}
A crystal lattice has a particular atom type appearing at coordinates (1, 2, –3), and repeating in three different directions at coordinates (–1, 2, 4), (2, 3, 5), and (4, –2, 1). What is a possible lattice basis for this crystal lattice?
See Answer
Answer: -1-12-24-(-3), 2-13-25-(-3), 4-1-2-21-(-3)=-207, 118, 3-44{“version”:”1.1″,”math”:”-1-12-24-(-3), 2-13-25-(-3), 4-1-2-21-(-3)=-207, 118, 3-44″}
Markov Chains
:
Suppose that Cosmo’s initial probability vector had for some reason been 001{“version”:”1.1″,”math”:”001″}, with all other probabilities the same. How would this affect Cosmo’s steady state vector?
See Answer
Answer: It wouldn’t, since when a steady state vector exists, which it does in this case, it is independent of the initial probability vector state.
For a Markov chain with stochastic matrix equal to the 2 by 2 matrix, P=0110{“version”:”1.1″,”math”:”P=0110″}, does a steady state vector exist? If so, what is it? If not, why not?
See Answer
Answer: Although this matrix does not satisfy the assumptions of the theorem, such a matrix may still have a steady state if it has a unique probability vector solution to the equation P-Ix=0{“version”:”1.1″,”math”:”P-Ix=0″}. In this case, P-I=-111-1{“version”:”1.1″,”math”:”P-I=-111-1″}, which has reduced row-echelon form, 1-100{“version”:”1.1″,”math”:”1-100″}. The solution to the homogeneous system is x=x211{“version”:”1.1″,”math”:”x=x211″}. Letting x2=1/2{“version”:”1.1″,”math”:”x2=1/2″} yields the unique probability vector solution, 1/21/2{“version”:”1.1″,”math”:”1/21/2″}, which is also the steady state vector.
Discuss!
Choose one of the applications, and one of the numbered Discussion Topics for that application to explore with your classmates in the Exploration Discussion in this Module.
Crystal Lattices – Discussion Topics
Using the lattice basis LB and the reciprocal lattice basis RLB given for graphite in Examples 3.1.1 and 3.1.2, you wish to determine the location of a point of interest within the lattice given by vRLB=abc{“version”:”1.1″,”math”:”vRLB=abc”}, where a, b,{“version”:”1.1″,”math”:”a, b,”}and c{“version”:”1.1″,”math”:”c”} are integers, not all zero, to be chosen by you and such that the vector is different from that chosen in any other posts on this topic . Give the coordinates of this point both with respect to the standard basis vectors for R3, and with respect to the lattice basis LB. Explain your steps carefully. Which basis of the three seems most natural for this point of interest? Be sure to include your chosen values for a, b,{“version”:”1.1″,”math”:”a, b,”}and c{“version”:”1.1″,”math”:”c”} in the Subject Heading of your post.
Diamond is a crystalline structure with a cubic unit cell. The lattice basis for diamond is given by LB=a2011, a2101, a2110, where a is a lattice constant. The reciprocal lattice basis for diamond is given by RLB=2πa-111, 2πa1-11, 2πa11-1. Using the appropriate 3 by 6 augmented matrices and row reduction techniques, find the change of coordinates matrices, PLB←RLB, PR3←RLB, and PLB←R3, {“version”:”1.1″,”math”:”PLB←RLB, PR3←RLB, and PLB←R3, ”}and show that one of these matrices may be obtained as a product of the other two.
Markov Chains – Discussion Topics
As indicated at the end of the Markov Chain discussion, realistically it would be more likely that Cosmo the cat would remember more than just one time period back. Assume that Cosmo recalls the last two time periods. To set up a Markov Chain in that setting, one needs to set up new states (probability vectors) that take the last two time periods from the original setting into account. So the state might be given by x=x1,1×1,2×1,3×2,1×2,2×2,3×3,1×3,2×3,3,where x(i, j) is the component of the probability vector that gives the probability that Cosmo will choose box i in this time period, and box j in the following time period. To set up a Markov chain, one possible stochastic matrix P is given by From: P= (1,1) (1,2) (1,3) (2,1) (2,2) (2,3) (3,1) (3,2) (3,3) 0.02000.21000.035000.49000.345000.27000.49000.445000.6950000.21000.4000.225000.345000.2000.125000.445000.4000.650000.035000.225000.05000.27000.125000.05000.695000.65000.9To:(1,1)(1,2)(1,3)(2,1)(2,2)(2,3)(3,1)(3,2)(3,3) where the (a,b), (b,c) entry corresponds to the probability that Cosmo will choose box c on the third try if she chooses box a on the first try and box b on the second try.Why do some of the entries have to be 0 in this new setting? Set up your own stochastic matrix for Cosmo with a two period memory like this one. The only entries that should be different are in the (1,1) column, the (2,2) column, and the (3,3) column. Get inside Cosmo’s head and decide what her probability is of choosing box 1 after having chosen it twice already, and similarly for boxes 2 and 3. This will give you the values for the (1,1), (1,1) entry, (2,2), (2,2) entry, and the (3,3), (3,3) entry. The remaining two nonzero values in each of those columns should be equal to each other, and the sum of the three nonzero entries in each column should sum to 1, since it is a stochastic matrix. For example, you might choose the value 0.01 (instead of 0.02) for the (1,1), (1,1) entry. In that case, the remaining two nonzero entries would both be equal to 0.99/2 = 0.495. For your stochastic matrix setup, what would Cosmo’s steady state vector q be? What is the interpretation of each of the 9 entries in q? What is Cosmo’s total probability of choosing box 1 in the long run, and how did you obtain it? Is this more realistic? Be sure to include your chosen values for the (1,1), (1,1) , (2,2), (2,2) , and (3,3), (3,3) entries in the Subject Heading of your main post, since they will need to be different from the values chosen in any other posts on this topic.
References
Crystal Lattices:
Chung, D. D. (2002). Review Graphite. Journal of Materials Science, 37(8), 1475-1489.
Averill, B., & Eldredge, P. (2007). Chemistry: Principles, patterns, and applications. San Francisco, CA: Pearson Benjamin Cummings.
Ungersböck, S.E. (1977, August 8). 3.4 Basic Properties of the Diamond Structure, Dissertation: Advanced Modeling of Strained CMOS Technology. Retrieved April 18, 2016, from
http://www.iue.tuwien.ac.at/phd/ungersboeck/node27.html
Markov Chains:
Khoury, J. (2003, October 28). Application to Markov Chains. Retrieved April 09, 2016, from
http://aix1.uottawa.ca/~jkhoury/markov.htm
Linear Algebra/Topic: Markov Chains. (2010, May 26). Retrieved April 09, 2016, from
https://en.wikibooks.org/wiki/Linear_Algebra/Topic:_Markov_Chains
Markov Chains, Math 22, Spring 2007. (2007). Retrieved April 9, 2016, from
https://math.dartmouth.edu/archive/m22s07/public_html/MarkovSlides
Introduction
Crystal Lattices
Markov Chains
Check Understanding
Discuss!
References
Module 03/M3_Exploration.html
Exploration of Vector Spaces and Subspaces
Introduction
Crystal Lattices
Markov Chains
Check Understanding
Discuss!
References
Vector spaces and subspaces
Outcomes
Recognize the general vector space framework with basic definitions and problems. (D)
Demonstrate how closely other vector spaces resemble Rn. (D)
Use vector space terminology to tie together important facts about rectangular matrices. (A, B, D)
Apply theory to Markov Chains models in biology, business, and engineering. (A, B, C, D, I)
In these Explorations, we will explore a few introductory applications of vector spaces and subspaces. Click on the tabs to explore each application, and read and watch all of the examples. Then, answer the questions on the Check Understanding tab.
The Discuss! tab contains a few discussion questions about the applications in this module, for use in the Exploration Discussion. The Exploration Discussion will provide you with an opportunity to expand further upon the ideas presented here.
Crystal Lattices
Note: If equations don’t render correctly when initially viewed, refresh your browser page and wait for it to load completely before visiting any of the tabs.
Some applications have a physically natural coordinate system that may be different from the standard coordinate system with standard unit basis vectors in Rn{“version”:”1.1″,”math”:”Rn”}.
In crystallography, a crystal lattice in R3{“version”:”1.1″,”math”:”R3″} is composed of unit cells, which are minimal parallelepipeds, that when translated by each of the lattice vectors along three adjacent edges, make up the crystal lattice, as shown.
Fig. 3.1.1 (generated from): Averill, B., & Eldredge, P. (2007). Chemistry: Principles, patterns, and applications. San Francisco, CA: Pearson Benjamin Cummings.
Figure 3.1.1: A three-dimensional unit cell with its three lattice vectors, and the resulting crystal lattice.
The arrangement of atoms within a unit cell is repeated throughout the crystal lattice. Therefore, knowledge about a crystal lattice can be reduced to knowledge about the unit cell and the arrangement of atoms within it. There are 14 basic types of unit cells comprising all known crystal lattice structures.
Scientists use the lattice vectors that make up the adjacent edges of a unit cell as their working basis of R3{“version”:”1.1″,”math”:”R3″}, forming a natural coordinate system. This gives the atomic arrangement within each unit cell its natural coordinates, whereas using the standard unit basis and resulting unit cubes might be useful for some calculations, but unintuitive for others, since they do not fit the crystal’s natural structure.
By using this natural coordinate system, scientists can combine atom configurations with formulas derived from the relevant type of unit cell to determine such things about a given material as: in what direction it is likely to fracture, where it is easiest to make a clean cut, in what direction it has the greatest tensile strength, or the greatest compressive resistance, and so on. These directions are naturally aligned with the directions of the lattice vectors, so using the lattice vectors as a basis is not only natural, but also a very useful choice.
Example 3.1.1: Finding a Coordinate Vector of an Atom Relative to a Lattice Basis for Graphite
A lattice basis for graphite is given by LB=a3/2-1/20, a3/21/20, c001{“version”:”1.1″,”math”:”LB=a3/2-1/20, a3/21/20, c001″}, where a=2.46 Ao{“version”:”1.1″,”math”:”a=2.46 Ao”} (angstroms), and c=6.71 Ao{“version”:”1.1″,”math”:”c=6.71 Ao”}.
An atom within a unit cell is located at v=-a/23-a/2c/2{“version”:”1.1″,”math”:”v=-a/23-a/2c/2″}. The following video shows how to find the coordinates of this atom relative to the lattice basis.
Video 3.1.1: Finding the Coordinates of an Atom Relative to a Lattice Basis
Your browser does not support the video tag.
Video by Dr. Lisa Korf: Finding the Coordinates of an Atom Relative to a Lattice Basis [© CCCOnline]
Video 3.1.1 Transcript: Finding the Coordinates of an Atom Relative to a Lattice Basis
It is also important to be able to convert from a lattice basis back to the standard unit basis or to another basis that may be useful for other computations. One such basis is called the reciprocal lattice basis.
Example 3.1.2: Converting from the Lattice Basis for Graphite to the Reciprocal Lattice basis
We begin again with the lattice basis for graphite given by LB=a3/2-1/20, a3/21/20, c001{“version”:”1.1″,”math”:”LB=a3/2-1/20, a3/21/20, c001″}.
This time, we are given a coordinate vector with respect to the lattice basis, vLB=-201{“version”:”1.1″,”math”:”vLB=-201″}, and we want to know the coordinates of v{“version”:”1.1″,”math”:”v”}
with respect to the reciprocal lattice basis given by RLB=2πa1/3-10, 2πa1/310, 2πc001.
Video 3.1.2: Converting Coordinates Between Two Lattice Bases
Your browser does not support the video tag.
Video by Dr. Lisa Korf: Converting the Coordinates Between Two Lattice Basis [© CCCOnline]
Video 3.1.2 Transcript: Converting Coordinates Between Two Lattice Bases
Crystallography has many areas of application, such as to metallurgy, solid state electronics, and geology. In all of these fields, it is important to be able to perform change of basis operations from one coordinate system to another.
Markov Chains
Consider a system that can be in exactly one of n possible states. Suppose that the system evolves in time such that the probability of being in a particular state at a given time period k depends only on the system’s state at the previous time period, k-1{“version”:”1.1″,”math”:”k-1″}.
For example, suppose you present your pet cat, Cosmo, with a Box Game, in which there are three boxes that she can open. Inside box 1 is a cat-chasing Chihuahua, while box 2 is empty, and box 3 contains an enticing mouse. If Cosmo opens box 1 and finds the Chihuahua, there is only a 2% probability that she will open that box the next time. If she opens box 2, finding it empty, there is a 20% chance she will open that box the next time, and if she opens box 3, there is a 90% chance she will open that box for a repeat game of cat and mouse. We’ll assume that Cosmo only remembers what happened the previous time she opened a box, but not any of the times prior to that.
Such a scenario can be modeled with a Markov Chain.
Definition 3.2.1. A probability vector is a vector of states whose entries sum to 1. A stochastic matrix is a square matrix with probability vectors as its columns, i.e. the entries of every column sum to one. A Markov Chain is a sequence of probability vectors x0, x1, …, xk, xk+1, …{“version”:”1.1″,”math”:”x0, x1, …, xk, xk+1, …”}, together with a stochastic matrix P, such thatx1=Px0, x2=Px1, …, xk+1=Pxk, ….The probability vectors represent a sequence of states, and the stochastic matrix P has entries Pij that give the transition probabilities of going from state j to state i.
Example 3.2.1: The Stochastic Matrix and Probability Vectors for the Cat and Box Game
If Cosmo the cat picks box 1 on her first try, then the probability that she will pick box 1 on her second try is only 2%, so the stochastic matrix P for this example would have entry P11=0.02. Assuming that the remaining probability of 0.98 is split evenly between the remaining boxes in that column, and similarly for the remaining columns, we can write down the stochastic matrix for Cosmo: From Box: P= 1 2 3 0.020.40.050.490.20.050.490.40.9To Box:123If Cosmo has a 1/3 probability of opening each of the three boxes at time period 0, her initial state vector would be x0=1/31/31/3{“version”:”1.1″,”math”:”x0=1/31/31/3″}. The following video computes the first 6 state vectors that result after Cosmo chooses which box to open at 6 consecutive time periods.
Video 3.2.1: Computing the First 6 State Vectors of a Markov Chain
Your browser does not support the video tag.
Video by Dr. Lisa Korf: Computing the First 6 State Vectors of a Markov Chain [© CCCOnline]
Video 3.2.1 Transcript: Computing the First 6 State Vectors of a Markov Chain
From the video, it appears that as Cosmo continues to choose boxes, her state vector starts to hone in on a single vector. This is not a coincidence, due to the following theorem.
Definition 3.2.2. A stochastic matrix P is regular if there is any power of P whose entries are all strictly positive.
Theorem 3.2.1. If P is a regular stochastic matrix, then there is a unique probability vector, q, such that Pq = q. The vector q is called the steady state vector for P. In addition, no matter what the initial probability vector x0 is, the sequence of probability vectors {xk} will converge to q as k tends to infinity.
Example 3.2.2: The Steady State Vector for the Cat and Box Game
To find a steady state vector q, we know from the theorem above that it must be the unique probability vector that solves the equation Px=x. So let’s use linear algebra to find Cosmo’s steady state vector!
Video 3.2.2: Finding the Steady State Vector
Your browser does not support the video tag.
Video by Dr. Lisa Korf: Finding the Steady State Vector [© CCCOnline]
Video 3.2.2 Transcript: Finding the Steady State Vector
We can now see that in the long run, Cosmo will pick box 1 with probability 40/481, or about 8.316% of the time, box 2 with probability 49/481, or about 10.187% of the time, and box 3 with probability 392/481, or about 81.497% of the time. Based on these results, do you think that in reality, Cosmo’s memory would really only go back one time period, or would she be more likely to avoid box 1, perhaps indicating a longer-term memory?
Image Attributions: Dog Silhouette: https://pixabay.com/en/dog-puppy-animal-black-silhouette-159481/ , Cat Silhouette: http://www.clipartbest.com/clipart-ndT88AKie , Mouse Silhouette: (mirrored image) http://www.rgbstock.com/bigphoto/mSdQGMK/silhouette+gerbil
check your understanding
Crystal Lattices:
Suppose that a crystal lattice has a lattice basis given by LB = { u, v, w }, and that an atom is located at x = u + 3v + w. What are the coordinates of the atom with respect to the lattice basis LB?
See Answer
Answer: xLB=131{“version”:”1.1″,”math”:”xLB=131″}
A crystal lattice has a particular atom type appearing at coordinates (1, 2, –3), and repeating in three different directions at coordinates (–1, 2, 4), (2, 3, 5), and (4, –2, 1). What is a possible lattice basis for this crystal lattice?
See Answer
Answer: -1-12-24-(-3), 2-13-25-(-3), 4-1-2-21-(-3)=-207, 118, 3-44{“version”:”1.1″,”math”:”-1-12-24-(-3), 2-13-25-(-3), 4-1-2-21-(-3)=-207, 118, 3-44″}
Markov Chains
:
Suppose that Cosmo’s initial probability vector had for some reason been 001{“version”:”1.1″,”math”:”001″}, with all other probabilities the same. How would this affect Cosmo’s steady state vector?
See Answer
Answer: It wouldn’t, since when a steady state vector exists, which it does in this case, it is independent of the initial probability vector state.
For a Markov chain with stochastic matrix equal to the 2 by 2 matrix, P=0110{“version”:”1.1″,”math”:”P=0110″}, does a steady state vector exist? If so, what is it? If not, why not?
See Answer
Answer: Although this matrix does not satisfy the assumptions of the theorem, such a matrix may still have a steady state if it has a unique probability vector solution to the equation P-Ix=0{“version”:”1.1″,”math”:”P-Ix=0″}. In this case, P-I=-111-1{“version”:”1.1″,”math”:”P-I=-111-1″}, which has reduced row-echelon form, 1-100{“version”:”1.1″,”math”:”1-100″}. The solution to the homogeneous system is x=x211{“version”:”1.1″,”math”:”x=x211″}. Letting x2=1/2{“version”:”1.1″,”math”:”x2=1/2″} yields the unique probability vector solution, 1/21/2{“version”:”1.1″,”math”:”1/21/2″}, which is also the steady state vector.
Discuss!
Choose one of the applications, and one of the numbered Discussion Topics for that application to explore with your classmates in the Exploration Discussion in this Module.
Crystal Lattices – Discussion Topics
Using the lattice basis LB and the reciprocal lattice basis RLB given for graphite in Examples 3.1.1 and 3.1.2, you wish to determine the location of a point of interest within the lattice given by vRLB=abc{“version”:”1.1″,”math”:”vRLB=abc”}, where a, b,{“version”:”1.1″,”math”:”a, b,”}and c{“version”:”1.1″,”math”:”c”} are integers, not all zero, to be chosen by you and such that the vector is different from that chosen in any other posts on this topic . Give the coordinates of this point both with respect to the standard basis vectors for R3, and with respect to the lattice basis LB. Explain your steps carefully. Which basis of the three seems most natural for this point of interest? Be sure to include your chosen values for a, b,{“version”:”1.1″,”math”:”a, b,”}and c{“version”:”1.1″,”math”:”c”} in the Subject Heading of your post.
Diamond is a crystalline structure with a cubic unit cell. The lattice basis for diamond is given by LB=a2011, a2101, a2110, where a is a lattice constant. The reciprocal lattice basis for diamond is given by RLB=2πa-111, 2πa1-11, 2πa11-1. Using the appropriate 3 by 6 augmented matrices and row reduction techniques, find the change of coordinates matrices, PLB←RLB, PR3←RLB, and PLB←R3, {“version”:”1.1″,”math”:”PLB←RLB, PR3←RLB, and PLB←R3, ”}and show that one of these matrices may be obtained as a product of the other two.
Markov Chains – Discussion Topics
As indicated at the end of the Markov Chain discussion, realistically it would be more likely that Cosmo the cat would remember more than just one time period back. Assume that Cosmo recalls the last two time periods. To set up a Markov Chain in that setting, one needs to set up new states (probability vectors) that take the last two time periods from the original setting into account. So the state might be given by x=x1,1×1,2×1,3×2,1×2,2×2,3×3,1×3,2×3,3,where x(i, j) is the component of the probability vector that gives the probability that Cosmo will choose box i in this time period, and box j in the following time period. To set up a Markov chain, one possible stochastic matrix P is given by From: P= (1,1) (1,2) (1,3) (2,1) (2,2) (2,3) (3,1) (3,2) (3,3) 0.02000.21000.035000.49000.345000.27000.49000.445000.6950000.21000.4000.225000.345000.2000.125000.445000.4000.650000.035000.225000.05000.27000.125000.05000.695000.65000.9To:(1,1)(1,2)(1,3)(2,1)(2,2)(2,3)(3,1)(3,2)(3,3) where the (a,b), (b,c) entry corresponds to the probability that Cosmo will choose box c on the third try if she chooses box a on the first try and box b on the second try.Why do some of the entries have to be 0 in this new setting? Set up your own stochastic matrix for Cosmo with a two period memory like this one. The only entries that should be different are in the (1,1) column, the (2,2) column, and the (3,3) column. Get inside Cosmo’s head and decide what her probability is of choosing box 1 after having chosen it twice already, and similarly for boxes 2 and 3. This will give you the values for the (1,1), (1,1) entry, (2,2), (2,2) entry, and the (3,3), (3,3) entry. The remaining two nonzero values in each of those columns should be equal to each other, and the sum of the three nonzero entries in each column should sum to 1, since it is a stochastic matrix. For example, you might choose the value 0.01 (instead of 0.02) for the (1,1), (1,1) entry. In that case, the remaining two nonzero entries would both be equal to 0.99/2 = 0.495. For your stochastic matrix setup, what would Cosmo’s steady state vector q be? What is the interpretation of each of the 9 entries in q? What is Cosmo’s total probability of choosing box 1 in the long run, and how did you obtain it? Is this more realistic? Be sure to include your chosen values for the (1,1), (1,1) , (2,2), (2,2) , and (3,3), (3,3) entries in the Subject Heading of your main post, since they will need to be different from the values chosen in any other posts on this topic.
References
Crystal Lattices:
Chung, D. D. (2002). Review Graphite. Journal of Materials Science, 37(8), 1475-1489.
Averill, B., & Eldredge, P. (2007). Chemistry: Principles, patterns, and applications. San Francisco, CA: Pearson Benjamin Cummings.
Ungersböck, S.E. (1977, August 8). 3.4 Basic Properties of the Diamond Structure, Dissertation: Advanced Modeling of Strained CMOS Technology. Retrieved April 18, 2016, from
http://www.iue.tuwien.ac.at/phd/ungersboeck/node27.html
Markov Chains:
Khoury, J. (2003, October 28). Application to Markov Chains. Retrieved April 09, 2016, from
http://aix1.uottawa.ca/~jkhoury/markov.htm
Linear Algebra/Topic: Markov Chains. (2010, May 26). Retrieved April 09, 2016, from
https://en.wikibooks.org/wiki/Linear_Algebra/Topic:_Markov_Chains
Markov Chains, Math 22, Spring 2007. (2007). Retrieved April 9, 2016, from
https://math.dartmouth.edu/archive/m22s07/public_html/MarkovSlides
Introduction
Crystal Lattices
Markov Chains
Check Understanding
Discuss!
References
Module 03/M3_Exploration.html
Exploration of Vector Spaces and Subspaces
Introduction
Crystal Lattices
Markov Chains
Check Understanding
Discuss!
References
Vector spaces and subspaces
Outcomes
Recognize the general vector space framework with basic definitions and problems. (D)
Demonstrate how closely other vector spaces resemble Rn. (D)
Use vector space terminology to tie together important facts about rectangular matrices. (A, B, D)
Apply theory to Markov Chains models in biology, business, and engineering. (A, B, C, D, I)
In these Explorations, we will explore a few introductory applications of vector spaces and subspaces. Click on the tabs to explore each application, and read and watch all of the examples. Then, answer the questions on the Check Understanding tab.
The Discuss! tab contains a few discussion questions about the applications in this module, for use in the Exploration Discussion. The Exploration Discussion will provide you with an opportunity to expand further upon the ideas presented here.
Crystal Lattices
Note: If equations don’t render correctly when initially viewed, refresh your browser page and wait for it to load completely before visiting any of the tabs.
Some applications have a physically natural coordinate system that may be different from the standard coordinate system with standard unit basis vectors in Rn{“version”:”1.1″,”math”:”Rn”}.
In crystallography, a crystal lattice in R3{“version”:”1.1″,”math”:”R3″} is composed of unit cells, which are minimal parallelepipeds, that when translated by each of the lattice vectors along three adjacent edges, make up the crystal lattice, as shown.
Fig. 3.1.1 (generated from): Averill, B., & Eldredge, P. (2007). Chemistry: Principles, patterns, and applications. San Francisco, CA: Pearson Benjamin Cummings.
Figure 3.1.1: A three-dimensional unit cell with its three lattice vectors, and the resulting crystal lattice.
The arrangement of atoms within a unit cell is repeated throughout the crystal lattice. Therefore, knowledge about a crystal lattice can be reduced to knowledge about the unit cell and the arrangement of atoms within it. There are 14 basic types of unit cells comprising all known crystal lattice structures.
Scientists use the lattice vectors that make up the adjacent edges of a unit cell as their working basis of R3{“version”:”1.1″,”math”:”R3″}, forming a natural coordinate system. This gives the atomic arrangement within each unit cell its natural coordinates, whereas using the standard unit basis and resulting unit cubes might be useful for some calculations, but unintuitive for others, since they do not fit the crystal’s natural structure.
By using this natural coordinate system, scientists can combine atom configurations with formulas derived from the relevant type of unit cell to determine such things about a given material as: in what direction it is likely to fracture, where it is easiest to make a clean cut, in what direction it has the greatest tensile strength, or the greatest compressive resistance, and so on. These directions are naturally aligned with the directions of the lattice vectors, so using the lattice vectors as a basis is not only natural, but also a very useful choice.
Example 3.1.1: Finding a Coordinate Vector of an Atom Relative to a Lattice Basis for Graphite
A lattice basis for graphite is given by LB=a3/2-1/20, a3/21/20, c001{“version”:”1.1″,”math”:”LB=a3/2-1/20, a3/21/20, c001″}, where a=2.46 Ao{“version”:”1.1″,”math”:”a=2.46 Ao”} (angstroms), and c=6.71 Ao{“version”:”1.1″,”math”:”c=6.71 Ao”}.
An atom within a unit cell is located at v=-a/23-a/2c/2{“version”:”1.1″,”math”:”v=-a/23-a/2c/2″}. The following video shows how to find the coordinates of this atom relative to the lattice basis.
Video 3.1.1: Finding the Coordinates of an Atom Relative to a Lattice Basis
Your browser does not support the video tag.
Video by Dr. Lisa Korf: Finding the Coordinates of an Atom Relative to a Lattice Basis [© CCCOnline]
Video 3.1.1 Transcript: Finding the Coordinates of an Atom Relative to a Lattice Basis
It is also important to be able to convert from a lattice basis back to the standard unit basis or to another basis that may be useful for other computations. One such basis is called the reciprocal lattice basis.
Example 3.1.2: Converting from the Lattice Basis for Graphite to the Reciprocal Lattice basis
We begin again with the lattice basis for graphite given by LB=a3/2-1/20, a3/21/20, c001{“version”:”1.1″,”math”:”LB=a3/2-1/20, a3/21/20, c001″}.
This time, we are given a coordinate vector with respect to the lattice basis, vLB=-201{“version”:”1.1″,”math”:”vLB=-201″}, and we want to know the coordinates of v{“version”:”1.1″,”math”:”v”}
with respect to the reciprocal lattice basis given by RLB=2πa1/3-10, 2πa1/310, 2πc001.
Video 3.1.2: Converting Coordinates Between Two Lattice Bases
Your browser does not support the video tag.
Video by Dr. Lisa Korf: Converting the Coordinates Between Two Lattice Basis [© CCCOnline]
Video 3.1.2 Transcript: Converting Coordinates Between Two Lattice Bases
Crystallography has many areas of application, such as to metallurgy, solid state electronics, and geology. In all of these fields, it is important to be able to perform change of basis operations from one coordinate system to another.
Markov Chains
Consider a system that can be in exactly one of n possible states. Suppose that the system evolves in time such that the probability of being in a particular state at a given time period k depends only on the system’s state at the previous time period, k-1{“version”:”1.1″,”math”:”k-1″}.
For example, suppose you present your pet cat, Cosmo, with a Box Game, in which there are three boxes that she can open. Inside box 1 is a cat-chasing Chihuahua, while box 2 is empty, and box 3 contains an enticing mouse. If Cosmo opens box 1 and finds the Chihuahua, there is only a 2% probability that she will open that box the next time. If she opens box 2, finding it empty, there is a 20% chance she will open that box the next time, and if she opens box 3, there is a 90% chance she will open that box for a repeat game of cat and mouse. We’ll assume that Cosmo only remembers what happened the previous time she opened a box, but not any of the times prior to that.
Such a scenario can be modeled with a Markov Chain.
Definition 3.2.1. A probability vector is a vector of states whose entries sum to 1. A stochastic matrix is a square matrix with probability vectors as its columns, i.e. the entries of every column sum to one. A Markov Chain is a sequence of probability vectors x0, x1, …, xk, xk+1, …{“version”:”1.1″,”math”:”x0, x1, …, xk, xk+1, …”}, together with a stochastic matrix P, such thatx1=Px0, x2=Px1, …, xk+1=Pxk, ….The probability vectors represent a sequence of states, and the stochastic matrix P has entries Pij that give the transition probabilities of going from state j to state i.
Example 3.2.1: The Stochastic Matrix and Probability Vectors for the Cat and Box Game
If Cosmo the cat picks box 1 on her first try, then the probability that she will pick box 1 on her second try is only 2%, so the stochastic matrix P for this example would have entry P11=0.02. Assuming that the remaining probability of 0.98 is split evenly between the remaining boxes in that column, and similarly for the remaining columns, we can write down the stochastic matrix for Cosmo: From Box: P= 1 2 3 0.020.40.050.490.20.050.490.40.9To Box:123If Cosmo has a 1/3 probability of opening each of the three boxes at time period 0, her initial state vector would be x0=1/31/31/3{“version”:”1.1″,”math”:”x0=1/31/31/3″}. The following video computes the first 6 state vectors that result after Cosmo chooses which box to open at 6 consecutive time periods.
Video 3.2.1: Computing the First 6 State Vectors of a Markov Chain
Your browser does not support the video tag.
Video by Dr. Lisa Korf: Computing the First 6 State Vectors of a Markov Chain [© CCCOnline]
Video 3.2.1 Transcript: Computing the First 6 State Vectors of a Markov Chain
From the video, it appears that as Cosmo continues to choose boxes, her state vector starts to hone in on a single vector. This is not a coincidence, due to the following theorem.
Definition 3.2.2. A stochastic matrix P is regular if there is any power of P whose entries are all strictly positive.
Theorem 3.2.1. If P is a regular stochastic matrix, then there is a unique probability vector, q, such that Pq = q. The vector q is called the steady state vector for P. In addition, no matter what the initial probability vector x0 is, the sequence of probability vectors {xk} will converge to q as k tends to infinity.
Example 3.2.2: The Steady State Vector for the Cat and Box Game
To find a steady state vector q, we know from the theorem above that it must be the unique probability vector that solves the equation Px=x. So let’s use linear algebra to find Cosmo’s steady state vector!
Video 3.2.2: Finding the Steady State Vector
Your browser does not support the video tag.
Video by Dr. Lisa Korf: Finding the Steady State Vector [© CCCOnline]
Video 3.2.2 Transcript: Finding the Steady State Vector
We can now see that in the long run, Cosmo will pick box 1 with probability 40/481, or about 8.316% of the time, box 2 with probability 49/481, or about 10.187% of the time, and box 3 with probability 392/481, or about 81.497% of the time. Based on these results, do you think that in reality, Cosmo’s memory would really only go back one time period, or would she be more likely to avoid box 1, perhaps indicating a longer-term memory?
Image Attributions: Dog Silhouette: https://pixabay.com/en/dog-puppy-animal-black-silhouette-159481/ , Cat Silhouette: http://www.clipartbest.com/clipart-ndT88AKie , Mouse Silhouette: (mirrored image) http://www.rgbstock.com/bigphoto/mSdQGMK/silhouette+gerbil
check your understanding
Crystal Lattices:
Suppose that a crystal lattice has a lattice basis given by LB = { u, v, w }, and that an atom is located at x = u + 3v + w. What are the coordinates of the atom with respect to the lattice basis LB?
See Answer
Answer: xLB=131{“version”:”1.1″,”math”:”xLB=131″}
A crystal lattice has a particular atom type appearing at coordinates (1, 2, –3), and repeating in three different directions at coordinates (–1, 2, 4), (2, 3, 5), and (4, –2, 1). What is a possible lattice basis for this crystal lattice?
See Answer
Answer: -1-12-24-(-3), 2-13-25-(-3), 4-1-2-21-(-3)=-207, 118, 3-44{“version”:”1.1″,”math”:”-1-12-24-(-3), 2-13-25-(-3), 4-1-2-21-(-3)=-207, 118, 3-44″}
Markov Chains
:
Suppose that Cosmo’s initial probability vector had for some reason been 001{“version”:”1.1″,”math”:”001″}, with all other probabilities the same. How would this affect Cosmo’s steady state vector?
See Answer
Answer: It wouldn’t, since when a steady state vector exists, which it does in this case, it is independent of the initial probability vector state.
For a Markov chain with stochastic matrix equal to the 2 by 2 matrix, P=0110{“version”:”1.1″,”math”:”P=0110″}, does a steady state vector exist? If so, what is it? If not, why not?
See Answer
Answer: Although this matrix does not satisfy the assumptions of the theorem, such a matrix may still have a steady state if it has a unique probability vector solution to the equation P-Ix=0{“version”:”1.1″,”math”:”P-Ix=0″}. In this case, P-I=-111-1{“version”:”1.1″,”math”:”P-I=-111-1″}, which has reduced row-echelon form, 1-100{“version”:”1.1″,”math”:”1-100″}. The solution to the homogeneous system is x=x211{“version”:”1.1″,”math”:”x=x211″}. Letting x2=1/2{“version”:”1.1″,”math”:”x2=1/2″} yields the unique probability vector solution, 1/21/2{“version”:”1.1″,”math”:”1/21/2″}, which is also the steady state vector.
Discuss!
Choose one of the applications, and one of the numbered Discussion Topics for that application to explore with your classmates in the Exploration Discussion in this Module.
Crystal Lattices – Discussion Topics
Using the lattice basis LB and the reciprocal lattice basis RLB given for graphite in Examples 3.1.1 and 3.1.2, you wish to determine the location of a point of interest within the lattice given by vRLB=abc{“version”:”1.1″,”math”:”vRLB=abc”}, where a, b,{“version”:”1.1″,”math”:”a, b,”}and c{“version”:”1.1″,”math”:”c”} are integers, not all zero, to be chosen by you and such that the vector is different from that chosen in any other posts on this topic . Give the coordinates of this point both with respect to the standard basis vectors for R3, and with respect to the lattice basis LB. Explain your steps carefully. Which basis of the three seems most natural for this point of interest? Be sure to include your chosen values for a, b,{“version”:”1.1″,”math”:”a, b,”}and c{“version”:”1.1″,”math”:”c”} in the Subject Heading of your post.
Diamond is a crystalline structure with a cubic unit cell. The lattice basis for diamond is given by LB=a2011, a2101, a2110, where a is a lattice constant. The reciprocal lattice basis for diamond is given by RLB=2πa-111, 2πa1-11, 2πa11-1. Using the appropriate 3 by 6 augmented matrices and row reduction techniques, find the change of coordinates matrices, PLB←RLB, PR3←RLB, and PLB←R3, {“version”:”1.1″,”math”:”PLB←RLB, PR3←RLB, and PLB←R3, ”}and show that one of these matrices may be obtained as a product of the other two.
Markov Chains – Discussion Topics
As indicated at the end of the Markov Chain discussion, realistically it would be more likely that Cosmo the cat would remember more than just one time period back. Assume that Cosmo recalls the last two time periods. To set up a Markov Chain in that setting, one needs to set up new states (probability vectors) that take the last two time periods from the original setting into account. So the state might be given by x=x1,1×1,2×1,3×2,1×2,2×2,3×3,1×3,2×3,3,where x(i, j) is the component of the probability vector that gives the probability that Cosmo will choose box i in this time period, and box j in the following time period. To set up a Markov chain, one possible stochastic matrix P is given by From: P= (1,1) (1,2) (1,3) (2,1) (2,2) (2,3) (3,1) (3,2) (3,3) 0.02000.21000.035000.49000.345000.27000.49000.445000.6950000.21000.4000.225000.345000.2000.125000.445000.4000.650000.035000.225000.05000.27000.125000.05000.695000.65000.9To:(1,1)(1,2)(1,3)(2,1)(2,2)(2,3)(3,1)(3,2)(3,3) where the (a,b), (b,c) entry corresponds to the probability that Cosmo will choose box c on the third try if she chooses box a on the first try and box b on the second try.Why do some of the entries have to be 0 in this new setting? Set up your own stochastic matrix for Cosmo with a two period memory like this one. The only entries that should be different are in the (1,1) column, the (2,2) column, and the (3,3) column. Get inside Cosmo’s head and decide what her probability is of choosing box 1 after having chosen it twice already, and similarly for boxes 2 and 3. This will give you the values for the (1,1), (1,1) entry, (2,2), (2,2) entry, and the (3,3), (3,3) entry. The remaining two nonzero values in each of those columns should be equal to each other, and the sum of the three nonzero entries in each column should sum to 1, since it is a stochastic matrix. For example, you might choose the value 0.01 (instead of 0.02) for the (1,1), (1,1) entry. In that case, the remaining two nonzero entries would both be equal to 0.99/2 = 0.495. For your stochastic matrix setup, what would Cosmo’s steady state vector q be? What is the interpretation of each of the 9 entries in q? What is Cosmo’s total probability of choosing box 1 in the long run, and how did you obtain it? Is this more realistic? Be sure to include your chosen values for the (1,1), (1,1) , (2,2), (2,2) , and (3,3), (3,3) entries in the Subject Heading of your main post, since they will need to be different from the values chosen in any other posts on this topic.
References
Crystal Lattices:
Chung, D. D. (2002). Review Graphite. Journal of Materials Science, 37(8), 1475-1489.
Averill, B., & Eldredge, P. (2007). Chemistry: Principles, patterns, and applications. San Francisco, CA: Pearson Benjamin Cummings.
Ungersböck, S.E. (1977, August 8). 3.4 Basic Properties of the Diamond Structure, Dissertation: Advanced Modeling of Strained CMOS Technology. Retrieved April 18, 2016, from
http://www.iue.tuwien.ac.at/phd/ungersboeck/node27.html
Markov Chains:
Khoury, J. (2003, October 28). Application to Markov Chains. Retrieved April 09, 2016, from
http://aix1.uottawa.ca/~jkhoury/markov.htm
Linear Algebra/Topic: Markov Chains. (2010, May 26). Retrieved April 09, 2016, from
https://en.wikibooks.org/wiki/Linear_Algebra/Topic:_Markov_Chains
Markov Chains, Math 22, Spring 2007. (2007). Retrieved April 9, 2016, from
https://math.dartmouth.edu/archive/m22s07/public_html/MarkovSlides
Introduction
Crystal Lattices
Markov Chains
Check Understanding
Discuss!
References
Place an order in 3 easy steps. Takes less than 5 mins.