• Employing the self-learning quantum Monte Carlo algorithm, we investigate the frustrated transverse-field triangle-lattice Ising model coupled to a Fermi surface. Without fermions, the spin degrees of freedom undergoes a second-order quantum phase transition between paramagnetic and clock-ordered 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, non-Fermi-liquid behavior are observed at the hot spots, and the QCP is rendered into a different universality with Hertz-Millis 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 laser-assisted experimental characterization of topological kagome lattices, aimed at elucidating their in-plane 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 zero-energy floppy edge modes predicted for ideal configurations morph into finite-frequency 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 low-frequency 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 differential-algebraic equations about any given contingency under a specific operating condition. An alternative approach promising for online simulation is to offline derive a semi-analytical 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 semi-analytical 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 error-rate 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 semi-analytical 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 39-bus system and the Polish 2383-bus system test the performance of the proposed semi-analytical approach and compare to existing methods. The results show that the SAS based approach has potentials for online simulations.
  • Einstein-Podolsky-Rosen (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 semi-analytical approach for online time-domain power system simulation. The approach applies the differential transformation method (DTM) to the power system differential equation model to offline derive a semi-analytical 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 multi-machine power systems. The optimal order of a DTM-based SAS is studied for the best time performance of simulation. The paper also designs a parallel computing strategy for power system simulation using the DTM-based SAS. Tests on the IEEE 10-machine 39-bus system demonstrate significant speedup of simulation using the proposed approach compared with the Runge-Kutta 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 momentum-space mesh, dubbed elective momentum ultra-size 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 low-energy (IR) physics in the vicinity of the so-called hot-spots, while other fermion modes irrelevant for universal properties are ignored. As shown in the manuscript, for any cutoff-independent 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 transverse-field 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.
  • For dynamic security assessment considering uncertainties in grid operations, this paper proposes an approach for time-domain simulation of a power system having stochastic loads. The proposed approach solves a stochastic differential equation model of the power system in a semi-analytical way using the Adomian decomposition method. The approach generates semi-analytical solutions expressing both deterministic and stochastic variables explicitly as symbolic variables so as to embed stochastic processes directly into the solutions for efficient simulation and analysis. The proposed approach is tested on the New England 10-machine 39-bus system with different levels of stochastic loads. The approach is also benchmarked with a traditional stochastic simulation approach based on the Euler-Maruyama 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 over-temperature 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 thermal-steady-state 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 system-wide analysis. In this regard, this paper derives an approximate analytical solution of transmission-line temperature evolution enabling efficient analysis on multiple operation states. Considering the uncertainties in environmental parameters, the region of over-temperature is constructed in the environmental parameter space to realize the over-temperature risk assessment in both the planning stage and real-time 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 140-bus 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 honeycomb lattice, we study superconductivity of singlet and triplet pairing on honeycomb lattice Hubbard model. We show that a superconducting state with coexisting spin-singlet and spin-triplet 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 closed-form analytical solutions for AC power flow equations do not exist in general. This paper proposes a multi-dimensional holomorphic embedding method (MDHEM) to obtain an explicit approximate analytical AC power-flow 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 high-dimensional 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 non-linear multivariate power series. Case studies implemented on the 4-bus test system and the IEEE 14-bus standard system confirm the effectiveness of the proposed method.
  • We construct a two-dimensional lattice model of fermions coupled to Ising ferromagnetic critical fluctuations. Using extensive sign-problem-free 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 non-Fermi-liquid behavior in the fermion propagators are observed at the QCP. The critical fluctuations at the QCP partially resemble Hertz-Millis-Moriya behavior. However, careful scaling analysis reveals that the QCP belongs to a different universality class, deviating from both (2+1)d Ising and Hertz-Millis-Moriya predictions.
  • It is well-known that in an m-machine power system where each machine is represented by a second-order differential equation, the Jacobian of the system equation contains (m-1) 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 center-of-inertia (COI) coordinate to analyze the rotor angle stability in a multi-machine power system.
  • This paper proposes an online steady-state voltage stability assessment scheme to evaluate the proximity to voltage collapse at each bus of a load area. Using a non-iterative 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 measurement-based equivalent for the external system. The HEM employs a power series to calculate an accurate Power-Voltage (P-V) curve at each load bus and accordingly evaluates the voltage stability margin considering load variations in the next period. An adaptive two-stage Pade approximants method is proposed to improve the convergence of the power series for accurate determination of the nose point on the P-V curve with moderate computational burden. The proposed method is illustrated in detail on a 4-bus test system and then demonstrated on a load area of the Northeast Power Coordinating Council (NPCC) 48-geneartor, 140-bus power system.
  • We report discovery of a topological Mott insulator in strongly-correlated Dirac semimetals. Such an interaction-driven 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 weather-related cascading outages. Based on the analysis of historical outage records and temperature-dependent 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 weather-dependent outage model, a two-stage 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 140-bus 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 non-iterative load flow method, called the holomorphic embedding method, may encounter the precision issue, i.e. nontrivial round-off errors caused by the limit of digits used in computation when calculating the power-voltage (P-V) curve for a heavily loaded power system. This letter proposes a multi-stage scheme to solve such a precision issue and calculate an accurate P-V curve. The scheme is verified on the New Eng-land 39-bus power system and benchmarked with the result from the traditional continuation power flow method.
  • Einstein-Podolsky-Rosen (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 one-way 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 two-qubit systems, is still lacking. In this work, we experimentally demonstrate one-way EPR steering with multimeasurement settings for a class of two-qubit states, which are still one-way steerable even with infinite settings. The steerability is quantified by the steering radius, which represents a necessary and sufficient steering criterion. The demonstrated one-way steering in the simplest bipartite quantum system is of fundamental interest and may provide potential applications in one-way 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 forward-backward 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 4-bus test system and the RTS-96 3-area test system verify the accuracy of risk gradient computation and effectiveness of the RR. And the performance of the IRR is demonstrated on RTS-96 3-area system and a 410-bus US-Canada 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 spin-momentum 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 fault-tolerant way. Majorana zero modes are experimentally the simplest realisation of anyons that can non-trivially 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 all-optical quantum simulator that give rise to a fault-tolerant qubit. We experimentally simulate their braiding and demonstrate both the non-Abelian character and the topological nature of the resulting geometric phase. We realise a full set of topological and non-topological gates that can arbitrarily rotate the encoded qubit. As an application, we implement the Deutsch-Jozsa algorithm exclusively by topological gates. Our experiment indicates the intriguing possibility of the experimental simulation of Majorana-based quantum computation with scalable technologies.
  • Interaction-driven topological phase transitions in Dirac semimetals are investigated by means of large-scale 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 two-dimensional 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 second-order Ising quantum phase transition and a first-order 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 first-order 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 self-consistent 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 10-12 K feature in data for these parameters, demonstrating a crossover from bulk dominated transport above this temperature to surface-dominated transport below this temperature. We find good agreement between our model and a collection of transport data from 4-40 K. This model can also be generalized to materials with similar band structure.
  • This paper studies the semi-analytic solution (SAS) of a power system's differential-algebraic equation. A SAS is a closed-form 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 two-stage SAS-based 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 fault-on period to determine the post-disturbance 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 SAS-based new approach is validated on the IEEE 10-machine, 39-bus 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 three-dimensional topological insulators---the Hopf insulators. In particular, we demonstrate that the momentum-space 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 time-of-flight 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 multi-oscillator 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 inter-modal terms but allow intra-modal 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 single-machine-infinite-bus (SMIB) systems and a proposed non-SMIB oscillator. Numerical studies on a 3-machine 9-bus system and New England 10-machine 39-bus 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 non-SMIB oscillators may keep more authentic dynamics of the original system than decoupled SMIB systems.