Abstract: Improving the efficiency and longevity of energy storage systems based on Li- and Na-ion rechargeable batteries presents a major challenge. The main problems are essentially capacity loss and limited cyclability. These effects are due to a hierarchy of factors spanning various length and time scales, interconnected in a complex manner. As a consequence, and in spite of several decades of research, a proper understanding of the ageing process has remained somewhat elusive. In recent years, however, combinations of advanced spectroscopy techniques and first-principles simulations have been applied with success to tackle this problem. In this Special Issue, we are pleased to present a selection of articles that, by precisely applying these methods, unravel key aspects of the reduction-oxidation reaction and intercalation processes. Furthermore, the approaches presented provide improvements to standard diagnostic and characterisation techniques, enabling the detection of possible Li-ion flow bottlenecks causing the degradation of capacity and cyclability.

Abstract: Condensation of spatially indirect excitons, with the electrons and holes confined in two separate layers, has recently been observed in two different double layer heterostructures. High transition temperatures were reported in a double Transition Metal Dichalcogenide (TMD) monolayer system. We briefly review electron-hole double layer systems that have been proposed as candidates for this interesting phenomenon. We investigate the double TMD system WSe2/hBN/MoSe2, using a mean-field approach that includes multiband effects due to the spin-orbit coupling and self-consistent screening of the electron-hole Coulomb interaction. We demonstrate that the transition temperature observed in the double TMD monolayers, which is remarkably high relative to the other systems, is the result of (i) the large electron and hole effective masses in TMDs, (ii) the large TMD band gaps, and (iii) the presence of multiple superfluid condensates in the TMD system. The net effect is that the superfluidity is strong across a wide range of densities, which leads to high transition temperatures that extend as high as TBKT=150 K.

Abstract: In this work, the magnetic properties of superparamagnetic iron oxide nanoparticles (SPIONs) submitted to an external magnetic field are studied using a Metropolis algorithm. The influence on the M(B) curves of the size distribution of the nanoparticles, of uniaxial anisotropy, and of dipolar interaction between the cores are examined, as well as the influence of drying the samples under a zero or non-zero magnetic field. It is shown that the anisotropy impacts the shape of the magnetization curves, which then deviate from a pure Langevin behaviour, whereas the dipolar interaction has no influence on the curves at 300 K for small particles (with a radius of 3 nm). The fitting of the magnetization curves of particles with magnetic anisotropy to a Langevin model (including a size distribution of the particles) can then lead to erroneous values of the distribution parameters. The simulation results are qualitatively compared to experimental results obtained for iron oxide nanoparticles (with a 3.21 nm median radius).

Abstract: When subjected to the interfacially induced Dzyaloshinskii-Moriya interaction, the ground state in thin ferromagnetic films with high perpendicular anisotropy is cycloidal. The period of this cycloidal state depends on the strength of the Dzyaloshinskii-Moriya interaction. In this work, we have studied the effect of confinement on the magnetic ground state and excited states, and we determined the phase diagram of thin strips and thin square platelets by means of micromagnetic calculations. We show that multiple cycloidal states with different periods can be stable in laterally confined films, where the period of the cycloids does not depend solely on the Dzyaloshinskii-Moriya interaction strength but also on the dimensions of the film. The more complex states comprising skyrmions are also found to be stable, though with higher energy.

Abstract: We study transport properties of a MoS2 monolayer in the presence of a perpendicular magnetic field B. We derive and discuss its band structure and take into account spin and valley Zeeman effects. Compared to a conventional two-dimensional electron gas, these effects lead to new quantum Hall plateaus and new peaks in the longitudinal resistivity as functions of the magnetic field. The field B leads to a significant enhancement of the spin splitting in the conduction band, to a beating of the Shubnikov-de Haas (SdH) oscillations in the low-field regime, and to their splitting in the high-field regime. The Zeeman fields suppress significantly the beating of the SdH oscillations in the low-field regime and strongly enhance their splitting at high fields. The spin and valley polarizations show a similar beating pattern at low fields and are clearly separated at high fields in which they attain a value higher than 90%.

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: We study ballistic transport of Dirac fermions in silicene through arrays of barriers, of width d, in the presence of an exchange field M and a tunable potential of height U or depth-U. The spin-and valley-resolved conductances as functions of U or M, exhibit resonances away from the Dirac point (DP) and close to it a pronounced dip that becomes a gap when a critical electric field E-z is applied. This gap widens by increasing the number of barriers and can be used to realize electric field-controlled switching of the current. The spin p(s) and valley p(v) polarizations of the current near the DP increase with Ez or M and can reach 100% for certain of their values. These field ranges widen significantly by increasing the number of barriers. Also, ps and pv oscillate nearly periodically with the separation between barriers or wells and can be inverted by reversing M.

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 structural properties of partially hydrogenated and fluorinated graphene with different percentages of H/F atoms are investigated using molecular dynamics simulations based on reactive force field (ReaxFF) potentials. We found that the roughness of graphene varies with the percentage (p) of H or F and in both cases is maximal around p = 50%. Similar results were obtained for partially oxidized graphene. The two-dimensional area size of partially fluorinated and hydrogenated graphene exhibits a local minimum around p = 35% coverage. The lattice thermal contraction in partially functionalized graphene is found to be one order of magnitude larger than that of fully covered graphene. We also show that the armchair structure for graphene oxide (similar to the structure of fully hydrogenated and fluorinated graphene) is unstable. Our results show that the structure of partially functionalized graphene changes nontrivially with the C : H and C : F ratio as well as with temperature.

Abstract: In 1976, Larkin and Ovchinnikov [Zh. Eksp. Teor. Fiz. 68, 1915 (1975) [Sov. Phys.–JETP 41, 960 (1976)]] predicted that vortex matter in superconductors driven by an electrical current can undergo an abrupt dynamic transition from a flux-flow regime to a more dissipative state at sufficiently high vortex velocities. Typically, this transition manifests itself as a large voltage jump at a particular current density, so-called instability current density J∗, which is smaller than the depairing current. By tuning the effective pinning strength in Al films, using an artificial periodic pinning array of triangular holes, we show that a unique and well-defined instability current density exists if the pinning is strong, whereas a series of multiple voltage transitions appear in the relatively weaker pinning regime. This behavior is consistent with time-dependent Ginzburg-Landau simulations, where the multiple-step transition can be unambiguously attributed to the progressive development of vortex chains and subsequently phase-slip lines. In addition, we explore experimentally the magnetic braking effects, caused by a thick Cu layer deposited on top of the superconductor, on the instabilities and the vortex ratchet effect.

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: The properties of a new class of phosphorene nanotubes (PNT) are investigated by performing first-principles calculations. We demonstrate that it is advantageous to use blue phosphorene in order to make small nanotubes and propose a way to create low-energy PNTs by the inclusion of defect lines. Five different types of defect lines are investigated and incorporated in various combinations. The resulting defect-induced faceted PNTs have negligible bending stresses which leads to a reduction in the formation energy with respect to round PNTs. Our armchair faceted PNTs have similar formation energies than the recently proposed multiphase faceted PNTs, but they have a larger variety of possible structures. Our zigzag faceted PNTs have lower formation energies than round tubes and multiphase faceted nanotubes. The electronic properties of the defect-induced faceted PNTs are determined by the defect lines which control the band gap and the shape of the electronic states at the band edges. These band gaps increase with the radius of the nanotubes and converge to those of isolated defect lines.

Abstract: Recent experimental studies showed that formation of a two-dimensional crystal structure of nitrogenated holey graphene (NHG) is possible. Similar to graphene, NHGs have an atomically thin and strong crystal structure. Using first-principles calculations, we investigate the structural, phononic, and thermal properties of monolayer NHG crystal. Our charge analysis reveals that the charged holey sites of NHG provide a reactive ground for further functionalization by adatoms or molecules. We also found that similar to graphene, the NHG structure has quite high-frequency phonon modes and the presence of nitrogen atoms leads to the emergence of additional vibrational modes. Our phonon analysis reveals the presence of three characteristic Raman-active modes of NHG. Furthermore, the analysis of constant-volume heat capacity showed that the NHG structure has a linear temperature dependence in the low-temperature region. The strong lattice structure and unique thermal properties of the NHG crystal structure are desirable in nanoscale device applications.

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: The electronic properties of monolayer phosphorene are exotic due to its puckered structure and large intrinsic direct band gap. We derive and discuss its band structure in the presence of a perpendicular magnetic field. Further, we evaluate the magneto-optical Hall and longitudinal optical conductivities as functions of temperature, magnetic field, and Fermi energy, and show that they are strongly influenced by the magnetic field. The imaginary part of the former and the real part of the latter exhibit regular interband oscillations as functions of the frequency omega in the range (h) over bar omega similar to 1.5-2 eV. Strong intraband responses in the latter and weak ones in the former occur at much lower frequencies. The magneto-optical response can be tuned in the microwave-to-terahertz and visible frequency ranges in contrast with a conventional two-dimensional electron gas or graphene in which the response is limited to the terahertz regime. This ability to isolate carriers in an anisotropic structure may make phosphorene a promising candidate for new optical devices.