Abstract: Different allotropes of phosphorene are possible of which black and blue phosphorus are the most stable. While blue phosphorus has isotropic properties, black phosphorus is strongly anisotropic in its electronic and optical properties due to its anisotropic crystal structure. In this work, we systematically investigated the lattice thermal properties of black and blue phosphorene by using first-principles calculations based on the quasiharmonic approximation approach. Similar to the optoelectronic and electronic properties, we predict that black phosphorene has highly anisotropic thermal properties, in contrast to the blue phase. The linear thermal expansion coefficients along the zigzag and armchair direction differ up to 20% in black phosphorene. The armchair direction of black phosphorene is more expandable as compared to the zigzag direction and the biaxial expansion of blue phosphorene under finite temperature. Our comparative analysis reveals that the inclusion of finite-temperature effects makes the blue phase thermodynamically more stable over the black phase above 135 K.

Abstract: Water confined between two graphene layers with a separation of a few A forms a layered two-dimensional ice structure. Using large scale molecular dynamics simulations with the adoptable ReaxFF interatomic potential we found that flat monolayer ice with a rhombic-square structure nucleates between the graphene layers which is nonpolar and nonferroelectric. We provide different energetic considerations and H-bonding results that explain the interlayer and intralayer properties of two-dimensional ice. The controversial AA stacking found experimentally [Algara-Siller et al., Nature (London) 519, 443 (2015)] is consistent with our minimum-energy crystal structure of bilayer ice. Furthermore, we predict that an odd number of layers of ice has the same lattice structure as monolayer ice, while an even number of ice layers exhibits the square ice AA stacking of bilayer ice.

Abstract: Employing density functional theory-based methods, we investigate monolayer and bilayer structures of hexagonal SnS2, which is a recently synthesized monolayer metal dichalcogenide. Comparison of the 1H and 1T phases of monolayer SnS2 confirms the ground state to be the 1T phase. In its bilayer structure we examine different stacking configurations of the two layers. It is found that the interlayer coupling in bilayer SnS2 is weaker than that of typical transition-metal dichalcogenides so that alternative stacking orders have similar structural parameters and they are separated with low energy barriers. A possible signature of the stacking order in the SnS2 bilayer has been sought in the calculated absorbance and reflectivity spectra. We also study the effects of the external electric field, charging, and loading pressure on the characteristic properties of bilayer SnS2. It is found that (i) the electric field increases the coupling between the layers at its preferred stacking order, so the barrier height increases, (ii) the bang gap value can be tuned by the external E field and under sufficient E field, the bilayer SnS2 can become a semimetal, (iii) the most favorable stacking order can be switched by charging, and (iv) a loading pressure exceeding 3 GPa changes the stacking order. The E-field tunable band gap and easily tunable stacking sequence of SnS2 layers make this 2D crystal structure a good candidate for field effect transistor and nanoscale lubricant applications.

Abstract: In a recent paper [Phys. Rev. B 89, 125418 (2014)], the authors argue that it is possible to map the electronic properties of twisted bilayer graphene to those of bilayer graphene in an in-plane magnetic field. However, their description of the low-energy dynamics of twisted bilayer graphene is restricted to the extended zone scheme and therefore neglects the effects of the superperiodic structure. If the energy spectrum is studied in the supercell Brillouin zone, we find that the comparison with an in-plane magnetic field fails because (i) the energy spectra of the two situations exhibit different symmetries and (ii) the low-energy spectra are very different.

Abstract: Simulations are carried out based on the dynamical mean-field theory (DMFT) in order to investigate the properties of correlated thin films for various values of the chemical potential, temperature, interaction strength, and applied transverse electric field. Application of a sufficiently strong field to a thin film at half filling leads to the appearance of conducting regions near the surfaces of the film, whereas in doped slabs the application of a field leads to a conductivity enhancement on one side of the film and a gradual transition to the insulating state on the opposite side. In addition to the inhomogeneous DMFT, a local density approximation (LDA) is considered in which the particle density n, quasiparticle residue Z, and spectral weight at the Fermi level A(ω=0) of each layer are approximated by a homogeneous bulk environment. A systematic comparison between the two approaches reveals that the less expensive LDA results are in good agreement with the DMFT approach, except close to the metal-to-insulator transition points and in the layers immediately at the film surfaces. LDA values for n are overall more reliable than those for Z and A(ω=0). The hysteretic behavior (memory effect) characteristic of the bulk doping driven Mott transition persists in the slab.

Abstract: Within a tight binding approach we investigate the energy levels of hexagonal and triangular bilayer graphene (BLG) quantum dots (QDs) with zigzag and armchair edges. We study AA- and AB-(Bernal) stacked BLG QDs and obtain the energy levels in both the absence and the presence of a perpendicular electric field (i.e., biased BLG QDs). Our results show that the size dependence of the energy levels is different from that of monolayer graphene QDs. The energy spectrum of AB-stacked BLG QDs with zigzag edges exhibits edge states which spread out into the opened energy gap in the presence of a perpendicular electric field. We found that the behavior of these edges states is different for the hexagonal and triangular geometries. In the case of AA-stacked BLG QDs, the electron and hole energy levels cross each other in both cases of armchair and zigzag edges as the dot size or the applied bias increases.

Abstract: Using density functional theory, we study the electronic properties of a graphene-silicene bilayer (GSB). A single layer of silicene binds to the graphene layer with adhesion energy of about 25 meV/atom. This adhesion energy between the two layers follows accurately the well-known -1/z(2) dispersion energy as found between two infinite parallel plates. In small flakes of GSB with hydrogenated edges, negative charge is transferred from the graphene layer to the silicene layer, producing a permanent and a switchable polar bilayer, while in an infinite GSB, the negative charge is transferred from the silicene layer to the graphene layer. The graphene-silicene bilayer is a good candidate for a nanocapacitor with piezoelectric capabilities. We found that the permanent dipole of the bilayer can be tuned by an external perpendicular electric field.

Abstract: Using the tight-binding approach, we investigate confined states in two different hybrid monolayer-bilayer systems: (i) a hexagonal monolayer area surrounded by bilayer graphene in the presence of a perpendicularly applied electric field and (ii) a hexagonal bilayer graphene dot surrounded by monolayer graphene. The dependence of the energy levels on dot size and external magnetic field is calculated. We find that the energy spectrum for quantum dots with zigzag edges consists of states inside the gap which range from dot-localized states, edge states, to mixed states coexisting together, whereas for dots with armchair edges, only dot-localized states are observed.

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: Finding novel atomically thin heterostructures and understanding their characteristic properties are critical for developing better nanoscale optoelectronic devices. In this study, we investigate the electronic and optical properties of a GaS-Ca(OH)(2) heterostructure using first-principle calculations. The band gap of the GaS-Ca(OH)(2) heterostructure is significantly reduced when compared to those of the isolated constituent layers. Our calculations showthat the GaS-Ca(OH)(2) heterostructure is a type-II heterojunction which can be used to separate photoinduced charge carriers where electrons are localized in GaS and holes in the Ca(OH)(2) layer. This leads to spatially indirect excitons which are important for solar energy and optoelectronic applications due to their long lifetime. By solving the Bethe-Salpeter equation on top of a single shot GW calculation (G(0)W(0)), the dielectric function and optical oscillator strength of the constituent monolayers and the heterostructure are obtained. The oscillator strength of the optical transition for the GaS monolayer is an order of magnitude larger than the Ca(OH)(2) monolayer. We also found that the calculated optical spectra of different stacking types of the heterostructure show dissimilarities, although their electronic structures are rather similar. This prediction can be used to determine the stacking type of ultrathin heterostructures.

Abstract: In order to study quantum effects in a two-dimensional crystal lattice of a free-standing monolayer graphene, we have performed both path-integral Monte Carlo (PIMC) and classical Monte Carlo (MC) simulations for temperatures up to 2000 K. The REBO potential is used for the interatomic interaction. The total energy, interatomic distance, root-mean-square displacement of the atom vibrations, and the free energy of the graphene layer are calculated. The obtained lattice vibrational energy per atom from the classical MC simulation is very close to the energy of a three-dimensional harmonic oscillator 3k(B)T. The PIMC simulation shows that quantum effects due to zero-point vibrations are significant for temperatures T < 1000 K. The quantum contribution to the lattice vibrational energy becomes larger than that of the classical lattice for T < 400 K. The lattice expansion due to the zero-point motion causes an increase of 0.53% in the lattice parameter. A minimum in the lattice parameter appears at T similar or equal to 500 K. Quantum effects on the atomic vibration amplitude of the graphene lattice and its free energy are investigated.

Abstract: The dynamics of driven pancake vortices in layered superconductors is studied using molecular-dynamics simulations. We found that, with increasing driving force, for strong interlayer coupling, the preexisted vortex lines either directly depin or first transform to two-dimensional (2D) pinned states before they are depinned, depending on the pinning strength. In a narrow region of pinning strengths, we found an interesting repinning process, which results in a negative differential resistance. For weak interlayer coupling, individually pinned pancake vortices first form disordered 2D flow and then transform to ordered three-dimensional (3D) flow with increasing driving force. However, for extremely strong pinning, the random pinning-induced thermal-like Langevin forces melt 3D vortex lines, which results in a persistent 2D flow in the fast-sliding regime. In the intermediate regime, the peak effect is found: With increasing driving force, the moving pancake vortices first crystallize to moving 3D vortex lines, and then these 3D vortex lines are melted, leading to the appearance of a reentrant 2D flow state. Our results are summarized in a dynamical phase diagram.

Abstract: We investigate silicene superlattices in the presence of a tunable barrier potential U, an exchange field M, and a perpendicular electric field E-z. The resulting miniband structure depends on the spin and valley indices and on the fields M and E-z. These fields determine the minigaps and also affect the additional Dirac points brought about by the periodic potential U. In addition, we consider diffusive transport and assess its dependence on the spin and valley indices as well as on temperature. The corresponding spin and valley polarizations strongly depend on the potential U and can be made almost 100% at very low temperatures at particular values of the Fermi energy.

Abstract: The electronic properties, carrier mobility, and strain response of TiS3 nanoribbons (TiS3 NRs) are investigated by first-principles calculations. We found that the electronic properties of TiS3 NRs strongly depend on the edge type (a or b). All a-TiS3 NRs are metallic with a magnetic ground state, while b-TiS3 NRs are direct band gap semiconductors. Interestingly, the size of the band gap and the band edge position are almost independent of the ribbon width. This feature promises a constant band gap in a b-TiS3 NR with rough edges, where the ribbon width differs in different regions. The maximum carrier mobility of b-TiS3 NRs is calculated by using the deformation potential theory combined with the effective mass approximation and is found to be of the order 10(3) cm(2) V-1 s(-1). The hole mobility of the b-TiS3 NRs is one order of magnitude lower, but it is enhanced compared to the monolayer case due to the reduction in hole effective mass. The band gap and the band edge position of b-TiS3 NRs are quite sensitive to applied strain. In addition we investigate the termination of ribbon edges by hydrogen atoms. Upon edge passivation, the metallic and magnetic features of a-TiS3 NRs remain unchanged, while the band gap of b-TiS3 NRs is increased significantly. The robust metallic and ferromagnetic nature of a-TiS3 NRs is an essential feature for spintronic device applications. The direct, width-independent, and strain-tunable band gap, as well as the high carrier mobility, of b-TiS3 NRs is of potential importance in many fields of nanoelectronics, such as field-effect devices, optoelectronic applications, and strain sensors.

Abstract: Anharmonic effects in an atomic monolayer thin crystal with honeycomb lattice structure are investigated by analytical and numerical lattice dynamical methods. Starting from a semiempirical model for anharmonic couplings of third and fourth orders, we study the in-plane and out-of-plane (flexural) mode components of the generalized wave vector dependent Gruneisen parameters, the thermal tension and the thermal expansion coefficients as a function of temperature and crystal size. From the resonances of the displacement-displacement correlation functions, we obtain the renormalization and decay rate of in-plane and flexural phonons as a function of temperature, wave vector, and crystal size in the classical and in the quantum regime. Quantitative results are presented for graphene. There, we find that the transition temperature T-alpha from negative to positive thermal expansion is lowered with smaller system size. Renormalization of the flexural mode has the opposite effect and leads to values of T-alpha approximate to 300 K for systems of macroscopic size. Extensive numerical analysis throughout the Brillouin zone explores various decay and scattering channels. The relative importance of normal and umklapp processes is investigated. The work is complementary to crystalline membrane theory and computational studies of anharmonic effects in two-dimensional crystals.