Abstract: Molecular dynamics simulations are used to study the effects of different cations on the permeation of charged polymers through flat capillaries with heights below 2 nm. Interestingly, we found that, despite being monovalent, Li+ , Na+ , and K+ cations have different effects on polymer permeation, which consequently affects their transmission speed throughout those capillaries. We attribute this phenomenon to the interplay of the cations' hydration free energies and the hydrodynamic drag in front of the polymer when it enters the capillary. Different alkali cations exhibit different surface versus bulk preferences in small clusters of water under the influence of an external electric field. This paper presents a tool to control the speed of charged polymers in confined spaces using cations.

Abstract: Predicting the mechanical and thermal properties of quasi-two-dimensional (2D) transition metal dichalco-genides (TMDs) is an essential task necessary for their implementation in device applications. Although rigorous density-functional-theory-based calculations are able to predict mechanical and electronic properties, mostly they are limited to zero temperature. Classical molecular dynamics facilitates the investigation of temperature-dependent properties, but its performance highly depends on the potential used for defining interactions between the atoms. In this study, we calculated temperature-dependent phonon properties of single-layer TMDs, namely, MoS2, MoSe2, WS2, and WSe2, by utilizing Stillinger-Weber-type potentials with optimized sets of parameters with respect to the first-principles results. The phonon lifetimes and contribution of each phonon mode in thermal conductivities in these monolayer crystals are systematically investigated by means of the spectralenergy-density method based on molecular dynamics simulations. The obtained results from this approach are in good agreement with previously available results from the Green-Kubo method. Moreover, detailed analysis of lattice thermal conductivity, including temperature-dependent mode decomposition through the entire Brillouin zone, shed more light on the thermal properties of these 2D crystals. The LA and TA acoustic branches contribute most to the lattice thermal conductivity, while ZA mode contribution is less because of the quadratic dispersion around the Brillouin zone center, particularly in MoSe2 due to the phonon anharmonicity, evident from the redshift, especially in optical modes, by increasing temperature. For all the considered 2D crystals, the phonon lifetime values are compelled by transition metal atoms, whereas the group velocity spectrum is dictated by chalcogen atoms. Overall, the lattice thermal conductivity is linearly proportional with inverse temperature.

Abstract: Recently, interlayer and intralayer excitons in transition metal dichalcogenide heterobilayers have been studied both experimentally and theoretically. In spite of a growing interest, these layer-resolved excitons in the presence of external stimuli, such as strain, remain not fully understood. Here, using density-functional theory calculations with many-body effects, we explore the excitonic properties of vertically stacked MoSe2/WSe2 heterobilayer in the presence of in-plane biaxial strain of up to 5%. We calculate the strain dependence of exciton absorption spectrum, oscillator strength, wave function, and binding energy by solving the Bethe-Salpeter equation on top of the standard GW approach. We identify the interlayer and intralayer excitons by analyzing their electron-hole weights and spatial wave functions. We show that with the increase in strain magnitude, the absorption spectrum of the interlayer and intralayer excitons is red-shifted and re-ordered, and the binding energies of these layer-resolved excitons decrease monotonically and almost linearly. We derive the sensitivity of exciton binding energy to the applied strain and find that the intralayer excitons are more sensitive to strain than the interlayer excitons. For instance, a sensitivity of -7.9 meV/% is derived for the intra-MoSe2-layer excitons, which is followed by -7.4 meV/% for the intra-WSe2-layer excitons, and by -4.2 meV/% for the interlayer excitons. Our results indicate that interlayer and intralayer excitons in vertically stacked MoSe2/WSe2 heterobilayer are efficiently tunable by in-plane biaxial strain.

Abstract: Structural and electronic properties of silicon and germanium monolayers with two different diamondoid crystal structures are detailed ab initio. Our results show that, despite Si and Ge being well-known indirect gap semiconductors in their bulk form, their adamantane and diamantane monolayers can exhibit optically active direct gap in the visible frequency range, with highly anisotropic effective masses, depending on the monolayer crystal structure. Moreover, we reveal that gaps in these materials are highly tunable with applied strain. These stable monolayer forms of Si and Ge are therefore expected to help bridging the gap between the fast growing area of opto-electronics in two-dimensional materials and the established silicon-based technologies.

Abstract: The diffusion of charged particles interacting through a repulsive Yukawa potential, exp(-r/lambda)/r, confined by a parabolic potential in the y direction and subjected to a periodic substrate potential in the x direction is investigated. Langevin dynamic simulations are used to investigate the effect of the particle density, the amplitude of the periodic substrate, and the range of the interparticle interaction potential on the diffusive behavior of the particles. We found that in general the diffusion is suppressed with increasing the amplitude of the periodic potential, but for specific values of the strength of the substrate potential a remarkable increase of the diffusion is found with increasing the periodic potential amplitude. In addition, we found a strong dependence of the diffusion on the specific arrangement of the particles, e. g., single-chain versus multichain configuration. For certain particle configurations, a reentrant behavior of the diffusion is found as a function of the substrate strength due to structural transitions in the ordering of the particles.

Abstract: The giant isotope effect on the ferro- and antiferroelectric transition temperature upon deuteration of hydrogen-bonded systems is well known experimentally since various decades. Yet, theoretically only recently a microscopic understanding of this effect has been achieved which, specifically, took into account the geometry of the O ... H ... O bond. The implications of this modeling are multiple as numerous hydrogen-bonded organic systems show the same effects as ferro- and antiferroelectrics, i.e., cooperative proton tunneling at a well-defined temperature.

Abstract: We investigate how the interplay of quantum confinement and level broadening caused by disorder affects superconducting correlations in ultra-small metallic grains. We use the electron-phonon interaction-induced electron mass renormalization and the reduced static-path approximation of the BCS formalism to calculate the critical temperature as a function of the grain size. We show how the strong electron-impurity scattering additionally smears the peak structure in the electronic density of states of a metallic grain and imposes additional limits on the critical temperature under strong quantum confinement.

Abstract: The effects of strain, induced by a Gaussian bump, on the magnetic field dependent transport properties of a graphene Hall bar are investigated. The numerical simulations are performed using both classical and quantum mechanical transport theory and we found that both approaches exhibit similar characteristic features. The effects of the Gaussian bump are manifested by a decrease of the bend resistance, RB, around zero-magnetic field and the occurrence of side-peaks in RB. These features are explained as a consequence of bump-assisted scattering of electrons towards different terminals of the Hall bar. Using these features we are able to give an estimate of the size of the bump. Additional oscillations in RB are found in the quantum description that are due to the population/depopulation of Landau levels. The bump has a minor influence on the Hall resistance even for very high values of the pseudo-magnetic field. When the bump is placed outside the center of the Hall bar valley polarized electrons can be collected in the leads.

Abstract: There is an increasing need to understand interfaces between two-dimensional materials to realize an energy efficient boundary with low contact resistance and small heat dissipation. In this respect, we investigated the impact of charge and substitutional atom doping on the electronic transport properties of the hybrid metallic-semiconducting lateral junctions, formed between metallic (1T and 1T(d)) and semiconducting (1H) phases of MoS2 by means of first-principles and non-equilibrium Green function formalism based calculations. Our results clearly revealed the strong influence of the type of interface and crystallographic orientation of the metallic phase on the transport properties of these systems. The Schottky barrier height, which is the dominant mechanism for contact resistance, was found to be as large as 0.63 eV and 1.19 eV for holes and electrons, respectively. We found that armchair interfaces are more conductive as compared to zigzag termination due to the presence of the metallic Mo zigzag chains that are directed along the transport direction. In order to manipulate these barrier heights we investigated the influence of electron doping of the metallic part (i.e. 1T(d) -MoS2). We observed that the Fermi level of the hybrid system moves towards the conduction band of semiconducting 1H-MoS2 due to filling of 4d-orbital of metallic MoS2, and thus the Schottky barrier for electrons decreases considerably. Besides electron doping, we also investigated the effect of substitutional doping of metallic MoS2 by replacing Mo atoms with either Re or Ta. Due to its valency, Re (Ta) behaves as a donor (acceptor) and reduces the Schottky barrier for electrons (holes). Since Re and Ta based transition metal dichalcogenides crystallize in either the 1T(d) or 1T phase, substitutional doping with these atom favors the stabilization of the 1T(d) phase of MoS2. Co-doping of hybrid structure results in an electronic structure, which facilities easy dissociation of excitons created in the 1H part.

Abstract: We present a theoretical investigation of the plasmon-polariton modes in grating coupled graphene inside a Fabry-Perot cavity. The cavity or photon modes of the device are determined by the Finite Difference Time Domain (FDTD) simulations and the corresponding plasmon-polariton modes are obtained by applying a many-body self-consistent field theory. We find that in such a device structure, the electric field strength of the incident electromagnetic (EM) field can be significantly enhanced near the edges of the grating strips. Thus, the strong coupling between the EM field and the plasmons in graphene can be achieved and the features of the plasmon-polariton oscillations in the structure can be observed. It is found that the frequencies of the plasmon-polariton modes are in the terahertz (THz) bandwidth and depend sensitively on electron density which can be tuned by applying a gate voltage. Moreover, the coupling between the cavity photons and the plasmons in graphene can be further enhanced by increasing the filling factor of the device. This work can help us to gain an in-depth understanding of the THz plasmonic properties of graphene-based structures.

Abstract: Using the anisotropic time-dependent Ginzburg-Landau theory we study the effect of ordered and disordered pinning on the time response of superconducting strips to an external current that switched on abruptly. The pinning centers result in a considerable delay of the response time of the system to such abrupt switching on of the current, whereas the output voltage is always larger when pinning is present. The resistive state in both cases are characterized either by dynamically stable phase-slip centers/lines or expanding in-time hot-spots, which are the main mechanisms for dissipation in current-carrying superconductors. We find that hot-spots are always initiated by the phase-slip state. However, the range of the applied current for the phase-slip state increases significantly when pinning is introduced. Qualitative changes are observed in the dynamics of the superconducting condensate in the presence of pinning.

Abstract: The melting of a binary system of charged particles confined in a quasi-one-dimensional parabolic channel is studied through Monte Carlo simulations. At zero temperature the particles are ordered in parallel chains. The melting is anisotropic and different melting temperatures are obtained according to the spatial direction, and the different kinds of particles present in the system. Melting is very different for the single-, two- and four-chain configurations. A temperature induced structural phase transition is found between two different four-chain ordered states which is absent in the mono-disperse system. In the mixed regime, where the two kinds of particles are only slightly different, melting is almost isotropic and a thermally induced homogeneous distribution of the distinct kinds of charges is observed.

Abstract: Recent progress in materials synthesis enabled fabrication of superconducting atomically flat single-crystalline metallic nanofilms with thicknesses down to a few monolayers. Interest in such nano-thin systems is attracted by the dimensional 3D-2D crossover in their coherent properties which occurs with decreasing the film thickness. The first fundamental aspect of this crossover is dictated by the Mermin-Wagner-Hohenberg theorem and concerns frustration of the long-range order due to superconductive fluctuations and the possibility to track its impact with an unprecedented level of control. The second important aspect is related to the Fabri-Perot modes of the electronic motion strongly bound in the direction perpendicular to the nanofilm. The formation of such modes results in a pronounced multiband structure that changes with the nanofilm thickness and affects both the mean-field behavior and superconductive fluctuations. Though the subject is very rich in physics, it is scarcely investigated to date. The main obstacle is that there are no manageable models to study a complex magnetic response in this case. Full microscopic consideration is rather time consuming, if practicable at all, while the standard Ginzburg-Landau theory is not applicable. In the present work we review the main achievements in the subject to date, and construct and justify an efficient multiband mean-field formalism which allows for numerical and even analytical treatment of nano-thin superconductors in applied magnetic fields.

Abstract: The valley splitting of the conduction-band states in a thin silicon-on-insulator layer is investigated using the 30-band k . p theory. The system composed of a few nm thick Si layer embedded within thick SiO2 layers is analyzed. The valley split states are found to cross periodically with increasing quantum well width, and therefore the energy splitting is an oscillatory function of the quantum well width, with period determined by the wave vector K-0 of the conduction band minimum. Because the valley split states are classified by parity, the optical transition between the ground hole state and one of those valley split conduction band states is forbidden. The oscillations in the valley splitting energy decrease with electric field and with smoothing of the composition profile between the well and the barrier by diffusion of oxygen from the SiO2 layers to the Si quantum well. Such a smoothing also leads to a decrease of the interband transition matrix elements. The obtained results are well parametrized by the effective two-valley model, but are found to disagree from previous 30-band calculations. This discrepancy could be traced back to the fact that the basis for the numerical solution of the eigenproblem must be restricted to the first Brillouin zone in order to obtain quantitatively correct results for the valley splitting.

Abstract: We compare classical versus quantum electron transport in recently investigated magnetic focusing devices (Bhandari et al 2016 Nano Lett. 16 1690) exposed to the perturbing potential of a scanning gate microscope (SGM). Using the Landauer-Buttiker formalism for a multi-terminal device, we calculate resistance maps that are obtained as the SGM tip is scanned over the sample. There are three unique regimes in which the scanning tip can operate (focusing, repelling, and mixed regime) which are investigated. Tip interacts mostly with electrons with cyclotron trajectories passing directly underneath it, leaving a trail of modified current density behind it. Other (indirect) trajectories become relevant when the tip is placed near the edges of the sample, and current is scattered between the tip and the edge. We point out that, in contrast to SGM experiments on gapped semiconductors, the STM tip can induce a pn junction in graphene, which improves contrast and resolution in SGM. We also discuss possible explanations for spatial asymmetry of experimentally measured resistance maps, and connect it with specific configurations of the measuring probes.