In a closed system, the total number of particles is fixed. We ask how much does this conservation law restricts the amount of entanglement that can be created. We derive a tight upper bound on the bipartite entanglement entropy in closed systems, and find what a maximally entangled state looks like in such a system. Finally, we illustrate numerically on an isolated system of one-dimensional fermionic gas, that the upper bound can be reached during its unitary evolution, when starting in a pure state that emulates a thermal state with high enough temperature. These results are in accordance with current experiments measuring Rényi-2 entanglement entropy, all of which employ particle-conserving Hamiltonian, where our bound acts as a loose bound, and will become especially important for bounding the amount of entanglement that can be spontaneously created, once a direct measurement of entanglement entropy becomes feasible.Entanglement is one of the most intriguing characteristics of quantum systems. It evolved from its perception as a mathematical artifact, as a result of EPR paradox [1], to becoming closely related and applicable to the fields of condensed matter [2-7], quantum information [8][9][10][11][12][13][14], quantum metrology [15][16][17][18][19][20], and quantum gravity [21][22][23][24][25].In the field of quantum information, entangled states are the backbone of quantum information protocols as they are considered a resource for tasks such as quantum teleportation [9,26], cryptography [8], and dense coding [27].In these quantum information protocols, more entanglement usually leads to a better performance. Therefore, it is important to set precise upper bounds on how much entanglement is in principle available in performing these tasks [28][29][30][31][32][33][34][35][36][37].As different tasks require different types of entangled states, numerous measures of entanglement have been introduced [38][39][40][41]. An important measure of entanglement is entanglement entropy [33,34,42]. It is defined as the von Neumann entropy of the reduced density matrix ρ A = tr B [ρ], whereρ denotes the density matrix of the composite system,This is a valuable measure as it draws a direct connection between density matrix and the amount of nonlocal correlations present in a given system. Entanglement entropy also gained significant attention in the past few decades due to the discovery of its geometric scaling in thermal state as well as ground states (famously known as the volume law [43] and the area law [44-46] respectively), and its use for characterizing quantum phase transition [2,[47][48][49].Despite its importance, this quantity has proven extremely difficult to probe experimentally, and related Rényi-2 entanglement entropy has been measured instead [50][51][52]. However, an experimental proposal has been put forward recently [53], opening new exciting possibilities.There exists a general bound on entanglement entropy. For a pure state of a bipartite system, it is straight forward to show that S ent ≡ S(ρ A ) = S...