Abstract: In a recent paper [Phys. Rev. B 89, 125418 (2014)], the authors argue that it is possible to map the electronic properties of twisted bilayer graphene to those of bilayer graphene in an in-plane magnetic field. However, their description of the low-energy dynamics of twisted bilayer graphene is restricted to the extended zone scheme and therefore neglects the effects of the superperiodic structure. If the energy spectrum is studied in the supercell Brillouin zone, we find that the comparison with an in-plane magnetic field fails because (i) the energy spectra of the two situations exhibit different symmetries and (ii) the low-energy spectra are very different.

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: We investigate the effects related to vortex-core deformations when vortices approach each other. As a result of these vortex-core deformations, the vortex-vortex interaction effectively acquires an attractive component leading to a variety of vortex patterns typical for systems with nonmonotonic repulsive-attractive interaction, such as stripes and labyrinths. The core deformations are anisotropic and can induce frustration in the vortex-vortex interaction. In turn, this frustration has an impact on the resulting vortex patterns, which are analyzed in the presence of additional random pinning, as a function of the pinning strength. This analysis can be applicable to vortices in multiband superconductors or to vortices in Bose-Einstein condensates.

Abstract: We propose that low-strain Si/O superlattices can be constructed by connecting reconstructed Si{001} surfaces by Si-O-Si bridges. Ab initio calculations show that our models are energetically more favorable than all the models proposed so far. The part of our Si/O superlattice model is experimentally accessible just by oxidizing a Si( 001) substrate. To complete our Si/O superlattice model, we propose a three-step method. We also explore the potential of our Si/O superlattice models for new materials used in future Si electronics. We find that the location of the channel where the carriers travel can be controlled between the interfaces and the Si layers by the insertion of O atoms into the Si-Si dimers. By revealing the origins of the interface electron and hole states, we find that similar interface states should be easily achieved for Si slabs and Si substrates. Interestingly, the interface electrons and holes have small effective masses in the direction parallel to the channel and large effective masses in the direction normal to the channel, which makes the Si/O superlattices attractive to be used for channel materials. We also find that the valley splitting of Si is enhanced by the formation of the Si/O/Si interfaces, which is ideal for developing Si-based qubits. Our findings open new perspectives to design and control the electronic properties of Si.

Abstract: Monte Carlo (MC) methods have a long-standing history as partners of molecular dynamics (MD) to simulate the evolution of materials at the atomic scale. Among these techniques, the uniform-acceptance force-bias Monte Carlo (UFMC) method [ G. Dereli Mol. Simul. 8 351 (1992)] has recently attracted attention [ M. Timonova et al. Phys. Rev. B 81 144107 (2010)] thanks to its apparent capacity of being able to simulate physical processes in a reduced number of iterations compared to classical MD methods. The origin of this efficiency remains, however, unclear. In this work we derive a UFMC method starting from basic thermodynamic principles, which leads to an intuitive and unambiguous formalism. The approach includes a statistically relevant time step per Monte Carlo iteration, showing a significant speed-up compared to MD simulations. This time-stamped force-bias Monte Carlo (tfMC) formalism is tested on both simple one-dimensional and three-dimensional systems. Both test-cases give excellent results in agreement with analytical solutions and literature reports. The inclusion of a time scale, the simplicity of the method, and the enhancement of the time step compared to classical MD methods make this method very appealing for studying the dynamics of many-particle systems.

Abstract: Recently, core-shell Ag-Co, Ag-Cu, and “onionlike” Cu-Co equilibrium configurations were predicted in the case of isolated face centered cubic (fcc) bimetallic clusters, and three shell onionlike configurations were predicted in the case of ternary metallic clusters with spherical and truncated octahedral morphologies. In the present paper, immiscible binary CuCo and ternary AgCuCo clusters with icosahedral structures are studied as functions of their size and composition. Clusters studied are formed by 13, 55, 147, 309, and 561 atoms corresponding to the five smallest possible closed shell icosahedral structures. An embedded atom model potential is used to describe their cohesion. Equilibrium configurations are investigated by means of Metropolis Monte Carlo free energy minimization in the (NPT) canonical ensemble. Most simulations are achieved at 10 and 300 K. The effect of temperature on segregation ordering is systematically investigated. Selected cases are used to identify the effect of size and composition on melting. In contrast with fcc clusters, homogeneous onionlike configurations of binary clusters are not predicted. When it is allowed by the composition, a complete outer shell is formed by Cu in binary Cu-Co clusters and by Ag in ternary Ag-Cu-Co clusters. Depending on temperature, Co may precipitate into decahedral groups under the Cu vertices of the icosahedra in binary clusters, while the Co-Cu configuration in ternary clusters drastically depends on the Ag coating. Despite the multicomponent character of the clusters and the immiscibility of the species forming them, for most compositions and sizes, equilibrium structures remain close to perfectly icosahedral at 10 K as well as at 300 K.

Abstract: Multifrequency electron paramagnetic resonance (EPR) and high resolution transmission electron microscopy (HRTEM) investigations were performed on small (2 nm) cubic ZnS nanocrystals (quantum dotsQDs) doped with 0.2% mol Mn2+, self-assembled into a mesoporous structure. The EPR data analysis shows that the substitutional Mn2+ ions are localized at Zn2+ sites subjected to a local axial lattice distortion, resulting in the observed zero-field-splitting parameter |D|=41×10−4 cm−1. The local distortion is attributed to the presence in the second shell of ligands of a stacking fault or twin, which alters the normal stacking sequence of the cubic structure. The HRTEM results confirm the presence of such extended planar defects in a large percentage of the investigated QDs, which makes possible the proposed substitutional Mn2+ impurity ions localization model. Based on these results it is suggested that the high doping levels of Mn2+ ions observed in cubic ZnS and possible in other II-VI semiconductor QDs prepared at low temperatures can be explained by the assistance of the extended lattice defects in the impurities incorporation.

Abstract: Dielectric properties of materials are crucial in describing the electromagnetic response of materials. As devices are becoming considerably smaller than the optical wavelength, the conventional measuring methods based on optical response are limited by their spatial resolution. Electron energy loss spectroscopy performed in a scanning transmission electron microscope is a good alternative to obtain the dielectric properties with excellent spatial resolution. Due to the overlap of diffraction discs in scanning transmission electron microscopy, it is difficult to apply conventional experimental settings to suppress retardation losses. In this contribution, a relativistic dielectric model for the loss function is presented which is used in a model based optimization scheme to estimate the complex dielectric function of a material. The method is applied to experiments on bulk diamond and SrTiO3 and shows a good agreement with optical reference data when retardation effects are included. Application of this technique to nanoparticles is possible but several theoretical assumptions made in the model of the loss function are violated and interpretation becomes problematic.

Abstract: Bi2Se3 is a three-dimensional topological insulator which has been extensively studied because it has a single Dirac cone on the surface, inside a relatively large bulk band gap. However, the effect of two-dimensional topological insulator Bi bilayers on the properties of Bi2Se3 and vice versa, has not been explored much. Bi bilayers are often present between the quintuple layers of Bi2Se3, since (Bi2)n(Bi2Se3)m form stable ground-state structures. Moreover, Bi2Se3 is a good substrate for growing ultrathin Bi bilayers. By first-principles techniques, we first show that there is no preferable surface termination by either Bi or Se. Next, we investigate the electronic structure of Bi bilayers on top of, or inside a Bi2Se3 slab. If the Bi bilayers are on top, we observe a charge transfer to the quintuple layers that increases the binding energy of the surface Dirac cones. The extra states, originating from the Bi bilayers, were declared to form a topological Dirac cone, but here we show that these are ordinary Rashba-split states. This result, together with the appearance of a new Dirac cone that is localized slightly deeper, might necessitate the reinterpretation of several experimental results. When the Bi bilayers are located inside the Bi2Se3 slab, they tend to split the slab into two topological insulators with clear surface states. Interface states can also be observed, but an energy gap persists because of strong coupling between the neighboring quintuple layers and the Bi bilayers.

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 use density functional theory within the local density approximation to investigate the structural, electronic, and optical properties of Al-doped MgO. The concentrations considered range from 6% to 56%. In the latter case, we also compare the optical properties of the amorphous and crystalline phases. We find that, overall, the electronic properties of the crystalline phases change qualitatively little with Al concentration. On the other hand, the changes in the electronic structure in the amorphous phase are more important, most notably because of deep impurity levels in the band gap that are absent in the crystalline phase. This leads to observable effects in, e.g., the optical absorption edge and in the refractive index. Thus, the latter can be used to characterize the crystalline to amorphous transition with Al doping level.