A Luneburg lens bends an Electromagnetic wave without a curved boundary. Instead its dielectric constant changes continuously from the centre to the edge. This forces different parts of the wavelength to travel at different speeds. Here, we solve the full 2D TM EM wave equation on an unstructured triangular mesh using the Finite Element Method (FEM).

Einstein Once Attempted to Build a Static Universe But It Could Not Stay Still Einstein once modified his own field equations because he believed the universe should remain perfectly static forever. The extra term he introduced, the cosmological constant Λ, balanced gravity against cosmic expansion. However, even this balance was unstable, and later observations of distant galaxies showed the universe was expanding, not holding still.

General Relativity Keeps Revealing Unexpected Spacetime Geometries That Go Far Beyond Everyday Intuition. In 1988, Kip Thorne and Michael Morris described a traversable Wormhole geometry that strangely connects two distant regions of Spacetime without an event horizon or central singularity. No traversable Wormhole has ever been observed. Within General Relativity, though, these solutions are mathematically valid if Spacetime is supported by exotic stress-energy that violates classical energy conditions. It's a solution Einstein’s field equations allow, but nature has not confirmed. #Einstein# #GeneralRelativity# #Wormholes# #Spacetime# #Physics# #Relativity#

What if Your Neural Network Was Forced to Obey Physics? Physics-Informed Neural Networks (PINNs) are neural networks trained to satisfy a differential equation by building the PDE residual directly into the loss. They emerged from a very practical problem...classical PDE pipelines can be brilliant, but they often demand heavy discretization work (meshes, stencils, stability tuning), and the method you build is usually tied to one geometry and one solver setup. A PINN flips the workflow by representing the solution itself as a smooth function uᵩ(x,t) and enforcing the physics everywhere you choose to sample the domain. People often meet PINNs in the least helpful way...via a flashy solution plot, and almost no explanation of what was enforced to get it. In this series we keep the enforcement visible. We pick a differential equation, represent the unknown solution as a flexible function, measure how well that function satisfies the equation across the domain, and train it to reduce that mismatch everywhere we sample. A normal neural net learns from labels...you give it inputs and target outputs. A PINN learns from a differential equation...you give it inputs (x,t) and it gets punished whenever its output fails the PDE. By punish we mean that the loss increases when the mismatch is large we reward it if the loss decreases as the mismatch gets smaller. The network isn’t replacing physics, it’s becoming a flexible function that is forced to satisfy the same calculus you’d impose on any candidate solution. The math breakdown: We start with a PDE we want to solve on a domain Ω. Write it as uₜ(x,t) + N(u(x,t), uₓ(x,t), uₓₓ(x,t), …) = 0 for (x,t) in Ω A PINN replaces the unknown function u with a neural network output uᵩ(x,t) Now define the physics residual by plugging uᵩ into the PDE rᵩ(x,t) = ∂uᵩ/∂t + N(uᵩ, ∂uᵩ/∂x, ∂²uᵩ/∂x², …) If uᵩ were an exact solution, we would have rᵩ(x,t) = 0 everywhere. We may also have data points (xᵢ,tᵢ,uᵢ) from measurements or a known initial condition. The training objective is just a weighted sum of squared errors L(ᵩ) = L_data(ᵩ) + λ L_phys(ᵩ) + L_bc/ic(ᵩ) with L_data(ᵩ) = meanᵢ |uᵩ(xᵢ,tᵢ) − uᵢ|² L_phys(ᵩ) = meanⱼ |rᵩ(xⱼ,tⱼ)|² where (xⱼ,tⱼ) are the collocation points in Ω L_bc/ic(ᵩ) = penalties enforcing boundary conditions and initial conditions The key technical step is that the derivatives inside rᵩ are computed by automatic differentiation ∂uᵩ/∂t, ∂uᵩ/∂x, ∂²uᵩ/∂x², … So we can differentiate the total loss L(ᵩ) with respect to ᵩ and train with gradient descent. This is the whole idea behind PINNs. Learn a function, but make the PDE part of the loss, so the network is trained to be a solution, not just a curve-fitter. In the render, the main 3D surface is the network’s current guess uᵩ(x,t), drawn as a living sheet over the (x,t) plane. Hovering above is the neural scaffold...a visible graph of feature nodes and connections. The bright tension threads are the physics residual rᵩ(x,t): each thread tethers a collocation bead on the sheet up to the scaffold, and it thickens and brightens exactly where |rᵩ| is large (color encodes the sign). As training runs, those threads go slack across the domain not because we hid the error, but because the network has actually been pushed toward rᵩ(x,t) ≈ 0. #PINNs# #PhysicsInformedNeuralNetworks# #ScientificMachineLearning# #PDE# #DifferentialEquations# #Optimization# #MachineLearning# #AppliedMath# #ComputationalPhysics#

Electromagnetic Y-Splitter (Finite Element Simulation) The field naturally splits into two waveguide branches, with energy distribution set by geometry and boundary conditions. Used in integrated photonics, optical communication routing, and interferometric sensing systems where controlled beam splitting without active components is required.

How Do You Measure Motion Smaller Than an Atom? The Michelson Interferometer measures incredibly tiny distance changes by converting mirror motion into shifting interference fringes. That sensitivity makes it useful in gravitational-wave detection, seismic sensing, vibration analysis, optical calibration, and other high-precision measurements.

Oscillatory Quadrupole Distortions A direct visualization of the h₊ and hₓ polarization modes of gravitational waves propagating through curved spacetime around a black hole.

Quantum Turbulence Can Organize Into Crystals. Take the same underlying Schrödinger-Poisson system as the previous two scenes, but now the wavefield is periodically driven and mechanically stirred, forcing quantized vortices to nucleate, interact, and eventually self-organize into ordered crystal-like lattices. The bright cyan and amber defects are genuine phase singularities with quantized circulation. What looks chaotic at first slowly develops long-range structure, as Floquet forcing and self-gravity push the condensate toward coherent vortex ordering #QuantumPhysics# #QuantumFluid# #VortexDynamics# #SchrodingerEquation# #ComputationalPhysics# #Mathematics# #Physics#

Quantum Matter Can Collapse Into Stellar-Like Structures. Take the same underlying system as the previous phase-helicoid scene we just posted, but now we stop looking at phase geometry and focus directly on how the density evolves under self-gravity. We note that self-attracting quantum waves do not always spread out. Under Schrödinger-Poisson dynamics, the density begins to cluster into bright gravitational condensations, forming turbulent filaments, rotating cores, and star-like structures driven entirely by the wavefunction’s own gravity. Result looks less like particles moving through space and more like Spacetime teaching a quantum fluid how to organize itself. #QuantumPhysics# #WaveFunction# #SchrodingerEquation# #Astrophysics# #ComputationalPhysics# #Mathematics# #Physics#

The Schrödinger Equation Can Fold Phase Into Geometry. Density is only half the story. Here, the quantum phase itself twists into moving helicoidal ribbons, while self-gravity from the Schrödinger-Poisson coupling bends and compresses the wavefield into glowing caustics and vortex singularities. Tiny white pearls mark places where the phase becomes undefined topological defects drifting through a self-generated gravitational landscape. #QuantumPhysics# #WaveFunction# #SchrodingerEquation# #ComputationalPhysics# #ScientificVisualization# #Mathematics#

The High-Dimensional Ball Is Full, but Its Probability Mass Lives Near the Boundary High-dimensional probability is weird. As the dimension grows, a uniform point in the ball does not usually land deep inside. Most of the mass gets pushed into a very thin shell near the boundary.