Abstract: We report on the crystal structure and magnetism of the iron-based oxyfluoride Bi4Fe5O13F, a material prototype of the Cairo pentagonal spin lattice. The crystal structure of Bi4Fe5O13F is determined by a combination of neutron diffraction, synchrotron x-ray diffraction, and transmission electron microscopy. It comprises layers of FeO6 octahedra and FeO4 tetrahedra forming deformed pentagonal units. The topology of these layers resembles a pentagonal least-perimeter tiling, which is known as the Cairo lattice. This topology gives rise to frustrated exchange couplings and underlies a sequence of magnetic transitions at T-1 = 62 K, T-2 = 71 K, and T-N = 178 K, as determined by thermodynamic measurements and neutron diffraction. Below T-1, Bi4Fe5O13F forms a fully ordered non-collinear antiferromagnetic structure, whereas the magnetic state between T-1 and T-N may be partially disordered according to the sizable increase in the magnetic entropy at T-1 and T-2. Bi4Fe5O13F reveals unanticipated magnetic transitions on the pentagonal Cairo spin lattice and calls for a further work on finite-temperature properties of this strongly frustrated spin model. DOI: 10.1103/PhysRevB.87.024423

Abstract: A procedure to derive the Ginzburg-Landau (GL) theory from the multiband BCS Hamiltonian is developed in a general case with an arbitrary number of bands and arbitrary interaction matrix. It combines the standard Gor'kov truncation and a subsequent reconstruction in order to match accuracies of the obtained terms. This reconstruction recovers the phenomenological GL theory as obtained from the Landau model of phase transitions but offers explicit microscopic expressions for the relevant parameters. Detailed calculations are presented for a three-band system treated as a prototype multiband superconductor. It is demonstrated that the symmetry in the coupling matrix may lead to the chiral ground state with the phase frustration, typical for systems with broken time-reversal symmetry. DOI: 10.1103/PhysRevB.87.134510

Abstract: Density-functional tight-binding and classical molecular dynamics simulations are used to investigate the structural deformations and melting of planar carbon nanoclusters C-N with N = 2-55. The minimum-energy configurations for different clusters are used as starting configurations for the study of the temperature effects on the bond breaking and rotation in carbon lines (N < 6), carbon rings (5 < N < 19), and graphene nanoflakes. The larger the rings (graphene nanoflakes) the higher the transition temperature (melting point) with ring-to-line (perfect-to-defective) transition structures. The melting point was obtained by using the bond energy, the Lindemann criteria, and the specific heat. We found that hydrogen-passivated graphene nanoflakes (CNHM) have a larger melting temperature with a much smaller dependence on size. The edges in the graphene nanoflakes exhibit several different metastable configurations (isomers) during heating before melting occurs. DOI: 10.1103/PhysRevB.87.134103

Abstract: Motivated by recent experimental observations of Tongay et al. [Nano Lett. 12, 5576 (2012)] we show how the electronic properties and Raman characteristics of single layer MoSe2 are affected by elastic biaxial strain. We found that with increasing strain: (1) the E' and E '' Raman peaks (E-2g and E-1g in bulk) exhibit significant redshifts (up to similar to 30 cm(-1)), (2) the position of the A'(1) peak remains at similar to 180 cm(-1) (A(1g) in bulk) and does not change considerably with further strain, (3) the dispersion of low energy flexural phonons crosses over from quadratic to linear, and (4) the electronic band structure undergoes a direct to indirect band gap crossover under similar to 3% biaxial tensile strain. Thus the application of strain appears to be a promising approach for a rapid and reversible tuning of the electronic, vibrational, and optical properties of single layer MoSe2 and similar MX2 dichalcogenides. DOI:10.1103/PhysRevB.87.125415

Abstract: We study Dirac plasmons and their damping in spatially separated N-layer graphene structures at finite doping and temperatures. The plasmon spectrum consists of one optical excitation with square-root dispersion and N – 1 acoustical excitations with linear dispersion, which are undamped at zero temperature and finite doping within a triangular energy region outside the electron-hole continuum. In the long-wavelength limit the energy and weight of the optical plasmon modes increase, respectively, as the square root and linearly with N in agreement with recent experimental findings. The energy and weight of the upper-lying acoustical branches also increase with N. This increase is strongest for the uppermost acoustical mode, and we find that its energy can exceed at some value of momentum the plasmon energy in an individual graphene sheet. Meanwhile, the energy of the low-lying acoustical branches decreases weakly with N as compared with the single acoustical mode in double-layer graphene structures. Our numerical calculations provide a detailed understanding of the overall behavior of the wave-vector dependence of the optical and acoustical multilayer plasmon modes and show how their dispersion and damping are modified as a function of temperature, interlayer spacing, and inlayer carrier density in (un)balanced graphene multilayer structures. DOI: 10.1103/PhysRevB.87.085401

Abstract: We investigate the magnetic interface states of graphene quantum dots that contain p-n junctions. Within a tight-binding approach, we consider rectangular quantum dots in the presence of a perpendicular magnetic field containing p-n as well as p-n-p and n-p-n junctions. The results show the interplay between the edge states associated with the zigzag terminations of the sample and the snake states that arise at the p-n junction due to the overlap between electron and hole states at the potential interface. Remarkable localized states are found at the crossing of the p-n junction with the zigzag edge having a dumb-bell-shaped electron distribution. The results are presented as a function of the junction parameters and the applied magnetic flux. DOI: 10.1103/PhysRevB.87.035426

Abstract: Neutral excitons in strained axially symmetric In(Ga)As/GaAs quantum dots with a ringlike shape are investigated. Similar to experimental self-assembled quantum rings, the analyzed quantum dots have volcano-like shapes. The continuum mechanical model is employed to determine the strain distribution, and the single-band envelope function approach is adopted to compute the electron states. The hole states are determined by the axially symmetric multiband Luttinger-Kohn Hamiltonian, and the exciton states are obtained from an exact diagonalization. We found that the presence of the inner layer covering the ring opening enhances the excitonic Aharonov-Bohm (AB) oscillations. The reason is that the hole becomes mainly localized in the inner part of the quantum dot due to strain, whereas the electron resides mainly inside the ring-shaped rim. Interestingly, larger AB oscillations are found in the analyzed quantum dot than in a fully opened quantum ring of the same width. Comparison with the unstrained ringlike quantum dot shows that the amplitude of the excitonic Aharonov-Bohm oscillations are almost doubled in the presence of strain. The computed oscillations of the exciton energy levels are comparable in magnitude to the oscillations measured in recent experiments. DOI: 10.1103/PhysRevB.87.085314

Abstract: Nonlinear dynamics of Josephson vortices (fluxons) in artificial stacks of superconducting-normal-superconducting Josephson junctions under simultaneously applied time-periodic ac and constant biasing dc currents is studied using the time dependent Ginzburg-Landau formalism with a Lawrence-Doniach extension. At zero external magnetic field and dc biasing current the resistive state of the system is characterized by periodic nucleation and annihilation of fluxon-antifluxon pairs, relative positions of which are determined by the state of neighboring junctions. Due to the mutual repulsive interaction, fluxons in different junctions move out of phase. Their collective motion can be synchronized by adding a small ac component to the biasing dc current. Coherent motion of fluxons is observed for a broad frequency range of the applied drive. In the coherent state the maximal output voltage, which is proportional to the number of junctions in the stack, is observed near the characteristic frequency of the system determined by the crossing of the fluxons across the sample. However, in this frequency range the dynamically synchronized state has an alternative-a less ordered state with smaller amplitude of the output voltage. Collective behavior of the junctions is strongly affected by the sloped sidewalls of the stack. Synchronization is observed only for weakly trapezoidal cross sections, whereas irregular motion of fluxons is observed for larger slopes of the sample edge.

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: In this work we have investigated the structural and electronic properties of SnO, which is built up from layers kept together by van der Waals (vdW) forces. The combination of a vdW functional within density functional theory (DFT) and quasiparticle band structure calculations within the GW approximation provides accurate values for the lattice parameters, atomic positions, and the electronic band structure including the fundamental (indirect) and the optical (direct) band gap without the need of experimental or empirical input. A systematic comparison is made between different levels of self-consistency within the GW approach {following the scheme of Shishkin et al. [Phys. Rev. B 75, 235102 (2007)]} and the results are compared with DFT and hybrid functional results. Furthermore, the effect of the vdW-corrected functional as a starting point for the GW calculation of the band gap has been investigated. Finally, we studied the effect of the vdW functional on the electron charge density.

Abstract: We theoretically demonstrate that the transport inefficiency recently found experimentally for branched-out mesoscopic networks can also be observed in a quantum ring of finite width with an attached central horizontal branch. This is done by investigating the time evolution of an electron wave packet in such a system. Our numerical results show that the conductivity of the ring does not necessary improve if one adds an extra channel. This ensures that there exists a quantum analog of the Braess paradox, originating from quantum scattering and interference.

Abstract: A parabolic quantum dot (QD) as realized by biasing nanostructured gates on bilayer graphene is investigated in the presence of electron-electron interaction. The energy spectrum and the phase diagram reveal unexpected transitions as a function of a magnetic field. For example, in contrast to semiconductor QDs, we find a valley transition rather than only the usual singlet-triplet transition in the ground state of the interacting system. The origin of these features can be traced to the valley degree of freedom in bilayer graphene. These transitions have important consequences for cyclotron resonance experiments.

Abstract: The electronic properties of a triaxially strained hexagonal graphene flake with either armchair or zigzag edges are investigated using molecular dynamics simulations and tight-binding calculations. We found that (i) the pseudomagnetic field in strained graphene flakes is not uniform neither in the center nor at the edge of zigzag terminated flakes, (ii) the pseudomagnetic field is almost zero in the center of armchair terminated flakes but increases dramatically near the edges, (iii) the pseudomagnetic field increases linearly with strain, for strains lower than 15% but increases nonlinearly beyond it, (iv) the local density of states in the center of the zigzag hexagon exhibits pseudo-Landau levels with broken sublattice symmetry in the zeroth pseudo-Landau level, and in addition there is a shift in the Dirac cone due to strain induced scalar potentials, and (v) there is size effect in pseudomagnetic field. This study provides a realistic model of the electronic properties of inhomogeneously strained graphene where the relaxation of the atomic positions is correctly included together with strain induced modifications of the hopping terms up to next-nearest neighbors.

Abstract: Based on the metastable electron-pair energy band in a two-dimensional (2D) periodic potential obtained previously by Hai and Castelano [J. Phys.: Condens. Matter 26, 115502 (2014)], we present in this work a Hamiltonian of many electrons consisting of single electrons and electron pairs in the 2D system. The electron-pair states are metastable of energies higher than those of the single-electron states at low electron density. We assume two different scenarios for the single-electron band. When it is considered as the lowest conduction band of a crystal, we compare the obtained Hamiltonian with the phenomenological model Hamiltonian of a boson-fermion mixture proposed by Friedberg and Lee [Phys. Rev. B 40, 6745 (1989)]. Single-electron-electron-pair and electron-pair-electron-pair interaction terms appear in our Hamiltonian and the interaction potentials can be determined from the electron-electron Coulomb interactions. When we consider the single-electron band as the highest valence band of a crystal, we show that holes in this valence band are important for stabilization of the electron-pair states in the system.