Loading video...

Video Failed to Load

Go Home

What Shape Does Chaos Take Before It Looks Random? Take one particle moving in the Hénon-Heiles Hamiltonian H = ½(pₓ² + pᵧ²) + ½(x² + y²) + λ(x²y - ⅓y³) Its motion is forced by Hamilton’s equations ẋ = pₓ ẏ = pᵧ ṗₓ = -x - 2λxy ṗᵧ...

mathelirium's profile picture

Mathelirium

36,596 subscribers

26,893 views • 1 month ago •via X (Twitter)

LIVE

Anya Rossi• Live Now

Private livecam show

0 Comments

No comments available

Comments from the original post will appear here

Related Videos

Phase Space Is the Real Stage of Dynamics. So, Stop Watching the Object and Watch the State. Ordinary Space tells you where something is. Phase Space tells you what state it is in. Position alone is not enough because you also need the variable that tells you where the system is trying to go. In the animation, the left panel shows the particle moving in Ordinary Space, while the right panel shows the same system as a moving point in Phase Space. For one degree of freedom, the state is (x, p) and as time evolves, that state traces a trajectory t ↦ (x(t), p(t)) Instead of thinking of a differential equation as just a formula, Phase Space lets you see it as a flow on the space of states. For this lecture we use ẋ = p ṗ = x − x³ with Hamiltonian H(x,p) = 0.5 p² + 0.25 x⁴ − 0.5 x² #PhaseSpace #DynamicalSystems #NonlinearDynamics #HamiltonianMechanics #Mathematics
0:30

Phase Space Is the Real Stage of Dynamics. So, Stop Watching the Object and Watch the State. Ordinary Space tells you where something is. Phase Space tells you what state it is in. Position alone is not enough because you also need the variable that tells you where the system is trying to go. In the animation, the left panel shows the particle moving in Ordinary Space, while the right panel shows the same system as a moving point in Phase Space. For one degree of freedom, the state is (x, p) and as time evolves, that state traces a trajectory t ↦ (x(t), p(t)) Instead of thinking of a differential equation as just a formula, Phase Space lets you see it as a flow on the space of states. For this lecture we use ẋ = p ṗ = x − x³ with Hamiltonian H(x,p) = 0.5 p² + 0.25 x⁴ − 0.5 x² #PhaseSpace #DynamicalSystems #NonlinearDynamics #HamiltonianMechanics #Mathematics

Mathelirium

30,340 views • 4 months ago

When the Lorenz Attractor Becomes a Cosmos Many of us have seen plots of the famous Lorenz system that famous butterfly-shaped chaos: ẋ = σ(y - x) ẏ = x(ρ - z) - y ż = xy - βz But here, we went one layer deeper. At every point (x, y, z), we computed the local Jacobian matrix, J(x, y, z) = [ -σ σ 0 ] [ ρ - z -1 - x ] [ y x - β ] Then we found its eigenvalues λ and plotted them in the complex plane. Re(λ) tells us whether nearby motion is locally growing or decaying. Im(λ) tells us whether that local motion is rotating. That gives us the spectral flow of the Lorenz system and not the attractor itself, but the hidden local linear geometry driving its instability. When I saw that flow, it looked strangely cosmic with two mirrored swarms curling like particles trapped in a mathematical whirlpool. So I mapped that spectral flow onto our living fabric, not as Einstein Spacetime this time, but as a creative geometric stage for chaos itself. #LorenzSystem #ChaosTheory #DynamicalSystems #Eigenvalues #SpectralFlow #ComplexPlane #NonlinearDynamics #Mathematics #Physics
0:40

When the Lorenz Attractor Becomes a Cosmos Many of us have seen plots of the famous Lorenz system that famous butterfly-shaped chaos: ẋ = σ(y - x) ẏ = x(ρ - z) - y ż = xy - βz But here, we went one layer deeper. At every point (x, y, z), we computed the local Jacobian matrix, J(x, y, z) = [ -σ σ 0 ] [ ρ - z -1 - x ] [ y x - β ] Then we found its eigenvalues λ and plotted them in the complex plane. Re(λ) tells us whether nearby motion is locally growing or decaying. Im(λ) tells us whether that local motion is rotating. That gives us the spectral flow of the Lorenz system and not the attractor itself, but the hidden local linear geometry driving its instability. When I saw that flow, it looked strangely cosmic with two mirrored swarms curling like particles trapped in a mathematical whirlpool. So I mapped that spectral flow onto our living fabric, not as Einstein Spacetime this time, but as a creative geometric stage for chaos itself. #LorenzSystem #ChaosTheory #DynamicalSystems #Eigenvalues #SpectralFlow #ComplexPlane #NonlinearDynamics #Mathematics #Physics

Mathelirium

22,460 views • 1 month ago

Warmup to Statistical Mechanics What Exactly is a Hamiltonian A System? In ordinary Mechanics, you might begin with position and velocity. Hamiltonian Mechanics rewrites the same motion in a different language. Instead of position and velocity, it uses position and momentum. We write the position variables as q and the momentum variables as p. Then the full state of the system at one instant is (q, p) That pair is one point in phase space. Why do we do this? Because in these variables, the equations of motion take a remarkably clean form. Everything is generated by one single function, the Hamiltonian H(q, p) and in the simplest cases this Hamiltonian is just the total energy written in terms of position and momentum. So if you know H, you know the dynamics. You might wonder, but how can one function generate motion? The rule is dqᵢ/dt = ∂H/∂pᵢ dpᵢ/dt = −∂H/∂qᵢ These are Hamilton’s equations. Now read them slowly 😄 The rate of change of position comes from differentiating H with respect to momentum. The rate of change of momentum comes from differentiating H with respect to position, with a minus sign. This constitutes the whole engine. A simple example makes this less abstract: Take one particle of mass m moving in a potential V(q). Then the Hamiltonian is H(q, p) = p²/(2m) + V(q) The first term is kinetic energy. The second term is potential energy. Now apply Hamilton’s equations. First, dq/dt = ∂H/∂p = p/m So momentum tells you how position changes. Second, dp/dt = −∂H/∂q = −dV/dq Thus, momentum changes because of force. If you now combine these two equations, you recover ordinary Newtonian mechanics. Since p = m dq/dt, we get m d²q/dt² = −dV/dq So, Hamiltonian mechanics is not a different theory. It is the same mechanics, written in a form that exposes its geometric structure much more clearly. The animation The full 3D surface is the Hamiltonian itself, the energy landscape H(q, p). The floor underneath is phase space, marked by energy contours and the local flow field. The bright moving point is one actual state (q(t), p(t)) evolving under Hamilton’s equations. Its trail shows that the motion is not arbitrary. It is guided everywhere by the geometry of the same single function H. The render is doing more than illustrating a particle moving, it is showing how one function organizes the whole phase-space motion. The math breakdown: Start with one degree of freedom. The state is described by position q and momentum p. So the system lives in a two-dimensional phase space with coordinates (q, p) Now choose a Hamiltonian H(q, p) Think of H as the energy function. In many standard systems, H(q, p) = kinetic energy + potential energy For a particle of mass m in a potential V(q), this becomes H(q, p) = p²/(2m) + V(q) Hamilton’s equations say dq/dt = ∂H/∂p dp/dt = −∂H/∂q Now substitute this specific H. First compute the p derivative: ∂H/∂p = ∂/∂p (p²/(2m) + V(q)) = p/m So dq/dt = p/m Now compute the q derivative: ∂H/∂q = ∂/∂q (p²/(2m) + V(q)) = dV/dq So dp/dt = −dV/dq These two first-order equations completely determine the motion. Now, connect this back to Newton’s law. From dq/dt = p/m we get p = m dq/dt Differentiate both sides with respect to time: dp/dt = m d²q/dt² But Hamilton’s second equation gives dp/dt = −dV/dq So , together they imply m d²q/dt² = −dV/dq This is exactly Newton’s second law for motion in the potential V(q). Thus, Hamilton’s equations do not replace mechanic, they reorganize it. #HamiltonianMechanics #PhaseSpace #ClassicalMechanics #MathematicalPhysics #DifferentialEquations #Mathematics #Physics
1:00

Warmup to Statistical Mechanics What Exactly is a Hamiltonian A System? In ordinary Mechanics, you might begin with position and velocity. Hamiltonian Mechanics rewrites the same motion in a different language. Instead of position and velocity, it uses position and momentum. We write the position variables as q and the momentum variables as p. Then the full state of the system at one instant is (q, p) That pair is one point in phase space. Why do we do this? Because in these variables, the equations of motion take a remarkably clean form. Everything is generated by one single function, the Hamiltonian H(q, p) and in the simplest cases this Hamiltonian is just the total energy written in terms of position and momentum. So if you know H, you know the dynamics. You might wonder, but how can one function generate motion? The rule is dqᵢ/dt = ∂H/∂pᵢ dpᵢ/dt = −∂H/∂qᵢ These are Hamilton’s equations. Now read them slowly 😄 The rate of change of position comes from differentiating H with respect to momentum. The rate of change of momentum comes from differentiating H with respect to position, with a minus sign. This constitutes the whole engine. A simple example makes this less abstract: Take one particle of mass m moving in a potential V(q). Then the Hamiltonian is H(q, p) = p²/(2m) + V(q) The first term is kinetic energy. The second term is potential energy. Now apply Hamilton’s equations. First, dq/dt = ∂H/∂p = p/m So momentum tells you how position changes. Second, dp/dt = −∂H/∂q = −dV/dq Thus, momentum changes because of force. If you now combine these two equations, you recover ordinary Newtonian mechanics. Since p = m dq/dt, we get m d²q/dt² = −dV/dq So, Hamiltonian mechanics is not a different theory. It is the same mechanics, written in a form that exposes its geometric structure much more clearly. The animation The full 3D surface is the Hamiltonian itself, the energy landscape H(q, p). The floor underneath is phase space, marked by energy contours and the local flow field. The bright moving point is one actual state (q(t), p(t)) evolving under Hamilton’s equations. Its trail shows that the motion is not arbitrary. It is guided everywhere by the geometry of the same single function H. The render is doing more than illustrating a particle moving, it is showing how one function organizes the whole phase-space motion. The math breakdown: Start with one degree of freedom. The state is described by position q and momentum p. So the system lives in a two-dimensional phase space with coordinates (q, p) Now choose a Hamiltonian H(q, p) Think of H as the energy function. In many standard systems, H(q, p) = kinetic energy + potential energy For a particle of mass m in a potential V(q), this becomes H(q, p) = p²/(2m) + V(q) Hamilton’s equations say dq/dt = ∂H/∂p dp/dt = −∂H/∂q Now substitute this specific H. First compute the p derivative: ∂H/∂p = ∂/∂p (p²/(2m) + V(q)) = p/m So dq/dt = p/m Now compute the q derivative: ∂H/∂q = ∂/∂q (p²/(2m) + V(q)) = dV/dq So dp/dt = −dV/dq These two first-order equations completely determine the motion. Now, connect this back to Newton’s law. From dq/dt = p/m we get p = m dq/dt Differentiate both sides with respect to time: dp/dt = m d²q/dt² But Hamilton’s second equation gives dp/dt = −dV/dq So , together they imply m d²q/dt² = −dV/dq This is exactly Newton’s second law for motion in the potential V(q). Thus, Hamilton’s equations do not replace mechanic, they reorganize it. #HamiltonianMechanics #PhaseSpace #ClassicalMechanics #MathematicalPhysics #DifferentialEquations #Mathematics #Physics

Mathelirium

50,635 views • 4 months ago

A Hamiltonian system is a way of describing motion where position and momentum evolve together as one coupled system. The plot shows an energy landscape in phase space, with the motion of the system traced directly on top of it and projected onto the underlying phase portrait. #HamiltonianSystems #PhaseSpace #PhysicsSimulation #DynamicalSystems #MathematicalPhysics
0:40

A Hamiltonian system is a way of describing motion where position and momentum evolve together as one coupled system. The plot shows an energy landscape in phase space, with the motion of the system traced directly on top of it and projected onto the underlying phase portrait. #HamiltonianSystems #PhaseSpace #PhysicsSimulation #DynamicalSystems #MathematicalPhysics

Mathelirium

35,799 views • 2 months ago

Normalize Working in Phase Space Ordinary space shows where the particle is. Phase space shows its complete state, position and momentum together, revealing the hidden geometry of its motion.
0:30

Normalize Working in Phase Space Ordinary space shows where the particle is. Phase space shows its complete state, position and momentum together, revealing the hidden geometry of its motion.

Mathelirium

103,032 views • 29 days ago

We know gravity pulls things together in space...but have you ever seen what it does in phase space? In this scene we're not watching a blob of matter in x–space. We're watching a thin probability sheet in (x, v).... position–velocity space... being pulled by its own gravity. In the collisionless limit, Liouville’s theorem says this sheet can stretch, fold, and swirl, but never tear or change area. Collapse in x is "paid for" by stretching in v. Numerically, this is a 1D cold-dark-matter toy universe: a sheet of particles evolved under the Vlasov-Poisson system. We solve Poisson on a periodic line, get the self-consistent gravitational field, then integrate Hamilton’s equations so the sheet folds into a multi-layer spiral. The glowing filaments are where the sheet keeps revisiting, carving out dense caustics in phase space. It's gravity not as "stuff falling down," but as area-preserving origami in phase space. #Mathelirium #VlasovPoisson #PhaseSpace #Liouville #DarkMatter #HamiltonianFlow #ProbabilityFlow
0:40

We know gravity pulls things together in space...but have you ever seen what it does in phase space? In this scene we're not watching a blob of matter in x–space. We're watching a thin probability sheet in (x, v).... position–velocity space... being pulled by its own gravity. In the collisionless limit, Liouville’s theorem says this sheet can stretch, fold, and swirl, but never tear or change area. Collapse in x is "paid for" by stretching in v. Numerically, this is a 1D cold-dark-matter toy universe: a sheet of particles evolved under the Vlasov-Poisson system. We solve Poisson on a periodic line, get the self-consistent gravitational field, then integrate Hamilton’s equations so the sheet folds into a multi-layer spiral. The glowing filaments are where the sheet keeps revisiting, carving out dense caustics in phase space. It's gravity not as "stuff falling down," but as area-preserving origami in phase space. #Mathelirium #VlasovPoisson #PhaseSpace #Liouville #DarkMatter #HamiltonianFlow #ProbabilityFlow

Mathelirium

47,636 views • 8 months ago

Quantum Mechanics Series Lecture 4 Lecture 1 established that ρ(x,t) = |ψ(x,t)|² behaves like a conserved probability density. Lecture 2 showed what drives that flow. We also saw that writing ψ = r exp(iθ) makes the probability current proportional to the phase gradient, making it clear that phase geometry literally steers the motion. Lecture 3 then showed that the centroid of that flow can move almost classically when the packet is tight and the external potential is smooth. However, that raises yet another question. If the centroid can look classical, why does the full wave still spread, bend, split, and interfere in ways no classical particle cloud would? This is because the wave is not driven only by the external potential. It is also driven by its own curvature. Write ψ(x,t) = r(x,t) exp(iθ(x,t)) with ρ = r². Then Schrödinger’s equation gives two coupled real equations. One is the continuity equation you already know. The other looks like a Hamilton-Jacobi equation, but with one extra term: Q = −(1/2m) ∇²r / r This is the so-called Quantum Potential. It depends entirely on how the amplitude bends across space. So, the wave is being shaped not only by V(x,t), but also by the geometry of its own envelope. In the animation, the upper surface is still |ψ| and its skin is still colored by arg(ψ). The glowing threads still trace the probability current. But now a second membrane hangs underneath. That lower membrane encodes the quantum potential Q itself. The porcelain bead marks the quantum centroid. The amber bead follows a classical centroid under the same external V. When those paths separate, the lower membrane tells you why. The difference is not magic but the extra term classical mechanics does not have. The math breakdown: Start from Schrödinger evolution in units with ħ = 1: i ∂ψ/∂t = [ −(1/2m) ∇² + V(x,t) ] ψ Write the state in polar form: ψ = r exp(iθ) Then ρ = |ψ|² = r² From the imaginary part, you recover probability conservation: ∂ρ/∂t + ∇·j = 0 with j = (1/m) Im(ψ* ∇ψ) = (ρ/m) ∇θ So the local velocity field is v = j / ρ = ∇θ / m Now take the real part of Schrödinger’s equation. That gives ∂θ/∂t + |∇θ|² / (2m) + V + Q = 0 where Q = −(1/2m) ∇²r / r This is the classical Hamilton-Jacobi equation with one extra term. That extra term is what makes quantum motion locally different from classical motion. Take a gradient of that phase equation and use v = ∇θ / m. Then the flow obeys an Euler-like equation: ∂v/∂t + (v·∇)v = −(1/m) ∇(V + Q) In other words, there are really two forces in the problem. One comes from the external potential V. The other comes from the wave’s own curvature through Q. That is why Ehrenfest is only approximate. The centroid can still satisfy d⟨x⟩/dt = ⟨p⟩/m d⟨p⟩/dt = −⟨∇V⟩ but the internal shape of the packet evolves under the combined influence of V and Q. When the packet stays broad and smooth, Q is gentle and the motion looks more classical. When the packet develops sharp curvature or interference structure, Q becomes strong and the classical picture breaks down. That is what this scene is designed to show live. #QuantumMechanics #Wavefunction #SchrodingerEquation #BornRule #ProbabilityCurrent #ContinuityEquation #Phase #EhrenfestTheorem #QuantumPotential #Madelung #HamiltonJacobi #MathematicalPhysics #Mathematics #Physics
1:00

Quantum Mechanics Series Lecture 4 Lecture 1 established that ρ(x,t) = |ψ(x,t)|² behaves like a conserved probability density. Lecture 2 showed what drives that flow. We also saw that writing ψ = r exp(iθ) makes the probability current proportional to the phase gradient, making it clear that phase geometry literally steers the motion. Lecture 3 then showed that the centroid of that flow can move almost classically when the packet is tight and the external potential is smooth. However, that raises yet another question. If the centroid can look classical, why does the full wave still spread, bend, split, and interfere in ways no classical particle cloud would? This is because the wave is not driven only by the external potential. It is also driven by its own curvature. Write ψ(x,t) = r(x,t) exp(iθ(x,t)) with ρ = r². Then Schrödinger’s equation gives two coupled real equations. One is the continuity equation you already know. The other looks like a Hamilton-Jacobi equation, but with one extra term: Q = −(1/2m) ∇²r / r This is the so-called Quantum Potential. It depends entirely on how the amplitude bends across space. So, the wave is being shaped not only by V(x,t), but also by the geometry of its own envelope. In the animation, the upper surface is still |ψ| and its skin is still colored by arg(ψ). The glowing threads still trace the probability current. But now a second membrane hangs underneath. That lower membrane encodes the quantum potential Q itself. The porcelain bead marks the quantum centroid. The amber bead follows a classical centroid under the same external V. When those paths separate, the lower membrane tells you why. The difference is not magic but the extra term classical mechanics does not have. The math breakdown: Start from Schrödinger evolution in units with ħ = 1: i ∂ψ/∂t = [ −(1/2m) ∇² + V(x,t) ] ψ Write the state in polar form: ψ = r exp(iθ) Then ρ = |ψ|² = r² From the imaginary part, you recover probability conservation: ∂ρ/∂t + ∇·j = 0 with j = (1/m) Im(ψ* ∇ψ) = (ρ/m) ∇θ So the local velocity field is v = j / ρ = ∇θ / m Now take the real part of Schrödinger’s equation. That gives ∂θ/∂t + |∇θ|² / (2m) + V + Q = 0 where Q = −(1/2m) ∇²r / r This is the classical Hamilton-Jacobi equation with one extra term. That extra term is what makes quantum motion locally different from classical motion. Take a gradient of that phase equation and use v = ∇θ / m. Then the flow obeys an Euler-like equation: ∂v/∂t + (v·∇)v = −(1/m) ∇(V + Q) In other words, there are really two forces in the problem. One comes from the external potential V. The other comes from the wave’s own curvature through Q. That is why Ehrenfest is only approximate. The centroid can still satisfy d⟨x⟩/dt = ⟨p⟩/m d⟨p⟩/dt = −⟨∇V⟩ but the internal shape of the packet evolves under the combined influence of V and Q. When the packet stays broad and smooth, Q is gentle and the motion looks more classical. When the packet develops sharp curvature or interference structure, Q becomes strong and the classical picture breaks down. That is what this scene is designed to show live. #QuantumMechanics #Wavefunction #SchrodingerEquation #BornRule #ProbabilityCurrent #ContinuityEquation #Phase #EhrenfestTheorem #QuantumPotential #Madelung #HamiltonJacobi #MathematicalPhysics #Mathematics #Physics

Mathelirium

20,456 views • 4 months ago

Ask anyone who’s taken a course in Ordinary Differential Equations (ODEs) what a solution to an ODE represents geometrically, and most of them won’t have a clean answer. When I first took ordinary differential equations, the pattern was always the same. Early on it turns into a speedrun of methods: separation of variables, integrating factors, variation of parameters, Bernoulli, exact equations. Then pretty quickly the course slides into hammer-picking. Spot the form, apply the recipe, move on. Too mechanical! And the real problem is what you don’t walk away with. You leave with a toolkit, but without a feel for what a differential equation even is, especially geometrically. That matters because in real modeling the equations you meet are rarely nice enough to reward memorised recipes. So you get trained to solve toy forms, while the actual subject stays blurry. The behavior. The flow. The shape of solutions. It wasn't until I watched the first lecture of Professor Arthur Mattuck that I realized I didn’t actually know what a solution to a differential equation represents geometrically. His point is almost embarrassingly simple. A first-order ODE is a slope field, and a solution is a curve that stays tangent to that field everywhere. The math breakdown: Write the ODE as dy/dx = f(x,y). At each point (x,y), attach a tiny line segment with slope f(x,y). A function y = y₁(x) is a solution exactly when its graph follows those slopes. At every x, the slope of the curve equals the slope prescribed by the field at the point on the curve. That’s the one line that ties both viewpoints together: y₁′(x) = f(x, y₁(x)). So solving the ODE and drawing an integral curve are the same statement in two languages. Once you see that, you stop obsessing over whether you can write y(x) in closed form. You start asking the questions that actually matter. Where do solutions flow. Where do they get trapped. Where do they blow up. Where does existence or uniqueness fail because the field isn’t even defined? That’s the perspective shift I wish every ODE course forces early. It’s also why I keep pairing math with animation. #DifferentialEquations #ODEs #VectorFields #AppliedMathematics #Mathematics #
20:13

Ask anyone who’s taken a course in Ordinary Differential Equations (ODEs) what a solution to an ODE represents geometrically, and most of them won’t have a clean answer. When I first took ordinary differential equations, the pattern was always the same. Early on it turns into a speedrun of methods: separation of variables, integrating factors, variation of parameters, Bernoulli, exact equations. Then pretty quickly the course slides into hammer-picking. Spot the form, apply the recipe, move on. Too mechanical! And the real problem is what you don’t walk away with. You leave with a toolkit, but without a feel for what a differential equation even is, especially geometrically. That matters because in real modeling the equations you meet are rarely nice enough to reward memorised recipes. So you get trained to solve toy forms, while the actual subject stays blurry. The behavior. The flow. The shape of solutions. It wasn't until I watched the first lecture of Professor Arthur Mattuck that I realized I didn’t actually know what a solution to a differential equation represents geometrically. His point is almost embarrassingly simple. A first-order ODE is a slope field, and a solution is a curve that stays tangent to that field everywhere. The math breakdown: Write the ODE as dy/dx = f(x,y). At each point (x,y), attach a tiny line segment with slope f(x,y). A function y = y₁(x) is a solution exactly when its graph follows those slopes. At every x, the slope of the curve equals the slope prescribed by the field at the point on the curve. That’s the one line that ties both viewpoints together: y₁′(x) = f(x, y₁(x)). So solving the ODE and drawing an integral curve are the same statement in two languages. Once you see that, you stop obsessing over whether you can write y(x) in closed form. You start asking the questions that actually matter. Where do solutions flow. Where do they get trapped. Where do they blow up. Where does existence or uniqueness fail because the field isn’t even defined? That’s the perspective shift I wish every ODE course forces early. It’s also why I keep pairing math with animation. #DifferentialEquations #ODEs #VectorFields #AppliedMathematics #Mathematics #

Mathelirium

40,739 views • 6 months ago

Erwin Schrödinger’s 1926 equation changed the game by turning Quantum Mechanics into wave dynamics. That move gave Physics a new way to think. Instead of forcing a particle onto one sharp path, it lets a complex wavefunction evolve in time, with its shape and phase holding the structure of the phenomenon. What makes this so striking is that Schrödinger’s equation does not start from vague mystery. It starts from a precise and daring idea The state of a particle is a complex field ψ(x,t), and the dynamics must push ψ forward in a way that preserves total probability.
0:33

Erwin Schrödinger’s 1926 equation changed the game by turning Quantum Mechanics into wave dynamics. That move gave Physics a new way to think. Instead of forcing a particle onto one sharp path, it lets a complex wavefunction evolve in time, with its shape and phase holding the structure of the phenomenon. What makes this so striking is that Schrödinger’s equation does not start from vague mystery. It starts from a precise and daring idea The state of a particle is a complex field ψ(x,t), and the dynamics must push ψ forward in a way that preserves total probability.

Mathelirium

54,730 views • 4 months ago

Lecture 2 on our Quantum Mechanics Series Schrödinger’s equation doesn’t start from mystery. It starts from a very specific bet…the state of a particle is a complex field ψ(x,t), and whatever dynamics we write down must move ψ forward in time in a way that preserves total probability. We ask a basic question…what equation should ψ satisfy so that |ψ|² behaves like a conserved density, the way mass density does in fluid flow? What is ψ? Think of ψ(x,t) as the amplitude assigned to “the particle is at position x at time t”. It’s not a probability. It’s the object you add first, and only at the end do you square p(x,t) = |ψ(x,t)|² Because ψ is complex, it has magnitude and phase. Write it in polar form ψ(x,t) = r(x,t) exp(i θ(x,t)) Then r² = |ψ|² is the density, and θ will end up controlling flow (the probability current). Where does Schrödinger’s equation come from? Start with two empirical inputs about waves and particles: E = ħ ω p = ħ k Here ħ (“h-bar”) is Planck’s constant divided by 2π. It’s the unit conversion factor between the wave description (frequency ω, wavevector k) and the particle description (energy E, momentum p). In units, ħ has units of joule-seconds, so multiplying ω (1/seconds) gives energy (joules), and multiplying k (1/meters) gives momentum (kg·m/s). It’s the number that tells you how much energy or momentum you get per unit frequency or wavenumber. A plane wave with angular frequency ω and wavevector k is ψ(x,t) = A exp(i(k·x − ω t)) Now notice what derivatives do to this wave: ∂ψ/∂t = −i ω ψ ∇ψ = i k ψ ∇²ψ = −|k|² ψ Multiply those identities by ħ: i ħ ∂ψ/∂t = ħ ω ψ = E ψ −i ħ ∇ψ = ħ k ψ = p ψ −ħ² ∇²ψ = ħ² |k|² ψ = p² ψ So for plane waves, the operators Ê = i ħ ∂/∂t p̂ = −i ħ ∇ act like energy and momentum! Now use the classical nonrelativistic energy relation: E = p²/(2m) + V(x) This is bookkeeping for a particle moving slow enough that relativity can be ignored. The term p²/(2m) is kinetic energy. If p = mv, then p²/(2m) = (m²v²)/(2m) = (1/2)mv². The term V(x) is potential energy. It depends on position because forces come from spatially varying energy. A slope in V pushes the particle. Examples: for a charged particle in an electric potential φ(x), V(x) = q φ(x). Near Earth, V(z) = mgz. The point is total energy equals kinetic plus potential. Turn that into an equation for ψ by replacing E and p with the operators above: Ê ψ = (p̂²/(2m) + V) ψ Compute p̂² = (−i ħ ∇)·(−i ħ ∇) = −ħ² ∇², so we get i ħ ∂ψ/∂t = ( −ħ²/(2m) ∇² + V(x) ) ψ That is the time-dependent Schrödinger equation. The derivation here is a controlled heuristic: we matched the plane-wave identities to the measured relations E = ħω and p = ħk, then imposed the same energy bookkeeping as classical mechanics. Why this equation is the right kind of rule If ψ is the state, we need a rule that preserves total probability: ∫ |ψ(x,t)|² dx = 1 Schrödinger evolution does. You can see it by deriving a continuity equation. Let ρ(x,t) = |ψ|² = ψ* ψ. Take a time derivative: ∂ρ/∂t = ψ* ∂ψ/∂t + ψ ∂ψ*/∂t Use Schrödinger and its complex conjugate: ∂ψ/∂t = (1/(i ħ)) ( −ħ²/(2m) ∇²ψ + Vψ ) ∂ψ*/∂t = (−1/(i ħ)) ( −ħ²/(2m) ∇²ψ* + Vψ* ) Plug in. The V terms cancel exactly, and what remains can be rearranged into a divergence: ∂ρ/∂t + ∇·j = 0 where the probability current is j = (ħ/(2mi)) ( ψ* ∇ψ − ψ ∇ψ* ) This is the best way to explain ehat ψ is: |ψ|² behaves like a conserved density, and the phase of ψ is what drives the current j. So in this series, ψ isn’t a slogan. It’s the object whose modulus squared is the density, whose phase generates flow, and whose time evolution is fixed (up to V) by matching wave relations to energy bookkeeping: i ħ ∂ψ/∂t = ( -ħ²/(2m) ∇² + V ) ψ #QuantumMechanics #SchrodingerEquation #WaveFunction #BornRule #Physics #MathematicalPhysics
0:40

Lecture 2 on our Quantum Mechanics Series Schrödinger’s equation doesn’t start from mystery. It starts from a very specific bet…the state of a particle is a complex field ψ(x,t), and whatever dynamics we write down must move ψ forward in time in a way that preserves total probability. We ask a basic question…what equation should ψ satisfy so that |ψ|² behaves like a conserved density, the way mass density does in fluid flow? What is ψ? Think of ψ(x,t) as the amplitude assigned to “the particle is at position x at time t”. It’s not a probability. It’s the object you add first, and only at the end do you square p(x,t) = |ψ(x,t)|² Because ψ is complex, it has magnitude and phase. Write it in polar form ψ(x,t) = r(x,t) exp(i θ(x,t)) Then r² = |ψ|² is the density, and θ will end up controlling flow (the probability current). Where does Schrödinger’s equation come from? Start with two empirical inputs about waves and particles: E = ħ ω p = ħ k Here ħ (“h-bar”) is Planck’s constant divided by 2π. It’s the unit conversion factor between the wave description (frequency ω, wavevector k) and the particle description (energy E, momentum p). In units, ħ has units of joule-seconds, so multiplying ω (1/seconds) gives energy (joules), and multiplying k (1/meters) gives momentum (kg·m/s). It’s the number that tells you how much energy or momentum you get per unit frequency or wavenumber. A plane wave with angular frequency ω and wavevector k is ψ(x,t) = A exp(i(k·x − ω t)) Now notice what derivatives do to this wave: ∂ψ/∂t = −i ω ψ ∇ψ = i k ψ ∇²ψ = −|k|² ψ Multiply those identities by ħ: i ħ ∂ψ/∂t = ħ ω ψ = E ψ −i ħ ∇ψ = ħ k ψ = p ψ −ħ² ∇²ψ = ħ² |k|² ψ = p² ψ So for plane waves, the operators Ê = i ħ ∂/∂t p̂ = −i ħ ∇ act like energy and momentum! Now use the classical nonrelativistic energy relation: E = p²/(2m) + V(x) This is bookkeeping for a particle moving slow enough that relativity can be ignored. The term p²/(2m) is kinetic energy. If p = mv, then p²/(2m) = (m²v²)/(2m) = (1/2)mv². The term V(x) is potential energy. It depends on position because forces come from spatially varying energy. A slope in V pushes the particle. Examples: for a charged particle in an electric potential φ(x), V(x) = q φ(x). Near Earth, V(z) = mgz. The point is total energy equals kinetic plus potential. Turn that into an equation for ψ by replacing E and p with the operators above: Ê ψ = (p̂²/(2m) + V) ψ Compute p̂² = (−i ħ ∇)·(−i ħ ∇) = −ħ² ∇², so we get i ħ ∂ψ/∂t = ( −ħ²/(2m) ∇² + V(x) ) ψ That is the time-dependent Schrödinger equation. The derivation here is a controlled heuristic: we matched the plane-wave identities to the measured relations E = ħω and p = ħk, then imposed the same energy bookkeeping as classical mechanics. Why this equation is the right kind of rule If ψ is the state, we need a rule that preserves total probability: ∫ |ψ(x,t)|² dx = 1 Schrödinger evolution does. You can see it by deriving a continuity equation. Let ρ(x,t) = |ψ|² = ψ* ψ. Take a time derivative: ∂ρ/∂t = ψ* ∂ψ/∂t + ψ ∂ψ*/∂t Use Schrödinger and its complex conjugate: ∂ψ/∂t = (1/(i ħ)) ( −ħ²/(2m) ∇²ψ + Vψ ) ∂ψ*/∂t = (−1/(i ħ)) ( −ħ²/(2m) ∇²ψ* + Vψ* ) Plug in. The V terms cancel exactly, and what remains can be rearranged into a divergence: ∂ρ/∂t + ∇·j = 0 where the probability current is j = (ħ/(2mi)) ( ψ* ∇ψ − ψ ∇ψ* ) This is the best way to explain ehat ψ is: |ψ|² behaves like a conserved density, and the phase of ψ is what drives the current j. So in this series, ψ isn’t a slogan. It’s the object whose modulus squared is the density, whose phase generates flow, and whose time evolution is fixed (up to V) by matching wave relations to energy bookkeeping: i ħ ∂ψ/∂t = ( -ħ²/(2m) ∇² + V ) ψ #QuantumMechanics #SchrodingerEquation #WaveFunction #BornRule #Physics #MathematicalPhysics

Mathelirium

40,835 views • 7 months ago

We know gravity pulls things together in space… but have you ever watched it work in phase space? In this scene we’re looking at a thin probability sheet in (x, v) position-velocity space evolving under its own self-gravity. In the collisionless limit, Liouville’s theorem says the sheet can stretch, fold, and swirl, but it can’t tear and it can’t change area. Any collapse in x is paid for by stretching in v. Numerically this is a 1D cold-dark-matter toy universe: a sheet of particles evolved under the Vlasov-Poisson system. We solve Poisson on a periodic line to get the self-consistent gravitational field, then integrate Hamilton’s equations so the sheet folds into a multi-layer spiral. The glowing filaments are the places the sheet keeps revisiting dense, sharp caustics carved into phase space. Gravity here isn’t stuff falling down. It’s area-preserving origami in (x, v). #Physics #VlasovPoisson #PhaseSpace #Liouville #DarkMatter #HamiltonianFlow #ProbabilityFlow #Mathematics
0:40

We know gravity pulls things together in space… but have you ever watched it work in phase space? In this scene we’re looking at a thin probability sheet in (x, v) position-velocity space evolving under its own self-gravity. In the collisionless limit, Liouville’s theorem says the sheet can stretch, fold, and swirl, but it can’t tear and it can’t change area. Any collapse in x is paid for by stretching in v. Numerically this is a 1D cold-dark-matter toy universe: a sheet of particles evolved under the Vlasov-Poisson system. We solve Poisson on a periodic line to get the self-consistent gravitational field, then integrate Hamilton’s equations so the sheet folds into a multi-layer spiral. The glowing filaments are the places the sheet keeps revisiting dense, sharp caustics carved into phase space. Gravity here isn’t stuff falling down. It’s area-preserving origami in (x, v). #Physics #VlasovPoisson #PhaseSpace #Liouville #DarkMatter #HamiltonianFlow #ProbabilityFlow #Mathematics

Mathelirium

16,437 views • 6 months ago

Ever Heard of a Metric in Differential Geometry? How About a Geodesic? These two ideas are central to how Mathematics and Physics describe Spacetime. A metric is how a surface measures distance and angle. A geodesic is the path that counts as "straight" after that measuring rule is chosen. Once the metric changes from point to point, straightness is no longer something you judge by eye. On a curved surface, the shortest paths are geodesics. They are chosen by the metric. Take a surface written as r(x,y) = (x, y, f(x,y)) Its metric is gᵢⱼ = ∂ᵢr · ∂ⱼr A geodesic q(t) = (x(t), y(t)) obeys q̈ᵏ + Γᵏᵢⱼ q̇ᶦq̇ʲ = 0 where Γᵏᵢⱼ are the Christoffel symbols of the surface metric. Launch many almost-parallel paths. At the start, they look like a clean beam until curvature takes over. Some paths focus, others separate and some slide around ridges. The surface is essentially the rulebook for motion. The shortest path is local. The geometry remembers everything. #DifferentialGeometry #Geodesics #RiemannianGeometry #MathematicalPhysics #Curvature #TensorCalculus #AppliedMathematics #Geometry #ScienceAnimation
0:33

Ever Heard of a Metric in Differential Geometry? How About a Geodesic? These two ideas are central to how Mathematics and Physics describe Spacetime. A metric is how a surface measures distance and angle. A geodesic is the path that counts as "straight" after that measuring rule is chosen. Once the metric changes from point to point, straightness is no longer something you judge by eye. On a curved surface, the shortest paths are geodesics. They are chosen by the metric. Take a surface written as r(x,y) = (x, y, f(x,y)) Its metric is gᵢⱼ = ∂ᵢr · ∂ⱼr A geodesic q(t) = (x(t), y(t)) obeys q̈ᵏ + Γᵏᵢⱼ q̇ᶦq̇ʲ = 0 where Γᵏᵢⱼ are the Christoffel symbols of the surface metric. Launch many almost-parallel paths. At the start, they look like a clean beam until curvature takes over. Some paths focus, others separate and some slide around ridges. The surface is essentially the rulebook for motion. The shortest path is local. The geometry remembers everything. #DifferentialGeometry #Geodesics #RiemannianGeometry #MathematicalPhysics #Curvature #TensorCalculus #AppliedMathematics #Geometry #ScienceAnimation

Mathelirium

22,170 views • 1 month ago

A Triple Spring-Mass Pendulum Just for The Fun of Physics We model the mass positions as r₁ = (x₁, y₁), r₂ = (x₂, y₂), r₃ = (x₃, y₃) with y positive downward in the simulation. The three spring segments are r₁ = r₁ − (0,0) r₁₂ = r₂ − r₁ r₂₃ = r₃ − r₂ Each spring applies a Hooke force along its own direction, F(r; k, L) = −k (|r| − L) (r / |r|) and each mass also feels gravity, Fgᵢ = (0, mᵢ g) So the equations of motion are m₁ r̈₁ = Fg₁ + F(r₁; k₁, L₁) − F(r₁₂; k₂, L₂) m₂ r̈₂ = Fg₂ + F(r₁₂; k₂, L₂) − F(r₂₃; k₃, L₃) m₃ r̈₃ = Fg₃ + F(r₂₃; k₃, L₃) We integrate these ODEs forward with RK4 using a small dt. Different scenes come from changing the mass ratios and from giving the third mass a small initial kick. The yellow trail is the path of the third mass. The inset phase portraits show the same motion in state space: (x₃, vₓ₃) and (−y₃, vᵧ₃) When those curves stay tight and loop cleanly, the motion is more regular. When they spread, fold, and smear, energy is being redistributed across the whole chain between gravity, spring stretch, and kinetic motion. #NonlinearDynamics #CoupledOscillators #SpringPendulum #ChaosInMotion #ComputationalPhysics #MathematicalVisualization
0:49

A Triple Spring-Mass Pendulum Just for The Fun of Physics We model the mass positions as r₁ = (x₁, y₁), r₂ = (x₂, y₂), r₃ = (x₃, y₃) with y positive downward in the simulation. The three spring segments are r₁ = r₁ − (0,0) r₁₂ = r₂ − r₁ r₂₃ = r₃ − r₂ Each spring applies a Hooke force along its own direction, F(r; k, L) = −k (|r| − L) (r / |r|) and each mass also feels gravity, Fgᵢ = (0, mᵢ g) So the equations of motion are m₁ r̈₁ = Fg₁ + F(r₁; k₁, L₁) − F(r₁₂; k₂, L₂) m₂ r̈₂ = Fg₂ + F(r₁₂; k₂, L₂) − F(r₂₃; k₃, L₃) m₃ r̈₃ = Fg₃ + F(r₂₃; k₃, L₃) We integrate these ODEs forward with RK4 using a small dt. Different scenes come from changing the mass ratios and from giving the third mass a small initial kick. The yellow trail is the path of the third mass. The inset phase portraits show the same motion in state space: (x₃, vₓ₃) and (−y₃, vᵧ₃) When those curves stay tight and loop cleanly, the motion is more regular. When they spread, fold, and smear, energy is being redistributed across the whole chain between gravity, spring stretch, and kinetic motion. #NonlinearDynamics #CoupledOscillators #SpringPendulum #ChaosInMotion #ComputationalPhysics #MathematicalVisualization

Mathelirium

36,155 views • 5 months ago

When I first took ordinary differential equations, the pattern was always the same. Week 1 turns into a speedrun of methods: separation of variables, integrating factors, variation of parameters, Bernoulli, exact equations… and by Week 2 or 3 the course has quietly degenerated into hammer-picking. Spot the form, apply the recipe, move on. Mechanical! Fuuuuck!😫😫😫😫 The problem is what you don’t walk away with. You leave with a toolkit, but without a feel for what a differential equation even is, especially geometrically. And that’s a big deal, because in real modeling the equations you meet are rarely nice enough to reward memorized recipes. So you end up trained to solve toy forms, while the actual subject...the behavior, the flow, the shape of solutions stays blurry. This is why I’m biased toward the old-timers. Their old-school way of doing things always surprises me:...they’ll spend time on one idea until it sticks, instead of sprinting through a syllabus checklist. One lecture from them and you start noticing a contrast. A lot of modern teaching feels like "finish the content,". You get marched through techniques, but you’re not left with a single thought that keeps bothering you later...the kind of thought that actually pushes you toward research-level curiosity. MIT OpenCourseWare’s Professor Arthur Mattuck did that to me in his very first ODE lecture. One lecture, and your whole relationship with dy/dx = f(x,y) changes. In this segment, Prof. Mattuck is basically saying: A first-order ODE is a slope field, and a solution is a curve that moves everywhere tangent to that field. The math breakdown Write the ODE as dy/dx = f(x,y). At each point (x,y) you attach a tiny line segment with slope f(x,y). A function y = y₁(x) is a solution exactly when its graph follows those slopes:. At every x, the slope of the curve equals the slope prescribed by the field at the point on the curve. That’s the single line that unifies both viewpoints: y₁′(x) = f(x, y₁(x)). So solving the ODE and drawing an integral curve are the same statement in two languages!👌🏻 Once you see that, you can stop obsessing over whether you can write y(x) in closed form. You can start asking the questions that matter: where do solutions flow, where do they get trapped, where do they blow up, and where does existence/uniqueness fail just because the field isn’t even defined? That’s the perspective shift I wish every ODE course forces early and it’s exactly why I keep pairing math with animation. #DifferentialEquations #ODEs #VectorFields #MathAnimation #Mathematics
20:13

When I first took ordinary differential equations, the pattern was always the same. Week 1 turns into a speedrun of methods: separation of variables, integrating factors, variation of parameters, Bernoulli, exact equations… and by Week 2 or 3 the course has quietly degenerated into hammer-picking. Spot the form, apply the recipe, move on. Mechanical! Fuuuuck!😫😫😫😫 The problem is what you don’t walk away with. You leave with a toolkit, but without a feel for what a differential equation even is, especially geometrically. And that’s a big deal, because in real modeling the equations you meet are rarely nice enough to reward memorized recipes. So you end up trained to solve toy forms, while the actual subject...the behavior, the flow, the shape of solutions stays blurry. This is why I’m biased toward the old-timers. Their old-school way of doing things always surprises me:...they’ll spend time on one idea until it sticks, instead of sprinting through a syllabus checklist. One lecture from them and you start noticing a contrast. A lot of modern teaching feels like "finish the content,". You get marched through techniques, but you’re not left with a single thought that keeps bothering you later...the kind of thought that actually pushes you toward research-level curiosity. MIT OpenCourseWare’s Professor Arthur Mattuck did that to me in his very first ODE lecture. One lecture, and your whole relationship with dy/dx = f(x,y) changes. In this segment, Prof. Mattuck is basically saying: A first-order ODE is a slope field, and a solution is a curve that moves everywhere tangent to that field. The math breakdown Write the ODE as dy/dx = f(x,y). At each point (x,y) you attach a tiny line segment with slope f(x,y). A function y = y₁(x) is a solution exactly when its graph follows those slopes:. At every x, the slope of the curve equals the slope prescribed by the field at the point on the curve. That’s the single line that unifies both viewpoints: y₁′(x) = f(x, y₁(x)). So solving the ODE and drawing an integral curve are the same statement in two languages!👌🏻 Once you see that, you can stop obsessing over whether you can write y(x) in closed form. You can start asking the questions that matter: where do solutions flow, where do they get trapped, where do they blow up, and where does existence/uniqueness fail just because the field isn’t even defined? That’s the perspective shift I wish every ODE course forces early and it’s exactly why I keep pairing math with animation. #DifferentialEquations #ODEs #VectorFields #MathAnimation #Mathematics

Mathelirium

53,338 views • 7 months ago

Russian Genius Found Structure in The Mess of Turbulence In 1941, Andrey Kolmogorov proposed a simple but powerful picture of turbulence: Energy does not disappear where it is injected but moves. Large eddies pass energy to smaller ones, those smaller ones feed even finer motion, and viscosity finally removes it at the tiniest scales. That cascade picture became one of the foundations of modern turbulence theory. The point is not just that turbulence looks messy. The point is that the mess has structure, and Kolmogorov’s great insight was to treat that structure statistically.
0:30

Russian Genius Found Structure in The Mess of Turbulence In 1941, Andrey Kolmogorov proposed a simple but powerful picture of turbulence: Energy does not disappear where it is injected but moves. Large eddies pass energy to smaller ones, those smaller ones feed even finer motion, and viscosity finally removes it at the tiniest scales. That cascade picture became one of the foundations of modern turbulence theory. The point is not just that turbulence looks messy. The point is that the mess has structure, and Kolmogorov’s great insight was to treat that structure statistically.

Mathelirium

58,963 views • 4 months ago

Why Does Quantum Mechanics Use a Complex Wavefunction? Schrödinger’s equation doesn’t start from mystery. It starts from a very specific bet. The state of a particle is a complex field ψ(x,t), and whatever time-evolution rule we choose has to move ψ forward while preserving total probability. So the basic question is simple. What equation should ψ satisfy so that |ψ|² behaves like a conserved density, the way mass density does in fluid flow? What is ψ? Think of ψ(x,t) as an amplitude attached to the statement the particle is at position x at time t. It’s not a probability. It’s the thing you add first, and only at the end do you square it: p(x,t) = |ψ(x,t)|² Because ψ is complex, it has magnitude and phase. Write it as ψ(x,t) = r(x,t) exp(i θ(x,t)) Then r² = |ψ|² is the density, and the phase θ ends up controlling the flow through the probability current. Where does Schrödinger’s equation come from? Start with two empirical inputs that tie waves to particles: E = ħ ω p = ħ k Here ħ is Planck’s constant divided by 2π. It’s the conversion factor between frequency and energy, and between wavenumber and momentum. A plane wave with angular frequency ω and wavevector k is ψ(x,t) = A exp(i(k·x − ωt)) Now watch what derivatives do to this wave: ∂ψ/∂t = −i ω ψ ∇ψ = i k ψ ∇²ψ = −|k|² ψ Multiply by ħ and you get: i ħ ∂ψ/∂t = ħ ω ψ = E ψ −i ħ ∇ψ = ħ k ψ = p ψ −ħ² ∇²ψ = ħ² |k|² ψ = p² ψ So for plane waves, the operators Ê = i ħ ∂/∂t p̂ = −i ħ ∇ act like energy and momentum. Now bring in the classical, nonrelativistic energy bookkeeping: E = p²/(2m) + V(x) Kinetic plus potential. That’s it. Turn it into an equation for ψ by replacing E and p with the operators above: Ê ψ = (p̂²/(2m) + V) ψ Since p̂² = (−i ħ ∇)·(−i ħ ∇) = −ħ² ∇², this becomes i ħ ∂ψ/∂t = ( −ħ²/(2m) ∇² + V(x) ) ψ That’s the time-dependent Schrödinger equation. This derivation is a controlled heuristic. Match the plane-wave identities to the measured relations E = ħω and p = ħk, then impose the same energy bookkeeping you trust in classical mechanics. Why this is the right kind of rule If ψ is the state, we need a rule that preserves total probability: ∫ |ψ(x,t)|² dx = 1 Schrödinger evolution does, and you can see it by deriving a continuity equation. Let ρ(x,t) = |ψ|² = ψ*ψ. Differentiate: ∂ρ/∂t = ψ* ∂ψ/∂t + ψ ∂ψ*/∂t Use Schrödinger and its complex conjugate. The potential terms cancel, and what’s left can be rearranged into ∂ρ/∂t + ∇·j = 0 with probability current j = (ħ/(2mi)) ( ψ* ∇ψ − ψ ∇ψ* ) That’s the cleanest way to say what ψ is. |ψ|² behaves like a conserved density, the phase drives a current, and the time evolution is fixed, up to V, by combining wave relations with energy bookkeeping: i ħ ∂ψ/∂t = ( −ħ²/(2m) ∇² + V ) ψ #QuantumMechanics #SchrodingerEquation #WaveFunction #BornRule #Physics #MathematicalPhysics
1:09

Why Does Quantum Mechanics Use a Complex Wavefunction? Schrödinger’s equation doesn’t start from mystery. It starts from a very specific bet. The state of a particle is a complex field ψ(x,t), and whatever time-evolution rule we choose has to move ψ forward while preserving total probability. So the basic question is simple. What equation should ψ satisfy so that |ψ|² behaves like a conserved density, the way mass density does in fluid flow? What is ψ? Think of ψ(x,t) as an amplitude attached to the statement the particle is at position x at time t. It’s not a probability. It’s the thing you add first, and only at the end do you square it: p(x,t) = |ψ(x,t)|² Because ψ is complex, it has magnitude and phase. Write it as ψ(x,t) = r(x,t) exp(i θ(x,t)) Then r² = |ψ|² is the density, and the phase θ ends up controlling the flow through the probability current. Where does Schrödinger’s equation come from? Start with two empirical inputs that tie waves to particles: E = ħ ω p = ħ k Here ħ is Planck’s constant divided by 2π. It’s the conversion factor between frequency and energy, and between wavenumber and momentum. A plane wave with angular frequency ω and wavevector k is ψ(x,t) = A exp(i(k·x − ωt)) Now watch what derivatives do to this wave: ∂ψ/∂t = −i ω ψ ∇ψ = i k ψ ∇²ψ = −|k|² ψ Multiply by ħ and you get: i ħ ∂ψ/∂t = ħ ω ψ = E ψ −i ħ ∇ψ = ħ k ψ = p ψ −ħ² ∇²ψ = ħ² |k|² ψ = p² ψ So for plane waves, the operators Ê = i ħ ∂/∂t p̂ = −i ħ ∇ act like energy and momentum. Now bring in the classical, nonrelativistic energy bookkeeping: E = p²/(2m) + V(x) Kinetic plus potential. That’s it. Turn it into an equation for ψ by replacing E and p with the operators above: Ê ψ = (p̂²/(2m) + V) ψ Since p̂² = (−i ħ ∇)·(−i ħ ∇) = −ħ² ∇², this becomes i ħ ∂ψ/∂t = ( −ħ²/(2m) ∇² + V(x) ) ψ That’s the time-dependent Schrödinger equation. This derivation is a controlled heuristic. Match the plane-wave identities to the measured relations E = ħω and p = ħk, then impose the same energy bookkeeping you trust in classical mechanics. Why this is the right kind of rule If ψ is the state, we need a rule that preserves total probability: ∫ |ψ(x,t)|² dx = 1 Schrödinger evolution does, and you can see it by deriving a continuity equation. Let ρ(x,t) = |ψ|² = ψ*ψ. Differentiate: ∂ρ/∂t = ψ* ∂ψ/∂t + ψ ∂ψ*/∂t Use Schrödinger and its complex conjugate. The potential terms cancel, and what’s left can be rearranged into ∂ρ/∂t + ∇·j = 0 with probability current j = (ħ/(2mi)) ( ψ* ∇ψ − ψ ∇ψ* ) That’s the cleanest way to say what ψ is. |ψ|² behaves like a conserved density, the phase drives a current, and the time evolution is fixed, up to V, by combining wave relations with energy bookkeeping: i ħ ∂ψ/∂t = ( −ħ²/(2m) ∇² + V ) ψ #QuantumMechanics #SchrodingerEquation #WaveFunction #BornRule #Physics #MathematicalPhysics

Mathelirium

20,781 views • 6 months ago

A physical system, its phase-space state, and the energy function that governs the flow, all linked together. The point of the animation is that these are not three different objects. They are three views of the same dynamics. The pendulum on the left is the physical motion. The curve on the floor is the evolution of its state. The translucent surface shows the energy landscape that organizes that motion.
0:30

A physical system, its phase-space state, and the energy function that governs the flow, all linked together. The point of the animation is that these are not three different objects. They are three views of the same dynamics. The pendulum on the left is the physical motion. The curve on the floor is the evolution of its state. The translucent surface shows the energy landscape that organizes that motion.

Mathelirium

19,865 views • 3 months ago

We just watched Professor Arthur Mattuck kick off MIT’s ODE course with the one interpretation that most differential equations classes somehow postpone: An ordinary differential equation isn’t primarily a method hunt. It’s a geometric rule. You write dy/dx = f(x,y), and that right-hand side is literally telling you the slope your solution curve must have at each point (x,y). So I built this animation as a visual companion to that first lecture. It draws the direction field (little line elements whose slope is f(x,y)) and then shows integral curves sliding through it...curves that are tangent to the field everywhere they go. Two quick examples from the animation: For dy/dx = −x/y, the slope field steers you onto circles x² + y² = R². You also see a subtle point that gets missed when everything is taught as y(x)...even when the curve exists smoothly, the graph y(x) may only exist on a limited interval (|x|<R for the upper semicircle). For dy/dx = 1 + x − y, the isoclines (curves where the slope is constant) make the global behavior obvious...trajectories get funneled into a corridor and become asymptotic to the special solution y=x. You learn qualitative behavior without solving it the traditional way. #DifferentialEquations #ODEs #MITOCW #VectorFields #MathAnimation #Mathematics
0:30

We just watched Professor Arthur Mattuck kick off MIT’s ODE course with the one interpretation that most differential equations classes somehow postpone: An ordinary differential equation isn’t primarily a method hunt. It’s a geometric rule. You write dy/dx = f(x,y), and that right-hand side is literally telling you the slope your solution curve must have at each point (x,y). So I built this animation as a visual companion to that first lecture. It draws the direction field (little line elements whose slope is f(x,y)) and then shows integral curves sliding through it...curves that are tangent to the field everywhere they go. Two quick examples from the animation: For dy/dx = −x/y, the slope field steers you onto circles x² + y² = R². You also see a subtle point that gets missed when everything is taught as y(x)...even when the curve exists smoothly, the graph y(x) may only exist on a limited interval (|x|<R for the upper semicircle). For dy/dx = 1 + x − y, the isoclines (curves where the slope is constant) make the global behavior obvious...trajectories get funneled into a corridor and become asymptotic to the special solution y=x. You learn qualitative behavior without solving it the traditional way. #DifferentialEquations #ODEs #MITOCW #VectorFields #MathAnimation #Mathematics

Mathelirium

22,440 views • 7 months ago

Can Spacetime Move Faster Than Light? In 1994, Miguel Alcubierre proposed a strange solution to Einstein’s field equations: A Spacetime geometry where a compact region could be transported without the object inside locally breaking the speed of light. The idea is not that a ship accelerates through space faster than light. The local interior remains comparatively calm. The violent part is the surrounding geometry. Space contracts in front of the bubble, expands behind it, and forms a moving transport shell that carries the interior region forward. A simplified form of the Alcubierre metric is ds² = −dt² + [dx − vₛ f(rₛ)dt]² + dy² + dz² where f(rₛ) shapes the warp bubble and vₛ describes the motion of its centre. The physics is still highly speculative, especially because the original construction requires exotic stress-energy, but the mathematical idea is powerful. #GeneralRelativity #AlcubierreDrive #WarpDrive #Spacetime #EinsteinEquations #Relativity #Physics #Astrophysics #Cosmology
0:56

Can Spacetime Move Faster Than Light? In 1994, Miguel Alcubierre proposed a strange solution to Einstein’s field equations: A Spacetime geometry where a compact region could be transported without the object inside locally breaking the speed of light. The idea is not that a ship accelerates through space faster than light. The local interior remains comparatively calm. The violent part is the surrounding geometry. Space contracts in front of the bubble, expands behind it, and forms a moving transport shell that carries the interior region forward. A simplified form of the Alcubierre metric is ds² = −dt² + [dx − vₛ f(rₛ)dt]² + dy² + dz² where f(rₛ) shapes the warp bubble and vₛ describes the motion of its centre. The physics is still highly speculative, especially because the original construction requires exotic stress-energy, but the mathematical idea is powerful. #GeneralRelativity #AlcubierreDrive #WarpDrive #Spacetime #EinsteinEquations #Relativity #Physics #Astrophysics #Cosmology

Mathelirium

37,986 views • 1 month ago

Primes Numbers Are Not Random Noise Prime numbers look scattered when you meet them one by one, but their disorder has structure. In this animation, every integer is placed on a growing spiral, primes ignite as sharp golden events, and the hidden oscillations come from the first Riemann zeta zeros. The scene is based on the von Mangoldt prime signal Λ(n) and the explicit formula for ψ(x), where primes appear as impulses and the zeta zeros behave like frequencies correcting the smooth drift of x. #PrimeNumbers #RiemannHypothesis #NumberTheory #Mathematics #MathAnimation
0:32

Primes Numbers Are Not Random Noise Prime numbers look scattered when you meet them one by one, but their disorder has structure. In this animation, every integer is placed on a growing spiral, primes ignite as sharp golden events, and the hidden oscillations come from the first Riemann zeta zeros. The scene is based on the von Mangoldt prime signal Λ(n) and the explicit formula for ψ(x), where primes appear as impulses and the zeta zeros behave like frequencies correcting the smooth drift of x. #PrimeNumbers #RiemannHypothesis #NumberTheory #Mathematics #MathAnimation

Mathelirium

16,018 views • 2 months ago