
Giant quantum oscillations of magnetothermal conductivity amounting to two
orders of magnitude of the estimation based on the WiedemannFranz 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 holetype Fermi pocket enclosing a Weyl node. A comparative analysis of
various potential scenarios suggests a significant electronphonon 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
chiralanomaly induced positive magnetoelectrical 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
manybody systems in both one (1D) and twodimensional (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 TrotterSuzuki 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 nonAbelian and Abelian
symmetries to greatly enhance their numerical performance. We employ these
cuttingedge techniques for finite temperature simulations to explore
lowtemperature 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 rodlike colloidal particles are known for spontaneously forming
both nematic and smectic liquid crystal phases, but their selfassembly was
typically exploited from the fundamental soft condensed matter physics
perspective. Here we demonstrate that depletion interactions, driven by
nonadsorbing polymers like dextran and surfactants, can be used to enrich
selforganization of photonupconversion nanorods into orientationally ordered
nematic and smecticlike membrane colloidal superstructures. We study
thermodynamic phase diagrams and demonstrate polarizationdependent photon
upconversion exhibited by the ensuing composites, which arises from the
superposition of unique properties of the solid nanostructures and the
longrange ordering enabled by liquid crystalline selforganization. Finally,
we discuss how our method of utilizing selfassembly 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 threedimensional
systems can be classified as accidental or essential. The former can be further
distinguished into typeI and typeII, 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 lowenergy band
structure of CaAgBifamily materials with a stuffed Wurtzite structure. Two
pairs of accidental Dirac points reside on the rotational axis, with one pair
being typeI and the other pair typeII; 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 pn junction with a double negative
refraction effect. We clarify their symmetry protections, and find both the
Diraccone and Fermi arc topological surface states.

In this work we propose a seriesexpansion thermal tensor network (SETTN)
approach for efficient simulations of quantum lattice models. This
continuoustime 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 Trottererror free. We
discover, through simulating XXZ spin chain and squarelattice 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 TrotterSuzuki renormalization group (RG)
approaches, and achieves an unprecedented standard of thermal RG simulations in
terms of accuracy and flexibility.

Lowdimensional 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 multitechnique
approach including thermal tensor network (TTN) simulations, firstprinciples
calculations, as well as magnetization measurements in experiments. Employing a
cuttingedge 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 singlecrystal CN specimens synthesized by us. Based on
firstprinciples 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 quantumcriticalityenhanced MCE
was uncovered, the TTN simulations are in good agreements with measured
isentropic lines in the subKelvin 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 monomerdimerloop
model, we find a secondorder phase transition between a trivial
monomercondensation and a loopcondensation (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 transfermatrix spectrum, as well as a
nonvanishing (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 monomerdimer 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 KosterlitzThouless (KT) transition
between the lowtemperature symmetry breaking and the hightemperature critical
phases; for the doped monomerdimer casewith finite chemical potential $\mu$,
we also find an orderdisorder phase transition which is of second order
instead. We use the boundary matrix product state approach to detect the KT and
secondorder phase transitions and obtain the phase diagrams $\nuT$ and
$\muT$ . Moreover, for the noninteracting monomerdimer model (setting $\mu =
\nu = 0$), we get an extraordinarily accurate determination of the free energy
per site (negative of the monomerdimer 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.

Longdistance collimation of fast electron beams generated by
lasermetallicwire targets has been observed in recent experiments, while the
mechanism behind this phenomenon remains unclear. In this work, we investigate
in detail the laserwire 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 longdistance 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 wirelength dependence of beam divergence in certain experiments, we
clarify that the role of the metallic wire is to balance the spacecharge
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 nT 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 twodimensional $\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, chargedensitywave and
stripe phases are identified and the phase diagram in $(V_1, V_2)$ plane is
obtained. When the nextnextnearest neighbor interaction is further included
to gapless Dirac equation, the topological phase expected in the meanfield
theory is absent. Our results are related to the possibility of dynamically
generating topological phase from the electronic correlations.