-
MCGAN: Enhancing GAN Training with Regression-Based Generator Loss
Authors:
Baoren Xiao,
Hao Ni,
Weixin Yang
Abstract:
Generative adversarial networks (GANs) have emerged as a powerful tool for generating high-fidelity data. However, the main bottleneck of existing approaches is the lack of supervision on the generator training, which often results in undamped oscillation and unsatisfactory performance. To address this issue, we propose an algorithm called Monte Carlo GAN (MCGAN). This approach, utilizing an innov…
▽ More
Generative adversarial networks (GANs) have emerged as a powerful tool for generating high-fidelity data. However, the main bottleneck of existing approaches is the lack of supervision on the generator training, which often results in undamped oscillation and unsatisfactory performance. To address this issue, we propose an algorithm called Monte Carlo GAN (MCGAN). This approach, utilizing an innovative generative loss function, termly the regression loss, reformulates the generator training as a regression task and enables the generator training by minimizing the mean squared error between the discriminator's output of real data and the expected discriminator of fake data. We demonstrate the desirable analytic properties of the regression loss, including discriminability and optimality, and show that our method requires a weaker condition on the discriminator for effective generator training. These properties justify the strength of this approach to improve the training stability while retaining the optimality of GAN by leveraging strong supervision of the regression loss. Numerical results on CIFAR-10 and CIFAR-100 datasets demonstrate that the proposed MCGAN significantly and consistently improves the existing state-of-the-art GAN models in terms of quality, accuracy, training stability, and learned latent space. Furthermore, the proposed algorithm exhibits great flexibility for integrating with a variety of backbone models to generate spatial images, temporal time-series, and spatio-temporal video data.
△ Less
Submitted 27 May, 2024;
originally announced May 2024.
-
High Rank Path Development: an approach of learning the filtration of stochastic processes
Authors:
Jiajie Tao,
Hao Ni,
Chong Liu
Abstract:
Since the weak convergence for stochastic processes does not account for the growth of information over time which is represented by the underlying filtration, a slightly erroneous stochastic model in weak topology may cause huge loss in multi-periods decision making problems. To address such discontinuities Aldous introduced the extended weak convergence, which can fully characterise all essentia…
▽ More
Since the weak convergence for stochastic processes does not account for the growth of information over time which is represented by the underlying filtration, a slightly erroneous stochastic model in weak topology may cause huge loss in multi-periods decision making problems. To address such discontinuities Aldous introduced the extended weak convergence, which can fully characterise all essential properties, including the filtration, of stochastic processes; however was considered to be hard to find efficient numerical implementations. In this paper, we introduce a novel metric called High Rank PCF Distance (HRPCFD) for extended weak convergence based on the high rank path development method from rough path theory, which also defines the characteristic function for measure-valued processes. We then show that such HRPCFD admits many favourable analytic properties which allows us to design an efficient algorithm for training HRPCFD from data and construct the HRPCF-GAN by using HRPCFD as the discriminator for conditional time series generation. Our numerical experiments on both hypothesis testing and generative modelling validate the out-performance of our approach compared with several state-of-the-art methods, highlighting its potential in broad applications of synthetic time series generation and in addressing classic financial and economic challenges, such as optimal stopping or utility maximisation problems.
△ Less
Submitted 23 May, 2024;
originally announced May 2024.
-
On a generalization of R. Chapman's "evil determinant"
Authors:
Li-Yuan Wang,
Hai-Liang Wu,
He-Xia Ni
Abstract:
Let $p$ be an odd prime and $x$ be an indeterminate. Recently, Z.-W. Sun proposed the following conjecture: $$\det\left[x+\left(\frac{j-i}{p}\right)\right]_{0\le i,j\le \frac{p-1}{2}}=\begin{cases} (\frac{2}{p})pb_px-a_p & \mbox{if}\ p\equiv 1\pmod4, 1 & \mbox{if}\ p\equiv 3\pmod4, \end{cases}$$ where $a_p$ and $b_p$ are rational numbers related to the fundamental unit and class number of the real…
▽ More
Let $p$ be an odd prime and $x$ be an indeterminate. Recently, Z.-W. Sun proposed the following conjecture: $$\det\left[x+\left(\frac{j-i}{p}\right)\right]_{0\le i,j\le \frac{p-1}{2}}=\begin{cases} (\frac{2}{p})pb_px-a_p & \mbox{if}\ p\equiv 1\pmod4, 1 & \mbox{if}\ p\equiv 3\pmod4, \end{cases}$$ where $a_p$ and $b_p$ are rational numbers related to the fundamental unit and class number of the real quadratic field $\mathbb{Q}(\sqrt{p})$. In this paper, we confirm the above conjecture of Sun based on Vsemirnov's decomposition of Chapman's "evil determinant".
△ Less
Submitted 3 May, 2024;
originally announced May 2024.
-
Quantum wave packet transforms with compact frequency support
Authors:
Hongkang Ni,
Lexing Ying
Abstract:
Different kinds of wave packet transforms are widely used for extracting multi-scale structures in signal processing tasks. This paper introduces the quantum circuit implementation of a broad class of wave packets, including Gabor atoms and wavelets, with compact frequency support. Our approach operates in the frequency space, involving reallocation and reshuffling of signals tailored for manipula…
▽ More
Different kinds of wave packet transforms are widely used for extracting multi-scale structures in signal processing tasks. This paper introduces the quantum circuit implementation of a broad class of wave packets, including Gabor atoms and wavelets, with compact frequency support. Our approach operates in the frequency space, involving reallocation and reshuffling of signals tailored for manipulation on quantum computers. The resulting implementation is different from the existing quantum algorithms for spatially compactly supported wavelets and can be readily extended to quantum transforms of other wave packets with compact frequency support.
△ Less
Submitted 3 May, 2024; v1 submitted 1 May, 2024;
originally announced May 2024.
-
On the determination of path signature from its unitary development
Authors:
Siran Li,
Zijiu Lyu,
Hao Ni,
Jiajie Tao
Abstract:
We establish an explicit, constructive approach to determine any element $X$ in the tensor algebra $\mathcal{T}\left(\mathbb{R}^d\right) = \bigoplus_{n=0}^\infty\left(\mathbb{R}^d\right)^{\otimes n}$ from its moment generating function. The only assumption is that $X$ has a nonzero radius of convergence, which relaxes the condition of having an infinite radius of convergence in the literature. The…
▽ More
We establish an explicit, constructive approach to determine any element $X$ in the tensor algebra $\mathcal{T}\left(\mathbb{R}^d\right) = \bigoplus_{n=0}^\infty\left(\mathbb{R}^d\right)^{\otimes n}$ from its moment generating function. The only assumption is that $X$ has a nonzero radius of convergence, which relaxes the condition of having an infinite radius of convergence in the literature. The key building block of our approach is tridiagonal antisymmetric matrices, whose sparsity offers a considerable advantage for dimension reduction in applications. In particular, specialising $X$ to the signature space of bounded $p$-variation paths in $\mathbb{R}^d$ with $1\leq p <2$, we show that the developments of such sparse matrices are sufficient to separate points over the space of signatures, which yields a refined answer to the "moment problem" concerning the signature. Based on the above theoretical investigations, we propose a new distance function for probability measures on the path space, termed as the "Restricted Path Characteristic Function Distance" (RPCFD), and validate its effectiveness via numerical experiments on hypothesis testing for examples of fractional Brownian motions.
△ Less
Submitted 29 April, 2024;
originally announced April 2024.
-
A PDE approach for solving the characteristic function of the generalised signature process
Authors:
Terry Lyons,
Hao Ni,
Jiajie Tao
Abstract:
The signature of a path, as a fundamental object in Rough path theory, serves as a generating function for non-commutative monomials on path space. It transforms the path into a grouplike element in the tensor algebra space, summarising the path faithfully up to a generalised form of re-parameterisation (a negligible equivalence class in this context). Our paper concerns stochastic processes and s…
▽ More
The signature of a path, as a fundamental object in Rough path theory, serves as a generating function for non-commutative monomials on path space. It transforms the path into a grouplike element in the tensor algebra space, summarising the path faithfully up to a generalised form of re-parameterisation (a negligible equivalence class in this context). Our paper concerns stochastic processes and studies the characteristic function of the path signature of the stochastic process. In contrast to the expected signature, it determines the law on the random signatures without any regularity condition. The computation of the characteristic function of the random signature offers potential applications in stochastic analysis and machine learning, where the expected signature plays an important role. In this paper, we focus on a time-homogeneous Itô diffusion process, and adopt a PDE approach to derive the characteristic function of its signature defined at any fixed time horizon. A key ingredient of our approach is the introduction of the generalised-signature process. This lifting enables us to establish the Feynman-Kac-type theorem for the characteristic function of the generalised-signature process by following the martingale approach. Moreover, as an application of our results, we present a novel derivation of the joint characteristic function of Brownian motion coupled with the Lévy area, leveraging the structure theorem of anti-symmetric matrices.
△ Less
Submitted 29 February, 2024; v1 submitted 4 January, 2024;
originally announced January 2024.
-
Quantum Hamiltonian Learning for the Fermi-Hubbard Model
Authors:
Hongkang Ni,
Haoya Li,
Lexing Ying
Abstract:
This work proposes a protocol for Fermionic Hamiltonian learning. For the Hubbard model defined on a bounded-degree graph, the Heisenberg-limited scaling is achieved while allowing for state preparation and measurement errors. To achieve $ε$-accurate estimation for all parameters, only $\tilde{\mathcal{O}}(ε^{-1})$ total evolution time is needed, and the constant factor is independent of the syste…
▽ More
This work proposes a protocol for Fermionic Hamiltonian learning. For the Hubbard model defined on a bounded-degree graph, the Heisenberg-limited scaling is achieved while allowing for state preparation and measurement errors. To achieve $ε$-accurate estimation for all parameters, only $\tilde{\mathcal{O}}(ε^{-1})$ total evolution time is needed, and the constant factor is independent of the system size. Moreover, our method only involves simple one or two-site Fermionic manipulations, which is desirable for experiment implementation.
△ Less
Submitted 1 May, 2024; v1 submitted 28 December, 2023;
originally announced December 2023.
-
Generative Modelling of Lévy Area for High Order SDE Simulation
Authors:
Andraž Jelinčič,
Jiajie Tao,
William F. Turner,
Thomas Cass,
James Foster,
Hao Ni
Abstract:
It is well known that, when numerically simulating solutions to SDEs, achieving a strong convergence rate better than O(\sqrt{h}) (where h is the step size) requires the use of certain iterated integrals of Brownian motion, commonly referred to as its "Lévy areas". However, these stochastic integrals are difficult to simulate due to their non-Gaussian nature and for a d-dimensional Brownian motion…
▽ More
It is well known that, when numerically simulating solutions to SDEs, achieving a strong convergence rate better than O(\sqrt{h}) (where h is the step size) requires the use of certain iterated integrals of Brownian motion, commonly referred to as its "Lévy areas". However, these stochastic integrals are difficult to simulate due to their non-Gaussian nature and for a d-dimensional Brownian motion with d > 2, no fast almost-exact sampling algorithm is known.
In this paper, we propose LévyGAN, a deep-learning-based model for generating approximate samples of Lévy area conditional on a Brownian increment. Due to our "Bridge-flipping" operation, the output samples match all joint and conditional odd moments exactly. Our generator employs a tailored GNN-inspired architecture, which enforces the correct dependency structure between the output distribution and the conditioning variable. Furthermore, we incorporate a mathematically principled characteristic-function based discriminator. Lastly, we introduce a novel training mechanism termed "Chen-training", which circumvents the need for expensive-to-generate training data-sets. This new training procedure is underpinned by our two main theoretical results.
For 4-dimensional Brownian motion, we show that LévyGAN exhibits state-of-the-art performance across several metrics which measure both the joint and marginal distributions. We conclude with a numerical experiment on the log-Heston model, a popular SDE in mathematical finance, demonstrating that high-quality synthetic Lévy area can lead to high order weak convergence and variance reduction when using multilevel Monte Carlo (MLMC).
△ Less
Submitted 4 August, 2023;
originally announced August 2023.
-
Robust iterative method for symmetric quantum signal processing in all parameter regimes
Authors:
Yulong Dong,
Lin Lin,
Hongkang Ni,
Jiasu Wang
Abstract:
This paper addresses the problem of solving nonlinear systems in the context of symmetric quantum signal processing (QSP), a powerful technique for implementing matrix functions on quantum computers. Symmetric QSP focuses on representing target polynomials as products of matrices in SU(2) that possess symmetry properties. We present a novel Newton's method tailored for efficiently solving the nonl…
▽ More
This paper addresses the problem of solving nonlinear systems in the context of symmetric quantum signal processing (QSP), a powerful technique for implementing matrix functions on quantum computers. Symmetric QSP focuses on representing target polynomials as products of matrices in SU(2) that possess symmetry properties. We present a novel Newton's method tailored for efficiently solving the nonlinear system involved in determining the phase factors within the symmetric QSP framework. Our method demonstrates rapid and robust convergence in all parameter regimes, including the challenging scenario with ill-conditioned Jacobian matrices, using standard double precision arithmetic operations. For instance, solving symmetric QSP for a highly oscillatory target function $α\cos(1000 x)$ (polynomial degree $\approx 1433$) takes $6$ iterations to converge to machine precision when $α=0.9$, and the number of iterations only increases to $18$ iterations when $α=1-10^{-9}$ with a highly ill-conditioned Jacobian matrix. Leveraging the matrix product states the structure of symmetric QSP, the computation of the Jacobian matrix incurs a computational cost comparable to a single function evaluation. Moreover, we introduce a reformulation of symmetric QSP using real-number arithmetics, further enhancing the method's efficiency. Extensive numerical tests validate the effectiveness and robustness of our approach, which has been implemented in the QSPPACK software package.
△ Less
Submitted 23 July, 2023;
originally announced July 2023.
-
Heisenberg-limited Hamiltonian learning for interacting bosons
Authors:
Haoya Li,
Yu Tong,
Hongkang Ni,
Tuvia Gefen,
Lexing Ying
Abstract:
We develop a protocol for learning a class of interacting bosonic Hamiltonians from dynamics with Heisenberg-limited scaling. For Hamiltonians with an underlying bounded-degree graph structure, we can learn all parameters with root mean squared error $ε$ using $\mathcal{O}(1/ε)$ total evolution time, which is independent of the system size, in a way that is robust against state-preparation and mea…
▽ More
We develop a protocol for learning a class of interacting bosonic Hamiltonians from dynamics with Heisenberg-limited scaling. For Hamiltonians with an underlying bounded-degree graph structure, we can learn all parameters with root mean squared error $ε$ using $\mathcal{O}(1/ε)$ total evolution time, which is independent of the system size, in a way that is robust against state-preparation and measurement error. In the protocol, we only use bosonic coherent states, beam splitters, phase shifters, and homodyne measurements, which are easy to implement on many experimental platforms. A key technique we develop is to apply random unitaries to enforce symmetry in the effective Hamiltonian, which may be of independent interest.
△ Less
Submitted 10 July, 2023;
originally announced July 2023.
-
A Neural RDE-based model for solving path-dependent PDEs
Authors:
Bowen Fang,
Hao Ni,
Yue Wu
Abstract:
The concept of the path-dependent partial differential equation (PPDE) was first introduced in the context of path-dependent derivatives in financial markets. Its semilinear form was later identified as a non-Markovian backward stochastic differential equation (BSDE). Compared to the classical PDE, the solution of a PPDE involves an infinite-dimensional spatial variable, making it challenging to a…
▽ More
The concept of the path-dependent partial differential equation (PPDE) was first introduced in the context of path-dependent derivatives in financial markets. Its semilinear form was later identified as a non-Markovian backward stochastic differential equation (BSDE). Compared to the classical PDE, the solution of a PPDE involves an infinite-dimensional spatial variable, making it challenging to approximate, if not impossible. In this paper, we propose a neural rough differential equation (NRDE)-based model to learn PPDEs, which effectively encodes the path information through the log-signature feature while capturing the fundamental dynamics. The proposed continuous-time model for the PPDE solution offers the benefits of efficient memory usage and the ability to scale with dimensionality. Several numerical experiments, provided to validate the performance of the proposed model in comparison to the strong baseline in the literature, are used to demonstrate its effectiveness.
△ Less
Submitted 1 June, 2023;
originally announced June 2023.
-
Small mass limit of expected signature for physical Brownian motion
Authors:
Siran Li,
Hao Ni,
Qianyu Zhu
Abstract:
The model of physical Brownian motion describes the dynamics of a Brownian particle experiencing both friction and external magnetic field. It was investigated as a physically meaningful approach to realising the standard ''mathematical'' Brownian motion, via sending the mass $m \to 0^+$ and performing natural scaling. In this paper, we are concerned with the singular limit analysis of a generalis…
▽ More
The model of physical Brownian motion describes the dynamics of a Brownian particle experiencing both friction and external magnetic field. It was investigated as a physically meaningful approach to realising the standard ''mathematical'' Brownian motion, via sending the mass $m \to 0^+$ and performing natural scaling. In this paper, we are concerned with the singular limit analysis of a generalised stochastic differential equation (SDE) model, motivated by and encompassing the physical Brownian motion. We show that the expected signature of the solution $\left\{P_t\right\}_{t \geq 0}$ for the generalised SDE converges in the limit $m \to 0^+$ to a nontrivial tensor, at each degree of tensors and on each compact time interval $[0,T]$. The solution $\left\{P_t\right\}_{t \geq 0}$, viewed as a rough path, generalises the momentum of particle in classical physical Brownian motion. This singular limit is identified through a delicate convergence analysis based on the graded PDE system for the expected signature of Itô diffusion processes. Explicit closed-form solutions exhibiting intriguing combinatorial patterns are obtained when the coefficient matrix $\mathbf{M}$ in our SDE, which generalises the stress tensor in physical Brownian motion, is diagonalisable. Our work appears among the very first endeavours to study the singular limit of expected signature of diffusion processes, in addition to the literature on the analysis of the expected signature of stochastic processes.
△ Less
Submitted 29 April, 2023;
originally announced May 2023.
-
A note on spike localization for line spectrum estimation
Authors:
Haoya Li,
Hongkang Ni,
Lexing Ying
Abstract:
This note considers the problem of approximating the locations of dominant spikes for a probability measure from noisy spectrum measurements under the condition of residue signal, significant noise level, and no minimum spectrum separation. We show that the simple procedure of thresholding the smoothed inverse Fourier transform allows for approximating the spike locations rather accurately.
This note considers the problem of approximating the locations of dominant spikes for a probability measure from noisy spectrum measurements under the condition of residue signal, significant noise level, and no minimum spectrum separation. We show that the simple procedure of thresholding the smoothed inverse Fourier transform allows for approximating the spike locations rather accurately.
△ Less
Submitted 13 March, 2023; v1 submitted 1 March, 2023;
originally announced March 2023.
-
On efficient quantum block encoding of pseudo-differential operators
Authors:
Haoya Li,
Hongkang Ni,
Lexing Ying
Abstract:
Block encoding lies at the core of many existing quantum algorithms. Meanwhile, efficient and explicit block encodings of dense operators are commonly acknowledged as a challenging problem. This paper presents a comprehensive study of the block encoding of a rich family of dense operators: the pseudo-differential operators (PDOs). First, a block encoding scheme for generic PDOs is developed. Then…
▽ More
Block encoding lies at the core of many existing quantum algorithms. Meanwhile, efficient and explicit block encodings of dense operators are commonly acknowledged as a challenging problem. This paper presents a comprehensive study of the block encoding of a rich family of dense operators: the pseudo-differential operators (PDOs). First, a block encoding scheme for generic PDOs is developed. Then we propose a more efficient scheme for PDOs with a separable structure. Finally, we demonstrate an explicit and efficient block encoding algorithm for PDOs with a dimension-wise fully separable structure. Complexity analysis is provided for all block encoding algorithms presented. The application of theoretical results is illustrated with worked examples, including the representation of variable coefficient elliptic operators and the computation of the inverse of elliptic operators without invoking quantum linear system algorithms (QLSAs).
△ Less
Submitted 31 May, 2023; v1 submitted 21 January, 2023;
originally announced January 2023.
-
Infinite quantum signal processing
Authors:
Yulong Dong,
Lin Lin,
Hongkang Ni,
Jiasu Wang
Abstract:
Quantum signal processing (QSP) represents a real scalar polynomial of degree $d$ using a product of unitary matrices of size $2\times 2$, parameterized by $(d+1)$ real numbers called the phase factors. This innovative representation of polynomials has a wide range of applications in quantum computation. When the polynomial of interest is obtained by truncating an infinite polynomial series, a nat…
▽ More
Quantum signal processing (QSP) represents a real scalar polynomial of degree $d$ using a product of unitary matrices of size $2\times 2$, parameterized by $(d+1)$ real numbers called the phase factors. This innovative representation of polynomials has a wide range of applications in quantum computation. When the polynomial of interest is obtained by truncating an infinite polynomial series, a natural question is whether the phase factors have a well defined limit as the degree $d\to \infty$. While the phase factors are generally not unique, we find that there exists a consistent choice of parameterization so that the limit is well defined in the $\ell^1$ space. This generalization of QSP, called the infinite quantum signal processing, can be used to represent a large class of non-polynomial functions. Our analysis reveals a surprising connection between the regularity of the target function and the decay properties of the phase factors. Our analysis also inspires a very simple and efficient algorithm to approximately compute the phase factors in the $\ell^1$ space. The algorithm uses only double precision arithmetic operations, and provably converges when the $\ell^1$ norm of the Chebyshev coefficients of the target function is upper bounded by a constant that is independent of $d$. This is also the first numerically stable algorithm for finding phase factors with provable performance guarantees in the limit $d\to \infty$.
△ Less
Submitted 21 September, 2022;
originally announced September 2022.
-
Some $q$-congruences involving central $q$-binomial coefficients
Authors:
He-Xia Ni
Abstract:
Suppose that $p$ is an odd prime and $m$ is an integer not divisible by $p$. Sun and Tauraso [Adv. in Appl. Math., 45(2010), 125--148] gave $\sum_{k=0}^{n-1}\binom{2k}{k+d}/m^k$ and $\sum_{k=0}^{n-1}\binom{2k}{k+d}/(km^k)$ modulo $p$ for all $d=0,1, \ldots n$ and $n= p^a$, where $a$ is a positive integer. In this paper, we present some $q$-analogues of these congruences in the cases $m=2, 4$ for a…
▽ More
Suppose that $p$ is an odd prime and $m$ is an integer not divisible by $p$. Sun and Tauraso [Adv. in Appl. Math., 45(2010), 125--148] gave $\sum_{k=0}^{n-1}\binom{2k}{k+d}/m^k$ and $\sum_{k=0}^{n-1}\binom{2k}{k+d}/(km^k)$ modulo $p$ for all $d=0,1, \ldots n$ and $n= p^a$, where $a$ is a positive integer. In this paper, we present some $q$-analogues of these congruences in the cases $m=2, 4$ for any positive integer $n$.
△ Less
Submitted 19 October, 2021;
originally announced October 2021.
-
Towards fast weak adversarial training to solve high dimensional parabolic partial differential equations using XNODE-WAN
Authors:
Paul Valsecchi Oliva,
Yue Wu,
Cuiyu He,
Hao Ni
Abstract:
Due to the curse of dimensionality, solving high dimensional parabolic partial differential equations (PDEs) has been a challenging problem for decades. Recently, a weak adversarial network (WAN) proposed in (Y.Zang et al., 2020) offered a flexible and computationally efficient approach to tackle this problem defined on arbitrary domains by leveraging the weak solution. WAN reformulates the PDE pr…
▽ More
Due to the curse of dimensionality, solving high dimensional parabolic partial differential equations (PDEs) has been a challenging problem for decades. Recently, a weak adversarial network (WAN) proposed in (Y.Zang et al., 2020) offered a flexible and computationally efficient approach to tackle this problem defined on arbitrary domains by leveraging the weak solution. WAN reformulates the PDE problem as a generative adversarial network, where the weak solution (primal network) and the test function (adversarial network) are parameterized by the multi-layer deep neural networks (DNNs). However, it is not yet clear whether DNNs are the most effective model for the parabolic PDE solutions as they do not take into account the fundamentally different roles played by time and spatial variables in the solution. To reinforce the difference, we design a novel so-called XNODE model for the primal network, which is built on the neural ODE (NODE) model with additional spatial dependency to incorporate the a priori information of the PDEs and serve as a universal and effective approximation to the solution. The proposed hybrid method (XNODE-WAN), by integrating the XNODE model within the WAN framework, leads to significant improvement in the performance and efficiency of training. Numerical results show that our method can reduce the training time to a fraction of that of the WAN model.
△ Less
Submitted 14 October, 2021;
originally announced October 2021.
-
$q$-Supercongruences from transformation formulas
Authors:
He-Xia Ni,
Li-Yuan Wang,
Hai-Liang Wu
Abstract:
Let $Φ_{n}(q)$ denote the $n$-th cyclotomic polynomial in $q$. Recently, Guo and Schlosser [Constr. Approx. 53 (2021), 155--200] put forward the following conjecture: for an odd integer $n>1$, \begin{align*} &\sum_{k=0}^{n-1}[8k-1]\frac{(q^{-1};q^4)_k^6(q^2;q^2)_{2k}}{(q^4;q^4)_k^6(q^{-1};q^2)_{2k}}q^{8k}\notag\\ &\quad\equiv\begin{cases}0 \pmod{[n]Φ_n(q)^2}, &\text{if }n\equiv 1\pmod{4},\\[5pt] 0…
▽ More
Let $Φ_{n}(q)$ denote the $n$-th cyclotomic polynomial in $q$. Recently, Guo and Schlosser [Constr. Approx. 53 (2021), 155--200] put forward the following conjecture: for an odd integer $n>1$, \begin{align*} &\sum_{k=0}^{n-1}[8k-1]\frac{(q^{-1};q^4)_k^6(q^2;q^2)_{2k}}{(q^4;q^4)_k^6(q^{-1};q^2)_{2k}}q^{8k}\notag\\ &\quad\equiv\begin{cases}0 \pmod{[n]Φ_n(q)^2}, &\text{if }n\equiv 1\pmod{4},\\[5pt] 0 \pmod{[n]},&\text{if }n\equiv 3\pmod{4}. \end{cases} \end{align*} Applying the `creative microscoping' method and several summation and transformation formulas for basic hypergeometric series and the Chinese remainder theorem for coprime polynomials, we confirm the above conjecture, as well as another similar $q$-supercongruence conjectured by Guo and Schlosser.
△ Less
Submitted 24 September, 2021;
originally announced September 2021.
-
A conjecture of Zhi-Wei Sun on determinants over finite fields
Authors:
Hai-Liang Wu,
Yue-Feng She,
He-Xia Ni
Abstract:
In this paper, we study certain determinants over finite fields. Let $\mathbb{F}_q$ be the finite field of $q$ elements and let $a_1,a_2,\cdots,a_{q-1}$ be all nonzero elements of $\mathbb{F}_q$. Let $T_q=\left[\frac{1}{a_i^2-a_ia_j+a_j^2}\right]_{1\le i,j\le q-1}$ be a matrix over $\mathbb{F}_q$. We obtain the explicit value of $\det T_q$. Also, as a consequence of our result, we confirm a conjec…
▽ More
In this paper, we study certain determinants over finite fields. Let $\mathbb{F}_q$ be the finite field of $q$ elements and let $a_1,a_2,\cdots,a_{q-1}$ be all nonzero elements of $\mathbb{F}_q$. Let $T_q=\left[\frac{1}{a_i^2-a_ia_j+a_j^2}\right]_{1\le i,j\le q-1}$ be a matrix over $\mathbb{F}_q$. We obtain the explicit value of $\det T_q$. Also, as a consequence of our result, we confirm a conjecture posed by Zhi-Wei Sun.
△ Less
Submitted 12 January, 2022; v1 submitted 24 August, 2021;
originally announced August 2021.
-
Expected signature of stopped Brownian motion on $d$-dimensional $C^{2, α}$-domains has finite radius of convergence everywhere: $2\leq d \leq 8$
Authors:
Siran Li,
Hao Ni
Abstract:
A fundamental question in rough path theory is whether the expected signature of a geometric rough path completely determines the law of signature. One sufficient condition is that the expected signature has infinite radius of convergence, which is satisfied by various stochastic processes on a fixed time interval, including the Brownian motion. In contrast, for the Brownian motion stopped upon th…
▽ More
A fundamental question in rough path theory is whether the expected signature of a geometric rough path completely determines the law of signature. One sufficient condition is that the expected signature has infinite radius of convergence, which is satisfied by various stochastic processes on a fixed time interval, including the Brownian motion. In contrast, for the Brownian motion stopped upon the first exit time from a bounded domain $Ω$, it is only known that the radius of convergence for the expected signature on sufficiently regular $Ω$ is strictly positive everywhere, and that the radius of convergence is finite at some point when $Ω$ is the $2$-dimensional unit disc ([1]).
In this paper, we prove that on any bounded $C^{2,α}$-domain $Ω\subset \mathbb{R}^d$ with $2\leq d \leq 8$, the expected signature of the stopped Brownian motion has finite radius of convergence everywhere. A key ingredient of our proof is the introduction of a "domain-averaging hyperbolic development" (see Definition 4.1), which allows us to symmetrize the PDE system for the hyperbolic development of expected signature by averaging over rotated domains.
△ Less
Submitted 7 October, 2021; v1 submitted 16 November, 2020;
originally announced November 2020.
-
Two $q$-supercongruences from Watson's transformation
Authors:
He-Xia Ni,
Li-Yuan Wang
Abstract:
Guo and Zudilin [Adv. Math. 346 (2019), 329--358] introduced a new method called `creative microscoping', to prove many $q$-supercongruences in a unified way. In this paper, we apply this method and Watson's ${}_8φ_7$ transformation formula to prove two
$q$-supercongruences, which were recently conjectured by Guo and Schlosser.
Guo and Zudilin [Adv. Math. 346 (2019), 329--358] introduced a new method called `creative microscoping', to prove many $q$-supercongruences in a unified way. In this paper, we apply this method and Watson's ${}_8φ_7$ transformation formula to prove two
$q$-supercongruences, which were recently conjectured by Guo and Schlosser.
△ Less
Submitted 29 September, 2020; v1 submitted 16 September, 2020;
originally announced September 2020.
-
A $q$-Dwork-type generalization of Rodriguez-Villegas' supercongruences
Authors:
He-Xia Ni
Abstract:
Guo and Zudilin [Adv. Math. 346 (2019), 329--358] developed an analytical method, called `creative microscoping', to prove many supercongruences by establishing their $q$-analogues. In this paper, we apply this method to give a $q$-Dwork-type generalization of Rodriguez-Villegas' supercongruences, which was recently conjectured by Guo and Zudilin.
Guo and Zudilin [Adv. Math. 346 (2019), 329--358] developed an analytical method, called `creative microscoping', to prove many supercongruences by establishing their $q$-analogues. In this paper, we apply this method to give a $q$-Dwork-type generalization of Rodriguez-Villegas' supercongruences, which was recently conjectured by Guo and Zudilin.
△ Less
Submitted 19 October, 2021; v1 submitted 6 August, 2020;
originally announced August 2020.
-
Supercongruences for Almkvist--Zudilin sequences
Authors:
Ji-Cai Liu,
He-Xia Ni
Abstract:
In this note, we prove two supercongruences involving Almkvist--Zudilin sequences, which were originally conjectured by Z.-H. Sun.
In this note, we prove two supercongruences involving Almkvist--Zudilin sequences, which were originally conjectured by Z.-H. Sun.
△ Less
Submitted 19 April, 2020; v1 submitted 16 April, 2020;
originally announced April 2020.
-
Some $q$-congruences arising from certain identities
Authors:
Chen Wang,
He-Xia Ni
Abstract:
In this paper, by constructing some identities, we prove some $q$-analogues of some congruences. For example, for any odd integer $n>1$, we show that \begin{gather*} \sum_{k=0}^{n-1} \frac{(q^{-1};q^2)_k}{(q;q)_k} q^k \equiv (-1)^{(n+1)/2} q^{(n^2-1)/4} - (1+q)[n] \pmod{Φ_n(q)^2},\\ \sum_{k=0}^{n-1}\frac{(q^3;q^2)_k}{(q;q)_k} q^k \equiv (-1)^{(n+1)/2} q^{(n^2-9)/4} + \frac{1+q}{q^2}[n]\pmod{Φ_n(q)…
▽ More
In this paper, by constructing some identities, we prove some $q$-analogues of some congruences. For example, for any odd integer $n>1$, we show that \begin{gather*} \sum_{k=0}^{n-1} \frac{(q^{-1};q^2)_k}{(q;q)_k} q^k \equiv (-1)^{(n+1)/2} q^{(n^2-1)/4} - (1+q)[n] \pmod{Φ_n(q)^2},\\ \sum_{k=0}^{n-1}\frac{(q^3;q^2)_k}{(q;q)_k} q^k \equiv (-1)^{(n+1)/2} q^{(n^2-9)/4} + \frac{1+q}{q^2}[n]\pmod{Φ_n(q)^2}, \end{gather*} where the $q$-Pochhanmmer symbol is defined by $(x;q)_0=1$ and $(x;q)_k = (1-x)(1-xq)\cdots(1-xq^{k-1})$ for $k\geq1$, the $q$-integer is defined by $[n]=1+q+\cdots+q^{n-1}$ and $Φ_n(q)$ is the $n$-th cyclotomic polynomial. The $q$-congruences above confirm some recent conjectures of Gu and Guo.
△ Less
Submitted 24 March, 2020;
originally announced March 2020.
-
The expected signature of Brownian motion stopped on the boundary of a circle has finite radius of convergence
Authors:
Horatio Boedihardjo,
Joscha Diehl,
Marc Mezzarobba,
Hao Ni
Abstract:
The expected signature is an analogue of the Laplace transform for rough paths. Chevyrev and Lyons showed that, under certain moment conditions, the expected signature determines the laws of signatures. Lyons and Ni posed the question of whether the expected signature of Brownian motion up to the exit time of a domain satisfies Chevyrev and Lyons' moment condition. We provide the first example whe…
▽ More
The expected signature is an analogue of the Laplace transform for rough paths. Chevyrev and Lyons showed that, under certain moment conditions, the expected signature determines the laws of signatures. Lyons and Ni posed the question of whether the expected signature of Brownian motion up to the exit time of a domain satisfies Chevyrev and Lyons' moment condition. We provide the first example where the answer is negative.
△ Less
Submitted 5 February, 2020; v1 submitted 30 May, 2019;
originally announced May 2019.
-
$ε$-Strong Simulation of Fractional Brownian Motion and Related Stochastic Differential Equations
Authors:
Yi Chen,
Jing Dong,
Hao Ni
Abstract:
Consider the fractional Brownian Motion (fBM) $B^H=\{B^H(t): t \in [0,1] \}$ with Hurst index $H\in (0,1)$. We construct a probability space supporting both $B^H$ and a fully simulatable process $\hat B_ε^H $ such that $$\sup_{t\in [0,1]}|B^H(t)-\hat B_ε^H(t)| \le ε$$ with probability one for any user specified error parameter $ε>0$. When $H>1/2$, we further enhance our error guarantee to the $α$-…
▽ More
Consider the fractional Brownian Motion (fBM) $B^H=\{B^H(t): t \in [0,1] \}$ with Hurst index $H\in (0,1)$. We construct a probability space supporting both $B^H$ and a fully simulatable process $\hat B_ε^H $ such that $$\sup_{t\in [0,1]}|B^H(t)-\hat B_ε^H(t)| \le ε$$ with probability one for any user specified error parameter $ε>0$. When $H>1/2$, we further enhance our error guarantee to the $α$-Hölder norm for any $α\in (1/2,H)$. This enables us to extend our algorithm to the simulation of fBM driven stochastic differential equations $Y=\{Y(t):t \in[0,1]\}$. Under mild regularity conditions on the drift and diffusion coefficients of $Y$, we construct a probability space supporting both $Y$ and a fully simulatable process $\hat Y_ε$ such that $$\sup_{t\in [0,1]}|Y(t)-\hat Y_ε^H(t)| \le ε$$ with probability one. Our algorithms enjoy the tolerance-enforcement feature, i.e., the error bounds can be updated sequentially. Thus, the algorithms can be readily combined with other simulation techniques like multilevel Monte Carlo to estimate expectation of functionals of fBMs efficiently.
△ Less
Submitted 20 February, 2019;
originally announced February 2019.
-
Divisibility of some binomial sums
Authors:
He-Xia Ni,
Hao Pan
Abstract:
With help of $q$-congruence, we prove the divisibility of some binomial sums. For example, for any integers $ρ,n\geq 2$, $$\sum_{k=0}^{n-1}(4k+1) \binom{2k}{k}^ρ\cdot (-4)^{ρ(n-1-k)} \equiv 0\pmod{2^{ρ-2}n\binom{2n}{n}}.$$
With help of $q$-congruence, we prove the divisibility of some binomial sums. For example, for any integers $ρ,n\geq 2$, $$\sum_{k=0}^{n-1}(4k+1) \binom{2k}{k}^ρ\cdot (-4)^{ρ(n-1-k)} \equiv 0\pmod{2^{ρ-2}n\binom{2n}{n}}.$$
△ Less
Submitted 9 August, 2018;
originally announced August 2018.
-
On the almost universality of $\lfloor x^2/a\rfloor+\lfloor y^2/b\rfloor+\lfloor z^2/c\rfloor$
Authors:
Hai-Liang Wu,
He-Xia Ni,
Hao Pan
Abstract:
In 2013, Farhi conjectured that for each $m\geq 3$, every natural number $n$ can be represented as $\lfloor x^2/m\rfloor+\lfloor y^2/m\rfloor+\lfloor z^2/m\rfloor$ with $x,y,z\in\Z$, where $\lfloor\cdot\rfloor$ denotes the floor function. Moreover, in 2015, Sun conjectured that every natural number $n$ can be written as $\lfloor x^2/a\rfloor+\lfloor y^2/b\rfloor+\lfloor z^2/c\rfloor$ with…
▽ More
In 2013, Farhi conjectured that for each $m\geq 3$, every natural number $n$ can be represented as $\lfloor x^2/m\rfloor+\lfloor y^2/m\rfloor+\lfloor z^2/m\rfloor$ with $x,y,z\in\Z$, where $\lfloor\cdot\rfloor$ denotes the floor function. Moreover, in 2015, Sun conjectured that every natural number $n$ can be written as $\lfloor x^2/a\rfloor+\lfloor y^2/b\rfloor+\lfloor z^2/c\rfloor$ with $x,y,z\in\Z$, where $a,b,c$ are integers and $(a,b,c)\neq (1,1,1),(2,2,2)$.
In this paper, with the help of congruence theta functions, we prove that for each $m\geq 3$, Farhi's conjecture is true for every sufficiently large integer $n$. And for $a,b,c\geq 5$ with $a,b,c$ are pairwisely co-prime, we also confirm Sun's conjecture for every sufficiently large integer $n$.
△ Less
Submitted 10 October, 2018; v1 submitted 26 June, 2018;
originally announced June 2018.
-
Signature inversion for monotone paths
Authors:
Jiawei Chang,
Nick Duffield,
Hao Ni,
Weijun Xu
Abstract:
The aim of this article is to provide a simple sampling procedure to reconstruct any monotone path from its signature. For every N, we sample a lattice path of N steps with weights given by the coefficient of the corresponding word in the signature. We show that these weights on lattice paths satisfy the large deviations principle. In particular, this implies that the probability of picking up a "…
▽ More
The aim of this article is to provide a simple sampling procedure to reconstruct any monotone path from its signature. For every N, we sample a lattice path of N steps with weights given by the coefficient of the corresponding word in the signature. We show that these weights on lattice paths satisfy the large deviations principle. In particular, this implies that the probability of picking up a "wrong" path is exponentially small in N. The argument relies on a probabilistic interpretation of the signature for monotone paths.
△ Less
Submitted 13 December, 2016;
originally announced December 2016.
-
Some symmetric $q$-congruences
Authors:
He-Xia Ni,
Hao Pan
Abstract:
We prove some symmetric $q$-congruences.
We prove some symmetric $q$-congruences.
△ Less
Submitted 5 January, 2016;
originally announced January 2016.
-
Uniqueness of signature for simple curves
Authors:
H. Boedihardjo,
H. Ni,
Z. Qian
Abstract:
We propose a topological approach to the problem of determining a curve from its iterated integrals. In particular, we prove that a family of terms in the signature series of a two dimensional closed curve with finite p variation, 1\leq p<2, are in fact moments of its winding number. This relation allows us to prove that the signature series of a class of simple non-smooth curves uniquely determin…
▽ More
We propose a topological approach to the problem of determining a curve from its iterated integrals. In particular, we prove that a family of terms in the signature series of a two dimensional closed curve with finite p variation, 1\leq p<2, are in fact moments of its winding number. This relation allows us to prove that the signature series of a class of simple non-smooth curves uniquely determine the curves. This implies that outside a Chordal SLE_κ null set, where 0<κ\leq 4, the signature series of curves uniquely determine the curves. Our calculations also enable us to express the Fourier transform of the n-point functions of SLE curves in terms of the expected signature of SLE curves. Although the techniques used in this article are deterministic, the results provide a platform for studying SLE curves through the signatures of their sample paths.
△ Less
Submitted 15 July, 2014; v1 submitted 2 April, 2013;
originally announced April 2013.
-
Concentration and exact convergence rates for expected Brownian signatures
Authors:
Hao Ni,
Weijun Xu
Abstract:
The signature of a $d$-dimensional Brownian motion is a sequence of iterated Stratonovich integrals along the Brownian paths, an object taking values in the tensor algebra over $\RR^{d}$. In this note, we derive the exact rate of convergence for the expected signatures of piecewise linear approximations to Brownian motion. The computation is based on the identification of the set of words whose co…
▽ More
The signature of a $d$-dimensional Brownian motion is a sequence of iterated Stratonovich integrals along the Brownian paths, an object taking values in the tensor algebra over $\RR^{d}$. In this note, we derive the exact rate of convergence for the expected signatures of piecewise linear approximations to Brownian motion. The computation is based on the identification of the set of words whose coefficients are of the leading order, and the convergence is concentrated on this subset of words. Moreover, under the choice of projective tensor norm, we give the explicit value of the leading term constant.
△ Less
Submitted 21 November, 2012; v1 submitted 5 August, 2012;
originally announced August 2012.
-
Expected signature of Brownian motion up to the first exit time from a bounded domain
Authors:
Terry Lyons,
Hao Ni
Abstract:
The signature of a path provides a top down description of the path in terms of its effects as a control [Differential Equations Driven by Rough Paths (2007) Springer]. The signature transforms a path into a group-like element in the tensor algebra and is an essential object in rough path theory. The expected signature of a stochastic process plays a similar role to that played by the characterist…
▽ More
The signature of a path provides a top down description of the path in terms of its effects as a control [Differential Equations Driven by Rough Paths (2007) Springer]. The signature transforms a path into a group-like element in the tensor algebra and is an essential object in rough path theory. The expected signature of a stochastic process plays a similar role to that played by the characteristic function of a random variable. In [Chevyrev (2013)], it is proved that under certain boundedness conditions, the expected value of a random signature already determines the law of this random signature. It becomes of great interest to be able to compute examples of expected signatures and obtain the upper bounds for the decay rates of expected signatures. For instance, the computation for Brownian motion on $[0,1]$ leads to the ``cubature on Wiener space'' methodology [Lyons and Victoir, Proc. R. Soc. Lond. Ser. A Math. Phys. Eng. Sci. 460 (2004) 169-198]. In this paper we fix a bounded domain $Γ$ in a Euclidean space $E$ and study the expected signature of a Brownian path starting at $z\inΓ$ and stopped at the first exit time from $Γ$. We denote this tensor series valued function by $Φ_Γ(z)$ and focus on the case $E=\mathbb{R}^d$. We show that $Φ_Γ(z)$ satisfies an elliptic PDE system and a boundary condition. The equations determining $Φ_Γ$ can be recursively solved; by an iterative application of Sobolev estimates we are able, under certain smoothness and boundedness condition of the domain $Γ$, to prove geometric bounds for the terms in $Φ_Γ(z)$. However, there is still a gap and we have not shown that $Φ_Γ(z)$ determines the law of the signature of this stopped Brownian motion even if $Γ$ is a unit ball.
△ Less
Submitted 17 November, 2015; v1 submitted 31 January, 2011;
originally announced January 2011.