One-Dimensional Magnetic Conduction Channels across Zigzag Graphene Nanoribbon/Hexagonal Boron Nitride Heterojunctions.

We examine the electronic structure of recently fabricated in-plane heterojunctions of zigzag graphene nanoribbons embedded in hexagonal boron nitride. We focus on hitherto unexplored interface configurations in which both edges of the nanoribbon are bonded to the same chemical species, either boron or nitrogen atoms. Using ab initio and mean-field Hubbard model calculations, we reveal the emergence of one-dimensional magnetic conducting channels at these interfaces. These channels originate from the energy shift of the magnetic interface states that is induced by charge transfer between the nanoribbon and hexagonal boron nitride. We further address the response of these heterojunctions to external electric and magnetic fields, demonstrating the tunability of energy and spin splittings in the electronic structure. Our findings establish that zigzag graphene nanoribbon/hexagonal boron nitride heterojunctions are a suitable platform for exploring and engineering spin transport in the atomically thin limit, with potential applications in integrated spintronic devices.

Dirac Half-Semimetallicity and Antiferromagnetism in Graphene Nanoribbon/Hexagonal Boron Nitride Heterojunctions.

Half-metals have been envisioned as active components in spintronic devices by virtue of their completely spin-polarized electrical currents. Actual materials hosting half-metallic phases, however, remain scarce. Here, we predict that recently fabricated heterojunctions of zigzag nanoribbons embedded in two-dimensional hexagonal boron nitride are half-semimetallic, featuring fully spin-polarized Dirac points at the Fermi level. The half-semimetallicity originates from the transfer of charges from hexagonal boron nitride to the embedded graphene nanoribbon. These charges give rise to opposite energy shifts of the states residing at the two edges, while preserving their intrinsic antiferromagnetic exchange coupling. Upon doping, an antiferromagnetic-to-ferrimagnetic phase transition occurs in these heterojunctions, with the sign of the excess charge controlling the spatial localization of the net magnetic moments. Our findings demonstrate that such heterojunctions realize tunable one-dimensional conducting channels of spin-polarized Dirac fermions seamlessly integrated into a two-dimensional insulator, thus holding promise for the development of carbon-based spintronics.

Electrons Surf Phason Waves in Moiré Bilayers.

We investigate the effect of thermal fluctuations on the atomic and electronic structure of a twisted MoSe2/WSe2 heterobilayer using a combination of classical molecular dynamics and ab initio density functional theory calculations. Our calculations reveal that thermally excited phason modes give rise to an almost rigid motion of the moiré lattice. Electrons and holes in low-energy states are localized in specific stacking regions of the moiré unit cell and follow the thermal motion of these regions. In other words, charge carriers surf phason waves that are excited at finite temperatures. We also show that such surfing survives in the presence of a substrate and frozen potential. This effect has potential implications for the design of charge and exciton transport devices based on moiré materials.

Unveiling and Manipulating Hidden Symmetries in Graphene Nanoribbons.

Armchair graphene nanoribbons are a highly promising class of semiconductors for all-carbon nanocircuitry. Here, we present a new perspective on their electronic structure from simple model Hamiltonians and ab initio calculations. We focus on a specific set of nanoribbons of width n=3p+2, where n is the number of carbon atoms across the nanoribbon axis and p is a positive integer. We demonstrate that the energy-gap opening in these nanoribbons originates from the breaking of a previously unidentified hidden symmetry by long-ranged hopping of π electrons and structural distortions occurring at the edges. This hidden symmetry can be restored or manipulated through the application of in-plane lattice strain, which enables continuous energy-gap tuning, the emergence of Dirac points at the Fermi level, and topological quantum phase transitions. Our work establishes an original interpretation of the semiconducting character of armchair graphene nanoribbons and offers guidelines for rationally designing their electronic structure.

Electrically Induced Dirac Fermions in Graphene Nanoribbons.

Graphene nanoribbons are widely regarded as promising building blocks for next-generation carbon-based devices. A critical issue to their prospective applications is whether their electronic structure can be externally controlled. Here, we combine simple model Hamiltonians with extensive first-principles calculations to investigate the response of armchair graphene nanoribbons to transverse electric fields. Such fields can be achieved either upon laterally gating the nanoribbon or incorporating ambipolar chemical codopants along the edges. We reveal that the field induces a semiconductor-to-semimetal transition with the semimetallic phase featuring zero-energy Dirac fermions that propagate along the armchair edges. The transition occurs at critical fields that scale inversely with the width of the nanoribbons. These findings are universal to group-IV honeycomb lattices, including silicene and germanene nanoribbons, irrespective of the type of edge termination. Overall, our results create new opportunities to electrically engineer Dirac semimetallic phases in otherwise semiconducting graphene-like nanoribbons.

OPTIMADE, an API for exchanging materials data.

The Open Databases Integration for Materials Design (OPTIMADE) consortium has designed a universal application programming interface (API) to make materials databases accessible and interoperable. We outline the first stable release of the specification, v1.0, which is already supported by many leading databases and several software packages. We illustrate the advantages of the OPTIMADE API through worked examples on each of the public materials databases that support the full API specification.

Accurate and Efficient Computation of Optical Absorption Spectra of Molecular Crystals: The Case of the Polymorphs of ROY.

When calculating the optical absorption spectra of molecular crystals from first principles, the influence of the crystalline environment on the excitations is of significant importance. For such systems, however, methods to describe the excitations accurately can be computationally prohibitive due to the relatively large system sizes involved. In this work, we demonstrate a method that allows optical absorption spectra to be computed both efficiently and at high accuracy. Our approach is based on the spectral warping method successfully applied to molecules in solvent. It involves calculating the absorption spectrum of a supercell of the full molecular crystal using semi-local time-dependent density functional theory (TDDFT), before warping the spectrum using a transformation derived from smaller-scale semi-local and hybrid TDDFT calculations on isolated dimers. We demonstrate the power of this method on three polymorphs of the well-known color polymorphic compound ROY and find that it outperforms both small-scale hybrid TDDFT dimer calculations and large-scale semi-local TDDFT supercell calculations, when compared to the experiment.

The CECAM electronic structure library and the modular software development paradigm.

First-principles electronic structure calculations are now accessible to a very large community of users across many disciplines, thanks to many successful software packages, some of which are described in this special issue. The traditional coding paradigm for such packages is monolithic, i.e., regardless of how modular its internal structure may be, the code is built independently from others, essentially from the compiler up, possibly with the exception of linear-algebra and message-passing libraries. This model has endured and been quite successful for decades. The successful evolution of the electronic structure methodology itself, however, has resulted in an increasing complexity and an ever longer list of features expected within all software packages, which implies a growing amount of replication between different packages, not only in the initial coding but, more importantly, every time a code needs to be re-engineered to adapt to the evolution of computer hardware architecture. The Electronic Structure Library (ESL) was initiated by CECAM (the European Centre for Atomic and Molecular Calculations) to catalyze a paradigm shift away from the monolithic model and promote modularization, with the ambition to extract common tasks from electronic structure codes and redesign them as open-source libraries available to everybody. Such libraries include "heavy-duty" ones that have the potential for a high degree of parallelization and adaptation to novel hardware within them, thereby separating the sophisticated computer science aspects of performance optimization and re-engineering from the computational science done by, e.g., physicists and chemists when implementing new ideas. We envisage that this modular paradigm will improve overall coding efficiency and enable specialists (whether they be computer scientists or computational scientists) to use their skills more effectively and will lead to a more dynamic evolution of software in the community as well as lower barriers to entry for new developers. The model comes with new challenges, though. The building and compilation of a code based on many interdependent libraries (and their versions) is a much more complex task than that of a code delivered in a single self-contained package. Here, we describe the state of the ESL, the different libraries it now contains, the short- and mid-term plans for further libraries, and the way the new challenges are faced. The ESL is a community initiative into which several pre-existing codes and their developers have contributed with their software and efforts, from which several codes are already benefiting, and which remains open to the community.

Correction: Atomistic QM/MM simulations of the strength of covalent interfaces in carbon nanotube-polymer composites.

Correction for 'Atomistic QM/MM simulations of the strength of covalent interfaces in carbon nanotube-polymer composites' by Jacek R. Golebiowski et al., Phys. Chem. Chem. Phys., 2020, 22, 12007-12014, DOI: 10.1039/d0cp01841d.

The ONETEP linear-scaling density functional theory program.

We present an overview of the onetep program for linear-scaling density functional theory (DFT) calculations with large basis set (plane-wave) accuracy on parallel computers. The DFT energy is computed from the density matrix, which is constructed from spatially localized orbitals we call Non-orthogonal Generalized Wannier Functions (NGWFs), expressed in terms of periodic sinc (psinc) functions. During the calculation, both the density matrix and the NGWFs are optimized with localization constraints. By taking advantage of localization, onetep is able to perform calculations including thousands of atoms with computational effort, which scales linearly with the number or atoms. The code has a large and diverse range of capabilities, explored in this paper, including different boundary conditions, various exchange-correlation functionals (with and without exact exchange), finite electronic temperature methods for metallic systems, methods for strongly correlated systems, molecular dynamics, vibrational calculations, time-dependent DFT, electronic transport, core loss spectroscopy, implicit solvation, quantum mechanical (QM)/molecular mechanical and QM-in-QM embedding, density of states calculations, distributed multipole analysis, and methods for partitioning charges and interactions between fragments. Calculations with onetep provide unique insights into large and complex systems that require an accurate atomic-level description, ranging from biomolecular to chemical, to materials, and to physical problems, as we show with a small selection of illustrative examples. onetep has always aimed to be at the cutting edge of method and software developments, and it serves as a platform for developing new methods of electronic structure simulation. We therefore conclude by describing some of the challenges and directions for its future developments and applications.

Atomistic QM/MM simulations of the strength of covalent interfaces in carbon nanotube-polymer composites.

We investigate the failure of carbon-nanotube/polymer composites by using a recently-developed hybrid quantum-mechanical/molecular-mechanical (QM/MM) approach to simulate nanotube pull-out from a cross-linked polyethene matrix. Our study focuses on the strength and failure modes of covalently-bonded nanotube-polymer interfaces based on amine, carbene and carboxyl functional groups and a [2+1] cycloaddition. We find that the choice of the functional group linking the polymer matrix to the nanotube determines the effective strength of the interface, which can be increased by up to 50% (up to the limit dictated by the strength of the polymer backbone itself) by choosing groups with higher interfacial binding energy. We rank the functional groups presented in this work based on the strength of the resulting interface and suggest broad guidelines for the rational design of nanotube functionalisation for nanotube-polymer composites.

Corrigendum: Recent progress in linear-scaling density functional calculations with plane waves and pseudopotentials: the ONETEP code (2008 J. Phys.: Condens. Matter 20 064209).

N/A.

Corrigendum: Density kernel optimization in the ONETEP code (2008 J. Phys.: Condens. Matter 20 294207).

N/A.

Corrigendum: Using ONETEP for accurate and efficient O(N) density functional calculations (2005 J. Phys.: Condens. Matter 17 5757).

N/A.

Combining Embedded Mean-Field Theory with Linear-Scaling Density-Functional Theory.

We demonstrate the capability of embedded mean-field theory (EMFT) within the linear-scaling density-functional-theory code ONETEP, which enables DFT-in-DFT quantum embedding calculations on systems containing thousands of atoms at a fraction of the cost of a full calculation. We perform simulations on a wide range of systems from molecules to complex nanostructures to demonstrate the performance of our implementation with respect to accuracy and efficiency. This work paves the way for the application of this class of quantum embedding method to large-scale systems that are beyond the reach of existing implementations.

Wannier90 as a community code: new features and applications.

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.

Negative thermal expansion in high pressure layered perovskite Ca2GeO4.

We report the high pressure synthesis of a layered perovskite Ca2GeO4 which is found to have the Ruddlesden-Popper structure with I41/acd symmetry. Consonant with our recent predictions [Ablitt et al., npj Comput. Mater., 2017, 3, 44], the phase displays pronounced uniaxial negative thermal expansion over a large temperature range. Negative thermal expansion that persists over a large temperature range is very unusual in the perovskite structure, and its occurrence in this instance can be understood to arise due to both soft lattice vibrations associated with a phase competition and the unusually compliant nature of this structure, which effectively couples thermal expansion in the layer plane to lattice contractions perpendicular to the layering direction via a "corkscrew" mechanism.

Multiscale simulations of critical interfacial failure in carbon nanotube-polymer composites.

Computational investigation of interfacial failure in composite materials is challenging because it is inherently multi-scale: the bond-breaking processes that occur at the covalently bonded interface and initiate failure involve quantum mechanical phenomena, yet the mechanisms by which external stresses are transferred through the matrix occur on length and time scales far in excess of anything that can be simulated quantum mechanically. In this work, we demonstrate and validate an adaptive quantum mechanics (QM)/molecular mechanics simulation method that can be used to address these issues and apply it to study critical failure at a covalently bonded carbon nanotube (CNT)-polymer interface. In this hybrid approach, the majority of the system is simulated with a classical forcefield, while areas of particular interest are identified on-the-fly and atomic forces in those regions are updated based on QM calculations. We demonstrate that the hybrid method results are in excellent agreement with fully QM benchmark simulations and offers qualitative insights missing from classical simulations. We use the hybrid approach to show how the chemical structure at the CNT-polymer interface determines its strength, and we propose candidate chemistries to guide further experimental work in this area.

Control of Uniaxial Negative Thermal Expansion in Layered Perovskites by Tuning Layer Thickness.

Uniaxial negative thermal expansion (NTE) is known to occur in low n members of the A n+1B n O3n+1 Ruddlesden-Popper (RP) layered perovskite series with a frozen rotation of BO6 octahedra about the layering axis. Previous work has shown that this NTE arises due to the combined effects of a close proximity to a transition to a competing phase, so called "symmetry trapping", and highly anisotropic elastic compliance specific to the symmetry of the NTE phase. We extend this analysis to the broader RP family (n = 1, 2, 3, 4, …, ∞), demonstrating that by changing the fraction of layer interface in the structure (i.e., the value of 1/n) one may control the anisotropic compliance that is necessary for the pronounced uniaxial NTE observed in these systems. More detailed analysis of how the components of the compliance matrix develop with 1/n allows us to identify different regimes, linking enhancements in compliance between these regimes to the crystallographic degrees of freedom in the structure. We further discuss how the perovskite layer thickness affects the frequencies of soft zone boundary modes with large negative Grüneisen parameters, associated with the aforementioned phase transition, that constitute the thermodynamic driving force for NTE. This new insight complements our previous work-showing that chemical control may be used to switch from positive to negative thermal expansion in these systems-since it makes the layer thickness, n, an additional design parameter that may be used to engineer layered perovskites with tuneable thermal expansion. In these respects, we predict that, with appropriate chemical substitution, the n = 1 phase will be the system in which the most pronounced NTE could be achieved.

Tuning electronic properties of transition-metal dichalcogenides via defect charge.

Defect engineering is a promising route for controlling the electronic properties of monolayer transition-metal dichalcogenide (TMD) materials. Here, we demonstrate that the electronic structure of MoS2 depends sensitively on the defect charge, both its sign and magnitude. In particular, we study shallow bound states induced by charged defects using large-scale tight-binding simulations with screened defect potentials and observe qualitative changes in the orbital character of the lowest lying impurity states as function of the impurity charge. To gain further insights, we analyze the competition of impurity states originating from different valleys of the TMD band structure using effective mass theory and find that impurity state binding energies are controlled by the effective mass of the corresponding valley, but with significant deviations from hydrogenic behaviour due to unconventional screening of the defect potential.