Abstract: We modeled the electron and hole states in Si/SiO2 quantum wells within a basis of standing waves using the 30-band k . p theory. The hard-wall confinement potential is assumed, and the influence of the peculiar band structure of bulk silicon on the quantum-well sub-bands is explored. Numerous spurious solutions in the conduction-band and valence-band energy spectra are found and are identified to be of two types: (1) spurious states which have large contributions of the bulk solutions with large wave vectors (the high-k spurious solutions) and (2) states which originate mainly from the spurious valley outside the Brillouin zone (the extravalley spurious solutions). An algorithm to remove all those nonphysical solutions from the electron and hole energy spectra is proposed. Furthermore, slow and oscillatory convergence of the hole energy levels with the number of basis functions is found and is explained by the peculiar band mixing and the confinement in the considered quantum well. We discovered that assuming the hard-wall potential leads to numerical instability of the hole states computation. Nonetheless, allowing the envelope functions to exponentially decay in a barrier of finite height is found to improve the accuracy of the computed hole states.

Abstract: Recently penta-graphene was proposed as a stable two-dimensional carbon allotrope consisting of a single layer of interconnected carbon pentagons [Zhang et al., PNAS, 2015, 112, 2372]. Its silicon counterpart, penta-silicene, however, is not stable. In this work, we show that multilayers of penta-silicene form stable materials with semiconducting or metallic properties, depending on the stacking mode. We demonstrate their dynamic stability through their phonon spectrum and using molecular dynamics. A particular type of bilayer penta-silicene is found to have lower energy than all of the known hexagonal silicene bilayers and forms therefore the most stable bilayer silicon material predicted so far. The electronic and mechanical properties of these new silicon allotropes are studied in detail and their behavior under strain is investigated. We demonstrate that strain can be used to tune its band gap.

Abstract: Nanoparticle charging in a capacitively coupled radio frequency discharge in argon is studied using a particle in cell Monte Carlo collisions method. The plasma parameters and dust potential were calculated self-consistently for different unmovable dust profiles. A new method for definition of the dust floating potential is proposed, based on the information about electron and ion energy distribution functions, obtained during the kinetic simulations. This approach provides an accurate balance of the electron and ion currents on the dust particle surface and allows us to precisely calculate the dust floating potential. A comparison of the obtained floating potentials with the results of the traditional orbital motion limit (OML) theory shows that in the presence of the ion resonant charge exchange collisions, even when the OML approximation is valid, its results are correct only in the region of a weak electric field, where the ion drift velocity is much smaller than the thermal one. With increasing ion drift velocity, the absolute value of the calculated dust potential becomes significantly smaller than the theory predicts. This is explained by a non-Maxwellian shape of the ion energy distribution function for the case of fast ion drift.

Abstract: A simplified quantum mechanical model is developed to investigate quantum transport features such as the electron concentration and the current flowing through a silicon nanowire metal-oxide-semiconductor field-effect transistor (MOSFET). In particular, the electron concentration is extracted from a self-consistent solution of the Schrödinger and Poisson equations as well as the ballistic Boltzmann equation which have been solved by exploiting a nonlinear variational principle within the framework of the generalized local density approximation. A suitable action functional has been minimized and details of the implementation and its numerical minimization are given. The current density and its related current-voltage characteristics are calculated from the one-dimensional ballistic steady-state Boltzmann transport equation which is solved analytically by using the method of characteristic curves. The straightforward implementation, the computational speed and the good qualitative behavior of the transport characteristics observed in our approach make it a promising simulation method for modeling quantum transport in nanowire MOSFETs.

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: 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: Dielectric properties of nano-confined water are important in several areas of science, i.e., it is relevant in the dielectric double layer that exists in practically all heterogeneous fluid-based systems. Molecular dynamics simulations are used to predict the in-plane dielectric properties of confined water in planar channels of width ranging from sub-nanometer to bulk. Because of suppressed rotational degrees of freedom near the confining walls, the dipole of the water molecules tends to be aligned parallel to the walls, which results in a strongly enhanced in-plane dielectric constant (epsilon (parallel to)) reaching values of about 120 for channels with height 8 angstrom < h < 10 angstrom. With the increase in the width of the channel, we predict that epsilon (parallel to) decreases nonlinearly and reaches the bulk value for h > 70 angstrom. A stratified continuum model is proposed that reproduces the h > 10 angstrom dependence of epsilon (parallel to). For sub-nanometer height channels, abnormal behavior of epsilon (parallel to) is found with two orders of magnitude reduction of epsilon (parallel to) around h similar to 7.5 angstrom, which is attributed to the formation of a particular ice phase that exhibits long-time (similar to mu s) stable ferroelectricity. This is of particular importance for the understanding of the influence of confined water on the functioning of biological systems.

Abstract: Hydrodynamic flow in two-dimensional electron systems has so far been probed only by dc transport and scanning gate microscopy measurements. In this work we discuss theoretically signatures of the hydrodynamic regime in near-field optical microscopy. We analyze the dispersion of acoustic plasmon modes in two-dimensional electron liquids using a nonlocal conductivity that takes into account the effects of (momentumconserving) electron-electron collisions, (momentum-relaxing) electron-phonon and electron-impurity collisions, and many-body interactions beyond the celebrated random phase approximation. We derive the dispersion and, most importantly, the damping of acoustic plasmon modes and their coupling to a near-field probe, identifying key experimental signatures of the crossover between collisionless and hydrodynamic regimes.

Abstract: Ab initio calculations within the density-functional theory formalism are performed to investigate the chemical functionalization of a graphene-like monolayer of siliconsilicenewith B, N, Al, or P atoms. The structural, electronic, magnetic, and vibrational properties are reported. The most preferable adsorption sites are found to be valley, bridge, valley and hill sites for B, N, Al, and P adatoms, respectively. All the relaxed systems with adsorbed/substituted atoms exhibit metallic behavior with strongly bonded B, N, Al, and P atoms accompanied by an appreciable electron transfer from silicene to the B, N, and P adatom/substituent. The Al atoms exhibit opposite charge transfer, with n-type doping of silicene and weaker bonding. The adatoms/substituents induce characteristic branches in the phonon spectrum of silicene, which can be probed by Raman measurements. Using molecular dynamics, we found that the systems under study are stable up to at least T=500 K. Our results demonstrate that silicene has a very reactive and functionalizable surface.

Abstract: Motivated by the recent realization of graphene sensors to detect individual gas molecules, we investigate the adsorption of H2O, NH3, CO, NO2, and NO on a graphene substrate using first-principles calculations. The optimal adsorption position and orientation of these molecules on the graphene surface is determined and the adsorption energies are calculated. Molecular doping, i.e., charge transfer between the molecules and the graphene surface, is discussed in light of the density of states and the molecular orbitals of the adsorbates. The efficiency of doping of the different molecules is determined and the influence of their magnetic moment is discussed.

Abstract: We investigate the adsorption process of small molecules on graphene through first-principles calculations and show the presence of two main charge transfer mechanisms. Which mechanism is the dominant one depends on the magnetic properties of the adsorbing molecules. We explain these mechanisms through the density of states of the system and the molecular orbitals of the adsorbates, and demonstrate the possible difficulties in calculating the charge transfer from first principles between a graphene sheet and a molecule. Our results are in good agreement with experiment.

Abstract: The adsorption and diffusion of Na, K, and Ca atoms on MXene/graphene heterostructures of MXene systems Sc2C(OH)(2), Ti2CO2, and V2CO2 are systematically investigated by using first-principles methods. We found that alkali metal intercalation is energetically favorable and thermally stable for Ti2CO2/graphene and V2CO2/graphene heterostructures but not for Sc2C(OH)(2). Diffusion kinetics calculations showed the advantage of MXene/graphene heterostructures over sole MXene systems as the energy barriers are halved for the considered alkali metals. Low energy barriers are found for Na and K ions, which are promising for fast charge/discharge rates. Calculated voltage profiles reveal that estimated high capacities can be fully achieved for Na ion in V2CO2/graphene and Ti2CO2/graphene heterostructures. Our results indicate that Ti2CO2/graphene and V2CO2/graphene electrode materials are very promising for Na ion battery applications. The former could be exploited for low voltage applications while the latter will be more appropriate for higher voltages.

Abstract: Van der Waals (vdW) heterostructures continue to attract intense interest as a route of designing materials with novel properties that cannot be found in nature. Unfortunately, this approach is currently limited to only a few layers that can be stacked on top of each other. Here, we report a bulk vdW material consisting of superconducting 1H TaS2 monolayers interlayered with 1T TaS2 monolayers displaying charge density waves (CDW). This bulk vdW heterostructure is created by phase transition of 1T-TaS2 to 6R at 800 degrees C in an inert atmosphere. Its superconducting transition (T-c) is found at 2.6 K, exceeding the T-c of the bulk 2H phase. Using first-principles calculations, we argue that the coexistence of superconductivity and CDW within 6R-TaS2 stems from amalgamation of the properties of adjacent 1H and 1T monolayers, where the former dominates the superconducting state and the latter the CDW behavior.

Abstract: The electron states in axially symmetric quantum wires are computed by means of the effective-mass Schrodinger equation, which is written in cylindrical coordinates phi, rho, and z. We show that a direct discretization of the Schrodinger equation by central finite differences leads to a non-symmetric Hamiltonian matrix. Because diagonalization of such matrices is more complex it is advantageous to transform it in a symmetric form. This can be done by the Liouville-like transformation proposed by Rizea et al. (2008), which replaces the wave function psi(rho) with the function F(rho) = psi(rho)root rho and transforms the Hamiltonian accordingly. Even though a symmetric Hamiltonian matrix is produced by this procedure, the computed wave functions are found to be inaccurate near the origin, and the accuracy of the energy levels is not very high. In order to improve on this, we devised a finite-difference scheme which discretizes the Schrodinger equation in the first step, and then applies the Liouville-like transformation to the difference equation. Such a procedure gives a symmetric Hamiltonian matrix, resulting in an accuracy comparable to the one obtained with the finite element method. The superior efficiency of the new finite-difference scheme (FDM) is demonstrated for a few p-dependent one-dimensional potentials which are usually employed to model the electron states in free-standing and core shell quantum wires. The new scheme is compared with the other FDM schemes for solving the effective-mass Schrodinger equation, and is found to deliver energy levels with much smaller numerical error for all the analyzed potentials. It also gives more accurate results than the scheme of Rizea et al., except for the ground state of an infinite rectangular potential in freestanding quantum wires. Moreover, the PT symmetry is invoked to explain similarities and differences between the considered FDM schemes. (C) 2015 Elsevier B.V. All rights reserved.