Abstract: Monte Carlo (MC) methods have a long-standing history as partners of molecular dynamics (MD) to simulate the evolution of materials at the atomic scale. Among these techniques, the uniform-acceptance force-bias Monte Carlo (UFMC) method [ G. Dereli Mol. Simul. 8 351 (1992)] has recently attracted attention [ M. Timonova et al. Phys. Rev. B 81 144107 (2010)] thanks to its apparent capacity of being able to simulate physical processes in a reduced number of iterations compared to classical MD methods. The origin of this efficiency remains, however, unclear. In this work we derive a UFMC method starting from basic thermodynamic principles, which leads to an intuitive and unambiguous formalism. The approach includes a statistically relevant time step per Monte Carlo iteration, showing a significant speed-up compared to MD simulations. This time-stamped force-bias Monte Carlo (tfMC) formalism is tested on both simple one-dimensional and three-dimensional systems. Both test-cases give excellent results in agreement with analytical solutions and literature reports. The inclusion of a time scale, the simplicity of the method, and the enhancement of the time step compared to classical MD methods make this method very appealing for studying the dynamics of many-particle systems.

Abstract: The superconducting state. of a thin superconducting disk with a hole is studied within the, nonlinear Ginzburg-Landau theory in which the demagnetization effect is accurately taken into account. We find that the flux through the hole is not quantized, the superconducting state is stabilized with increasing size of the hole for fixed radius of the disk, and a transition to a multivortex state is found if the disk is sufficiently large. Breaking the circular symmetry through a non-central-location of the hole in the disk favors the multivortex state.

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: Recently, core-shell Ag-Co, Ag-Cu, and “onionlike” Cu-Co equilibrium configurations were predicted in the case of isolated face centered cubic (fcc) bimetallic clusters, and three shell onionlike configurations were predicted in the case of ternary metallic clusters with spherical and truncated octahedral morphologies. In the present paper, immiscible binary CuCo and ternary AgCuCo clusters with icosahedral structures are studied as functions of their size and composition. Clusters studied are formed by 13, 55, 147, 309, and 561 atoms corresponding to the five smallest possible closed shell icosahedral structures. An embedded atom model potential is used to describe their cohesion. Equilibrium configurations are investigated by means of Metropolis Monte Carlo free energy minimization in the (NPT) canonical ensemble. Most simulations are achieved at 10 and 300 K. The effect of temperature on segregation ordering is systematically investigated. Selected cases are used to identify the effect of size and composition on melting. In contrast with fcc clusters, homogeneous onionlike configurations of binary clusters are not predicted. When it is allowed by the composition, a complete outer shell is formed by Cu in binary Cu-Co clusters and by Ag in ternary Ag-Cu-Co clusters. Depending on temperature, Co may precipitate into decahedral groups under the Cu vertices of the icosahedra in binary clusters, while the Co-Cu configuration in ternary clusters drastically depends on the Ag coating. Despite the multicomponent character of the clusters and the immiscibility of the species forming them, for most compositions and sizes, equilibrium structures remain close to perfectly icosahedral at 10 K as well as at 300 K.

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: We study the effect of structural inversion asymmetry, induced by the presence of substrates or by external electric fields, on charge transport in magnetic topological ultra-thin films. We consider general orientations of the magnetic impurities. Our results are based on the Boltzmann formalism along with a modified relaxation time scheme. We show that the structural inversion asymmetry enhances the charge transport anisotropy induced by the magnetic impurities and when only one conduction subband contributes to the charge transport a dissipationless charge current is accessible. We demonstrate how a substrate or gate voltage can control the effect of the magnetic impurities on the charge transport, and how this depends on the orientation of the magnetic impurities.