Abstract: It is known that a perpendicular electric field applied to multilayers of graphene modifies the electronic structure near the K point and may induce an energy gap in the electronic spectrum which is tunable by the gate voltage. Here we consider a system of graphene multilayers in the presence of a positively charged top and a negatively charged back gate to control independently the density of electrons on the graphene layers and the Fermi energy of the system. The band structure of three- and four-layer graphene systems in the presence of the top and back gates is obtained using a tight-binding approach. A self-consistent Hartree approximation is used to calculate the induced charges on the different graphene layers. We predict that for opposite and equal charges on the top and bottom layers an energy gap is opened at the Fermi level. For an even number of layers this gap is larger than in the case of an odd number of graphene layers. We find that the circular asymmetry of the spectrum, which is a consequence of the trigonal warping, changes the size of the induced electronic gap, even when the total density of the induced electrons on the graphene layers is low.

Abstract: We derive an exact solution to the problem of spin snake states induced in a nonhomogeneous magnetic field by a combined action of the Rashba spin-orbit and Zeeman fields. The electron spin behavior as a function of the cyclotron orbit center position and an external homogeneous magnetic field was obtained. It is shown that in an antisymmetric magnetic field the electron spin in the snake states has only an in-plane projection, perpendicular to the magnetic interface, which vanishes at large positive momenta. Applying an external homogeneous magnetic field adds a finite out-of-plane spin component and simultaneously gaps out the spectral branches, which results in regular beating patterns of the spin current components.

Abstract: The electronic properties of boron nitride nanoribbons (BNNRs) doped with a line of carbon atoms are investigated using density functional calculations. By replacing a line of alternating B and N atoms with carbons, three different configurations are possible depending on the type of the atoms which bond to the carbons. We found very different electronic properties for these configurations: (i) the NCB arrangement is strongly polarized with a large dipole moment having an unexpected direction, (ii) the BCB and NCN arrangements are nonpolar with zero dipole moment, (iii) the doping by a carbon line reduces the band gap regardless of the local arrangement of the borons and the nitrogens around the carbon line, and (iv) the polarization and energy gap of the carbon-doped BNNRs can be tuned by an electric field applied parallel to the carbon line. Similar effects were found when either an armchair or zigzag line of carbon was introduced.

Abstract: The dynamics of fullerene-cubane (C60⋅C8H8) cocrystal is studied combining experimental [x-ray diffuse scattering, quasielastic and inelastic neutron scattering (INS)] and simulation (molecular dynamics) investigations. Neutron scattering gives direct evidence of the free rotation of fullerenes and of the libration of cubanes in the high-temperature phase, validating the rotor-stator description of this molecular system. X-ray diffuse scattering shows that orientational disorder survives the order/disorder transition in the low-temperature phase, although the loss of fullerene isotropic rotational diffusion is featured by the appearance of a 2.2 meV mode in the INS spectra. The coupling between INS and simulations allows identifying a degeneracy lift of the cubane librations in the low temperature phase, which is used as a tool for probing the environment of cubane in this phase and for getting further insights into the phase transition mechanism.

Abstract: We performed a detailed investigation of the oxidation of the technologically relevant GaAs(001)-beta 2(2x4) surface via density functional calculations. The purpose is to gain insights on the atomistic mechanisms and local bondings that underlie the degradation of the surface properties once exposed to oxygen. The study comprises the adsorption of single O atoms, through the sampling of several adsorption sites, and the subsequent formation of the O adsorbate at increasing coverage by taking into account multiple-atom adsorption. Based on the evaluation of the energetics and the structural properties of the atomistic models generated, the results here reported delineate a consistent picture of the initial stage of the surface oxidation: (i) at low coverage, in the limit of single O insertions, oxygen is incorporated on the surface forming a twofold-bridging Ga-O-As bond; (ii) at increasing coverage, as multiple O atoms are involved, this is accompanied by the formation of a threefold-coordinated bond (with two Ga and one As atoms); (iii) the latter has important implications regarding the electronic properties of the adsorbate since this O bonding may result in the formation of As dangling bonds. Moreover, a clear trend of increased energy gain for the incorporation of neighboring O atoms compared to single O insertions indicates that the formation of oxide clusters is favored over a regime of uniform oxidation. Our findings provide a detailed description of the O bonding and stress the importance of modeling the adsorption of multiple O atoms for an accurate description of the surface oxidation.

Abstract: A data-acquisition and data-processing method is proposed that aims at minimizing the effect of retardation on the Kramers-Kronig analysis of valence electron energy-loss spectra. This method is applied to diamond, which, due to its high dielectric constant, is a material that shows strong retardation effects and thus is a challenging material to be studied by valence electron energy-loss spectroscopy. The results obtained show a significant improvement but still show small discrepancies with respect to optical data, which are most likely due to the residual retardation contributions and the fact that nonzero momentum transfers are measured.

Abstract: Using atomistic simulations we investigate the morphological properties of graphene deposited on top of a nanostructured substrate. Sinusoidally corrugated surfaces, steps, elongated trenches, one-dimensional and cubic barriers, spherical bubbles, Gaussian bumps, and Gaussian depressions are considered as support structures for graphene. The graphene-substrate interaction is governed by van der Waals forces and the profile of the graphene layer is determined by minimizing the energy using molecular dynamics simulations. Based on the obtained optimum configurations, we found that (i) for graphene placed over sinusoidally corrugated substrates with corrugation wavelengths longer than 2 nm, the graphene sheet follows the substrate pattern while for supported graphene it is always suspended across the peaks of the substrate, (ii) the conformation of graphene to the substrate topography is enhanced when increasing the energy parameter in the van der Waals model, (iii) the adhesion of graphene into the trenches depends on the width of the trench and on the graphene's orientation, i. e., in contrast to a small-width (3 nm) nanoribbon with armchair edges, the one with zigzag edges follows the substrate profile, (iv) atomic-scale graphene follows a Gaussian bump substrate but not the substrate with a Gaussian depression, and (v) the adhesion energy due to van der Waals interaction varies in the range [0.1-0.4] J/m(2).

Abstract: Correlated, or extended, impurities play an important role in the transport properties of dirty metals. Here, we examine, in the framework of a tight-binding lattice, the transmission of a single electron through an array of correlated impurities. In particular we show that particles transmit through an impurity array in identical fashion, regardless of the direction of traversal. The demonstration of this fact is straightforward in the continuum limit, but requires a detailed proof for the discrete lattice. We also briefly demonstrate and discuss the time evolution of these scattering states, to delineate regions (in time and space) where the aforementioned symmetry is violated.

Abstract: The influence of a perpendicular magnetic field on the binding energy and structural properties of excitons, trions, and biexcitons in monolayers of semiconducting transition metal dichalcogenides (TMDs) is investigated. The stochastic variational method (SVM) with a correlated Gaussian basis is used to calculate the different properties of these few-particle systems. In addition, we present a simplified variational approach which supports the SVM results for excitons as a function of magnetic field. The exciton diamagnetic shift is compared with recent experimental results, and we extend this concept to trions and biexcitons. The effect of a local potential fluctuation, which we model by a circular potential well, on the binding energy of trions and biexcitons is investigated and found to significantly increase the binding of those excitonic complexes.

Abstract: In a recent paper [Phys. Rev. B 99, 155124 (2019)], the spectrum of a regularized Coulomb charge was studied in gapped pseudospin-1 systems generated by an alpha – T-3 lattice. The electronic spectrum was studied as a function of the impurity strength Z alpha. However, the results and conclusions on the behavior of the flatband states as a function of the impurity strength are incomplete. In this Comment, I argue that because of the dispersionless nature of the flatband, the states spread out under the influence of a charged impurity forming a continuous band of states. I support my arguments with explicit numerical calculations which show the emergence of a continuum of states.

Abstract: Recently, an extremely low lattice thermal conductivity value has been reported for the alkali-based telluride material BaIn2Te4. The value is comparable with low-thermal conductivity metal chalcogenides, and the glass limit is highly intriguing. Therefore, to shed light on this issue, we performed first-principles phonon thermal transport calculations. We predicted highly anisotropic lattice thermal conductivity along different directions via the solution of the linearized phonon Boltzmann transport equation. More importantly, we determined several different factors as the main sources of the predicted ultralow lattice thermal conductivity of this crystal, such as the strong interactions between low-frequency optical phonons and acoustic phonons, small phonon group velocities, and lattice anharmonicity indicated by large negative mode Gruneisen parameters. Along with thermal transport calculations, we also investigated the electronic transport properties by accurately calculating the scattering mechanisms, namely the acoustic deformation potential, ionized impurity, and polar optical scatterings. The inclusion of spin-orbit coupling (SOC) for electronic structure is found to strongly affect the p-type Seebeck coefficients. Finally, we calculated the thermoelectric properties accurately, and the optimal ZT value of p-type doping, which originated from high Seebeck coefficients, was predicted to exceed unity after 700 K and have a direction averaged value of 1.63 (1.76 in the y-direction) at 1000 K around 2 x 1020 cm-3 hole concentration. For n-type doping, a ZT around 3.2 x 1019 cm-3 concentration was predicted to be a direction-averaged value of 1.40 (1.76 in the z-direction) at 1000 K, mostly originating from its high electron mobility. With the experimental evidence of high thermal stability, we showed that the BaIn2Te4 compound has the potential to be a promising mid- to high-temperature thermoelectric material for both p-type and n-type systems with appropriate doping.

Abstract: We study the structural phase transitions in confined systems of strongly interacting particles. We consider infinite quasi-one-dimensional systems with different pairwise repulsive interactions in the presence of an external confinement following a power law. Within the framework of Landaus theory, we find the necessary conditions to observe continuous transitions and demonstrate that the only allowed continuous transition is between the single- and the double-chain configurations and that it only takes place when the confinement is parabolic. We determine analytically the behavior of the system at the transition point and calculate the critical exponents. Furthermore, we perform Monte Carlo simulations and find a perfect agreement between theory and numerics.

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.