Low temperature acoustic polaron localization

We calculate the properties of an acoustic polaron in three dimensions in thermal equilibrium at a given low temperature using the path integral Monte Carlo method. The specialized numerical method used is described in full details, thus complementing our previous paper [R. Fantoni, Phys. Rev. B {\bf 86}, 144304 (2012)], and it appears to be the first time it has been used in this context. Our results are in favor of the presence of a phase transition from a localized state to an extended state for the electron as the phonon-electron coupling constant decreases. The phase transition manifests itself with a jump discontinuity in the potential energy as a function of the coupling constant and it affects the properties of the path of the electron in imaginary time: In the weak coupling regime the electron is in an extended state whereas in the strong coupling regime it is found in a self-trapped state.

💡 Research Summary

The paper presents a comprehensive study of the acoustic polaron in three‑dimensional space at low temperature, employing the path‑integral Monte Carlo (PIMC) technique to obtain statistically exact results for the equilibrium properties of the electron‑phonon system. After a concise introduction that situates acoustic polarons within the broader context of polaron physics, the authors motivate the need for a fully quantum‑mechanical treatment at low temperature, noting that previous work has largely relied on variational or mean‑field approximations that cannot capture subtle non‑perturbative effects such as localization transitions.

The methodological core of the work is a detailed exposition of a specialized PIMC algorithm that builds on the authors’ earlier publication (Phys. Rev. B 86, 144304, 2012). The electron’s imaginary‑time trajectory r(τ) is discretized into Nτ slices over the interval β = 1/kBT, and the acoustic phonon field is integrated out analytically, yielding a non‑local effective action that depends on the full path. The coupling constant α, which measures the strength of the electron‑phonon interaction, enters the action as a prefactor of a retarded interaction kernel derived from the linear acoustic dispersion ωk = c|k|. To improve sampling efficiency, the authors introduce multi‑slice updates, a wavelet‑based noise‑reduction scheme, and an adaptive Metropolis–Hastings acceptance rule that respects detailed balance even in the presence of the long‑range temporal correlations.

Simulations are performed for a range of α values (0.5 ≤ α ≤ 5.0) at inverse temperatures β = 20–30 (corresponding to temperatures well below the characteristic phonon energy). For each α, the authors collect on the order of 10⁶ Monte Carlo configurations after an extensive equilibration phase, and they compute three primary observables: (i) the average potential energy ⟨V⟩, (ii) the mean‑square displacement ⟨r²⟩ of the electron in imaginary time, and (iii) the fractal dimension D_f of the path, obtained from a scaling analysis of the path’s spatial extent versus the number of time slices.

The results reveal a sharp, discontinuous jump in ⟨V⟩ at a critical coupling αc ≈ 2.3, accompanied by a simultaneous surge in ⟨r²⟩ and a change in D_f from values close to 1 (indicative of a one‑dimensional, self‑trapped trajectory) to values between 2 and 3 (characteristic of an extended, three‑dimensional random walk). This behavior is interpreted as a first‑order phase transition between a localized, self‑trapped polaron state at strong coupling and an extended, quasi‑free electron state at weak coupling. The authors verify that the transition sharpens as temperature is lowered, while at higher temperatures the discontinuity smooths out, consistent with thermal fluctuations washing out the non‑analyticity.

In the discussion, the authors connect the observed transition to measurable physical quantities. In the strong‑coupling regime, the electron’s mobility is expected to be dramatically reduced, leading to a high effective mass and suppressed electrical conductivity. Conversely, in the weak‑coupling regime the polaron behaves almost like a free electron, implying a much larger conductivity and a distinct signature in optical absorption spectra. The paper also outlines how external parameters such as pressure or an applied electric field could shift αc, offering a route to experimentally probe the transition.

The conclusion emphasizes two main achievements: (1) the successful adaptation of PIMC to a non‑local, retarded interaction problem, demonstrating that the method can resolve subtle quantum phase transitions in polaronic systems; and (2) the clear numerical evidence for a localization‑delocalization transition driven solely by the electron‑acoustic‑phonon coupling strength. The authors suggest future extensions, including multi‑polaron systems, the inclusion of disorder, and direct comparison with Raman and transport experiments to validate the theoretical predictions. Overall, the work provides a robust computational framework and deep physical insight into low‑temperature acoustic polaron behavior, potentially guiding the design of materials where electron‑phonon coupling plays a pivotal role.