En esta tesis estudiaremos diversos aspectos del entrelazamiento entre diferentes particiones de sistemas de muchos cuerpos, tanto para el estado fundamental como para estados térmicos. Consideraremos en general aquellos sistemas definidos por Hamiltonianos cuadráticos en operadores locales y en particular el caso de redes de espines. Con este fin se desarrollará en primer lugar una generalización de la aproximación denominada CSPA de la función de partición de sistemas cuánticos compuestos, así como del formalismo RPA, al que se reduce el primero en situaciones donde el sistema admite soluciones de campo medio estables. A partir de la función de partición aproximada se estimará la concurrencia de pares para redes de espines 1/2 invariantes ante translaciones. Se mostrará también que en redes de espín 1/2 invariantes traslacionales la función de partición en la aproximación RPA puede evaluarse en forma analítica. En segundo lugar se presentará un mapeo bosónico consistente con RPA, que nos permitirá mapear el estado fundamental de un sistema general a un estado en el sistema de bosones que en primera aproximación pertenece a la clase de los estados Gaussianos. Esto nos permitirá estimar el entrelazamiento entre cualquier par de partes del sistema original en términos del entrelazamiento en estados Gaussianos, que puede ser evaluado exactamente. Luego se mostrarán resultados concretos en redes generales de espín s invariantes ante transformaciones de “paridad de espín z”. Veremos además como este tratamiento nos permite determinar las condiciones exactas para la existencia de un campo factorizante del sistema, esto es, configuraciones particulares del campo magnético externo para las que el sistema admite un (o en general varios) estado(s) fundamental(es) completamente separables. En los ejemplos estudiados, mostraremos cómo en el caso en que estos estados presenten ruptura espontánea de simetría el sistema desarrolla correlaciones cuánticas de largo alcance en la vecindad de ese punto, independientemente del alcance de las interacciones entre las partes. Finalmente se analizará en detalle el comportamiento exacto del sistema en esas condiciones. / In this work we will study several aspects of the entanglement of different partitions of a many body system, both for the ground state as well as for the thermal equilibrium states. We will consider in general those defined by Hamiltonians which are quadratic in local observables and in particular, the case of spin arrays. With this aim, we will first develop a generalization of the CSPA method for evaluating the partition function of a general composite quantum system, as well as of the RPA formalism, to which reduces the former in situations where the system admits stable mean field solutions. From these approximations to the partition function, the pairwise entanglement in spin 1/2 translational-invariant arrays will be estimated. We will also show that in this kind of systems the RPA partition function can be evaluated in a fully analytical way. Secondly, we will show a bosonic map consistent with RPA, which allows us to map the ground state of a general composite quantum system to a boson state which, in first approximation, belongs to the class of Gaussian states. This allows us to estimate the entanglement between any pair of parts of the original system in terms of the entanglement of Gaussian-states, which can be evaluated exactly. Later we will show explicit results for general, “Z Spin-Parity” invariant spin s arrays. We will also see that this treatment allows us to determine the exact conditions for the existence of a factorizing field in the system, i.e., particular configurations of the external magnetic field for which the system admits a (in general, several) fully separable ground state(s). In the analyzed examples, we show how the system develops, when these states present a spontaneous breaking of a Hamiltonian symmetry, long range quantum correlations in its immediate vicinity, independently of the interaction range. Finally, we will discuss in detail the exact behavior of the system in these conditions.