The set of Gaussian Type Orbitals g(n1, n2, n3) of order (n + 1)(n + 2)/2 and of common n ≡ n1 + n2 + n3 ≤ 7, common center and exponential, is customized to define a set of 2n + 1 linear combinations tn,m (−n ≤ m ≤ n) such that each tn,m depends on the azimuthal and polar angle of the spherical coordinate system like the real or imaginary part of the associated Spherical Harmonic. (Results cover both Hermite and Cartesian Gaussian Type Orbitals.) Overlap, kinetic energy and Coulomb energy matrix elements are presented for generalized basis functions of the type r s tn,m (s = 0, 2, 4 . . .). In addition, normalization integrals |g(n1, n2, n3)|d 3 r are calculated up to n = 7 and normalization integrals |r s tn,m|d 3 r up to n = 5. Gaussian Type Orbitals (GTO's) are widely used construction elements of basis sets in quantum chemistry. They became highly successful owing to fairly simple analytical representations of key integrals 1,2,3,4,5 that pay off in terms of speed when the matrix elements are calculated. This work deals with the conversion of Cartesian or Hermite GTO's which have a product representation in Cartesian coordinates into harmonic oscillator functions which have a product representation in spherical coordinates. The latter are orthogonal with respect to two quantum indices -the total number of integrals to be calculated is reduced by the overlap and kinetic energy integrals over products of orbitals with the same center.Sec. II introduces part of the notation. Sec. III lists those linear combinations of Gaussians defined in Cartesian coordinates that obtain the orthogonal property and are found to be the harmonic oscillator eigenfunctions. To extend coverage of the functional space, the oscillator functions are generalized in Sec. IV. Section V computes two-center overlap, kinetic energy and Coulomb integrals of those without recourse to the expansions in Cartesian coordinates. Sec. VI tabulates some of their products in support of 2-particle Coulomb Integrals. Absolute norms of orbitals play some role in Fitting Function techniques and are given in Sec. VII for the types of Gaussians discussed before.(1) be a primitive Hermite GTO (HGTO) with exponent α and quantum number n ≡ n 1 + n 2 + n 3 centered at the origin. H ni are Hermite Polynomials, and r 2 ≡ x 2 + y 2 + z 2 . It is normalized as 56(n + 1)(n + 2)/2 different g(n 1 , n 2 , n 3 ) exist for a given n. If n ≥ 2, they build an overcomplete set of states compared to only 2n + 1 eigenstates Y m n of the angular momentum operator.shall denote Spherical Harmonics in spherical coordinates, 78 andgeneralized Legendre Polynomials. Their real-valued counterparts areand(0 < m ≤ l), all normalized to unity, with parity (−) l , and orthogonal, The main result of this work is support to use of GTO's in systems with heavy or highly polarized atoms by reduction of the overcomplete sets to sets of 2n + 1 linearly independent combinations of GTO's.