• Giant quantum oscillations of magneto-thermal conductivity amounting to two orders of magnitude of the estimation based on the Wiedemann-Franz law have been observed in the prototypical Weyl semimetal TaAs. The characteristic oscillation frequency ($F$ $\approx$ 7 T) agrees well with that confirmed for a small hole-type Fermi pocket enclosing a Weyl node. A comparative analysis of various potential scenarios suggests a significant electron-phonon coupling that strongly modulates the phonon mean free path through Landau quantization of the electronic density of states to be at the heart. Resembling the chiral-anomaly induced positive magneto-electrical conductivity, an increase of the thermal conductivity in parallel magnetic field has also been observed. Our findings pose the question whether these are characteristic also for other recently discovered topological electronic materials, calling for more intensive investigations along this line.
  • In this work, we exponentially speed up thermal simulations of quantum many-body systems in both one- (1D) and two-dimensional (2D) models using matrix product operators (MPOs). Instead of evolving the density operator {$\hat{\rho}=e^{-\beta \hat{H}}$} linearly in inverse temperature $\beta \equiv 1/T$ as conventional Trotter-Suzuki methods do, we cool down the system {by iteratively projecting the MPO representation of $\hat{\rho}$ to itself, i.e., doubling $\beta$} in the process of imaginary time evolution. This exponential tensor renormalization group (XTRG) scheme, compared to linear evolution schemes, reaches low temperatures much faster, and thus not only saves computational time but also merits better accuracy due to significantly fewer truncation steps. For similar reasons, we also find that the series expansion thermal tensor network (SETTN) approach benefits in both efficiency and precision, from the logarithmic temperature scale setup. For both thermal algorithms, XTRG as well as SETTN, we fully implement non-Abelian and Abelian symmetries to greatly enhance their numerical performance. We employ these cutting-edge techniques for finite temperature simulations to explore low-temperature thermal states of both 1D and 2D Heisenberg models. The entanglement properties, as well as the renormalization group flow of entanglement spectra in MPOs, are discussed, where logarithmic entropies ($\sim \ln{\beta}$) are shown in both spin chains and square lattice models with gapless tower of states.
  • Monodisperse rod-like colloidal particles are known for spontaneously forming both nematic and smectic liquid crystal phases, but their self-assembly was typically exploited from the fundamental soft condensed matter physics perspective. Here we demonstrate that depletion interactions, driven by non-adsorbing polymers like dextran and surfactants, can be used to enrich self-organization of photon-upconversion nanorods into orientationally ordered nematic and smectic-like membrane colloidal superstructures. We study thermodynamic phase diagrams and demonstrate polarization-dependent photon upconversion exhibited by the ensuing composites, which arises from the superposition of unique properties of the solid nanostructures and the long-range ordering enabled by liquid crystalline self-organization. Finally, we discuss how our method of utilizing self-assembly due to the steric and electrostatic interactions, along with attractive depletion forces, can enable technological uses of lyotropic colloidal liquid crystals and mesostructured composite materials enabled by them, even when they are formed by anisotropic nanoparticles with relatively small aspect ratios.
  • Based on their formation mechanisms, Dirac points in three-dimensional systems can be classified as accidental or essential. The former can be further distinguished into type-I and type-II, depending on whether the Dirac cone spectrum is completely tipped over along certain direction. Here, we predict the coexistence of all three kinds of Dirac points in the low-energy band structure of CaAgBi-family materials with a stuffed Wurtzite structure. Two pairs of accidental Dirac points reside on the rotational axis, with one pair being type-I and the other pair type-II; while another essential Dirac point is pinned at the high symmetry point on the Brillouin zone boundary. Due to broken inversion symmetry, the band degeneracy around accidental Dirac points is completely lifted except along the rotational axis, which may enable the splitting of chiral carriers at a ballistic p-n junction with a double negative refraction effect. We clarify their symmetry protections, and find both the Dirac-cone and Fermi arc topological surface states.
  • In this work we propose a series-expansion thermal tensor network (SETTN) approach for efficient simulations of quantum lattice models. This continuous-time SETTN method is based on the numerically exact Taylor series expansion of equilibrium density operator $e^{-\beta H}$ (with $H$ the total Hamiltonian and $\beta$ the imaginary time), and is thus Trotter-error free. We discover, through simulating XXZ spin chain and square-lattice quantum Ising models, that not only the Hamiltonian $H$, but also its powers $H^n$, can be efficiently expressed as matrix product operators, which enables us to calculate with high precision the equilibrium and dynamical properties of quantum lattice models at finite temperatures. Our SETTN method provides an alternative to conventional Trotter-Suzuki renormalization group (RG) approaches, and achieves an unprecedented standard of thermal RG simulations in terms of accuracy and flexibility.
  • Low-dimensional quantum magnets, due to the existence of abundant exotic quantum phases therein and experimental feasibilities in laboratories, continues intriguing people in condensed matter physics. In this work, a comprehensive study of Cu(NO$_3$)$_2$ $\cdot$ 2.5H$_2$O (copper nitrate hemipentahydrate, CN), a spin chain material, is performed with multi-technique approach including thermal tensor network (TTN) simulations, first-principles calculations, as well as magnetization measurements in experiments. Employing a cutting-edge TTN method developed in the present work, we determine the couplings $J=5.13$ K, $\alpha=0.23(1)$ and Land\'e factors $g_{\parallel}=2.31$, $g_{\perp}=2.14$ in an alternating Heisenberg antiferromagnetic chain model, with which one can fit strikingly well the magnetothermodynamic properties. Part of the fitted experimental data are measured on the single-crystal CN specimens synthesized by us. Based on first-principles calculations, we reveal explicitly the spin chain scenario in CN by displaying the calculated electron density distributions, from which the distinct superexchange paths are visualized. On top of that, we investigated the magnetocaloric effect (MCE) in CN by calculating its isentropes and magnetic Gr\"ueisen parameter (GP). Prominent quantum-criticality-enhanced MCE was uncovered, the TTN simulations are in good agreements with measured isentropic lines in the sub-Kelvin region. We propose that CN is potentially a very promising quantum critical coolant, due to the remarkably enhanced MCE near both critical fields of moderate strengths as 2.87 and 4.08 T, respectively.
  • In this work, we investigate the classical loop models doped with monomers and dimers on a square lattice, whose partition function can be expressed as a tensor network (TN). In the thermodynamic limit, we use the boundary matrix product state technique to contract the partition function TN, and determine the thermodynamic properties with high accuracy. In this monomer-dimer-loop model, we find a second-order phase transition between a trivial monomer-condensation and a loop-condensation (LC) phases, which can not be distinguished by any local order parameter, while nevertheless the two phases have distinct topological properties. In the LC phase, we find two degenerate dominating eigenvalues in the transfer-matrix spectrum, as well as a non-vanishing (nonlocal) string order parameter, both of which identify the \textit{topological ergodicity breaking} in the LC phase and can serve as the order parameter for detecting the phase transitions.
  • Using the tensor network approach, we investigate the monomer-dimer models on a checkerboard lattice, in which there are interactions (with strength $\nu$) between the parallel dimers on half of the plaquettes. For the fully packed interacting dimer model, we observe a Kosterlitz-Thouless (KT) transition between the lowtemperature symmetry breaking and the high-temperature critical phases; for the doped monomer-dimer casewith finite chemical potential $\mu$, we also find an order-disorder phase transition which is of second order instead. We use the boundary matrix product state approach to detect the KT and second-order phase transitions and obtain the phase diagrams $\nu-T$ and $\mu-T$ . Moreover, for the noninteracting monomer-dimer model (setting $\mu = \nu = 0$), we get an extraordinarily accurate determination of the free energy per site (negative of the monomer-dimer constant $h_2$) as $f=-0.662\, 798\, 972\, 833\, 746$ with the dimer density $n=0.638\, 123\, 109\, 228\, 547$, both of 15 correct digits.
  • Long-distance collimation of fast electron beams generated by laser-metallic-wire targets has been observed in recent experiments, while the mechanism behind this phenomenon remains unclear. In this work, we investigate in detail the laser-wire interaction processes with a simplified model and Classical Trajectory Monte Carlo simulations, and demonstrate the significance of the self magnetic fields of the beams in the long-distance collimation. Good agreements of simulated image plate patterns with various experiments and detailed analysis of electron trajectories show that the self magnetic fields provide restoring force that is critical for the beam collimation. By studying the wire-length dependence of beam divergence in certain experiments, we clarify that the role of the metallic wire is to balance the space-charge effect and thus maintain the collimation.
  • The possibility of spontaneous magnetization due to the "asymmetry in mass" of charge carriers in a system is investigated. Analysis shows that when the masses of positive and negative charge carriers are identical, no magnetization is predicted. However, if the masses of two species are different, spontaneous magnetic field would appear, either due to the equipartition of magnetic energy or due to fluctuations together with a feedback mechanism. The conditions for magnetization to occur are also obtained, in the form of n-T phase diagram. The theory proposed here, if confirmed by future observations and/or experiments, would provide a new insight on the origin of magnetic fields in the universe.
  • Based on the Dirac equations in the two-dimensional $\pi-$ flux model, we study the interaction effects both in nontrivial gapped and gapless Dirac equations with numerical exact diagonalization method. In the presence of the nearest and next nearest neighbor interactions: for nontrivial gapped Dirac equation, the topological phase is robust and persists in a finite region of the phase diagram; while for gapless Dirac equation, charge-density-wave and stripe phases are identified and the phase diagram in $(V_1, V_2)$ plane is obtained. When the next-next-nearest neighbor interaction is further included to gapless Dirac equation, the topological phase expected in the mean-field theory is absent. Our results are related to the possibility of dynamically generating topological phase from the electronic correlations.