-
The Landau--Lifshitz--Bloch equation: Unique existence and finite element approximation
Authors:
Kim-Ngan Le,
Agus L. Soenjaya,
Thanh Tran
Abstract:
The Landau--Lifshitz--Bloch equation (LLBE) describes the evolution of magnetic spin field in a ferromagnet at high temperatures. We consider a viscous (pseudo-parabolic) regularisation of the LLBE for temperatures higher than the Curie temperature, which we call the $ε$-LLBE. Variants of the $ε$-LLBE are applicable to model pattern formation, phase transition, and heat conduction for non-simple m…
▽ More
The Landau--Lifshitz--Bloch equation (LLBE) describes the evolution of magnetic spin field in a ferromagnet at high temperatures. We consider a viscous (pseudo-parabolic) regularisation of the LLBE for temperatures higher than the Curie temperature, which we call the $ε$-LLBE. Variants of the $ε$-LLBE are applicable to model pattern formation, phase transition, and heat conduction for non-simple materials, among other things. In this paper, we show well-posedness of the $ε$-LLBE and the convergence of the solution $\boldsymbol{u}^ε$ of the regularised equation to the solution $\boldsymbol{u}$ of the LLBE as $ε\to 0^+$. As a by-product of our analysis, we show the existence and uniqueness of regular solution to the LLBE for temperatures higher than the Curie temperature. Furthermore, we propose a linear fully discrete conforming finite element scheme to approximate the solution of the $ε$-LLBE. Error analysis is performed to show unconditional stability and optimal uniform-in-time convergence rate for the schemes. Several numerical simulations corroborate our theoretical results.
△ Less
Submitted 9 June, 2024;
originally announced June 2024.
-
Algorithmic aspects of left-orderings of solvable Baumslag--Solitar groups via its dynamical realization
Authors:
Meng-Che "Turbo" Ho,
Khanh Le,
Dino Rossegger
Abstract:
We answer a question of Calderoni and Clay by showing that the conjugation equivalence relation of left orderings of the Baumslag-Solitar groups $\mathrm{BS}(1,n)$ is hyperfinite for any $n$. Our proof relies on a classification of $\mathrm{BS}(1,n)$'s left-orderings via its one-dimensional dynamical realizations. We furthermore use the effectiveness of the dynamical realizations of…
▽ More
We answer a question of Calderoni and Clay by showing that the conjugation equivalence relation of left orderings of the Baumslag-Solitar groups $\mathrm{BS}(1,n)$ is hyperfinite for any $n$. Our proof relies on a classification of $\mathrm{BS}(1,n)$'s left-orderings via its one-dimensional dynamical realizations. We furthermore use the effectiveness of the dynamical realizations of $\mathrm{BS}(1,n)$ to study algorithmic properties of the left-orderings on $\mathrm{BS}(1,n)$.
△ Less
Submitted 14 May, 2024;
originally announced May 2024.
-
Quantitative Propagation of Chaos for Singular Interacting Particle Systems Driven by Fractional Brownian Motion
Authors:
Lucio Galeati,
Khoa Lê,
Avi Mayorcas
Abstract:
We consider interacting systems particle driven by i.i.d. fractional Brownian motions, subject to irregular, possibly distributional, pairwise interactions. We show propagation of chaos and mean field convergence to the law of the associated McKean--Vlasov equation, as the number of particles $N\to\infty$, with quantitative sharp rates of order $N^{-1/2}$. Our results hold for a wide class of poss…
▽ More
We consider interacting systems particle driven by i.i.d. fractional Brownian motions, subject to irregular, possibly distributional, pairwise interactions. We show propagation of chaos and mean field convergence to the law of the associated McKean--Vlasov equation, as the number of particles $N\to\infty$, with quantitative sharp rates of order $N^{-1/2}$. Our results hold for a wide class of possibly time-dependent interactions, which are only assumed to satisfy a Besov-type regularity, related to the Hurst parameter $H\in (0,+\infty)\setminus \mathbb{N}$ of the driving noises. In particular, as $H$ decreases to $0$, interaction kernels of arbitrary singularity can be considered, a phenomenon frequently observed in regularization by noise results. Our proofs rely on a combinations of Sznitman's direct comparison argument with stochastic sewing techniques.
△ Less
Submitted 8 March, 2024;
originally announced March 2024.
-
Zagier-Hoffman's conjectures in positive characteristic II
Authors:
Bo-Hae Im,
Hojin Kim,
Khac Nhuan Le,
Tuan Ngo Dac,
Lan Huong Pham
Abstract:
Zagier-Hoffman's conjectures predict the dimension and a basis for the $\mathbb Q$-vector spaces spanned by $N$th cyclotomic multiple zeta values (MZV's) of fixed weight where $N$ is a natural number.
For $N=1$ (MZV's case), half of these conjectures have been solved by the work of Terasoma, Deligne-Goncharov and Brown with the help of Zagier's identity. The other half are completely open. For…
▽ More
Zagier-Hoffman's conjectures predict the dimension and a basis for the $\mathbb Q$-vector spaces spanned by $N$th cyclotomic multiple zeta values (MZV's) of fixed weight where $N$ is a natural number.
For $N=1$ (MZV's case), half of these conjectures have been solved by the work of Terasoma, Deligne-Goncharov and Brown with the help of Zagier's identity. The other half are completely open. For $N=2$ (alternating MZV's case) and $N=3,4,8$, Deligne-Goncharov and Deligne solved the same half of these conjectures for $N$th-cyclotomic MZV's. For other values of $N$, no sharp upper bound on the dimension is known.
In this paper we completely establish, for all $N$, Zagier-Hoffman's conjectures for $N$th cyclotomic multiple zeta values in positive characteristic. By working with the tower of all cyclotomic extensions, we present a proof that is uniform on $N$ and give an effective algorithm to express any cyclotomic multiple zeta value in the chosen basis. This generalizes all previous work on these conjectures for MZV's and alternating MZV's in positive characteristic.
△ Less
Submitted 18 February, 2024;
originally announced February 2024.
-
Sliding Mode Observers for Set-valued Lur'e Systems with Uncertainties Beyond Observational Range
Authors:
Samir Adly,
Jun Huang,
Ba Khiet Le
Abstract:
In this paper, we introduce a new sliding mode observer for Lur'e set-valued dynamical systems, particularly addressing challenges posed by uncertainties not within the standard range of observation. Traditionally, most of Luenberger-like observers and sliding mode observer have been designed only for uncertainties in the range of observation. Central to our approach is the treatment of the uncert…
▽ More
In this paper, we introduce a new sliding mode observer for Lur'e set-valued dynamical systems, particularly addressing challenges posed by uncertainties not within the standard range of observation. Traditionally, most of Luenberger-like observers and sliding mode observer have been designed only for uncertainties in the range of observation. Central to our approach is the treatment of the uncertainty term which we decompose into two components: the first part in the observation subspace and the second part in its complemented subspace. We establish that when the second part converges to zero, an exact sliding mode observer for the system can be obtained. In scenarios where this convergence does not occur, our methodology allows for the estimation of errors between the actual state and the observer state. This leads to a practical interval estimation technique, valuable in situations where part of the uncertainty lies outside the observable range. Finally, we show that our observer is also a T- observer as well as a strong H-infinity observer.
△ Less
Submitted 25 April, 2024; v1 submitted 8 February, 2024;
originally announced February 2024.
-
R-Continuity with Applications to Convergence Analysis of Tikhonov Regularization and DC Programming
Authors:
Ba Khiet Le
Abstract:
In the paper, we study the convergence analysis of Tikhonov regularization in finding a zero of a maximal monotone operator using the notion of R-continuity. Applications to convex minimization and DC programming are provided.
In the paper, we study the convergence analysis of Tikhonov regularization in finding a zero of a maximal monotone operator using the notion of R-continuity. Applications to convex minimization and DC programming are provided.
△ Less
Submitted 21 November, 2023;
originally announced November 2023.
-
State-Dependent Sweeping Processes: Asymptotic Behavior and Algorithmic Approaches
Authors:
Samir Adly,
Monica G. Cojocaru,
Ba Khiet Le
Abstract:
In this paper, we investigate the asymptotic properties of a particular class of state-dependent sweeping processes. While extensive research has been conducted on the existence and uniqueness of solutions for sweeping processes, there is a scarcity of studies addressing their behavior in the limit of large time. Additionally, we introduce novel algorithms designed for the resolution of quasi-vari…
▽ More
In this paper, we investigate the asymptotic properties of a particular class of state-dependent sweeping processes. While extensive research has been conducted on the existence and uniqueness of solutions for sweeping processes, there is a scarcity of studies addressing their behavior in the limit of large time. Additionally, we introduce novel algorithms designed for the resolution of quasi-variational inequalities. As a result, we introduce a new derivative-free algorithm to find zeros of nonsmooth Lipschitz continuous mappings with a linear convergence rate. This algorithm can be effectively used in nonsmooth and nonconvex optimization problems that do not possess necessarily second-order differentiability conditions of the data.
△ Less
Submitted 21 November, 2023;
originally announced November 2023.
-
A new twist on modular links from an old perspective
Authors:
Khanh Le
Abstract:
We show that the complement of arithmetic modular links found in arXiv:2307.09409 is homeomorphic to the complement of augmented chainlinks. In particular, these link complements arise as n-fold cyclic covers of the Whitehead link complement.
We show that the complement of arithmetic modular links found in arXiv:2307.09409 is homeomorphic to the complement of augmented chainlinks. In particular, these link complements arise as n-fold cyclic covers of the Whitehead link complement.
△ Less
Submitted 20 November, 2023;
originally announced November 2023.
-
Asymptotically accurate and locking-free finite element implementation of first order shear deformation theory for plates
Authors:
Khanh Chau Le,
Hoang Giang Bui
Abstract:
A formulation of the asymptotically exact first-order shear deformation theory for linear-elastic homogeneous plates in the rescaled coordinates and rotation angles is considered. This allows the development of its asymptotically accurate and shear-locking-free finite element implementation. As applications, numerical simulations are performed for circular and rectangular plates, showing complete…
▽ More
A formulation of the asymptotically exact first-order shear deformation theory for linear-elastic homogeneous plates in the rescaled coordinates and rotation angles is considered. This allows the development of its asymptotically accurate and shear-locking-free finite element implementation. As applications, numerical simulations are performed for circular and rectangular plates, showing complete agreement between the analytical solution and the numerical solutions based on two-dimensional theory and three-dimensional elasticity theory.
△ Less
Submitted 16 April, 2024; v1 submitted 30 October, 2023;
originally announced October 2023.
-
A central limit theorem for the Euler method for SDEs with irregular drifts
Authors:
Konstantinos Dareiotis,
Máté Gerencsér,
Khoa Lê
Abstract:
The goal of this article is to establish a central limit theorem for the Euler-Maruyama scheme approximating multidimensional SDEs with elliptic Brownian diffusion, under very mild regularity requirements on the drift coefficients. When the drift is Hölder continuous, we show that the limiting law of the rescaled fluctuations around the true solution is characterised as the unique solution of a hy…
▽ More
The goal of this article is to establish a central limit theorem for the Euler-Maruyama scheme approximating multidimensional SDEs with elliptic Brownian diffusion, under very mild regularity requirements on the drift coefficients. When the drift is Hölder continuous, we show that the limiting law of the rescaled fluctuations around the true solution is characterised as the unique solution of a hybrid Young-Itô differential equation. When the drift has positive Sobolev regularity, this limit is characterised by the solution of a transformed SDE. Our result is an extension of the results of Jacod-Kurtz-Protter (1991, 1998) in which SDEs with differentiable coefficients were considered. To compensate for the lack of regularity of the drifts, we utilize the regularisation effect from the non-degenerate noise.
△ Less
Submitted 28 September, 2023;
originally announced September 2023.
-
A class of space-time discretizations for the stochastic $p$-Stokes system
Authors:
Kim-Ngan Le,
Jörn Wichmann
Abstract:
The main objective of the present paper is to construct a new class of space-time discretizations for the stochastic $p$-Stokes system and analyze its stability and convergence properties.
We derive regularity results for the approximation that are similar to the natural regularity of solutions. One of the key arguments relies on discrete extrapolation that allows to relate lower moments of disc…
▽ More
The main objective of the present paper is to construct a new class of space-time discretizations for the stochastic $p$-Stokes system and analyze its stability and convergence properties.
We derive regularity results for the approximation that are similar to the natural regularity of solutions. One of the key arguments relies on discrete extrapolation that allows to relate lower moments of discrete maximal processes.
We show that, if the generic spatial discretization is constraint conforming, then the velocity approximation satisfies a best-approximation property in the natural distance.
Moreover, we present an example such that the resulting velocity approximation converges with rate $1/2$ in time and $1$ in space towards the (unknown) target velocity with respect to the natural distance.
△ Less
Submitted 25 July, 2023;
originally announced July 2023.
-
On a Conjecture of Gezmis and Pellarin
Authors:
Khac Nhuan Le,
Kien Huu Nguyen
Abstract:
In 2022, Gezmis and Pellarin introduced and studied the concept of trivial multiple zeta values, along with a map from the vector space spanned by these values to the vector space spanned by Thakur's multiple zeta values. Their construction allows us to generate some linear relations among the latter values using the former. In our work, we determine the structure of the kernel of the aforemention…
▽ More
In 2022, Gezmis and Pellarin introduced and studied the concept of trivial multiple zeta values, along with a map from the vector space spanned by these values to the vector space spanned by Thakur's multiple zeta values. Their construction allows us to generate some linear relations among the latter values using the former. In our work, we determine the structure of the kernel of the aforementioned map. As a consequence, we give an answer to a conjecture proposed by Gezmis and Pellarin regarding the injectivity of this specific map.
△ Less
Submitted 22 June, 2023;
originally announced June 2023.
-
Numerical analysis of the stochastic Stefan problem
Authors:
Jerome Droniou,
Muhammad Awais Khan,
Kim Ngan Le
Abstract:
The gradient discretisation method (GDM) -- a generic framework encompassing many numerical methods -- is studied for a general stochastic Stefan problem with multiplicative noise. The convergence of the numerical solutions is proved by compactness method using discrete functional analysis tools, Skorohod theorem and the martingale representation theorem. The generic convergence results establishe…
▽ More
The gradient discretisation method (GDM) -- a generic framework encompassing many numerical methods -- is studied for a general stochastic Stefan problem with multiplicative noise. The convergence of the numerical solutions is proved by compactness method using discrete functional analysis tools, Skorohod theorem and the martingale representation theorem. The generic convergence results established in the GDM framework are applicable to a range of different numerical methods, including for example mass-lumped finite elements, but also some finite volume methods, mimetic methods, lowest-order virtual element methods, etc. Theoretical results are complemented by numerical tests based on two methods that fit in GDM framework.
△ Less
Submitted 26 June, 2023; v1 submitted 22 June, 2023;
originally announced June 2023.
-
Path-by-path uniqueness for stochastic differential equations under Krylov-Röckner condition
Authors:
Lukas Anzeletti,
Khoa Lê,
Chengcheng Ling
Abstract:
We show that any stochastic differential equation (SDE) driven by Brownian motion with drift satisfying the Krylov-Röckner condition has exactly one solution in an ordinary sense for almost every trajectory of the Brownian motion. Additionally, we show that such SDE is strongly complete, i.e. for almost every trajectory of the Brownian motion, the family of solutions with different initial data fo…
▽ More
We show that any stochastic differential equation (SDE) driven by Brownian motion with drift satisfying the Krylov-Röckner condition has exactly one solution in an ordinary sense for almost every trajectory of the Brownian motion. Additionally, we show that such SDE is strongly complete, i.e. for almost every trajectory of the Brownian motion, the family of solutions with different initial data forms a continuous semiflow for all nonnegative times.
△ Less
Submitted 13 April, 2023;
originally announced April 2023.
-
Hopf algebras and alternating multiple zeta values in positive characteristic
Authors:
Bo-Hae Im,
Hojin Kim,
Khac Nhuan Le,
Tuan Ngo Dac,
Lan Huong Pham
Abstract:
In \cite{IKLNDP23} we presented a systematic study of algebra structures of multiple zeta values in positive characteristic introduced by Thakur as analogues of classical multiple zeta values of Euler. In this paper we construct algebra and Hopf algebra structures of alternating multiple zeta values introduced by Harada, extending our previous work. Our results could be considered as an analogue o…
▽ More
In \cite{IKLNDP23} we presented a systematic study of algebra structures of multiple zeta values in positive characteristic introduced by Thakur as analogues of classical multiple zeta values of Euler. In this paper we construct algebra and Hopf algebra structures of alternating multiple zeta values introduced by Harada, extending our previous work. Our results could be considered as an analogue of those of Hoffman \cite{Hof00} and Racinet \cite{Rac02} in the classical setting. The proof is based on two new ingredients: the first one is a direct and explicit construction of the shuffle Hopf algebra structure, and the second one is the notion of horizontal maps.
△ Less
Submitted 5 April, 2023;
originally announced April 2023.
-
Sliding Mode Observer for Set-valued Lur'e Systems and Chattering Removing
Authors:
Samir Adly,
Ba Khiet Le
Abstract:
In this paper, we study a sliding mode observer for a class of set-valued Lur'e systems subject to uncertainties. We show that our approach has obvious advantages than the existing Luenberger-like observers. Furthermore, we provide an effective continuous approximation to eliminate the chattering effect in the sliding mode technique.
In this paper, we study a sliding mode observer for a class of set-valued Lur'e systems subject to uncertainties. We show that our approach has obvious advantages than the existing Luenberger-like observers. Furthermore, we provide an effective continuous approximation to eliminate the chattering effect in the sliding mode technique.
△ Less
Submitted 11 June, 2023; v1 submitted 30 March, 2023;
originally announced March 2023.
-
Stochastic equations with singular drift driven by fractional Brownian motion
Authors:
Oleg Butkovsky,
Khoa Lê,
Leonid Mytnik
Abstract:
We consider stochastic differential equation $$ d X_t=b(X_t) dt +d W_t^H, $$ where the drift $b$ is either a measure or an integrable function, and $W^H$ is a $d$-dimensional fractional Brownian motion with Hurst parameter $H\in(0,1)$, $d\in\mathbb{N}$. For the case where $b\in L_p(\mathbb{R}^d)$, $p\in[1,\infty]$ we show weak existence of solutions to this equation under the condition…
▽ More
We consider stochastic differential equation $$ d X_t=b(X_t) dt +d W_t^H, $$ where the drift $b$ is either a measure or an integrable function, and $W^H$ is a $d$-dimensional fractional Brownian motion with Hurst parameter $H\in(0,1)$, $d\in\mathbb{N}$. For the case where $b\in L_p(\mathbb{R}^d)$, $p\in[1,\infty]$ we show weak existence of solutions to this equation under the condition $$ \frac{d}p<\frac1H-1, $$ which is an extension of the Krylov-Röckner condition (2005) to the fractional case. We construct a counter-example showing optimality of this condition. If $b$ is a Radon measure, particularly the delta measure, we prove weak existence of solutions to this equation under the optimal condition $H<\frac1{d+1}$. We also show strong well-posedness of solutions to this equation under certain conditions. To establish these results, we utilize the stochastic sewing technique and develop a new version of the stochastic sewing lemma.
△ Less
Submitted 23 February, 2023;
originally announced February 2023.
-
Hopf algebras and multiple zeta values in positive characteristic
Authors:
Bo-Hae Im,
Hojin Kim,
Khac Nhuan Le,
Tuan Ngo Dac,
Lan Huong Pham
Abstract:
Multiples zeta values (MZV's for short) in positive characteristic were introduced by Thakur as analogues of classical multiple zeta values of Euler. In this paper we give a systematic study of algebraic structures of MZV's in positive characteristic. We construct both the stuffle algebra and the shuffle algebra of these MZV's and equip them with algebra and Hopf algebra structures. In particular,…
▽ More
Multiples zeta values (MZV's for short) in positive characteristic were introduced by Thakur as analogues of classical multiple zeta values of Euler. In this paper we give a systematic study of algebraic structures of MZV's in positive characteristic. We construct both the stuffle algebra and the shuffle algebra of these MZV's and equip them with algebra and Hopf algebra structures. In particular, we completely solve a problem suggested by Deligne and Thakur \cite{Del17} in 2017 and establish Shi's conjectures \cite{Shi18}. The construction of the stuffle algebra is based on our recent work \cite{IKLNDP22}.
△ Less
Submitted 14 January, 2023;
originally announced January 2023.
-
Numerical method and Error estimate for stochastic Landau--Lifshitz--Bloch equation
Authors:
Beniamin Goldys,
Chunxi Jiao,
Kim-Ngan Le
Abstract:
We study numerical methods for solving a system of quasilinear stochastic partial differential equations known as the stochastic Landau-Lifshitz-Bloch (LLB) equation on a bounded domain in $\mathbb R^d$ for $d=1,2$. Our main results are estimates of the rate of convergence of the Finite Element Method to the solutions of stochastic LLB. To overcome the lack of regularity of the solution in the cas…
▽ More
We study numerical methods for solving a system of quasilinear stochastic partial differential equations known as the stochastic Landau-Lifshitz-Bloch (LLB) equation on a bounded domain in $\mathbb R^d$ for $d=1,2$. Our main results are estimates of the rate of convergence of the Finite Element Method to the solutions of stochastic LLB. To overcome the lack of regularity of the solution in the case $d=2$, we propose a Finite Element scheme for a regularised version of the equation. We then obtain error estimates of numerical solutions and for the solution of the regularised equation as well as the rate of convergence of this solution to the solution of the stochastic LLB equation. As a consequence, the convergence in probability of the approximate solutions to the solution of the stochastic LLB equation is derived. To the best of our knowledge this is the first result on error estimates for a system of stochastic quasilinear partial differential equations. A stronger result is obtained in the case $d=1$ due to a new regularity result for the LLB equation which allows us to avoid regularisation.
△ Less
Submitted 21 December, 2022;
originally announced December 2022.
-
Maximal inequalities and weighted BMO processes
Authors:
Khoa Lê
Abstract:
For a general adapted integrable right-continuous with left limits (RCLL) process $(X_t)_{t\in[0,τ]}$ taking values in a metric space $(\mathcal E,d)$, we show (among other things) that for every $m\in(1,\infty)$…
▽ More
For a general adapted integrable right-continuous with left limits (RCLL) process $(X_t)_{t\in[0,τ]}$ taking values in a metric space $(\mathcal E,d)$, we show (among other things) that for every $m\in(1,\infty)$
$$
\frac{m-1}{2m-1}\|\sup_{t\in[0,τ]}\mathbb{E}(d(X_{t-},X_τ)|\mathcal F_t)\|_m\le \|\sup_{t\in[0,τ]}d(X_0,X_t)\|_m\le c\frac{m^2}{m-1} \|\sup_{t\in[0,τ]}\mathbb{E}(d(X_{t-},X_τ)|\mathcal F_t)\|_m
$$
with a universal constant $c$. This is a probabilistic version of Fefferman--Stein estimate for the sharp maximal functions. While the former inequality is derived easily from Doob's martingale inequality, the later inequality is a consequence of John--Nirenberg inequalities for weighted BMO processes, which are obtained in this note. We explain how John--Nirenberg inequalities can be utilized to obtain inequalities for martingales, both old and new alike in a unified way.
△ Less
Submitted 19 December, 2022; v1 submitted 28 November, 2022;
originally announced November 2022.
-
Quantitative John--Nirenberg inequality for stochastic processes of bounded mean oscillation
Authors:
Khoa Lê
Abstract:
Stroock and Varadhan in 1997 and Geiss in 2005 independently introduced stochastic processes with bounded mean oscillation (BMO) and established their exponential integrability with some unspecified exponential constant. This result is an analogue of the John--Nirenberg inequality for functions of bounded mean oscillation. In this work, we quantify the size of the exponential constant by the modul…
▽ More
Stroock and Varadhan in 1997 and Geiss in 2005 independently introduced stochastic processes with bounded mean oscillation (BMO) and established their exponential integrability with some unspecified exponential constant. This result is an analogue of the John--Nirenberg inequality for functions of bounded mean oscillation. In this work, we quantify the size of the exponential constant by the modulus of mean oscillation. Some new applications of BMO processes in rough stochastic differential equations, numerical approximations and regularization by noise are discussed.
△ Less
Submitted 11 November, 2022; v1 submitted 27 October, 2022;
originally announced October 2022.
-
Linear recurrences of order at most two in nontrivial small divisors and large divisors
Authors:
Hung Viet Chu,
Kevin Huu Le,
Steven J. Miller,
Yuan Qiu,
Liyang Shen
Abstract:
For each positive integer $N$, define $$S'_N \ =\ \{1 < d < \sqrt{N}: d|N\}\mbox{ and }L'_N \ =\ \{\sqrt{N} < d < N : d|N\}.$$ Recently, Chentouf characterized all positive integers $N$ such that the set of small divisors $\{d\le \sqrt{N}: d|N\}$ satisfies a linear recurrence of order at most two. We nontrivially extend the result by excluding the trivial divisor $1$ from consideration, which dram…
▽ More
For each positive integer $N$, define $$S'_N \ =\ \{1 < d < \sqrt{N}: d|N\}\mbox{ and }L'_N \ =\ \{\sqrt{N} < d < N : d|N\}.$$ Recently, Chentouf characterized all positive integers $N$ such that the set of small divisors $\{d\le \sqrt{N}: d|N\}$ satisfies a linear recurrence of order at most two. We nontrivially extend the result by excluding the trivial divisor $1$ from consideration, which dramatically increases the analysis complexity. Our first result characterizes all positive integers $N$ such that $S'_N$ satisfies a linear recurrence of order at most two. Moreover, our second result characterizes all positive $N$ such that $L'_N$ satisfies a linear recurrence of order at most two, thus extending considerably a recent result that characterizes $N$ with $L'_N$ being in an arithmetic progression.
△ Less
Submitted 1 October, 2022;
originally announced October 2022.
-
Existence, uniqueness and regularity of solutions to the stochastic Landau-Lifschitz-Slonczewski equation
Authors:
Beniamin Goldys,
Chunxi Jiao,
Kim Ngan Le
Abstract:
In this paper we are concerned with the stochastic Landau-Lifshitz-Slonczewski equation (LLS) that describes magnetisation of an inifnite nanowire evolving under current driven spin torque. The current brings into the system a multiplicative gradient noise that appears as a transport term in the equation. We prove the existence, uniqueness and regularity of pathwise solutions to this equation.
In this paper we are concerned with the stochastic Landau-Lifshitz-Slonczewski equation (LLS) that describes magnetisation of an inifnite nanowire evolving under current driven spin torque. The current brings into the system a multiplicative gradient noise that appears as a transport term in the equation. We prove the existence, uniqueness and regularity of pathwise solutions to this equation.
△ Less
Submitted 22 August, 2022;
originally announced August 2022.
-
Geodesic surfaces in the complement of knots with small crossing number
Authors:
Khanh Le,
Rebekah Palmer
Abstract:
In this article, we investigate the problem of counting totally geodesic surfaces in the complement of hyperbolic knots with at most 9 crossings. Adapting previous counting techniques of boundary slope and intersection, we establish uniqueness of a totally geodesic surface for the knots $7_4$ and $9_{35}$. Extending an obstruction to the existence of totally geodesic surfaces due to Calegari, we s…
▽ More
In this article, we investigate the problem of counting totally geodesic surfaces in the complement of hyperbolic knots with at most 9 crossings. Adapting previous counting techniques of boundary slope and intersection, we establish uniqueness of a totally geodesic surface for the knots $7_4$ and $9_{35}$. Extending an obstruction to the existence of totally geodesic surfaces due to Calegari, we show that there is no totally geodesic surface in the complement of 47 knots.
△ Less
Submitted 16 March, 2023; v1 submitted 29 July, 2022;
originally announced July 2022.
-
Zagier-Hoffman's conjectures in positive characteristic
Authors:
Bo-Hae Im,
Hojin Kim,
Khac Nhuan Le,
Tuan Ngo Dac,
Lan Huong Pham
Abstract:
Multiples zeta values and alternating multiple zeta values in positive characteristic were introduced by Thakur and Harada as analogues of classical multiple zeta values of Euler and Euler sums. In this paper we determine all linear relations among alternating multiple zeta values and settle the main goals of these theories. As a consequence we completely establish Zagier-Hoffman's conjectures in…
▽ More
Multiples zeta values and alternating multiple zeta values in positive characteristic were introduced by Thakur and Harada as analogues of classical multiple zeta values of Euler and Euler sums. In this paper we determine all linear relations among alternating multiple zeta values and settle the main goals of these theories. As a consequence we completely establish Zagier-Hoffman's conjectures in positive characteristic formulated by Todd and Thakur which predict the dimension and an explicit basis of the span of multiple zeta values of Thakur of fixed weight.
△ Less
Submitted 10 June, 2024; v1 submitted 14 May, 2022;
originally announced May 2022.
-
The Allen-Cahn equation with generic initial datum
Authors:
Martin Hairer,
Khoa Lê,
Tommaso Rosati
Abstract:
We consider the Allen-Cahn equation $\partial_t u- Δu=u-u^3$ with a rapidly mixing Gaussian field as initial condition. We show that provided that the amplitude of the initial condition is not too large, the equation generates fronts described by nodal sets of the Bargmann-Fock Gaussian field, which then evolve according to mean curvature flow.
We consider the Allen-Cahn equation $\partial_t u- Δu=u-u^3$ with a rapidly mixing Gaussian field as initial condition. We show that provided that the amplitude of the initial condition is not too large, the equation generates fronts described by nodal sets of the Bargmann-Fock Gaussian field, which then evolve according to mean curvature flow.
△ Less
Submitted 20 January, 2022;
originally announced January 2022.
-
Multi-valued variational inequalities for variable exponent double phase problems: comparison and extremality results
Authors:
Siegfried Carl,
Vy Khoi Le,
Patrick Winkert
Abstract:
We prove existence and comparison results for multi-valued variational inequalities in a bounded domain $Ω$ of the form \begin{equation*} u\in K\,:\, 0 \in Au+\partial I_K(u)+\mathcal{F}(u)+\mathcal{F}_Γ(u)\quad\text{in }W^{1,\mathcal{H}}(Ω)^*, \end{equation*} where $A\colon W^{1, \mathcal{H}}(Ω) \to W^{1, \mathcal{H}}(Ω)^*$ given by \begin{equation*} Au:=-\text{div}\left(|\nabla u|^{p(x)-2} \nabl…
▽ More
We prove existence and comparison results for multi-valued variational inequalities in a bounded domain $Ω$ of the form \begin{equation*} u\in K\,:\, 0 \in Au+\partial I_K(u)+\mathcal{F}(u)+\mathcal{F}_Γ(u)\quad\text{in }W^{1,\mathcal{H}}(Ω)^*, \end{equation*} where $A\colon W^{1, \mathcal{H}}(Ω) \to W^{1, \mathcal{H}}(Ω)^*$ given by \begin{equation*} Au:=-\text{div}\left(|\nabla u|^{p(x)-2} \nabla u+ μ(x) |\nabla u|^{q(x)-2} \nabla u\right) \end{equation*} for $u \in W^{1, \mathcal{H}}(Ω)$, is the double phase operator with variable exponents and $W^{1, \mathcal{H}}(Ω)$ is the associated Musielak-Orlicz Sobolev space. First, an existence result is proved under some weak coercivity condition. Our main focus aims at the treatment of the problem under consideration when coercivity fails. To this end we establish the method of sub-supersolution for the multi-valued variational inequality in the space $W^{1, \mathcal{H}}(Ω)$ based on appropriately defined sub- and supersolutions, which yields the existence of solutions within an ordered interval of sub-supersolution. Moreover, the existence of extremal solutions will be shown provided the closed, convex subset $K$ of $W^{1, \mathcal{H}}(Ω)$ satisfies a lattice condition. As an application of the sub-supersolution method we are able to show that a class of generalized variational-hemivariational inequalities with a leading double phase operator are included as a special case of the multi-valued variational inequality considered here. Based on a fixed point argument, we also study the case when the corresponding Nemytskij operators $\mathcal{F}, \mathcal{F}_Γ$ need not be continuous. At the end, we give a nontrivial example of the construction of sub- and supersolutions related to the problem above.
△ Less
Submitted 30 March, 2023; v1 submitted 8 January, 2022;
originally announced January 2022.
-
Taming singular stochastic differential equations: A numerical method
Authors:
Khoa Lê,
Chengcheng Ling
Abstract:
We consider a generic and explicit tamed Euler--Maruyama scheme for multidimensional time-inhomogeneous stochastic differential equations with multiplicative Brownian noise. The diffusion coefficient is uniformly elliptic, Hölder continuous and weakly differentiable in the spatial variables while the drift satisfies the Ladyzhenskaya--Prodi--Serrin condition, as considered by Krylov and Röckner (2…
▽ More
We consider a generic and explicit tamed Euler--Maruyama scheme for multidimensional time-inhomogeneous stochastic differential equations with multiplicative Brownian noise. The diffusion coefficient is uniformly elliptic, Hölder continuous and weakly differentiable in the spatial variables while the drift satisfies the Ladyzhenskaya--Prodi--Serrin condition, as considered by Krylov and Röckner (2005). In the discrete scheme, the drift is tamed by replacing it by an approximation. A strong rate of convergence of the scheme is provided in terms of the approximation error of the drift in a suitable and possibly very weak topology. A few examples of approximating drifts are discussed in detail. The parameters of the approximating drifts can vary and be fine-tuned to achieve the standard 1/2-strong convergence rate with a logarithmic factor. The result is then applied to provide numerical solutions for stochastic transport equations with singular vector fields satisfying the aforementioned condition.
△ Less
Submitted 14 June, 2022; v1 submitted 4 October, 2021;
originally announced October 2021.
-
Left orderability for surgeries on the $[1,1,2,2,2j]$ two-bridge knots
Authors:
Khanh Le
Abstract:
Let $M$ be a $\mathbb{Q}$-homology solid torus. In this paper, we give a cohomological criterion for the existence of an interval of left-orderable Dehn surgeries on $M$. We apply this criterion to prove that the two-bridge knot that corresponds to the continued fraction $[1,1,2,2,2j]$ for $j\geq 1$ admits an interval of left-orderable Dehn surgeries. This family of two-bridge knots gives some pos…
▽ More
Let $M$ be a $\mathbb{Q}$-homology solid torus. In this paper, we give a cohomological criterion for the existence of an interval of left-orderable Dehn surgeries on $M$. We apply this criterion to prove that the two-bridge knot that corresponds to the continued fraction $[1,1,2,2,2j]$ for $j\geq 1$ admits an interval of left-orderable Dehn surgeries. This family of two-bridge knots gives some positive evidence for a question of Xinghua Gao.
△ Less
Submitted 13 September, 2021;
originally announced September 2021.
-
Entropic Gromov-Wasserstein between Gaussian Distributions
Authors:
Khang Le,
Dung Le,
Huy Nguyen,
Dat Do,
Tung Pham,
Nhat Ho
Abstract:
We study the entropic Gromov-Wasserstein and its unbalanced version between (unbalanced) Gaussian distributions with different dimensions. When the metric is the inner product, which we refer to as inner product Gromov-Wasserstein (IGW), we demonstrate that the optimal transportation plans of entropic IGW and its unbalanced variant are (unbalanced) Gaussian distributions. Via an application of von…
▽ More
We study the entropic Gromov-Wasserstein and its unbalanced version between (unbalanced) Gaussian distributions with different dimensions. When the metric is the inner product, which we refer to as inner product Gromov-Wasserstein (IGW), we demonstrate that the optimal transportation plans of entropic IGW and its unbalanced variant are (unbalanced) Gaussian distributions. Via an application of von Neumann's trace inequality, we obtain closed-form expressions for the entropic IGW between these Gaussian distributions. Finally, we consider an entropic inner product Gromov-Wasserstein barycenter of multiple Gaussian distributions. We prove that the barycenter is a Gaussian distribution when the entropic regularization parameter is small. We further derive a closed-form expression for the covariance matrix of the barycenter.
△ Less
Submitted 24 February, 2022; v1 submitted 24 August, 2021;
originally announced August 2021.
-
On Multimarginal Partial Optimal Transport: Equivalent Forms and Computational Complexity
Authors:
Khang Le,
Huy Nguyen,
Tung Pham,
Nhat Ho
Abstract:
We study the multi-marginal partial optimal transport (POT) problem between $m$ discrete (unbalanced) measures with at most $n$ supports. We first prove that we can obtain two equivalence forms of the multimarginal POT problem in terms of the multimarginal optimal transport problem via novel extensions of cost tensor. The first equivalence form is derived under the assumptions that the total masse…
▽ More
We study the multi-marginal partial optimal transport (POT) problem between $m$ discrete (unbalanced) measures with at most $n$ supports. We first prove that we can obtain two equivalence forms of the multimarginal POT problem in terms of the multimarginal optimal transport problem via novel extensions of cost tensor. The first equivalence form is derived under the assumptions that the total masses of each measure are sufficiently close while the second equivalence form does not require any conditions on these masses but at the price of more sophisticated extended cost tensor. Our proof techniques for obtaining these equivalence forms rely on novel procedures of moving mass in graph theory to push transportation plan into appropriate regions. Finally, based on the equivalence forms, we develop optimization algorithm, named ApproxMPOT algorithm, that builds upon the Sinkhorn algorithm for solving the entropic regularized multimarginal optimal transport. We demonstrate that the ApproxMPOT algorithm can approximate the optimal value of multimarginal POT problem with a computational complexity upper bound of the order $\tilde{\mathcal{O}}(m^3(n+1)^{m}/ \varepsilon^2)$ where $\varepsilon > 0$ stands for the desired tolerance.
△ Less
Submitted 24 February, 2022; v1 submitted 18 August, 2021;
originally announced August 2021.
-
Rough stochastic differential equations
Authors:
Peter K. Friz,
Antoine Hocquet,
Khoa Lê
Abstract:
We build a hybrid theory of rough stochastic analysis which seamlessly combines the advantages of both Itô's stochastic and Lyons' rough differential equations. This gives a direct and intrinsic understanding of multidimensional diffusion with Brownian noise $(B,\tilde B)$ $$
dY_t(ω)=b(ω,t,Y_t(ω))dt+σ(ω,t,Y_t(ω))dB_t(ω)+f({ω,t},Y_t(ω))d\tilde B_t \,,
$$ in the annealed form, when conditioned o…
▽ More
We build a hybrid theory of rough stochastic analysis which seamlessly combines the advantages of both Itô's stochastic and Lyons' rough differential equations. This gives a direct and intrinsic understanding of multidimensional diffusion with Brownian noise $(B,\tilde B)$ $$
dY_t(ω)=b(ω,t,Y_t(ω))dt+σ(ω,t,Y_t(ω))dB_t(ω)+f({ω,t},Y_t(ω))d\tilde B_t \,,
$$ in the annealed form, when conditioned on its environmental noise $\tilde B$. This situation arises naturally e.g. in filtering theory, for Feynman--Kac representations of solutions to stochastic partial differential equations, in Lions--Souganidis' theory of pathwise stochastic control, and for McKean--Vlasov stochastic differential equations with common noise. In fact, we establish well-posedness of rough stochastic differential equations, with $\tilde B$ replaced by a genuine rough path. As consequence, the `annealed' process $Y$ is a locally Lipschitz function of its environmental noise in rough path metrics. There is also interest in taking $\tilde B=\tilde B^H$, a fractional Brownian motion which fits our theory for $H>1/3$. Our assumptions for $b,σ$ agree with those from Itô theory, those for $f$ with rough paths theory, including an extension of Davie's critical regularity result for deterministic rough differential equations. A major role in our analysis is played by a new scale of stochastic controlled rough paths spaces, related to a $(L^m,L^n)$-variant of stochastic sewing.
△ Less
Submitted 28 March, 2023; v1 submitted 18 June, 2021;
originally announced June 2021.
-
Stochastic sewing in Banach space
Authors:
Khoa Lê
Abstract:
A stochastic sewing lemma which is applicable for processes taking values in Banach spaces is introduced. Applications to additive functionals of fractional Brownian motion of distributional type are discussed.
A stochastic sewing lemma which is applicable for processes taking values in Banach spaces is introduced. Applications to additive functionals of fractional Brownian motion of distributional type are discussed.
△ Less
Submitted 9 January, 2022; v1 submitted 19 May, 2021;
originally announced May 2021.
-
On Robust Optimal Transport: Computational Complexity and Barycenter Computation
Authors:
Khang Le,
Huy Nguyen,
Quang Nguyen,
Tung Pham,
Hung Bui,
Nhat Ho
Abstract:
We consider robust variants of the standard optimal transport, named robust optimal transport, where marginal constraints are relaxed via Kullback-Leibler divergence. We show that Sinkhorn-based algorithms can approximate the optimal cost of robust optimal transport in $\widetilde{\mathcal{O}}(\frac{n^2}{\varepsilon})$ time, in which $n$ is the number of supports of the probability distributions a…
▽ More
We consider robust variants of the standard optimal transport, named robust optimal transport, where marginal constraints are relaxed via Kullback-Leibler divergence. We show that Sinkhorn-based algorithms can approximate the optimal cost of robust optimal transport in $\widetilde{\mathcal{O}}(\frac{n^2}{\varepsilon})$ time, in which $n$ is the number of supports of the probability distributions and $\varepsilon$ is the desired error. Furthermore, we investigate a fixed-support robust barycenter problem between $m$ discrete probability distributions with at most $n$ number of supports and develop an approximating algorithm based on iterative Bregman projections (IBP). For the specific case $m = 2$, we show that this algorithm can approximate the optimal barycenter value in $\widetilde{\mathcal{O}}(\frac{mn^2}{\varepsilon})$ time, thus being better than the previous complexity $\widetilde{\mathcal{O}}(\frac{mn^2}{\varepsilon^2})$ of the IBP algorithm for approximating the Wasserstein barycenter.
△ Less
Submitted 27 October, 2021; v1 submitted 12 February, 2021;
originally announced February 2021.
-
Quantifying a convergence theorem of Gyöngy and Krylov
Authors:
Konstantinos Dareiotis,
Máté Gerencsér,
Khoa Lê
Abstract:
We derive sharp strong convergence rates for the Euler-Maruyama scheme approximating multidimensional SDEs with multiplicative noise without imposing any regularity condition on the drift coefficient. In case the noise is additive, we show that Sobolev regularity can be leveraged to obtain improved rate: drifts with regularity of order $α\in (0,1)$ lead to rate $(1+α)/2$.
We derive sharp strong convergence rates for the Euler-Maruyama scheme approximating multidimensional SDEs with multiplicative noise without imposing any regularity condition on the drift coefficient. In case the noise is additive, we show that Sobolev regularity can be leveraged to obtain improved rate: drifts with regularity of order $α\in (0,1)$ lead to rate $(1+α)/2$.
△ Less
Submitted 21 August, 2022; v1 submitted 28 January, 2021;
originally announced January 2021.
-
Well-posedness of stochastic heat equation with distributional drift and skew stochastic heat equation
Authors:
Siva Athreya,
Oleg Butkovsky,
Khoa Lê,
Leonid Mytnik
Abstract:
We study stochastic reaction--diffusion equation $$ \partial_tu_t(x)=\frac12 \partial^2_{xx}u_t(x)+b(u_t(x))+\dot{W}_{t}(x), \quad t>0,\, x\in D $$ where $b$ is a generalized function in the Besov space $\mathcal{B}^β_{q,\infty}({\mathbb R})$, $D\subset{\mathbb R}$ and $\dot W$ is a space-time white noise on ${\mathbb R}_+\times D$. We introduce a notion of a solution to this equation and obtain e…
▽ More
We study stochastic reaction--diffusion equation $$ \partial_tu_t(x)=\frac12 \partial^2_{xx}u_t(x)+b(u_t(x))+\dot{W}_{t}(x), \quad t>0,\, x\in D $$ where $b$ is a generalized function in the Besov space $\mathcal{B}^β_{q,\infty}({\mathbb R})$, $D\subset{\mathbb R}$ and $\dot W$ is a space-time white noise on ${\mathbb R}_+\times D$. We introduce a notion of a solution to this equation and obtain existence and uniqueness of a strong solution whenever $β-1/q\ge-1$, $β>-1$ and $q\in[1,\infty]$. This class includes equations with $b$ being measures, in particular, $b=δ_0$ which corresponds to the skewed stochastic heat equation. For $β-1/q > -3/2$, we obtain existence of a weak solution. Our results extend the work of Bass and Chen (2001) to the framework of stochastic partial differential equations and generalizes the results of Gyöngy and Pardoux (1993) to distributional drifts. To establish these results, we exploit the regularization effect of the white noise through a new strategy based on the stochastic sewing lemma introduced in Lê~(2020).
△ Less
Submitted 11 February, 2022; v1 submitted 26 November, 2020;
originally announced November 2020.
-
Totally geodesic surfaces in twist knot complements
Authors:
Khanh Le,
Rebekah Palmer
Abstract:
In this article, we give explicit examples of infinitely many non-commensurable (non-arithmetic) hyperbolic $3$-manifolds admitting exactly $k$ totally geodesic surfaces for any positive integer $k$, answering a question of Bader, Fisher, Miller and Stover. The construction comes from a family of twist knot complements and their dihedral covers. The case $k=1$ arises from the uniqueness of an imme…
▽ More
In this article, we give explicit examples of infinitely many non-commensurable (non-arithmetic) hyperbolic $3$-manifolds admitting exactly $k$ totally geodesic surfaces for any positive integer $k$, answering a question of Bader, Fisher, Miller and Stover. The construction comes from a family of twist knot complements and their dihedral covers. The case $k=1$ arises from the uniqueness of an immersed totally geodesic thrice-punctured sphere, answering a question of Reid. Applying the proof techniques of the main result, we explicitly construct non-elementary maximal Fuchsian subgroups of infinite covolume within twist knot groups, and we also show that no twist knot complement with odd prime half twists is right-angled in the sense of Champanerkar, Kofman, and Purcell.
△ Less
Submitted 28 January, 2022; v1 submitted 9 September, 2020;
originally announced September 2020.
-
What are GT-shadows?
Authors:
Vasily A. Dolgushev,
Khanh Q. Le,
Aidan A. Lorenz
Abstract:
Let $B_4$ (resp. $PB_4$) be the braid group (resp. the pure braid group) on 4 strands and $NFI_{PB_4}(B_4)$ be the poset whose objects are finite index normal subgroups of $B_4$ that are contained in $PB_4$. In this paper, we introduce GT-shadows which may be thought of as "approximations" to elements of the profinite version $\widehat{GT}$ of the Grothendieck-Teichmueller group (see V. Drinfeld,…
▽ More
Let $B_4$ (resp. $PB_4$) be the braid group (resp. the pure braid group) on 4 strands and $NFI_{PB_4}(B_4)$ be the poset whose objects are finite index normal subgroups of $B_4$ that are contained in $PB_4$. In this paper, we introduce GT-shadows which may be thought of as "approximations" to elements of the profinite version $\widehat{GT}$ of the Grothendieck-Teichmueller group (see V. Drinfeld, Algebra i Analiz, 1990). We prove that GT-shadows form a groupoid whose objects are elements of $NFI_{PB_4}(B_4)$. We show that GT-shadows coming from elements of $\widehat{GT}$ satisfy various additional properties and we investigate these properties. We establish an explicit link between GT-shadows and the group $\widehat{GT}$ (see Theorem 3.8). We also present selected results of computer experiments on GT-shadows. In the appendix of this paper, we give a complete description of GT-shadows in the Abelian setting. We also prove that, in the Abelian setting, every GT-shadow comes from an element of $\widehat{GT}$. Objects very similar to GT-shadows were introduced in a paper by D. Harbater and L. Schneps in 1997. A variation of the concept of $GT$-shadows for the coarse version of $\widehat{GT}$ was studied in papers by P. Guillot.
△ Less
Submitted 31 October, 2021; v1 submitted 31 July, 2020;
originally announced August 2020.
-
Design and convergence analysis of numerical methods for stochastic evolution equations with Leray-Lions operator
Authors:
Jerome Droniou,
Beniamin Goldys,
Kim-Ngan Le
Abstract:
*The gradient discretisation method (GDM) is a generic framework, covering many classical methods (Finite Elements, Finite Volumes, Discontinuous Galerkin, etc.), for designing and analysing numerical schemes for diffusion models. In this paper, we study the GDM for a general stochastic evolution problem based on a Leray--Lions type operator. The problem contains the stochastic $p$-Laplace equatio…
▽ More
*The gradient discretisation method (GDM) is a generic framework, covering many classical methods (Finite Elements, Finite Volumes, Discontinuous Galerkin, etc.), for designing and analysing numerical schemes for diffusion models. In this paper, we study the GDM for a general stochastic evolution problem based on a Leray--Lions type operator. The problem contains the stochastic $p$-Laplace equation as a particular case. The convergence of the Gradient Scheme (GS) solutions is proved by using Discrete Functional Analysis techniques, Skorohod theorem and the Kolmogorov test. In particular, we provide an independent proof of the existence of weak martingale solutions for the problem. In this way, we lay foundations and provide techniques for proving convergence of the GS approximating stochastic partial differential equations.
△ Less
Submitted 29 December, 2020; v1 submitted 24 April, 2020;
originally announced April 2020.
-
On Unbalanced Optimal Transport: An Analysis of Sinkhorn Algorithm
Authors:
Khiem Pham,
Khang Le,
Nhat Ho,
Tung Pham,
Hung Bui
Abstract:
We provide a computational complexity analysis for the Sinkhorn algorithm that solves the entropic regularized Unbalanced Optimal Transport (UOT) problem between two measures of possibly different masses with at most $n$ components. We show that the complexity of the Sinkhorn algorithm for finding an $\varepsilon$-approximate solution to the UOT problem is of order…
▽ More
We provide a computational complexity analysis for the Sinkhorn algorithm that solves the entropic regularized Unbalanced Optimal Transport (UOT) problem between two measures of possibly different masses with at most $n$ components. We show that the complexity of the Sinkhorn algorithm for finding an $\varepsilon$-approximate solution to the UOT problem is of order $\widetilde{\mathcal{O}}(n^2/ \varepsilon)$, which is near-linear time. To the best of our knowledge, this complexity is better than the complexity of the Sinkhorn algorithm for solving the Optimal Transport (OT) problem, which is of order $\widetilde{\mathcal{O}}(n^2/\varepsilon^2)$. Our proof technique is based on the geometric convergence of the Sinkhorn updates to the optimal dual solution of the entropic regularized UOT problem and some properties of the primal solution. It is also different from the proof for the complexity of the Sinkhorn algorithm for approximating the OT problem since the UOT solution does not have to meet the marginal constraints.
△ Less
Submitted 18 November, 2020; v1 submitted 9 February, 2020;
originally announced February 2020.
-
L-space knots with tunnel number >1 by experiment
Authors:
Chris Anderson,
Kenneth L. Baker,
Xinghua Gao,
Marc Kegel,
Khanh Le,
Kyle Miller,
Sinem Onaran,
Geoffrey Sangston,
Samuel Tripp,
Adam Wood,
Ana Wright
Abstract:
In Dunfield's catalog of the hyperbolic manifolds in the SnapPy census which are complements of L-space knots in $S^3$, we determine that $22$ have tunnel number $2$ while the remaining all have tunnel number $1$. Notably, these $22$ manifolds contain $9$ asymmetric L-space knot complements. Furthermore, using SnapPy and KLO we find presentations of these $22$ knots as closures of positive braids…
▽ More
In Dunfield's catalog of the hyperbolic manifolds in the SnapPy census which are complements of L-space knots in $S^3$, we determine that $22$ have tunnel number $2$ while the remaining all have tunnel number $1$. Notably, these $22$ manifolds contain $9$ asymmetric L-space knot complements. Furthermore, using SnapPy and KLO we find presentations of these $22$ knots as closures of positive braids that realize the Morton-Franks-Williams bound on braid index. The smallest of these has genus $12$ and braid index $4$.
△ Less
Submitted 18 January, 2021; v1 submitted 2 September, 2019;
originally announced September 2019.
-
The gradient discretisation method for slow and fast diffusion porous media equations
Authors:
Jerome Droniou,
Kim-Ngan Le
Abstract:
The gradient discretisation method (GDM) is a generic framework for designing and analysing numerical schemes for diffusion models. In this paper, we study the GDM for the porous medium equation, including fast diffusion and slow diffusion models, and a concentration-dependent diffusion tensor. Using discrete functional analysis techniques, we establish a strong $L^2$-convergence of the approximat…
▽ More
The gradient discretisation method (GDM) is a generic framework for designing and analysing numerical schemes for diffusion models. In this paper, we study the GDM for the porous medium equation, including fast diffusion and slow diffusion models, and a concentration-dependent diffusion tensor. Using discrete functional analysis techniques, we establish a strong $L^2$-convergence of the approximate gradients and a uniform-in-time convergence for the approximate solution, without assuming non-physical regularity assumptions on the data or continuous solution. Being established in the generic GDM framework, these results apply to a variety of numerical methods, such as finite volume, (mass-lumped) finite elements, etc. The theoretical results are illustrated, in both fast and slow diffusion regimes, by numerical tests based on two methods that fit the GDM framework: mass-lumped conforming $\mathbb{P}_1$ finite elements and the Hybrid Mimetic Mixed method.
△ Less
Submitted 31 March, 2020; v1 submitted 5 May, 2019;
originally announced May 2019.
-
Well-posedness and nonsmooth Lyapunov pairs for state-dependent maximal monotone differential inclusions
Authors:
Ba Khiet Le
Abstract:
In this paper, we introduce for the first time a class of state-dependent maximal monotone differential inclusions. Then the existence and uniqueness of solutions are obtained by using an implicit discretization scheme and a kind of hypo-monotonicity assumption respectively. In addition, a characterization for nonsmooth Lyapunov pairs associated with such systems is provided. Our result can be app…
▽ More
In this paper, we introduce for the first time a class of state-dependent maximal monotone differential inclusions. Then the existence and uniqueness of solutions are obtained by using an implicit discretization scheme and a kind of hypo-monotonicity assumption respectively. In addition, a characterization for nonsmooth Lyapunov pairs associated with such systems is provided. Our result can be applied to study state-dependent sweeping processes and Lur'e dynamical systems. It is new even the involved maximal monotone operators depend only on the time.
△ Less
Submitted 5 March, 2019;
originally announced March 2019.
-
Lur'e dynamical systems with state-dependent set-valued feedback
Authors:
Ba Khiet Le
Abstract:
Using a new implicit discretization scheme, we study in this paper the existence and uniqueness of strong solutions for a class of Lur'e dynamical systems where the set-valued feedback depends on both time and state. This work is a generalization of \cite{abc} where the time-dependent set-valued feedback is considered to acquire only weak solutions. Obviously, strong solutions and implicit discret…
▽ More
Using a new implicit discretization scheme, we study in this paper the existence and uniqueness of strong solutions for a class of Lur'e dynamical systems where the set-valued feedback depends on both time and state. This work is a generalization of \cite{abc} where the time-dependent set-valued feedback is considered to acquire only weak solutions. Obviously, strong solutions and implicit discretization scheme are nice properties, especially for numerical simulation. We also provide some conditions such that the solutions are exponentially attractive. The obtained results can be used to study the time-varying Lur'e systems with errors in data. Our result is new even the set-valued feedback depends only on the time.
△ Less
Submitted 5 March, 2019;
originally announced March 2019.
-
A semidiscrete finite element approximation of a time-fractional Fokker-Planck equation with nonsmooth initial data
Authors:
Kim Ngan Le,
William McLean,
Kassem Mustapha
Abstract:
We present a new stability and convergence analysis for the spatial discretization of a time-fractional Fokker--Planck equation in a convex polyhedral domain, using continuous, piecewise-linear, finite elements. The forcing may depend on time as well as on the spatial variables, and the initial data may have low regularity. Our analysis uses a novel sequence of energy arguments in combination with…
▽ More
We present a new stability and convergence analysis for the spatial discretization of a time-fractional Fokker--Planck equation in a convex polyhedral domain, using continuous, piecewise-linear, finite elements. The forcing may depend on time as well as on the spatial variables, and the initial data may have low regularity. Our analysis uses a novel sequence of energy arguments in combination with a generalized Gronwall inequality. Although this theory covers only the spatial discretization, we present numerical experiments with a fully discrete scheme employing a very small time step, and observe results consistent with the predicted convergence behavior.
△ Less
Submitted 8 February, 2019;
originally announced February 2019.
-
Existence, uniqueness and regularity of the solution of the time-fractional Fokker-Planck equation with general forcing
Authors:
Kim-Ngan Le,
William McLean,
Martin Stynes
Abstract:
A time-fractional Fokker-Planck initial-boundary value problem is considered, with differential operator $u_t-\nabla\cdot(\partial_t^{1-α}κ_α\nabla u-\textbf{F}\partial_t^{1-α}u)$, where $0<α<1$. The forcing function $\textbf{F} = \textbf{F}(t,x)$, which is more difficult to analyse than the case $\textbf{F}=\textbf{F}(x)$ investigated previously by other authors. The spatial domain…
▽ More
A time-fractional Fokker-Planck initial-boundary value problem is considered, with differential operator $u_t-\nabla\cdot(\partial_t^{1-α}κ_α\nabla u-\textbf{F}\partial_t^{1-α}u)$, where $0<α<1$. The forcing function $\textbf{F} = \textbf{F}(t,x)$, which is more difficult to analyse than the case $\textbf{F}=\textbf{F}(x)$ investigated previously by other authors. The spatial domain $Ω\subset\mathbb{R}^d$, where $d\ge 1$, has a smooth boundary. Existence, uniqueness and regularity of a mild solution $u$ is proved under the hypothesis that the initial data $u_0$ lies in $L^2(Ω)$. For $1/2<α<1$ and $u_0\in H^2(Ω)\cap H_0^1(Ω)$, it is shown that $u$ becomes a classical solution of the problem. Estimates of time derivatives of the classical solution are derived---these are known to be needed in numerical analyses of this problem.
△ Less
Submitted 7 February, 2019;
originally announced February 2019.
-
Existence of a unique solution and invariant measures for the stochastic Landau--Lifshitz--Bloch equation
Authors:
Zdzislaw Brzeźniak,
Beniamin Goldys,
Kim Ngan Le
Abstract:
The Landau--Lifshitz--Bloch equation perturbed by a space-dependent noise was proposed in Garanin 1991 as a model for evolution of spins in ferromagnatic materials at the full range of temperatures, including the temperatures higher than the Curie temperature. In the case of a ferromagnet filling a bounded domain $D\subset \mathbb R^d$, $d=1,2,3$, we show the existence of strong (in the sense of P…
▽ More
The Landau--Lifshitz--Bloch equation perturbed by a space-dependent noise was proposed in Garanin 1991 as a model for evolution of spins in ferromagnatic materials at the full range of temperatures, including the temperatures higher than the Curie temperature. In the case of a ferromagnet filling a bounded domain $D\subset \mathbb R^d$, $d=1,2,3$, we show the existence of strong (in the sense of PDEs) martingale solutions. Furthermore, in cases $d=1,2$ we prove uniqueness of pathwise solutions and the existence of invariant measures.
△ Less
Submitted 9 January, 2019;
originally announced January 2019.
-
A stochastic sewing lemma and applications
Authors:
Khoa Lê
Abstract:
We introduce a stochastic version of Gubinelli's sewing lemma, providing a sufficient condition for the convergence in moments of some random Riemann sums. Compared with the deterministic sewing lemma, adaptiveness is required and the regularity restriction is improved by a half. The limiting process exhibits a Doob-Meyer-type decomposition. Relations with Itô calculus are established. To illustra…
▽ More
We introduce a stochastic version of Gubinelli's sewing lemma, providing a sufficient condition for the convergence in moments of some random Riemann sums. Compared with the deterministic sewing lemma, adaptiveness is required and the regularity restriction is improved by a half. The limiting process exhibits a Doob-Meyer-type decomposition. Relations with Itô calculus are established. To illustrate further potential applications, we use the stochastic sewing lemma in studying stochastic differential equations driven by Brownian motions or fractional Brownian motions with irregulardrifts.
△ Less
Submitted 13 March, 2020; v1 submitted 24 October, 2018;
originally announced October 2018.
-
Joint Hölder continuity of parabolic Anderson model
Authors:
Yaozhong Hu,
Khoa Lê
Abstract:
We show that the random field solution to the parabolic Anderson equation $(\partial_t-\frac12 Δ)u=u\diamond \dot{W}$ is jointly Hölder continuous in space and time.
We show that the random field solution to the parabolic Anderson equation $(\partial_t-\frac12 Δ)u=u\diamond \dot{W}$ is jointly Hölder continuous in space and time.
△ Less
Submitted 26 September, 2018;
originally announced September 2018.
-
A combined GDM--ELLAM--MMOC scheme for advection dominated PDEs
Authors:
Hanz Martin Cheng,
Jérôme Droniou,
Kim-Ngan Le
Abstract:
We propose a combination of the Eulerian Lagrangian Localised Adjoint Method (ELLAM) and the Modified Method of Characteristics (MMOC) for time-dependent advection-domina\-ted PDEs. The combined scheme, so-called GEM scheme, takes advantages of both ELLAM scheme (mass conservation) and MMOC scheme (easier computations), while at the same time avoids their disadvantages (respectively, harder tracki…
▽ More
We propose a combination of the Eulerian Lagrangian Localised Adjoint Method (ELLAM) and the Modified Method of Characteristics (MMOC) for time-dependent advection-domina\-ted PDEs. The combined scheme, so-called GEM scheme, takes advantages of both ELLAM scheme (mass conservation) and MMOC scheme (easier computations), while at the same time avoids their disadvantages (respectively, harder tracking around the injection regions, and loss of mass).
We present a precise analysis of mass conservation properties for these three schemes, and after achieving global mass balance, an adjustment yielding local volume conservation is then proposed. Numerical results for all three schemes are then compared, illustrating the advantages of the GEM scheme.
A convergence result of the MMOC scheme, motivated by our previous work on the convergence of ELLAM schemes, is provided, which can be extended to obtain the convergence of GEM scheme.
△ Less
Submitted 13 October, 2021; v1 submitted 15 May, 2018;
originally announced May 2018.