A research team has developed a new quantum algorithm capable of solving nonlinear scalar conservation laws, a type of partial differential equation (PDE) fundamental in physics and engineering. This advance is significant because nonlinear dynamics are inherently challenging to simulate directly with unitary quantum algorithms. The proposed method transforms the nonlinear equation into a linear Liouville equation, which is then discretized and embedded into a unitary evolution through a process of "Schrödingerization."
The algorithm also includes quantum procedures for estimating relevant observables from the evolved state. Error bounds and gate-complexity estimates have been established for the complete algorithm. The resulting complexity comparison suggests a quantum advantage for observable estimation in sufficiently high spatial dimensions, under standard assumptions on state preparation and oracle access. This implies the algorithm could outperform classical methods in certain computational scenarios.
The strategy builds upon the level-set formulation, a technique that allows addressing the complexity of nonlinear equations. Numerical experiments validate the method's accuracy, its applicability to multidimensional problems, and its predicted scaling. This work opens new avenues for applying quantum computing to complex problems in fluid dynamics, plasma physics, and other areas where nonlinear conservation laws are crucial.