Abstract An Approximation of Löwdin's Natural Orbitals for Molecules with a Green's Function Method The many-body-pertubation theorie of the single-particle Green's function is used to get an approximate first-order density matrix. Slightly modified SCF-orbitals form the basis for the expansion. The mass-operator in Dyson's equation is considered up to second order in the Perturbation. In the present form the method is only applicable to ground states with closed shells. The ground states of the molecules LiH and NH3 serve as examples to demonstrate the usefulness of the directly calculated natural orbitals for application in the C I-method. The natural orbitals give a much better convergence of the C I-expansion than the SCF-orbitals do.