How to determine which p-type GTO to choose in HF method? I am trying to understand the working principle of basis sets in quantum chemistry. As far as I understand in the Hartree-Fock method, the basis set is used to calculate the matrix representation of Fock operator in each iteration in the SCF procedure. Usually each basis function in a chosen basis set is represented as a linear combination of several primitive Gaussian type orbital (GTO). I want to use cc-PVDZ basis for C atom which I can take from the EMSL Basis Set Exchange Library. For example in the Gaussian format:


!  cc-pVDZ  EMSL  Basis Set Exchange Library   10/19/17 2:40 AM
! Elements                             References
! --------                             ----------
! H     : T.H. Dunning, Jr. J. Chem. Phys. 90, 1007 (1989).
! He    : D.E. Woon and T.H. Dunning, Jr. J. Chem. Phys. 100, 2975 (1994).
! Li - Ne: T.H. Dunning, Jr. J. Chem. Phys. 90, 1007 (1989).
! Na - Mg: D.E. Woon and T.H. Dunning, Jr.  (to be published)
! Al - Ar: D.E. Woon and T.H. Dunning, Jr.  J. Chem. Phys. 98, 1358 (1993).
! Sc - Zn: N.B. Balabanov and K.A. Peterson, J. Chem. Phys. 123, 064107 (2005),
! N.B. Balabanov and K.A. Peterson, J. Chem. Phys. 125, 074110 (2006)
! Ca     : J. Koput and K.A. Peterson, J. Phys. Chem. A, 106, 9595 (2002).
! 



****
C     0 
S   8   1.00
   6665.0000000              0.0006920        
   1000.0000000              0.0053290        
    228.0000000              0.0270770        
     64.7100000              0.1017180        
     21.0600000              0.2747400        
      7.4950000              0.4485640        
      2.7970000              0.2850740        
      0.5215000              0.0152040        
S   8   1.00
   6665.0000000             -0.0001460        
   1000.0000000             -0.0011540        
    228.0000000             -0.0057250        
     64.7100000             -0.0233120        
     21.0600000             -0.0639550        
      7.4950000             -0.1499810        
      2.7970000             -0.1272620        
      0.5215000              0.5445290        
S   1   1.00
      0.1596000              1.0000000        
P   3   1.00
      9.4390000              0.0381090        
      2.0020000              0.2094800        
      0.5456000              0.5085570        
P   1   1.00
      0.1517000              1.0000000        
D   1   1.00
      0.5500000              1.0000000        
****



In the lines such as S   8   1.00 I know that the letter S is the orbital type, 8 is the number of GTO in this orbital's expansion, but what is 1.00 for?
The GTO for a p-type orbital has the form $x^l y^m z^n \exp(-ar^2)$ where the powers satisfy $l+m+n=1$ and thus there are three possible sets of $\{l,m,n\}$. But for instance in that link for lines with p-type orbitals, which possibilities shall I pick? The same question for d-type orbitals.
Have a look in the Gaussian manual if something is unclear.

If you want to use an extern basis set in Gaussian you have to respect the gaussian input conventions for the basis set, which you can find here (gaussian.com/gen/). Since you are using the EMSL Basis Set Exchange the server gives you a file that satisfies these conventions. Like:

Type  NGauss  Sc
$α_1\quad d_{1μ}$
$α_2\quad d_{2μ}$
…
$α_N\quad d_{Nμ}$  

Type stands for the shell, then you have the number of primitives with NGauss and scaling factor Sc. Afterwards you put your exponents and coeffients of the primitive gaussian functions which then will be contracted to the respective contracted gaussian functions. You might notice that the scaling factor is set to unity for almost all basis sets—I checked some Dunning, Pople and Ahlrichs basis sets and always found it set to unity, which seems resonable. In the TURBOMOLE format the scaling factor is omitted totally, so it seems to be a Gaussian specific hack.

Now to your second question. You are specifing the shell, not a single function, a shell contains always all spherical harmonics belonging to the given quantum number or all cartesian gaussian functions. The spherical harmonics are the solution of the Schrödinger equation for an electron on a surface. You have generally $2l+1$ spherical harmonics where $l$ is the azimudal quantum number specified as 0/s, 1/p, 2/d, 3/f, 4/g, 5/h and so on. There exists also the possibilty to use cartesian gaussians with $(l+1)(l+2)/2$ degenerated functions.

Why am I telling you that? You always choose a whole shell with a azimudal quantum number, how many function belong to this shell depends on the quantum chemistry code you are using. ASAIK, Gaussian supports both.

In case of a p-shell it does not matter, both expressions evaluate to three p-functions (as expected). So you get three p-functions with one p-shell.

The basis set you have chosen would be evaluted as:

(9s4p1d) → [3s2p1d]

Having nine primitive s type gaussians (note that Dunning basis sets use a general contraction scheme), four primitive p type gaussians and one primitive d type gaussian which form a basis set out of three s shells, two p shells and one d shell with 14 contracted spherical harmonic gaussian functions or 15 contracted cartesian gaussian functions for carbon.