Critical behavior of a cubic-lattice 3D Ising model for systems with quenched disorder

A Monte Carlo method is applied to simulate the static critical behavior of a cubic-lattice 3D Ising model for systems with quenched disorder. Numerical results are presented for the spin concentrations of p = 1.0, 0.95, 0.9, 0.8, 0.6 on L × L × L lattices with L = 20–60 under periodic boundary conditions. The critical temperature is determined by the Binder cumulant method. A finite-size scaling technique is used to calculate the static critical exponents α, β, γ, and ν (for specific heat, susceptibility, magnetization, and correlation length, respectively) in the range of p under study. Universality classes of critical behavior are discussed for three-dimensional diluted systems.     

Investigation of the critical behavior of the three-dimensional weakly diluted Potts model

The static critical behavior of the three-dimensional weakly diluted Potts model with the state q = 3 on a simple cubic lattice has been investigated by the Monte Carlo method using the Wolff single-cluster algorithm. It is shown that at the spin concentrations p = 0.9 and 0.8 a second-order phase transition is observed in the three-dimensional weakly diluted Potts model with the state q = 3. On the basis of the finite-size scaling theory, we calculated the static critical exponents of the specific heat α, susceptibility γ, magnetization β, and the correlation-length exponent v.     

Static critical behavior of 3D frustrated Heisenberg model on stacked triangular lattice with variable interlayer exchange coupling

Several Monte Carlo algorithms are used to examine the critical behavior of the 3D frustrated Heisenberg model on stacked triangular lattice with variable interlayer exchange coupling for values of the interlayer-to-intralayer exchange ratio R = |′/J| in the interval between 0.01 and 1.0. A finite-size scaling technique is used to calculate the static magnetic and chiral critical exponents α (specific heat), γ and γ_k (susceptibility), β and β_k(magnetization), ν and ν_k(correlation length), and the Fisher exponent η. It is shown that 3D frustrated Heisenberg models on stacked triangular lattice with R > 0.05 constitute a new universality class of critical behavior. At lower R, a crossover from 3D to 2D critical behavior is observed.     

Investigation of the critical properties of the three-dimensional weakly diluted potts model

The problem of the type of the phase transition in the three-dimensional weakly diluted Potts model with the number of spin states q= 3 has been investigated by the Monte Carlo method. The temperature dependences of the Binder cumulants, energy, magnetization, specific heat, and susceptibility have been calculated. It is found that the second-order phase transition occurs in a system at the spin concentration p = 0.9. The critical exponents of the magnetization (β), specific heat (α), and susceptibility (γ) and the critical correlation-length exponent v were calculated on the basis of the finite-size scaling theory at p = 0.9.     

Effects of Superaging and Percolation Crossover on the Nonequilibrium Critical Behavior of the Two-Dimensional Disordered Ising Model

A Monte Carlo study of the specific features of the nonequilibrium critical behavior has been performed for the two-dimensional “pure” and structurally disordered Ising models in the course of their evolution from the low-temperature initial state at spin concentrations p = 1.0, 0.9, and 0.8. It is shown for the first time that the pinning of domain walls by structural defects leads to the anomalously strong slowing down in the evolution of the autocorrelation function characterized by the superaging effect with exponents μ = 6.25(5) and μ = 6.75(5) for the model with the spin concentrations p = 0.9 and 0.8, respectively. The pure model exhibits the conventional aging with the exponent μ = 1. It is found that the superaging effects in structurally disordered systems lead to vanishing of the limiting fluctuation?dissipation ratio X^∞, whereas X^∞ = 0.751(24) for the pure model.     

A comparison on absorption coefficients for secondary electron emission obtained from two different formulas

In this work, the absorption coefficients for secondary electron emißsion, α and β, that appeared respectively in the two different formulas, \(\delta (E_p ) = k\int_0^\infty {\left( {\frac{{dE}}{{dz}}} \right)E_p \exp ( - \alpha z)dz} \) and \(\delta (E_p ) = k\int_0^\infty {\left( {\frac{{dE}}{{dz}}} \right)E_p \exp ( - \alpha z)dz} \), were derived with a standard deviation rate analysis method based on a Monte Carlo simulated secondary electron yield, δ(E_p). Both the energy dissipation in depth for primary electrons, \(\left( {dE/dz} \right)E_p \), and the depth distribution for the number of secondary electrons including cascade electrons, n(z, E_p), were obtained by the same Monte Carlo method, in which the discrete inelastic scattering model was employed. The calculation results for Cu and Mg show that the n(z, E_p)-curve differs significantly from the \(\left( {dE/dz} \right)E_p \)-curve, and thus as well as a from b, for varied incidence angles (0°–80°) and low-energy primary electrons (up to 3 keV). The absorption coefficient β-values derived from a realistic depth distribution of cascade secondary electrons, n(z, E_p), then describe more accurately the nature of attenuation behavior of secondary electrons than a-values that associated with the approximate formula.     

Critical behavior of a stochastic anisotropic Bak–Sneppen model

In this paper we present our study on the critical behavior of a stochastic anisotropic Bak–Sneppen (saBS) model, in which a parameter α is introduced to describe the interaction strength among nearest species. We estimate the threshold fitness f _c and the critical exponent τ _r by numerically integrating a master equation for the distribution of avalanche spatial sizes. Other critical exponents are then evaluated from previously known scaling relations. The numerical results are in good agreement with the counterparts yielded by the Monte Carlo simulations. Our results indicate that all saBS models with nonzero interaction strength exhibit self-organized criticality, and fall into the same universality class, by sharing the universal critical exponents.     

Scaling exponents of forced polymer translocation through a nanopore

We investigate several properties of a translocating homopolymer through a thin pore driven by an external field present inside the pore only using Langevin Dynamics (LD) simulations in three dimensions (3D). Motivated by several recent theoretical and numerical studies that are apparently at odds with each other, we estimate the exponents describing the scaling with chain length (Nof the average translocation time \(\ensuremath \langle\tau\rangle\) , the average velocity of the center of mass \(\ensuremath \langle v_{{\rm CM}}\rangle\) , and the effective radius of gyration \(\ensuremath \langle {R}_g\rangle\) during the translocation process defined as \(\ensuremath \langle\tau\rangle \sim N^{\alpha}\) , \(\ensuremath \langle v_{{\rm CM}} \rangle \sim N^{-\delta}\) , and \(\ensuremath {R}_g \sim N^{\bar{\nu}}\) respectively, and the exponent of the translocation coordinate (s -coordinate) as a function of the translocation time \(\ensuremath \langle s^2(t)\rangle\sim t^{\beta}\) . We find \(\ensuremath \alpha=1.36 \pm 0.01\) , \(\ensuremath \beta=1.60 \pm 0.01\) for \(\ensuremath \langle s^2(t)\rangle\sim \tau^{\beta}\) and \(\ensuremath \bar{\beta}=1.44 \pm 0.02\) for \(\ensuremath \langle\Delta s^2(t)\rangle\sim\tau^{\bar{\beta}}\) , \(\ensuremath \delta=0.81 \pm 0.04\) , and \(\ensuremath \bar{\nu}\simeq\nu=0.59 \pm 0.01\) , where \( \nu\) is the equilibrium Flory exponent in 3D. Therefore, we find that \(\ensuremath \langle\tau\rangle\sim N^{1.36}\) is consistent with the estimate of \(\ensuremath \langle\tau\rangle\sim\langle R_g \rangle/\langle v_{{\rm CM}} \rangle\) . However, as observed previously in Monte Carlo (MC) calculations by Kantor and Kardar (Y. Kantor, M. Kardar, Phys. Rev. E 69, 021806 (2004)) we also find the exponent α = 1.36 ± 0.01 < 1 + ν. Further, we find that the parallel and perpendicular components of the gyration radii, where one considers the “cis” and “trans” parts of the chain separately, exhibit distinct out-of-equilibrium effects. We also discuss the dependence of the effective exponents on the pore geometry for the range of N studied here.     

Universality in the partially anisotropic three-dimensional Ising lattice

Using transfer-matrix extended phenomenological renormalization-group methods, we study the critical properties of the spin-1/2 Ising model on a simple-cubic lattice with partly anisotropic coupling strengths \(\mathop J\limits^ \to = (J',J',J)\). The universality of both fundamental critical exponents y _t and y _h is confirmed. It is shown that the critical finite-size scaling amplitude ratios \(U = A_{\chi ^{(4)} } A_\kappa /A_\chi ^2 ,Y_1 = A_{\kappa ''} /A_\chi\), and \(Y_2 = A_{\kappa ^{(4)} } /A_{\chi ^{(4)} }\) are independent of the lattice anisotropy parameter Δ=J′/J. For the Y² invariant of the three-dimensional Ising universality class, we give the first quantitative estimate Y²≈2.013 (shape L×L×∞, periodic boundary conditions in both transverse directions).     

Suche nach Strahlungsprozessen zweiter Ordnung beim Zerfall von Ag109m

The decay of an excited state by the emission of twoγ-quanta (γ γ-transitions) or two conversion electrons (e e-transitions) or oneγ-quantum and one conversion electron (γ e-transitions) is expected as a second order radiation process. The decay of Ag^109m was examined for such events using a special arrangement of two NaJ-scintillation counters in coincidence. The energies of coincident quanta were displayed on the two axes of an “X-Y”-Oscilloscope respectively. For the ratio ofγ γ-transitions to one-quantum transitions an upper limit of\(\frac{{W_{\gamma \gamma } }}{{W_\gamma }} \leqq 1,9 \cdot 10^{ - 5} \) was obtained. Furthermore theγ-spectrum in coincidence withK X-rays was studied. From these measurementse e- andγ e-transition rates can be calculated for the case ofK shell conversion. The results obtained are:

$$\frac{{W_{^e K^e K} }}{{W_\gamma }} = \left( {8,1_{ - 1,7}^{ + 0,6} } \right) \cdot 10^{ - 3} and\frac{{W_{\gamma ^e K} }}{{W_\gamma }}< 1,5 \cdot 10^{ - 3} .$$     

The convergence behaviour of expansions in the classical collision theory

In the classical collision theory the scattering angle? depends on the impact parameterb and on the kinetic energyE _r of the relative motion. This angle?(b, E _r ) is expanded for two limiting cases: 1. Expansion in powers of the potentialV(r)/E _r (momentum approximation). 2. Expansion in powers of the impact parameterb (central collision approximation). The radius of convergence of the series depends onb andE _r . It will be given for the following potentialsV(r):

$$A\left( {\frac{a}{r}} \right)^\mu ;Ae^{ - \frac{r}{a}} ;A\frac{a}{r}e^{ - \frac{r}{a}} ;A\left( {\frac{a}{r}} \right)^2 e^{ - \left( {\frac{r}{a}} \right)^2 } .$$     

Die relative Isotopenlage des Os186

Using three different samples of Os¹⁹², Os¹⁹⁰ and Os¹⁸⁶ enriched to 98.68%, 95.46%, and 61.27%, respectively, the isotopic shiftsΔν(Os¹⁹²?Os¹⁹⁰) andΔν(Os¹⁹²?Os¹⁸⁶) have been measured in four Os I-lines with a recording Fabry-Perot spectrometer. The relative isotopic position of Os¹⁸⁶ has been found to be

$$\frac{{\Delta \nu (Os^{192} - Os^{186} )}}{{\Delta \nu (Os^{192} - Os^{190} )}} = 3.52 \pm 0.03.$$     

Gleichzeitige Messung von Hyperfeinstruktur,Starkeffekt und Zeemaneffekt des23Na19F mit einer Molekülstrahl-Resonanzapparatur

An electric molecular beam resonance spectrometer has been used to measure simultaneously the Zeeman- and Stark-effect splitting of the hyperfine structure of²³Na¹⁹F. Electric four pole lenses served as focusing and refocusing fields of the spectrometer. A homogenous magnetic field (Zeeman field) was superimposed to the electric field (Stark field) in the transition region of the apparatus. The observed (Δm _J=±1)-transitions were induced electrically. Completely resolved spectra of NaF in theJ=1 rotational state have been measured in several vibrational states. The obtained quantities are: The electric dipolmomentμ _el of the molecule forv=0, 1 and 2, the rotational magnetic dipolmomentμ _J forv=0, 1, the difference of the magnetic shielding (σ _⊥-σ _‖) by the electrons of both nuclei as well as the difference of the molecular susceptibility (ξ _⊥-ξ _‖), the spin rotational constantsc _F andc _Na, the scalar and the tensor part of the molecular spin-spin interaction, the quadrupol interactione q Q forv=0, 1 and 2. The numerical values are

$$\begin{gathered} \mu _{\mathfrak{e}1} = 8,152(6) deb \hfill \\ \frac{{\mu _{\mathfrak{e}1} (v = 1)}}{{\mu _{\mathfrak{e}1} (v = 0)}} = 1,007985 (7) \hfill \\ \frac{{\mu _{\mathfrak{e}1} (v = 2)}}{{\mu _{\mathfrak{e}1} (v = 1)}} = 1,00798 (5) \hfill \\ \mu _J = - 2,89(3)10^{ - 6} \mu _B \hfill \\ \frac{{\mu _J (v = 0)}}{{\mu _J (v = 1)}} = 1,020 (13) \hfill \\ (\sigma _ \bot - \sigma _\parallel )_{Na} = - 51(12) \cdot 10^{ - 5} \hfill \\ (\sigma _ \bot - \sigma _\parallel )_F = - 51(12) \cdot 10^{ - 6} \hfill \\ (\xi _ \bot - \xi _\parallel ) = - 1,59(120)10^{ - 30} erg/Gau\beta ^2 \hfill \\ {}^CNa/^h = 1,7 (2)kHz \hfill \\ {}^CF/^h = 2,2 (2)kHz \hfill \\ {}^dT/^h = 3,7 (2)kHz \hfill \\ {}^dS/^h = 0,2 (2)kHz \hfill \\ eq Q/h = - 8,4393 (19)MHz \hfill \\ \frac{{eq Q(v = 0)}}{{eq Q(v = 1)}} = 1,0134 (2) \hfill \\ \frac{{eq Q(v = 1)}}{{eq Q(v = 2)}} = 1,0135 (2) \hfill \\ \end{gathered} $$     

Two Phases of the Non-Commutative Quantum Mechanics with the Generalized Uncertainty Relations

We consider the quantum mechanics on the noncommutative plane with the generalized uncertainty relations \({\Delta } x_{1} {\Delta } x_{2} \ge \frac {\theta }{2}, {\Delta } p_{1} {\Delta } p_{2} \ge \frac {\bar {\theta }}{2}, {\Delta } x_{i} {\Delta } p_{i} \ge \frac {\hbar }{2}, {\Delta } x_{1} {\Delta } p_{2} \ge \frac {\eta }{2}\). We show that the model has two essentially different phases which is determined by \(\kappa = 1 + \frac {1}{\hbar ^{2} } (\eta ^{2} - \theta \bar {\theta })\). We construct a operator \(\hat {\pi }_{i}\) commuting with \(\hat {x}_{j} \) and discuss the harmonic oscillator model in two dimensional non-commutative space for three case κ > 0, κ = 0, κ < 0. Finally, we discuss the thermodynamics of a particle whose hamiltonian is related to the harmonic oscillator model in two dimensional non-commutative space.     

Entanglement Involved in Time Evolution of Two-Mode Squeezed State in Single-Mode Diffusion Channel

We derive the evolution law of an initial two-mode squeezed vacuum state \( \text {sech}^{2}\lambda e^{a^{\dag }b^{\dagger }\tanh \lambda }\left \vert 00\right \rangle \left \langle 00\right \vert e^{ab\tanh \lambda }\) (a pure state) passing through an a-mode diffusion channel described by the master equation

$$\frac{d\rho \left( t\right) }{dt}=-\kappa \left[ a^{\dagger}a\rho \left( t\right) -a^{\dagger}\rho \left( t\right) a-a\rho \left( t\right) a^{\dagger}+\rho \left( t\right) aa^{\dagger}\right] , $$

since the two-mode squeezed state is simultaneously an entangled state, the final state which emerges from this channel is a two-mode mixed state. Performing partial trace over the b-mode of ρ(t) yields a new chaotic field, \(\rho _{a}\left (t\right ) =\frac {\text {sech}^{2}\lambda }{1+\kappa t \text {sech}^{2}\lambda }:\exp \left [ \frac {- \text {sech}^{2}\lambda }{1+\kappa t\text {sech}^{2}\lambda }a^{\dagger }a \right ] :,\) which exhibits higher temperature and more photon numbers, showing the diffusion effect. Besides, measuring a-mode of ρ(t) to find n photons will result in the collapse of the two-mode system to a new Laguerre polynomial-weighted chaotic state in b-mode, which also exhibits entanglement.

     

Parabolic Anderson Model in a Dynamic Random Environment: Random Conductances

The parabolic Anderson model is defined as the partial differential equation ? u(x, t)/? t = κ Δ u(x, t) + ξ(x, t)u(x, t), x ∈ ? ^d , t ≥ 0, where κ ∈ [0, ∞) is the diffusion constant, Δ is the discrete Laplacian, and ξ is a dynamic random environment that drives the equation. The initial condition u(x, 0) = u ₀(x), x ∈ ? ^d , is typically taken to be non-negative and bounded. The solution of the parabolic Anderson equation describes the evolution of a field of particles performing independent simple random walks with binary branching: particles jump at rate 2d κ, split into two at rate ξ ∨ 0, and die at rate (?ξ) ∨ 0. In earlier work we looked at the Lyapunov exponents

$$ \lambda _{p}(\kappa ) = \lim\limits _{t\to \infty } \frac {1}{t} \log \mathbb {E} ([u(0,t)]^{p})^{1/p}, \quad p \in \mathbb{N} , \qquad \lambda _{0}(\kappa ) = \lim\limits _{t\to \infty } \frac {1}{t}\log u(0,t). $$

For the former we derived quantitative results on the κ-dependence for four choices of ξ : space-time white noise, independent simple random walks, the exclusion process and the voter model. For the latter we obtained qualitative results under certain space-time mixing conditions on ξ. In the present paper we investigate what happens when κΔ is replaced by Δ^??, where ?? = {??(x, y) : x, y ∈ ? _d , x ～ y} is a collection of random conductances between neighbouring sites replacing the constant conductances κ in the homogeneous model. We show that the associated annealed Lyapunov exponents λ _p (??), p ∈ ?, are given by the formula

$$ \lambda _{p}(\mathcal{K} ) = \text{sup} \{\lambda _{p}(\kappa ) : \, \kappa \in \text{Supp} (\mathcal{K} )\}, $$

where, for a fixed realisation of ??, Supp(??) is the set of values taken by the ??-field. We also show that for the associated quenched Lyapunov exponent λ ₀(??) this formula only provides a lower bound, and we conjecture that an upper bound holds when Supp(??) is replaced by its convex hull. Our proof is valid for three classes of reversible ξ, and for all ?? satisfying a certain clustering property, namely, there are arbitrarily large balls where ?? is almost constant and close to any value in Supp(??). What our result says is that the annealed Lyapunov exponents are controlled by those pockets of ?? where the conductances are close to the value that maximises the growth in the homogeneous setting. In contrast our conjecture says that the quenched Lyapunov exponent is controlled by a mixture of pockets of ?? where the conductances are nearly constant. Our proof is based on variational representations and confinement arguments.

     

Thermodynamics of the frustrated J<Subscript>1</Subscript>-J<Subscript>2</Subscript> Heisenberg ferromagnet on the body-centered cubic lattice with arbitrary spin

We use the spin-rotation-invariant Green’s function method as well as thehigh-temperature expansion to discuss the thermodynamic properties of the frustratedspin-S J ₁-J ₂ Heisenbergmagnet on the body-centered cubic lattice. We consider ferromagnetic nearest-neighborbonds J ₁<0 and antiferromagnetic next-nearest-neighbor bonds J ₂ ≥ 0 andarbitrary spin S. We find that the transition point\hbox{$J_2^c$}J2cbetween the ferromagnetic ground state and theantiferromagnetic one is nearly independent of the spin S, i.e., it is very closeto the classical transition point\hbox{$J_2^{c,{\rm clas}}= \frac{2}{3}|J_1|$}J2c,clas=23|J1|. At finite temperatures we focus on the parameterregime\hbox{$J_2<J_2^c$}J2<J2cwith a ferromagnetic ground-state. We calculate theCurie temperature T _C (S, J ₂)and derive an empirical formula describing the influence of the frustration parameterJ ₂ and spin S on T _C . We find that theCurie temperature monotonically decreases with increasing frustration J ₂, where veryclose to\hbox{$J_2^{c,{\rm clas}}$}J2c,clasthe T _C (J ₂)-curveexhibits a fast decay which is well described by a logarithmic term\hbox{$1/\textrm{log}(\frac{2}{3}|J_1|-J_{2})$}1/log(23|J1|?J2). To characterize the magnetic ordering below and aboveT _C , we calculate thespin-spin correlation functions ?S ₀ S _R ?, the spontaneous magnetization, the uniform static susceptibilityχ ₀ as well as the correlation lengthξ.Moreover, we discuss the specific heat C _V and the temperaturedependence of the excitation spectrum. As approaching the transition point\hbox{$J_2^c$}J2csome unusual features were found, such as negativespin-spin correlations at temperatures above T _C even though theground state is ferromagnetic or an increase of the spin stiffness with growingtemperature.     

Rules of correspondence in atomic physics

The Smirnov method of analytic continuation (B.M. Smirnov, Sov. Phys. JETP 20, 345 (1964)) has been justified and developed for atomic physics. It has been shown that the polarizability of alkali atoms α, their van der Waals interaction constant C ₆, and the oscillator strength of the transition to the first P state f ₀₁ are related to the parameter 〈r ²〉 and gap in the spectrum \(\frac{3}{2}\frac{f}{\Delta } \approx \frac{3}{2}\alpha \Delta \approx {\left( {3{C_6}\Delta } \right)^{{\raise0.7ex\hbox{$1$} \!\mathord{\left/ {\vphantom {1 2}}\right.\kern-\nulldelimiterspace}\!\lower0.7ex\hbox{$2$}}}} \approx \left\langle {{r^2}} \right\rangle \). The average square of the coordinate of the valence electron 〈r ²〉 in the first approximation has a hydrogen dependence \({J_1} = \frac{1}{{2{v^2}}}.\) on the filling factor ν, which is defined in terms of the first ionization potential: xxxxxxxxx     

Dynamic scaling behaviors of the restricted-solid-on-solid model on honeycomb and square-octagon lattice substrates

The dynamic scaling behaviors of the restricted-solid-on-solid (RSOS) model on two new types of substrate, which are honeycomb and square-octagon lattice substrates, are studied by means of Kinetic Monte Carlo simulations. The growth exponent β and the roughness exponent α defined, respectively, by the surface width via W ~ t ^β and the saturated width via W _sat ~ L ^α , L being the system size, were obtained by a power-counting analysis. Our simulation results show that the Family-Vicsek scaling is still satisfied. However, the structures of the substrates indeed affect the dynamic behavior of the growth model. The values of the roughness exponents fall between regular and fractal lattices. Deeper analysis show that the coordination number of the substrates play an crucial role.     

Spherical Model on a Cayley Tree: Large Deviations

We study the spherical model of a ferromagnet on a Cayley tree and show that in the case of empty boundary conditions a ferromagnetic phase transition takes place at the critical temperature \(T_\mathrm{c} =\frac{6\sqrt{2}}{5}J\), where J is the interaction strength. For any temperature the equilibrium magnetization, \(m_n\), tends to zero in the thermodynamic limit, and the true order parameter is the renormalized magnetization \(r_n=n^{3/2}m_n\), where n is the number of generations in the Cayley tree. Below \(T_\mathrm{c}\), the equilibrium values of the order parameter are given by \(\pm \rho ^*\), where

$$\begin{aligned} \rho ^*=\frac{2\pi }{(\sqrt{2}-1)^2}\sqrt{1-\frac{T}{T_\mathrm{c}}}. \end{aligned}$$

One more notable temperature in the model is the penetration temperature

$$\begin{aligned} T_\mathrm{p}=\frac{J}{W_\mathrm{Cayley}(3/2)}\left( 1-\frac{1}{\sqrt{2}}\left( \frac{h}{2J}\right) ^2\right) . \end{aligned}$$

Below \(T_\mathrm{p}\) the influence of homogeneous boundary field of magnitude h penetrates throughout the tree. The main new technical result of the paper is a complete set of orthonormal eigenvectors for the discrete Laplace operator on a Cayley tree.