Abstract: Starting from the single-particle Dirac Hamiltonian for charge carriers in monolayer transition metal dichalcogenides (TMDs), we construct a four-band Hamiltonian describing interlayer excitons consisting of an electron in one TMD layer and a hole in the other TMD layer. An expression for the electron-hole interaction potential is derived, taking into account the effect of the dielectric environment above, below, and between the two TMD layers as well as polarization effects in the transition metal layer and in the chalcogen layers of the TMD layers. We calculate the interlayer exciton binding energy and average in-plane interparticle distance for different TMD heterostructures. The effect of different dielectric environments on the exciton binding energy is investigated and a remarkable dependence on the dielectric constant of the barrier between the two layers is found, resulting from competing effects as a function of the in-plane and out-of-plane dielectric constants of the barrier. The polarization effects in the chalcogen layers, which in general reduce the exciton binding energy, can lead to an increase in binding energy in the presence of strong substrate effects by screening the substrate. The excitonic absorbance spectrum is calculated and we show that the interlayer exciton peak depends linearly on a perpendicular electric field, which agrees with recent experimental results.

Abstract: We explore the electronic and the optical properties of monolayer TiS3, which shows in-plane anisotropy and is composed of a chain-like structure along one of the lattice directions. Together with its robust direct band gap, which changes very slightly with stacking order and with the thickness of the sample, the anisotropic physical properties of TiS3 make the material very attractive for various device applications. In this study, we present a detailed investigation on the effect of the crystal anisotropy on the excitons and the trions of the TiS3 monolayer. We use many-body perturbation theory to calculate the absorption spectrum of anisotropic TiS3 monolayer by solving the Bethe-Salpeter equation. In parallel, we implement and use a Wannier-Mott model for the excitons that takes into account the anisotropic effective masses and Coulomb screening, which are obtained from ab initio calculations. This model is then extended for the investigation of trion states of monolayer TiS3. Our calculations indicate that the absorption spectrum of monolayer TiS3 drastically depends on the polarization of the incoming light, which excites different excitons with distinct binding energies. In addition, the binding energies of positively and the negatively charged trions are observed to be distinct and they exhibit an anisotropic probability density distribution.

Abstract: Edge functionalization of graphene nanoribbons with nitrogen atoms for various adatom configurations at armchair and zigzag edges are investigated. We provide comprehensive information on the electronic and magnetic properties and investigate the stability of the various systems. Two types of rippling of the nanoribbons, namely edge and bulk rippling depending on the sign of edge stress induced at the edge, are found. They are found to play the decisive role for the stability of the structures. We also propose a type of edge decoration in which every third nitrogen adatom at the zigzag edges is replaced by an oxygen atom. In this way, the electron count is compatible with a full aromatic structure, leading to additional stability and a disappearance of magnetism that is usually associated with zigzag nanoribbons.

Abstract: Massless Dirac fermions (MDFs) emerge as quasiparticles in various novel materials such as graphene and topological insulators, and they exhibit several intriguing properties, of which Veselago focusing is an outstanding example with a lot of possible applications. However, up to now Veselago focusing merely occurred in p-n junction devices based on the isotropic MDF, which lacks the tunability needed for realistic applications. Here, motivated by the emergence of novel Dirac materials, we investigate the propagation behaviors of anisotropic MDFs in such a p-n junction structure. By projecting the Hamiltonian of the anisotropic MDF to that of the isotropic MDF and deriving an exact analytical expression for the propagator, precise Veselago focusing is demonstrated without the need for mirror symmetry of the electron source and its focusing image. We show a tunable focusing position that can be used in a device to probe masked atom-scale defects. This study provides an innovative concept to realize Veselago focusing relevant for potential applications, and it paves the way for the design of novel electron optics devices by exploiting the anisotropic MDF.

Abstract: Excitonic superfluidity in double phosphorene monolayers is investigated using the BCS mean-field equations. Highly anisotropic superfluidity is predicted where we found that the maximum superfluid gap is in the Bose-Einstein condensate (BEC) regime along the armchair direction and in the BCS-BEC crossover regime along the zigzag direction. We estimate the highest Kosterlitz-Thouless transition temperature with maximum value up to similar to 90 K with onset carrier densities as high as 4 x 10(12) cm(-2). This transition temperature is significantly larger than what is found in double electron-hole few-layers graphene. Our results can guide experimental research toward the realization of anisotropic condensate states in electron-hole phosphorene monolayers.

Abstract: Water inside a nanocapillary exhibits unconventional structural and dynamical behavior due to its ordered structure. The confining walls, density, and lateral pressures control profoundly the microscopic structure of trapped water. Here we study the electrostriction of confined water subjected to pressures of the order of GPa for two different setups: (i) a graphene nanochannel containing a constant number of water molecules independent of the height of the channel, (ii) an open nanochannel where water molecules can be exchanged with those in a reservoir. For the former case, a square-rhombic structure of confined water is formed when the height of the channel is d = 6.5 angstrom having a density of rho = 1.42 g cm(-3). By increasing the height of the channel, a transition from a flat to a buckled state occurs, whereas the density rapidly decreases and reaches the bulk density for d congruent to 8.5 angstrom. When a perpendicular electric field is applied, the water structure and the lateral pressure change. For strong electric fields (similar to 1 V/angstrom), the square-rhombic structure is destroyed. For an open setup, a solid phase of confined water consisting of an imperfect square-rhombic structure is formed. By applying a perpendicular field, the density and phase of confined water change. However, the density and pressure inside the channel decrease as compared to the first setup. Our study is closely related to recent experiments on confined water, and it reveals the sensitivity of the microscopic structure of confined water to the size of the channel, the external electric field, and the experimental setup.

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: By performing density functional theory-based ab initio calculations, Raman-active phonon modes of single-layer two-dimensional (2D) materials and the effect of in-plane biaxial strain on the peak frequencies and corresponding activities of the Raman-active modes are calculated. Our findings confirm the Raman spectrum of the unstrained 2D crystals and provide expected variations in the Raman-active modes of the crystals under in-plane biaxial strain. The results are summarized as follows: (i) frequencies of the phonon modes soften (harden) under applied tensile (compressive) strains; (ii) the response of the Raman activities to applied strain for the in-plane and out-of-plane vibrational modes have opposite trends, thus, the built-in strains in the materials can be monitored by tracking the relative activities of those modes; (iii) in particular, the A peak in single-layer Si and Ge disappears under a critical tensile strain; (iv) especially in mono-and diatomic single layers, the shift of the peak frequencies is a stronger indication of the strain rather than the change in Raman activities; (v) Raman-active modes of single-layer ReX2 (X = S, Se) are almost irresponsive to the applied strain. Strain-induced modifications in the Raman spectrum of 2D materials in terms of the peak positions and the relative Raman activities of the modes could be a convenient tool for characterization.

Abstract: Using the tight-binding approach, we investigate the energy spectrum of square, triangular, and hexagonal MoS2 quantum dots (QDs) in the presence of a perpendicular magnetic field. Novel edge states emerge in MoS2 QDs, which are distributed over the whole edge which we call ring states. The ring states are robust in the presence of spin-orbit coupling (SOC). The corresponding energy levels of the ring states oscillate as a function of the perpendicular magnetic field which are related to Aharonov-Bohm oscillations. Oscillations in the magnetic field dependence of the energy levels and the peaks in the magneto-optical spectrum emerge (disappear) as the ring states are formed (collapsed). The period and the amplitude of the oscillation decrease with the size of the MoS2 QDs.

Abstract: Defects are almost inevitably present in realistic materials and defective materials are expected to exhibit very different properties than their nondefective (perfect) counterparts. Here, using a combination of the tight-binding approach and the scattering matrix formalism, we investigate the electronic transport properties of defective phosphorene nanoribbons (PNRs) containing atomic vacancies. We find that for both armchair PNRs (APNRs) and zigzag PNRs (ZPNRs), single vacancies can create quasilocalized states, which can affect their conductance. With increasing vacancy concentration, three different transport regimes are identified: ballistic, diffusive, and Anderson localized ones. In particular, ZPNRs that are known to be metallic due to the presence of edge states become semiconducting: edge conductance vanishes and transport gap appears due to Anderson localization. Moreover, we find that for a fixed vacancy concentration, both APNRs and ZPNRs of narrower width and/or longer length are more sensitive to vacancy disorder than their wider and/or shorter counterparts, and that for the same ribbon length and width, ZPNRs are more sensitive to vacancy disorder than APNRs.

Abstract: <script type='text/javascript'>document.write(unpmarked('We show that delaminations in bilayer graphene (BLG) with electrostatically induced interlayer symmetry can provide one with ballistic channels for electrons with energies inside the electrostatically induced BLG gap. These channels are formed by a combination of valley-polarized evanescent states propagating along the delamination edges (which persist in the presence of a strong magnetic field) and standing waves bouncing between them inside the delaminated region (in a strong magnetic field, these transform into Landau levels in the monolayers). For inverted stackings in BLGs on the left and right of the delamination (AB-2ML-BA or BA-2ML-AB, where 2ML indicates two decoupled monolayers of graphene), the lowest-energy ballistic channels are gapless, have linear dispersion, and appear to be weakly topologically protected. When BLG stackings on both sides of the delamination are the same (AB-2ML-AB or BA-2ML-BA), the lowest-energy ballistic channels are gapped, with a gap epsilon(g) scaling as epsilon(g) alpha W-1 with delamination width and epsilon(g) alpha delta(-1) with the on-layer energy difference in the delaminated part of the structure. Depending on the width, delaminations may also support several \u0022higher-energy\u0022 waveguide modes. Our results are based on both the analytical study of the wave matching of Dirac states and tight-binding model calculations, and we analyze in detail the dependence of the delamination spectrum on the electrostatic conditions in the structure, such as the vertical displacement field.'));

Abstract: <script type='text/javascript'>document.write(unpmarked('We present an advanced method to study spin fluctuations in superconductors quantitatively and entirely from first principles. This method can be generally applied to materials where electron-phonon coupling and spin fluctuations coexist. We employ it here to examine the recently synthesized superconductor iron tetraboride (FeB4) with experimental T-c similar to 2.4 K [H. Gou et al., Phys. Rev. Lett, 111, 157002 (2013)]. We prove that FeB4 is particularly prone to ferromagnetic spin fluctuations due to the presence of iron, resulting in a large Stoner interaction strength, I = 1.5 eV, as calculated from first principles. The other important factor is its Fermi surface that consists of three separate sheets, among which two are nested ellipsoids. The resulting susceptibility has a ferromagnetic peak around q = 0, from which we calculated the repulsive interaction between Cooper pair electrons using the random phase approximation. Subsequently, we combined the electron-phonon interaction calculated from first principles with the spin fluctuation interaction in fully anisotropic Eliashberg theory calculations. We show that the resulting superconducting gap spectrum is conventional, yet very strongly depleted due to coupling to the spin fluctuations. The critical temperature decreases from T-c = 41 K, if they are not taken into account, to T-c = 1.7 K, in good agreement with the experimental value.'));

Abstract: Starting from the shell structure in atoms and the significant correlation within electron pairs, we distinguish the exchange-correlation effects between two electrons of opposite spins occupying the same orbital from the average correlation among many electrons in a crystal. In the periodic potential of the crystal with lattice constant larger than the effective Bohr radius of the valence electrons, these correlated electron pairs can form a metastable energy band above the corresponding single-electron band separated by an energy gap. In order to determine if these metastable electron pairs can be stabilized, we calculate the many-electron exchange-correlation renormalization and the polaron correction to the two-band system with single electrons and electron pairs. We find that the electron-phonon interaction is essential to counterbalance the Coulomb repulsion and to stabilize the electron pairs. The interplay of the electron-electron and electron-phonon interactions, manifested in the exchange-correlation energies, polaron effects, and screening, is responsible for the formation of electron pairs (bipolarons) that are located on the Fermi surface of the single-electron band.

Abstract: Environmental parameters such as temperature, relative humidity, visible light, UV radiation and pollution influence the deterioration rate of heritage items. To judge on the environmental appropriateness for heritage conservation, it is therefore important to monitor the environment. Often, an incomplete set of environmental parameters is measured, or sporadic or time-averaged measurements are performed. As a result, a wide range of undesirable situations and hazards remain unnoticed. This might lead to an underestimation of environmental dangers (i.e., inaccurate judgement) or to inappropriate mitigation measures (i.e., inaccurate decision making). We present an innovative and user-friendly monitoring device that simultaneously and continuously measures (1) environmental parameters and (2) material behavior. An extended combination of off-the-shelf sensors for temperature, relative humidity, air speed, CO2, NO2, O-3 and particulate matter are connected to a multipurpose datalogger. In-house developed sensors for the shrinkage and expansion behavior of wood, as well as sensors for metal corrosion rates are connected to the same datalogger. Such extended monitoring shows the identification of a wider range of undesirable situations, and it facilitates the search for correlations between such situations and the sources that cause them, i.e., the hazards.

Abstract: By means of replica exchange molecular dynamics simulations we investigate how the length of a silk-like, alternating diblock oligopeptide influences its secondary and quaternary structure. We carry out simulations for two protein sizes consisting of three and five blocks, and study the stability of a single protein, a dimer, a trimer and a tetramer. Initial configurations of our simulations are β-roll and β-sheet structures. We find that for the triblock the secondary and quaternary structures upto and including the tetramer are unstable: the proteins melt into random coil structures and the aggregates disassemble either completely or partially. We attribute this to the competition between conformational entropy of the proteins and the formation of hydrogen bonds and hydrophobic interactions between proteins. This is confirmed by our simulations on the pentablock proteins, where we find that, as the number of monomers in the aggregate increases, individual monomers form more hydrogen bonds whereas their solvent accessible surface area decreases. For the pentablock β-sheet protein, the monomer and the dimer melt as well, although for the β-roll protein only the monomer melts. For both trimers and tetramers remain stable. Apparently, for these the entropy loss of forming β-rolls and β-sheets is compensated for in the free-energy gain due to the hydrogen-bonding and hydrophobic interactions. We also find that the middle monomers in the trimers and tetramers are conformationally much more stable than the ones on the top and the bottom. Interestingly, the latter are more stable on the tetramer than on the trimer, suggesting that as the number of monomers increases protein-protein interactions cooperatively stabilize the assembly.
According to our simulations, the β-roll and β-sheet aggregates must be approximately equally
stable.