Abstract: Using the Ginzburg-Landau (GL) theory, we calculate the stability of sample symmetry-induced vortex-antivortex molecules in a mesoscopic superconducting bilayer exposed to a homogeneous magnetic field. We demonstrate the conditions under which the two condensates cooperatively broaden the field-temperature stability range of the composite (joint) vortex-antivortex state. In cases when such broadening is not achieved, a reentrance of the vortex-antivortex state is found at lower temperatures. In a large portion of the phase diagram noncomposite states are possible, in which the antivortex is present in only one of the layers. In this case, we demonstrate that the vortex-antivortex molecule in one of the layers can be pinned and enlarged by interaction with a vortex molecule in the other. Using analogies in the respective GL formalisms, we map our findings for the bilayer onto mesoscopic two-band superconductors.

Abstract: ZnAl2O4 (gahnite) is a ceramic which is considered a possible transparent conducting oxide (TCO) due to its wide band gap and transparency for UV. Defects play an important role in controlling the conductivity of a TCO material along with the dopant, which is the main source of conductivity in an otherwise insulating oxide. A comprehensive first-principles density functional theory study for point defects in ZnAl2O4 spinel is presented using the Heyd, Scuseria, and Ernzerhof hybrid functional (HSE06) to overcome the band gap problem. We have investigated the formation energies of intrinsic defects which include the Zn, Al, and O vacancy and the antisite defects: Zn at the Al site (ZnAl) and Al at the Zn site (AlZn). The antisite defect AlZn has the lowest formation energy and acts as a shallow donor, indicating possible n-type conductivity in ZnAl2O4 spinel by Al doping.

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: The cyclotron resonance energies, the corresponding oscillator strengths, and the cyclotron absorption spectrum for trilayer graphene are calculated for both ABA and ABC stacking. A gate potential across the stacked layers leads to (1) a reduction of the transition energies, (2) a lifting of the degeneracy of the zero Landau level, and (3) the removal of the electron-hole symmetry.

Abstract: Richardson approach provides an exact solution of the pairing Hamiltonian. This Hamiltonian is characterized by the electron-hole pairing symmetry, which is however hidden in Richardson equations. By analyzing this symmetry and using an additional conjecture, fulfilled in solvable limits, we suggest a simple expression of the ground state energy for an equally-spaced energy-level model, which is applicable along the whole crossover from the superconducting state to the pairing fluctuation regime. Solving Richardson equations numerically, we demonstrate a good accuracy of our expression.

Abstract: The drag of massless fermions in graphene double-layer structures is investigated over a wide range of temperatures and interlayer separations. We show that the inhomogeneity of the dielectric background in such graphene structures, for experimentally relevant parameters, results in a significant enhancement of the drag resistivity. At intermediate temperatures the dynamical screening via plasmon-mediated drag enhances the drag resistivity and results in an upturn in its behavior at large interlayer separations. In a range of interlayer separations, corresponding to the crossover from strong to weak coupling of graphene layers, we find that the decrease of the drag resistivity with interlayer spacing is approximately quadratic. This dependence weakens below this range of interlayer spacing while for larger separations we find a cubic (quartic) dependence at intermediate (low) temperatures.