Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userA hierarchy of self-consistent stochastic boundary conditions for Ising lattice simulations
151
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
We describe a hierarchy of stochastic boundary conditions (SBCs) that can be used to systematically eliminate finite size effects in Monte Carlo simulations of Ising lattices. For an Ising model on a $100 times 100$ square lattice, we measured the specific heat, the magnetic susceptibility, and the spin-spin correlation using SBCs of the two lowest orders, to show that they compare favourably against periodic boundary conditions (PBC) simulations and analytical results. To demonstrate how versatile the SBCs are, we then simulated an Ising lattice with a magnetized boundary, and another with an open boundary, measuring the magnetization, magnetic susceptibility, and longitudinal and transverse spin-spin correlations as a function of distance from the boundary.
rate research
Read More
We present the complete set of stochastic Verlet-type algorithms that can provide correct statistical measures for both configurational and kinetic sampling in discrete-time Langevin systems. The approach is a brute-force general representation of the Verlet-algorithm with free parameter coefficients that are determined by requiring correct Boltzmann sampling for linear systems, regardless of time step. The result is a set of statistically correct methods given by one free functional parameter, which can be interpreted as the one-time-step velocity attenuation factor. We define the statistical characteristics of both true on-site $v^n$ and true half-step $u^{n+frac{1}{2}}$ velocities, and use these definitions for each statistically correct Stormer-Verlet method to find a unique associated half-step velocity expression, which yields correct kinetic Maxwell-Boltzmann statistics for linear systems. It is shown that no other similar, statistically correct on-site velocity exists. We further discuss the use and features of finite-difference velocity definitions that are neither true on-site, nor true half-step. The set of methods is written in convenient and conventional stochastic Verlet forms that lend themselves to direct implementation for, e.g., Molecular Dynamics applications. We highlight a few specific examples, and validate the algorithms through comprehensive Langevin simulations of both simple nonlinear oscillators and complex Molecular Dynamics.
Multifluid simulations of plasma sheaths are increasingly used to model a wide variety of problems in plasma physics ranging from global magnetospheric flows around celestial bodies to plasma-wall interactions in thrusters and fusion devices. For multifluid problems, accurate boundary conditions to model an absorbing wall that resolves a classical sheath remains an open research area. This work justifies the use of vacuum boundary conditions for absorbing walls to show comparable accuracy between a multifluid sheath and lower moments of a continuum-kinetic sheath.
We present a generalization of Blochs theorem to finite-range lattice systems of independent fermions, in which translation symmetry is broken only by arbitrary boundary conditions, by providing exact, analytic expressions for all energy eigenvalues and eigenstates. By transforming the single-particle Hamiltonian into a corner-modified banded block-Toeplitz matrix, a key step is a bipartition of the lattice, which splits the eigenvalue problem into a system of bulk and boundary equations. The eigensystem inherits most of its solutions from an auxiliary, infinite translation-invariant Hamiltonian that allows for non-unitary representations of translation symmetry. A reformulation of the boundary equation in terms of a boundary matrix ensures compatibility with the boundary conditions, and determines the allowed energy eigenstates. We show how the boundary matrix captures the interplay between bulk and boundary properties, leading to efficient indicators of bulk-boundary correspondence. Remarkable consequences of our generalized Bloch theorem are the engineering of Hamiltonians that host perfectly localized, robust zero-energy edge modes, and the predicted emergence, e.g. in Kitaevs chain, of localized excitations whose amplitudes decay exponentially with a power-law prefactor. We further show how the theorem yields diagonalization algorithms for the class of Hamiltonians under consideration, and use the proposed bulk-boundary indicator to characterize the topological response of a multi-band time-reversal invariant s-wave superconductor under twisted boundary conditions, showing how a fractional Josephson effect can occur without a fermionic parity switch. Finally, we establish connections to the transfer matrix method and demonstrate, using the paradigmatic Kitaevs chain example, that a non-diagonalizable transfer matrix signals the presence of solutions with a power-law prefactor.
Phonon lifetime calculations from first principles usually rely on time consuming molecular dynamics calculations, or density functional perturbation theory (DFPT) where the zero temperature crystal structure is assumed to be dynamically stable. Here a new and effective method for calculating phonon lifetimes from first principles is presented, not limited to crystal structures stable at 0 K, and potentially much more effective than most corresponding molecular dynamics calculations. The method is based on the recently developed self consistent lattice dynamical method and is here tested by calculating the bcc phase phonon lifetimes of Li, Na, Ti and Zr, as representative examples.
Using the parallel tempering algorithm and GPU accelerated techniques, we have performed large-scale Monte Carlo simulations of the Ising model on a square lattice with antiferromagnetic (repulsive) nearest-neighbor(NN) and next-nearest-neighbor(NNN) interactions of the same strength and subject to a uniform magnetic field. Both transitions from the (2x1) and row-shifted (2x2) ordered phases to the paramagnetic phase are continuous. From our data analysis, reentrance behavior of the (2x1) critical line and a bicritical point which separates the two ordered phases at T=0 are confirmed. Based on the critical exponents we obtained along the phase boundary, Suzukis weak universality seems to hold.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
