Electrostatic Interaction of Point Charges in Three-Layer Structures: The Classical Model

• Published in: Condensed matter Vol. 4; no. 2; p. 44
• Main Authors: Gabovich, Alexander M., Voitenko, Alexander I.
• Format: Journal Article
• Language: English
• Published: Basel MDPI AG 01.06.2019
• Subjects: Approximation; Configurations; Electron-hole interaction; electrostatic interaction; Energy; Excitons; Graphene; Green’s function method; Heterostructures; Interfaces; Mathematical analysis; Phase transitions; polarization forces; three-layer structures
• DOI: 10.3390/condmat4020044
• ISSN: 2410-3896

Electrostatic interaction energy W between two point charges in a three-layer plane system was calculated on the basis of the Green’s function method in the classical model of constant dielectric permittivities for all media involved. A regular method for the calculation of W ( Z , Z ′ , R ) , where Z and Z ′ are the charge coordinates normal to the interfaces, and R the lateral (along the interfaces) distance between the charges, was proposed. The method consists in substituting the evaluation of integrals of rapidly oscillating functions over the semi-infinite interval by constructing an analytical series of inverse radical functions to a required accuracy. Simple finite-term analytical approximations of the dependence W ( Z , Z ′ , R ) were proposed. Two especially important particular cases of charge configurations were analyzed in more detail: (i) both charges are in the same medium and Z = Z ′ ; and (ii) the charges are located at different interfaces across the slab. It was demonstrated that the W dependence on the charge–charge distance S = R 2 + Z − Z ′ 2 differs from the classical Coulombic one W ∼ S − 1 . This phenomenon occurs due to the appearance of polarization charges at both interfaces, which ascribes a many-body character to the problem from the outset. The results obtained testify, in particular, that the electron–hole interaction in heterostructures leading to the exciton formation is different in the intra-slab and across-slab charge configurations, which is usually overlooked in specific calculations related to the subject concerned. Our consideration clearly demonstrates the origin, the character, and the consequences of the actual difference. The often used Rytova–Keldysh approximation was analyzed. The cause of its relative success was explained, and the applicability limits were determined.