Entanglement scaling in critical two-dimensional fermionic and bosonic systems

We relate the reduced density matrices of quadratic bosonic and fermionic models to their Green's function matrices in a unified way and calculate the scaling of bipartite entanglement of finite systems in an infinite universe exactly. For critical fermionic 2D systems at T=0, two regimes of scaling are identified: generically, we find a logarithmic correction to the area law with a prefactor dependence on the chemical potential that confirms earlier predictions based on the Widom conjecture. If, however, the Fermi surface of the critical system is zero-dimensional, we find an area law with a sublogarithmic correction. For a critical bosonic 2D array of coupled oscillators at T=0, our results show that entanglement follows the area law without corrections.