Loading video...

Video Failed to Load

Go Home

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...

mathelirium's profile picture

Mathelirium

39,450 subscribers

30,447 views • 4 months ago •via X (Twitter)

EducationScience & Technology
LIVE

Anya Rossi• Live Now

Private livecam show

0 Comments

No comments available

Comments from the original post will appear here

Related Videos

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 • 1 month 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

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,681 views • 4 months ago

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 ṗᵧ = -y - λ(x² - y²) The saddle energy is Eₛ = 1/(6λ²) Below that energy, the particle remains trapped. But a tiny packet of nearby initial conditions stretches into filaments, folds through the potential, and begins to look almost random. However, it is not random. Hamiltonian flow preserves the hidden phase-space structure. What looks like chaos in ordinary space is still a deterministic loom of energy conserved, equations obeyed, geometry stretched without being erased. #HamiltonianMechanics #ChaosTheory #PhaseSpace #DynamicalSystems #MathematicalPhysics #NonlinearDynamics #HenonHeiles #SymplecticGeometry #ScienceAnimation
0:58

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 ṗᵧ = -y - λ(x² - y²) The saddle energy is Eₛ = 1/(6λ²) Below that energy, the particle remains trapped. But a tiny packet of nearby initial conditions stretches into filaments, folds through the potential, and begins to look almost random. However, it is not random. Hamiltonian flow preserves the hidden phase-space structure. What looks like chaos in ordinary space is still a deterministic loom of energy conserved, equations obeyed, geometry stretched without being erased. #HamiltonianMechanics #ChaosTheory #PhaseSpace #DynamicalSystems #MathematicalPhysics #NonlinearDynamics #HenonHeiles #SymplecticGeometry #ScienceAnimation

Mathelirium

26,893 views • 2 months 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,651 views • 9 months ago

Visualization of the harmonic oscillator. Showing how x(t) = A cos(ωt + φ) is the horizontal projection of uniform circular motion. The spring–mass system and the closed ellipse in (x, v) phase space represent the same evolving state, while kinetic and elastic potential energy continuously exchange as the total energy remains constant.
0:45

Visualization of the harmonic oscillator. Showing how x(t) = A cos(ωt + φ) is the horizontal projection of uniform circular motion. The spring–mass system and the closed ellipse in (x, v) phase space represent the same evolving state, while kinetic and elastic potential energy continuously exchange as the total energy remains constant.

Philosophy Of Physics

10,497 views • 27 days ago

Do you actually know what convex optimization is in the geometric, guarantee-theoretic sense or have you only met it through solvers and loss curves? Convexity is rare comfort in optimization...there are no spurious local minima, no surprise traps, and inequalities you can use like tools instead of prayers. So, what is this convexity? Let x = (x₁, x₂) and let f(x) be convex. Plot the surface z = f(x). Pick a contact point x₀. The local slope is the gradient p = ∇f(x₀). That p is exactly the data that defines the supporting plane: z = f(x₀) + p · (x − x₀). Thus, f is said to be convex because for every x, f(x) ≥ f(x₀) + p · (x − x₀). So the plane at x₀ can slide under the surface, but it never slices through it. Not near the point...everywhere. Now for here is the interesting part: The slope becomes a coordinate system! Rewrite the same plane as z = p · x − b, where b is the offset. Because the plane passes through (x₀, f(x₀)), the offset is forced to be b = p · x₀ − f(x₀). And that number isn’t just geometry trivia. It’s the convex conjugate: f*(p) = sup over x ( p · x − f(x) ). At a differentiable contact point, the supporting plane touches f tightly enough that the supremum is achieved at x₀, giving the identity f*(p) = p · x₀ − f(x₀) when p = ∇f(x₀). So one moving contact point gives two linked readouts: primal position x₀ dual position (slope) p = ∇f(x₀) dual offset f*(p) One surface. Two worlds. #ConvexOptimization #Optimization #MachineLearning #SignalProcessing #AppliedMath #Engineering
0:50

Do you actually know what convex optimization is in the geometric, guarantee-theoretic sense or have you only met it through solvers and loss curves? Convexity is rare comfort in optimization...there are no spurious local minima, no surprise traps, and inequalities you can use like tools instead of prayers. So, what is this convexity? Let x = (x₁, x₂) and let f(x) be convex. Plot the surface z = f(x). Pick a contact point x₀. The local slope is the gradient p = ∇f(x₀). That p is exactly the data that defines the supporting plane: z = f(x₀) + p · (x − x₀). Thus, f is said to be convex because for every x, f(x) ≥ f(x₀) + p · (x − x₀). So the plane at x₀ can slide under the surface, but it never slices through it. Not near the point...everywhere. Now for here is the interesting part: The slope becomes a coordinate system! Rewrite the same plane as z = p · x − b, where b is the offset. Because the plane passes through (x₀, f(x₀)), the offset is forced to be b = p · x₀ − f(x₀). And that number isn’t just geometry trivia. It’s the convex conjugate: f*(p) = sup over x ( p · x − f(x) ). At a differentiable contact point, the supporting plane touches f tightly enough that the supremum is achieved at x₀, giving the identity f*(p) = p · x₀ − f(x₀) when p = ∇f(x₀). So one moving contact point gives two linked readouts: primal position x₀ dual position (slope) p = ∇f(x₀) dual offset f*(p) One surface. Two worlds. #ConvexOptimization #Optimization #MachineLearning #SignalProcessing #AppliedMath #Engineering

Mathelirium

38,506 views • 7 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

Quantum mechanics has a reputation for being mystical mainly because people skip the rules and jump to interpretations. In this lecture series, we’re doing the opposite. We start from the rules, follow the algebra, and let the picture be the calculation. Classical Probability Theory combines alternatives by adding their probabilities. Quantum Theory combines them one step earlier…add complex amplitudes first, then square at the end. That swap in order is everything. Expand |a₁ + a₂|² and you don’t just get |a₁|² + |a₂|²…you get a cross-term, 2 Re(a₁ a₂*). Its sign is set by phase, so the same two contributions can reinforce or cancel. Interference is just the algebra of squaring a sum. In the 3D render, the surface height is proportional to |a(x)| (so peaks become bright bands after squaring), while the surface skin is colored by the local phase arg(a(x)). As the phase knob φ(t) is swept on path 2, the cross-term oscillates, and you literally watch the interference ridges slide across the screen. We model a detector screen with coordinates x in R² (think x = (x,y)). A quantum state assigns a complex amplitude a(x). The rule for outcomes is p(x) = |a(x)|² Now the key situation: two coherent alternatives contribute to the same outcome x. Let their amplitudes be a₁(x) and a₂(x). Quantum says a(x) = a₁(x) + a₂(x) So the probability density becomes p(x) = |a₁(x) + a₂(x)|² Expand it (this is the whole episode): p(x) = (a₁ + a₂)(a₁* + a₂*) = |a₁|² + |a₂|² + a₁ a₂* + a₁* a₂ = |a₁|² + |a₂|² + 2 Re(a₁ a₂*) That last term is the interference term. It can be positive or negative. To see phase explicitly, write each contribution in polar form: a₁(x) = r₁(x) exp(i θ₁(x)) a₂(x) = r₂(x) exp(i θ₂(x)) Then a₁ a₂* = r₁ r₂ exp(i(θ₁ − θ₂)) So the cross-term is 2 Re(a₁ a₂*) = 2 r₁ r₂ cos(θ₁(x) − θ₂(x)) That’s the fringe engine: p(x) = r₁² + r₂² + 2 r₁ r₂ cos(Δθ(x)) Now the phase knob we animate: Add a controllable phase shift φ to path 2: a₂(x) → a₂(x) exp(i φ) Then Δθ(x) → Δθ(x) − φ, so p(x; φ) = r₁² + r₂² + 2 r₁ r₂ cos(Δθ(x) − φ) As φ changes smoothly, the bright/dark pattern slides continuously. Same setup, same geometry, same magnitudes r₁,r₂, only phase changed. #QuantumMechanics #WaveInterference #ComplexAmplitudes #DoubleSlit #Physics #Mathematics
1:00

Quantum mechanics has a reputation for being mystical mainly because people skip the rules and jump to interpretations. In this lecture series, we’re doing the opposite. We start from the rules, follow the algebra, and let the picture be the calculation. Classical Probability Theory combines alternatives by adding their probabilities. Quantum Theory combines them one step earlier…add complex amplitudes first, then square at the end. That swap in order is everything. Expand |a₁ + a₂|² and you don’t just get |a₁|² + |a₂|²…you get a cross-term, 2 Re(a₁ a₂*). Its sign is set by phase, so the same two contributions can reinforce or cancel. Interference is just the algebra of squaring a sum. In the 3D render, the surface height is proportional to |a(x)| (so peaks become bright bands after squaring), while the surface skin is colored by the local phase arg(a(x)). As the phase knob φ(t) is swept on path 2, the cross-term oscillates, and you literally watch the interference ridges slide across the screen. We model a detector screen with coordinates x in R² (think x = (x,y)). A quantum state assigns a complex amplitude a(x). The rule for outcomes is p(x) = |a(x)|² Now the key situation: two coherent alternatives contribute to the same outcome x. Let their amplitudes be a₁(x) and a₂(x). Quantum says a(x) = a₁(x) + a₂(x) So the probability density becomes p(x) = |a₁(x) + a₂(x)|² Expand it (this is the whole episode): p(x) = (a₁ + a₂)(a₁* + a₂*) = |a₁|² + |a₂|² + a₁ a₂* + a₁* a₂ = |a₁|² + |a₂|² + 2 Re(a₁ a₂*) That last term is the interference term. It can be positive or negative. To see phase explicitly, write each contribution in polar form: a₁(x) = r₁(x) exp(i θ₁(x)) a₂(x) = r₂(x) exp(i θ₂(x)) Then a₁ a₂* = r₁ r₂ exp(i(θ₁ − θ₂)) So the cross-term is 2 Re(a₁ a₂*) = 2 r₁ r₂ cos(θ₁(x) − θ₂(x)) That’s the fringe engine: p(x) = r₁² + r₂² + 2 r₁ r₂ cos(Δθ(x)) Now the phase knob we animate: Add a controllable phase shift φ to path 2: a₂(x) → a₂(x) exp(i φ) Then Δθ(x) → Δθ(x) − φ, so p(x; φ) = r₁² + r₂² + 2 r₁ r₂ cos(Δθ(x) − φ) As φ changes smoothly, the bright/dark pattern slides continuously. Same setup, same geometry, same magnitudes r₁,r₂, only phase changed. #QuantumMechanics #WaveInterference #ComplexAmplitudes #DoubleSlit #Physics #Mathematics

Mathelirium

81,608 views • 7 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

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 • 7 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 • 5 months ago

Lecture 3 of our Quantum Mechanics series. Lecture 2 gave us the one clean privilege quantum theory offers: treat ψ(x,t) as the state and ρ(x,t) = |ψ(x,t)|² as probability, because Schrödinger evolution forces ρ to obey a continuity equation. Lecture 3 is what that continuity equation is really telling you. If ρ behaves like a fluid, then the only question that matters is: What is the velocity field? Write ψ(x,t) = r(x,t) exp(i θ(x,t)). The magnitude r sets how much probability is sitting there. The phase θ sets where it tries to go. When you unpack the current j = Im(ψ* ∇ψ), it collapses to j = (ρ/m) ∇θ, which means the flow lines you draw are literally contours of phase geometry. Then the constraint that makes the picture bite: ψ has to be single-valued, so θ can’t wind by an arbitrary amount. Around any closed loop the total phase change must be 2π n, with n an integer. That’s why vortices aren’t features you add...they’re defects the math permits, in quantized units. In the render you see both layers at once...the 3D surface shows |ψ| breathing while the phase skin slides, and the 2D panel exposes the engine...current lines steering around discrete vortex charges. The math breakdown We write the state as a complex field ψ(x,t) on the plane (x in R²). The Born rule defines the probability density ρ(x,t) = |ψ(x,t)|² Schrödinger evolution (ħ = 1 units) is i ∂ψ/∂t = [ −(1/2m) ∇² + V(x,t) ] ψ Now derive conservation of probability. Start with ρ = ψ*ψ: ∂ρ/∂t = ψ* (∂ψ/∂t) + ψ (∂ψ*/∂t) Use Schrödinger and its complex conjugate: ∂ψ/∂t = (1/i) [ −(1/2m) ∇²ψ + Vψ ] ∂ψ*/∂t = (−1/i) [ −(1/2m) ∇²ψ* + Vψ* ] Substitute. The V terms cancel, and the remaining terms rearrange into the continuity equation ∂ρ/∂t + ∇·j = 0 with probability current j = (1/2mi) ( ψ* ∇ψ − ψ ∇ψ* ) = (1/m) Im(ψ* ∇ψ) So "probability density" really behaves like a conserved fluid density with flux j. Now expose the phase mechanism. Write ψ in polar form ψ(x,t) = r(x,t) exp(i θ(x,t)) Compute the gradient ∇ψ = exp(iθ) (∇r + i r ∇θ) Then ψ* ∇ψ = r (∇r + i r ∇θ) Taking the imaginary part gives Im(ψ* ∇ψ) = r² ∇θ = ρ ∇θ So the current becomes j = (ρ/m) ∇θ That’s the steering-wheel statement: Phase gradient sets the flow direction and speed (modulated by density and m). Finally, quantized vortices. Because ψ must be single-valued, going around any closed loop must return the same complex value. That forces the phase winding to be an integer multiple of 2π: ∮ ∇θ · dl = 2π n with n in Z n is the vortex charge. Vortex cores sit where ρ ≈ 0 (phase is undefined), and the current streamlines circulate around them. #QuantumMechanics #Wavefunction #SchrodingerEquation #BornRule #ProbabilityCurrent #ContinuityEquation #Phase #Vortices #TopologicalDefects #ComplexAnalysis #MathematicalPhysics #Mathematics #Physics
0:34

Lecture 3 of our Quantum Mechanics series. Lecture 2 gave us the one clean privilege quantum theory offers: treat ψ(x,t) as the state and ρ(x,t) = |ψ(x,t)|² as probability, because Schrödinger evolution forces ρ to obey a continuity equation. Lecture 3 is what that continuity equation is really telling you. If ρ behaves like a fluid, then the only question that matters is: What is the velocity field? Write ψ(x,t) = r(x,t) exp(i θ(x,t)). The magnitude r sets how much probability is sitting there. The phase θ sets where it tries to go. When you unpack the current j = Im(ψ* ∇ψ), it collapses to j = (ρ/m) ∇θ, which means the flow lines you draw are literally contours of phase geometry. Then the constraint that makes the picture bite: ψ has to be single-valued, so θ can’t wind by an arbitrary amount. Around any closed loop the total phase change must be 2π n, with n an integer. That’s why vortices aren’t features you add...they’re defects the math permits, in quantized units. In the render you see both layers at once...the 3D surface shows |ψ| breathing while the phase skin slides, and the 2D panel exposes the engine...current lines steering around discrete vortex charges. The math breakdown We write the state as a complex field ψ(x,t) on the plane (x in R²). The Born rule defines the probability density ρ(x,t) = |ψ(x,t)|² Schrödinger evolution (ħ = 1 units) is i ∂ψ/∂t = [ −(1/2m) ∇² + V(x,t) ] ψ Now derive conservation of probability. Start with ρ = ψ*ψ: ∂ρ/∂t = ψ* (∂ψ/∂t) + ψ (∂ψ*/∂t) Use Schrödinger and its complex conjugate: ∂ψ/∂t = (1/i) [ −(1/2m) ∇²ψ + Vψ ] ∂ψ*/∂t = (−1/i) [ −(1/2m) ∇²ψ* + Vψ* ] Substitute. The V terms cancel, and the remaining terms rearrange into the continuity equation ∂ρ/∂t + ∇·j = 0 with probability current j = (1/2mi) ( ψ* ∇ψ − ψ ∇ψ* ) = (1/m) Im(ψ* ∇ψ) So "probability density" really behaves like a conserved fluid density with flux j. Now expose the phase mechanism. Write ψ in polar form ψ(x,t) = r(x,t) exp(i θ(x,t)) Compute the gradient ∇ψ = exp(iθ) (∇r + i r ∇θ) Then ψ* ∇ψ = r (∇r + i r ∇θ) Taking the imaginary part gives Im(ψ* ∇ψ) = r² ∇θ = ρ ∇θ So the current becomes j = (ρ/m) ∇θ That’s the steering-wheel statement: Phase gradient sets the flow direction and speed (modulated by density and m). Finally, quantized vortices. Because ψ must be single-valued, going around any closed loop must return the same complex value. That forces the phase winding to be an integer multiple of 2π: ∮ ∇θ · dl = 2π n with n in Z n is the vortex charge. Vortex cores sit where ρ ≈ 0 (phase is undefined), and the current streamlines circulate around them. #QuantumMechanics #Wavefunction #SchrodingerEquation #BornRule #ProbabilityCurrent #ContinuityEquation #Phase #Vortices #TopologicalDefects #ComplexAnalysis #MathematicalPhysics #Mathematics #Physics

Mathelirium

37,998 views • 7 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

A quick explanation of the color you see in my space photography. Space is just as colorful as anything else, so why the misconceptions?
1:03

A quick explanation of the color you see in my space photography. Space is just as colorful as anything else, so why the misconceptions?

Andrew McCarthy

84,658 views • 2 years ago

Final Lecture of our Statistical Mechanics Series. Lecture 2 showed how we move from the Full Phase-Space Density ρ(q₁, …, qₙ, p₁, …, pₙ, t) to smaller statistical objects by integrating out variables we do not want to keep. That gives reduced descriptions like the One-Particle Density f₁(q₁,p₁,t) and the Two-Particle Density f₂(q₁,p₁,q₂,p₂,t) This was the simplification. Now comes the catch. If the Full Density obeys Liouville’s Equation, the reduced densities do not evolve independently. The equation for one level depends on the next one and this is referred to as the BBGKY hierarchy. The One-Particle Density depends on the Two-Particle Density. The Two-Particle Density depends on the Three-Particle Density. And the chain keeps going. That happens because particles interact. Once one particle feels the rest, one-particle information is no longer enough. Correlations enter, and the lower level is fed from above. If the full Hamiltonian is H = Σᵢ pᵢ²/(2m) + Σᵢ U(qᵢ) + (1/2) Σᵢ Σⱼ≠ᵢ Φ(qᵢ − qⱼ) then reducing the full density does not make the interaction terms disappear. It leaves behind coupling to higher-order reduced densities. So, schematically, ∂f₁/∂t + transport of one particle = interaction term involving f₂ and more generally the heirarchy is such that ∂fₛ/∂t + s-particle transport = interaction term involving fₛ₊₁ Therefore, Lecture 3 is really about the price of reduction. We simplify the description, but the information we remove comes back as coupling to higher-order correlations. So, how do you actually compute anything if every level depends on the next one? This is the so-called Closure Problem. To make the hierarchy usable, you need an extra assumption that cuts the chain. You replace the exact higher-order object by an approximation in terms of lower-order ones. The most basic example is a factorized closure at the pair level, where the exact correlated Two-Particle Density is replaced schematically by a product of One-Particle Densities: f₂(q₁,p₁,q₂,p₂,t) ≈ f₁(q₁,p₁,t) f₁(q₂,p₂,t) That approximation is not exact. It throws away part of the correlation structure. But it gives you something the raw hierarchy does not... a closed equation for the lower-level description. That is why closure matters so much. Without it, the hierarchy is exact but open. With it, the theory becomes approximate but usable. Thus, the combined point of this final Statistical Mechanics post is simple. First, reduced descriptions are not closed because interactions generate correlations across levels. Second, if you want a workable Kinetic Theory, you must close the hierarchy by approximating those higher-order correlations. It is the bridge from formal many-body mechanics to equations people can actually solve. In the render, that is exactly the story you are seeing. The first part shows the hierarchy itself: one reduced level feeding the next, with lower descriptions inheriting structure from higher ones. The second part shows the closure step where the exact correlated pair level is replaced by a factorized ansatz, and that approximation gives back a closed one-particle description. That is, the animation moves from dependence to approximation, and from approximation to solvability. #StatisticalMechanics #BBGKY #ClosureProblem #KineticTheory #PhaseSpace #ReducedDistribution #HamiltonianMechanics #MathematicalPhysics #Mathematics #Physics
1:35

Final Lecture of our Statistical Mechanics Series. Lecture 2 showed how we move from the Full Phase-Space Density ρ(q₁, …, qₙ, p₁, …, pₙ, t) to smaller statistical objects by integrating out variables we do not want to keep. That gives reduced descriptions like the One-Particle Density f₁(q₁,p₁,t) and the Two-Particle Density f₂(q₁,p₁,q₂,p₂,t) This was the simplification. Now comes the catch. If the Full Density obeys Liouville’s Equation, the reduced densities do not evolve independently. The equation for one level depends on the next one and this is referred to as the BBGKY hierarchy. The One-Particle Density depends on the Two-Particle Density. The Two-Particle Density depends on the Three-Particle Density. And the chain keeps going. That happens because particles interact. Once one particle feels the rest, one-particle information is no longer enough. Correlations enter, and the lower level is fed from above. If the full Hamiltonian is H = Σᵢ pᵢ²/(2m) + Σᵢ U(qᵢ) + (1/2) Σᵢ Σⱼ≠ᵢ Φ(qᵢ − qⱼ) then reducing the full density does not make the interaction terms disappear. It leaves behind coupling to higher-order reduced densities. So, schematically, ∂f₁/∂t + transport of one particle = interaction term involving f₂ and more generally the heirarchy is such that ∂fₛ/∂t + s-particle transport = interaction term involving fₛ₊₁ Therefore, Lecture 3 is really about the price of reduction. We simplify the description, but the information we remove comes back as coupling to higher-order correlations. So, how do you actually compute anything if every level depends on the next one? This is the so-called Closure Problem. To make the hierarchy usable, you need an extra assumption that cuts the chain. You replace the exact higher-order object by an approximation in terms of lower-order ones. The most basic example is a factorized closure at the pair level, where the exact correlated Two-Particle Density is replaced schematically by a product of One-Particle Densities: f₂(q₁,p₁,q₂,p₂,t) ≈ f₁(q₁,p₁,t) f₁(q₂,p₂,t) That approximation is not exact. It throws away part of the correlation structure. But it gives you something the raw hierarchy does not... a closed equation for the lower-level description. That is why closure matters so much. Without it, the hierarchy is exact but open. With it, the theory becomes approximate but usable. Thus, the combined point of this final Statistical Mechanics post is simple. First, reduced descriptions are not closed because interactions generate correlations across levels. Second, if you want a workable Kinetic Theory, you must close the hierarchy by approximating those higher-order correlations. It is the bridge from formal many-body mechanics to equations people can actually solve. In the render, that is exactly the story you are seeing. The first part shows the hierarchy itself: one reduced level feeding the next, with lower descriptions inheriting structure from higher ones. The second part shows the closure step where the exact correlated pair level is replaced by a factorized ansatz, and that approximation gives back a closed one-particle description. That is, the animation moves from dependence to approximation, and from approximation to solvability. #StatisticalMechanics #BBGKY #ClosureProblem #KineticTheory #PhaseSpace #ReducedDistribution #HamiltonianMechanics #MathematicalPhysics #Mathematics #Physics

Mathelirium

10,970 views • 4 months ago

Quantum Probability Surface Goes Negative A Schrödinger-Cat state is built from two separated coherent wavepackets, |ψ⟩=N(|α⟩+|−α⟩), and its Wigner function reveals the part of Quantum Mechanics that ordinary probability cannot show. The two bright lobes behave like classical possibilities, while the interference between them carves negative valleys into phase space, where W(x,p)<0. In this animation, the surface height is the signed Wigner function. So, gold rises above zero, cyan sinks below it, and the blades between them expose the nonclassical structure of the wavefunction. #QuantumMechanics #WignerFunction #SchrodingerCat #Physics #QuantumArt #ScienceVisualization
0:40

Quantum Probability Surface Goes Negative A Schrödinger-Cat state is built from two separated coherent wavepackets, |ψ⟩=N(|α⟩+|−α⟩), and its Wigner function reveals the part of Quantum Mechanics that ordinary probability cannot show. The two bright lobes behave like classical possibilities, while the interference between them carves negative valleys into phase space, where W(x,p)<0. In this animation, the surface height is the signed Wigner function. So, gold rises above zero, cyan sinks below it, and the blades between them expose the nonclassical structure of the wavefunction. #QuantumMechanics #WignerFunction #SchrodingerCat #Physics #QuantumArt #ScienceVisualization

Mathelirium

12,603 views • 2 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,841 views • 6 months ago

“When you look at the Earth from space, you realize that our planet is a beautiful, interconnected system. We are all in this together.” - Mae Jemison, Space Shuttle NASA astronaut This time-lapse taken aboard the International Space Station shows a blanket of stars surrounding our Earth. The view is one that is known to create a moment of pause in appreciation of the beauty of our shared home. #WorldViewWednesday
0:21

“When you look at the Earth from space, you realize that our planet is a beautiful, interconnected system. We are all in this together.” - Mae Jemison, Space Shuttle NASA astronaut This time-lapse taken aboard the International Space Station shows a blanket of stars surrounding our Earth. The view is one that is known to create a moment of pause in appreciation of the beauty of our shared home. #WorldViewWednesday

NASA's Johnson Space Center

61,946 views • 1 year ago

The Trap in Every Mathematics Lecture If you’ve taken a lot of math courses, you start to recognize a pattern. There’s a moment where the lecturer is warming up with the obvious stuff...add matrices entrywise, scale by α, do the row-column product...and you’re thinking, alright… where is this going? Then you relax. You stop resisting. And right there, they slip in one line that changes how you see the whole subject. When Benedict Gross says "matrices represent linear operators,"he’s telling you to stop treating a matrix as a rectangle of numbers and start treating it as an action. A linear operator is a function T: Rⁿ → Rⁿ that respects two rules: T(u+v)=T(u)+T(v) and T(αu)=αT(u). Once you pick a basis, T is completely determined by where it sends the basis vectors e₁,…,eₙ. Put T(e₁),…,T(eₙ) into columns and you get a matrix A. That is what "A represents T" means...A is the coordinate portrait of the transformation. Now the punchline that makes matrix multiplication feel inevitable. If B represents S and A represents T, then doing S first and then T is the composition T∘S. In coordinates that becomes A(Bx)=(AB)x. So multiplying matrices is really composing transformations. That’s why multiplication is usually not commutative: T∘S is generally not the same transformation as S∘T, and the matrices inherit that noncommutativity. This explains half of Linear Algebra because it tells you what the course is really about...functions that move vectors around, not grids of numbers. A matrix is just the written form of that function once you choose coordinates. Then the rules stop feeling random Multiplying matrices means doing one move and then another, an inverse means you can undo the move, eigenvectors are directions that don’t get turned, and changing basis is just describing the same move in a different language. That one idea makes a lot of linear algebra click. #LinearAlgebra #Matrices #GroupTheory #GLn #MathLectures #Mathematics
0:46

The Trap in Every Mathematics Lecture If you’ve taken a lot of math courses, you start to recognize a pattern. There’s a moment where the lecturer is warming up with the obvious stuff...add matrices entrywise, scale by α, do the row-column product...and you’re thinking, alright… where is this going? Then you relax. You stop resisting. And right there, they slip in one line that changes how you see the whole subject. When Benedict Gross says "matrices represent linear operators,"he’s telling you to stop treating a matrix as a rectangle of numbers and start treating it as an action. A linear operator is a function T: Rⁿ → Rⁿ that respects two rules: T(u+v)=T(u)+T(v) and T(αu)=αT(u). Once you pick a basis, T is completely determined by where it sends the basis vectors e₁,…,eₙ. Put T(e₁),…,T(eₙ) into columns and you get a matrix A. That is what "A represents T" means...A is the coordinate portrait of the transformation. Now the punchline that makes matrix multiplication feel inevitable. If B represents S and A represents T, then doing S first and then T is the composition T∘S. In coordinates that becomes A(Bx)=(AB)x. So multiplying matrices is really composing transformations. That’s why multiplication is usually not commutative: T∘S is generally not the same transformation as S∘T, and the matrices inherit that noncommutativity. This explains half of Linear Algebra because it tells you what the course is really about...functions that move vectors around, not grids of numbers. A matrix is just the written form of that function once you choose coordinates. Then the rules stop feeling random Multiplying matrices means doing one move and then another, an inverse means you can undo the move, eigenvectors are directions that don’t get turned, and changing basis is just describing the same move in a different language. That one idea makes a lot of linear algebra click. #LinearAlgebra #Matrices #GroupTheory #GLn #MathLectures #Mathematics

Mathelirium

66,892 views • 7 months ago