We investigate a random field of mutually dependent random variables ("spins"), indexed by a finite one-dimensional lattice, called in physical sciences the one-dimensional Ising model, in which the random variables can take only ±1 values (see the text for a precise definition). One of the couplings, termed a "bond," that describes the mutual influence of two adjacent random variables is altered—it does not equal the others, thereby introducing a single "defect" bond. This defect bond represents a localised perturbation within an otherwise uniform system. Utilising the recurrence relations of Chebyshev polynomials and the bijective map between the number of spins and the polynomial index, we present a new method for calculations and systematically explore, using it, the system’s properties across different chain lengths and boundary conditions. As an application, we derive analytical expressions for the dependence of the average values of the random variables on their position within the chain, which we refer to as the "local magnetisation profile". From the findings related to the system with a defect bond, we present a novel result for this profile under free (Dirichlet) boundary conditions and re-derive the corresponding result for antiperiodic boundary conditions.
We investigate a random field of mutually dependent random variables ("spins"), indexed by a finite one-dimensional lattice, called in physical sciences the one-dimensional Ising model, in which the random variables can take only ±1 values (see the text for a precise definition). One of the couplings, termed a "bond," that describes the mutual influence of two adjacent random variables is altered—it does not equal the others, thereby introducing a single "defect" bond. This defect bond represents a localised perturbation within an otherwise uniform system. Utilising the recurrence relations of Chebyshev polynomials and the bijective map between the number of spins and the polynomial index, we present a new method for calculations and systematically explore, using it, the system’s properties across different chain lengths and boundary conditions. As an application, we derive analytical expressions for the dependence of the average values of the random variables on their position within the chain, which we refer to as the "local magnetisation profile". From the findings related to the system with a defect bond, we present a novel result for this profile under free (Dirichlet) boundary conditions and re-derive the corresponding result for antiperiodic boundary conditions.
For studying of the finite-size behavior of the Ising model under different boundary conditions we propose an alternative to the standard transfer matrix technique approach based on Abel\`{e}s theorem and Chebyshev polynomials. Using it we easily reproduce the known for periodic boundary conditions results for Lee-Yang zeros, the exact position space renormalization group transformation, etc., and extend them deriving new results for antiperiodic boundary conditions. Note that in the latter case one has a nontrivial order-parameter profile, which we also calculate, where the average value of a given spin depends on the distance from the seam with opposite bond in the system. It is interesting to stress that under both boundary conditions, the one-dimensional case exhibits Schottky anomaly.
The one-dimensional Ising model with general spin S and either open or periodic boundary conditions has been formulated as a difference equation of order 2S+1 with non-zero magnetic field B. This difference equation has the characteristic equation of the transfer matrix, V, as its indicial equation. When B=0, for the open chain, the order of the difference equation is reduced to [S+1]. It is shown that the correlation functions for the open chain satisfy similar difference equations. This result explicitly displays the lack of translational invariance of the open chain. For several cases it is shown that the partition function can be represented as the difference of two Chebyshev polynomials of the second kind.
The quantum wetting transition in the one-dimensional transverse-field Ising model with random bonds is studied. Opposite boundary fields are applied at the two ends of the Ising chain. The interface between two domains with oppositely oriented magnetization is localized near the boundary or in the middle of the chain. The wetting transition refers to the phase transition of localization and delocalization of the interface. The tuning parameter hL is the boundary field at one end, e.g., the left end. At the transition point hL=hw, the wetting transition occurs. First, we study the wetting transition with one defect bond. It is a first-order transition, and the interface jumps from the left end of the chain to the defect bond at the transition point. Second, we study the wetting transition with two defect bonds and find that the weaker defect bond dominates the phase transition. The wetting transition is still first order, and the interface is localized at the left end of the chain or at the weaker defect bond. Lastly, the random bond case is studied. The random bonds have a rectangular distribution. The wetting transition is still first order. The finite-size effects are studied. The statistics of the transition points and the energy gaps are obtained. We study lattices with sizes N=100, 150, 200, 250, 300, and 350. The deviation of the phase transition points hw does not decrease as the lattice size increases. It is argued that the transition point is sample dependent, i.e., the variation in the transition point hw does not approach zero, even at the thermodynamic limit.
Following the discovery of high Tc superconductivity in the copper oxides, there has been a great deal of interest in the RVB wave function proposed by Anderson [1]. As a warm-up exercise we have considered a valence-bond wave function for the one dimensional spin-1/2 Heisenberg chain. The main virtue of our work is to propose a new variational singlet wavefunction which is almost analytically tractable by a transfer-matrix technique. We have obtained the ground state energy for finite as well as infinite chains, in good agreement with exact results. Correlation functions, excited states, and the effects of other interactions (e.g., spin-Peierls) are also accessible within this scheme [2]. Since the ground state of the chain is known to be a singlet (Lieb & Mattis [3]), we write the appropriate wave function as a superposition of valence-bond singlets, [Formula: see text], where | k > is a spin configuration obtained by pairing all spins into singlet pairs, in a way which is common in valence-bond calculations of large molecules. As in that case, each configuration, | k >, can be represented by a Rümer diagram, with directed bonds connecting each pair of spins on the chain. The c k 's are variational co-efficients, the form of which is determined as follows: Each singlet configuration (Rümer diagram) is divided into "zones", a "zone" corresponding to the region between two consecutive sites. Each zone is indexed by its distance from the end of the chain and by the number of bonds crossing it. Our procedure assigns a variational parameter, x ij , to the j th zone, when crossed by i bonds. The resulting wavefunction for an N-site chain is written as [Formula: see text] where m ij(k) equals 1 when zone j is crossed by i bonds and zero otherwise. To make the calculation tractable we reduce the number of variational parameters by disallowing configurations with bonds connecting any two sites separated by more than 2M lattice points. (For simplicity, we have limited ourselves to M=3, but the scheme can be used for any M). With the simple ansatz, matrix elements can be calculated by a transfer-matrix method. To understand the transfer-matrix method note that since only local zone parameters appear in the description of each state | k >, matrix elements and overlaps, [Formula: see text] and < k | k '>, are completely specified by a small number of "local states" associated with each zone. Within a given zone a local state is defined by (i) the number of bonds crossing the zone and (ii), by whether the bonds originate from the initial (| k >) or final (| k '>) state. It is then easy to see that "local states" of consecutive zones are connected by a 15 × 15 transfer matrix (for the case M=3). Furthermore, the overlap matrix element can be written as a product of transfer-matrices associated with each zone of the chain. When calculating matrix elements of the Hamiltonian, an additional matrix, U , must be defined, to represent the particular zone involving the two spins connected by the Heisenberg interaction. The relevant details as well as the comparison with exact results will be given elsewhere. We are planning to ultimately apply this method to the two dimensional case, and hope to include the effects of holes.
We develop a recently proposed mapping of the two-dimensional Ising model with random exchange (RBIM) via the transfer matrix, to a network model for a disordered system of noninteracting fermions. The RBIM transforms in this way to a localization problem belonging to one of a set of nonstandard symmetry classes, known as class D; the transition between paramagnet and ferromagnet is equivalent to a delocalization transition between an insulator and a quantum Hall conductor. We establish the mapping as an exact and efficient tool for numerical analysis: using it, the computational effort required to study a system of width M is proportional to M3, and not exponential in M as with conventional algorithms. We show how the approach may be used to calculate for the RBIM the free energy, typical correlation lengths in quasi-one dimension for both the spin and the disorder operators, and the even powers of spin-spin correlation functions and their disorder averages. We examine in detail the square-lattice, nearest-neighbor $\ifmmode\pm\else\textpm\fi{}J$ RBIM, in which bonds are independently antiferromagnetic with probability p, and ferromagnetic with probability $1\ensuremath{-}p.$ Studying temperatures T>0.4J, we obtain precise coordinates in the $p\ensuremath{-}T$ plane for points on the phase boundary between ferromagnet and paramagnet, and for the multicritical (Nishimori) point. We demonstrate scaling flow towards the pure Ising fixed point at small p, and determine critical exponents at the multicritical point.
The one-dimensional Ising model with arbitrary nearest neighbor interactions and a magnetic field is formulated in terms of the spin correlation function relations derived in the preceding article I. The local magnetization and the two-spin correlation functions are shown to satisfy a set of coupled algebraic summation equations. The solution of this set is given in terms of a power series in hyberbolic tangents of the magnetic field strength, for chains of arbitrary length. For cases in which the interaction constants are random variables, explicit expressions are derived for the configurational averaged magnetization and the two-spin correlation functions to the lowest order in the magnetic field.