Monte Carlo determination of the low-energy constants for a two-dimensional spin-1 Heisenberg model with spatial anisotropy

The low-energy constants, namely the spin stiffness ρ _s , the staggered magnetization density ? _s per area, and the spinwave velocity c of the two-dimensional (2D) spin-1 Heisenberg model on the square and rectangular lattices are determined using the first principles Monte Carlo method. In particular, the studied models have different antiferromagnetic couplings J ₁ and J ₂ in the spatial 1- and 2-directions, respectively. For each considered J ₂∕J ₁, the aspect ratio of the corresponding linear box sizes L ₂∕L ₁ used in the simulations is adjusted so that the squares of the two spatial winding numbers take the same values. In addition, the relevant finite-volume and -temperature predictions from magnon chiral perturbation theory are employed in extracting the numerical values of these low-energy constants. Our results of ρ _s1 are in quantitative agreement with those obtained by the series expansion method over a broad range of J ₂∕J ₁. This in turn provides convincing numerical evidence for the quantitative correctness of our approach. The ? _s and c presented here for the spatially anisotropic models are new and can be used as benchmarks for future related studies.     

Anisotropic<Emphasis Type="Italic">s</Emphasis>-wave or<Emphasis Type="Italic">d</Emphasis>-wave pairing? Analysis of experiments on ion irradiation of YBa<Subscript>2</Subscript>Cu<Subscript>3</Subscript>O<Subscript>7-x</Subscript> films

The effect of ion irradiation on the superconducting transition temperatureT _c and resistivityρ _ab (T) of YBa₂Cu₃O_7-x films with different oxygen content (initial temperatureT _c0≈90 K and 60 K) is studied experimentally. The dependenciesT _c /T _c0 on residual resistivityρ _o are obtained in very wide range 0.2<T _c /T _c0 <1 andρ _o μΩ·cm. The critical values ofρ _o , corresponding to the vanishing of superconductivity, are found to be an order of magnitude larger then those predicted by theory ford-wave pairing. At 0.5÷0.6<T _c /T _c0<1 the experimental data are in close agreement with theoretical dependencies, obtained for the anisotropics-wave superconductor within the BCS-framework.     

Dependence of structure factor and correlation energy on the width of electron wires

The structure factor and correlation energy of a quantum wire of thickness b ? a _B are studied in random phase approximation (RPA) and for the less investigated region r _s < 1. Using the single-loop approximation, analytical expressions of the structure factor are obtained. The exact expressions for the exchange energy are also derived for a cylindrical and harmonic wire. The correlation energy in RPA is found to be represented by ? _c (b, r _s ) = α(r _s )/b + β(r _s ) ln(b) + η(r _s ), for small b and high densities. For a pragmatic width of the wire, the correlation energy is in agreement with the quantum Monte Carlo simulation data.     

Role of critical fluctuations in the formation of a skyrmion lattice in MnSi

The region in the H–T phase diagram near the critical temperature (T _c ) of the cubic helicoidal MnSi magnet is comprehensively studied by small-angle neutron diffraction. Magnetic field H is applied along the [111] axis. The experimental geometry is chosen to simultaneously observe the following three different magnetic states of the system: (a) critical fluctuations of a spin spiral with randomly orientated wavevector k _f , (b) conical structure with k _c ∥ H, and (c) hexagonal skyrmion lattice with k_sk ⊥ H. Both states (conical structure, and skyrmion lattice) are shown to exist above critical temperature T _c = 29 K against the background of the critical fluctuations of a spin spiral. The conical lattice is present up to the temperatures where fluctuation correlation length ξ becomes comparable with pitch of spiral d _s . The skyrmion lattice is localized near T _c and is related to the fluctuations of a spiral with correlation length ξ ≈ 2d _s , and the propagation vector is normal to the field (k_sk ⊥ H). These spiral fluctuations are assumed to be the defects that stabilize the skyrmion lattice and promote its formation.     

Topologic distance in the Lucena network

The Lucena network (LN) is the dual of a multifractal partition of the square. We analyzethe relation between the typical topologic distance l and the number ofvertices Nof the LN. The multifractal partition has one parameter ρ which controls thegeometrical asymmetry of the multifractal. In the limit of ρ → 1 the blocks of thepartition are squared, the connections amont the blocks are short range, the LN is moreregular and the relation l ∝ √N is observed. For the limit ρ → 0 the blocks arestrongly asymmetric, long range connections appear, and the topologic distance followsl ∝(log?N)^α, a weak smallworld phenomenon. For any network size we calculate analytically the size of the minimumdistance l_min (ρ → 0) and the maximaldistance l_max (ρ → 1). The distance in theweak small world regime is calculated using the number of vertices inside a radius oflength land taking into account the network average connectivity and the exponent α.     

Probing the sign-changeable interaction between dark energy and dark matter with current observations

We consider the models of vacuum energy interacting with cold dark matter in this study, in which the coupling can change sigh during the cosmological evolution. We parameterize the running coupling b by the form b(a) = b_0 a + b_e(1-a), where at the earlytime the coupling is given by a constant b_e and today the coupling is described by another constant b_0. We explore six specific models with(i) Q = b(a)H_0ρ_0,(ii) Q = b(a)H_0ρ_(de),(iii) Q = b(a)H_0ρ_c,(iv) Q = b(a)Hρ_0,(v) Q = b(a)Hρ_(de), and(vi) Q = b(a)Hρ_c.The current observational data sets we use to constrain the models include the JLA compilation of type Ia supernova data, the Planck 2015 distance priors data of cosmic microwave background observation, the baryon acoustic oscillations measurements,and the Hubble constant direct measurement. We find that, for all the models, we have b_0 0 and b_e 0 at around the 1σ level,and b_0 and b_e are in extremely strong anti-correlation. Our results show that the coupling changes sign during the evolution at about the 1σ level, i.e., the energy transfer is from dark matter to dark energy when dark matter dominates the universe and the energy transfer is from dark energy to dark matter when dark energy dominates the universe.     

A Monte Carlo study of the first-order transition in a Heisenberg FCC antiferromagnet

A nearest neighbor Heisenberg antiferromagnet on a face-centered cubic lattice is studied by extensive Monte Carlo simulations in zero magnetic field. The parallel tempering algorithm is utilized, which allows one to overcome a slow relaxation of the magnetic order parameter and to fully equilibrate moderately sized clusters with up to N ? 7 × 10³ spins. By collecting energy and order parameter histograms on clusters with up to N ? 2 × 10⁴ sites, we accurately locate the first-order transition point at T _c = 0.4459(1)J.     

Spin-precession-induced currents through ferromagnet/nanomagnet/superconductor junctions

The spin-precession-induced current through ferromagnet/nanomagnet/superconductorjunctions is investigated by using the nonequilibrium Green’s function method. It is foundthat the charge current I _c for the spinprecession frequency ω less than the energy gap Δ onlyarises from the equal-spin Andreev reflection, which is independent of the spinpolarization p of the ferromagnetic lead, while that forω > Δ mainly originates from the quasiparticle’scontribution. Both equal-spin AR and quasiparticle scattering processes contribute to thespin current I _s and the quasiparticlescattering process plays a dominant role. WhileI _c forω < Δ can be enhanced by the spin polarizationp, I _s decreases withp. These features may be of interest for ongoing experiments in thefield of molecular spintronics.     

Duality-mediated critical amplitude ratios for the (2 + 1)-dimensional S = 1XY model

The phase transition for the (2 + 1)-dimensional spin-S = 1XY model was investigated numerically. Because of the boson-vortex duality, the spin stiffness ρ _s in the ordered phase and the vortex-condensate stiffness ρ _v in the disordered phase should have a close relationship. We employed the exact diagonalization method, which yields the excitation gap directly. As a result, we estimate the amplitude ratios ρ _s,v /Δ (Δ: Mott insulator gap) by means of the scaling analyses for the finite-size cluster with N ≤ 22 spins. The ratio ρ _s /ρ _v admits a quantitative measure of deviation from selfduality.     

Ferrimagnetic nanoparticles for self-controlled magnetic hyperthermia

Based on the Heisenberg model including single-site uniaxial anisotropy and using aGreen’s function technique we studied the influence of size and composition effects on theCurie temperature T _C , saturationmagnetization M _S and coercivityH _C of spherical nanoparticles with astructural formulaM e _1?x Zn _x Fe₂O₄,Me = Ni, Cu, Co, Mn. It is shown that for x = 0.4–0.5and d = 10–20 nm these nanoparticles have aT _C  = 315 K and are suitable for aself-controlled magnetic hyperthermia.     

Quantifying Correlations via the Wigner-Yanase-Dyson Skew Information

In this paper, based on a discussion about the Wigner-Yanase-Dyson (WYD) skew information, the measure F_a,α(ρ_ab) for correlations in terms of the WYD skew information is introduced and discussed. The following conclusions are obtained. For a classical-quantum state ρ_ab, F_a,α(ρ_ab)=0 if and only if ρ_ab is a product state; F_a,α(ρ_ab) is locally unitary invariant and convex on the set of states with the fixed marginal ρ_a; F_a,α(ρ_ab) decreases under local random unitary operation on H_b; For a quantum-classical state ρ_ab, F_a,α(ρ_ab) decreases under local operation on H_b; Lastly, F_a,α(ρ_ab) is computed for the pure states and the Bell-diagonal states, respectively.     

Long-range spin correlations in a honeycomb spin model with a magnetic field

We consider the spin-1/2 model on the honeycomb lattice [A. Kitaev, Ann. Phys. 321, 2 (2006)] in the presence of a weak magnetic field h _α ? J. Such a perturbation treated in the lowest nonvanishing order over h _α leads [K.S. Tikhonov, M.V. Feigel’man, and A.Yu. Kitaev, Phys. Rev. Lett. 106, 067203 (2011)] to a powerlaw decay of irreducible spin correlations 《s ^z (t, r)s ^z (0, 0)》 ∝ h _z ² f(t, r), where f(t, r) ∝ [max(t, Jr)]^–4. We have studied the effects of the next order of perturbation in h _z and found an additional term of the order h _z ⁴ in the correlation function 《s ^z (t, r)s ^z (0, 0)》 which scales as h _z ⁴ cosγ/r ³ at Jt? r, where γ is the polar angle in the 2D plane. We demonstrate that such a contribution can be understood as a result of a perturbation of the effective Majorana Hamiltonian by the weak imaginary vector potential A _x ∝ i h _z ² .     

<Emphasis Type="Italic">K</Emphasis><Superscript>±</Superscript><Emphasis Type="Italic">p</Emphasis> elastic scattering in the QCD inspired eikonalized model

Based on our previous study of the QCD inspired eikonalized model for describing vector meson photoproduction, pp, and \(\bar p\) p elastic scattering at high energies, we apply the mode to high energy K ^± p elastic scattering. The total cross section σ _tot(s), differential cross section dσ/dt, the ratio of the real part to imaginary part of the forward scattering amplitude ρ(s), and nuclear slope parameter function β(s) are calculated in the model. Our results show that the theoretical prediction for σ _tot(s) is in a good agreement with the experimental data within error bars of the data. For the other theoretical predictions there are no data to test the predictive power of the model. We need the corresponding experimental data to examinate the validity of our QCD inspired eikonalized model. However, our calculations clearly show that the Odderon exchange in the process makes a significant contribution to the observable of ρ(s) and β(s). Therefore, we may conclude that there is a good opportunity to find the QCD Odderon in the K ^± p elastic scattering at high energies.     

Schur Polynomials and The Yang-Baxter Equation

We describe a parametrized Yang-Baxter equation with nonabelian parameter group. That is, we show that there is an injective map \({g \mapsto R (g)}\) from \({ \rm{GL}(2, \mathbb{C}) \times \rm{GL}(1, \mathbb{C})}\) to End \({(V \otimes V)}\) , where V is a two-dimensional vector space such that if \({g, h \in G}\) then R ₁₂(g)R ₁₃(gh) R ₂₃(h) = R ₂₃(h) R ₁₃(gh)R ₁₂(g). Here R _{i j} denotes R applied to the i, j components of \({V \otimes V \otimes V}\) . The image of this map consists of matrices whose nonzero coefficients a ₁, a ₂, b ₁, b ₂, c ₁, c ₂ are the Boltzmann weights for the non-field-free six-vertex model, constrained to satisfy a ₁ a ₂ + b ₁ b ₂ ? c ₁ c ₂ = 0. This is the exact center of the disordered regime, and is contained within the free fermionic eight-vertex models of Fan and Wu. As an application, we show that with boundary conditions corresponding to integer partitions λ, the six-vertex model is exactly solvable and equal to a Schur polynomial s _λ times a deformation of the Weyl denominator. This generalizes and gives a new proof of results of Tokuyama and Hamel and King.     

Revisiting the interacting model of new agegraphic dark energy

In this paper,a new version of the interacting model of new agegraphic dark energy(INADE)is proposed and analyzed in detail.The interaction between dark energy and dark matter is reconsidered.The interaction term Q=bH0ραdeρ1αdm is adopted,which abandons the Hubble expansion rate H and involves bothρde andρdm.Moreover,the new initial condition for the agegraphic dark energy is used,which solves the problem of accommodating baryon matter and radiation in the model.The solution of the model can be given using an iterative algorithm.A concrete example for the calculation of the model is given.Furthermore,the model is constrained by using the combined Planck data(Planck+BAO+SNIa+H0)and the combined WMAP-9 data(WMAP+BAO+SNIa+H0).Three typical cases are considered:(A)Q=bH0ρde,(B)Q=bH0√ρdeρdm,and(C)Q=bH0ρdm,which correspond toα=1,1/2,and 0,respectively.The departures of the models from theΛCDM model are measured by the BIC and AIC values.It is shown that the INADE model is better than the NADE model in the fit,and the INADE(A)model is the best in fitting data among the three cases.     

Dependence of the magnitude and direction of the persistent current on the magnetic flux in superconducting rings

The obtained periodic magnetic-field dependences I _c+(Φ/Φ₀) and I _c?(Φ/Φ₀) of the critical current measured in opposite directions on asymmetric superconducting aluminum rings has made it possible to explain previously observed quantum oscillations of dc voltage as a result of alternating current rectification. It was found that a higher rectification efficiency of both single rings and ring systems is caused by hysteresis of the current-voltage characteristics. The asymmetry of current-voltage characteristics providing the rectification effect is due to the relative shifts of the magnetic dependences I _c?(Φ/Φ₀) = I _c+(Φ/Φ₀ + Δ?) of the critical current measured in opposite directions. This shift means that the position of I _c+(Φ/Φ₀) and I _c?(Φ/Φ₀) minima does not correspond to n + 0.5 magnetic flux Φ quanta, which is in direct contradiction to measured Little-Parks resistance oscillations. Despite this contradiction, the amplitude I _{c, an}(Φ/Φ₀) = I _c+(Φ/Φ₀) ? I _c?(Φ/Φ₀) of critical current anisotropy oscillations and its variations with temperature correspond to the expected amplitude of persistent current oscillations and its variations with temperature.     

Multistage radiation-stimulated changes in the microhardness of silicon single crystals exposed to low-intensity β irradiation

Radiation-stimulated and postradiation changes in the microhardness of silicon single crystals exposed to irradiation with a low-intensity flux of β particles (I = 9 × 10⁵ cm^?2 s^?1, W = 0.20 + 0.93 MeV) are studied. It is established that the inversion of the radiation-induced plastic effect occurs at a characteristic irradiation time τ_c = 75 min; i.e., irradiation of silicon single crystals for a time τ < τ_c leads to nonmonotonic reversible hardening, whereas nonmonotonic reversible softening is observed under irradiation for a time τ > τ_c. It is demonstrated that there exists a correlation between the nonmonotonic dependences of the microhardness and the concentration of electrically active defects at acceptor levels with energies E _c ? 0.11 eV, E _c ? 0.13 eV, and E _c ? 0.18 eV on the irradiation time.     

Understanding cooperative behavior in structurally disordered populations

The effects of an inhomogeneous competing environment on the extent of cooperation are studied within the context of a site-diluted evolutionary snowdrift game on a square lattice, with the occupied sites representing the players, both numerically and analytically. The frequency of cooperation ?_C generally shows a non-monotonic dependence on the fraction of occupied sites ρ, for different values of the payoff parameter r. Slightly diluting a lattice leads to a lower cooperation for small and high values of r. For a range of r, however, dilution leads to an enhanced cooperation. An analytic treatment is developed for?_C^I + ?_C^II, with?_C^I emphasizing the importance of the small clusters of players especially for?_C^II from the other players is shown to be inadequate. A local configuration approximation (LCA) that treats the local competing configurations as the variables and amounts to include spatial correlation up to the neighborhood of a player’s neighbors is developed. Results of ?_C (ρ) and the number of different local configurations from LCA are in good agreement with simulation results. A transparent physical picture of the dynamics stemming from LCA is also presented. The theoretical approach provides a framework that can be readily applied to competing agent-based models in structurally ordered and disordered populations.     

Structural Phase Transition of Perovskite-Type N(CH<Subscript>3</Subscript>)<Subscript>4</Subscript>CdBr<Subscript>3</Subscript> Studied by MAS NMR and Static NMR

The structural geometry change in the perovskite-type N(CH₃)₄CdBr₃ single crystal near the phase transition temperature of T _C = 390 K was investigated using magic angle spinning nuclear magnetic resonance techniques. For ¹H and ¹³C nuclei, the temperature dependences of their chemical shift, spectral intensity, and spin–lattice relaxation time (T _1ρ) in the rotating frame were obtained and analyzed. While the chemical shift and T _1ρ of ¹H showed change near T _C, those of ¹³C did not. In addition, the ¹¹³Cd spin–lattice relaxation time T ₁ in the laboratory frame near T _C show no evidence of anomalous change near the phase transition temperature, which coincides with the measured changes in the ¹H T _1ρ. The driving force for this phase transition was connected to the ¹H in the CH₃ groups.     

Role of initial correlation in coarsening of a ferromagnet

We study the dynamics of ordering in ferromagnets via Monte Carlo simulations of theIsing model, employing the Glauber spin-flip mechanism, in space dimensionsd = 2 and3, on square and simplecubic lattices. Results for the persistence probability and the domain growth arediscussed for quenches to various temperatures (T _f ) below the criticalone (T _c ), from differentinitial temperatures T _i ≥T _c . In long timelimit, for T _i >T _c ,the persistence probability exhibits power-law decay with exponents θ ? 0.22 and? 0.18 in d = 2 and 3, respectively. For finite T _i , the early timebehavior is a different power-law whose life-time diverges and exponent decreases asT _i →T _c . The two steps areconnected via power-law as a function of domain length and the crossover to the secondstep occurs when this characteristic length exceeds the equilibrium correlation length atT =T _i . T _i =T _c is expected toprovide a new universality class for which we obtain θ ≡θ _c ? 0.035 ind = 2 and?0.105 in d = 3. The time dependenceof the average domain size ?, however, is observed to be rather insensitive tothe choice of T _i .