Abstract: Using a numerical tight-binding approach based on the Chebyshev–Bogoliubov–de Gennes method we describe Josephson junctions made of multilayer graphene contacted by top superconducting gates. Both Bernal (ABA) and rhombohedral (ABC) stacking are considered and we find that the type of stacking has a strong effect on the proximity effect and the supercurrent flow. For both cases the pair amplitude shows a polarization between dimer and nondimer atoms, being more pronounced for rhombohedral stacking. Even though the proximity effect in nondimer sites is enhanced when compared to single-layer graphene, we find that the supercurrent is suppressed. The spatial distribution of the supercurrent shows that for Bernal stacking the current flows only in the topmost layers while for rhombohedral stacking the current flows throughout the whole structure.

Abstract: The magneto- and Hall resistance of a locally gated cross shaped graphene Hall bar is calculated. The edge of the top gate is placed diagonally across the center of the Hall cross. Four-probe resistance is calculated using the Landauer-Büttiker formalism, while the transmission coefficients are obtained using the non-equilibrium Green's function approach. The interplay between transport due to edge channels and snake states is investigated. When two edge channels are occupied, we predict oscillations in the Hall and the bend resistance as function of the magnetic field, which are a consequence of quantum interference between the occupied snake states.

Abstract: We investigate the magnetic field dependence of the Hall and the bend resistances for a ballistic Hall bar structure containing a pn-junction sculptured from a bilayer of graphene. The electric response is obtained using the billiard model, and we investigate the cases of bilayer graphene with and without a band gap. Two different conduction regimes are possible: (i) both sides of the junction have the same carrier type and (ii) one side of the junction is n-type while the other one is p-type. The first case shows Hall plateau-like features in the Hall resistance that fade away as the band gap opens. The second case exhibits a bend resistance that is asymmetric in magnetic field as a consequence of snake states along the pn-interface, where the maximum is shifted away from zero magnetic field.

Abstract: Transverse magnetic focusing properties of graphene using a ballistic four terminal structure are investigated. The electric response is obtained using the semiclassical billiard model. The transmission exhibits pronounced peaks as a consequence of skipping orbits at the edge of the structure. When we add a pn-junction between the two probes, snake states along the pn-interface appear. Injected electrons are guided by the pn-interface to one of the leads depending on the value of the applied magnetic field. Oscillations in the resistance are found depending on the amount of particles that end up in each lead.

Abstract: Hydrogen adsorption on boron and nitrogen doped graphene is investigated in detail by means of first-principles calculations. A comprehensive study is performed of the structural, electronic, and magnetic properties of chemisorbed hydrogen atoms and atom pairs near the dopant sites. The main effect of the substitutional atoms is charge doping which is found to greatly affect the adsorption process by increasing the binding energy at the sites closest to the substitutional species. It is also found that doping does not induce magnetism despite the odd number of electrons per atom introduced by the foreign species, and that it quenches the paramagnetic response of chemisorbed H atoms on graphene. Overall, the effects are similar for B and N doping, with only minor differences in the adsorption energetics due to different sizes of the dopant atoms and the accompanying lattice distortions.

Abstract: The properties of a new class of phosphorene nanotubes (PNT) are investigated by performing first-principles calculations. We demonstrate that it is advantageous to use blue phosphorene in order to make small nanotubes and propose a way to create low-energy PNTs by the inclusion of defect lines. Five different types of defect lines are investigated and incorporated in various combinations. The resulting defect-induced faceted PNTs have negligible bending stresses which leads to a reduction in the formation energy with respect to round PNTs. Our armchair faceted PNTs have similar formation energies than the recently proposed multiphase faceted PNTs, but they have a larger variety of possible structures. Our zigzag faceted PNTs have lower formation energies than round tubes and multiphase faceted nanotubes. The electronic properties of the defect-induced faceted PNTs are determined by the defect lines which control the band gap and the shape of the electronic states at the band edges. These band gaps increase with the radius of the nanotubes and converge to those of isolated defect lines.

Abstract: Thermal fluctuations play an increasingly important role in micromagnetic research relevant for various biomedical and other technological applications. Until now, it was deemed necessary to use a time stepping algorithm with a fixed time step in order to perform micromagnetic simulations at nonzero temperatures. However, Berkov and Gorn have shown in [D. Berkov and N. Gorn, J. Phys.: Condens. Matter,14, L281, 2002] that the drift term which generally appears when solving stochastic differential equations can only influence the length of the magnetization. This quantity is however fixed in the case of the stochastic Landau-Lifshitz-Gilbert equation. In this paper, we exploit this fact to straightforwardly extend existing high order solvers with an adaptive time stepping algorithm. We implemented the presented methods in the freely available GPU-accelerated micromagnetic software package MuMax3 and used it to extensively validate the presented methods. Next to the advantage of having control over the error tolerance, we report a twenty fold speedup without a loss of accuracy, when using the presented methods as compared to the hereto best practice of using Heun’s solver with a small fixed time step.

Abstract: Magnetic skyrmions are topological spin configurations in materials with chiral Dzyaloshinskii-Moriya interaction (DMI), that are potentially useful for storing or processing information. To date, DMI has been found in few bulk materials, but can also be induced in atomically thin magnetic films in contact with surfaces with large spin-orbit interactions. Recent experiments have reported that isolated magnetic skyrmions can be stabilized even near room temperature in few-atom-thick magnetic layers sandwiched between materials that provide asymmetric spin-orbit coupling. Here we present the minimum-energy path analysis of three distinct mechanisms for the skyrmion collapse, based on ab initio input and the performed atomic-spin simulations. We focus on the stability of a skyrmion in three atomic layers of Co, either epitaxial on the Pt(111) surface or within a hybrid multilayer where DMI nontrivially varies per monolayer due to competition between different symmetry breaking from two sides of the Co film. In laterally finite systems, their constrained geometry causes poor thermal stability of the skyrmion toward collapse at the boundary, which we show to be resolved by designing the high-DMI structure within an extended film with lower or no DMI.

Abstract: Van der Waals (vdW) heterostructures continue to attract intense interest as a route of designing materials with novel properties that cannot be found in nature. Unfortunately, this approach is currently limited to only a few layers that can be stacked on top of each other. Here, we report a bulk vdW material consisting of superconducting 1H TaS2 monolayers interlayered with 1T TaS2 monolayers displaying charge density waves (CDW). This bulk vdW heterostructure is created by phase transition of 1T-TaS2 to 6R at 800 degrees C in an inert atmosphere. Its superconducting transition (T-c) is found at 2.6 K, exceeding the T-c of the bulk 2H phase. Using first-principles calculations, we argue that the coexistence of superconductivity and CDW within 6R-TaS2 stems from amalgamation of the properties of adjacent 1H and 1T monolayers, where the former dominates the superconducting state and the latter the CDW behavior.

Abstract: The pyroantimonate pigments Naples yellow and lead tin antimonate yellow are recognized as some of the most stable synthetic yellow pigments in the history of art. However, this exceptional lightfastness is in contrast with experimental evidence suggesting that this class of mixed oxides is of semiconducting nature. In this study the electronic structure and light-induced behavior of the lead pyroantimonate pigments were determined by means of a combined multifaceted analytical and computational approach (photoelectrochemical measurements, UV-vis diffuse reflectance spectroscopy, STEM-EDS, STEM-HAADF, and density functional theory calculations). The results demonstrate both the semiconducting nature and the lightfastness of these pigments. Poor optical absorption and minority carrier mobility are the main properties responsible for the observed stability. In addition, novel fundamental insights into the role played by Na atoms in the stabilization of the otherwise intrinsically unstable Pb2Sb2O7 pyrochlore were obtained.

Abstract: Layered Li-rich oxides, demonstrating both cationic and anionic redox chemistry being used as positive electrodes for Li-ion batteries,have raised interest due to their high specific discharge capacities exceeding 250 mAh/g. However, irreversible structural transformations triggered by anionic redox chemistry result in pronounced voltagefade (i.e., lowering the specific energy by a gradual decay of discharge potential) upon extended galvanostatic cycling. Activating or suppressing oxygen anionic redox through structural stabilization induced by redox-inactivecation substitution is a well-known strategy. However, less emphasishas been put on the correlation between substitution degree and theactivation/suppression of the anionic redox. In this work, Ti4+-substituted Li2MnO3 was synthesizedvia a facile solution-gel method. Ti4+ is selected as adopant as it contains no partially filled d-orbitals. Our study revealedthat the layered “honeycomb-ordered” C2/m structure is preserved when increasing the Ticontent to x = 0.2 in the Li2Mn1-x Ti (x) O-3 solidsolution, as shown by electron diffraction and aberration-correctedscanning transmission electron microscopy. Galvanostatic cycling hintsat a delayed oxygen release, due to an improved reversibility of theanionic redox, during the first 10 charge-discharge cyclesfor the x = 0.2 composition compared to the parentmaterial (x = 0), followed by pronounced oxygen redoxactivity afterward. The latter originates from a low activation energybarrier toward O-O dimer formation and Mn migration in Li2Mn0.8Ti0.2O3, as deducedfrom first-principles molecular dynamics (MD) simulations for the“charged” state. Upon lowering the Ti substitution to x = 0.05, the structural stability was drastically improvedbased on our MD analysis, stressing the importance of carefully optimizingthe substitution degree to achieve the best electrochemical performance.

Abstract: Chalcopyrite semiconductors are of considerable interest for application as absorber layers in thin-film photovoltaic cells. When growing films of these compounds, however, they are often found to contain CuAu-like domains, a metastable phase of chalcopyrite. It has been reported that for CuInS2, the presence of the CuAu-like phase improves the short circuit current of the chalcopyrite-based photovoltaic cell. We investigate the thermodynamic stability of both phases for a selected list of I-III-VI2 materials using a first-principles density functional theory approach. For the CuIn-VI2 compounds, the difference in formation energy between the chalcopyrite and CuAu-like phase is found to be close to 2 meV per atom, indicating a high likelihood of the presence of CuAu-like domains. Next, we calculate the spectroscopic limited maximum efficiency (SLME) of the CuAu-like phase and compare the results with those of the corresponding chalcopyrite phase. We identify several candidates with a high efficiency, such as CuAu-like CuInS2, for which we obtain an SLME of 29% at a thickness of 500 nm. We observe that the SLME can have values above the Shockley-Queisser (SQ) limit, and show that this can occur because the SQ limit assumes the absorptivity to be a step function, thus overestimating the radiative recombination in the detailed balance approach. This means that it is possible to find higher theoretical efficiencies within this framework simply by calculating the J-V characteristic with an absorption spectrum. Finally, we expand our SLME analysis to indirect band gap absorbers by studying silicon, and find that the SLME quickly overestimates the reverse saturation current of indirect band gap materials, drastically lowering their calculated efficiency.

Abstract: Cu-based chalcogenides are promising materials for thin-film solar cells with more than 20% measured
cell efficiency. Using first-principles calculations based on density functional theory, the
optoelectronic properties of a group of Cu-based chalcogenides Cu2-II-IV-VI4 is studied. They are
then screened with the aim of identifying potential absorber materials for photovoltaic applications.
The spectroscopic limited maximum efficiency (SLME) introduced by Yu and Zunger [Phys. Rev.
Lett. 108, 068701 (2012)] is used as a metric for the screening. After constructing the currentvoltage
curve, the SLME is calculated from the maximum power output. The role of the nature of
the band gap, direct or indirect, and also of the absorptivity of the studied materials on the maximum
theoretical power conversion efficiency is studied. Our results show that Cu2II-GeSe4 with
II¼ Cd and Hg, and Cu2-II-SnS4 with II ¼ Cd, Hg, and Zn have a higher theoretical efficiency
compared with the materials currently used as absorber layer.

Abstract: Lithium-rich layered oxides such as Li₂MnO₃have shown great potential as cathodes in Li-ion batteries, mainly because of their large capacities. However, these materials still suffer from structural degradation as the battery is cycled, reducing the average voltage and capacity of the cell. The voltage fade is believed to be related to the migration of transition metals into the lithium layer, linked to the formation of O-O dimers with a short bond length, which in turn is driven by the presence of oxygen holes due to the participation of oxygen in the redox process. We investigate the formation of O-O dimers for partially charged O1-Li₂MnO₃using a first-principles density functional theory approach by calculating the reaction energy and kinetic barriers for dimer formation. Next, we perform similar calculations for partially charged O1-Li₂IrO₃, a Li-rich material for which the voltage fade was not observed during cycling. When we compare the stability of the oxygen framework, we conclude that the formation of O-O dimers is both thermodynamically and kinetically viable for O1-Li_0.5MnO₃. For O1-Li_0.5IrO₃, we observe that the oxygen lattice is much more stable, either returning to its original state when perturbed, or resulting in a structure with an O-O dimer that is much higher in energy. This can be explained by the mixed redox process for Li₂IrO₃, which is also shown from the calculated magnetic moments. The lack of O-O dimer formation in O1-Li_0.5IrO₃provides valuable insight as to why Li₂IrO₃does not demonstrate a voltage fade as the battery is cycled, which can be used to design Li-rich battery cathodes with an improved cycling performance.

Abstract: We investigate the role of transition metal atoms of group V-b (V, Nb, Ta) and VI-b (Cr, Mo, W) as n- or p-type dopants in anatase TiO$_2$ using thermodynamic
principles and density functional theory with the Heyd-Scuseria-Ernzerhof HSE06 hybrid functional. The HSE06 functional provides a realistic value for the band gap, which ensures a correct classification of dopants as shallow or deep donors or acceptors. Defect formation energies and thermodynamic transition levels are calculated taking into account the constraints imposed by the stability of TiO$_2$ and the solubility limit of the impurities.
Nb, Ta, W and Mo are identified as shallow donors. Although W provides two electrons, Nb and Ta show a considerable lower formation energy, in particular under O-poor conditions. Mo donates in principle one electron, but under specific conditions can turn into a double donor. V impurities are deep donors and Cr
shows up as an amphoteric defect, thereby acting as an electron trapping center in n-type TiO$_2$ especially under O-rich conditions. A comparison with the available experimental data yields excellent agreement.