Dynamical Renormalization Group Research Articles

Collective dynamics of dividing chemotactic cells.

The large scale behavior of a population of cells that grow and interact through the concentration field of the chemicals they secrete is studied using dynamical renormalization group methods. The combination of the effective long-range chemotactic interaction and lack of number conservation leads to a rich variety of phase behavior in the system, which includes a sharp transition from a phase that has moderate (or controlled) growth and regulated chemical interactions to a phase with strong (or uncontrolled) growth and no chemical interactions. The transition point has nontrivial critical exponents. Our results might help shed light on the interplay between chemical signaling and growth in tissues and colonies, and in particular on the challenging problem of cancer metastasis.

Skewness in (1+1)-dimensional Kardar-Parisi-Zhang-type growth.

We use the (1+1)-dimensional Kardar-Parisi-Zhang equation driven by a Gaussian white noise and employ the dynamic renormalization-group of Yakhot and Orszag without rescaling [J. Sci. Comput. 1, 3 (1986)]. Hence we calculate the second- and third-order moments of height distribution using the diagrammatic method in the large-scale and long-time limits. The moments so calculated lead to the value S=0.3237 for the skewness. This value is comparable with numerical and experimental estimates.

Chebyshev expansion for impurity models using matrix product states

We improve a recently developed expansion technique for calculating real frequency spectral functions of any one-dimensional model with short-range interactions, by postprocessing computed Chebyshev moments with linear prediction. This can be achieved at virtually no cost and, in sharp contrast to existing methods based on the dampening of the moments, improves the spectral resolution rather than lowering it. We validate the method for the exactly solvable resonating level model and the single impurity Anderson model. It is capable of resolving sharp Kondo resonances, as well as peaks within the Hubbard bands when employed as an impurity solver for dynamical mean-field theory (DMFT). Our method works at zero temperature and allows for arbitrary discretization of the bath spectrum. It achieves similar precision as the dynamical density matrix renormalization group (DDMRG), at lower cost. We also propose an alternative expansion, of 1-exp(-tau H) instead of the usual H, which opens the possibility of using established methods for the time evolution of matrix product states to calculate spectral functions directly.

Dynamical diffusion and renormalization group equation for the Fermi velocity in doped graphene

The aim of this work is to study the electron transport in graphene with impurities by introducing a generalization of linear response theory for linear dispersion relations and spinor wave functions. Current response and density response functions are derived and computed in the Boltzmann limit showing that in the former case a minimum conductivity appears in the no-disorder limit. In turn, from the generalization of both functions, an exact relation can be obtained that relates both. Combining this result with the relation given by the continuity equation it is possible to obtain general functional behavior of the diffusion pole. Finally, a dynamical diffusion is computed in the quasistatic limit using the definition of relaxation function. A lower cutoff must be introduced to regularize infrared divergences which allow us to obtain a full renormalization group equation for the Fermi velocity, which is solved up to order O(ℏ2).

Dynamical correlation functions for the one-dimensional Bose-Hubbard insulator

We calculate the dynamical current and kinetic-energy correlation functions for the first Mott lobe of the one-dimensional Bose-Hubbard model. We employ the strong-coupling expansion up to sixth order in $x=t/U$ and the dynamical density-matrix renormalization group method on rings with 64 sites. The correlation functions are finite above the single-particle gap with a square-root onset, as also found from field theory close to the Mott transition. The correlation functions display a featureless superposition of the primary and tertiary Hubbard bands. We find very good agreement between all methods in the interaction and frequency regimes where they are applicable.

Monte Carlo renormalization group of dilute 3D Ising dynamics

Simulation of critical relaxation of the magnetization in a three-dimensional Ising model diluted by quenched nonmagnetic impurity atoms is carried out. The cubic systems with linear size L = 144 and spin concentrations p = 0.95, 0.8, 0.6, and 0.5 are studied. The dynamical critical exponent z is determined by the Monte Carlo method combined with the dynamical renormalization group method. The obtained results are: z(p): z(0.95) = 2.19 ± 0.07, z(0.8) = 2.29 ± 0.06, z(0.6) = 2.72 ± 0.08, and z(0.5) = 2.75 ± 0.08.

Blind deconvolution of density-matrix renormalization-group spectra

We present a numerical method for calculating piecewise smooth spectral functions of correlated quantum systems in the thermodynamic limit from the spectra of finite systems computed using the dynamical or correction-vector density-matrix renormalization group method. The key idea is to consider this problem as a blind deconvolution with an unknown kernel which causes both a broadening and finite-size corrections of the spectrum. In practice, the method reduces to a least-square optimization under non-linear constraints which enforce the positivity and piecewise smoothness of spectral functions. The method is demonstrated on the single-particle density of states of one-dimensional paramagnetic Mott insulators represented by the half-filled Hubbard model on an open chain. Our results confirm that the density of states has a step-like shape but no square-root singularity at the spectrum onset.

Strong anisotropy in two-dimensional surfaces with generic scale invariance: nonlinear effects.

We expand a previous study [Phys. Rev. E 86, 051611 (2012)] on the conditions for occurrence of strong anisotropy in the scaling properties of two-dimensional surfaces displaying generic scale invariance. In that study, a natural scaling ansatz was proposed for strongly anisotropic systems, which arises naturally when analyzing data from, e.g., thin-film production experiments. The ansatz was tested in Gaussian (linear) models of surface dynamics and in nonlinear models, like the Hwa-Kardar (HK) equation [Phys. Rev. Lett. 62, 1813 (1989)], which are susceptible of accurate approximations through the former. In contrast, here we analyze nonlinear equations for which such approximations fail. Working within generically scale-invariant situations, and as representative case studies, we formulate and study a generalization of the HK equation for conserved dynamics and reconsider well-known systems, such as the conserved and the nonconserved anisotropic Kardar-Parisi-Zhang equations. Through the combined use of dynamic renormalization group analysis and direct numerical simulations, we conclude that the occurrence of strong anisotropy in two-dimensional surfaces requires dynamics to be conserved. We find that, moreover, strong anisotropy is not generic in parameter space but requires, rather, specific forms of the terms appearing in the equation of motion, whose justification needs detailed information on the dynamical process that is being modeled in each particular case.

Invited review: KPZ. Recent developments via a variational formulation

Recently, a variational approach has been introduced for the paradigmatic Kardar-Parisi-Zhang (KPZ) equation. Here we review that approach, together with the functional Taylor expansion that the KPZ nonequilibrium potential (NEP) admits. Such expansion becomes naturally truncated at third order, giving rise to a nonlinear stochastic partial differential equation to be regarded as a gradient-flow counterpart to the KPZ equation. A dynamic renormalization group analysis at one-loop order of this new mesoscopic model yields the KPZ scaling relation $\alpha + z = 2$, as a consequence of the exact cancelation of the different contributions to vertex renormalization. This result is quite remarkable, considering the lower degree of symmetry of this equation, which is in particular not Galilean invariant. In addition, this scheme is exploited to inquire about the dynamical behavior of the KPZ equation through a path-integral approach. Each of these aspects offers novel points of view and sheds light on particular aspects of the dynamics of the KPZ equation. Received: 25 June 2013, Accepted: 10 December 2013; Reviewed by: F. Reis, Instituto de Fisica, Univ. Fed. Fluminense, Brazil; Edited by: R. Dickman; DOI: http://dx.doi.org/10.4279/PIP.050010Cite as: H S Wio, R R Deza, C Escudero, J A Revelli, Papers in Physics 5, 050010 (2013)

The δN formula is the dynamical renormalization group

We derive the `separate universe' methodfor the inflationary bispectrum,beginning directly froma field-theory calculation.We work to tree-level in quantum effects but to all orders in theslow-roll expansion,with masses accommodated perturbatively.Our method provides a systematic basis to account fornovel sources of time-dependencein inflationary correlation functions,and has immediate applications.First, we use our result to obtain the correct matching prescriptionbetween the `quantum' and `classical' parts of theseparate universe computation.Second, we elaborate on theapplication of this method in situations where its validity is not clear.As a by-product of our calculation we give the leading slow-rollcorrections to the three-point function of field fluctuations onspatially flat hypersurfaces in a canonical, multiple-field model.

Model Evidence of an Anisotropic Chirald+id-Wave Pairing State for the Water-IntercalatedNaxCoO2·yH2OSuperconductor

Since its discovery, the superconducting phase in water-intercalated sodium cobaltates NaxCoO2 yH2O (x~0.3, y~1.3) has posed fundamental challenges in terms of experimental investigation and theoretical understanding. By a combined dynamical mean-field and renormalization group approach, we find an anisotropic chiral d+id wave state as a consequence of multi-orbital effects, Fermi surface topology, and magnetic fluctuations. It naturally explains the singlet property and close-to-nodal gap features of the superconducting phase as indicated by experiments.

Nonequilibrium photo dynamics of low-dimensional strongly correlated electron systems

One of the outstanding contemporary challenges in condensed matter physics is to understand the dynamics of interacting quantum systems exposed to an external perturbation. We theoretically examine nonequilibrium photo dynamics and its interplay of charge, spin, and lattice degrees of freedom on a Hubbard-Holstein chain in one dimension and a t-J-Holstein square lattice in two dimensions. In the chain, performing dynamical density-matrix renormalization group calculations, we find that many phonons generated dynamically after photo irradiation in Mott insulators cause initial relaxation process. On the other hand, in the square lattice with model parameters as relevant for cuprates, a Lanczos-type exact diagonalization calculation shows that the majority of absorbed energy flows into spin subsystem rather than phonon subsystem.

A numerical renormalisation group method for the analysis of critical spreading activity in spiking neural networks

In recent years criticality has been used as a mechanism for understanding the presence of long-range correlations and intermittence in macroscopic neural activity. Little work has been done to demonstrate the usefulness of this paradigm for understanding the statistics and dynamics of neural networks. The explanation for self-similarity and universality is grounded in the renormalisation group (RG) theory of critical phenomena [1]. Until recently [2], this theory has not been leveraged for the inspection and understanding of neural network models or in neural data analysis. Here, we present a range of statistics derived from a dynamic renormalisation group theory [3], shown to elucidate the nature of self-organised criticality in sandpile models, which we have operationalised for the numeric data analysis of spike trains. The relationship between universality via self-similarity, the power-law exponents and a set of renormalised parameters is explored in the case of two automata models to demonstrate the application of the RG. In particular, we focus on the temporal profile of avalanches, the branching ratio and the avalanche size distribution and how these can be related to the renormalised parameters we study in these systems. Next, the same statistics are used to investigate the universality of a simplified integrate-and-fire model. Finally we use these methods to assess the critical state on a number of generic neural network models which are used to study the balanced state of neocortex.

Many-Body Localization in One Dimension as a Dynamical Renormalization Group Fixed Point

We formulate a dynamical real space renormalization group (RG) approach to describe the time evolution of a random spin-1/2 chain, or interacting fermions, initialized in a state with fixed particle positions. Within this approach we identify a many-body localized state of the chain as a dynamical infinite randomness fixed point. Near this fixed point our method becomes asymptotically exact, allowing analytic calculation of time dependent quantities. In particular, we explain the striking universal features in the growth of the entanglement seen in recent numerical simulations: unbounded logarithmic growth delayed by a time inversely proportional to the interaction strength. This is in striking contrast to the much slower entropy growth as loglogt found for noninteracting fermions with bond disorder. Nonetheless, even the interacting system does not thermalize in the long time limit. We attribute this to an infinite set of approximate integrals of motion revealed in the course of the RG flow, which become asymptotically exact conservation laws at the fixed point. Hence we identify the many-body localized state with an emergent generalized Gibbs ensemble.

Fixing the fixed-point system—Applying Dynamic Renormalization Group to systems with long-range interactions

In this paper a mode of using the Dynamic Renormalization Group (DRG) method is suggested in order to cope with inconsistent results obtained when applying it to a continuous family of one-dimensional nonlocal models. The key observation is that the correct fixed-point dynamical system has to be identified during the analysis in order to account for all the relevant terms that are generated under renormalization. This is well established for static problems, however poorly implemented in dynamical ones. An application of this approach to a nonlocal extension of the Kardar-Parisi-Zhang equation resolves certain problems in one-dimension. Namely, obviously problematic predictions are eliminated and the existing exact analytic results are recovered.

Dynamic density-density correlations in interacting Bose gases on optical lattices

In order to identify possible experimental signatures of the superfluid to Mott-insulator quantum phase transition we calculate the charge structure factor S(k, ω) for the one-dimensional Bose-Hubbard model using the dynamical density-matrix renormalisation group (DDMRG) technique. Particularly we analyse the behaviour of S(k, ω) by varying – at zero temperature–the Coulomb interaction strength within the first Mott lobe. For strong interactions, in the Mott-insulator phase, we demonstrate that the DDMRG results are well reproduced by a strong-coupling expansion, just as the quasi-particle dispersion. In the super-fluid phase we determine the linear excitation spectrum near k = 0. In one dimension, the amplitude mode is absent which mean-field theory suggests for higher dimensions.

Local density of states of the one-dimensional spinless fermion model

We investigate the local density of states of the one-dimensional half-filled spinless fermion model with nearest-neighbour hopping t > 0 and interaction V in its Luttinger liquid phase −2t < V ≤ 2t. The bulk density of states and the local density of states in open chains are calculated over the full band width ∼4t with an energy resolution ≤0.08t using the dynamical density-matrix renormalization group (DDMRG) method. We also perform DDMRG simulations with a resolution of 0.01t around the Fermi energy to reveal the power-law behaviour D(ϵ) ∼ |ϵ − ϵF|α predicted by the Luttinger liquid theory for bulk and boundary density of states. The exponents α are determined using a finite-size scaling analysis of DDMRG data for lattices with up to 3200 sites. The results agree with the exact exponents given by the Luttinger liquid theory combined with the Bethe Ansatz solution. The crossover from boundary to bulk density of states is analysed. We have found that boundary effects can be seen in the local density of states at all energies, even far away from the chain edges.

Discrete antiferromagnetic spin-wave excitations in the giant ferric wheel Fe18

The low-temperature elementary spin excitations in the AFM molecular wheel Fe18 were studied experimentally by inelastic neutron scattering and theoretically by modern numerical methods, such as dynamical density matrix renormalization group or quantum Monte Carlo techniques, and analytical spin-wave theory calculations. Fe18 involves eighteen spin-5/2 Fe(III) ions with a Hilbert space dimension of 10^14, constituting a physical system that is situated in a region between microscopic and macroscopic. The combined experimental and theoretical approach allowed us to characterize and discuss the magnetic properties of Fe18 in great detail. It is demonstrated that physical concepts such as the rotational-band or L&E-band concepts developed for smaller rings are still applicable. In particular, the higher-lying low-temperature elementary spin excitations in Fe18 or AFM wheels in general are of discrete antiferromagnetic spin-wave character.

Characterization of Mott-insulating and superfluid phases in the one-dimensional Bose-Hubbard model

We use strong-coupling perturbation theory, the variational cluster approach (VCA), and the dynamical density-matrix renormalization group (DDMRG) method to investigate static and dynamical properties of the one-dimensional Bose--Hubbard model in both the Mott-insulating and superfluid phases. From the von Neumann entanglement entropy we determine the central charge and the transition points for the first two Mott lobes. Our DMRG results for the ground-state energy, momentum distribution function, boson correlation function decay, Mott gap, and single particle-spectral function are reproduced very well by the strong-coupling expansion to fifth order, and by VCA with clusters up to 12 sites as long as the ratio between the hopping amplitude and on-site repulsion, t/U, is smaller than 0.15 and 0.25, respectively. In addition, in the superfluid phase VCA captures well the ground-state energy and the sound velocity of the linear phonon modes. This comparison provides an authoritative estimate for the range of applicability of these methods. In strong-coupling theory for the Mott phase, the dynamical structure factor is obtained from the solution of an effective single-particle problem with an attractive potential. The resulting resonances show up as double-peak structure close to the Brillouin zone boundary. These high-energy features also appear in the superfluid phase which is characterized by a pronounced phonon mode at small momenta and energies, as predicted by Bogoliubov and field theory. In one dimension, there are no traces of an amplitude mode in the dynamical single-particle and two-particle correlation functions.

Transport coefficients ofO(N)scalar field theories close to the critical point

We investigate the critical dynamics of $O(N)$-symmetric scalar field theories to determine the critical exponents of transport coefficients as a second-order phase transition is approached from the symmetric phase. A set of stochastic equations of motion for the slow modes is formulated, and the long-wavelength dynamics is examined for an arbitrary number of field components $N$ in the framework of the dynamical renormalization group within the $\ensuremath{\epsilon}$ expansion. We find that for a single component scalar field theory, $N=1$, the system reduces to the model C of critical dynamics, whereas for $N&gt;1$, the model G is effectively restored owing to dominance of $O(N)$-symmetric charge fluctuations. In both cases, the shear viscosity remains finite in the critical region. On the other hand, we find that the bulk viscosity diverges as the correlation length squared, for $N=1$, while it remains finite for $N&gt;1$.