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: We investigate the self-assembly of vortices in a type-II superconducting disk subjected to highly nonuniform confining potentials produced by inhomogeneous magnetic textures. Using a series of numerical experiments performed within the Ginzburg–Landau theory, we show that vortices can arrange spontaneously in highly nonuniform, defect-free crystals, reminiscent of conformal lattices, even though the strict conditions for the conformal crystal are not fulfilled. These results contradict continuum-limit theory, which predicts that the order of a nonuniform crystal is unavoidably frustrated by the presence of topological defects. By testing different cooling routes of the superconductor, we observed several different self-assembled configurations, each of which corresponding to one in a set of allowed conformal transformations, which depends on the magnetic and thermal histories of the system.

Abstract: A unified theory of phonon dispersions and piezoelectricity in bulk and multilayers of hexagonal boron nitride (h-BN) is derived. The dynamical matrix is calculated on the basis of an empirical force constant model of intralayer valence and interlayer van der Waals interactions. Coulomb interactions are calculated by Ewalds method, adapted for the three-dimensional (3D) and the multilayer case. The deformation of the ionic charge distribution with long-wave lattice displacements is taken into account. Special attention is devoted to the nonanalytic long-range Coulomb contribution to the dynamical matrix which is different for the 3D crystal and the multilayer case. Consequently there is a splitting of the transverse optical (TO) and longitudinal optical (LO) phonon branches of E1u symmetry and a discontinuity of the A2u branch at the Γ point in 3D h-BN. No such splitting and discontinuity at Γ are present in multilayer crystals with a finite number N of layers. There a diverging bundle of N overbending optical phonon branches emerges from Γ. Borns long-wave theory is applied and extended for the study of piezoelectricity in layered crystals. While 3D h-BN and h-BN multilayers with an even number of layers (symmetry D6h) are not piezoelectric, multilayers with an uneven number of Nu layers (symmetry D3h) are piezoelectric; the piezoelectric coefficient e1,11 is inversely proportional to Nu.

Abstract: We study magnetic flux interacting with arrays of pinning sites (APSs) placed on vertices of hyperbolic tesselations (HTs). We show that, due to the gradient in the density of pinning sites, HT APSs are capable of trapping vortices for a broad range of applied magnetic fluxes. Thus, the penetration of magnetic field in HT APSs is essentially different from the usual scenario predicted by the Bean model. We demonstrate that, due to the enhanced asymmetry of the surface barrier for vortex entry and exit, this HT APS could be used as a “capacitor” to store magnetic flux.