Abstract: The energy spectrum and eigenstates of single-layer black phosphorus nanoribbons in the presence of a perpendicular magnetic field and an in-plane transverse electric field are investigated by means of a tight-binding method, and the effect of different types of edges is examined analytically. A description based on a continuum model is proposed using an expansion of the tight-binding model in the long-wavelength limit. Thewave functions corresponding to the flatband part of the spectrum are obtained analytically and are shown to agree well with the numerical results from the tight-binding method for both narrow (10 nm) and wide (100 nm) nanoribbons. Analytical expressions for the critical magnetic field at which Landau levels are formed and the ranges of wave numbers in the dispersionless flatband segments in the energy spectra are derived. We examine the evolution of the Landau levels when an in-plane lateral electric field is applied, and we determine analytically how the edge states shift withmagnetic field. For wider nanoribbons, the conductance is shown to have a characteristic staircase shape in combined magnetic and electric fields. Some of the stairs in zigzag and skewed armchair nanoribbons originate from edge states that are found in the band gap.

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: The influence of a perpendicular magnetic field on the binding energy and structural properties of excitons, trions, and biexcitons in monolayers of semiconducting transition metal dichalcogenides (TMDs) is investigated. The stochastic variational method (SVM) with a correlated Gaussian basis is used to calculate the different properties of these few-particle systems. In addition, we present a simplified variational approach which supports the SVM results for excitons as a function of magnetic field. The exciton diamagnetic shift is compared with recent experimental results, and we extend this concept to trions and biexcitons. The effect of a local potential fluctuation, which we model by a circular potential well, on the binding energy of trions and biexcitons is investigated and found to significantly increase the binding of those excitonic complexes.

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 magnetic-field dependence of the energy spectrum, wave function, binding energy, and oscillator strength of exciton states confined in a circular graphene quantum dot (CGQD) is obtained within the configuration interaction method. We predict that (i) excitonic effects are very significant in the CGQD as a consequence of a combination of geometric confinement, magnetic confinement, and reduced screening; (ii) two types of excitons (intravalley and intervalley) are present in the CGQD because of the valley degree of freedom in graphene; (iii) the intravalley and intervalley exciton states display different magnetic-field dependencies due to the different electron-hole symmetries of the single-particle energy spectra; (iv) with increasing magnetic field, the exciton ground state in the CGQD undergoes an intravalley to intervalley transition accompanied by a change of angular momentum; (v) the exciton binding energy does not increase monotonically with the magnetic field due to the competition between geometric and magnetic confinements; and (vi) the optical transitions of the intervalley and intravalley excitons can be tuned by the magnetic field, and valley-dependent excitonic transitions can be realized in a CGQD.

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: The modification of the quantum states in a Dice lattice due to a Coulomb impurity are investigated. The energy-band structure of a pristine Dice lattice consists of a Dirac cone and a flat band at the Dirac point. We use the tight-binding formalism and find that the flat band states transform into a set of discrete bound states whose electron density is localized on a ring around the impurity mainly on two of the three sublattices. Its energy is proportional to the strength of the Coulomb impurity. Beyond a critical strength of the Coulomb potential atomic collapse states appear that have some similarity with those found in graphene with the difference that the flat band states contribute with an additional ringlike electron density that is spatially decoupled from the atomic collapse part. At large value of the strength of the Coulomb impurity the flat band bound states anticross with the atomic collapse states.

Abstract: The possibility of an inhomogeneous charge density wave phase is investigated in a system of two coupled electron and hole monolayers separated by a hexagonal boron nitride insulating layer. The charge-density-wave state is induced through the assumption of negative compressibility of electron/hole gases in a Coulomb drag configuration between the electron and hole sheets. Under equilibrium conditions, we derive analytical expressions for the density oscillation along the zigzag and armchair directions. We find that the density modulation not only depends on the sign of the compressibility but also on the anisotropy of the low-energy bands. Our results are applicable to any two-dimensional system with anisotropic parabolic bands, characterized by different effective masses. For equal effective masses, i.e., isotropic energy bands, our results agree with Hroblak et al. [Phys. Rev. B 96, 075422 (2017)]. Our numerical results are applied to phosphorene.

Abstract: The recent development in the field of two-dimensional magnetic materials urges reliable theoretical methodology for determination of magnetic properties. Among the available methods, ab initio four-state energy mapping based on density functional theory stands out as a powerful technique to calculate the magnetic exchange interaction in the Heisenberg spin model. Although the required formulas were explained in earlier works, the considered Hamiltonian in those studies always corresponded to the specific case that the off-diagonal part of J matrix is antisymmetric, which may be misleading in other cases. Therefore, using the most general form of the Heisenberg spin Hamiltonian, we here derive the generic formulas. With a proper choice of four different magnetic states, a single formula governs all elements of the exchange interaction matrix for any considered pair of spin sites.

Abstract: The report on square ice sandwiched between two graphene layers by Algara-Siller et al. [Nature (London) 519, 443 (2015)] has generated a large interest in this system. By applying high lateral pressure on nanoconfined water, we found that monolayer ice is transformed to bilayer ice when the two graphene layers are separated by H = 6,7 angstrom. It was also found that three layers of a denser phase of ice with smaller lattice constant are formed if we start from bilayer ice and apply a lateral pressure of about 0.7 GPa with H = 8,9 angstrom. The lattice constant (2.5-2.6 angstrom) in both transitions is found to be smaller than those typical for the known phases of ice and water, i.e., 2.8 angstrom. We validate these results using ab initio calculations and find good agreement between ab initio O-O distance and those obtained from classical molecular dynamics simulations. The reversibility of the mentioned transitions is confirmed by decompressing the systems.

Abstract: The spectrum of a Coulomb impurity in bilayer graphene is investigated as function of the strength of a perpendicular magnetic field for different values of the angular quantum number m and for different values of the gate voltage. We point out fundamental differences between the results from the two-band and four-band model. The supercritical instability and fall-to-center phenomena are investigated in the presence of a magnetic field. We find that in the four-band model the fall-to-center phenomenon occurs as in monolayer graphene, while this is not the case in the two-band model. We find that in a magnetic field the supercritical instability manifests itself as a series of anticrossings in the hole part of the spectrum for states coming from the low-energy band. However, we also find very distinct anticrossings in the electron part of the spectrum that continue into the hole part, which are related to the higher energy band of the four-band model. At these anticrossings, we find a very sharp peak in the probability density close to the impurity, reminiscent for the fall-to-center phenomenon. In this paper, these peculiar and interesting effects are studied for different magnetic field, interlayer coupling, and bias potential strengths.

Abstract: The structural and mechanical properties of N-doped graphene (NG) are investigated using reactive force field (ReaxFF) potentials in large-scale molecular dynamics simulations. We found that ripples, which are induced by the dopants, change the roughness of NG, which depends on the number of dopants and their local arrangement. For any doping ratio N/C, the NG becomes ferroelectric with a net dipole moment. The formation energy increases nonlinearly with N/C ratio, while the Young's modulus, tensile strength, and intrinsic strain decrease with the number of dopants. Our results for the structural deformation and the thermoelectricity of the NG sheet are in good agreement with recent experiments and ab initio calculations.

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: 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 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.