اشترك بالحزمة الذهبية واحصل على وصول غير محدود شمرا أكاديميا
تسجيل مستخدم جديدSimulating graphene impurities using the worm algorithm
53
0
0.0
(
0
)
تأليف
Marcin Szyniszewski
اسأل ChatGPT حول البحث
ﻻ يوجد ملخص باللغة العربية
The two-dimensional Ising model is studied by performing computer simulations with one of the Monte Carlo algorithms - the worm algorithm. The critical temperature T_C of the phase transition is calculated by the usage of the critical exponents and the results are compared to the analytical result, giving a very high accuracy. We also show that the magnetic ordering of impurities distributed on a graphene sheet is possible, by simulating the properly constructed model using the worm algorithm. The value of T_C is estimated. Furthermore, the dependence of T_C on the interaction constants is explored. We outline how one can proceed in investigating this relation in the future.
قيم البحث
اقرأ أيضاً
We observe a localized phase of ultracold bosonic quantum gases in a 3-dimensional optical lattice induced by a small contribution of fermionic atoms acting as impurities in a Fermi-Bose quantum gas mixture. In particular we study the dependence of t
his transition on the fermionic 40K impurity concentration by a comparison to the corresponding superfluid to Mott insulator transition in a pure bosonic 87Rb gas and find a significant shift in the transition parameter. The observed shift is larger than expected based on a mean-field argument, which is a strong indication that disorder-related effects play a significant role.
Quantum spin systems with strong geometric restrictions give rise to rich quantum phases such as valence bond solids and spin liquid states. However, the geometric restrictions often hamper the application of sophisticated numerical approaches. Based
on the stochastic series expansion method, we develop an efficient and exact quantum Monte Carlo sweeping cluster algorithm which automatically satisfies the geometrical restrictions. Here we use the quantum dimer model as a benchmark to demonstrate the reliability and power of this algorithm. Comparing to existing numerical methods, we can obtain higher accuracy results for a wider parameter region and much more substantial system sizes.
Analyzing in detail the first corrections to the scaling hypothesis, we develop accelerated methods for the determination of critical points from finite size data. The output of these procedures are sequences of pseudo-critical points which rapidly c
onverge towards the true critical points. In fact more rapidly than previously existing methods like the Phenomenological Renormalization Group approach. Our methods are valid in any spatial dimensionality and both for quantum or classical statistical systems. Having at disposal fast converging sequences, allows to draw conclusions on the basis of shorter system sizes, and can be extremely important in particularly hard cases like two-dimensional quantum systems with frustrations or when the sign problem occurs. We test the effectiveness of our methods both analytically on the basis of the one-dimensional XY model, and numerically at phase transitions occurring in non integrable spin models. In particular, we show how a new Homogeneity Condition Method is able to locate the onset of the Berezinskii-Kosterlitz-Thouless transition making only use of ground-state quantities on relatively small systems.
The simplest model of DNA mechanics describes the double helix as a continuous rod with twist and bend elasticity. Recent work has discussed the relevance of a little-studied coupling $G$ between twisting and bending, known to arise from the groove a
symmetry of the DNA double helix. Here, the effect of $G$ on the statistical mechanics of long DNA molecules subject to applied forces and torques is investigated. We present a perturbative calculation of the effective torsional stiffness $C_text{eff}$ for small twist-bend coupling. We find that the bare $G$ is screened by thermal fluctuations, in the sense that the low-force, long-molecule effective free energy is that of a model with $G=0$, but with long-wavelength bending and twisting rigidities that are shifted by $G$-dependent amounts. Using results for torsional and bending rigidities for freely-fluctuating DNA, we show how our perturbative results can be extended to a non-perturbative regime. These results are in excellent agreement with numerical calculations for Monte Carlo triad and molecular dynamics oxDNA models, characterized by different degrees of coarse-graining, validating the perturbative and non-perturbative analyses. While our theory is in generally-good quantitative agreement with experiment, the predicted torsional stiffness does systematically deviate from experimental data, suggesting that there are as-yet-uncharacterized aspects of DNA twisting-stretching mechanics relevant to low-force, long-molecule mechanical response, which are not captured by widely-used coarse-grained models.
We study the distribution of lengths and other statistical properties of worms constructed by Monte Carlo worm algorithms in the power-law three-sublattice ordered phase of frustrated triangular and kagome lattice Ising antiferromagnets. Viewing each
step of the worm construction as a position increment (step) of a random walker, we demonstrate that the persistence exponent $theta$ and the dynamical exponent $z$ of this random walk depend only on the universal power-law exponents of the underlying critical phase, and not on the details of the worm algorithm or the microscopic Hamiltonian. Further, we argue that the detailed balance criterion obeyed by such worm algorithms and the power-law correlations of the underlying equilibrium system together give rise to two related properties of this random walk: First, the steps of the walk are expected to be power-law correlated in time. Second, the position distribution of the walker relative to its starting point is given by the equilibrium position distribution of a particle in an attractive logarithmic central potential of strength $eta_m$, where $eta_m$ is the universal power-law exponent of the equilibrium defect-antidefect correlation function of the underlying spin system. We derive a scaling relation, $z = (2-eta_m)/(1-theta)$, that allows us to express the dynamical exponent $z(eta_m)$ of this process in terms of its persistence exponent $theta(eta_m)$. Our measurements of $z(eta_m)$ and $theta(eta_m)$ are consistent with this relation over a range of values of the universal equilibrium exponent $eta_m$, and yield subdiffusive ($z>2$) values of $z$ in the entire range. Thus we demonstrate that the worms represent a discrete-time realization of a fractional Brownian motion characterized by these properties.
سجل دخول لتتمكن من نشر تعليقات
التعليقات
جاري جلب التعليقات
سجل دخول لتتمكن من متابعة معايير البحث التي قمت باختيارها
