Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userSimulation of 1+1 dimensional surface growth and lattices gases using GPUs
208
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
Restricted solid on solid surface growth models can be mapped onto binary lattice gases. We show that efficient simulation algorithms can be realized on GPUs either by CUDA or by OpenCL programming. We consider a deposition/evaporation model following Kardar-Parisi-Zhang growth in 1+1 dimensions related to the Asymmetric Simple Exclusion Process and show that for sizes, that fit into the shared memory of GPUs one can achieve the maximum parallelization speedup ~ x100 for a Quadro FX 5800 graphics card with respect to a single CPU of 2.67 GHz). This permits us to study the effect of quenched columnar disorder, requiring extremely long simulation times. We compare the CUDA realization with an OpenCL implementation designed for processor clusters via MPI. A two-lane traffic model with randomized turning points is also realized and the dynamical behavior has been investigated.
rate research
Read More
We implement several symplectic integrators, which are based on two part splitting, for studying the chaotic behavior of one- and two-dimensional disordered Klein-Gordon lattices with many degrees of freedom and investigate their numerical performance. For this purpose, we perform extensive numerical simulations by considering many different initial energy excitations and following the evolution of the created wave packets in the various dynamical regimes exhibited by these models. We compare the efficiency of the considered integrators by checking their ability to correctly reproduce several features of the wave packets propagation, like the characteristics of the created energy distribution and the time evolution of the systems maximum Lyapunov exponent estimator. Among the tested integrators the fourth order $ABA864$ scheme cite{BCFLMM13} showed the best performance as it needed the least CPU time for capturing the correct dynamical behavior of all considered cases when a moderate accuracy in conserving the systems total energy value was required. Among the higher order schemes used to achieve a better accuracy in the energy conservation, the sixth order scheme $s11ABA82_6$ exhibited the best performance.
With Monte Carlo methods, we investigate the universality class of the depinning transition in the two-dimensional Ising model with quenched random fields. Based on the short-time dynamic approach, we accurately determine the depinning transition field and both static and dynamic critical exponents. The critical exponents vary significantly with the form and strength of the random fields, but exhibit independence on the updating schemes of the Monte Carlo algorithm. From the roughness exponents $zeta, zeta_{loc}$ and $zeta_s$, one may judge that the depinning transition of the random-field Ising model belongs to the new dynamic universality class with $zeta eq zeta_{loc} eq zeta_s$ and $zeta_{loc} eq 1$. The crossover from the second-order phase transition to the first-order one is observed for the uniform distribution of the random fields, but it is not present for the Gaussian distribution.
We show how U(1) lattice gauge theories display key signatures of ergodicity breaking in the presence of a random charge background. Contrary to the widely studied case of spin models, in the presence of Coulomb interactions, the spectral properties of such lattice gauge theories are very weakly affected by finite-volume effects. This allows to draw a sharp boundary for the ergodic regime, and thus the breakdown of quantum chaos for sufficiently strong gauge couplings, at the system sizes accessible via exact diagonalization. Our conclusions are independent on the value of a background topological angle, and are contrasted with a gauge theory with truncated Hilbert space, where instead we observe very strong finite-volume effects akin to those observed in spin chains.
The characterizing feature of a many-body localized phase is the existence of an extensive set of quasi-local conserved quantities with an exponentially localized support. This structure endows the system with the signature logarithmic in time entanglement growth between spatial partitions. This feature differentiates the phase from Anderson localization, in a non-interacting model. Experimentally measuring the entanglement between large partitions of an interacting many-body system requires highly non-local measurements which are currently beyond the reach of experimental technology. In this work we demonstrate that the defining structure of many-body localization can be detected by the dynamics of a simple quantity from quantum information known as the total correlations which is connected to the local entropies. Central to our finding is the necessity to propagate specific initial states, drawn from the Hamiltonian unbiased basis (HUB). The dynamics of the local entropies and total correlations requires only local measurements in space and therefore is potentially experimentally accessible in a range of platforms.
Differentiable programming has emerged as a key programming paradigm empowering rapid developments of deep learning while its applications to important computational methods such as Monte Carlo remain largely unexplored. Here we present the general theory enabling infinite-order automatic differentiation on expectations computed by Monte Carlo with unnormalized probability distributions, which we call automatic differentiable Monte Carlo (ADMC). By implementing ADMC algorithms on computational graphs, one can also leverage state-of-the-art machine learning frameworks and techniques to traditional Monte Carlo applications in statistics and physics. We illustrate the versatility of ADMC by showing some applications: fast search of phase transitions and accurately finding ground states of interacting many-body models in two dimensions. ADMC paves a promising way to innovate Monte Carlo in various aspects to achieve higher accuracy and efficiency, e.g. easing or solving the sign problem of quantum many-body models through ADMC.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
