Abstract-Artificial boundary conditions, which can be identified as Robin boundary conditions positioned at a complex space coordinate, are introduced in order to obtain pertinent approximations for the Green's functions in grounded layered media. These artificial boundary conditions include perfectly matched layers backed by perfectly electric or magnetic conductors. As a first result, we obtain analytical expressions for the differences of Green's functions subject to different boundary conditions. Since weighted sums of Green's functions are again Green's functions, the need arises to solve an optimization problem, in the sense of obtaining the optimal weighted mixture of Green's functions, as compared to the exact Green's function. Comprehensive eigenexpansions for the Green's functions are given in the general case, and a few examples illustrate the goodness of fit between the approximate Green's functions and the exact Green's function.