Source author record

Ke Ye

Ke Ye 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

18works
17topics
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

18 published item(s)

preprint2026arXiv

From Text to Simulation: A Multi-Agent LLM Workflow for Automated Chemical Process Design

Process simulation is a critical cornerstone of chemical engineering design. Current automated chemical design methodologies focus mainly on various representations of process flow diagrams. However, transforming these diagrams into executable simulation flowsheets remains a time-consuming and labor-intensive endeavor, requiring extensive manual parameter configuration within simulation software. In this work, we propose a novel multi-agent workflow that leverages the semantic understanding capabilities of large language models(LLMs) and enables iterative interactions with chemical process simulation software, achieving end-to-end automated simulation from textual process specifications to computationally validated software configurations for design enhancement. Our approach integrates four specialized agents responsible for task understanding, topology generation, parameter configuration, and evaluation analysis, respectively, coupled with Enhanced Monte Carlo Tree Search to accurately interpret semantics and robustly generate configurations. Evaluated on Simona, a large-scale process description dataset, our method achieves a 31.1% improvement in the simulation convergence rate compared to state-of-the-art baselines and reduces the design time by 89. 0% compared to the expert manual design. This work demonstrates the potential of AI-assisted chemical process design, which bridges the gap between conceptual design and practical implementation. Our workflow is applicable to diverse process-oriented industries, including pharmaceuticals, petrochemicals, food processing, and manufacturing, offering a generalizable solution for automated process design.

preprint2022arXiv

Skew-sparse matrix multiplication

Based on the observation that $\mathbb{Q}^{(p-1) \times (p-1)}$ is isomorphic to a quotient skew polynomial ring, we propose a new method for $(p-1)\times (p-1)$ matrix multiplication over $\mathbb{Q}$, where $p$ is a prime number. The main feature of our method is the acceleration for matrix multiplication if the product is skew-sparse. Based on the new method, we design a deterministic algorithm with complexity $O(T^{ω-2} p^2)$, where $T\le p-1$ is a parameter determined by the skew-sparsity of input matrices and $ω$ is the asymptotic exponent of matrix multiplication. Moreover, by introducing randomness, we also propose a probabilistic algorithm with complexity $O^\thicksim(t^{ω-2}p^2+p^2\log\frac{1}ν)$, where $t\le p-1$ is the skew-sparsity of the product and $ν$ is the probability parameter.

preprint2022arXiv

When geometry meets optimization theory: partially orthogonal tensors

Due to the multi-linearity of tensors, most algorithms for tensor optimization problems are designed based on the block coordinate descent method. Such algorithms are widely employed by practitioners for their implementability and effectiveness. However, these algorithms usually suffer from the lack of theoretical guarantee of global convergence and analysis of convergence rate. In this paper, we propose a block coordinate descent type algorithm for the low rank partially orthogonal tensor approximation problem and analyse its convergence behaviour. To achieve this, we carefully investigate the variety of low rank partially orthogonal tensors and its geometric properties related to the parameter space, which enable us to locate KKT points of the concerned optimization problem. With the aid of these geometric properties, we prove without any assumption that: (1) Our algorithm converges globally to a KKT point; (2) For any given tensor, the algorithm exhibits an overall sublinear convergence with an explicit rate which is sharper than the usual $O(1/k)$ for first order methods in nonconvex optimization; {(3)} For a generic tensor, our algorithm converges $R$-linearly.

preprint2020arXiv

Symmetric Tensor Decompositions On Varieties

This paper discusses the problem of symmetric tensor decomposition on a given variety $X$: decomposing a symmetric tensor into the sum of tensor powers of vectors contained in $X$. In this paper, we first study geometric and algebraic properties of such decomposable tensors, which are crucial to the practical computations of such decompositions. For a given tensor, we also develop a criterion for the existence of a symmetric decomposition on $X$. Secondly and most importantly, we propose a method for computing symmetric tensor decompositions on an arbitrary $X$. As a specific application, Vandermonde decompositions for nonsymmetric tensors can be computed by the proposed algorithm.

preprint2019arXiv

Linear Convergence of an Alternating Polar Decomposition Method for Low Rank Orthogonal Tensor Approximations

Low rank orthogonal tensor approximation (LROTA) is an important problem in tensor computations and their applications. A classical and widely used algorithm is the alternating polar decomposition method (APD). In this article, an improved version iAPD of the classical APD is proposed. For the first time, all the following four fundamental properties are established for iAPD: (i) the algorithm converges globally and the whole sequence converges to a KKT point without any assumption; (ii) it exhibits an overall sublinear convergence with an explicit rate which is sharper than the usual $O(1/k)$ for first order methods in optimization; (iii) more importantly, it converges $R$-linearly for a generic tensor without any assumption; (iv) for almost all LROTA problems, iAPD reduces to APD after finitely many iterations if it converges to a local minimizer.

preprint2017arXiv

Argyres-Douglas matter and S-duality: Part II

We study S-duality of Argyres-Douglas theories obtained by compactification of 6d (2,0) theories of ADE type on a sphere with irregular punctures. The weakly coupled descriptions are given by the degeneration limit of auxiliary Riemann sphere with marked points, among which three punctured sphere represents isolated superconformal theories. We also discuss twisted irregular punctures and their S-duality.

preprint2016arXiv

A 3d-3d appetizer

We test the 3d-3d correspondence for theories that are labelled by Lens spaces. We find a full agreement between the index of the 3d ${\cal N}=2$ "Lens space theory" $T[L(p,1)]$ and the partition function of complex Chern-Simons theory on $L(p,1)$. In particular, for $p=1$, we show how the familiar $S^3$ partition function of Chern-Simons theory arises from the index of a free theory. For large $p$, we find that the index of $T[L(p,1)]$ becomes a constant independent of $p$. In addition, we study $T[L(p,1)]$ on the squashed three-sphere $S^3_b$. This enables us to see clearly, at the level of partition function, to what extent $G_\mathbb{C}$ complex Chern-Simons theory can be thought of as two copies of Chern-Simons theory with compact gauge group $G$.

preprint2016arXiv

Algorithms for structured matrix-vector product of optimal bilinear complexity

We present explicit algorithms for computing structured matrix-vector products that are optimal in the sense of Strassen, i.e., using a provably minimum number of multiplications. These structures include Toeplitz/Hankel/circulant, symmetric, Toeplitz-plus-Hankel, sparse, and multilevel structures. The last category include \textsc{bttb}, \textsc{bhhb}, \textsc{bccb} but also any arbitrarily complicated nested structures built out of other structures.

preprint2016arXiv

Fast structured matrix computations: tensor rank and Cohn--Umans method

We discuss a generalization of the Cohn-Umans method, a potent technique developed for studying the bilinear complexity of matrix multiplication by embedding matrices into an appropriate group algebra. We investigate how the Cohn-Umans method may be used for bilinear operations other than matrix multiplication, with algebras other than group algebras, and we relate it to Strassen's tensor rank approach, the traditional framework for investigating bilinear complexity. To demonstrate the utility of the generalized method, we apply it to find the fastest algorithms for forming structured matrix-vector product, the basic operation underlying iterative algorithms for structured matrices. The structures we study include Toeplitz, Hankel, circulant, symmetric, skew-symmetric, f-circulant, block-Toeplitz-Toeplitz-block, triangular Toeplitz matrices, Toeplitz-plus-Hankel, sparse/banded/triangular. Except for the case of skew-symmetric matrices, for which we have only upper bounds, the algorithms derived using the generalized Cohn-Umans method in all other instances are the fastest possible in the sense of having minimum bilinear complexity. We also apply this framework to a few other bilinear operations including matrix-matrix, commutator, simultaneous matrix products, and briefly discuss the relation between tensor nuclear norm and numerical stability.

preprint2016arXiv

Inverse tensor eigenvalue problem

A tensor $\mathcal T\in \mathbb T(\mathbb C^n,m+1)$, the space of tensors of order $m+1$ and dimension $n$ with complex entries, has $nm^{n-1}$ eigenvalues (counted with algebraic multiplicities). The inverse eigenvalue problem for tensors is a generalization of that for matrices. Namely, given a multiset $S\in \mathbb C^{nm^{n-1}}/\mathfrak{S}(nm^{n-1})$ of total multiplicity $nm^{n-1}$, is there a tensor in $\mathbb T(\mathbb C^n,m+1)$ such that the multiset of eigenvalues of $\mathcal{T}$ is exact $S$? The solvability of the inverse eigenvalue problem for tensors is studied in this paper. With tools from algebraic geometry, it is proved that the necessary and sufficient condition for this inverse problem to be generically solvable is $m=1,\ \text{or }n=2,\ \text{or }(n,m)=(3,2),\ (4,2),\ (3,3)$.

preprint2016arXiv

New classes of matrix decompositions

The idea of decomposing a matrix into a product of structured matrices such as triangular, orthogonal, diagonal matrices is a milestone of numerical computations. In this paper, we describe six new classes of matrix decompositions, extending our work in arXiv:1307.5132. We prove that every $n\times n$ matrix is a product of finitely many bidiagonal, skew symmetric (when n is even), generic, companion matrices and generalized Vandermonde matrices, respectively. We also prove that a generic $n\times n$ centrosymmetric matrix is a product of finitely many symmetric Toeplitz (resp. persymmetric Hankel) matrices. We determine an upper bound of the number of structured matrices needed to decompose a matrix for each case.

preprint2016arXiv

Schubert varieties and distances between subspaces of different dimensions

We resolve a basic problem on subspace distances that often arises in applications: How can the usual Grassmann distance between equidimensional subspaces be extended to subspaces of different dimensions? We show that a natural solution is given by the distance of a point to a Schubert variety within the Grassmannian. This distance reduces to the Grassmann distance when the subspaces are equidimensional and does not depend on any embedding into a larger ambient space. Furthermore, it has a concrete expression involving principal angles, and is efficiently computable in numerically stable ways. Our results are largely independent of the Grassmann distance --- if desired, it may be substituted by any other common distances between subspaces. Our approach depends on a concrete algebraic geometric view of the Grassmannian that parallels the differential geometric perspective that is well-established in applied and computational mathematics.

preprint2014arXiv

Every matrix is a product of Toeplitz matrices

We show that every n-by-n matrix is generically a product of [n/2] + 1 Toeplitz matrices and always a product of at most 2n+5 Toeplitz matrices. The same result holds true if the word "Toeplitz" is replaced by "Hankel", and the generic bound [n/2] + 1 is sharp. We will see that these decompositions into Toeplitz or Hankel factors are unusual: We may not in general replace the subspace of Toeplitz or Hankel matrices by an arbitrary (2n-1)-dimensional subspace of n-by-n matrices. Furthermore such decompositions do not exist if we require the factors to be symmetric Toeplitz, persymmetric Hankel, or circulant matrices, even if we allow an infinite number of factors. Lastly, we discuss how the Toeplitz and Hankel decompositions of a generic matrix may be computed by either (i) solving a system of linear and quadratic equations if the number of factors is required to be [n/2] + 1, or (ii) Gaussian elimination in O(n^3) time if the number of factors is allowed to be 2n.

preprint2014arXiv

Multiplicities of eigenvalues of tensors

We study in this article multiplicities of eigenvalues of tensors. There are two natural multiplicities associated to an eigenvalue $λ$ of a tensor: algebraic multiplicity $\operatorname{am}(λ)$ and geometric multiplicity $\operatorname{gm}(λ)$. The former is the multiplicity of the eigenvalue as a root of the characteristic polynomial, and the latter is the dimension of the eigenvariety (i.e., the set of eigenvectors) corresponding to the eigenvalue. We show that the algebraic multiplicity could change along the orbit of tensors by the orthogonal linear group action, while the geometric multiplicity of the zero eigenvalue is invariant under this action, which is the main difficulty to study their relationships. However, we show that for a generic tensor, every eigenvalue has a unique (up to scaling) eigenvector, and both the algebraic multiplicity and geometric multiplicity are one. In general, we suggest for an $m$-th order $n$-dimensional tensor the relationship \[ \operatorname{am}(λ)\geq \operatorname{gm}(λ)(m-1)^{\operatorname{gm}(λ)-1}. \] We show that it is true for serveral cases, especially when the eigenvariety contains a linear subspace of dimension $\operatorname{gm}(λ)$ in coordinate form. As both multiplicities are invariants under the orthogonal linear group action in the matrix counterpart, this generalizes the classical result for a matrix: the algebraic mutliplicity is not smaller than the geometric multiplicity.

preprint2013arXiv

Probing W+W-gamma Production and Anomalous Quartic Gauge Boson Couplings at the CERN LHC

Triple gauge boson associated production at the LHC serves as an interesting channel to test the robustness of the Standard Model. Any deviation from its SM prediction may indicate possible existence of relevant new physics, e.g., anomalous quartic gauge boson couplings. In this paper, a Monte-Carlo feasibility study of measuring WWA production with pure leptonic decays and probing anomalous quartic gauge-boson (e.g., WWAA) couplings, is presented in detail for the first time, with parton shower and detector simulation effects taken into account. Our results show that at the sqrt{s} = 14 TeV LHC with an integrated luminosity of 100 (30) fb-1, one can reach a significance of 9 (5) sigma to observe the SM WWA production, and can constrain at the 95% CL the anomalous WWAA coupling parameters, e.g., a_{0,c}^W/Λ^2 (see Ref.[15] for their definitions), at 10^{-5} GeV^{-2}, respectively.

preprint2013arXiv

The CERN LHC Sensitivity on measuring WZGamma Production and Anomalous WWZGamma Coupling

In this paper we present for the first time a detailed Monte Carlo study of measuring WZGamma production with pure leptonic decays and probing anomalous quartic gauge-boson WWZGamma couplings at the sqrt(s) = 14 TeV LHC, with parton shower and detector simulation effects taken into account. We find that with an integrated luminosity of 100 fb-1 and proper selection cuts, the Standard Model WZGamma signal significance can be improved to as much as 3 sigma. After reviewing previous parametrization on anomalous WWZGamma couplings (see e.g.an/Lambda^2 or k2m/Lambda^2 as shown in Ref. [17]), we propose a more general parametrization scheme with 4 free inputs leading only to genuine WWZGamma aQGC couplings. Finally, our numerical results show that one can reach constraints at 95% confidence level of -5.7 * 10^{-5} GeV^{-2} < k2m/Lambda^2 < 5.5 * 10^{-5} GeV^{-2} and -2.2 x 10^{-5} GeV^{-2} < an/Lambda^2 < 2.4 * 10^{-5} GeV^{-2}, which are more stringent than LEP's results by three orders of magnitude.

preprint2012arXiv

On the geometry of tensor network states

We answer a question of L. Grasedyck that arose in quantum information theory, showing that the limit of tensors in a space of tensor network states need not be a tensor network state. We also give geometric descriptions of spaces of tensor networks states corresponding to trees and loops. Grasedyck's question has a surprising connection to the area of Geometric Complexity Theory, in that the result is equivalent to the statement that the boundary of the Mulmuley-Sohoni type variety associated to matrix multiplication is strictly larger than the projections and re-labelings of matrix multiplication. Tensor Network States are also related to graphical models in algebraic statistics.

preprint2011arXiv

The Stabilizer Of Immanants

We describe immanants as trivial modules of the symmetric group and show that any homogeneous polynomial of degree n on the space of n by n matrices preserved up to scalar by left and right action by diagonal matrices and conjugation by permutation matrices is a linear combination of immanants. we prove that the identity component of the stabilizer of any immanant (except determinant, permanent) is the expected one. We also prove that for n>5 the stabilizer of the immanant of any non-symmetric partition (except determinant and permanent) is again the expected one.