Abstract: The formation of lattice defects and their effect on the electronic properties of single-layer FeCl2 are investigated by means of first-principles calculations. Among the vacancy defects, namely mono-, di-, and three-Cl vacancies and mono-Fe vacancy, the formation of mono-Cl vacancy is the most preferable. Comparison of two different antisite defects reveals that the formation of the Fe-antisite defect is energetically preferable to the Cl-antisite defect. While a single Cl vacancy leads to a 1 mu(B) decrease in the total magnetic moment of the host lattice, each Fe vacant site reduces the magnetic moment by 4 mu(B). However, adsorption of an excess Cl atom on the surface changes the electronic structure to a ferromagnetic metal or to a ferromagnetic semiconductor depending on the adsorption site without changing the ferromagnetic state of the host lattice. Both Cl-antisite and Fe-antisite defected domains change the magnetic moment of the host lattice by -1 mu(B) and +3 mu(B), respectively. The electronic ground state of defected structures reveals that (i) single-layer FeCl2 exhibits half-metallicity under the formation of vacancy and Cl-antisite defects; (ii) ferromagnetic metallicity is obtained when a single Cl atom is adsorbed on upper-Cl and Fe sites, respectively; and (iii) ferromagnetic semiconducting behavior is found when a Cl atom is adsorbed on a lower-Cl site or a Fe-antisite defect is formed. Simulated scanning electron microscope images show that atomic-scale identification of defect types is possible from their electronic charge density. Further investigation of the periodically Fe-defected structures reveals that the formation of the single-layer FeCl3 phase, which is a dynamically stable antiferromagnetic semiconductor, is possible. Our comprehensive analysis on defects in single-layer FeCl2 will complement forthcoming experimental observations.

Abstract: Atomic collapse in graphene nanoribbons behaves in a fundamentally different way as compared to monolayer graphene due to the presence of multiple energy bands and the effect of edges. For armchair nanoribbons we find that bound states gradually transform into atomic collapse states with increasing impurity charge. This is very different in zigzag nanoribbons where multiple quasi-one-dimensional bound states are found that originates from the zero-energy zigzag edge states. They are a consequence of the flat band and the electron distribution of these bound states exhibits two peaks. The lowest-energy edge state transforms from a bound state into an atomic collapse resonance and shows a distinct relocalization from the edge to the impurity position with increasing impurity charge.

Abstract: Electrons in alpha-T-3 lattices behave as condensed-matter analogies of integer-spin Dirac fermions. The three atoms making up the unit cell bestow the energy spectrum with an additional energy band that is completely flat, providing unique electronic properties. The interatomic hopping term, alpha, is known to strongly affect the electronic spectrum of the two-dimensional (2D) lattice, allowing it to continuously morph from graphenelike responses to the behavior of fermions in a dice lattice. For pristine lattice structures the energy bands are gapless, but small deviations in the atomic equivalence of the three sublattices will introduce gaps in the spectrum. It is unknown how these affect transport and electronic properties such as the energy spectrum of superlattice minibands. Here we investigate the dependency of these properties on the parameter a accounting for different symmetry-breaking terms, and we show how it affects band-gap formation. Furthermore, we find that superlattices can force band gaps to close and shift in energy. Our results demonstrate that alpha-T-3 superlattices provide a versatile material for 2D band-gap engineering purposes.

Abstract: We examine from first principles the flexoelectric properties of phosphorene nanoribbons under mechanical bending along armchair and zigzag directions. In both cases we find that the radial polarization depends linearly on the strain gradient. The flexoelectricity along the armchair direction is over 40% larger than along the zigzag direction. The obtained flexoelectric coefficients of phosphorene are four orders of magnitude larger than those of graphene and comparable to transition metal dichalcogenides. Analysis of charge density shows that the flexoelectricity mainly arises from the pz orbitals of phosphorus atoms. The electron mobilities in bent phosphorene can be enhanced by over 60% along the armchair direction, which is significantly higher than previous reports of mobility tuned by uniaxial strain. Our results indicate phosphorene is a candidate for a two-dimensional material applicable in flexible-electronic devices.

Abstract: Structural symmetry breaking in two-dimensional materials can lead to superior physical properties and introduce an additional degree of piezoelectricity. In the present paper, we propose three structural phases (1H, 1T, and 1T') of Janus WXO (X = S, Se, and Te) monolayers and investigate their vibrational, thermal, elastic, piezoelectric, and electronic properties by using first-principles methods. Phonon spectra analysis reveals that while the 1H phase is dynamically stable, the 1T phase exhibits imaginary frequencies and transforms to the distorted 1T' phase. Ab initio molecular dynamics simulations confirm that 1H- and 1T'-WXO monolayers are thermally stable even at high temperatures without any significant structural deformations. Different from binary systems, additional Raman active modes appear upon the formation of Janus monolayers. Although the mechanical properties of 1H-WXO are found to be isotropic, they are orientation dependent for 1T'-WXO. It is also shown that 1H-WXO monolayers are indirect band-gap semiconductors and the band gap narrows down the chalcogen group. Except 1T'-WSO, 1T'-WXO monolayers have a narrow band gap correlated with the Peierls distortion. The effect of spin-orbit coupling on the band structure is also examined for both phases and the alteration in the band gap is estimated. The versatile mechanical and electronic properties of Janus WXO monolayers together with their large piezoelectric response imply that these systems are interesting for several nanoelectronic applications.

Abstract: By using density functional theory (DFT) based first-principles calculations, electronic, vibrational, piezo-electric, and optical properties of monoclinic Janus single-layer Ga2TeX (X = S or Se) are investigated. The dynamical, mechanical, and thermal stability of the proposed Janus single layers are verified by means of phonon bands, stiffness tensor, and quantum molecular dynamics simulations. The calculated vibrational spectrum reveals the either pure or coupled optical phonon branches arising from Ga-Te and Ga-X atoms. In addition to the in-plane anisotropy, single-layer Janus Ga2TeX exhibits additional out-of-plane asymmetry, which leads to important consequences for its electronic and optical properties. Electronic band dispersions indicate the direct band-gap semiconducting nature of the constructed Janus structures with energy band gaps falling into visible spectrum. Moreover, while orientation-dependent linear-elastic properties of Janus single layers indicate their strong anisotropy, the calculated in-plane stiffness values reveal the ultrasoft nature of the structures. In addition, predicted piezoelectric coefficients show that while there is a strong in-plane anisotropy between piezoelectric constants along armchair (AC) and zigzag (ZZ) directions, there exists a tiny polarization along the out-of-plane direction as a result of the formation of Janus structure. The optical response to electromagnetic radiation has been also analyzed through density functional theory by considering the independent-particle approximation. Finally, the optical spectra of Janus Ga2TeX structures is investigated and it showed a shift from the ultraviolet region to the visible region. The fact that the spectrum is between these regions will allow it to be used in solar energy and many nanoelectronics applications. The predicted monoclinic single-layer Janus Ga2TeX are relevant for promising applications in optoelectronics, optical dichroism, and anisotropic nanoelasticity.

Abstract: The effect of the spatially varying spin current on a skyrmion is numerically investigated. It is shown that an inhomogeneous current density induces an elongation of the skyrmion. This elongation can be controlled using current pulses of different strength and duration. Long current pulses lead to a splitting that forms two replicas of the initial skyrmion while for short pulses the elongated skyrmion relaxes back to its initial circular state through rotation in the MHz-GHz frequency range. The frequency is dependent on the strength of the damping coefficient.

Abstract: Thermal expansion behavior of two-dimensional (2D) nitrides and graphene were studied by ab initio molecular dynamics (MD) simulations as well as quasiharmonic approximation (QHA). Anharmonicity of the acoustic phonon modes are related to the unusual negative thermal expansion (NTE) behavior of the nitrides. Our results also hint that direct ab initio MD simulations are a more elaborate method to investigate thermal expansion behavior of 2D materials than the QHA. Nevertheless, giant NTE coefficients are found for h-GaN and h-AlN within the covered temperature range 100-600 K regardless of the chosen computational method. This unusual NTE of 2D nitrides is reasoned with the out-of-plane oscillations related to the rippling behavior of the monolayers.

Abstract: The interplay between a magnetic field and the Coulomb potential from a charged vacancy on the electron states in graphene is investigated within the tight-binding model. The Coulomb potential removes locally Landau level degeneracy, while the vacancy introduces a satellite level next to the normal Landau level. These satellite levels are found throughout the positive-energy region, but in the negative-energy region, they turn into atomic collapse resonances. Crossings between Landau levels with different angular quantum number m are found. Unlike the point impurity system in which an anticrossing occurs between Landau levels of the same m, in this work anticrossing is found between the normal Landau level and the vacancy-induced level. The atomic collapse resonance hybridizes with the Landau levels. The charge at which the lowest Landau level m = -1, N = 1 crosses E = 0 increases with enhancing magnetic field. A Landau level scaling anomaly occurs when the charge is larger than the critical charge beta 0.6 and this critical charge is independent of the magnetic field.

Abstract: The McMillan-Ginzburg-Landau (MGL) model for charge density waves (CDW) is employed in a systematic phenomenological study of the different phases that have been probed in recent experiments involving transition metal dichalcogenides. We implemented an efficient imaginary time evolution method to solve the MGL equations, which enabled us to investigate the role of different coupling parameters on the CDW patterns and to perform calculations with different energy functionals that lead to several experimentally observed singularities in the CDW phase profiles. In particular, by choosing the appropriate energy functionals, we were able to obtain phases that go beyond the well-known periodic phase slips (discommensurations), exhibiting also topological defects (i.e., vortex-antivortex pairs), domain walls where the CDW order parameter is suppressed, and even CDW with broken rotational symmetry. Finally, we briefly discuss the effect of these different CDW phases on the profile and critical temperature of the competing superconducting state.

Abstract: By using the full 3D generalized time-dependent Ginzbug-Landau equation, we study a long superconducting film of finite width and thickness under an applied transport current. We show that, for sufficiently large thickness, the vortices and the antivortices become curved before they annihilate each other. As they approach the center of the sample, their ends combine, producing a single closed vortex. We also determine the critical values of the thickness for which the closed vortex sets in for different values of the Ginzburg-Ladau parameter. Finally, we propose a model of how to detect a closed vortex experimentally.

Abstract: Flat bands in the energy spectrum have attracted a lot of attention in recent years because of their unique properties and promising applications. Special arrangement of impurities on monolayer graphene are proposed to generate multiflat bands in the electronic band structure. In addition to the single midgap states in the spectrum of graphene with low hydrogen density, we found closely spaced bands around the Fermi level with increasing impurity density, which are similar to discrete lines in the spectrum of quantum dots, as well as the unusual Landau-level energy spectrum of graphene in the presence of a strong magnetic field. The presence of flat bands crucially depends on whether or not there are odd or even electrons of H(F) atoms bound to graphene. Interestingly, we found that a fully hydrogenated (fluoridated) of a hexagon of graphene sheet with six hydrogen (fluorine) atoms sitting on top and bottom in consecutive order exhibits Dirac cones in the electronic band structure with a 20% smaller Fermi velocity as compared to the pristine graphene. Functionalizing graphene introduces various C-C bond lengths resulting in nonuniform strains. Such a nonuniform strain may induce a giant pseudomagnetic field in the system, resulting in quantum Hall effect.

Abstract: A wave-packet time evolution method, based on the split-operator technique, is developed to investigate the scattering of quasiparticles at a normal-superconductor interface of arbitrary profile and shape. As a practical application, we consider a system where low-energy electrons can be described as Dirac particles, which is the case for most two-dimensional materials, such as graphene and transition-metal dichalcogenides. However, the method is easily adapted for other cases such as electrons in few-layer black phosphorus or any Schrodinger quasiparticles within the effective mass approximation in semiconductors. We employ the method to revisit Andreev reflection in mono-, bi-, and trilayer graphene, where specular-and retro-reflection cases are observed for electrons scattered by a steplike superconducting region. The effect of opening a zero-gap channel across the superconducting region on the electron and hole scattering is also addressed, as an example of the versatility of the technique proposed here.

Abstract: The superconducting properties of two recently synthesized ternary silicides with unit formula X₄CoSi (X = Nb, Ta) are investigated through ab initio calculations combined with Eliashberg theory. Interestingly, their crystal structure comprises interlocking honeycomb networks of Nb/Ta atoms. Nb₄CoSi is found to harbor better conditions for phonon-mediated superconductivity, as it possesses a higher density of states at the Fermi level, fostering stronger electron-phonon coupling. The superconducting critical temperatures (T_c) follow the same trend, with Nb₄CoSi having a twice higher value than Ta₄CoSi. Furthermore, the calculated T_c values (5.9 K vs 3.1 K) agree excellently with the experimentally obtained ones, establishing superconductivity in this new materials class as mediated by the electron-phonon coupling. Furthermore, my calculations show that the superconducting properties of these compounds do not simply correlate with the parameters of their honeycomb networks, contrary to proposals raised in the literature. Rather, their complete fermiology and phonon spectrum should be taken into account in order to explain their respective superconducting properties.

Abstract: We theoretically investigate the binding energy and electron-hole (e-h) overlap of excitonic states confined at the interface between two-dimensional materials with type-II band alignment, i.e., with lowest conduction and highest valence band edges placed in different materials, arranged in a side-by-side planar heterostructure. We propose a variational procedure within the effective mass approximation to calculate the exciton ground state and apply our model to a monolayer MoS2/WS2 heterostructure. The role of nonabrupt interfaces between the materials is accounted for in our model by assuming a WxMo1-xS2 alloy around the interfacial region. Our results demonstrate that (i) interface-bound excitons are energetically favorable only for small interface thickness and/or for systems under high dielectric screening by the materials surrounding the monolayer, and that (ii) the interface exciton binding energy and its e-h overlap are controllable by the interface width and dielectric environment.