
Employing the selflearning quantum Monte Carlo algorithm, we investigate the
frustrated transversefield trianglelattice Ising model coupled to a Fermi
surface. Without fermions, the spin degrees of freedom undergoes a secondorder
quantum phase transition between paramagnetic and clockordered phases. This
quantum critical point (QCP) has an emergent U(1) symmetry and thus belongs to
the (2+1)D XY universality class. In the presence of fermions, spin
fluctuations introduce effective interactions among fermions and distort the
bare Fermi surface towards an interacting one with hot spots and Fermi pockets.
Near the QCP, nonFermiliquid behavior are observed at the hot spots, and the
QCP is rendered into a different universality with HertzMillis type exponents.
The detailed properties of this QCP and possibly related experimental systems
are also discussed.

Although topological mechanical metamaterials have been extensively studied
from a theoretical perspective, their experimental characterization has been
lagging. To address this shortcoming, we present a systematic laserassisted
experimental characterization of topological kagome lattices, aimed at
elucidating their inplane phononic and topological characteristics. We
specifically explore the continuum elasticity limit, which is established when
the ideal hinges that appear in the theoretical models are replaced by
ligaments capable of supporting bending deformation, as observed for instance
in realistic physical lattices. We reveal how the zeroenergy floppy edge modes
predicted for ideal configurations morph into finitefrequency phonon modes
that localize at the edges. By probing the lattices with carefully designed
excitation signals, we are able to extract and characterize all the features of
a complex lowfrequency acoustic regime in which bulk modes and topological
edge modes overlap and entangle in response. The experiments provide
unequivocal evidence of the existence of strong asymmetric wave transport
regimes at finite frequencies.

Time domain simulation is the basis of dynamic security assessment for power
systems. Traditionally, numerical integration methods are adopted by simulation
software to solve nonlinear power system differentialalgebraic equations about
any given contingency under a specific operating condition. An alternative
approach promising for online simulation is to offline derive a semianalytical
solution (SAS) and then online evaluate the SAS over consecutive time windows
regarding the operating condition and contingency until obtaining the
simulation result over a desired period. This paper proposes a general
semianalytical approach that derives and evaluates an SAS in the form of power
series in time to approximate the solutions of power system differential
equations. An errorrate upper bound of the SAS is also proposed to guarantee
the reliable use of adaptive time windows for evaluation of the SAS. A dynamic
bus method is proposed to extend the semianalytical approach for solving
general power system DAEs by efficiently linking the SASs for dynamic
components through the numerical solution of the network algebraic equations.
Case studies performed on the New England 39bus system and the Polish 2383bus
system test the performance of the proposed semianalytical approach and
compare to existing methods. The results show that the SAS based approach has
potentials for online simulations.

EinsteinPodolskyRosen (EPR) steering describes a quantum nonlocal
phenomenon in which one party can nonlocally affect the other's state through
local measurements. It reveals an additional concept of quantum nonlocality,
which stands between quantum entanglement and Bell nonlocality. Recently, a
quantum information task named as subchannel discrimination (SD) provides a
necessary and sufficient characterization of EPR steering. The success
probability of SD using steerable states is higher than using any unsteerable
states, even when they are entangled. However, the detailed construction of
such subchannels and the experimental realization of the corresponding task are
still technologically challenging. In this work, we designed a feasible
collection of subchannels for a quantum channel and experimentally demonstrated
the corresponding SD task where the probabilities of correct discrimination are
clearly enhanced by exploiting steerable states. Our results provide a concrete
example to operationally demonstrate EPR steering and shine a new light on the
potential application of EPR steering.

This paper proposes a new semianalytical approach for online timedomain
power system simulation. The approach applies the differential transformation
method (DTM) to the power system differential equation model to offline derive
a semianalytical solution (SAS) having symbolic variables about time, the
initial state and system conditions. When simulation is online needed for a
contingency under the current system condition, the SAS can be evaluated in
real time to generate simulation results. Compared to the Adomian decomposition
method in obtaining a power system SAS, an SAS derived by the DTM adopts a
recursive form to avoid generating and storing its complete symbolic
expression, which makes both derivation and evaluation of the SAS more
efficient especially for multimachine power systems. The optimal order of a
DTMbased SAS is studied for the best time performance of simulation. The paper
also designs a parallel computing strategy for power system simulation using
the DTMbased SAS. Tests on the IEEE 10machine 39bus system demonstrate
significant speedup of simulation using the proposed approach compared with the
RungeKutta method.

One bottleneck of quantum Monte Carlo (QMC) simulation of strongly correlated
electron systems lies at the scaling relation of computational complexity with
respect to the system sizes. For generic lattice model of interacting fermions,
the best methodology at hand still scales with $\beta N^3$ where $\beta$ is the
inverse temperature and $N$ is the system size. Such scaling behavior has
greatly hampered the accessibility of the universal infrared (IR) physics of
many interesting correlated electron models at (2+1)D, let alone (3+1)D. To
reduce the computational complexity, we develop a new QMC method with
inhomogeneous momentumspace mesh, dubbed elective momentum ultrasize quantum
Monte Carlo (EQMC) method. Instead of treating all fermionic excitations on an
equal footing as in conventional QMC methods, by converting the fermion
determinant into the momentum space, our method focuses on fermion modes that
are directly associated with lowenergy (IR) physics in the vicinity of the
socalled hotspots, while other fermion modes irrelevant for universal
properties are ignored. As shown in the manuscript, for any cutoffindependent
quantities, e.g. scaling exponents, this method can achieve the same level of
accuracy with orders of magnitude increase in computational efficiency. We
demonstrate this method with a model of antiferromagnetic itinerant quantum
critical point, realized via coupling itinerant fermions with a frustrated
transversefield Ising model on a triangle lattice. The system size of $48
\times 48 \times 48$ ($L\times L\times\beta$, more than 4 times larger than
previous investigations) are comfortably accessed with EQMC. With much larger
system sizes, the scaling exponents are unveiled with unprecedentedly high
accuracy, and this result sheds new light on the open debate about the nature
and the universality class of itinerant quantum critical points.

Considering increasing distributed energy resources and responsive loads in
smart grid, this paper proposes a stochastic simulation approach for stability
analysis of a power system having stochastic loads. The proposed approach
solves a stochastic, nonlinear differential equation model of the system in an
analytical way by the Adomian decomposition method and generates
semianalytical solutions that express both deterministic and stochastic state
variables explicitly as symbolic variables so as to embed stochastic processes
directly into the solutions for efficient stability analysis with
uncertainties. The proposed approach is tested on the New England 10machine
39bus system with different penetration levels of stochastic loads. The
approach is also benchmarked with a traditional stochastic simulation approach
based on the EulerMaruyama method. The results show that the new approach has
better time performance and a comparable accuracy.

Transmission lines are vital components in power systems. Tripping of
transmission lines caused by overtemperature is a major threat to the security
of system operations, so it is necessary to efficiently simulate line
temperature under both normal operation conditions and foreseen fault
conditions. Existing methods based on thermalsteadystate analyses cannot
reflect transient temperature evolution, and thus cannot provide timing
information needed for taking remedial actions. Moreover, conventional
numerical method requires huge computational efforts and barricades systemwide
analysis. In this regard, this paper derives an approximate analytical solution
of transmissionline temperature evolution enabling efficient analysis on
multiple operation states. Considering the uncertainties in environmental
parameters, the region of overtemperature is constructed in the environmental
parameter space to realize the overtemperature risk assessment in both the
planning stage and realtime operations. A test on a typical conductor model
verifies the accuracy of the approximate analytical solution. Based on the
analytical solution and numerical weather prediction (NWP) data, an efficient
simulation method for temperature evolution of transmission systems under
multiple operation states is proposed. As demonstrated on an NPCC 140bus
system, it achieves over 1000 times of efficiency enhancement, verifying its
potentials in online risk assessment and decision support.

Motivated by the recent numerical simulations for doped $t$$J$ model on the
honeycomb lattice, we study superconductivity of singlet and triplet pairing on
the honeycomb lattice Hubbard model. We show that a superconducting state with
coexisting spinsinglet and spintriplet pairings is induced by the
antiferromagnetic order near half filling. The superconducting state we obtain
has a topological phase transition that separates a topologically trivial state
and a nontrivial state with Chern number two. Possible experimental realization
of such a topological superconductivity is also discussed.

It is well known that closedform analytical solutions for AC power flow
equations do not exist in general. This paper proposes a multidimensional
holomorphic embedding method (MDHEM) to obtain an explicit approximate
analytical AC powerflow solution by finding a physical germ solution and
arbitrarily embedding each power, each load or groups of loads with respective
scales. Based on the MDHEM, the complete approximate analytical solutions to
the power flow equations in the highdimensional space become achievable, since
the voltage vector of each bus can be explicitly expressed by a convergent
multivariate power series of all the loads. Unlike the traditional iterative
methods for power flow calculation and inaccurate sensitivity analysis method
for voltage control, the algebraic variables of a power system in all operating
conditions can be prepared offline and evaluated online by only plugging in the
values of any operating conditions into the scales of the nonlinear
multivariate power series. Case studies implemented on the 4bus test system
and the IEEE 14bus standard system confirm the effectiveness of the proposed
method.

We construct a twodimensional lattice model of fermions coupled to Ising
ferromagnetic critical fluctuations. Using extensive signproblemfree quantum
Monte Carlo simulations, we show that the model realizes a continuous itinerant
quantum phase transition. In comparison with other similar itinerant quantum
critical points (QCPs), our QCP shows much weaker superconductivity tendency
with no superconducting state down to the lowest temperature investigated,
hence making the system a good platform for the exploration of quantum critical
fluctuations. Remarkably, clear signatures of nonFermiliquid behavior in the
fermion propagators are observed at the QCP. The critical fluctuations at the
QCP partially resemble HertzMillisMoriya behavior. However, careful scaling
analysis reveals that the QCP belongs to a different universality class,
deviating from both (2+1)d Ising and HertzMillisMoriya predictions.

It is wellknown that in an mmachine power system where each machine is
represented by a secondorder differential equation, the Jacobian of the system
equation contains (m1) pairs of conjugate eigenvalues and two real
eigenvalues, including at least one zero. This letter proves that under the
uniform damping condition, the dynamics associated with the two real
eigenvalues do not have any impact on the dynamics associated with those
complex eigenvalues. This conclusion is important to justify the use of the
relative motions or centerofinertia (COI) coordinate to analyze the rotor
angle stability in a multimachine power system.

This paper proposes an online steadystate voltage stability assessment
scheme to evaluate the proximity to voltage collapse at each bus of a load
area. Using a noniterative holomorphic embedding method (HEM) with a proposed
physical germ solution, an accurate loading limit at each load bus can be
calculated based on online state estimation on the entire load area and a
measurementbased equivalent for the external system. The HEM employs a power
series to calculate an accurate PowerVoltage (PV) curve at each load bus and
accordingly evaluates the voltage stability margin considering load variations
in the next period. An adaptive twostage Pade approximants method is proposed
to improve the convergence of the power series for accurate determination of
the nose point on the PV curve with moderate computational burden. The
proposed method is illustrated in detail on a 4bus test system and then
demonstrated on a load area of the Northeast Power Coordinating Council (NPCC)
48geneartor, 140bus power system.

We report discovery of a topological Mott insulator in stronglycorrelated
Dirac semimetals. Such an interactiondriven topological state has been
theoretically proposed but not yet observed with unbiased large scale numerical
simulations. In our model, interactions between electrons are mediated by Ising
spins in a transverse field. The results indicate that the topological mass
term is dynamically generated and the resulting quantum phase transition
belongs to the (2+1)D $N=8$ chiral Ising universality class. These conclusions
stem from large scale sign free quantum Monte Carlo simulations.

Weather and environmental factors are verified to have played significant
roles in historical major cascading outages and blackouts. Therefore, in the
simulation and risk assessment of cascading outages in power systems, it is
necessary to consider the weather and environmental effects. This paper
proposes a method for the risk assessment of weatherrelated cascading outages.
Based on the analysis of historical outage records and temperaturedependent
physical outage mechanisms of transmission lines, an outage rate model
considering weather condition and conductor temperature is proposed, and the
analytical form of outage probability of lines are derived. With the
weatherdependent outage model, a twostage risk assessment method based on
Markovian tree (MT) search is proposed, which consists of offline full
assessment, and online efficient update of risk assessment results and
continued MT search using updated NWP data. The test cases on NPCC 140bus test
system model in winter and summer scenarios verify the advantages of the
proposed risk assessment method in both accuracy and efficiency.

The recently proposed noniterative load flow method, called the holomorphic
embedding method, may encounter the precision issue, i.e. nontrivial roundoff
errors caused by the limit of digits used in computation when calculating the
powervoltage (PV) curve for a heavily loaded power system. This letter
proposes a multistage scheme to solve such a precision issue and calculate an
accurate PV curve. The scheme is verified on the New England 39bus power
system and benchmarked with the result from the traditional continuation power
flow method.

EinsteinPodolskyRosen (EPR) steering describes the ability of one party to
remotely affect another's state through local measurements. One of the most
distinguishable properties of EPR steering is its asymmetric aspect. Steering
can work in one direction but fail in the opposite direction. This type of
oneway steering, which is different from the symmetry concepts of entanglement
and Bell nonlocality, has garnered much interest. However, an experimental
demonstration of genuine EPR steering in the simplest scenario, i.e., one that
employs twoqubit systems, is still lacking. In this work, we experimentally
demonstrate oneway EPR steering with multimeasurement settings for a class of
twoqubit states, which are still oneway steerable even with infinite
settings. The steerability is quantified by the steering radius, which
represents a necessary and sufficient steering criterion. The demonstrated
oneway steering in the simplest bipartite quantum system is of fundamental
interest and may provide potential applications in oneway quantum information
tasks.

Since cascading outages are major threats to power system operations, it is
of great significance to reduce the risk of potential cascading outages. In
this paper, a method for reduction of cascading outage risk based on Markovian
tree (MT) search is proposed. Based on the MT expansion of the cascading outage
risk, the risk gradient is computed with a forwardbackward tree search scheme.
The computation of risk gradient is incorporated into the procedure of risk
assessment based on MT search, which is efficient with good convergence. Then
the an optimization model for risk reduction (RR) is formulated by using risk
gradient, which minimizes the cost of control while effectively reduces the
cascading outage risk. Moreover, to overcome the limitation of linearization,
an iterative risk reduction (IRR) algorithm is further developed. Test results
on a 4bus test system and the RTS96 3area test system verify the accuracy of
risk gradient computation and effectiveness of the RR. And the performance of
the IRR is demonstrated on RTS96 3area system and a 410bus USCanada
northeast system model. The results show that the subsequent cascade risk and
the total risk are reduced by 93.6% and 54.5%, respectively.

Weyl semimetal (WSM) is expected to be an ideal spintronic material owing to
its spin currents carried by the bulk and surface states with spinmomentum
locking. A photocurrent generation in noncentrosymmetric WSM was also predicted
owing to its broken inversion symmetry and linear energy dispersion which are
unique to Weyl systems. In our recent measurements, the circular photogalvanic
effect (CPGE) has been demonstrated in WSM of TaAs. The CPGE voltage is
proportional to the helicity of the incident light and reverses its direction
on changing the radiation helicity from left handed to right handed, a
periodical oscillation therefore appears following with the alteration of
optical polarization. We herein attribute the CPGE to the asymmetric optical
excitation of the Weyl cone, which could result in an asymmetric distribution
of photoexcited carriers in momentum space according to an optical spin
selection rule.

Topological quantum computation aims to employ anyonic quasiparticles with
exotic braiding statistics to encode and manipulate quantum information in a
faulttolerant way. Majorana zero modes are experimentally the simplest
realisation of anyons that can nontrivially process quantum information.
However, their braiding evolutions, necessary for realising topological gates,
still remain beyond current technologies. Here we report the experimental
encoding of four Majorana zero modes in an alloptical quantum simulator that
give rise to a faulttolerant qubit. We experimentally simulate their braiding
and demonstrate both the nonAbelian character and the topological nature of
the resulting geometric phase. We realise a full set of topological and
nontopological gates that can arbitrarily rotate the encoded qubit. As an
application, we implement the DeutschJozsa algorithm exclusively by
topological gates. Our experiment indicates the intriguing possibility of the
experimental simulation of Majoranabased quantum computation with scalable
technologies.

Interactiondriven topological phase transitions in Dirac semimetals are
investigated by means of largescale quantum Monte Carlo (QMC) simulations. The
interaction among Dirac fermions is introduced by coupling them to Ising spins
that realize the quantum dynamics of the twodimensional transverse field Ising
model. The ground state phase diagram, in which the tuning parameters are the
transverse field and the coupling between fermion and Ising spins, is
determined. At weak and intermediate coupling, a secondorder Ising quantum
phase transition and a firstorder topological phase transition between two
topologically distinct Dirac semimetals are observed. Interestingly, at the
latter, the Dirac points smear out to form nodal lines in the Brillouin zone,
and collective bosonic fluctuations with SO(4) symmetry are strongly enhanced.
At strong coupling, these two phase boundaries merge into a firstorder
transition.

We present a new model to explain the difference between the transport and
spectroscopy gaps in samarium hexaboride (SmB$_6$), which has been a mystery
for some time. We propose that SmB$_6$ can be modeled as an intrinsic
semiconductor with a depletion length that diverges at cryogenic temperatures.
In this model, we find a selfconsistent solution to Poisson's equation in the
bulk, with boundary conditions based on Fermi energy pinning due to surface
charges. The solution yields band bending in the bulk; this explains the
difference between the two gaps because spectroscopic methods measure the gap
near the surface, while transport measures the average over the bulk. We also
connect the model to transport parameters, including the Hall coefficient and
thermopower, using semiclassical transport theory. The divergence of the
depletion length additionally explains the 1012 K feature in data for these
parameters, demonstrating a crossover from bulk dominated transport above this
temperature to surfacedominated transport below this temperature. We find good
agreement between our model and a collection of transport data from 440 K.
This model can also be generalized to materials with similar band structure.

This paper studies the semianalytic solution (SAS) of a power system's
differentialalgebraic equation. A SAS is a closedform function of symbolic
variables including time, the initial state and the parameters on system
operating conditions, and hence able to directly give trajectories on system
state variables, which are accurate for at least a certain time window. A
twostage SASbased approach for fast transient stability simulation is
proposed, which offline derives the SAS by the Adomian Decomposition Method and
online evaluates the SAS for each of sequential time windows until making up a
desired simulation period. When applied to fault simulation, the new approach
employs numerical integration only for the faulton period to determine the
postdisturbance initial state of the SAS. The paper further analyzes the
maximum length of a time window for a SAS to keep its accuracy, and
accordingly, introduces a divergence indicator for adaptive time windows. The
proposed SASbased new approach is validated on the IEEE 10machine, 39bus
system.

Knots and links are fascinating and intricate topological objects. Their
influence spans from DNA and molecular chemistry to vortices in superfluid
helium, defects in liquid crystals and cosmic strings in the early universe.
Here, we find that knotted structures also exist in a peculiar class of
threedimensional topological insulatorsthe Hopf insulators. In particular,
we demonstrate that the momentumspace spin textures of Hopf insulators are
twisted in a nontrivial way, which implies the presence of various knot and
link structures. We further illustrate that the knots and nontrivial spin
textures can be probed via standard timeofflight images in cold atoms as
preimage contours of spin orientations in stereographic coordinates. The
extracted Hopf invariants, knots, and links are validated to be robust to
typical experimental imperfections. Our work establishes the existence of
knotted structures in Hopf insulators, which may have potential applications in
spintronics and quantum information processing.

Many natural and manmade dynamical systems that are modeled as large
nonlinear multioscillator systems like power systems are hard to analyze. For
such a system, we propose a nonlinear modal decoupling (NMD) approach inversely
constructing as many decoupled nonlinear oscillators as the system oscillation
modes so that individual decoupled oscillators can easily be analyzed to infer
dynamics and stability of the original system. The NMD follows a similar idea
to the normal form except that we eliminate intermodal terms but allow
intramodal terms of desired nonlinearities in decoupled systems, so decoupled
systems can flexibly be shaped into desired forms of nonlinear oscillators. The
NMD is then applied to power systems towards two types of nonlinear
oscillators, i.e. the singlemachineinfinitebus (SMIB) systems and a proposed
nonSMIB oscillator. Numerical studies on a 3machine 9bus system and New
England 10machine 39bus system show that (i) decoupled oscillators keep a
majority of the original system modal nonlinearities and the NMD provides a
bigger validity region than the normal form, and (ii) decoupled nonSMIB
oscillators may keep more authentic dynamics of the original system than
decoupled SMIB systems.