2025-12-25 | | Total: 6
Quantum algorithms to integrate nonlinear PDEs governing flow problems are challenging to discover but critical to enhancing the practical usefulness of quantum computing. We present here a near-optimal, robust, and end-to-end quantum algorithm to solve time-dependent, dissipative, and nonlinear PDEs. We embed the PDEs in a truncated, high dimensional linear space on the basis of quantum homotopy analysis. The linearized system is discretized and integrated using finite-difference methods that use a compact quantum algorithm. The present approach can adapt its input to the nature of nonlinearity and underlying physics. The complexity estimates improve existing approaches in terms of scaling of matrix operator norms, condition number, simulation time, and accuracy. We provide a general embedding strategy, bounds on stability criteria, accuracy, gate counts and query complexity. A physically motivated measure of nonlinearity is connected to a parameter that is similar to the flow Reynolds number $Re_{\textrm{H}}$, whose inverse marks the allowed integration window, for given accuracy and complexity. We illustrate the embedding scheme with numerical simulations of a one-dimensional Burgers problem. This work shows the potential of the hybrid quantum algorithm for simulating practical and nonlinear phenomena on near-term and fault-tolerant quantum devices.
Within the framework of many-body perturbation theory based on Green's functions, the $GW$ approximation has emerged as a pivotal method for computing quasiparticle energies and excitation spectra. However, its high computational cost and steep scaling present significant challenges for applications to large molecular systems. In this work, we extend the block tensor decomposition (BTD) algorithm, recently developed in our previous work [J. Chem. Phys. 163, 174109 (2025)] for low-rank tensor compression, to enable a formally $O(N^3)$-scaling $GW$ algorithm. By integrating BTD with an imaginary-time $GW$ formalism and introducing a real space screening strategy for the polarizability, we achieve an observed scaling of approximately $O(N^2)$ in test systems. Key parameters of the algorithm are optimized on the S66 dataset using the JADE algorithm, ensuring a balanced compromise between accuracy and efficiency. Our BTD-based random phase approximation also exhibits $O(N^2)$ scaling, and eigenvalue-self-consistent $GW$ calculations become feasible for systems with over 3000 basis functions. This work establishes BTD as an efficient and scalable approach for large-scale $GW$ calculations in molecular systems.
We study the emergence of a velocity-dip phenomenon in turbulent mixed convection in open Poiseuille-Rayleigh-Bénard (PRB) channels with a free-slip upper boundary. Three-dimensional direct numerical simulations (DNS) are performed for Rayleigh numbers in the range $10^5 \leq Ra \leq 10^8$, at a fixed Prandtl number $Pr = 0.71$ and a bulk Reynolds number $Re_b = 2850$. In the shear-dominated regime, the flow is characterised by small-scale structures such as near-wall streaks. As buoyancy becomes comparable to shear, streamwise-oriented large-scale rolls emerge and span the full channel height. At higher Rayleigh numbers, buoyancy dominates and the rolls fragment, giving rise to a convection-cell-dominated regime. Short-time-averaged flow fields show that streamwise rolls transport low-speed fluid from the bottom wall towards the upper boundary, forming laterally extended low-speed regions, while roll fragmentation induces upstream low-speed regions near the upper boundary. Both mechanisms locally reduce the near-surface mean velocity, leading to a velocity dip in which the maximum mean streamwise velocity is located below the upper boundary. Consistent with the mean momentum budget, the near-surface region exhibits a large-scale Reynolds shear stress that exceeds the local total shear stress, implying a negative viscous contribution and a reversal of the mean velocity gradient. To model this behaviour, we propose a model based on a balance between buoyancy and shear production with dissipation, incorporating a linear wall-normal profile for the Reynolds shear stress, a wall-normal-independent buoyancy-production term, and a decomposition of the dissipation into shear-induced and buoyancy-induced contributions. Our model accurately reproduces the DNS mean velocity profiles across the explored $Ra$ range.
An ab initio approach is presented for studying the collective excitations in excitonic insulators, charge/spin density waves and superconductors. We derive the Bethe-Salpeter-Equation for the particle-hole excitations in the quasiparticle representation, from which the collective excited states are solved and the corresponding order parameter fluctuations are computed. This method is demonstrated numerically for the excitonic insulating phases of the biased WSe2-MoSe2 bilayer. It reveals the gapless phase-mode, the subgap Bardasis-Schrieffer modes and the above-gap scattering states. Our work paves the way for quantitative predictions of excited state phenomena from first-principles calculations in electronic systems with spontaneous symmetry breaking.
The n body problem, fundamental to astrophysics, simulates the motion of n bodies acting under the effect of their own mutual gravitational interactions. Traditional machine learning models that are used for predicting and forecasting trajectories are often data intensive black box models, which ignore the physical laws, thereby lacking interpretability. Whereas Scientific Machine Learning ( Scientific ML ) directly embeds the known physical laws into the machine learning framework. Through robust modelling in the Julia programming language, our method uses the Scientific ML frameworks: Neural ordinary differential equations (NODEs) and Universal differential equations (UDEs) to predict and forecast the system dynamics. In addition, an essential component of our analysis involves determining the forecasting breakdown point, which is the smallest possible amount of training data our models need to predict future, unseen data accurately. We employ synthetically created noisy data to simulate real-world observational limitations. Our findings indicate that the UDE model is much more data efficient, needing only 20% of data for a correct forecast, whereas the Neural ODE requires 90%.
Flow Gym is a toolkit for research and deployment of flow-field quantification methods inspired by OpenAI Gym and Stable-Baselines3. It uses SynthPix as synthetic image generation engine and provides a unified interface for the testing, deployment and training of (learning-based) algorithms for flow-field quantification from a number of consecutive images of tracer particles. It also contains a growing number of integrations of existing algorithms and stable (re-)implementations in JAX.