While tremendous research has revealed that symmetry enriches topological phases of matter, more general principles that protect topological phases have yet to be explored. In this Letter, we elucidate the roles of subspaces in free-fermionic topological phases. A subspace property for Hamiltonians enables us to define new topological invariants. It results in peculiar topological boundary phenomena, i.e., the emergence of an unpaired zero mode or zero-winding skin modes, characterizing subspace-protected topological phases. We establish and demonstrate the bulk-boundary correspondence in subspace-protected topological phases. We further discuss the interplay of the subspace property and internal symmetries. Toward application, we also propose possible platforms possessing the subspace property.

We propose a real-space renormalization group algorithm for accurately coarse-graining two-dimensional tensor networks. The central innovation of our method lies in utilizing variational boundary tensors as a globally optimized environment for the entire system. Based on this optimized environment, we construct renormalization projectors that significantly enhance accuracy. By leveraging the canonical form of tensors, our algorithm maintains the same computational complexity as the original tensor renormalization group (TRG) method, yet achieves higher accuracy than existing approaches that do not incorporate entanglement filtering. Our work offers a practical pathway for extending TRG methods to higher dimensions while keeping computational costs manageable.

Given any symmetry acting on a $d$-dimensional quantum field theory, there is an associated $(d+1)$-dimensional topological field theory known as the Symmetry TFT (SymTFT). The SymTFT is useful for decoupling the universal quantities of quantum field theories, such as their generalized global symmetries and 't Hooft anomalies, from their dynamics. In this work, we explore the SymTFT for theories with Kramers-Wannier-like duality symmetry in both $(1+1)$d and $(3+1)$d quantum field theories. After constructing the SymTFT, we use it to reproduce the non-invertible fusion rules of duality defects, and along the way we generalize the concept of duality defects to \textit{higher} duality defects. We also apply the SymTFT to the problem of distinguishing intrinsically versus non-intrinsically non-invertible duality defects in $(1+1)$d.

We propose a new tensor renormalization group algorithm, Anisotropic Tensor Renormalization Group (ATRG), for lattice models in arbitrary dimensions. The proposed method shares the same versatility with the Higher-Order Tensor Renormalization Group (HOTRG) algorithm, i.e., it preserves the lattice topology after the renormalization. In comparison with HOTRG, both of the computation cost and the memory footprint of our method are drastically reduced, especially in higher dimensions, by renormalizing tensors in an anisotropic way after the singular value decomposition. We demonstrate the ability of ATRG for the square lattice and the simple cubic lattice Ising models. Although the accuracy of the present method degrades when compared with HOTRG of the same bond dimension, the accuracy with fixed computation time is improved greatly due to the drastic reduction of the computation cost.

Multi-controlled Pauli gates are typical high-level qubit operations that appear in the quantum circuits of various quantum algorithms. We find multi-controlled Pauli gate decompositions with smaller CNOT-count or $T$-depth while keeping the currently known minimum $T$-count. For example, for the CCCZ gate, we find decompositions with CNOT-count 7 or $T$-depth 2 while keeping the $T$-count at the currently known minimum of 6. The discovery of these efficient decompositions improves the computational efficiency of many quantum algorithms. What led to this discovery is the systematic procedure for constructing multi-controlled Pauli gate decompositions. This procedure not only deepens our theoretical understanding of quantum gate decomposition but also leads to more efficient decompositions that have yet to be discovered.

It is known that the 't Hooft anomalies of invertible global symmetries can be characterized by an invertible TQFT in one higher dimension. The analogous statement remains to be understood for non-invertible symmetries. In this note we discuss how the linking invariants in a non-invertible TQFT known as the Symmetry TFT (SymTFT) can be used as a diagnostic for 't Hooft anomalies of non-invertible symmetries. When the non-invertible symmetry is non-intrinsically non-invertible, and hence the SymTFT is a Dijkgraaf-Witten model, the linking invariants can be computed explicitly. We illustrate this proposal through the examples of the abelian Higgs model in 2d, as well as adjoint QCD and $\mathcal{N}=4$ super Yang-Mills in 4d. We also comment on how the 't Hooft anomalies of non-invertible symmetries impose new constraints on the dynamics.

We propose a setup to directly measure the anyonic statistical angle on a single edge of a fractional quantum Hall system, without requiring independent knowledge of non-universal parameters. We consider a Laughlin edge state bent into a closed loop geometry, where tunneling processes are controllably induced between the endpoints of the loop. To illustrate the underlying physical mechanism, we compute the time-dependent current generated by the injection of multiple anyons, and show that its behavior exhibits distinctive features governed by the anyonic statistical angle. The measured current reflects quantum interference effects due to the time-resolved braiding of anyons at the junction. To establish experimental relevance, we introduce a protocol where anyons are probabilistically injected upstream of the loop via a quantum point contact (QPC) source. Unlike in Fabry-Perot interferometers, where phase jumps occur spontaneously due to stochastic quasi-particle motion, here the phase jumps are deliberately induced by source injections. These events imprint measurable signatures in the cross-correlation noise, enabling a controlled statistical analysis of the braiding phase. We further show that, by varying the magnetic field while remaining within the same fractional quantum Hall plateau, the statistical angle can be extracted without relying on the knowledge of other non-universal system parameters. Our results provide a minimal and accessible platform for probing anyonic statistics using a single chiral edge.

The $(1+1)$-dimensional two-color lattice QCD is studied with the Grassmann tensor renormalization group. We construct tensor network representations of theories with the staggered fermion and the Wilson fermion and show that Grassmann tensor networks can describe both cases with the same bond dimension. We also propose an efficient initial tensor compression scheme to gauge degrees of freedom. We compute the number density, chiral condensate, and diquark condensate at finite density, employing the staggered fermions. For the theory with Wilson fermion, a critical point in the negative mass region is identified by inspecting the pseudoscalar condensate and the conformal field theory data.

We classify two-dimensional purely chiral conformal field theories which are defined on two-dimensional surfaces equipped with spin structure and have central charge less than or equal to 16, and discuss their duality webs. This result can be used to confirm that the list of non-supersymmetric ten-dimensional heterotic string theories found in the late 1980s is complete and does not miss any exotic example.

We explore non-invertible symmetries in two-dimensional lattice models with subsystem $\mathbb Z_2$ symmetry. We introduce a subsystem $\mathbb Z_2$-gauging procedure, called the subsystem Kramers-Wannier transformation, which generalizes the ordinary Kramers-Wannier transformation. The corresponding duality operators and defects are constructed by gaugings on the whole or half of the Hilbert space. By gauging twice, we derive fusion rules of duality operators and defects, which enriches ordinary Ising fusion rules with subsystem features. Subsystem Kramers-Wannier duality defects are mobile in both spatial directions, unlike the defects of invertible subsystem symmetries. We finally comment on the anomaly of the subsystem Kramers-Wannier duality symmetry, and discuss its subtleties.

Through coarse-graining, tensor network representations of a two-dimensional critical lattice model flow to a universal four-leg tensor, corresponding to a conformal field theory (CFT) fixed-point. We computed explicit elements of the critical fixed-point tensor, which we identify as the CFT four-point function. This allows us to directly extract the operator product expansion coefficients of the CFT from these tensor elements. Combined with the scaling dimensions obtained from the transfer matrix, we determine the complete set of the CFT data from the fixed-point tensor for any critical unitary lattice model.

In this short paper, we argue that the chiral central charge $c_-$ of a (2+1)d topological ordered state is sometimes strongly constrained by 't Hooft anomaly of anti-unitary global symmetry. For example, if a (2+1)d fermionic TQFT has a time reversal anomaly with $T^2=(-1)^F$ labeled as $\nu\in\mathbb{Z}_{16}$, the TQFT must have $c_-=1/4$ mod $1/2$ for odd $\nu$, while $c_-=0$ mod $1/2$ for even $\nu$. This generalizes the fact that the bosonic TQFT with $T$ anomaly in a particular class must carry $c_-=4$ mod $8$ to fermionic cases. We also study such a constraint for fermionic TQFT with $U(1)\times CT$ symmetry, which is regarded as a gapped surface of the topological superconductor in class AIII.

We study phase transition in the ferromagnetic Potts model with invisible states that are added as redundant states by mean-field calculation and Monte Carlo simulation. Invisible states affect the entropy and the free energy, although they do not contribute to the internal energy. The internal energy and the number of degenerated ground states do not change, if invisible states are introduced into the standard Potts model. A second-order phase transition takes place at finite temperature in the standard $q$-state ferromagnetic Potts model on two-dimensional lattice for $q=2,3$, and 4. However, our present model on two-dimensional lattice undergoes a first-order phase transition with spontaneous $q$-fold symmetry breaking ($q=2,3$, and 4) due to entropy effect of invisible states. We believe that our present model is a fundamental model for analysis of a first-order phase transition with spontaneous discrete symmetry breaking.

Superconductivity in the heavy-fermion compound CeCu2Si2 is a prototypical example of Cooper pairs formed by strongly correlated electrons. For more than 30 years, it has been believed to arise from nodal d-wave pairing mediated by a magnetic glue. Here, we report a detailed study of the specific heat and magnetization at low temperatures for a high-quality single crystal. Unexpectedly, the specific-heat measurements exhibit exponential decay with a two-gap feature in its temperature dependence, along with a linear dependence as a function of magnetic field and the absence of oscillations in the field angle, reminiscent of multiband full-gap superconductivity. In addition, we find anomalous behavior at high fields, attributed to a strong Pauli paramagnetic effect. A low quasiparticle density of states at low energies with a multiband Fermi-surface topology would open a new door into electron pairing in CeCu2Si2.

Kohn-Sham density functional theory is the base of modern computational approaches to electronic structures. Their accuracy vitally relies on the exchange-correlation energy functional, which encapsulates electron-electron interaction beyond the classical one. The functional provides a way to obtain the density and energy without solving the many-body equation and can, in principle, be determined to reproduce the exact ones universally. However, the past approaches are dependent on the theoretical development, which limits the possibility of the functional to human's intuition. Here, we demonstrate a systematic way to machine-learn a functional from a database, without complicated assumptions. The density and energy are related with a flexible feed-forward neural network, which is trained to reproduce accurate dataset, and the KS-DFT is solved by taking the functional derivatives with the back-propagation technique. Surprisingly, a trial functional, trained for just a few molecules, has been shown to be applicable to hundreds of molecular systems with comparable accuracy to the standard functionals. Also, by adding the nodes connected to the hidden layers, a non-local term is straightforwardly included to improve accuracy, which has been hitherto impractically difficult. Utilizing the strategy of rapidly advancing machine learning techniques, this novel approach is expected to enrich the DFT framework by constructing a functional just from a database for materials conventionally difficult to calculate accurately.

The open source ALPS (Algorithms and Libraries for Physics Simulations) project provides a collection of physics libraries and applications, with a focus on simulations of lattice models and strongly correlated electron systems. The libraries provide a convenient set of well-documented and reusable components for developing condensed matter physics simulation codes, and the applications strive to make commonly used and proven computational algorithms available to a non-expert community. In this paper we present an update of the core ALPS libraries. We present in particular new Monte Carlo libraries and new Green's function libraries.

The local Z_N quantized Berry phase for the SU(N) antiferromagnetic Heisenberg spin model is formulated. This quantity, which is a generalization of the local Z_2 Berry phase for SU(2) symmetry, has a direct correspondence to the number of singlet pairs spanning on a particular bond, and is effective as a tool to characterize and classify the various symmetry protected topological phases of one-dimensional SU(N) spin systems. We extend the path-integral quantum Monte Carlo method for the Z_2 Berry phase in order to calculate the ZN Berry phase numerically. We demonstrate our method by calculating the Z_4 Berry phase for the bond-alternating SU(4) antiferromagnetic Heisenberg chain, which is represented by the Young diagram of four columns.

We present an \textit{ab initio} framework to calculate anharmonic phonon frequency and phonon lifetime that is applicable to severely anharmonic systems. We employ self-consistent phonon (SCPH) theory with microscopic anharmonic force constants, which are extracted from density-functional calculations using the least absolute shrinkage and selection operator technique. We apply the method to the high-temperature phase of SrTiO$_{3}$ and obtain well-defined phonon quasiparticles that are free from imaginary frequencies. Here we show that the anharmonic phonon frequency of the antiferrodistortive mode depends significantly on the system size near the critical temperature of the cubic-to-tetragonal phase transition. By applying perturbation theory to the SCPH result, phonon lifetimes are calculated for cubic SrTiO$_{3}$, which are then employed to predict lattice thermal conductivity using the Boltzmann transport equation within the relaxation-time approximation. The presented methodology is efficient and accurate, paving the way toward a reliable description of thermodynamic, dynamic, and transport properties of systems with severe anharmonicity, including thermoelectric, ferroelectric, and superconducting materials.

X-ray-induced carrier dynamics in silicon and gallium arsenide were investigated through intensity variations of transmitted terahertz (THz) pulses in the pico to microsecond time scale with X-ray free-electron laser and synchrotron radiation. We observed a steep reduction in THz transmission with a picosecond scale due to the X-ray-induced carrier generation, followed by a recovery on a nano to microsecond scale caused by the recombination of carriers. The rapid response in the former process is applicable to a direct determination of temporal overlap between THz and X-ray pulses for THz pump-X-ray probe experiments with an accuracy of a few picoseconds.

A pure quantum state can fully describe thermal equilibrium as long as one focuses on local observables. Thermodynamic entropy can also be recovered as the entanglement entropy of small subsystems. When the size of the subsystem increases, however, quantum correlations break the correspondence and cause a correction to this simple volume-law. To elucidate the size dependence of the entanglement entropy is of essential importance in linking quantum physics with thermodynamics, and in addressing recent experiments in ultra-cold atoms. Here we derive an analytic formula of the entanglement entropy for a class of pure states called cTPQ states representing thermal equilibrium. We further find that our formula applies universally to any sufficiently scrambled pure states representing thermal equilibrium, i.e., general energy eigenstates of non-integrable models and states after quantum quenches. Our universal formula can be exploited as a diagnostic of chaotic systems; we can distinguish integrable models from chaotic ones and detect many-body localization with high accuracy.