Dispersion, static correlation, and delocalisation errors in density functional theory: An electrostatic theorem perspective J. Chem. Phys. 135, 164110 (2011) Evaluation of coupling terms between intra-and intermolecular vibrations in coarse-grained normal-mode analysis: Does a stronger acid make a stiffer hydrogen bond? J. Chem. Phys. 135, 154111 (2011) Calculating dispersion interactions using maximally localized Wannier functions J. Chem. Phys. 135, 154105 (2011) A theoretical study of Ne3 using hyperspherical coordinates and a slow variable discretization approach J. Chem. Phys. 135, 134312 (2011) Additional information on J. Chem. Phys. The frozen-density embedding ͑FDE͒ scheme ͓Wesolowski and Warshel, J. Phys. Chem. 97, 8050 ͑1993͔͒ relies on the use of approximations for the kinetic-energy component v T ͓ 1 , 2 ͔ of the embedding potential. While with approximations derived from generalized-gradient approximation kinetic-energy density functional weak interactions between subsystems such as hydrogen bonds can be described rather accurately, these approximations break down for bonds with a covalent character. Thus, to be able to directly apply the FDE scheme to subsystems connected by covalent bonds, improved approximations to v T are needed. As a first step toward this goal, we have implemented a method for the numerical calculation of accurate references for v T . We present accurate embedding potentials for a selected set of model systems, in which the subsystems are connected by hydrogen bonds of various strength ͑water dimer and F-H-F − ͒, a coordination bond ͑ammonia borane͒, and a prototypical covalent bond ͑ethane͒. These accurate potentials are analyzed and compared to those obtained from popular kinetic-energy density functionals.