Search

Interlayer rotation angle couples strongly to the electronic states of twisted van der Waals layers. However, not every angle is energetically favorable. Recent experiments on rotation-tunable electronics reveal the existence of a discrete set of angles at which the rotation-tunable electronics assume the most stable configurations. Nevertheless, a quantitative map for locating these intrinsically preferred twist angles in twisted bilayer system has not been available, posing challenges for the on-demand design of twisted electronics that are intrinsically stable at desired twist angles. Here we reveal a simple mapping between intrinsically preferred twist angles and geometry of the twisted bilayer system, in the form of geometric scaling laws for a wide range of intrinsically preferred twist angles as a function of only geometric parameters of the rotating flake on a supporting layer. We reveal these scaling laws for triangular and hexagonal flakes since they frequently appear in chemical vapor deposition growth. We also present a general method for handling arbitrary flake geometry. Such dimensionless scaling laws possess universality for all kinds of two-dimensional material bilayer systems, providing abundant opportunities for the on-demand design of intrinsic "twistronics". For example, the set of increasing magic-sizes that intrinsically prefers zero-approaching sequence of multiple magic-angles in bilayer graphene system can be revealed.

Germanium manganese compounds exhibit a variety of stable and metastable phases with different stoichiometry. These materials entail interesting electronic, magnetic and thermal properties both in their bulk form and as heterostructures. Here we develop and validate a transferable machine learning potential, based on the high-dimensional neural network formalism, to enable the study of Mn$_x$Ge$_y$ materials over a wide range of compositions. We show that a neural network potential fitted on a minimal training set reproduces successfully the structural and vibrational properties and the thermal conductivity of systems with different local chemical environments, and it can be used to predict phononic effects in nanoscale heterostructures., Comment: 29 pages, 9 figures

We report on the structural and transport properties of the smallest dislocation loop in graphene, known as a flower defect. First, by means of advanced experimental imaging techniques, we deduce how flower defects are formed during recrystallization of chemical vapor deposited graphene. We propose that the flower defects arise from a bulge type mechanism in which the flower domains are the grains left over by dynamic recrystallisation. Next, in order to evaluate the use of such defects as possible building blocks for all-graphene electronics, we combine multiscale modeling tools to investigate the structure and the electron and phonon transport properties of large monolayer graphene samples with a random distribution of flower defects. For large enough flower densities, we find that electron transport is strongly suppressed while, surprisingly, hole transport remains almost unaffected. These results suggest possible applications of flowered graphene for electron energy filtering. For the same defect densities, phonon transport is reduced by orders of magnitude as elastic scattering by defects becomes dominant. Heat transport by flexural phonons, key in graphene, is largely suppressed even for very low concentrations., Comment: 12 pages, 7 figures

The structure of the CiCs complex in silicon has long been the subject of debate. Numerous theoretical and experimental studies have attempted to shed light on the properties of these defects that are at the origin of the light emitting G-center. These defects are relevant for applications in lasing, and it would be advantageous to control their formation and concentration in bulk silicon. It is therefore essential to understand their structural and electronic properties. In this paper, we present the structural, electronic, and optical properties of four possible configurations of the CiCs complex in bulk silicon, namely the A-, B-, C-, and D-forms. The configurations were studied by density functional theory (DFT) and many-body perturbation theory (MBPT). Our results suggest that the C-form was misinterpreted as a B-form in some experiments. Our optical investigation also tends to exclude any contribution of A- and B-forms to light emission. Taken together, our results suggest that the C-form could play an important role in heavily carbon-doped silicon.

In this letter, we study the stability of the domain model for lithium intercalated graphite in stages III and II by means of Density Functional Theory and Kinetic Lattice Monte Carlo simulations. We find that the domain model is either thermodynamically or kinetically stable when compared to the standard model in stages III and II. The existence of domains in the intercalation sequence is well supported by recent high resolution transmission electron microscope observations in lithiated graphite. Moreover, we predict that such domain staging sequences leads to a wide range of diffusivity as reported in experiments.

By means of Density Functional Theory calculations we evaluate several lithium carbonate - graphite interface models as a prototype of the Solid Electrolyte Interphase capping layer on graphite anodes in lithium-ion batteries. It is found that only an (a,b)-oriented Li2CO3 slab promotes tight binding with graphite. Such mutual organization of the components combines their structural features and reproduces coordination environment of ions, resulting in an adhesive energy of 116 meV/{\AA}2 between graphite and lithium carbonate. This model also presents a high potential affinity with bulk. The corresponding charge distribution at such interface induces an electric potential gradient, such a gradient having been experimentally observed. We regard the mentioned criteria as the key descriptors of the interface stability and recommend them as the principal assessments for such interface study. In addition, we evaluate the impact of lithiated graphite on the stability of the model interface and study the generation of different point defects as mediators for Li interface transport. It is found that Li diffusion is mainly provided by interstitials. The induced potential gradient fundamentally assists the intercalation up to lithiation ratio of 70%., Comment: 24 pages, 2 figures

The structure of the ultimately-thin crystalline allotrope of silicon oxide, prepared onto a ruthenium surface, is unveiled down to atomic scale with chemical sensitivity, thanks to high resolution scanning tunneling microscopy and first principle calculations. An ordered oxygen lattice is imaged which coexists with the two-dimensional monolayer oxide. This coexistence signals a displacive transformation from an oxygen reconstructed-Ru(0001) to silicon oxide, along which latterally-shifted domains form, each with equivalent and degenerate epitaxial relationships with the substrate. The unavoidable character of defects at boundaries between these domains appeals for the development of alternative methods capable of producing single-crystalline two-dimensional oxides.

We have investigated the variation in the magnetization of highly ordered pyrolytic graphite (HOPG) after neutron irradiation, which introduces defects in the bulk sample and consequently gives rise to a large magnetic signal. We observe strong paramagnetism in HOPG, increasing with the neutron fluence. We correlate the induced paramagnetism with structural defects by comparison with density-functional theory calculations. In addition to the in-plane vacancies, the trans-planar defects also contribute to the magnetization. The lack of any magnetic order between the local moments is possibly due to the absence of hydrogen/nitrogen chemisorption, or the magnetic order cannot be established at all in the bulk form., Comment: 10 pages, 8 figures

The efficiency of minimum-energy configuration searching algorithms is closely linked to the energy landscape structure of complex systems. Here we characterize this structure by following the time evolution of two systems, vacancy aggregation in Fe and energy relaxation in ion-bombarded c-Si, using the kinetic Activation-Relaxation Technique (k-ART), an off-lattice kinetic Monte Carlo (KMC) method, and the well-known Bell-Evans-Polanyi (BEP) principle. We also compare the efficiency of two methods for handling non-diffusive flickering states -- an exact solution and a Tabu-like approach that blocks already visited states. Comparing these various simulations allow us to confirm that the BEP principle does not hold for complex system since forward and reverse energy barriers are completely uncorrelated. This means that following the lowest available energy barrier, even after removing the flickering states, leads to rapid trapping: relaxing complex systems requires crossing high-energy barriers in order to access new energy basins, in agreement with the recently proposed replenish-and-relax model [B\'eland et al., PRL 111, 105502 (2013)] This can be done by forcing the system through these barriers with Tabu-like methods. Interestingly, we find that following the fundamental kinetics of a system, though standard KMC approach, is at least as efficient as these brute-force methods while providing the correct kinetics information., Comment: 10 pages, 15 figures

We investigate Ge mixing at the Si(001) surface and characterize the $2\times N$ Si(001) reconstruction by means of hybrid quantum and molecular mechanics calculations (QM/MM). Avoiding fake elastic dampening, this scheme allows to correctly take into account long range deformation induced by reconstruted and defective surfaces. We focus in particular on the dimer vacancy line (DVL) and its interaction with Ge adatoms. We first show that calculated formation energies for these defects are highly dependent on the choice of chemical potential and that the latter must be chosen carefully. Characterizing the effect of the DVL on the deformation field, we also find that the DVL favors Ge segregation in the fourth layer close to the DVL. Using the activation-relaxation technique (ART nouveau) and QM/MM, we show that a complex diffusion path permits the substitution of the Ge atom in the fourth layer, with barriers compatible with mixing observed at intermediate temperature., Comment: 11 pages, 7 figures, 3 tables

In this letter, we study enhanced-retarded diffusion of oxygen in doped silicon by means of first principle calculations. We found that the migration of oxygen dimers can not be significantly affected by strain, doping type or rate. We attribute the enhanced oxygen diffusion in p-doped silicon to reduced monomer migration energy, and the retarded oxygen diffusion in Sb-doped to oxygen trapping close to a dopant site. The two proposed kinetic and thermodynamic mechanisms can appear at the same time and might lead to contradictory experimental results. Such mechanisms can be involved in the light induced degradation phenomenon in solar grade silicon.

We present novel structural motifs for boron-carbon nano-cages of the stochiometries B12C48 and B12C50, based on first principle calculations. These configurations are distinct from those proposed so far by the fact that the boron atoms are not isolated and distributed over the entire surface of the cages, but rather aggregate at one location to form a patch. Our putative ground state of B12C48 is 1.8 eV lower in energy than the previously proposed ground state and violates all the suggested empirical rules for constructing low energy fullerenes. The B12C50 configuration is energetically even more favorable than B12C48, showing that structures derived from the C60 Buckminsterfullerene are not necessarily magic sizes for heterofullerenes.

So far, no boron fullerenes were synthesized: more compact sp3-bonded clusters are energetically preferred. To circumvent this, metallic clusters have been suggested by Pochet et al. [Phys. Rev. B 83, 081403(R) (2011)] as "seeds" for a possible synthesis which would topologically protect the sp2 sector of the configuration space. In this paper, we identify a basic pentagonal unit which allows a balance between the release of strain and the self-doping rule. We formulate a guiding principle for the stability of boron fullerenes, which takes the form of an isolated filled pentagon rule (IFPR). The role of metallic clusters is then reexamined. It is shown that the interplay of the IFPR and the seed-induced doping breaks polymorphism and its related problems: it can effectively select between different isomers and reduce the reactivity of the boron shells. The balance between self and exterior doping represents the best strategy for boron buckyball synthesis., Comment: 15 pages, 3 figures; submitted

X-ray absorption (XAS) and x-ray magnetic circular dichroism (XMCD) spectra at the L$_{2,3}$ edges of Mn in (Ge,Mn) compounds have been measured and are compared to the results of first principles calculation. Early \textit{ab initio} studies show that the Density Functional Theory (DFT) can very well describe the valence band electronic properties but fails to reproduce a characteristic change of sign in the L$_{3}$ XMCD spectrum of Mn in Ge$_3$Mn$_5$, which is observed in experiments. In this work we demonstrate that this disagreement is partially related to an underestimation of the exchange splitting of Mn 2$p$ core states within the local density approximation. It is shown that the change in sign experimentally observed is reproduced if the exchange splitting is accurately calculated within the Hartree-Fock approximation, while the final states can be still described by the DFT. This approach is further used to calculate the XMCD in different (Ge,Mn) compounds. It demonstrates that the agreement between experimental and theoretical spectra can be improved by combining state of the art calculations for the core and valence states respectively., Comment: 8 pages

Magnetism in 2D atomic sheets has attracted considerable interest as its existence could allow the development of electronic and spintronic devices. The existence of magnetism is not sufficient for devices, however, as states must be addressable and modifiable through the application of an external drive. We show that defects in hexagonal boron nitride present a strong interplay between the the N-N distance in the edge and the magnetic moments of the defects. By stress-induced geometry modifications, we change the ground state magnetic moment of the defects. This control is made possible by the triangular shape of the defects as well as the strong spin localisation in the magnetic state., Comment: 10 pages, 3 figures, published in APL

We use extensive first principle simulations to show the major role played by interfaces in the mechanism of phase separation observed in semiconductor multifunctional materials. We make an analogy with the precipitation sequence observed in over-saturated AlCu alloys, and replace the Guinier-Preston zones in this new context. A new class of materials, the $\alpha$ phases, is proposed to understand the formation of the coherent precipitates observed in the GeMn system. The interplay between formation and interface energies is analyzed for these phases and for the structures usually considered in the literature. The existence of the alpha phases is assessed with both theoretical and experimental arguments.

Using the minima hopping global geometry optimization method on the density functional potential energy surface we study medium size and large boron clusters. Even though for isolated medium size clusters the ground state is a cage like structure they are unstable against external perturbations such as contact with other clusters. The energy landscape of larger boron clusters is glass like and has a large number of structures which are lower in energy than the cages. This is in contrast to carbon and boron nitride systems which can be clearly identified as structure seekers in our minima hopping runs. The differences in the potential energy landscape explain why carbon and boron nitride systems are found in nature whereas pure boron fullerenes have not been found., Comment: 5 pages 5 figures