Abstract: Themean-field Hubbard model is used to investigate the formation of the antiferromagnetic phase in hexagonal graphene rings with inner zigzag edges. The outer edge of the ring was taken to be either zigzag or armchair, and we found that both types of structures can have a larger antiferromagnetic interaction as compared with hexagonal dots. This difference could be partially ascribed to the larger number of zigzag edges per unit area in rings than in dots. Furthermore, edge states localized on the inner ring edge are found to hybridize differently than the edge states of dots, which results in important differences in the magnetism of graphene rings and dots. The largest staggered magnetization is found when the outer edge has a zigzag shape. However, narrow rings with armchair outer edge are found to have larger staggered magnetization than zigzag hexagons. The edge defects are shown to have the least effect on magnetization when the outer ring edge is armchair shaped. DOI: 10.1103/PhysRevB.87.085434

Abstract: Thermal fluctuations of single layer hydrogenated graphene (graphane) are investigated using large scale atomistic simulations. By analyzing the mean square value of the height fluctuations < h(2)> and the height-height correlation function H(q) for different system sizes and temperatures, we show that hydrogenated graphene is an unrippled system in contrast to graphene. The height fluctuations are bounded, which is confirmed by a H(q) tending to a constant in the long wavelength limit instead of showing the characteristic scaling law q(4-eta)(eta similar or equal to 0.85) predicted by membrane theory. This unexpected behavior persists up to temperatures of at least 900 K and is a consequence of the fact that in graphane the thermal energy can be accommodated by in-plane bending modes, i.e., modes involving C-C-C bond angles in the buckled carbon layer, instead of leading to significant out-of-plane fluctuations that occur in graphene.

Abstract: Thin ferromagnetic films with an interfacially induced DMI exhibit nontrivial asymmetric dispersion relations that lead to unique and useful magnonic properties. Here we derive an analytical expression for the magnon propagation angle within the micromagnetic framework and show how the dispersion relation can be approximated with a comprehensible geometrical interpretation in the k space of the propagation of spin waves. We further explore the refraction of spin waves at DMI interfaces in heterochiral magnetic films, after deriving a generalized Snell's law tunable by an in-plane magnetic field, that yields analytical expressions for critical incident angles. The found asymmetric Brewster angles at interfaces of regions with different DMI strengths, adjustable by magnetic field, support the conclusion that heterochiral ferromagnetic structures are an ideal platform for versatile spin-wave guides.

Abstract: This work provides values of electron scattering 002 structure factors for InxGa1-xAs as a function of the In concentration x=0 to 1. These results allow accurate compositional analysis of pseudomorphically grown InxGa1-xAs/GaAs layers by transmission electron microscopy methods relying on the chemical sensitivity of the (002) beam. The calculations go beyond the limits of the isolated atom approximation, because they take into account charge redistribution effects between atomic sites in the crystal, strain, and static atomic displacements. The computations were performed by the full potential linearized augmented plane-wave method using a generalized gradient approximation for the exchange and correlation part of the potential. The calculations of strained InxGa1-xAs correspond to the strain state in specimens with large, small, and intermediate thickness in the electron beam direction. Additionally, the effect of static atomic displacements is taken into account. All results are listed in a parameterized form. The calculated 002 structure factor vanishes at an In concentration of 16.4%. This value is in a good agreement with previously reported experimental measurements. Hence, our results are a significant improvement with respect to the isolated atom approximation which is conventionally applied in transmission electron microscopy simulations, and which predicts a value of 22.5%.

Abstract: To date, several experimental groups reported measurements of the thickness dependence of T-c of atomically uniform single-crystalline Pb nanofilms. The reported amplitude of the T-c oscillations varies significantly from one experiment to another. Here we propose that the reason for this unresolved issue is an interplay of the quantum-size variations in the single-electron density of states with thickness-dependent oscillations in the phonon-mediated electron-electron coupling. Such oscillations in the coupling depend on the substrate material, the quality of the interface, the protection cover, and other details of the fabrication process, changing from one experiment to another. This explains why the available data do not exhibit one-voice consistency about the amplitude of the T-c oscillations. Our analyses are based on a numerical solution of the Bogoliubov-de Gennes equations for a superconducting slab.

Abstract: To unravel the effects of phonon confinement, the influence of size and morphology on the atomic vibrations is investigated in Sn nanoislands and cluster-assembled films. Nuclear resonant inelastic x-ray scattering is used to probe the phonon densities of states of the Sn nanostructures which show significant broadening of the features compared to bulk phonon behavior. Supported by ab initio calculations, the broadening is attributed to phonon scattering and can be described within the damped harmonic oscillator model. Contrary to the expectations based on previous research, the appearance of high-energy modes above the cutoff energy is not observed. From the thermodynamic properties extracted from the phonon densities of states, it was found that grain boundary Sn atoms are bound by weaker forces than bulk Sn atoms.

Abstract: Transition metal dichalcogenides (TMDs) are considered as promising candidates for next generation of electronic and optoelectronic devices. To make use of these materials, for instance in field effect transistor applications, it is mandatory to know the detailed properties of contacts of such TMDs with metal electrodes. Here, we investigate the role of the contact structure on the electronic and transport properties of metal-MoSe2 interfaces. Two different contact types, namely face and edge contacts, are studied. We consider both low (Sc) and high (Au) work function metals in order to thoroughly elucidate the role of the metal work function and the type of metal. First principles plane wave calculations and transport calculations based on nonequilibrium Green's function formalism reveal that the contact type has a large impact on the electronic and transport properties of metal-MoSe2 interfaces. For the Sc electrode, the Schottky barrier heights are around 0.25 eV for face contact and bigger than 0.6 eV for edge contact. For the Au case, we calculate very similar barrier heights for both contact types with an average value of 0.5 eV. Furthermore, while the face contact is found to be highly advantageous as compared to the edge contact for the Sc electrode, the latter contact becomes a better choice for the Au electrode. Our findings provide guidelines for the fabrication of TMD-based devices.

Abstract: Transport properties of a superconducting NbN nanowire are studied experimentally and theoretically. Different attached leads (superconducting contacts) allowed us to measure current-voltage (I-V) characteristics of different segments of the wire independently. The experimental results show that with increasing the length of the segment the number of jumps in the I-V curve increases indicating an increasing number of phase-slip phenomena. The system shows a clear hysteresis in the direction of the current sweep, the size of which depends on the length of the superconducting segment. The interpretation of the experimental results is supported by theoretical simulations that are based on the time-dependent Ginzburg-Landau theory, the heat equation has been included in the Ginzbur-Landau theory.

Abstract: Tunneling of holes through a trilayer structure made of two diluted magnetic semiconductors, (Ga,Mn)As, separated by a thin layer of nonmagnetic AlAs is investigated. The problem is treated within the 6x6 Luttinger-Kohn model for valence bands with the split-off band included. The influence of the spin-orbit coupling is pronounced as the spin-splitting Delta(ex) is comparable with the split-off Delta(SO) splitting. It is assumed that direct tunneling is the dominant mechanism due to the high quality of the tunnel junctions. Our theoretical results predict the correct order of magnitude for the tunneling magnetoresistance ratio, but various other effects, such as scattering on impurities and defects, should be included in order to realize a quantitative agreement with experiment.

Abstract: Two-dimensional graphene-based van der Waals heterostructures have received considerable interest because of their intriguing characteristics compared with the constituent single-layer two-dimensional materials. Here, we investigate the interfacial characteristics, Schottky contact, and optical performance of graphene/Ga2SSe van der Waals (vdW) heterostructure using first-principles calculations. The effects of stacking patterns, electric gating, and interlayer coupling on the interfacial properties of graphene/Ga2SSe heterostructures are also examined. Our results demonstrate that the Dirac cone of graphene is well preserved at the F point in all stacking patterns due to the weak vdW interactions, which keep the heterostructures feasible such that they can be obtained in further experiments. Moreover, depending on the stacking patterns, a small band gap of about 13-17 meV opens in graphene and has a high carrier mobility, indicating that the graphene/Ga2SSe heterostructures are potential candidates for future high-speed nanoelectronic applications. In the ground state, the graphene/Ga2SSe heterostructures form an n-type Schottky contact. The transformation from an n-type to a p-type Schottky contact or to an Ohmic contact can be forced by electric gating or by varying the interlayer coupling. Our findings could provide physical guidance for designing controllable Schottky nanodevices with high electronic and optical performances.

Abstract: Typical Raman spectra of transition-metal dichalcogenides (TMDs) display two prominent peaks, E-2g and A(1g), that are well separated from each other. We find that these modes are degenerate in bulk WSe2 yielding one single Raman peak in contrast to other TMDs. As the dimensionality is lowered, the observed peak splits in two. In contrast, our ab initio calculations predict that the degeneracy is retained even for WSe2 monolayers. Interestingly, for minuscule biaxial strain, the degeneracy is preserved, but once the crystal symmetry is broken by a small uniaxial strain, the degeneracy is lifted. Our calculated phonon dispersion for uniaxially strained WSe2 shows a good match to the measured Raman spectrum, which suggests that uniaxial strain exists in WSe2 flakes, possibly induced during the sample preparation and/or as a result of the interaction between WSe2 and the substrate. Furthermore, we find that WSe2 undergoes an indirect-to-direct band-gap transition from bulk to monolayers, which is ubiquitous for semiconducting TMDs. These results not only allow us to understand the vibrational and electronic properties of WSe2, but also point to effects of the interaction between the monolayer TMDs and the substrate on the vibrational and electronic properties. DOI: 10.1103/PhysRevB.87.165409

Abstract: Unidirectional chiral states are predicted in single layer graphene which originate from the breaking of the sublattice symmetry due to an asymmetric mass potential. The latter can be created experimentally using boron-nitride (BN) substrates with a line defect (B-B or N-N) that changes the induced mass potential in graphene. Solving the Dirac-Weyl equation, the obtained energy spectrum is compared with the one calculated using ab initio density functional calculations. We found that these one-dimensional chiral states are very robust and they can even exist in the presence of a small gap between the mass regions. In the latter case additional bound states are found that are topologically different from those chiral states.

Abstract: Using a particle-based model we examine the depinning motion of solitons in skyrmion chains in quasi -onedimensional (1D) and two-dimensional (2D) systems containing embedded 1D interfaces. The solitons take the form of a particle or hole in a commensurate chain of skyrmions. Under an applied drive, just above a critical depinning threshold, the soliton moves with a skyrmion Hall angle of zero. For higher drives, the entire chain depins, and in a 2D system we observe that both the solitons and chain move at zero skyrmion Hall angle and then transition to a finite skyrmion Hall angle as the drive increases. In a 2D system with a 1D interface that is at an angle to the driving direction, there can be a reversal of the sign of the skyrmion Hall angle from positive to negative. Our results suggest that solitons in skyrmion systems could be used as information carriers in racetrack geometries that would avoid the drawbacks of finite skyrmion Hall angles. The soliton states become mobile at significantly lower drives than the depinning transition of the skyrmion chains themselves.

Abstract: Using an effective one-dimensional cluster expansion in combination with first-principles electronic structure calculations we have studied the energetics and electronic properties of Sb-Te layered systems. For a Te concentration between 0 and 60 at. % an almost continuous series of metastable structures is obtained consisting of consecutive Sb bilayers next to consecutive Sb2Te3 units, with the general formula (Sb-2)(n)(Sb2Te3)(m) (n, m = 1,2, ... ). Between 60 and 100 at.% no stable structures are found. We account explicitly for the weak van derWaals bonding between Sb bilayers and Sb2Te3 units by using a recently developed functional, which strongly improves the interlayer bonding distances. At T = 0 K, no evidence is found for the existence of two separate single-phase regions delta and gamma and a two-phase region delta + gamma. Metastable compounds with a Te concentration between 0 and 40 at. % are semimetallic, whereas compounds with a Te concentration between 50 and 60 at. % are semiconducting. Compounds with an odd number of Sb layers are metallic and have a much higher formation energy than those with an even number of consecutive Sb layers, thereby favoring the formation of Sb bilayers.

Abstract: Using atomistic simulations we investigate the morphological properties of graphene deposited on top of a nanostructured substrate. Sinusoidally corrugated surfaces, steps, elongated trenches, one-dimensional and cubic barriers, spherical bubbles, Gaussian bumps, and Gaussian depressions are considered as support structures for graphene. The graphene-substrate interaction is governed by van der Waals forces and the profile of the graphene layer is determined by minimizing the energy using molecular dynamics simulations. Based on the obtained optimum configurations, we found that (i) for graphene placed over sinusoidally corrugated substrates with corrugation wavelengths longer than 2 nm, the graphene sheet follows the substrate pattern while for supported graphene it is always suspended across the peaks of the substrate, (ii) the conformation of graphene to the substrate topography is enhanced when increasing the energy parameter in the van der Waals model, (iii) the adhesion of graphene into the trenches depends on the width of the trench and on the graphene's orientation, i. e., in contrast to a small-width (3 nm) nanoribbon with armchair edges, the one with zigzag edges follows the substrate profile, (iv) atomic-scale graphene follows a Gaussian bump substrate but not the substrate with a Gaussian depression, and (v) the adhesion energy due to van der Waals interaction varies in the range [0.1-0.4] J/m(2).