Abstract: Optical and electronic properties of black phosphorus strongly depend on the number of layers and type of stacking. Using first-principles calculations within the framework of density functional theory, we investigate the electronic properties of bilayer black phosphorus with an interlayer twist angle of 90 degrees. These calculations are complemented with a simple (k) over right arrow . (p) over right arrow model which is able to capture most of the low energy features and is valid for arbitrary twist angles. The electronic spectrum of 90 degrees twisted bilayer black phosphorus is found to be x-y isotropic in contrast to the monolayer. However x-y anisotropy, and a partial return to monolayer-like behavior, particularly in the valence band, can be induced by an external out-of-plane electric field. Moreover, the preferred hole effective mass can be rotated by 90 degrees simply by changing the direction of the applied electric field. In particular, a +0.4 (-0.4) V angstrom(1) out-of-plane electric field results in a similar to 60% increase in the hole effective mass along the y (x) axis and enhances the m(y)*/m(x)* (m(x)*/m(y)*) ratio as much as by a factor of 40. Our DFT and (k) over right arrow . (p) over right arrow simulations clearly indicate that the twist angle in combination with an appropriate gate voltage is a novel way to tune the electronic and optical properties of bilayer phosphorus and it gives us a new degree of freedom to engineer the properties of black phosphorus based devices.

Abstract: Quantized magnetic vortices driven by electric current determine key electromagnetic properties of superconductors. While the dynamic behavior of slow vortices has been thoroughly investigated, the physics of ultrafast vortices under strong currents remains largely unexplored. Here, we use a nanoscale scanning superconducting quantum interference device to image vortices penetrating into a superconducting Pb film at rates of tens of GHz and moving with velocities of up to tens of km/s, which are not only much larger than the speed of sound but also exceed the pair-breaking speed limit of superconducting condensate. These experiments reveal formation of mesoscopic vortex channels which undergo cascades of bifurcations as the current and magnetic field increase. Our numerical simulations predict metamorphosis of fast Abrikosov vortices into mixed Abrikosov-Josephson vortices at even higher velocities. This work offers an insight into the fundamental physics of dynamic vortex states of superconductors at high current densities, crucial for many applications.

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: Starting from first principles, we show the formation and evolution of superconducting gaps in MgB2 at its ultrathin limit. Atomically thin MgB2 is distinctly different from bulk MgB2 in that surface states become comparable in electronic density to the bulklike sigma and pi bands. Combining the ab initio electron-phonon coupling with the anisotropic Eliashberg equations, we showthat monolayer MgB2 develops three distinct superconducting gaps, on completely separate parts of the Fermi surface due to the emergent surface contribution. These gaps hybridize nontrivially with every extra monolayer added to the film owing to the opening of additional coupling channels. Furthermore, we reveal that the three-gap superconductivity in monolayer MgB2 is robust over the entire temperature range that stretches up to a considerably high critical temperature of 20 K. The latter can be boosted to >50K under biaxial tensile strain of similar to 4%, which is an enhancement that is stronger than in any other graphene-related superconductor known to date.

Abstract: By making use of circularly polarized light and electrostatic gating, monolayer molybdenum disulfide (ML – MoS2) can form a platform supporting multiple types of charge carriers. They can be discriminated by their spin, valley index, or whether they are electrons or holes. We investigate the collective properties of those charge carriers and are able to identify distinct plasmon modes. We analyze the corresponding dispersion relation, lifetime, and oscillator strength, and calculate the phase relation between the oscillations in the different components of the plasmon modes. All platforms in ML-MoS2 support a long-wavelength root q plasmon branch at zero kelvins. In addition to this, for an n-component system, n-1 distinct plasmon modes appear as acoustic modes with linear dispersion in the long-wavelength limit. These modes correspond to out-of-phase oscillations in the different fermion liquids and have, although being damped, a relatively long lifetime. Additionally, we also find distinct modes at large wave vectors that are more strongly damped by intraband processes.

Abstract: Starting from a Hamiltonian with anharmonic coupling between in-plane acoustic displacements and outof-plane (flexural) modes, we derived coupled equations of motion for in-plane displacements correlations and flexural mode density fluctuations. Linear response theory and time-dependent thermal Green's functions techniques are applied in order to obtain different response functions. As external perturbations we allow for stresses and thermal heat sources. The displacement correlations are described by a Dyson equation where the flexural density distribution enters as an additional perturbation. The flexural density distribution satisfies a kinetic equation where the in-plane lattice displacements act as a perturbation. In the hydrodynamic limit this system of coupled equations is at the basis of a unified description of elastic and thermal phenomena, such as isothermal versus adiabatic sound motion and thermal conductivity versus second sound. The general theory is formulated in view of application to graphene, two-dimensional h-BN, and 2H-transition metal dichalcogenides and oxides.

Abstract: We propose a tight-binding (TB) model, that includes spin-orbit coupling (SOC), to describe the electronic properties of methyl-substituted germanane (GeCH3). This model gives an electronic spectrum in agreement with first principle results close to the Fermi level. Using the Z(2) formalism, we show that a topological phase transition from a normal insulator (NI) to a quantum spin Hall (QSH) phase occurs at 11.6% biaxial tensile strain. The sensitivity of the electronic properties of this system on strain, in particular its transition to the topological insulating phase, makes it very attractive for applications in strain sensors and other microelectronic applications.

Abstract: The Wigner-Liouville equation is reformulated using a spectral decomposition of the classical force field instead of the potential energy. The latter is shown to simplify the Wigner-Liouville kernel both conceptually and numerically as the spectral force Wigner-Liouville equation avoids the numerical evaluation of the highly oscillatory Wigner kernel which is nonlocal in both position and momentum. The quantum mechanical evolution is instead governed by a term local in space and non-local in momentum, where the non locality in momentum has only a limited range. An interpretation of the time evolution in terms of two processes is presented; a classical evolution under the influence of the averaged driving field, and a probability-preserving quantum-mechanical generation and annihilation term. Using the inherent stability and reduced complexity, a direct deterministic numerical implementation using Chebyshev and Fourier pseudo-spectral methods is detailed. For the purpose of illustration, we present results for the time evolution of a one-dimensional resonant tunneling diode driven out of equilibrium. (C) 2017 Elsevier Inc. All rights reserved.

Abstract: Using the tight-binding approach, we investigate the electronic properties of bilayer phosphorene (BLP) quantum dots (QDs) in the presence of perpendicular electric and magnetic fields. Since BLP consists of two coupled phosphorene layers, it is of interest to examine the layer-dependent electronic properties of BLP QDs, such as the electronic distributions over the two layers and the so-produced layer-polarization features, and to see how these properties are affected by the magnetic field and the bias potential. We find that in the absence of a bias potential only edge states are layer polarized while the bulk states are not, and the layer-polarization degree (LPD) of the unbiased edge states increases with increasing magnetic field. However, in the presence of a bias potential both the edge and bulk states are layer polarized, and the LPD of the bulk (edge) states depends strongly (weakly) on the interplay of the bias potential and the interlayer coupling. At high magnetic fields, applying a bias potential renders the bulk electrons in a BLP QD to be mainly distributed over the top or bottom layer, resulting in layer-polarized bulk Landau levels (LLs). In the presence of a large bias potential that can drive a semiconductor-to-semimetal transition in BLP, these bulk LLs exhibit different magnetic-field dependences, i.e., the zeroth LLs exhibit a linearlike dependence on the magnetic field while the other LLs exhibit a square-root-like dependence.