Abstract: Today, further downscaling of mobile electronic devices poses serious problems, such as energy consumption and local heat dissipation. In this context, spin wave majority gates made of very thin ferromagnetic films may offer a viable alternative. However, similar downscaling of magnetic thin films eventually enforces the latter to operate as quasi-2D magnets, the magnetic properties of which are not yet fully understood, especially those related to anisotropies and external magnetic fields in arbitrary directions. To this end, we have investigated the behaviour of an easy-plane and easy-axis anisotropic ferromagnet-both in two and three dimensions-subjected to a uniform magnetic field, applied along an arbitrary direction. In this paper, a spin-1/2 Heisenberg Hamiltonian with anisotropic exchange interactions is solved using double-time temperature-dependent Green's functions and the Tyablikov decoupling approximation. We determine various magnetic properties such as the Curie temperature and the magnetization as a function of temperature and the applied magnetic field, discussing the impact of the system's dimensionality and the type of anisotropy. The magnetic reorientation transition taking place in anisotropic Heisenberg ferromagnets is studied in detail. Importantly, spontaneous magnetization is found to be absent for easy-plane 2D spin systems with short range interactions.

Abstract: Transport properties of a two-dimensional electron gas (2DEG) are studied in the presence of a normal magnetic field B, of a weak one-dimensional (1D) periodic potential modulation V(x) = V(0)cos(Kx), and of the Rashba spin-orbit interaction (SOI) of strength a. For V(x) = 0 the SOI mixes the up and down spin states of neighboring Landau levels into two, unequally spaced energy branches. For V(x) not equal 0 these levels broaden into bands and their bandwidths oscillate with B. The n-th level bandwidth of each series vanishes at different values of B. Relative to the ID-modulated 2DEG without SOI and one flat-band condition, there are two flat-band conditions that depend on a and the transport coefficients can change considerably. For weak a the Weiss oscillations show beating patterns while for strong a the Shubnikov-de Haas ones axe split in two.

Abstract: Twisted graphene layers produce a moire pattern (MP) structure with a predetermined wavelength for a given twist angle. However, predicting the membrane corrugation amplitude for any angle other than pure AB-stacked or AA-stacked graphene is impossible using first-principles density functional theory (DFT) due to the large supercell. Here, within elasticity theory, we define the MP structure as the minimum-energy configuration, thereby leaving the height amplitude as the only unknown parameter. The latter is determined from DFT calculations for AB-and AA-stacked bilayer graphene in order to eliminate all fitting parameters. Excellent agreement with scanning tunneling microscopy results across multiple substrates is reported as a function of twist angle.

Abstract: Unidirectional chiral states are predicted in single layer graphene which originate from the breaking of the sublattice symmetry due to an asymmetric mass potential. The latter can be created experimentally using boron-nitride (BN) substrates with a line defect (B-B or N-N) that changes the induced mass potential in graphene. Solving the Dirac-Weyl equation, the obtained energy spectrum is compared with the one calculated using ab initio density functional calculations. We found that these one-dimensional chiral states are very robust and they can even exist in the presence of a small gap between the mass regions. In the latter case additional bound states are found that are topologically different from those chiral states.

Abstract: Using a particle based model, we investigate the skyrmion dynamical behavior in a channel where the upper wall contains divots of one depth and the lower wall contains divots of a different depth. Under an applied driving force, skyrmions in the channels move with a finite skyrmion Hall angle that deflects them toward the upper wall for -x direction driving and the lower wall for +x direction driving. When the upper divots have zero height, the skyrmions are deflected against the flat upper wall for -x direction driving and the skyrmion velocity depends linearly on the drive. For +x direction driving, the skyrmions are pushed against the lower divots and become trapped, giving reduced velocities and a nonlinear velocity-force response. When there are shallow divots on the upper wall and deep divots on the lower wall, skyrmions get trapped for both driving directions; however, due to the divot depth difference, skyrmions move more easily under -x direction driving, and become strongly trapped for +x direction driving. The preferred -x direction motion produces what we call a Magnus diode effect since it vanishes in the limit of zero Magnus force, unlike the diode effects observed for asymmetric sawtooth potentials. We show that the transport curves can exhibit a series of jumps or dips, negative differential conductivity, and reentrant pinning due to collective trapping events. We also discuss how our results relate to recent continuum modeling on a similar skyrmion diode system.

Abstract: Using an effective one-dimensional cluster expansion in combination with first-principles electronic structure calculations we have studied the energetics and electronic properties of Sb-Te layered systems. For a Te concentration between 0 and 60 at. % an almost continuous series of metastable structures is obtained consisting of consecutive Sb bilayers next to consecutive Sb2Te3 units, with the general formula (Sb-2)(n)(Sb2Te3)(m) (n, m = 1,2, ... ). Between 60 and 100 at.% no stable structures are found. We account explicitly for the weak van derWaals bonding between Sb bilayers and Sb2Te3 units by using a recently developed functional, which strongly improves the interlayer bonding distances. At T = 0 K, no evidence is found for the existence of two separate single-phase regions delta and gamma and a two-phase region delta + gamma. Metastable compounds with a Te concentration between 0 and 40 at. % are semimetallic, whereas compounds with a Te concentration between 50 and 60 at. % are semiconducting. Compounds with an odd number of Sb layers are metallic and have a much higher formation energy than those with an even number of consecutive Sb layers, thereby favoring the formation of Sb bilayers.

Abstract: Using an extension of the Gutzwiller approximation for an inhomogeneous system, we study the two-band Hubbard model with unequal band widths for a slab geometry. The aim is to investigate the mutual effect of individual bands on the spatial distribution of quasi-particle weight and charge density, especially near the surface of the slab. The main effect of the difference in band width is the presence of two different length scales corresponding to the quasi-particle profile of each band. This is enhanced in the vicinity of the critical interaction of the narrow band where an orbitally selective Mott transition occurs and a surface dead layer forms for the narrow band. For the doped case, two different regimes of charge transfer between the surface and the bulk of the slab are revealed. The charge transfer from surface/ center to center/ surface depends on both the doping level and the average relative charge accumulated in each band. Such effects could also be of importance when describing the accumulation of charges at the interface between structures made of multi-band strongly correlated materials.

Abstract: Using atomistic simulations we investigate the thermodynamical properties of a single atomic layer of hexagonal boron nitride (h-BN). The thermal induced ripples, heat capacity, and thermal lattice expansion of large scale h-BN sheets are determined and compared to those found for graphene (GE) for temperatures up to 1000 K. By analyzing the mean-square height fluctuations < h(2)> and the height-height correlation function H(q) we found that the h-BN sheet is a less stiff material as compared to graphene. The bending rigidity of h-BN (i) is about 16% smaller than the one of GE at room temperature (300 K), and (ii) increases with temperature as in GE. The difference in stiffness between h-BN and GE results in unequal responses to external uniaxial and shear stress and different buckling transitions. In contrast to a GE sheet, the buckling transition of a h-BN sheet depends strongly on the direction of the applied compression. The molar heat capacity, thermal-expansion coefficient, and Gruneisen parameter are estimated to be 25.2 J mol(-1) K-1, 7.2 x 10(-6) K-1, and 0.89, respectively.

Abstract: Using density functional theory, we study the electronic properties of a graphene-silicene bilayer (GSB). A single layer of silicene binds to the graphene layer with adhesion energy of about 25 meV/atom. This adhesion energy between the two layers follows accurately the well-known -1/z(2) dispersion energy as found between two infinite parallel plates. In small flakes of GSB with hydrogenated edges, negative charge is transferred from the graphene layer to the silicene layer, producing a permanent and a switchable polar bilayer, while in an infinite GSB, the negative charge is transferred from the silicene layer to the graphene layer. The graphene-silicene bilayer is a good candidate for a nanocapacitor with piezoelectric capabilities. We found that the permanent dipole of the bilayer can be tuned by an external perpendicular electric field.

Abstract: Using first-principles calculations and reactive force field molecular dynamics simulations, we study the structural properties and dynamics of a fluorine (F) atom, either adsorbed on the surface of single layer graphene (F/GE) or between the layers of AB stacked bilayer graphene (F@ bilayer graphene). It is found that the diffusion of the F atom is very different in those cases, and that the mobility of the F atom increases by about an order of magnitude when inserted between two graphene layers. The obtained diffusion constant for F/GE is twice larger than that experimentally found for gold adatom and theoretically found for C-60 molecule on graphene. Our study provides important physical insights into the dynamics of fluorine atoms between and on graphene layers and explains the mechanism behind the separation of graphite layers due to intercalation of F atoms.

Abstract: Using highly efficient GPU-based simulations of the tight-binding Bogoliubov-de Gennes equations we solve self-consistently for the pair correlation in rhombohedral (ABC) and Bernal (ABA) multilayer graphene by considering a finite intrinsic s-wave pairing potential. We find that the two different stacking configurations have opposite bulk/surface behavior for the order parameter. Surface superconductivity is robust for ABC stacked multilayer graphene even at very low pairing potentials for which the bulk order parameter vanishes, in agreement with a recent analytical approach. In contrast, for Bernal stacked multilayer graphene, we find that the order parameter is always suppressed at the surface and that there exists a critical value for the pairing potential below which no superconducting order is achieved. We considered different doping scenarios and find that homogeneous doping strongly suppresses surface superconductivity while nonhomogeneous field-induced doping has a much weaker effect on the superconducting order parameter. For multilayer structures with hybrid stacking (ABC and ABA) we find that when the thickness of each region is small (few layers), high-temperature surface superconductivity survives throughout the bulk due to the proximity effect between ABC/ABA interfaces where the order parameter is enhanced. DOI: 10.1103/PhysRevB.87.134509

Abstract: Using improved Wigner-Brillouin perturbation theory we study resonant electron-phonon interaction in a semiconductor quantum dot. We predict pinning of the excited energy levels to the ground state level plus one optical phonon as a function of the strength of the confinement potential. This effect should be observable through optical spectroscopic measurements.

Abstract: Using relativistic density-functional calculations, we examine the magneto-crystalline anisotropy and exchange properties of transition-metal atoms adsorbed on 2D-GaAs. We show that single Mn and Mo atom (Co and Os) strongly bind on 2D-GaAs, and induce local out-of-plane (in-plane) magnetic anisotropy. When a pair of TM atoms is adsorbed on 2D-GaAs in a close range from each other, magnetisation properties change (become tunable) with respect to concentrations and ordering of the adatoms. In all cases, we reveal presence of strong Dzyaloshinskii–Moriya interaction. These results indicate novel pathways towards two-dimensional chiral magnetic materials by design, tailored for desired applications in magneto-electronics.

Abstract: Using the Ginzburg-Landau theory extended to the next-to-leading order, we determine numerically the healing lengths of the two order parameters at the two-gap superconductor/normal metal interface. We demonstrate on several examples that those can be different even in the strict domain of applicability of the Ginzburg-Landau theory. This justifies the use of this theory to describe relevant physics of two-gap superconductors, distinguishing them from their single-gap counterparts. The calculational degree of complexity increases only slightly with respect to the conventional Ginzburg-Landau expansion, thus the extended Ginzburg-Landau model remains numerically far less demanding compared to the full microscopic approaches.

Abstract: Using the tight-binding approach, we calculate the electronic structure of triangular and hexagonal MoS2 quantum dots. Due to the orbital asymmetry we show that it is possible to form quantum dots with the same shape but having different electronic properties. The electronic states of triangular and hexagonal quantum dots are explored, as well as the local and total density of states and the convergence towards the bulk spectrum with dot size is investigated. Our calculations show that: (1) edge states appear in the band gap, (2) that there are a larger number of electronic states in the conduction band as compared to the valence band, and (3) the relative number of edge states decreases with increasing dot size.