Source author record

Ying-Jer Kao

Ying-Jer Kao appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

26works
18topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

26 published item(s)

preprint2022arXiv

J1-J2 fractal studied by multi-recursion tensor-network method

We generalize a tensor-network algorithm to study thermodynamic properties of self-similar spin lattices constructed on a square-lattice frame with two types of couplings, $J_{1}^{}$ and $J_{2}^{}$, chosen to transform a regular square lattice ($J_{1}^{} = J_{2}^{}$) onto a fractal lattice if decreasing $J_{2}^{}$ to zero (the fractal fully reconstructs when $J_{2}^{} = 0$). We modified the Higher-Order Tensor Renormalization Group (HOTRG) algorithm for this purpose. Single-site measurements are performed by means of so-called impurity tensors. So far, only a single local tensor and uniform extension-contraction relations have been considered in HOTRG. We introduce ten independent local tensors, each being extended and contracted by fifteen different recursion relations. We applied the Ising model to the $J_{1}^{}-J_{2}^{}$ planar fractal whose Hausdorff dimension at $J_{2}^{} = 0$ is $d^{(H)} = \ln 12 / \ln 4 \approx 1.792$. The generalized tensor-network algorithm is applicable to a wide range of fractal patterns and is suitable for models without translational invariance.

preprint2022arXiv

Phase structure of the CP(1) model in the presence of a topological $θ$-term

We numerically study the phase structure of the CP(1) model in the presence of a topological $θ$-term, a regime afflicted by the sign problem for conventional lattice Monte Carlo simulations. Using a bond-weighted tensor renormalization group method, we compute the free energy for inverse couplings ranging from $0\leq β\leq 1.1$ and find a CP-violating, first-order phase transition at $θ=π$. In contrast to previous findings, our numerical results provide no evidence for a critical coupling $β_c<1.1$ above which a second-order phase transition emerges at $θ=π$ and/or the first-order transition line bifurcates at $θ\neqπ$. If such a critical coupling exists, as suggested by Haldane's conjecture, our study indicates that is larger than $β_c>1.1$.

preprint2021arXiv

Two-wire Junction of Inequivalent Tomonaga-Luttinger Liquids

We develop two novel numerical schemes to study the conductance of the two-wire junction of inequivalent Tomonaga-Luttinger Liquids. In the first scheme we use the static current-current correlation function across the junction to extract the linear conductance through a relation that is derived via the bosonization method. In the second scheme we apply a bias and evaluate the time-dependent current across the junction to obtain the current-voltage characteristic. The conductance is then extracted from the small bias result within the linear response regime. Both schemes are based on the infinite size matrix product state to minimize the finite-size effects. Due to the lack of the translational invariance, we focus on a finite-size window containing the junction. For time-independent calculations, we use infinite boundary conditions to evaluate the correlations within the window. For time-dependent calculations, we use the window technique to evaluate the local currents within the window. The numerical results obtained by both schemes show excellent agreement with the analytical predictions.

preprint2021arXiv

Variational Quantum Reinforcement Learning via Evolutionary Optimization

Recent advance in classical reinforcement learning (RL) and quantum computation (QC) points to a promising direction of performing RL on a quantum computer. However, potential applications in quantum RL are limited by the number of qubits available in the modern quantum devices. Here we present two frameworks of deep quantum RL tasks using a gradient-free evolution optimization: First, we apply the amplitude encoding scheme to the Cart-Pole problem; Second, we propose a hybrid framework where the quantum RL agents are equipped with hybrid tensor network-variational quantum circuit (TN-VQC) architecture to handle inputs with dimensions exceeding the number of qubits. This allows us to perform quantum RL on the MiniGrid environment with 147-dimensional inputs. We demonstrate the quantum advantage of parameter saving using the amplitude encoding. The hybrid TN-VQC architecture provides a natural way to perform efficient compression of the input dimension, enabling further quantum RL applications on noisy intermediate-scale quantum devices.

preprint2019arXiv

Emergent Snake Magnetic Domains in Canted Kagome Ice

We study the two-dimensional kagome-ice model derived from a pyrochlore lattice with second- and third-neighbor interactions. The canted moments align along the local $\langle 111 \rangle$ axes of the pyrochlore and respond to both in-plane and out-of-plane external fields. We find that the combination of further-neighbor interactions together with the external fields introduces a rich phase diagram with different spin textures. Close to the phase boundaries, metastable $\textit{"snake"}$ domains emerge with extremely long relaxation time. Our kinetic Monte Carlo analysis of the magnetic-field quench process from saturated state shows unusually slow dynamics. Despite that the interior spins are almost frozen in snake domains, the spins on the edge are free to fluctuate locally, leading to frequent creation and annihilation of monopole-anti-monopole bound states. Once the domains are formed, these excitations are localized and can hardly propagate due to the energy barrier of snakes. The emergence of such snake domains may shed light on the experimental observation of dipolar spin ice under tilted fields, and provide a new strategy to manipulate both spin and charge textures in artificial spin ice.

preprint2016arXiv

Field-Induced Ordering in Dipolar Spin Ice

We present numerical studies of dipolar spin ice in the presence of a magnetic field slightly tilted away from the [111] axis. We find a first-order transition from a kagome ice to a $\mathbf{q}=$X state when the external field is tilted toward the $[11\bar{2}]$ direction. This is consistent with the anomalous critical scattering previously observed in the neutron scattering experiment on the spin ice material $\textrm{Ho}_{2}\textrm{Ti}_{2}\textrm{O}_{7}$ in a tilted field [Nat. Phys. 3, 566 (2007)]. We show that this ordering originates from the antiferromagnetic alignment of spin chains on the kagome planes. The residual entropy of the kagome ice is fully recovered. Our result captures the features observed in the experiments and points to the importance of the dipolar interaction in determining ordered states in the spin ice materials. We place our results in the context of recent susceptibility measurements on $\textrm{Dy}_{2}\textrm{Ti}_{2}\textrm{O}_{7}$, showing two features for a [111] field.

preprint2015arXiv

Variational Monte Carlo simulations using tensor-product projected states

We propose an efficient numerical method, which combines the advantages of recently developed tensor-network based methods and standard trial wave functions, to study the ground state properties of quantum many-body systems. In this approach, we apply a projector in the form of a tensor-product operator to an input wave function, such as a Jastrow-type or Hartree-Fock wave function, and optimize the tensor elements via variational Monte Carlo. The entanglement already contained in the input wave function can considerably reduce the bond dimensions compared to the regular tensor-product state representation. In particular, this allows us to also represent states that do not obey the area law of entanglement entropy. In addition, for fermionic systems, the fermion sign structure can be encoded in the input wave function. We show that the optimized states provide good approximations of the ground-state energy and correlation functions in the cases of two-dimensional bosonic and fermonic systems.

preprint2014arXiv

Inverse Melting of an Electronic Liquid Crystal

Inverse melting refers to the rare thermodynamic phenomenon in which a solid melts into a liquid upon cooling, a transition that can occur only when the ordered (solid) phase has more entropy than the disordered (liquid) phase, and that has so far only been observed in a handful of systems. Here we report the first experimental observation for the inverse melting of an electronic liquid crystalline order in strontium-doped lanthanum nickelate, a compound isostructural with the superconducting cuprates, with a hole doping concentration of 1/3. Using x-ray scattering, we demonstrate that the isotropic charge modulation is driven to nematic order by fluctuating spins and shows an inverse melting transition. Using a phenomenological Landau theory, we show that this inverse melting transition is due to the interlayer coupling between the charge and spin orders. This discovery points to the importance of the interlayer correlations in the system, and provides a new perspective to study the intricate nature of the electronic liquid crystal phases in strongly correlated electronic systems, including possibly the Cu- and Fe-based high-Tc superconductors.

preprint2014arXiv

Long-time dynamics of quantum chains: transfer-matrix renormalization group and entanglement of the maximal eigenvector

By using a different quantum-to-classical mapping from the Trotter-Suzuki decomposition, we identify the entanglement structure of the maximal eigenvectors for the associated quantum transfer matrix. This observation provides a deeper insight into the problem of linear growth of the entanglement entropy in time evolution using conventional methods. Based on this observation, we propose a general method for arbitrary temperatures using the biorthonormal transfer-matrix renormalization group. Our method exhibits a competitive accuracy with a much cheaper computational cost in comparison with two recent proposed methods for long-time dynamics based on a folding algorithm [Phys. Rev. Lett. 102, 240603 (2009)] and a modified time-dependent density-matrix renormalization group [Phys. Rev. Lett. 108, 227206 (2012)].

preprint2014arXiv

Quantum Impurity in Luttinger Liquid: Universal Conductance with Entanglement Renormalization

We study numerically the universal conductance of Luttinger liquids wire with a single impurity via the Muti-scale Entanglement Renormalization Ansatz (MERA). The scale invariant MERA provides an efficient way to extract scaling operators and scaling dimensions for both the bulk and the boundary conformal field theories. By utilizing the key relationship between the conductance tensor and ground-state correlation function, the universal conductance can be evaluated within the framework of the boundary MERA. We construct the boundary MERA to compute the correlation functions and scaling dimensions for the Kane-Fisher fixed points by modeling the single impurity as a junction (weak link) of two interacting wires. We show that the universal behavior of the junction can be easily identified within the MERA and argue that the boundary MERA framework has tremendous potential to classify the fixed points in general multi-wire junctions.

preprint2013arXiv

Detection of symmetry-protected topological phases in one dimension with multiscale entanglement renormalization

Symmetry-protected topological (SPT) phases are short-range entangled quantum phases with symmetry, which have gapped excitations in the bulk and gapless modes at the edge. In this paper, we study the SPT phases in the spin-1 Heisenberg chain with a single-ion anisotropy D, which has a quantum phase transition between a Haldane phase and a large-D phase. Using symmetric multiscale entanglement renormalization ansatz (MERA) tensor networks, we study the nonlocal order parameters for time-reversal and inversion symmetry. For the inversion symmetric MERA, we propose a brick-and-rope representation that gives a geometrical interpretation of inversion symmetric tensors. Finally, we study the symmetric renormalization group (RG) flow of the inversion symmetric string-order parameter, and show that entanglement renormalization with symmetric tensors produces proper behavior of the RG fixed-points.

preprint2013arXiv

Finite-size scaling method for the Berezinskii-Kosterlitz-Thouless transition

We test an improved finite-size scaling method for reliably extracting the critical temperature $T_{\rm BKT}$ of a Berezinskii-Kosterlitz-Thouless (BKT) transition. Using known single-parameter logarithmic corrections to the spin stiffness $ρ_s$ at $T_{\rm BKT}$ in combination with the Kosterlitz-Nelson relation between the transition temperature and the stiffness, $ρ_s(T_{\rm BKT})=2T_{\rm BKT}/π$, we define a size dependent transition temperature $T_{\rm BKT}(L_1,L_2)$ based on a pair of system sizes $L_1,L_2$, e.g., $L_2=2L_1$. We use Monte Carlo data for the standard two-dimensional classical XY model to demonstrate that this quantity is well behaved and can be reliably extrapolated to the thermodynamic limit using the next expected logarithmic correction beyond the ones included in defining $T_{\rm BKT}(L_1,L_2)$. For the Monte Carlo calculations we use GPU (graphical processing unit) computing to obtain high-precision data for $L$ up to 512. We find that the sub-leading logarithmic corrections have significant effects on the extrapolation. Our result $T_{\rm BKT}=0.8935(1)$ is several error bars above the previously best estimates of the transition temperature; $T_{\rm BKT} \approx 0.8929$. If only the leading log-correction is used, the result is, however, consistent with the lower value, suggesting that previous works have underestimated $T_{\rm BKT}$ because of neglect of sub-leading logarithms. Our method is easy to implement in practice and should be applicable to generic BKT transitions.

preprint2013arXiv

Half-Magnetization Plateau of a Dipolar Spin Ice in a [100] Field

We report here numerical results of the low-temperature behavior of a dipolar spin ice in a magnetic field along the [100] direction. Tuning the magnetic field, the system exhibit a half-magnetization plateau at low temperature. This half-polarized phase should correspond to a quantum solid phase in an effective 2D quantum boson model, and the transition from the Coulomb phase with a power-law correlation to this state can be regarded as a superfluid to a quantum solid transition. We discuss possible experimental signatures of this half-polarized state.

preprint2012arXiv

Accurate computation of low-temperature thermodynamics for quantum spin chains

We apply the biorthonormal transfer-matrix renormalization group (BTMRG) [Phys. Rev. E 83, 036702 (2011)] to study low-temperature properties of quantum spin chains. Simulation on isotropic Heisenberg spin-1/2 chain demonstrates that the BTMRG outperforms the conventional transfer-matrix renormalization group (TMRG) by successfully accessing far lower temperature unreachable by conventional TMRG, while retaining the same level of accuracy. The power of the method is further illustrated by the calculation of the low-temperature specific heat for a frustrated spin chain.

preprint2012arXiv

Higgs transition from a magnetic Coulomb liquid to a ferromagnet in Yb_2Ti_2O_7

In a class of frustrated magnets known as spin ice, magnetic monopoles emerge as classical defects and interact via the magnetic Coulomb law. With quantum-mechanical interactions, these magnetic charges are carried by fractionalised bosonic quasi-particles, spinons, which can undergo Bose-Einstein condensation through a first-order transition via the Higgs mechanism. Here, we report evidence of a Higgs transition from a magnetic Coulomb liquid to a ferromagnet in single-crystal Yb2Ti2O7. Polarised neutron-scattering experiments show that the diffuse [111]-rod scattering and pinch-point features which develop on cooling are suddenly suppressed below T_C ~ 0.21 K, where magnetic Bragg peaks and a full depolarisation of the neutron spins are observed with thermal hysteresis, indicating a first-order ferromagnetic transition. Our results are explained on the basis of a quantum spin-ice model, whose high-temperature phase is effectively described as a magnetic Coulomb liquid, while the ground state shows a nearly collinear ferromagnetism with gapped spin excitations.

preprint2012arXiv

High-precision Monte Carlo study of the three-dimensional XY model on GPU

We perform large-scale Monte Carlo simulations of the classical XY model on a three-dimensional $L\times L \times L$ cubic lattice using the graphics processing unit (GPU). By the combination of Metropolis single-spin flip, over-relaxation and parallel-tempering methods, we simulate systems up to L=160. Performing the finite-size scaling analysis, we obtain estimates of the critical exponents for the three-dimensional XY universality class: $α=-0.01293(48)$ and $ν=0.67098(16)$. Our estimate for the correlation-length exponent $ν$, in contrast to previous theoretical estimates, agrees with the most recent experimental estimate $ν_{\rm exp}=0.6709(1)$ at the superfluid transition of $^4$He in a microgravity environment.

preprint2012arXiv

Spin-1/2 J1-J2 Heisenberg antiferromagnet on a square lattice: a plaquette renormalized tensor network study

We apply the plaquette renormalization scheme of tensor network states [Phys. Rev. E, 83, 056703 (2011)] to study the spin-1/2 frustrated Heisenberg J1-J2 model on a L*L square lattice with L = 8, 16 and 32. By treating tensor elements as variational parameters, we obtain the ground states for different J2/J1 values, and investigate the staggered magnetic order parameters, the nearest-neighbor spin-spin correlations and the plaquette order parameters. In addition to the well-known Néel-order and collinear-order at low and high J2/J1, we observe a plaquette order at intermediate J2/J1. A continuous transition between the Néel order and the plaquette order near J2 = 0.40J1 is observed. The collinear order emerges at J2 = 0.62J1 through a first-order phase transition.

preprint2012arXiv

Tuning the disorder of superglasses

We study the interplay of superfluidity, glassy and magnetic orders of hardcore bosons with random frustrating interactions. In the classical limit, this model reduces to a $\pm J$ Edwards-Anderson Ising model with concentration $p$ of the ferromagnetic bonds, which hosts a glassy-ferromagnetic transition at a critical concentration $p_c^{\rm cl}\sim 0.77$ on a 3D cubic lattice. Our quantum Monte Carlo simulation results using the worm algorithm show that quantum fluctuations stabilize the coexistence of superfluidity and glassy order ("superglass"), and shift the (super)glassy-ferromagnetic transition to $p_c> p_c^{\rm cl}$. In contrast, antiferromagnetic order coexists with superfluidity to form a supersolid, and the transition to the glassy phase occurs at a higher $p$.

preprint2011arXiv

Length and temperature dependent crossover of charge transport across molecular junctions

We study the electronic transport in a molecular junction in which each site is coupled to a local phonon bath using the non-equilibrium Green's function method. We observe the length period of the oscillatory conductance in odd-numbered chains depends strongly on the applied bias, and the oscillatory behavior is smeared out for the bias voltage near the phonon energy. In addition, a crossover from tunneling to thermally activated hopping transport as the length of the molecule increases is found for the phonon-free case. In the presence of electron-phonon interaction, hopping transport is dominant and a transition from the thermally suppressed to assisted conduction is observed.

preprint2011arXiv

Plaquette Renormalization Scheme for Tensor Network States

We present a method for contracting a square-lattice tensor network in two dimensions, based on auxiliary tensors accomplishing successive truncations (renormalization) of 8-index tensors for 2 by 2 plaquettes into 4-index tensors. The scheme is variational, and thus the tensors can be optimized by minimizing the energy. Test results for the quantum phase transition of the transverse-field Ising model confirm that even the smallest possible tensors (two values for each tensor index at each renormalization level) produce much better results than the simple product (mean-field) state.

preprint2010arXiv

Quantum order by disorder in a semiclassical spin ice

We study the S>1 nearest-neighbor Heisenberg model with a ferromagnetic interaction J and a large non-collinear <111> easy-axis anisotropy D on a pyrochlore lattice. For a finite D>>|J|, the low-energy physics is described by a < 111 > Ising model with second- and third-neighbor exchange interactions. The extensive degeneracy of the ground state manifold in the infinite anisotropic limit is lifted by the transverse quantum fluctuations, and a q=0 ordered state is selected via the quantum order by disorder machanism, through a first-order phase transition at low temperature.

preprint2010arXiv

Symmetry breaking and criticality in tensor-product states

We discuss variationally optimized matrix-product states for the transverse-field Ising chain, using D*D matrices with small D=2-10. For finite system size N there are energy minimums for symmetric as well as symmetry-broken states, which cross each other at a field value hc(N,D); thus the transition is first-order. A continuous transition develops as N->infinity. The asymptotic critical behavior is then always of mean-field type (the magnetization exponent beta=1/2), but a window of field strengths where true Ising scaling holds (beta=1/8) emerges with increasing D. We also demonstrate asymptotic mean-field behavior for infinite-size two-dimensional tensor-product (iPEPS) states with small tensors.

preprint2010arXiv

Towards A Holographic Model of D-Wave Superconductors

The holographic model for S-wave high T_c superconductors developed by Hartnoll, Herzog and Horowitz is generalized to describe D-wave superconductors. The 3+1 dimensional gravitational theory consists a symmetric, traceless second-rank tensor field and a U(1) gauge field in the background of the AdS black hole. Below T_c the tensor field which carries the U(1) charge undergoes the Higgs mechanism and breaks the U(1) symmetry of the boundary theory spontaneously. The phase transition characterized by the D-wave condensate is second order with the mean field critical exponent beta = 1/2. As expected, the AC conductivity is isotropic below T_c and the system becomes superconducting in the DC limit but has no hard gap.

preprint2009arXiv

Peak-Dip-Hump from Holographic Superconductivity

We study the fermionic spectral function in a holographic superconductor model. At zero temperature, the black hole has zero horizon and hence the entropy of the system is zero after the back reaction of the condensate is taken into account. We find the system exhibits the famous peak-dip-hump lineshape with a sharp low-energy peak followed by a dip then a hump at higher energies. This feature is widely observed in the spectrum of several high-T_c superconductors. We also find a linear relation between the gap in the fermionic spectrum and the condensate, indicating the condensate is formed by fermion pairing.

preprint2007arXiv

A short-loop algorithm for quantum Monte Carlo simulations

We present an algorithmic framework for a variant of the quantum Monte Carlo operator-loop algorithm, where non-local cluster updates are constructed in a way that makes each individual loop smaller. The algorithm is designed to increase simulation efficiency in cases where conventional loops become very large, do not close altogether, or otherwise behave poorly. We demonstrate and characterize some aspects of the short-loop on a square lattice spin-1/2 XXZ model where, remarkably, a significant increase in simulation efficiency is observed in some parameter regimes. The simplicity of the model provides a prototype for the use of short-loops on more complicated quantum systems.

preprint2000arXiv

Commensurate and Incommensurate Structure of the Neutron Cross Section in LaSrCuO and YBaCuO

We study the evolution of the d-wave neutron cross-section with variable frequency ωand fixed T (below and above Tc) in two different cuprate families. The evolution from incommensurate to commensurate to incommensurate peaks is rather generic within an RPA-like scheme. This behavior seems to be in reasonable accord with experiments, and may help distinguish between this and the "stripe" scenario.