Entanglement scaling and criticality of quantum many-body systems in canonical quantization picture using tensor network

Simulating strongly-correlated quantum systems in continuous space belongs to the most challenging and long-concerned issues in quantum physics. This work investigates the quantum entanglement and criticality of the ground-state wave-functions of infinitely-many coupled quantum oscillators (iCQOs). The essential task involves solving a set of partial differential equations (Schrödinger equations in the canonical quantization picture) with infinitely-many variables, which currently lacks valid methods. By extending the imaginary-time evolution algorithm with translationally-invariant functional tensor network, we simulate the ground state of iCQOs with the presence of two- and three-body couplings. We determine the range of coupling strengths where there exists a real ground-state energy (dubbed as physical region). With two-body couplings, we reveal the logarithmic scaling law of entanglement entropy (EE) and the polynomial scaling law of correlation length against the virtual bond dimension $\chi$ at the dividing point of physical and non-physical regions. These two scaling behaviors are signatures of criticality, according to the previous results in quantum lattice models, but were not reported in continuous-space quantum systems. The scaling coefficients result in a central charge $c=1$, indicating the presence of free boson conformal field theory (CFT). With three-body couplings where there exist no analytical nor numerical results, we show the breakdown of CFT at the dividing point even with an extremely small strength of three-body terms. Our work uncovers the scaling behaviors of EE in the continuous-space quantum many-body systems, providing numerical evidences for the efficiency of tensor networks in representing the continuous-space quantum many-body ground states in the thermodynamic limit.