"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"# The Schrödinger Equation\n",
"\n",
">The underlying physical laws necessary for the mathematical theory of a large part of physics and the whole of chemistry are thus completely known, and the difficulty is only that the exact application of these laws leads to equations much too complicated to be soluble. It therefore becomes desirable that approximate practical methods of applying quantum mechanics should be developed, which can lead to an explanation of the main features of complex atomic systems without too much computation.\n",
" - Paul A. M. Dirac (1929)\n",
"\n",
"\n",
"## 🥅 Learning Objectives\n",
"- Review of Plane Wave(functions)\n",
"- \"Derivation\" of the Time-Dependent Schrödinger Equation (TDSE)\n",
"- The Hamiltonian\n",
"- The Time-Independent Schrödinger Equation (TISE)\n",
"- Separation of Variables Method"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"## Motivation for the Schrödinger Equation\n",
"Just like Newton’s equations or Maxwell’s equations, the Schrödinger equation is a law of nature: it cannot be “derived.” Instead, it must be postulated, and then tested by experiment. However, just as Newton’s equations can be justified insofar as they encapsulate the “laws of motion” formulated by Galileo and Newton, and Maxwell’s equations can be justified insofar as they encapsulate the “laws of electromagnetism” formulated by Gauss, Faraday, Coulomb, etc., Schrödinger's equation can be justified by observing that it encapsulates the “laws” of Planck and De Broglie, namely that\n",
"$$ E = h \\nu = \\hbar \\omega $$\n",
"$$ p = \\tfrac{h}{\\lambda} = \\tfrac{h \\nu}{c} = \\hbar k = h \\tilde{\\nu} $$\n",
"\n",
"### Angular Frequency and Wavenumber\n",
"In these equations I have introduced several new symbols, mostly related to the fact it is often convenient to use angular frequency, \n",
"$$ \\omega = 2 \\pi \\nu $$\n",
"which naturally pairs with \n",
"$$ \\hbar = \\tfrac{h}{2 \\pi} $$\n",
"Similarly, it is often convenient, especially in spectroscopic studies, to use wavenumber,\n",
"$$ \\tilde{\\nu} = \\tfrac{1}{\\lambda} = \\lambda^{-1} $$\n",
"or its angular analogue\n",
"$$ k = \\tfrac{2 \\pi}{\\lambda} = 2 \\pi \\tilde{nu} $$\n",
"Sometimes it is also useful to consider the period of the wave, \n",
"$$ \\text{T} = \\nu^{-1} $$\n",
"Recall that the wavelength, frequency, period, and wavenumber of a wave can be related to the speed of light by relations like:\n",
"$$ c = \\lambda \\nu = \\tfrac{\\lambda}{\\text{T}} = \\tilde{\\nu} \\nu $$\n",
"$$ c = \\tfrac{\\omega}{k} $$"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"### The Wavefunction\n",
"![1D Plane Wave](https://github.com/PaulWAyers/IntroQChem/blob/main/linkedFiles/planewave1D.gif?raw=true? \"Visualization of a the real part of a 1-dimensional travelling plane wave, from Szilágyi, András (2019): EMANIM: Interactive visualization of electromagnetic waves. Web application available at URL https://emanim.szialab.org\")\n",
"\n",
"The use of angular frequency and wave number is especially useful when dealing with waves, because the need to explicitly include factors of $2 \\pi$ to describe the periodic oscillations is removed. A one-dimensional plane wave, moving in the $x$ direction, with amplitude $A$, wavelength $\\lambda = \\tfrac{2 \\pi}{k}$, and frequency $\\nu = \\tfrac{\\omega}{2 \\pi}$ is described by\n",
"$$ \\Psi(x,t) = A e^{i (kx - \\omega t)} = A\\left(\\cos(kx-\\omega t) + i \\sin(kx - \\omega t) \\right) $$\n",
"Since $\\Psi(x,t)$ is the function that describes a wave, we call it a \"wavefunction.\" \n",
"\n",
"### Expression for the Total Energy\n",
"Based on electromagnetism, we believe this equation should give us the essential insight we need to describe the \"waviness of particles.\" We need, however, a second equation that will allow us to describe the \"particulateness of waves.\" Since the energy of a wave is proportional to $\\omega$ (Planck-Einstein relation) and the momentum of a wave is proportional to $k$ (De Broglie relation), an equation that relates the energy to the momentum is an obvious candidate. For example,\n",
"\\begin{align}\n",
"\\text{energy} &= \\text{kinetic energy} + \\text{potential energy} \\\\\n",
"E &= T + V(x,t) \\\\\n",
"E &= \\frac{p^2}{2m} + V(x,t) \n",
"\\end{align}\n"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"### The Time-Dependent Schrödinger Equation (TDSE)\n",
"We need to find a way to combine the key equations from the previous two sections. To do this, note that\n",
"$$ \\frac{\\partial \\Psi(x,t)}{\\partial t} = (-i \\omega) A e^{i (kx - \\omega t)} = (-i \\omega) \\Psi(x,t) $$\n",
"We can link this to Planck's equation for the energy of a photon by multiplying both sides by $i \\hbar$, \n",
"\\begin{align}\n",
"i \\hbar \\frac{\\partial \\Psi(x,t)}{\\partial t} &= (i \\hbar) (-i \\omega) \\Psi(x,t) = \\hbar \\omega \\Psi(x,t) \\\\\n",
"&= E \\Psi(x,t) \n",
"\\end{align}\n",
"Similarly, differentiating with respect to position,\n",
"\\begin{align}\n",
"\\frac{\\partial \\Psi(x,t)}{\\partial x} &= (i k) A e^{i (kx - \\omega t)} \\\\\n",
"\\frac{\\partial^2 \\Psi(x,t)}{\\partial x^2} &= (i k)^2 A e^{i (kx - \\omega t)} = -k^2 \\Psi(x,t) \n",
"\\end{align}\n",
"Multiplying both sides by $-\\hbar^2$ provides a link to the De Broglie relation,\n",
"$$\n",
"-\\hbar^2 \\frac{\\partial^2 \\Psi(x,t)}{\\partial x^2} = (\\hbar k)^2 \\Psi(x,t) = p^2 \\Psi(x,t)\n",
"$$\n",
"\n",
"Now, let's consider our expression for the energy,\n",
"\\begin{align}\n",
"E &= \\frac{p^2}{2m} + V(x,t) \\\\\n",
"E \\Psi(x,t) &= \\frac{p^2}{2m} \\Psi(x,t) + V(x,t)\\Psi(x,t)\n",
"\\end{align}\n",
"Substituting in the equations for the energy and the momentum-squared in terms of wavefunction derivatives, we have:\n",
"$$ i \\hbar \\frac{\\partial \\Psi(x,t)}{\\partial t} = -\\frac{\\hbar^2}{2m}\\frac{\\partial^2 \\Psi(x,t)}{\\partial x^2} + V(x,t)\\Psi(x,t) $$\n",
"This is the [*time-dependent Schrödinger equation*](https://en.wikipedia.org/wiki/Schr%C3%B6dinger_equation#Time-dependent_equation), which we often abbreviate as the TDSE.\n",
"\n",
"The biggest limitations of the Schrödinger equation is that it was derived using the nonrelativistic expression for the kinetic energy of a particle, $T = \\tfrac{p^2}{2m}$. For this reason, the Schrödinger equation is not valid for particles that are moving at a substantial fraction of the speed of light (where relativistic effects are important) or for light itself. Quantum mechanics can be extended to such cases, and relativistic quantum mechanics is quite important for the heavier atoms in the periodic table. (Typically relativity is neglected for atoms up to Zinc (Z=30) or Krypton (Z=36), and the included past that. However there are a few molecules containing lighter elements—e.g., the Sulfur dimer—where there are chemically (or at least spectroscopically) interesting relativistic effects.) The full relativistic treatment of quantum mechanics leads to “quantum electrodynamics,” which is a rich and interesting subject that we, alas, will not have time to cover. "
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"### The Hamiltonian\n",
"An equation like Schrödinger’s may seem mysterious to us now, but at the time it was derived it seemed very natural. A top-rate physicist like Schrödinger had a rich background in classical physics, and he would have understood that there were many results, from classical physics, that suggested that the equations for radiation and for matter could be combined. One of the strongest such suggestions came from two principles: “Fermat’s principle of least time” and the “Hamilton-Lagrange principle of least action.” Fermat’s principle says that a ray of light goes from point A to point B in the fastest possible way, and can be used to derive all of classical optics. The Hamilton-Lagrange principle says that a particle goes from point A to point B in a way that makes minimizes its action (the product of its momentum and its change in position) along the path. Schrödinger would have noticed that Planck’s constant had units of action, and then quickly ascertained that laws of quantum mechanics might be governed by an equation similar to the Hamilton-Jacobi equation for the action, \n",
"$$\\frac{1}{2m}\\left(\\frac{\\partial S}{\\partial x}\\right)^2 + V(x,t) = \\frac{\\partial S}{\\partial t} $$\n",
"or, equivalently,\n",
"$$ S \\left(-\\frac{1}{2m} \\frac{\\partial^2S}{\\partial x^2} + V(x,t) \\right) = \\frac{\\partial S}{\\partial t} - \\frac{1}{2m}\\frac{\\partial^2S^2}{\\partial x^2}$$\n",
"The second term on the right-hand-side of this equation ordinarily integrates to zero (owing to the divergence theorem).\n",
"\n",
"Similarly, as hinted at in the previous set of notes, the Schrödinger equation can be deduced from the (quantum-mechanical) principle of least (really stationary) action. Indeed, the stationary-action principle is more fundamental than the Schrödinger equation, and is the basis for quantum electrodynamics, which is the relativistic theory for matter and electromagnetic radiation that undergirds classical mechanics and classical electromagnetism (in the $\\hbar \\rightarrow 0$ limit) and the Schrödinger equation (for nonrelativistic matter, without concern for how matter and radiation interact). Feynman and Schwinger won the Nobel Prize for quantum electrodynamics, but its origins go back to Schrödinger's and Dirac's observations about the links between the (least) action principle and the Schrödinger equation.\n",
"\n",
"The Schrödinger equation is usually written in the form:\n",
"$$ \\color{blue}{ \\left( -\\frac{\\hbar^2}{2m} \\frac{\\partial^2}{\\partial x^2}+V(x,t)\\right)}\\Psi(x,t) = i \\hbar \\frac{\\partial \\Psi(x,t)}{\\partial t} $$\n",
"The term in blue is the [Hamiltonian](https://en.wikipedia.org/wiki/Hamiltonian_(quantum_mechanics)) operator, $\\hat{H}(x,t)$; it is the operator for the energy in quantum mechanics. In classical mechanics, the [Hamiltonian](https://en.wikipedia.org/wiki/Hamiltonian_mechanics) is exactly the equation that we started with when we wrote a (classical) expression for the energy of matter, $H(x,t) = \\tfrac{p^2}{2m} + V(x,t)$. \n",
"\n",
"Using the Hamiltonian, we can then write the TDSE in a compact form, \n",
"$$ \\color{blue}{\\hat{H}(x,t)}\\Psi(x,t) = i \\hbar \\frac{\\partial \\Psi(x,t)}{\\partial t} $$\n",
"As in classical mechanics, the Hamiltonian is the operator for the energy. In quantum mechanics, every observable quantity corresponds to an operator. For example, the kinetic energy operator is\n",
"$$ \\hat{T} = -\\frac{\\hbar^2}{2m} \\frac{\\partial^2}{\\partial x^2} $$\n",
"Moreover, the values of the properties of any quantum-mechanical operator correspond to its eigenvalues,\n",
"$$ (\\text{operator}) \\Psi = (\\text{eigenvalue}) \\Psi $$\n",
"This suggests that there should be another equation with the form\n",
"$$ \\hat{H}\\Psi = E \\Psi $$\n",
"This equation is called the [time-independent Schrödinger equation](https://en.wikipedia.org/wiki/Schr%C3%B6dinger_equation#Time-independent_equation) (TISE).\n",
"\n",
"My father, an organic chemist, always says he knew that quantum mechanics would be problematic the first time he saw the Schrödinger equation. He still refers to $\\Psi$ as \"Satan's pitchfork\" and $\\hat{H}$ as the \"hell operator.\" Most people, however, call $\\hat{H}$ the Hamiltonian and $\\Psi$ the wavefunction."
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"### The Time-Independent Schrödinger Equation (TISE)\n",
"In much of chemistry, the Hamiltonian operator is not time-dependent. For examples, the electrons in a molecule feel a fixed (time-independent) attraction to the nuclei in the molecule and, to a first approximation, the (very small) motion of the nuclei can be neglected. Similarly, in the absence of a time-varying external field, a molecule as a whole (containing both electrons and nuclei) does not feel a time-varying potential. \n",
"\n",
"To derive the time-independent Schrödinger equation, we insert the time-independent Hamiltonian, \n",
"$$\\hat{H}(x) = -\\frac{\\hbar^2}{2m} \\frac{\\partial^2}{\\partial x^2} +V(x)$$\n",
"into the time-dependent Schrödinger equation,\n",
"$$ \\hat{H}(x)\\Psi(x,t) = i \\hbar \\frac{\\partial \\Psi(x,t)}{\\partial t} $$\n",
"and then rearrange the equation into the form\n",
"$$ \\left( \\hat{H}(x) - i\\hbar \\frac{\\partial}{\\partial t} \\right) \\Psi(x,t) = 0 $$\n",
"This is a differential equation that can be solved by the general technique of [separation of variables](https://en.wikipedia.org/wiki/Separation_of_variables). Specifically,\n",
"\n",
"> [**Separation of Variables**](https://en.wikipedia.org/wiki/Separation_of_variables#Partial_differential_equations) Given an operator that is a sum of two operators that depend on different (sets of) variables,\n",
"$$ \\hat{A}(x,y) = \\hat{A}_x(x) + \\hat{A}_y(y) $$\n",
"The solution to the eigenvalue equation\n",
"$$ \\hat{A}(x,y) \\Psi(x,y) = a \\Psi(x,y) $$\n",
"is \n",
"$$\\Psi(x,y) = \\Psi_x(x) + \\Psi_y(y) $$\n",
"$$ a = a_x + a_y $$\n",
"where \n",
"$$\\hat{A}_x(x) \\Psi_x(x) = a_x \\Psi_x(x) $$\n",
"$$\\hat{A}_y(y) \\Psi_y(y) = a_y \\Psi_y(y) $$\n",
"\n",
"\n",
"This tells us that for a time-independent Hamiltonian, the wavefunction has the form $\\Psi(x,t) = \\psi(x) \\phi(t)$, where\n",
"$$\\hat{H}(x) \\psi(x) = E \\psi(x) $$\n",
"$$- i\\hbar \\frac{\\partial}{\\partial t} \\phi(t) = -E\\phi(t) $$\n",
"In the first equation, we have chosen the eigenvalue of the TDSE to be the energy since we know that the Hamiltonian is the quantum-mechanical operator for the energy. The second equation is a separable differential equation, which can be solved as follows:\n",
"\\begin{align}\n",
"-i\\hbar \\frac{\\partial}{\\partial t} \\phi(t) &= -E\\phi(t) \\\\\n",
"\\frac{1}{\\phi} d\\phi &= \\frac{i}{\\hbar}E dt \\\\ \n",
"\\int \\frac{1}{\\phi} d\\phi = \\int \\frac{i}{\\hbar}E dt \\\\\n",
"\\ln \\phi(t) &= \\frac{iEt}{\\hbar} + \\text{constant of integration}\\\\\n",
"\\phi(t) &\\propto e^{\\frac{iEt}{\\hbar}}\n",
"\\end{align}\n",
"Therefore,\n",
"$$ \\Psi(x,t) = \\psi(x) \\phi(t) \\propto \\psi(x) e^{\\frac{iEt}{\\hbar}} $$\n",
"The constant of proportionality is not important here, because the time-dependent and time-independent Schrödinger equations are satisfied even if $\\Psi(x)$ and $\\phi(t)$ are multiplied by constants.\n"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"## Summary\n",
"When solving the time-independent Schrödinger equation, we need to find all the possible values of E that solve this equation. There are almost always very many different values of E that work and, for this reason, we usually label the solutions to the equation, which are the eigenfunctions (eigenvectors) and eigenvalues of the Hamiltonian operator, accordingly:\n",
"$$ \\hat{H}(x)\\psi_k(x) = E_k \\psi_k(x) $$\n",
"This course is largely about the “art” of solving the TISE\n",
"\n",
"Sometimes in chemistry we cannot use the time-independent Schrödinger equation, and need to use the time-dependent Schrödinger equation instead. This is particular important in spectroscopy, where we want to understand how a molecule responds to a time-dependent external electromagnetic field.\n"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"## 🪞 Self-Reflection\n",
"The Hamiltonian operator enters the Schrödinger equation in a special way, and in this sense is a \"more important\" operator than other quantum-mechanical operators. Why do you think this is the case? "
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"## 🤔 Thought-Provoking Questions\n",
"- Once upon a time, one of the professors at McMaster was asked \"why is the quantum-mechanical momentum operator associated with the derivative\" and, unable to answer the student, sent the student to me for clarification. Can you justify why the quantum mechanical momentum operator is $\\hat{p} = -i \\hbar \\nabla$?\n",
"- Once upon a time, a person interviewing for a job as a quantum chemist was asked how it was possible that the kinetic energy of a quantum system was always positive, given that the quantum-mechanical operator for the kinetic energy was \"the negative of something squared\", namely, $\\hat{T} = \\tfrac{-1}{2m}\\nabla^2$. How would you answer this question?\n",
"- Suppose that one has an electron that is tethered to the origin by a spring, so that the force the electrons feels towards the origin, $x=0$, increases proportionally to its distance, $F = -kx$. What would the Hamiltonian be?\n",
"- What are the eigenfunctions of the kinetic-energy operator?\n",
"- What are the eigenfunctions of the momentum operator?\n",
"- Show, mathematically, that the eigenvectors and eigenvalues of a separable operator, $\\hat{A}(x,y) = \\hat{A}(x) + \\hat{A}(y)$ are given by the above expressions."
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"## 🔁 Recapitulation\n",
"- Write the quantum-mechanical operator for the energy, kinetic-energy, potential energy, and momentum.\n",
"- Write the time-independent and time-dependent Schrödinger equations.\n",
"- What is the form of the time-dependent wavefunction that solves the TDSE when the Hamiltonian is time-independent.\n",
"\n",
"## 🔮 Next Up...\n",
"- Study a particle confined to a box with infinite sides.\n",
"- Explore the postulates of quantum mechanics."
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"## 📚 References\n",
"My favorite sources for this material are:\n",
"- R. Eisberg and R. Resnick, Quantum Physics of Atoms, Molecules, Solids, Nuclei, and Particles (Wiley, New York, 1974)\n",
"- R. Dumont, [An Emergent Reality, Part 2: Quantum Mechanics](https://github.com/PaulWAyers/IntroQChem/blob/main/documents/DumontBook.pdf?raw=true) (Chapters 1 and 2).\n",
"- Also see my (pdf) class [notes](https://github.com/PaulWAyers/IntroQChem/blob/main/documents/IntroQM.pdf?raw=true).\n",
"\n",
"There are also some excellent wikipedia articles:\n",
"- [Plane Wave](https://en.wikipedia.org/wiki/Sinusoidal_plane_wave)\n",
"- [Wavelength-Frequency Relation](https://en.wikipedia.org/wiki/Wavelength)\n",
"- Schrödinger Equation\n",
" - [Schrödinger Equation](http://en.wikipedia.org/wiki/Schrodinger_equation)\n",
" - [Theoretical Justification of S.E.](http://en.wikipedia.org/wiki/Theoretical_and_experimental_justification_for_the_Schr%C3%B6dinger_equation).\n",
" - [Momentum Operator](http://en.wikipedia.org/wiki/Momentum_operator)"
]
}
],
"metadata": {
"kernelspec": {
"display_name": "Python 3",
"language": "python",
"name": "python3"
},
"language_info": {
"codemirror_mode": {
"name": "ipython",
"version": 3
},
"file_extension": ".py",
"mimetype": "text/x-python",
"name": "python",
"nbconvert_exporter": "python",
"pygments_lexer": "ipython3",
"version": "3.8.3"
},
"latex_envs": {
"LaTeX_envs_menu_present": true,
"autoclose": false,
"autocomplete": true,
"bibliofile": "biblio.bib",
"cite_by": "apalike",
"current_citInitial": 1,
"eqLabelWithNumbers": true,
"eqNumInitial": 1,
"hotkeys": {
"equation": "Ctrl-E",
"itemize": "Ctrl-I"
},
"labels_anchors": false,
"latex_user_defs": false,
"report_style_numbering": false,
"user_envs_cfg": false
},
"toc": {
"base_numbering": "2",
"nav_menu": {},
"number_sections": true,
"sideBar": true,
"skip_h1_title": false,
"title_cell": "Table of Contents",
"title_sidebar": "Contents",
"toc_cell": true,
"toc_position": {},
"toc_section_display": true,
"toc_window_display": true
}
},
"nbformat": 4,
"nbformat_minor": 4
}