The quantitative description of the quantum entanglement between a qubit and its environment is considered. Specifically, for the ground state of the spin-boson model, the entropy of entanglement of the spin is calculated as a function of α, the strength of the ohmic coupling to the environment, and ε, the level asymmetry. This is done by a numerical renormalization group treatment of the related anisotropic Kondo model. For ε = 0, the entanglement increases monotonically with α, until it becomes maximal for α → 1 − . For fixed ε > 0, the entanglement is a maximum as a function of α for a value, α = αM < 1.Due to the promise of quantum computation there is currently considerable interest in the relationship between entanglement, decoherence, entropy, and measurement. Motivated by quantum information theory several authors have recently investigated entanglement in quantum many-body systems [1,2,3,4]. It is often stated that decoherence or a measurement causes a system to become entangled with its environment. The purpose of this paper is to make these ideas quantitative by a study of the simplest possible model, the spin-boson model [5,6]. This describes a qubit (two-level system) interacting with an infinite collection of harmonic oscillators that model the environment responsible for decoherence and dissipation. Specifically, we show how the entanglement between a superposition state of the qubit and the environment changes as the coupling between the qubit and environment increases. One interesting result is that we find that the qubit becomes maximally entangled with the environment when the coupling α approaches a particular finite value (α → 1 − ). Furthermore, at this value the model undergoes a quantum phase transition, which is consistent with recent observations that often entanglement is largest near quantum critical points [1,2,3,4].The spin-boson model. The Hamiltonian is [5, 6]where ∆ is the bare tunneling amplitude between the two quantum mechanical states ↑ and ↓, ε is the level asymmetry (or bias), ω i are the frequencies of the oscillators and λ i the strength with which they couple to the two quantum mechanical states. The effect of the oscillator bath is completely determined by the spectral function * Electronic address: tac@tkm.physik.uni-karlsruhe.de † Electronic address: mckenzie@physics.uq.edu.au J(ω), defined below [5]. We will only consider the ohmic case, where it is has a linear dependence on frequencyfor ω ≪ ω c , and α is the dimensionless dissipation strength. The cutoff frequency, ω c ≫ ∆. This model can describe the decoherence of Josephson junction qubits, such as those recently realized experimentally [7], due to voltage fluctuations in the electronic circuit [8], and α can be expressed in terms of resistances and capacitances in the circuit and so this is an experimentally tunable parameter. Recent results show it is possible to construct devices, with α ≪ 1, the regime required for quantum computation. However, when modeling measurements one has α ∼ 1. The dynamical properti...