Thermodynamic studies of the two dimensional Falicov-Kimball model on a triangular lattice

Thermodynamic properties of the spinless Falicov-Kimball model are studied on a triangular lattice using numerical diagonalization technique with Monte-Carlo simulation algorithm. Discontinuous metal-insulator transition is observed at finite temperature. Unlike the case of square lattice, here we observe that the finite temperature effect is not able to smear out the discontinuous metal-insulator transition seen in the ground state. Calculation of specific heat (C _v ) shows single and double peak structures for different values of parameters like on-site correlation strength (U), f-electron energy (E _f ) and temperature.     

Breakdown of universality in directed spiral percolation

Directed spiral percolation (DSP), percolation under both directional and rotational constraints, is studied on the triangular lattice in two dimensions (2D). The results are compared with that of the 2D square lattice. Clusters generated in this model are generally rarefied and have chiral dangling ends on both the square and triangular lattices. It is found that the clusters are more compact and less anisotropic on the triangular lattice than on the square lattice. The elongation of the clusters is in a different direction than the imposed directional constraint on both the lattices. The values of some of the critical exponents and fractal dimension are found considerably different on the two lattices. The DSP model then exhibits a breakdown of universality in 2D between the square and triangular lattices. The values of the critical exponents obtained for the triangular lattice are not only different from that of the square lattice but also different form other percolation models.Received: 12 March 2004, Published online: 23 July 2004PACS: 02.50.-r Probability theory, stochastic processes, and statistics - 64.60.-i General studies of phase transitions - 72.80.Tm Composite materials     

Dynamic phase transitions in the kinetic spin-1 Blume-Capel model on the Bethe lattice

The stationary states of the kinetic spin-1 Blume-Capel (BC) model on the Bethe lattice are analyzed in detail in terms of recursion relations. The model is described using a Glauber-type stochastic dynamics in the presence of a time-dependent oscillating external magnetic field (h) and crystal field (D) interactions. The dynamic order parameter, the hysteresis loop area and the dynamic correlation are calculated. It is found that the magnetization oscillates around nonzero values at low temperatures (T) for the ferromagnetic (F) phase while it only oscillates around zero values at high temperatures for the paramagnetic (P) phase. There are regions of the phase space where the two solutions coexist. The dynamic phase diagrams are obtained on the (kT/J,h/J) and (kT/J,D/J) planes for the coordination number q=4. In addition to second-order and first-order phase transitions, dynamical tricritical points and triple points are also observed.     

Three-dimensional array structures associated with Richtmyer-Meshkov and Rayleigh-Taylor instability

A boundary separating adjacent gas or liquid media is frequently unstable. Richtmyer-Meshkov and Rayleigh-Taylor instability cause the growth of intricate structures on such boundaries. All the lattice symmetries [rectangular (pmm2), square (p4mm), hexagonal (p6mm), and triangular (p3m1) lattices] which are of interest in connection with the instability of the surface of a fluid are studied for the first time. They are obtained from initial disturbances consisting of one (planar case, two-dimensional flow), two (rectangular cells), or three (hexagons and triangles) harmonic waves. It is shown that the dynamic system undergoes a transition during development from an initial, weakly disturbed state to a limiting or asymptotic stationary state (stationary point). The stability of these points (stationary states) is investigated. It is shown that the stationary states are stable toward large-scale disturbances both in the case of Richtmyer-Meshkov instability and in the case of Rayleigh-Taylor instability. It is discovered that the symmetry increases as the system evolves in certain cases. In one example the initial Richtmyer-Meshkov or Rayleigh-Taylor disturbance is a sum of two waves perpendicular to one another with equal wave numbers, but unequal amplitudes: a ₁(t=0)≠a ₂(t=0). Then, during evolution, the flow has p2 symmetry (rotation relative to the vertical axis by 180°), which goes over to p4 symmetry (rotation by 90°) at t→∞, since the amplitudes equalize in the stationary state: a ₁(t=∞)=a ₂(t=∞). It is shown that the hexagonal and triangular arrays are complementary. Upon time inversion (t→−t), “rephasing” occurs, and the bubbles of a hexagonal array transform into jets of a triangular array and vice versa. Zh. éksp. Teor. Fiz. 116, 908–939 (September 1999)     

Odd-frequency superconductivity on quasi-one-dimensional triangular lattice in the Hubbard model

We study the superconductivity in the Hubbard model on quasi-one-dimensional triangular lattice using random phase approximation (RPA). We find that odd-frequency spin-singlet p-wave pairing can be realized on isosceles quasi-one-dimensional triangular lattice.     

Dynamic phase transitions and the effects of frequency of oscillating magnetic field on the dynamic phase diagrams in the bilayer honeycomb lattice with AB stacking geometry

We study some dynamic properties of the bilayer honeycomb lattice with AB stacking geometry in the presence of a time-dependent oscillating external magnetic field. We employ the Glauber transition rates to construct the mean-field dynamical equations. First, we obtain dynamic phases in the system and observe the paramagnetic (p), ferromagnetic (f), compensated (c) antiferromagnetic (af), surface ferromagnetic (sf) and mixed (m) phases. Besides, coexistence phase regions also exist in the system. Second, we investigate the thermal behavior of the dynamic order parameters. From these study, the natures (first- or second-order) of the transitions are characterized and the dynamic phase transition (DPT) points are presented. The DPTs are obtained and the dynamic phase diagrams (DPD) are constructed plane of the temperature versus the amplitude of the magnetic field. We investigate the effect of the frequency of the oscillating external magnetic field on the DPD.     

Ferromagnetism in the Hubbard model: instability of the Nagaoka state on the triangular,honeycomb and kagome lattices

In order to analyse the lattice dependence of ferromagnetism in the two-dimensional Hubbard model we investigate the instability of the fully polarised ferromagnetic ground state (Nagaoka state) on the triangular, honeycomb and kagome lattices. We mainly focus on the local instability, applying single spin flip variational wave functions which include majority spin correlation effects. The question of global instability and phase separation is addressed in the framework of Hartree-Fock theory. We find a strong tendency towards Nagaoka ferromagnetism on the non-bipartite lattices (triangular, kagome) for more than half filling. For the triangular lattice we find the Nagaoka state to be unstable above a critical density of n = 1.887 at U = ∞, thereby significantly improving former variational results. For the kagome lattice the region where ferromagnetism prevails in the phase diagram widely exceeds the flat band regime. Our results even allow the stability of the Nagaoka state in a small region below half filling. In the case of the bipartite honeycomb lattice several disconnected regions are left for a possible Nagaoka ground state.     

Monte Carlo investigation of phase transitions in the frustrated Heisenberg model on a triangular lattice

The critical properties and phase transitions of the three-dimensional frustrated antiferromagnetic Heisenberg model on a triangular lattice have been investigated using the Monte Carlo method with a replica algorithm. The critical temperature has been determined and the character of the phase transitions has been analyzed using the method of fourth-order Binder cumulants. A second-order phase transition has been found in the three-dimensional frustrated Heisenberg model on a triangular lattice. The static magnetic and chiral critical exponents of the heat capacity α, the susceptibility γ and γ _k , the magnetization β and β _k , the correlation length ν and ν _k , as well as the Fisher exponents η and η _k , have been calculated in terms of the finite-size scaling theory. It has been demonstrated that the three-dimensional frustrated antiferromagnetic Heisenberg model on a triangular lattice forms a new universality class of the critical behavior.     

光子晶体波导中的孤子传输及其延迟特性研究

研究了正方晶格和三角晶格空气背景硅介质柱光子晶体线缺陷波导导模左带隙边缘处的亮孤子脉冲传播特性及其慢光延迟特性. 采用平面波展开法仿真分析了波导相邻两行介质柱大小r₁和r₂以及波导宽度D对孤子脉冲传输所需峰值功率P₀和延迟时间T_s的影响,总结了其变化规律. 通过调整波导结构得到了正方晶格和三角晶格优化波导结构,优化后,正方晶格结构波导P₀减小了81.17%,T_s增加了66.32%;三角晶格结构波导P₀减小了73.7%,T_s增加了67.63%,实现了孤子传输性能的大幅度优化. 关键词： 光子晶体波导 光孤子 峰值功率 延迟时间     

Effective-field theory on the kinetic Ising model

As an analytical method, the effective-field theory (EFT) is used to study the dynamical response of the kinetic Ising model in the presence of a sinusoidal oscillating field. The effective-field equations of motion of the average magnetization are given for the square lattice (Z=4) and the simple cubic lattice (Z=6), respectively. The dynamic order parameter, the hysteresis loop area and the dynamic correlation are calculated. In the field amplitude h₀/ZJ-temperature T/ZJ plane, the phase boundary separating the dynamic ordered and the disordered phase has been drawn, and the dynamical tricritical point has been observed. We also make the compare results of EFT with that given by using the mean field theory (MFT).     

Dynamic phase transitions and dynamic phase diagrams in the kinetic spin-5/2 Blume-Capel model in an oscillating external magnetic field: Effective-field theory and the Glauber-type stochastic dynamics approach

Using an effective field theory with correlations, we study a kinetic spin-5/2 Blume-Capel model with bilinear exchange interaction and single-ion crystal field on a square lattice. The effective-field dynamic equation is derived by employing the Glauber transition rates. First, the phases in the kinetic system are obtained by solving this dynamic equation. Then, the thermal behavior of the dynamic magnetization, the hysteresis loop area and correlation are investigated in order to characterize the nature of the dynamic transitions and to obtain dynamic phase transition temperatures. Finally, we present the phase diagrams in two planes, namely (T/zJ, h₀/zJ) and (T/zJ, D/zJ), where T absolute temperature, h₀, the amplitude of the oscillating field, D, crystal field interaction or single-ion anisotropy constant and z denotes the nearest-neighbor sites of the central site. The phase diagrams exhibit four fundamental phases and ten mixed phases which are composed of binary, ternary and tetrad combination of fundamental phases, depending on the crystal field interaction parameter. Moreover, the phase diagrams contain a dynamic tricritical point (T), a double critical end point (B), a multicritical point (A) and zero-temperature critical point (Z).     

Mixed-Spin Ising Model in an Oscillating Magnetic Field and Compensation Temperature

Godoy et al. (Phys. Rev. B 69, 054428, 2004) presented a study of the magnetic properties of a mixed spin (1/2,1), Ising ferrimagnetic model on a hexagonal lattice without an oscillating magnetic field. They employed dynamic mean-field calculations and Monte Carlo simulations to find the compensation point of the model and to present the phase diagrams. It has been found that the N-type compensation temperature appears only when the intrasublattice interaction between spins in the σ sublattice is ferromagnetic. Moreover, the system only undergoes a second-order phase transition. In this work, we extend the study a dynamic compensation temperature of a mixed spin-1/2 and spin-1 Ising ferrimagnetic system on a hexagonal lattice in the presence of oscillating magnetic field within the framework of dynamic mean-field calculations. We find that the system displays the N-type compensation temperature. We also calculate dynamic phase diagrams in which contain the paramagnetic, ferrimagnetic, nonmagnetic fundamental phases and two different mixed phases, depending on the interaction parameters and oscillating magnetic field. The system also exhibits tricritical and reentrant behaviors.     

Order and disorder lines in systems with competing interactions. III. Exact results from stochastic crystal growth

The methods presented in the first two articles of this series are simplified and generalized by growing stationary stochastic crystals from a given Ansatz layer. On the disorder trajectory the free energy, correlation functions, and multicritical points are calculated explicitly for a large class of models with competing interactions, including the staggered eight-vertex model, the general sixteen-vertex model, theq-state Potts model on a triangular lattice, a generalZ(q) model, and restricted spin glass models in two dimensions.     

Magnetic properties of <Superscript>3</Superscript>He on the recursive lattices

We considered the Heisenberg model on the recursive lattices with multi-spin interaction in a strong magnetic field as an approximation of the two-dimensional kagome lattice, as well as hexagonal recursive lattices as an approximation of triangular lattice, for solid ³He. In a strong magnetic field it is possible to approximate the Heisenberg model with the Izing one. Using dynamic approach, we obtain exact recursion relations for partition functions. Diagrams of the magnetization versus external magnetic field with different spin-exchange parameters and temperatures are presented. Magnetization plateaux, bifurcation points, and doublings are obtained.     

Dynamic hysteresis for Potts spin system: a Monte Carlo simulation

A Monte Carlo approach based on the Q-state Potts model is developed to describe and simulate the dynamic hysteresis of Potts spin lattice against periodic time-varying external field E. The dynamic responses of the hysteresis loops against frequency P of applied field and domain size R are studied. It is revealed that the hysteresis loops for the system energy W, polarization P, and domain wall fraction @ depend considerably on the frequency and domain size. The remnant polarization P_r shows a single-peak pattern as a function of P. The P-E loop exhibits thin rhomb and fat rhomb patterns at low P, whereas a tip-smoothed rhomb and roughly elliptical pattern is observed at high P. The loop area can be scaled with Q,P^1/3 at low P. The frequency dependency of the dynamic hysteresis is explained in terms of a simplified phenomenological model. At very small domain size, the dynamic hysteresis is significantly compressed, predicting the polarization weakening effect at small domain size.     

Hysteretic behavior of quadrupolar ordering in a 2D magnetic spin˗1 Ising nanoparticle

Pair approximation technique has been used to study the quadrupolar ordering properties for the 2D magnetic spin-1 Ising nanoparticles, consisting of core and surface parts with an interface coupling. Similar to the magnetic hysteresis, it has been shown that the quadrupole order parameter (Q) as a function of single-ion anisotropy (D) has a ‘hysteresis’ character which is generally called quadrupole (or QD) hysteresis. The observed QD loops strongly depend on temperature (T), external magnetic field (H) and biquadratic exchange interaction (K). Shifted QD loops with an asymmetry are found when H and K values are increased. These behaviors are also discussed in relation to other theoretical findings.     

Dynamical phase transition in two-dimensional fully frustrated Josephson junction arrays with resistively shunted junction dynamics

The dynamical phase transitions in two-dimensional fully frustrated Josephson junction arrays at zero temperature are investigated numerically with the resistively shunted junction model through the fluctuating twist boundary condition. The model is subjected to a driving current with nonzero orthogonal components i _x , i _y parallel to both axes of the square lattice. We find a roughly lattice size independent phase diagram with three dynamical phases: a pinned vortex lattice phase, a moving vortex lattice phase and a moving plastic phase. The phase diagram shows a direct transition from the pinned vortex to the moving vortex phase and the separation of the pinned vortex and the moving plastic phases. The time-dependent voltages v _x and v _y are periodic in the moving vortex lattice phase. But they are aperiodic in the moving plastic phase, resulting in non-monotonic characteristics and hysteresis in the current-voltage curves. It is found that the characteristic frequency is twice the time-averaged voltage in the moving vortex phase and around the time-averaged voltage in the plastic flow regime.Received: 29 May 2003, Published online: 2 October 2003PACS: 64.60.Ht Dynamic critical phenomena - 74.25.Sv Critical currents - 74.25.Fy Transport properties     

Accurate Expansions of Internal Energy and Specific Heat of Critical Two-Dimensional Ising Model with Free Boundaries

The bond-propagation algorithm for the specific heat of the two dimensional Ising model is developed and that for the internal energy is completed. Using these algorithms, we study the critical internal energy and specific heat of the model on the square lattice and triangular lattice with free boundaries. Comparing with previous works (Phys Rev E 86:041149, 2012; Phys Rev E 87:022124, 2013), we reach much higher accuracy ( $10^{-28}$ ) of the internal energy and specific heat, compared to the accuracy $10^{-11}$ of the internal energy and $10^{-9}$ of the specific heat reached in the previous works. This leads to much more accurate estimations of the edge and corner terms. The exact values of all edge and corner terms are therefore conjectured. The accurate forms of finite-size scaling for the internal energy and specific heat are determined for the rectangle-shaped square lattice with various aspect ratios and various shaped triangular lattice. For the rectangle-shaped square and triangular lattices and the triangle-shaped triangular lattice, there is no logarithmic correction terms of order higher than $1/S$ , with $S$ the area of the system. For the triangular lattice in rhombus, trapezoid and hexagonal shapes, there exist logarithmic correction terms of order higher than $1/S$ for the internal energy, and logarithmic correction terms of all orders for the specific heat.     

Nonintersecting string model and graphical approach: Equivalence with a Potts model

Using a graphical method we establish the exact equivalence of the partition function of aq-state nonintersecting string (NIS) model on an arbitrary planar, even-valenced, lattice with that of a q²-state Potts model on a related lattice. The NIS model considered in this paper is one in which the vertex weights are expressible as sums of those of basic vertex types, and the resulting Potts model generally has multispin interactions. For the square and Kagomé lattices this leads to the equivalence of a staggered NIS model with Potts models with anisotropic pair interactions, indicating that these NIS models have a first-order transition forq > 2. For the triangular lattice the NIS model turns out to be the five-vertex model of Wu and Lin and it relates to a Potts model with two- and three-site interactions. The most general model we discuss is an oriented NIS model which contains the six-vertex model and the NIS models of Stroganov and Schultz as special cases.     

Percolation of the Loss of Tension in an Infinite Triangular Lattice

We introduce a new class of bootstrap percolation models where the local rules are of a geometric nature as opposed to simple counts of standard bootstrap percolation. Our geometric bootstrap percolation comes from rigidity theory and convex geometry. We outline two percolation models: a Poisson model and a lattice model. Our Poisson model describes how defects--holes is one of the possible interpretations of these defects--imposed on a tensed membrane result in a redistribution or loss of tension in this membrane; the lattice model is motivated by applications of Hooke spring networks to problems in material sciences. An analysis of the Poisson model is given by Menshikov et al. ⁽⁴⁾ In the discrete set-up we consider regular and generic triangular lattices on the plane where each bond is removed with probability 1–p. The problem of the existence of tension on such lattice is solved by reducing it to a bootstrap percolation model where the set of local rules follows from the geometry of stresses. We show that both regular and perturbed lattices cannot support tension for any p<1. Moreover, the complete relaxation of tension--as defined in Section 4--occurs in a finite time almost surely. Furthermore, we underline striking similarities in the properties of the Poisson and lattice models.