{
"cells": [
{
"cell_type": "markdown",
"metadata": {},
"source": [
"# One Degree-of-Freedom (DoF) Hamiltonian Bifurcation of Equilibria"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"(ADD INTRODUCTORY LANGUAGE ABOUT PROBLEM DEVELOPMENT)\n",
"We will now consider two examples of bifurcation of equilibria in two dimensional Hamiltonian system; in particular, the Hamiltonian saddle-node and Hamiltonian pitchfork bifurcations. "
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"## Hamiltonian saddle-node bifurcation"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"We consider the Hamiltonian:\n",
"\n",
"\\begin{equation}\n",
"H (q, p) = \\frac{p^2}{2} - \\lambda q + \\frac{q^3}{3}, \\quad (q, p) \\in \\mathbb{R}^2.\n",
"\\label{eq:hamApp13}\n",
"\\end{equation}\n",
"\n",
"\n",
"where $\\lambda$ is considered to be a parameter that can be varied. From this Hamiltonian, we derive Hamilton's equations:\n",
"\n",
"\\begin{eqnarray}\n",
"\\dot{q} & = & \\frac{\\partial H}{\\partial p} = p, \\nonumber \\\\\n",
"\\dot{p} & = & -\\frac{\\partial H}{\\partial q} =\\lambda - q^2.\n",
"\\label{eq:hamApp14}\n",
"\\end{eqnarray}"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"### Revealing the Phase Space Structures and their implications for Reaction Dynamics"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"The fixed points for \\eqref{eq:hamApp14} are:\n",
"\n",
"\\begin{equation}\n",
"(q, p) = (\\pm\\sqrt{\\lambda}, 0),\n",
"\\end{equation}\n",
"\n",
"\n",
"from which it follows that there are no fixed points for $\\lambda <0$, one fixed point for $\\lambda =0$, and two fixed points for $\\lambda >0$. This is the scenario for a saddle-node bifurcation. \n",
"\n",
"Next we examine the stability of the fixed points. The Jacobian of \\eqref{eq:hamApp14} is given by:\n",
"\n",
"\\begin{equation}\n",
"J =\\left(\n",
"\\begin{array}{cc} \n",
"0 & 1\\\\\n",
"-2 q & 0\n",
"\\end{array}\n",
"\\right).\n",
"\\label{eq:hamApp15}\n",
"\\end{equation}\n",
"\n",
"\n",
"The eigenvalues of this matrix are:\n",
"\n",
"\\begin{equation}\n",
"\\Lambda_{1, 2} = \\pm \\sqrt{-2q}.\n",
"\\end{equation}\n",
"\n",
"\n",
"Hence $(q, p) = (-\\sqrt{\\lambda}, 0)$ is a saddle, $(q, p) = (\\sqrt{\\lambda}, 0)$ is a center, and $(q, p) = (0, 0)$ has two zero eigenvalues. The phase portraits are shown in Fig. [fig:1](#fig:appC_fig3)."
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"\n",
"\n",
"\n",
"