Larsen's Comment is based on the erroneous assumption that the second-order term from expansions of approximate variational calculations in powers of the coupling constant should agree with direct second-order perturbative calculations.
Direct evaluation of Feynman's functional integral for the polaron maps it into a standard type of many-body theory. The perturbation expansion is shown to be well controlled and does not require any renormalization. This perturbation expansion has the pleasing feature that nth-order terms in the polaron coupling constant are obtained from nth-order perturbation theory, whereas the traditional perturbation formulation in terms of creation and annihilation operators requires a 2nth-order perturbation calculation to generate the same results. In contrast to previous beliefs, the method is shown to enable the calculation of the polaron effective mass and relaxation time directly from the functional-integral expression. The ground-state energy is computed to second order in the polaron coupling constant to illustrate the use of the method and to establish agreement with previous perturbation treatments. In addition, we provide the first perturbation calculation of the polaron effective mass through second order in the coupling constant. Pad\'e approximants are shown to provide a vehicle for analytically extending the perturbation results into the region of intermediate and strong couplings, and the procedure yields reasonable agreement with variational calculations for values of the polaron coupling constant up to about $\ensuremath{\alpha}=10$. In addition, conditions are derived under which the original polaron action, used in the functional integral, reduces to some of the best variational trial actions utilized in the literature.
A refined approach where the wave function is gradually improved is developed in the one-dimensional (1D) large-Frohlich-polaron problem. The results for the ground-state energy, the effective mass, and the average number of virtual phonons are obtained for a wide range of the coupling constant It, and are in agreement with those of the Feynman path-integral formalism. In addition, in the weak-coupling limit, the expansions of these observables in powers of the coupling constant are exactly calculated up to the alpha 3 -term. They are E=- alpha -0.06066 alpha 2 -0.00844 alpha 3 , m*=1+0.5 alpha +0.19194 alpha 2 +0.06912 alpha 3 , and N=0.5 alpha +0.12132 alpha 2 +0.029540 alpha 3 , which agree with the known results of the fourth-order perturbation theory.