Abstract: The determination of the layer number ${\cal N}$ in nanoscale thin layered crystals is a challenging problem of technological relevance. In addition to innovative experimental techniques, a thorough knowledge of the underlying lattice dynamics is required. Starting from phenomenological atomic interaction potentials we have carried out an analytical study of the low-frequency optical phonon dispersions in layered crystals. At the gamma point of the two-dimensional Brillouin zone the optical phonon frequencies correspond to rigid-plane shearing and compression modes. We have investigated graphene multilayers (GML) and hexagonal boron-nitride multilayers (BNML). The frequencies show a characteristic dependence on ${\cal N}$. The results which are represented in the form of fan diagrams are very similar for both materials. Due to charge neutrality within layers Coulomb forces play no role, only van der Waals forces between nearest neighbor layers are relevant. The theoretical results agree with recent low-frequency Raman results on rigid-layer modes [Tan et al., Nature Mater. 11, 294 (2012)] in GML and double-resonant Raman scattering data on rigid-layer compression modes [Herziger et al., Phys. Rev. B 85, 235447 (2012)] in GML. (C) 2012 WILEY-VCH Verlag GmbH Co. KGaA, Weinheim

Abstract: Starting from an empirical force constant model of valence interactions and calculating by Ewald's method the electrostatic force constants, we derive the dynamical matrix for a monolayer and for multilayer systems of hexagonal boron nitride (h-BN). Solution of the secular problem leads to the corresponding phonon dispersion relations. The interplay between valence forces and Coulomb forces is discussed. A comparison with previous results on graphene and graphene multilayers is made. Our spectra on the h-BN monolayer are rather similar to previous ab initio theory results. Comparison is also made with Raman and infrared experimental results.

Abstract: Born's long wave method is used to study the elastic properties of graphite and graphene. Starting from an empirical force constant model derived from full inplane phonon dispersions of graphite [Mohr et al., Phys. Rev. B 76, 035439 (2007)] we calculate the tension coefficients of graphene. Extending the model by interplanar interactions, we calculate the elastic constants of graphite. The agreement of our theoretical values with inelastic x-ray scattering results on elastic constants of graphite [Bosak et al., Phys. Rev. B 75, 153408 (2007)] is very satisfactory.

Abstract: The thermal expansion alpha(T) in layered crystals is of fundamental and technological interest. As suggested by I. M. Lifshitz in 1952, in thin solid films (crystalline membranes) a negative contribution to alpha(T) is due to anharmonic couplings between in-plane stretching modes and out-of-plane bending (flexural modes). Genuine in-plane anharmonicities give a positive contribution to alpha(T). The competition between these two effects can lead to a change of sign (crossover) from a negative value of alpha(T) in a temperature (T) range T <= T-alpha to a positive value of alpha(T) for T > T-alpha in layered crystals. Here, we present an analytical lattice dynamical theory of these phenomena for a two-dimensional (2D) hexagonal crystal. We start from a Hamiltonian that comprises anharmonic terms of third and fourth order in the lattice displacements. The in-plane and out-of-plane contributions to the thermal expansion are studied as functions of T for crystals of different sizes. Besides, renormalization of the flexural mode frequencies plays a crucial role in determining the crossover temperature T-alpha. Numerical examples are given for graphene where the anharmonic couplings are determined from experiments. The theory is applicable to other layer crystals wherever the anharmonic couplings are known. (C) 2015 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim

Abstract: Piezo- and flexoelectricity are manifestations of electromechanical coupling in solids with potential applications in nanoscale materials. Naumov etal. [Phys. Rev. Lett. 102, 217601 (2009)] have shown by first principles calculations that a monolayer BN sheet becomes macroscopically polarized in-plane when in a corrugated state. Here, we investigate the interplay of layer corrugation and in-plane polarization by atomistic lattice dynamics. We treat the coupling between static flexural modes and in-plane atomic ion displacements as an anharmonic effect, similar to the membrane effect that is at the origin of negative thermal expansion in layered crystals. We have derived analytical expressions for the corrugation-induced static in-plane strains and the optical displacements with the resulting polarization response functions. Beyond h-BN, the theory applies to transition metal dichalcogenides and dioxides. Numerical calculations show that the effects are considerably stronger for 2D h-BN than for 2H-MoS2.

Abstract: First‐principles calculations, within the framework of density functional theory, are used to investigate the structural, electronic, optical, and thermoelectric properties of monolayer PbTe. The effect of layer thickness, electric field, strain, and vacancy defects on the electronic and magnetic properties is systematically studied. The results show that the bandgap decreases as the layer thickness increases from monolayer to bulk. With application of an electric field on bilayer PbTe, the bandgap decreases from 70 meV (0.2 V Å⁻¹) to 50 meV (1 V Å⁻¹) when including spin–orbit coupling (SOC). Application of uniaxial strain induces a direct‐to‐indirect bandgap transition for strain greater than +6%. In addition, the bandgap decreases under compressive biaxial strain (with SOC). The effect of vacancy defects on the electronic properties of PbTe is also investigated. Such vacancy defects turn PbTe into a ferromagnetic metal (single vacancy Pb) with a magnetic moment of 1.3 μB, and into an indirect semiconductor with bandgap of 1.2 eV (single Te vacancy) and 1.5 eV (double Pb + Te vacancy). In addition, with change of the Te vacancy concentration, a bandgap of 0.38 eV (5.55%), 0.43 eV (8.33%), and 0.46 eV (11.11%) is predicted.

Abstract: We study the effect of mechanical deformations on the elec- tronic properties of hexagonal flakes of bilayer graphene. The behavior of electrons induced by triaxial strain can be de- scribed by an effective pseudo-magnetic field which is homo- geneous in the center of the flake. We find that in-plane strain, applied to both layers equally, can break the layer symmetry leading to different behavior in the top and bottom layers of graphene. At low energy, just one of the layers feels
the pseudo-magnetic field: the zero-energy pseudo-Landau level is missing in the second layer, thus creating a gap be- tween the lowest non-zero levels. While the layer asymmetry is most significant at zero energy, interaction with the edges of the flake extends the effect to higher pseudo-Landau lev- els. The behavior of the top and bottom layers may be re- versed by rotating the triaxial strain by 60°.

Abstract: The DFT-1/2 method is a band-gap correction with GW precision at a density functional theory (DFT) computational cost. The method was also extended to correct the gap between defect levels, allowing for the calculation of optical transitions. However, this method fails when the atomic character of the occupied and unoccupied defect levels is similar as we illustrate by two examples, the tetrahedral hydrogen interstitial and the negatively charged vacancy in diamond. We solve this problem by decoupling the effect of the occupied and unoccupied defect levels and call this the decoupled DFT-1/2 method for defects.

Abstract: We propose a simple quantum mechanical model describing the time dependent diffusion current between two fermion reservoirs that were initially disconnected and characterized by different densities or chemical potentials. The exact, analytical solution of the model yields the transient behavior of the coupled fermion systems evolving to a final steady state, whereas the long-time behavior is determined by a power law rather than by exponential decay. Similar results are obtained for the entropy production which is proportional to the diffusion current. (C) 2014 Elsevier B.V. All rights reserved.

Abstract: Knowing the exact number of particles N, and taking this knowledge into account, the quantum canonical ensemble imposes a constraint on the occupation number operators. The constraint particularly hampers the systematic calculation of the partition function and any relevant thermodynamic expectation value for arbitrary but fixed N. On the other hand, fixing only the average number of particles, one may remove the above constraint and simply factorize the traces in Fock space into traces over single-particle states. As is well known, that would be the strategy of the grand-canonical ensemble which, however, comes with an additional Lagrange multiplier to impose the average number of particles. The appearance of this multiplier can be avoided by invoking a projection operator that enables a constraint-free computation of the partition function and its derived quantities in the canonical ensemble, at the price of an angular or contour integration. Introduced in the recent past to handle various issues related to particle-number projected statistics, the projection operator approach proves beneficial to a wide variety of problems in condensed matter physics for which the canonical ensemble offers a natural and appropriate environment. In this light, we present a systematic treatment of the canonical ensemble that embeds the projection operator into the formalism of second quantization while explicitly fixing N, the very number of particles rather than the average. Being applicable to both bosonic and fermionic systems in arbitrary dimensions, transparent integral representations are provided for the partition function Z(N) and the Helmholtz free energy F-N as well as for two- and four-point correlation functions. The chemical potential is not a Lagrange multiplier regulating the average particle number but can be extracted from FN+1 – F-N, as illustrated for a two-dimensional fermion gas. (C) 2017 Elsevier B.V. All rights reserved.

Abstract: Recently, we have used a projection operator to fix the number of particles in a second quantization approach in order to deal with the canonical ensemble. Having been applied earlier to handle various problems in nuclear physics that involve fixed particle numbers, the projector formalism was extended to grant access as well to quantum-statistical averages in condensed matter physics, such as particle densities and correlation functions. In this light, the occupation numbers of the subsequent single-particle energy eigenstates are key quantities to be examined. The goal of this paper is (1) to provide a sound extension of the projector formalism directly addressing the occupation numbers as well as the chemical potential, and (2) to demonstrate how the emerging problems related to numerical instability for fermions can be resolved to obtain the canonical statistical quantities for both fermions and bosons. (C) 2018 Elsevier B.V. All rights reserved.

Abstract: Cyclotron resonance (CR) of the electrons accumulated at sheets with heavy Si doping in InSb were observed using far infrared radiation. The angular dependence of the CR follows closely the 1/cos theta behaviour with some small deviations at high angles attributed to coupling between subbands. From the effective mass of the lowest subband, which is found to be 0.027m(o), the bottom of the lowest subband was determined to lie 125 meV below the Fermi level.