Table of contents

Volume 32

Number 16, 17 April 2020

Previous issue Next issue

Buy this issue in print

Letter

16LT01

and

The quantum correction to electrical conductivity is studied on the basis of two-dimensional Wolff Hamiltonian, which is an effective model for a spin–orbit coupled (SOC) lattice system. It is shown that weak anti-localization (WAL) arises in SOC lattices, although its mechanism and properties are different from the conventional WAL in normal metals with SOC impurities. The interband SOC effect induces the contribution from the interband singlet Cooperon, which plays a crucial role for WAL in the SOC lattice. It is also shown that there is a crossover from WAL to weak localization in SOC lattices when the Fermi energy or band gap changes. The implications of the present results to Bi–Sb alloys and PbTe under pressure are discussed.

Special Issue Paper

164001
The following article is Open access

, and

Emerging Leaders 2019 Emerging Leaders 2019

We study, by means of an exact analytical solution, the motion of a spheroidal, axisymmetric squirmer in an unbounded fluid, as well as the low Reynolds number hydrodynamic flow associated to it. In contrast to the case of a spherical squirmer—for which, e.g. the velocity of the squirmer and the magnitude of the stresslet associated with the flow induced by the squirmer are respectively determined by the amplitudes of the first two slip ('squirming') modes—for the spheroidal squirmer each squirming mode either contributes to the velocity, or contributes to the stresslet. The results are straightforwardly extended to the self-phoresis of axisymmetric, spheroidal, chemically active particles in the case when the phoretic slip approximation holds.

Papers

Surfaces and interfaces

165001
The following article is Open access

, , , and

Thin films of rhombohedral Sb2Se3 with thicknesses from 1 to 5 quintuple layers (QL) were grown on Bi2Se3/Si(1 1 1) substrate. The electronic band structure of the grown films and the Sb2Se3/Bi2Se3 interface were studied using angle-resolved photoemission spectroscopy. It was found that while Sb2Se3 has an electronic band structure generally similar to that of Bi2Se3, there is no fingerprints of band inversion in it. Instead, the one-QL-thick Sb2Se3 films show direct band gap of about 80 meV. With growing film thickness, the Fermi level of the Sb2Se3 films gradually shifts by 200 meV for 5 QL-thick film revealing the band bending of the Sb2Se3/Bi2Se3 hetero-junction.

Nanostructures and nanoelectronics

165301
The following article is Open access

, and

We theoretically study a silicon triple quantum dot (TQD) system coupled to a superconducting microwave resonator. The response signal of an injected probe signal can be used to extract information about the level structure by measuring the transmission and phase shift of the output field. This information can further be used to gain knowledge about the valley splittings and valley phases in the individual dots. Since relevant valley states are typically split by several , a finite temperature or an applied external bias voltage is required to populate energetically excited states. The theoretical methods in this paper include a capacitor model to fit experimental charging energies, an extended Hubbard model to describe the tunneling dynamics, a rate equation model to find the occupation probabilities, and an input–output model to determine the response signal of the resonator.

165302

, and

The control of geometric structure is a key aspect in the interplay between theoretical predictions and experimental realization in the science and applications of nanomaterials. This is particularly important in one-dimensional structures such as nanoribbons, in which the edge morphology dictates most of the electronic behavior in low energy scale. In the present work we demonstrate by means of first principles calculations that the oxidation of few-layer antimonene may lead to an atomic restructuring with formation of ordered multilayer zig-zag nanoribbons. The widths are uniquely determined by the number of layers of the initial structure, allowing the synthesis of ultranarrow ribbons and chains. We also show that the process may be extended to other compounds based on group V elements, such as arsenene. The characterization of the electronic structure of the resulting ribbons shows an important effect of stacking on band gaps and on modulation of electronic behavior.

Structure, dynamics and phase transitions

165401

, , and

By introducing the next-nearest-neighboring (NNN) intersite coupling, we investigate the eigenenergies of the -symmetric non-Hermitian Su–Schrieffer–Heeger (SSH) model with two conjugated imaginary potentials at the end sites. It is found that with the strengthening of NNN coupling, the particle-hole symmetry is destroyed. As a result, the bonding band is first narrowed and then undergoes the top-bottom reversal followed by the its width's increase, whereas the antibonding band is widened monotonously. In this process, the topological state extends into the topologically-trivial region, and its energy departs from the energy zero point, accompanied by the emergence of one new topological state in this region. All these results give rise to the complication of the topological properties and the manner of -symmetry breaking. It can be concluded that the NNN coupling takes important effects to the change of the topological properties of the non-Hermitian SSH system.

Electronic structure

165501

, and

We address the problem of hybridization between topological surface states and a non-topological flat bulk band. Our model, being a mixture of three-dimensional Bernevig–Hughes–Zhang and two-dimensional pseudospin-1 Hamiltonian, allows explicit treatment of the topological surface state evolution by continuously changing the hybridization between the inverted bands and an additional 'parasitic' flat band in the bulk. We show that the hybridization with a flat band lying below the edge of the conduction band converts the initial Dirac-like surface states into a branch below and one above the flat band. Our results univocally demonstrate that the upper branch of the topological surface states is formed by Dyakonov–Khaetskii surface states, known for HgTe since the 1980s. Additionally we explore an evolution of the surface states and the arising of Fermi arcs in Dirac semimetals when the flat band crosses the conduction band.

165502

, and

Vacancies in diamond can introduce defect bands in the band gap of diamond and lead to the sub-band gap absorption in the visible and infrared regions. At sufficiently high concentrations of vacancy, the defect bands overlap each other to form a single partially filled intermediate band (IB) in the band gap and the sub-band gap absorption in the infrared region is especially strong. Along with the decreasing of the vacancy concentrations, the IB splits into two separate bands and the sub-band gap absorption decreases sharply. The computed results are important for thoroughly understanding of the optical properties of black diamond.

Correlated electron systems

165601

, , , and

The non-equilibrium electronic transport through a nanoscale device composed of a single quantum dot between two metallic contacts is studied theoretically within the framework of the Keldysh formalism. The quantum dot consists of a single energy level subject to an applied magnetic field. Correlations due to the Coulomb repulsion between electrons on the dot are treated with a Green's function decoupling scheme which, although similar to the Hubbard-I approximation, captures some of the dynamics beyond. The scheme is exact in the so-called atomic limit, defined by vanishing tunneling between contacts and dot, and in the non-interacting limit, where the on-dot Coulomb repulsion is zero. Explicit analytic solutions, valid for arbitrary magnetic fields, are obtained for two important setups: (i) the stationary regime, with constant voltage bias between the leads, and (ii) the time-dependent regime for metallic leads with constant density of states of infinite width. In these regimes, the current through the dot is evaluated numerically for various parameter sets and its main features interpreted in terms of the underlying physical processes. The results are compared to the non-crossing approximation (NCA) and diagrammatic non-equilibrium quantum Monte-Carlo (QMC) where available.

Physics of materials

165701

, , , , and

Owing to the various ways of chemical bonding, carbon can form abundant allotropes with different frameworks, which harbor rich mechanical and electronic properties. Taking the cage-like isomer of C16 cluster as a building block, we design a new low-density carbon allotrope, which has tetragonal symmetry (I4/mmm) and a 56-atom unit cell, hence termed as T-C56. Our first-principles calculations reveal that T-C56 is not only energetically, dynamically, thermally (above 1800 K) and mechanically stable, but even more stable than the experimentally synthesized C20-sc and T-carbon. Remarkably, although the framework of T-C56 is low density (2.72 g cm−3), it exhibits novel superhard properties with a Vickers hardness of 48.71 GPa. The obtained Yang's modulus and ideal strength show that T-C56 is mechanically anisotropic. In particular, our analysis of electronic and optical properties suggest that T-C56 is a transparent indirect semiconductor with a wide bandgap of 3.18 eV (HSE06). These findings highlight a distinct carbon allotrope, having promising applications for optical and aerospace devices due to its light, transparent, superhard features.

165702
The following article is Open access

, , , , , , , , and

Tungsten tetraboride (WB4)-based solid solutions represent one of the most promising superhard metal candidates; however, their underlying hardening mechanisms have not yet been fully understood. Here, we explore the lattice compressibility of WB4 binary solid solutions with different manganese (Mn) concentrations using high-pressure x-ray diffraction (XRD) up to 52 GPa. Under initial compression, the lattices of low and high Mn-doped WB4 alloys (i.e. W0.96Mn0.04B4 and W0.84Mn0.16B4) are shown to be more and less compressible than pure WB4, respectively. Then, a c-axis softening is found to occur above 39 GPa in WB4, consistent with previous results. However, an anomalous sudden a-axis stiffening is revealed at ~36 GPa in W0.96Mn0.04B4, along with suppression of c-axis softening observed in WB4. Furthermore, upon Mn addition, a simultaneous stiffening of a- and c-axes is demonstrated in W0.84Mn0.16B4 at ~37 GPa. Speculation on the possible relationship between this anomalous stiffening and the combined effects of valence-electron concentration (VEC) and atomic size mismatch is also included to understand the origin of the nearly identical hardness enhancement in those two solid solutions compared to WB4. Our findings emphasize the importance of accurate bonding and structure manipulation via solute atoms to best optimize the hardness of WB4 solid solutions.

Magnetism

165801
The following article is Open access

, , , , , , , , , et al

The ability to create atomically perfect, epitaxial heterostructures of correlated complex perovskite oxides using state-of-art thin film deposition techniques has generated new physical phenomena at engineered interfaces. Here we report on the impact of growth kinetics on the magnetic structure and exchange coupling at the interface in heterostructures combining layers of antiferromagnetic La1/3Sr2/3FeO3 (LSFO) and ferromagnetic La2/3Sr1/3MnO3 (LSMO) on (0 0 1)-oriented SrTiO3 (STO) substrates. Two growth orders are investigated, (a) LSMO/LSFO/STO(0 0 1) and (b) LSFO/LSMO/STO(0 0 1), where the LSFO layer is grown by molecular beam epitaxy and the LSMO layer by high oxygen pressure sputtering. The interface has been investigated using electron microscopy and polarized neutron reflectometry. Interdiffusion over seven monolayers is observed in LSMO/LSFO (a) with an almost 50% reduction in magnetization at the interface and showing no exchange coupling. However, the exchange bias effect ( mT at 10 K) could be realized when the interface is atomically sharp, as in LSFO/LSMO (b). Our study therefore reveals that, even for well ordered and lattice-matched structures, the kinetics involved in the growth processes drastically influences the interface quality with a strong correlation to the magnetic properties.

165802

, , and

We have studied a Fe-based di-nuclear molecular complex having the chemical formula [{Fe(bpp)(NCS)2}2('-bipy)]·2MeOH (where bpp  =  -bis(pyrazol-3-yl) pyridine and '-bipy  =  '-bipyridine, 1) using density functional theory and model Hamiltonian approach. Our study provides insight to the pressure driven spin-crossover (SCO) phenomena observed experimentally in these systems. Upon increasing the pressure, the spin state of Fe(II) cation gradually changes from a high spin state (S  =2) to a low spin (LS) state (S  =0) accompanied by volume contraction. The gradual increase in pressure shrinks Fe–N bond length and also causes angular deviation of the FeN6 octahedron leading to full conversion to the LS state without global structural phase transition. We have carried out exact diagonalization study of an effective single site Hamiltonian and confirmed the importance of intramolecular interaction for SCO phenomena. We have investigated the cooperativity of the observed SCO phenomena. We have also studied the effect of Co doping on the spin state of Fe and find that the spin state of Fe has a subtle dependency on the concentration of dopant atoms. Excess Co doping pave the way towards the possibility of an intermediate spin state for Fe and can give rise to a bistable spin transition process.

165803

, , , and

We report the coexistence of the Kondo effect and spin glass behavior in Fe-doped NbS2 single crystals. The FexNbS2 shows the resistance minimum and negative magnetoresistance due to the Kondo effect, and exhibits no superconducting behavior at low temperatures. The resistance curve follows a numerical renormalization-group theory using the Kondo temperature K for x  =  0.01 as evidence of Kondo effect. Scanning tunneling microscope/spectroscopy (STM/STS) revealed the presence of Fe atoms near sulfur atoms and asymmetric spectra. The magnetic susceptibility exhibits a feature of spin glass. The static critical exponents determined by the universal scaling of the nonlinear part of the susceptibility suggest a three-dimensional Heisenberg spin glass. The doped-Fe atoms in the intra- and inter-layers revealed by the x-ray result can realize the coexistence of the Kondo effect and spin glass.

165804

and

Recently, it has been shown that two dimensional frustrated mixed-spin systems with anisotropic exchange interactions display supersolid phases in their ground state phase diagrams even in the absence of long-range interactions. In this paper, using cluster mean field theory, we investigate the effects of thermal fluctuations on the ground state phases of this kind of system and show that various thermal solids and thermal insulators emerge around the ground state solid and Mott insulating phases. We also study the thermodynamic properties and magnetocaloric effect of these systems and demonstrate that at low temperatures, a large cooling rate is seen in the vicinity of the solid-supersolid, solid-superfluid and Mott insulator-superfluid critical points, with the large accumulation of the entropy and the minimums of the isentropes. Our results show the sign change of the magnetocaloric parameter inside the solids and the Mott insulator, which is a characteristic of ordered phases.

165805

, and

We have formulated a twist operator argument for the geometrically frustrated quantum spin systems on the kagome and triangular lattices, thereby extending the application of the Lieb–Schultz–Mattis and Oshikawa–Yamanaka–Affleck theorems to these systems. The equivalent large gauge transformation for the geometrically frustrated lattice differs from that for non-frustrated systems due to the existence of multiple sublattices in the unit cell and non-orthogonal basis vectors. Our study for the S  =  1/2 kagome Heisenberg antiferromagnet at zero external magnetic field gives a criterion for the existence of a two-fold degenerate ground state with a finite excitation gap and fractionalized excitations. At finite field, we predict various plateaux at fractional magnetisation, in analogy with integer and fractional quantum Hall states of the primary sequence. These plateaux correspond to gapped quantum liquid ground states with a fixed number of singlets and spinons in the unit cell. A similar analysis for the triangular lattice predicts a single fractional magnetization plateau at 1/3. Our results are in broad agreement with numerical and experimental studies.

165806

, and

First principles calculations reveal that the antiferromagnetic double layer (AFMD) along the magnetization direction of 〈0 0 1〉 is the magnetic ground state of face-centered-cubic (FCC) structure of iron (Fe) due to its lowest total energy among all the studied magnetic states. This magnetic ground state of AFMD-〈0 0 1〉 is fundamentally due to a stronger chemical bonding in terms of electronic structure, and could be confirmed by its mechanical properties for the first time. Calculations also show that lattice constant has an important effect to determine the magnetic state of FCC Fe, and a transition of magnetic state between AFMD-〈0 0 1〉 and ferromagnetic happens at the critical lattice constant of 3.666 Å. The derived results are in good agreement with experimental and theoretical observations, and could clarify the controversy regarding magnetic ground state of FCC Fe in the literature.

Computational and experimental methods

165901

, , , and

In this work we describe the polarization and depolarization process of Eu2+–VC dipoles in Eu2+ doped solid solutions of KCl and KBr considering in detail the different jumping paths of the cation vacancy within individual arrangements of minority anions on the first two anion shells surrounding the dipolar complex. Experimental depolarization curves of binary crystals show in contrast to the single depolarization peak in unitary crystals, up to at least two additional depolarization peaks. These two additional peaks are associated with jumps perturbed by one, and two minority anions flanking the jumping path. In crystals for which Cl is the minority anion, the formation of the two additional bands is more prominent than for those where Br is the minority anion, suggesting that the divalent cation dopant has more affinity towards the Cl minority anion than towards the Br minority anion. Fittings of our model to the experimental data in which we also contemplate different affinities with the divalent dopant cation for Cl and Br minority ions, confirm this. In contrast to an earlier work, in which it had been suggested that only the first anion shell contributes to the perturbation of the dipole in order to reconcile the experimentally observed peak intensities and the supposed affinity between the dipole and the Cl minority anion, the results of our work suggest that both shells contribute in a similar way. This discrepancy is discussed in terms of the different approaches to the description of depolarization process.

165902
The following article is Open access

, , , , , , , , , et al

Wannier90 is an open-source computer program for calculating maximally-localised Wannier functions (MLWFs) from a set of Bloch states. It is interfaced to many widely used electronic-structure codes thanks to its independence from the basis sets representing these Bloch states. In the past few years the development of Wannier90 has transitioned to a community-driven model; this has resulted in a number of new developments that have been recently released in Wannier90 v3.0. In this article we describe these new functionalities, that include the implementation of new features for wannierisation and disentanglement (symmetry-adapted Wannier functions, selectively-localised Wannier functions, selected columns of the density matrix) and the ability to calculate new properties (shift currents and Berry-curvature dipole, and a new interface to many-body perturbation theory); performance improvements, including parallelisation of the core code; enhancements in functionality (support for spinor-valued Wannier functions, more accurate methods to interpolate quantities in the Brillouin zone); improved usability (improved plotting routines, integration with high-throughput automation frameworks), as well as the implementation of modern software engineering practices (unit testing, continuous integration, and automatic source-code documentation). These new features, capabilities, and code development model aim to further sustain and expand the community uptake and range of applicability, that nowadays spans complex and accurate dielectric, electronic, magnetic, optical, topological and transport properties of materials.