Entanglement entropy dynamics of Heisenberg chains

By means of the time dependent density matrix renormalization group algorithm we study the zero-temperature dynamics of the Von Neumann entropy of a block of spins in a Heisenberg chain after a sudden quench in the anisotropy parameter. In the absence of any disorder the block entropy increases linearly with time and then saturates. We analyse the velocity of propagation of the entanglement as a function of the initial and final anisotropies and compare our results, wherever possible, with those obtained by means of conformal field theory. In the disordered case we find a slower ( logarithmic) evolution which may signal the onset of entanglement localization.

Phase diagram of spin-1 bosons on one-dimensional lattices

Spinor Bose condensates loaded in optical lattices have a rich phase diagram characterized by different magnetic order. Here we apply the density matrix renormalization group to accurately determine the phase diagram for spin-1 bosons loaded on a one-dimensional lattice. The Mott lobes present an even or odd asymmetry associated to the boson filling. We show that for odd fillings the insulating phase is always in a dimerized state. The results obtained in this work are also relevant for the determination of the ground state phase diagram of the S=1 Heisenberg model with biquadratic interaction.

Bilinear-biquadratic spin-1 chain undergoing quadratic Zeeman effect

The Heisenberg model for spin-1 bosons in one dimension presents many different quantum phases, including the famous topological Haldane phase. Here we study the robustness of such phases in front of a SU(2) symmetry-breaking field as well as the emergence of unique phases. Previous studies have analyzed the effect of such uniaxial anisotropy in some restricted relevant points of the phase diagram. Here we extend those studies and present the complete phase diagram of the spin-1 chain with uniaxial anisotropy. To this aim, we employ the density-matrix renormalization group together with analytical approaches. The complete phase diagram can be realized using ultracold spinor gases in the Mott insulator regime under a quadratic Zeeman effect.

Strong electronic correlation in the hydrogen chain: A variational Monte Carlo study

In this paper, we report a fully ab initio variational Monte Carlo study of the linear and periodic chain of hydrogen atoms, a prototype system providing the simplest example of strong electronic correlation in low dimensions. In particular, we prove that numerical accuracy comparable to that of benchmark density-matrix renormalization-group calculations can be achieved by using a highly correlated Jastrow-antisymmetrized geminal power variational wave function. Furthermore, by using the so-called "modern theory of polarization" and by studying the spin-spin and dimer-dimer correlations functions, we have characterized in detail the crossover between the weakly and strongly correlated regimes of this atomic chain. Our results show that variational Monte Carlo provides an accurate and flexible alternative to highly correlated methods of quantum chemistry which, at variance with these methods, can be also applied to a strongly correlated solid in low dimensions close to a crossover or a phase transition.

Scaling of the entanglement spectrum near quantum phase transitions

The entanglement spectrum describing quantum correlations in many-body systems has been recently recognized as a key tool to characterize different quantum phases, including topological ones. Here we derive its analytically scaling properties in the vicinity of some integrable quantum phase transitions and extend our studies also to nonintegrable quantum phase transitions in one-dimensional spin models numerically. Our analysis shows that, in all studied cases, the scaling of the difference between the two largest nondegenerate Schmidt eigenvalues yields with good accuracy critical points and mass scaling exponents.

Full characterization of the quantum linear-zigzag transition in atomic chains

A string of repulsively interacting particles exhibits a phase transition to a zigzag structure, by reducing the transverse trap potential or the interparticle distance. Based on the emergent symmetry Z2 it has been argued that this instability is a quantum phase transition, which can be mapped to an Ising model in transverse field. An extensive Density Matrix Renormalization Group analysis is performed, resulting in an high-precision evaluation of the critical exponents and of the central charge of the system, confirming that the quantum linear-zigzag transition belongs to the critical Ising model universality class. Quantum corrections to the classical phase diagram are computed, and the range of experimental parameters where quantum effects play a role is provided. These results show that structural instabilities of one-dimensional interacting atomic arrays can simulate quantum critical phenomena typical of ferromagnetic systems.

Analytical Study of Light Propagation in Highly Nonlinear Media

The present dissertation is devoted to the construction of exact and approximate analytical solutions of the problem of light propagation in highly nonlinear media. It is demonstrated that for many experimental conditions, the problem can be studied under the geometrical optics approximation with a sufficient accuracy. Based on the renormalization group symmetry analysis, exact analytical solutions of the eikonal equations with a higher order refractive index are constructed. A new analytical approach to the construction of approximate solutions is suggested. Based on it, approximate solutions for various boundary conditions, nonlinear refractive indices and dimensions are constructed. Exact analytical expressions for the nonlinear self-focusing positions are deduced. On the basis of the obtained solutions a general rule for the single filament intensity is derived; it is demonstrated that the scaling law (the functional dependence of the self-focusing position on the peak beam intensity) is defined by a form of the nonlinear refractive index but not the beam shape at the boundary. Comparisons of the obtained solutions with results of experiments and numerical simulations are discussed.

Limitations of small x resummation methods from F2data

We discuss several methods of calculating the DIS structure functions F2(x,Q2) based on BFKL-type small x resummations. Taking into account new HERA data ranging down to small xand low Q2, the pure leading order BFKL-based approach is excluded. Other methods based on high energy factorization are closer to conventional renormalization group equations. Despite several difficulties and ambiguities in combining the renormalization group equations with small x resummed terms, we find that a fit to the current data is hardly feasible, since the data in the low Q2 region are not as steep as the BFKL formalism predicts. Thus we conclude that deviations from the (successful) renormalization group approach towards summing up logarithms in 1/x are disfavoured by experiment.

Small x resummations confronted with F2(x, Q2) data

It has been observed recently that a consistent LO BFKL gluon evolution leads to a steep growth of F2(x, Q2) for x → 0 almost independently of Q2. We show that current data from the DESY HERA collider are precise enough to finally rule out a pure BFKL behaviour in the accessible small x region. Several attempts have been made by other groups to treat the BFKL type small x resummations instead as additions to the conventional anomalous dimensions of the successful renormalization group “Altarelli-Parisi” equations. We demonstrate that all presently available F2 data, in particular at lower values of Q2, can not be described using the presently known NLO (two-loop consistent) small x resummations. Finally we comment on the common reason for the failure of these BFKL inspired methods which result, in general, in too steep >x-dependencies as x → 0.

Balitskiǐ-Fadin-Kuraev-Lipatov equation versus data from DESY HERA

The BFKL equation and the kT-factorization theorem are used to obtain predictions for F2 in the small Bjo/rken-x region over a wide range of Q2. The dependence on the parameters, especially on those concerning the infrared region, is discussed. After a background fit to recent experimental data obtained at DESY HERA and at Fermilab (E665 experiment) we find that the predicted, almost Q2 independent BFKL slope λ≳0.5 appears to be too steep at lower Q2 values. Thus there seems to be a chance that future HERA data can distinguish between pure BFKL and conventional field theoretic renormalization group approaches. © 1995 The American Physical Society.

Hierarchical spherical model from a geometric point of view

A continuous version of the hierarchical spherical model at dimension d=4 is investigated. Two limit distributions of the block spin variable X(gamma), normalized with exponents gamma = d + 2 and gamma=d at and above the critical temperature, are established. These results are proven by solving certain evolution equations corresponding to the renormalization group (RG) transformation of the O(N) hierarchical spin model of block size L(d) in the limit L down arrow 1 and N ->infinity. Starting far away from the stationary Gaussian fixed point the trajectories of these dynamical system pass through two different regimes with distinguishable crossover behavior. An interpretation of this trajectories is given by the geometric theory of functions which describe precisely the motion of the Lee-Yang zeroes. The large-N limit of RG transformation with L(d) fixed equal to 2, at the criticality, has recently been investigated in both weak and strong (coupling) regimes by Watanabe (J. Stat. Phys. 115:1669-1713, 2004) . Although our analysis deals only with N = infinity case, it complements various aspects of that work.

Influence of super-ohmic dissipation on a disordered quantum critical point

We investigate the combined influence of quenched randomness and dissipation on a quantum critical point with O(N) order-parameter symmetry. Utilizing a strong-disorder renormalization group, we determine the critical behavior in one space dimension exactly. For super-ohmic dissipation, we find a Kosterlitz-Thouless type transition with conventional (power-law) dynamical scaling. The dynamical critical exponent depends on the spectral density of the dissipative baths. We also discuss the Griffiths singularities, and we determine observables.

Thermal dependence of the zero-bias conductance through a nanostructure

We show that the conductance of a quantum wire side-coupled to a quantum dot, with a gate potential favoring the formation of a dot magnetic moment, is a universal function of the temperature. Universality prevails even if the currents through the dot and the wire interfere. We apply this result to the experimental data of Sato et al. (Phys. Rev. Lett., 95 (2005) 066801). Copyright (C) EPLA, 2009

Estudo das propriedades críticas do processo epidêmico por par com difusão de pares

The pair contact process - PCP is a nonequilibrium stochastic model which, like the basic contact process - CP, exhibits a phase transition to an absorbing state. While the absorbing state CP corresponds to a unique configuration (empty lattice), the PCP process infinitely many. Numerical and theoretical studies, nevertheless, indicate that the PCP belongs to the same universality class as the CP (direct percolation class), but with anomalies in the critical spreading dynamics. An infinite number of absorbing configurations arise in the PCP because all process (creation and annihilation) require a nearest-neighbor pair of particles. The diffusive pair contact process - PCPD) was proposed by Grassberger in 1982. But the interest in the problem follows its rediscovery by the Langevin description. On the basis of numerical results and renormalization group arguments, Carlon, Henkel and Schollwöck (2001), suggested that certain critical exponents in the PCPD had values similar to those of the party-conserving - PC class. On the other hand, Hinrichsen (2001), reported simulation results inconsistent with the PC class, and proposed that the PCPD belongs to a new universality class. The controversy regarding the universality of the PCPD remains unresolved. In the PCPD, a nearest-neighbor pair of particles is necessary for the process of creation and annihilation, but the particles to diffuse individually. In this work we study the PCPD with diffusion of pair, in which isolated particles cannot move; a nearest-neighbor pair diffuses as a unit. Using quasistationary simulation, we determined with good precision the critical point and critical exponents for three values of the diffusive probability: D=0.5 and D=0.1. For D=0.5: PC=0.89007(3), β/v=0.252(9), z=1.573(1), =1.10(2), m=1.1758(24). For D=0.1: PC=0.9172(1), β/v=0.252(9), z=1.579(11), =1.11(4), m=1.173(4)

Fases e criticalidade no modelo ashkin - teller de três cores

The usual Ashkin-Teller (AT) model is obtained as a superposition of two Ising models coupled through a four-spin interaction term. In two dimension the AT model displays a line of fixed points along which the exponents vary continuously. On this line the model becomes soluble via a mapping onto the Baxter model. Such richness of multicritical behavior led Grest and Widom to introduce the N-color Ashkin-Teller model (N-AT). Those authors made an extensive analysis of the model thus introduced both in the isotropic as well as in the anisotropic cases by several analytical and computational methods. In the present work we define a more general version of the 3-color Ashkin-Teller model by introducing a 6-spin interaction term. We investigate the corresponding symmetry structure presented by our model in conjunction with an analysis of possible phase diagrams obtained by real space renormalization group techniques. The phase diagram are obtained at finite temperature in the region where the ferromagnetic behavior is predominant. Through the use of the transmissivities concepts we obtain the recursion relations in some periodical as well as aperiodic hierarchical lattices. In a first analysis we initially consider the two-color Ashkin-Teller model in order to obtain some results with could be used as a guide to our main purpose. In the anisotropic case the model was previously studied on the Wheatstone bridge by Claudionor Bezerra in his Master Degree dissertation. By using more appropriated computational resources we obtained isomorphic critical surfaces described in Bezerra's work but not properly identified. Besides, we also analyzed the isotropic version in an aperiodic hierarchical lattice, and we showed how the geometric fluctuations are affected by such aperiodicity and its consequences in the corresponding critical behavior. Those analysis were carried out by the use of appropriated definitions of transmissivities. Finally, we considered the modified 3-AT model with a 6-spin couplings. With the inclusion of such term the model becomes more attractive from the symmetry point of view. For some hierarchical lattices we derived general recursion relations in the anisotropic version of the model (3-AAT), from which case we can obtain the corresponding equations for the isotropic version (3-IAT). The 3-IAT was studied extensively in the whole region where the ferromagnetic couplings are dominant. The fixed points and the respective critical exponents were determined. By analyzing the attraction basins of such fixed points we were able to find the three-parameter phase diagram (temperature £ 4-spin coupling £ 6-spin coupling). We could identify fixed points corresponding to the universality class of Ising and 4- and 8-state Potts model. We also obtained a fixed point which seems to be a sort of reminiscence of a 6-state Potts fixed point as well as a possible indication of the existence of a Baxter line. Some unstable fixed points which do not belong to any aforementioned q-state Potts universality class was also found