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: The strain-induced pseudomagnetic field in supported graphene deposited on top of a nanostructured substrate is investigated by using atomistic simulations. A step, an elongated trench, a one-dimensional barrier, a spherical bubble, a Gaussian bump, and a Gaussian depression are considered as support structures for graphene. From the obtained optimum configurations we found very strong induced pseudomagnetic fields which can reach up to similar to 1000 T due to the strain-induced deformations in the supported graphene. Different magnetic confinements with controllable geometries are found by tuning the pattern of the substrate. The resulting induced magnetic fields for graphene on top of a step, barrier, and trench are calculated. In contrast to the step and trench the middle part of graphene on top of a barrier has zero pseudomagnetic field. This study provides a theoretical background for designing magnetic structures in graphene by nanostructuring substrates. We found that altering the radial symmetry of the deformation changes the sixfold symmetry of the induced pseudomagnetic field.

Abstract: The Stark effect is one of the most efficient mechanisms to manipulate many-body states in nanostructured systems. In mono- and few-layer transition metal dichalcogenides, it has been successfully induced by optical and electric field means. Here, we tune the optical emission energies and dissociate excitonic states in MoSe2 monolayers employing the 220 MHz in-plane piezoelectric field carried by surface acoustic waves. We transfer the monolayers to high dielectric constant piezoelectric substrates, where the neutral exciton binding energy is reduced, allowing us to efficiently quench (above 90%) and red-shift the excitonic optical emissions. A model for the acoustically induced Stark effect yields neutral exciton and trion in-plane polarizabilities of 530 and 630 x 10(-5) meV/(kV/cm)(2), respectively, which are considerably larger than those reported for monolayers encapsulated in hexagonal boron nitride. Large in-plane polarizabilities are an attractive ingredient to manipulate and modulate multiexciton interactions in two-dimensional semiconductor nanostructures for optoelectronic applications.

Abstract: The split-operator technique for wave packet propagation in quantum systems is expanded here to the case of propagatingwave functions describing Schrodinger particles, namely, charge carriers in semiconductor nanostructures within the effective mass approximation, in the presence of Zeeman effect, as well as of Rashba and Dresselhaus spin-orbit interactions. We also demonstrate that simple modifications to the expanded technique allow us to calculate the time evolution of wave packets describing Dirac particles, which are relevant for the study of transport properties in graphene.

Abstract: The spin interaction and the effective g factor of a magnetic exciton (ME) are investigated theoretically in a diluted magnetic semiconductor (DMS) quantum dot (QD), including the Coulomb interaction and the sp-d exchange interaction. At low magnetic field, the ME energy decreases rapidly with increasing magnetic field and saturates at high magnetic field for high Mn concentration. The ground state of the ME exhibits an interesting crossing behavior between sigma(+)-ME and sigma(-)-ME for low Mn concentration. The g(ex) factor of the ME in a DMS QD displays a monotonic decrease with increasing magnetic field and can be tuned to zero by an external magnetic field. (C) 2003 American Institute of Physics.

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 spatial dependence of the van der Waals (vdW) energy between graphene and hexagonal boron-nitride (h-BN) is investigated using atomistic simulations. The van der Waals energy between graphene and h-BN shows a hexagonal superlattice structure identical to the observed Moire pattern in the local density of states, which depends on the lattice mismatch and misorientation angle between graphene and h-BN. Our results provide atomistic features of the weak van der Waals interaction between graphene and BN which are in agreement with experiment and provide an analytical expression for the size of the spatial variation of the weak van der Waals interaction. We also found that the A-B-lattice symmetry of graphene is broken along the armchair direction. (C) 2014 AIP Publishing LLC.

Abstract: The short exciton lifetime and strong exciton-exciton interaction in transition metal dichalcogenides limit the efficiency of exciton emission. Here, the authors show that exciton-exciton interaction in monolayer WS2 can be screened using proximal metal plates, leading to an improved quantum yield. Monolayer transition metal dichalcogenides (1L-TMDs) have tremendous potential as atomically thin, direct bandgap semiconductors that can be used as convenient building blocks for quantum photonic devices. However, the short exciton lifetime due to the defect traps and the strong exciton-exciton interaction in TMDs has significantly limited the efficiency of exciton emission from this class of materials. Here, we show that exciton-exciton interaction in 1L-WS2 can be effectively screened using an ultra-flat Au film substrate separated by multilayers of hexagonal boron nitride. Under this geometry, induced dipolar exciton-exciton interaction becomes quadrupole-quadrupole interaction because of effective image dipoles formed within the metal. The suppressed exciton-exciton interaction leads to a significantly improved quantum yield by an order of magnitude, which is also accompanied by a reduction in the exciton-exciton annihilation (EEA) rate, as confirmed by time-resolved optical measurements. A theoretical model accounting for the screening of the dipole-dipole interaction is in a good agreement with the dependence of EEA on exciton densities. Our results suggest that fundamental EEA processes in the TMD can be engineered through proximal metallic screening, which represents a practical approach towards high-efficiency 2D light emitters.

Abstract: The self rotation of a graphene flake over graphite is controlled by the size, initial misalignment and temperature. Using both ab initio calculations and molecular dynamics simulations, we investigate annealing effects on the self rotation of a graphene flake on a graphene substrate. The energy barriers for rotation and drift of a graphene flake over graphene is found to be smaller than 25 meV/atom which is comparable to thermal energy. We found that small flakes (of about similar to 4 nm) are more sensitive to temperature and initial misorientation angles than larger one (beyond 10 nm). The initial stacking configuration of the flake is found to be important for its dynamics and time evolution of misalignment. Large flakes, which are initially in the AA-or AB-stacking state with small misorientation angle, rotate and end up in the AB-stacking configuration. However small flakes can they stay in an incommensurate state specially when the initial misorientation angle is larger than 2 degrees. Our results are in agreement with recent experiments.

Abstract: The scattering of two-dimensional (2D) massless Dirac electrons is investigated in the presence of a random array of circular mass barriers. The inverse momentum relaxation time and the Hall factor are calculated and used to obtain parallel and perpendicular resistivity components within linear transport theory. We found a nonzero perpendicular resistivity component which has opposite sign for electrons in the different K and K′ valleys. This property can be used for valley filter purposes. The total cross section for scattering on penetrable barriers exhibits resonances due to the presence of quasibound states in the barriers that show up as sharp gaps in the cross section while for Schrödinger electrons they appear as peaks.

Abstract: The role of dielectric screening of electron-hole interaction in van der Waals heterostructures is theoretically investigated. A comparison between models available in the literature for describing these interactions is made and the limitations of these approaches are discussed. A simple numerical solution of Poisson's equation for a stack of dielectric slabs based on a transfer matrix method is developed, enabling the calculation of the electron-hole interaction potential at very low computational cost and with reasonable accuracy. Using different potential models, direct and indirect exciton binding energies in these systems are calculated within Wannier-Mott theory, and a comparison of theoretical results with recent experiments on excitons in two-dimensional materials is discussed.

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 reduction of dimensionality is a very effective way to achieve appealing properties in two-dimensional materials (2DMs). First-principles calculations can greatly facilitate the prediction of 2DM properties and find possible approaches to enhance their performance. We employed first-principles calculations to gain insight into the impact of different types of point defects (vacancies and substitutional dopants) on the electronic and magnetic properties of a ZnSnN2 (ZSN) monolayer. We show that Zn, Sn, and N + Zn vacancy-defected structures are p-type conducting, while the defected ZSN with a N vacancy is n-type conducting. For substitutional dopants, we found that all doped structures are thermally and energetically stable. The most stable structure is found to be B-doping at the Zn site. The highest work function value (5.0 eV) has been obtained for Be substitution at the Sn site. Li-doping (at the Zn site) and Be-doping (at the Sn site) are p-type conducting, while B-doping (at the Zn site) is n-type conducting. We found that the considered ZSN monolayer-based structures with point defects are magnetic, except those with the N vacancy defects and Be-doped structures. The ab initio molecular dynamics simulations confirm that all substitutionally doped and defected structures are thermally stable. Thus, our results highlight the possibility of tuning the magnetism in ZnSnN2 monolayers through defect engineering.

Abstract: The recently discovered structural phase transition in Polymerized KC60 at about 50K leads to a doubling of the unit cell volume and is accompanied by a metal-insulator transition. Here, we show that the ((a) over right arrow + (c) over right arrow, (b) over right arrow, (a) over right arrow – (c) over right arrow) superstructure results from orientational charge density waves along the polymer chains and correlated displacements of the K+ ions. The presented model can also account for the metal-insulator transition. The effect is specific for the space group Pmnn of KC60 and is absent in both Rb- and CsC60 (space group 12/m), in agreement with the present experimental knowledge of these compounds.

Abstract: The recent discovery of axion states in materials such as antiferromagnetic topological insulators has boosted investigations of the magnetoelectric response in topological insulators and their promise towards realizing dissipationless topological electronics. In this paper, we develop a tight-binding methodology to explore the emergence of axion states in Bi2Se3 in proximity to magnetic insulators on the top and bottom surfaces. The topological protection of the surface states is lifted by a time-reversal-breaking perturbation due to the proximity of a magnetic insulator, and a gap is opened on the surfaces, giving rise to half-quantized Hall conductance and a zero Hall plateau-evidencing an axion insulator state. We developed a real-space tight-binding Hamiltonian for Bi2Se3 using first-principles data. Transport properties of the system were obtained within the Landauer-Buttiker formalism, and we discuss the creation of axion states through Hall conductance and a zero Hall plateau at the surfaces, as a function of proximitized magnetization and corresponding potentials at the surfaces, as well as the thickness of the topological insulator.