We introduce a general method for optimizing real-space renormalization-group transformations to study the critical properties of a classical system. The scheme is based on minimizing the Kullback-Leibler divergence between the distribution of the system and the normalized normalizing factor of the transformation parametrized by a restricted Boltzmann machine. We compute the thermal critical exponent of the two-dimensional Ising model using the trained optimal projector and obtain a very accurate thermal critical exponent $y_t=1.0001(11)$ after the first step of the transformation.
Discrete amorphous materials are best described in terms of arbitrary networks which can be embedded in three dimensional space. Investigating the thermodynamic equilibrium as well as non-equilibrium behavior of such materials around second order phase transitions call for special techniques. We set up a renormalization group scheme by expanding an arbitrary scalar field living on the nodes of an arbitrary network, in terms of the eigenvectors of the normalized graph Laplacian. The renormalization transformation involves, as usual, the integration over the more rapidly varying components of the field, corresponding to eigenvectors with larger eigenvalues, and then rescaling. The critical exponents depend on the particular graph through the spectral density of the eigenvalues.
We show a way to perform the canonical renormalization group (RG) prescription in tensor space: write down the tensor RG equation, linearize it around a fixed-point tensor, and diagonalize the resulting linearized RG equation to obtain scaling dimensions. The tensor RG methods have had a great success in producing accurate free energy compared with the conventional real-space RG schemes. However, the above-mentioned canonical procedure has not been implemented for general tensor-network-based RG schemes. We extend the success of the tensor methods further to extraction of scaling dimensions through the canonical RG prescription, without explicitly using the conformal field theory. This approach is benchmarked in the context of the Ising models in 1D and 2D. Based on a pure RG argument, the proposed method has potential applications to 3D systems, where the existing bread-and-butter method is inapplicable.
We employ an adaptation of a strong-disorder renormalization-group technique in order to analyze the ferro-paramagnetic quantum phase transition of Ising chains with aperiodic but deterministic couplings under the action of a transverse field. In the presence of marginal or relevant geometric fluctuations induced by aperiodicity, for which the critical behavior is expected to depart from the Onsager universality class, we derive analytical and asymptotically exact expressions for various critical exponents (including the correlation-length and the magnetization exponents, which are not easily obtainable by other methods), and shed light onto the nature of the ground state structures in the neighborhood of the critical point. The main results obtained by this approach are confirmed by finite-size scaling analyses of numerical calculations based on the free-fermion method.
We revisit perturbative RG analysis in the replicated Landau-Ginzburg description of the Random Field Ising Model near the upper critical dimension 6. Working in a field basis with manifest vicinity to a weakly-coupled Parisi-Sourlas supersymmetric fixed point (Cardy, 1985), we look for interactions which may destabilize the SUSY RG flow and lead to the loss of dimensional reduction. This problem is reduced to studying the anomalous dimensions of leaders -- lowest dimension parts of $S_n$-invariant perturbations in the Cardy basis. Leader operators are classified as non-susy-writable, susy-writable or susy-null depending on their symmetry. Susy-writable leaders are additionally classified as belonging to superprimary multiplets transforming in particular $textrm{OSp}(d | 2)$ representations. We enumerate all leaders up to 6d dimension $Delta = 12$, and compute their perturbative anomalous dimensions (up to two loops). We thus identify two perturbations (with susy-null and non-susy-writable leaders) becoming relevant below a critical dimension $d_c approx 4.2$ - $4.7$. This supports the scenario that the SUSY fixed point exists for all $3 < d leq 6$, but becomes unstable for $d < d_c$.
Efficient sampling of complex high-dimensional probability densities is a central task in computational science. Machine Learning techniques based on autoregressive neural networks have been recently shown to provide good approximations of probability distributions of interest in physics. In this work, we propose a systematic way to remove the intrinsic bias associated with these variational approximations, combining it with Markov-chain Monte Carlo in an automatic scheme to efficiently generate cluster updates, which is particularly useful for models for which no efficient cluster update scheme is known. Our approach is based on symmetry-enforced cluster updates building on the neural-network representation of conditional probabilities. We demonstrate that such finite-cluster updates are crucial to circumvent ergodicity problems associated with global neural updates. We test our method for first- and second-order phase transitions in classical spin systems, proving in particular its viability for critical systems, or in the presence of metastable states.