“…The dimensions of the hBN sheet were 30.05 and 30.37 Å in the x - and y -directions, respectively, and a vacuum space of 60 Å in the z -direction was used to eliminate periodic image interactions. The induced dipoles on the boron and nitrogen atoms were modeled using the classical Drude oscillator model with Thole damping, while the vibration of these atoms was modeled using the DREIDING force field. , The expressions for the induced dipole moments on the boron and nitrogen atoms in the classical Drude oscillator model are given by μ i ,B ind = q D,B d i ,B and μ i ,N ind = q D,B d i ,N , respectively, where d is the DC–DP separation distance. The variations of the induced dipole moments, per unit cell of hBN, resulting from the external electric field applied along the y -axis and the z -axis, as a function of the MD simulation time, are shown in Figure a,b, respectively.…”