-
Hecke Equivariance of Divisor Lifting with respect to Sesquiharmonic Maass Forms
Authors:
Daeyeol Jeon,
Soon-Yi Kang,
Chang Heon Kim
Abstract:
We investigate the properties of Hecke operator for sesquiharmonic Maass forms. We begin by proving Hecke equivariance of the divisor lifting with respect to sesquiharmonic Mass functions, which maps an integral weight meromorphic modular form to the holomorphic part of the Fourier expansion of a weight 2 sesquiharmonic Maass form. Using this Hecke equivariance, we show that the sesquiharmonic Maa…
▽ More
We investigate the properties of Hecke operator for sesquiharmonic Maass forms. We begin by proving Hecke equivariance of the divisor lifting with respect to sesquiharmonic Mass functions, which maps an integral weight meromorphic modular form to the holomorphic part of the Fourier expansion of a weight 2 sesquiharmonic Maass form. Using this Hecke equivariance, we show that the sesquiharmonic Maass functions, whose images under the hyperbolic Laplace operator are the Faber polynomials $J_n$ of the $j$-function, form a Hecke system analogous to $J_n$. By combining the Hecke equivariance of the divisor lifting with that of the Borcherds isomorphism, we extend Matsusaka's finding on the twisted traces of sesquiharmonic Maass functions.
△ Less
Submitted 31 July, 2024;
originally announced July 2024.
-
Multiscale Modeling Framework using Element-based Galerkin Methods for Moist Atmospheric Limited-Area Simulations
Authors:
Soonpil Kang,
James F. Kelly,
Anthony P. Austin,
Francis X. Giraldo
Abstract:
This paper presents a multiscale modeling framework (MMF) to model moist atmospheric limited-area weather. The MMF resolves large-scale convection using a coarse grid while simultaneously resolving local features through numerous fine local grids and coupling them seamlessly. Both large- and small-scale processes are modeled using the compressible Navier-Stokes equations within the Nonhydrostatic…
▽ More
This paper presents a multiscale modeling framework (MMF) to model moist atmospheric limited-area weather. The MMF resolves large-scale convection using a coarse grid while simultaneously resolving local features through numerous fine local grids and coupling them seamlessly. Both large- and small-scale processes are modeled using the compressible Navier-Stokes equations within the Nonhydrostatic Unified Model of the Atmosphere (NUMA), and they are discretized using a continuous element-based Galerkin method (spectral elements) with high-order basis functions. Consequently, the large-scale and small-scale models share the same dynamical core but have the flexibility to be adjusted individually. The proposed MMF method is tested in 2D and 3D idealized limited-area weather problems involving storm clouds produced by squall line and supercell simulations. The MMF numerical results showed enhanced representation of cloud processes compared to the coarse model.
△ Less
Submitted 11 May, 2024;
originally announced July 2024.
-
Categorification of quantum Borcherds-Bozec algebras
Authors:
Seok-Jin Kang,
Young Rock Kim,
Bolun Tong
Abstract:
We categorify the quantum Borcherds-Bozec algebras by constructing their associated Khovanov-Lauda-Rouquier algebras.
We categorify the quantum Borcherds-Bozec algebras by constructing their associated Khovanov-Lauda-Rouquier algebras.
△ Less
Submitted 26 June, 2024;
originally announced June 2024.
-
Fast and Certifiable Trajectory Optimization
Authors:
Shucheng Kang,
Xiaoyang Xu,
Jay Sarva,
Ling Liang,
Heng Yang
Abstract:
We propose semidefinite trajectory optimization (STROM), a framework that computes fast and certifiably optimal solutions for nonconvex trajectory optimization problems defined by polynomial objectives and constraints. STROM employs sparse second-order Lasserre's hierarchy to generate semidefinite program (SDP) relaxations of trajectory optimization. Different from existing tools (e.g., YALMIP and…
▽ More
We propose semidefinite trajectory optimization (STROM), a framework that computes fast and certifiably optimal solutions for nonconvex trajectory optimization problems defined by polynomial objectives and constraints. STROM employs sparse second-order Lasserre's hierarchy to generate semidefinite program (SDP) relaxations of trajectory optimization. Different from existing tools (e.g., YALMIP and SOSTOOLS in Matlab), STROM generates chain-like multiple-block SDPs with only positive semidefinite (PSD) variables. Moreover, STROM does so two orders of magnitude faster. Underpinning STROM is cuADMM, the first ADMM-based SDP solver implemented in CUDA and runs in GPUs. cuADMM builds upon the symmetric Gauss-Seidel ADMM algorithm and leverages GPU parallelization to speedup solving sparse linear systems and projecting onto PSD cones. In five trajectory optimization problems (inverted pendulum, cart-pole, vehicle landing, flying robot, and car back-in), cuADMM computes optimal trajectories (with certified suboptimality below 1%) in minutes (when other solvers take hours or run out of memory) and seconds (when others take minutes). Further, when warmstarted by data-driven initialization in the inverted pendulum problem, cuADMM delivers real-time performance: providing certifiably optimal trajectories in 0.66 seconds despite the SDP has 49,500 variables and 47,351 constraints.
△ Less
Submitted 11 June, 2024; v1 submitted 9 June, 2024;
originally announced June 2024.
-
Distributionally Risk-Receptive and Robust Multistage Stochastic Integer Programs and Two-player Interdiction Games with(out) Decision-Dependent Uncertainty
Authors:
Sumin Kang,
Manish Bansal
Abstract:
In this paper, we study distributionally risk-receptive and distributionally robust (or risk-averse) multistage stochastic mixed-integer programs (denoted by DRR- and DRA-MSIPs). These frameworks are useful for optimization problems under uncertainty where the focus is on analyzing outcomes based on multiple decision-makers' differing perspectives, such as interdiction problems that are attacker-d…
▽ More
In this paper, we study distributionally risk-receptive and distributionally robust (or risk-averse) multistage stochastic mixed-integer programs (denoted by DRR- and DRA-MSIPs). These frameworks are useful for optimization problems under uncertainty where the focus is on analyzing outcomes based on multiple decision-makers' differing perspectives, such as interdiction problems that are attacker-defender games having non-cooperative players. We present cutting plane-based and reformulation-based approaches for solving DRR- and DRA-MSIPs without and with decision-dependent uncertainty to optimality. We show that these approaches are finitely convergent with probability one. Furthermore, we introduce generalizations of DRR- and DRA-MSIPs by presenting multistage stochastic disjunctive programs and algorithms for solving them. To assess the performance of the algorithms for DRR- and DRA-MSIPs, we consider instances of distributionally ambiguous multistage maximum flow and facility location interdiction problems that are important in their own right. Based on our computational results, we observe that the cutting plane-based approaches are 2800% and 2410% (on average) faster than the reformulation-based approaches for the foregoing instances with distributional risk-aversion and risk-receptiveness, respectively. Additionally, we conducted out-of-sample tests to showcase the significance of the DRR framework in revealing network vulnerabilities and also in mitigating the impact of data corruption.
△ Less
Submitted 7 June, 2024;
originally announced June 2024.
-
Cables of the figure-eight knot via real Frøyshov invariants
Authors:
Sungkyung Kang,
JungHwan Park,
Masaki Taniguchi
Abstract:
We prove that the $(2n,1)$-cable of the figure-eight knot is not smoothly slice when $n$ is odd, by using the real Seiberg-Witten Frøyshov invariant of Konno-Miyazawa-Taniguchi. For the computation, we develop an $O(2)$-equivariant version of the lattice homotopy type, originally introduced by Dai-Sasahira-Stoffregen. This enables us to compute the real Seiberg-Witten Floer homotopy type for a cer…
▽ More
We prove that the $(2n,1)$-cable of the figure-eight knot is not smoothly slice when $n$ is odd, by using the real Seiberg-Witten Frøyshov invariant of Konno-Miyazawa-Taniguchi. For the computation, we develop an $O(2)$-equivariant version of the lattice homotopy type, originally introduced by Dai-Sasahira-Stoffregen. This enables us to compute the real Seiberg-Witten Floer homotopy type for a certain class of knots. Additionally, we present some computations of Miyazawa's real framed Seiberg-Witten invariant for 2-knots.
△ Less
Submitted 15 May, 2024;
originally announced May 2024.
-
Perfect basis theory for quantum Borcherds-Bozec algebras
Authors:
Zhaobing Fan,
Shaolong Han,
Seok-Jin Kang,
Young Rock Kim
Abstract:
In this paper, we develop the perfect basis theory for quantum Borcherds-Bozec algebras $U_{q}(\mathfrak g)$ and their irreducible highest weight modules $V(λ)$. We show that the lower perfect graph (resp. upper perfect graph) of every lower perfect basis (resp. upper perfect basis) of $U_{q}^{-}(\mathfrak g)$ (resp. $V(λ)$) is isomorphic to the crystal $B(\infty)$ (resp. $B(λ)$).
In this paper, we develop the perfect basis theory for quantum Borcherds-Bozec algebras $U_{q}(\mathfrak g)$ and their irreducible highest weight modules $V(λ)$. We show that the lower perfect graph (resp. upper perfect graph) of every lower perfect basis (resp. upper perfect basis) of $U_{q}^{-}(\mathfrak g)$ (resp. $V(λ)$) is isomorphic to the crystal $B(\infty)$ (resp. $B(λ)$).
△ Less
Submitted 9 May, 2024;
originally announced May 2024.
-
Invariant splitting principles for the Lipshitz--Ozsváth--Thurston correspondence
Authors:
Gary Guth,
Sungkyung Kang
Abstract:
We prove that the Lipshitz-Ozsváth-Thurston correspondence between extended type D structures of knot complements and $\mathbb{F}[U, V]/(UV)$ knot Floer complexes can be arranged so that $ι_K$-invariant splittings of knot Floer chain complexes correspond to $ι_{S^3 \setminus K}$-invariant splittings of bordered Floer homology of knot complements. For patterns satisfying the satellite extension pro…
▽ More
We prove that the Lipshitz-Ozsváth-Thurston correspondence between extended type D structures of knot complements and $\mathbb{F}[U, V]/(UV)$ knot Floer complexes can be arranged so that $ι_K$-invariant splittings of knot Floer chain complexes correspond to $ι_{S^3 \setminus K}$-invariant splittings of bordered Floer homology of knot complements. For patterns satisfying the satellite extension property, which include cabling patterns, this provides a novel way to compute the involutive knot Floer homology of satellites from that of their companions. As a topological application, we show that our results can be applied to construct infinitely many examples of exotic pairs of contractible 4-manifolds which remain exotic after one stabilization. Along the way, we also establish first order naturality of bordered Floer homology.
△ Less
Submitted 9 April, 2024;
originally announced April 2024.
-
A new Young wall realization of $B(λ)$ and $B(\infty)$
Authors:
Zhaobing Fan,
Shaolong Han,
Seok-Jin Kang,
Young Rock Kim
Abstract:
Using new combinatorics of Young walls, we give a new construction of the arbitrary level highest weight crystal $B(λ)$ for the quantum affine algebras of types $A^{(2)}_{2n}$, $D^{(2)}_{n+1}$, $A^{(2)}_{2n-1}$, $D^{(1)}_n$, $B^{(1)}_n$ and $C^{(1)}_n$. We show that the crystal consisting of reduced Young walls is isomorphic to the crystal $B(λ)$. Moreover, we provide a new realization of the crys…
▽ More
Using new combinatorics of Young walls, we give a new construction of the arbitrary level highest weight crystal $B(λ)$ for the quantum affine algebras of types $A^{(2)}_{2n}$, $D^{(2)}_{n+1}$, $A^{(2)}_{2n-1}$, $D^{(1)}_n$, $B^{(1)}_n$ and $C^{(1)}_n$. We show that the crystal consisting of reduced Young walls is isomorphic to the crystal $B(λ)$. Moreover, we provide a new realization of the crystal $B(\infty)$ in terms of reduced virtual Young walls and reduced extended Young walls.
△ Less
Submitted 16 March, 2024;
originally announced March 2024.
-
Texture Edge detection by Patch consensus (TEP)
Authors:
Guangyu Cui,
Sung Ha Kang
Abstract:
We propose Texture Edge detection using Patch consensus (TEP) which is a training-free method to detect the boundary of texture. We propose a new simple way to identify the texture edge location, using the consensus of segmented local patch information. While on the boundary, even using local patch information, the distinction between textures are typically not clear, but using neighbor consensus…
▽ More
We propose Texture Edge detection using Patch consensus (TEP) which is a training-free method to detect the boundary of texture. We propose a new simple way to identify the texture edge location, using the consensus of segmented local patch information. While on the boundary, even using local patch information, the distinction between textures are typically not clear, but using neighbor consensus give a clear idea of the boundary. We utilize local patch, and its response against neighboring regions, to emphasize the similarities and the differences across different textures. The step of segmentation of response further emphasizes the edge location, and the neighborhood voting gives consensus and stabilize the edge detection. We analyze texture as a stationary process to give insight into the patch width parameter verses the quality of edge detection. We derive the necessary condition for textures to be distinguished, and analyze the patch width with respect to the scale of textures. Various experiments are presented to validate the proposed model.
△ Less
Submitted 16 March, 2024;
originally announced March 2024.
-
Young wall models for the level 1 highest weight and Fock space crystals of $U_q(E_6^{(2)})$ and $U_q(F_4^{(1)})$
Authors:
Shaolong Han,
Yuanfeng Jin,
Seok-Jin Kang,
Duncan Laurie
Abstract:
In this paper we construct Young wall models for the level $1$ highest weight and Fock space crystals of quantum affine algebras in types $E_6^{(2)}$ and $F_4^{(1)}$. Our starting point in each case is a combinatorial realization for a certain level $1$ perfect crystal in terms of Young columns. Then using energy functions and affine energy functions we define the notions of reduced and proper You…
▽ More
In this paper we construct Young wall models for the level $1$ highest weight and Fock space crystals of quantum affine algebras in types $E_6^{(2)}$ and $F_4^{(1)}$. Our starting point in each case is a combinatorial realization for a certain level $1$ perfect crystal in terms of Young columns. Then using energy functions and affine energy functions we define the notions of reduced and proper Young walls, which model the highest weight and Fock space crystals respectively.
△ Less
Submitted 24 February, 2024;
originally announced February 2024.
-
Sparse Polynomial Optimization with Unbounded Sets
Authors:
Lei Huang,
Shucheng Kang,
Jie Wang,
Heng Yang
Abstract:
This paper considers sparse polynomial optimization with unbounded sets. When the problem possesses correlative sparsity, we propose a sparse homogenized Moment-SOS hierarchy with perturbations to solve it. The new hierarchy introduces one extra auxiliary variable for each variable clique according to the correlative sparsity pattern. Under the running intersection property, we prove that this hie…
▽ More
This paper considers sparse polynomial optimization with unbounded sets. When the problem possesses correlative sparsity, we propose a sparse homogenized Moment-SOS hierarchy with perturbations to solve it. The new hierarchy introduces one extra auxiliary variable for each variable clique according to the correlative sparsity pattern. Under the running intersection property, we prove that this hierarchy has asymptotic convergence. Furthermore, we provide two alternative sparse hierarchies to remove perturbations while preserving asymptotic convergence. As byproducts, new Positivstellensätze are obtained for sparse positive polynomials on unbounded sets. Extensive numerical experiments demonstrate the power of our approach in solving sparse polynomial optimization problems on unbounded sets with up to thousands of variables. Finally, we apply our approach to tackle two trajectory optimization problems (block-moving with minimum work and optimal control of Van der Pol).
△ Less
Submitted 28 January, 2024;
originally announced January 2024.
-
Highest weight modules over Borcherds-Bozec superalgebras and their character formula
Authors:
Zhaobing Fan,
Jiaqi Huang,
Seok-Jin Kang,
Yong-Su Shin
Abstract:
We present and prove the Weyl-Kac type character formula for the irreducible highest weight modules over Borcherds-Bozec superalgebras with dominant integral highest weights.
We present and prove the Weyl-Kac type character formula for the irreducible highest weight modules over Borcherds-Bozec superalgebras with dominant integral highest weights.
△ Less
Submitted 28 January, 2024;
originally announced January 2024.
-
Monostatic imaging of an extended target with MCMC sampling
Authors:
Jiho Hong,
Sangwoo Kang,
Mikyoung Lim
Abstract:
We consider the imaging of a planar extended target from far-field data under a monostatic measurement configuration, in which the data is measured by a single moving transducer, as frequently encountered in practical application. In this paper, we develop a Bayesian approach to recover the shape of the extended target with MCMC sampling, where a new shape basis selection is proposed based on the…
▽ More
We consider the imaging of a planar extended target from far-field data under a monostatic measurement configuration, in which the data is measured by a single moving transducer, as frequently encountered in practical application. In this paper, we develop a Bayesian approach to recover the shape of the extended target with MCMC sampling, where a new shape basis selection is proposed based on the shape derivative analysis for the measurement data. In order to optimize the center and radius of the initial disk, we use the monostatic sampling method for the center and the explicit scattered field expression for disks for the radius. Numerical simulations are presented to validate the proposed method.
△ Less
Submitted 8 December, 2023;
originally announced January 2024.
-
Improving reconstructions in nanotomography for homogeneous materials via mathematical optimization
Authors:
Sebastian Kreuz,
Benjamin Apeleo Zubiri,
Silvan Englisch,
Sung-Gyu Kang,
Rajaprakash Ramachandramoorthy,
Erdmann Spiecker,
Frauke Liers,
Jan Rolfes
Abstract:
Compressed sensing is an image reconstruction technique to achieve high-quality results from limited amount of data. In order to achieve this, it utilizes prior knowledge about the samples that shall be reconstructed. Focusing on image reconstruction in nanotomography, this work proposes enhancements by including additional problem-specific knowledge. In more detail, we propose further classes of…
▽ More
Compressed sensing is an image reconstruction technique to achieve high-quality results from limited amount of data. In order to achieve this, it utilizes prior knowledge about the samples that shall be reconstructed. Focusing on image reconstruction in nanotomography, this work proposes enhancements by including additional problem-specific knowledge. In more detail, we propose further classes of algebraic inequalities that are added to the compressed sensing model. The first consists in a valid upper bound on the pixel brightness. It only exploits general information about the projections and is thus applicable to a broad range of reconstruction problems. The second class is applicable whenever the sample material is of roughly homogeneous composition. The model favors a constant density and penalizes deviations from it. The resulting mathematical optimization models are algorithmically tractable and can be solved to global optimality by state-of-the-art available implementations of interior point methods. In order to evaluate the novel models, obtained results are compared to existing image reconstruction methods, tested on simulated and experimental data sets. The experimental data comprise one 360° electron tomography tilt series of a macroporous zeolite particle and one absorption contrast nano X-ray computed tomography (nano-CT) data set of a copper microlattice structure. The enriched models are optimized quickly and show improved reconstruction quality, outperforming the existing models. Promisingly, our approach yields superior reconstruction results, particularly when information about the samples is available for a small number of tilt angles only
△ Less
Submitted 9 December, 2023;
originally announced December 2023.
-
Fourier Features for Identifying Differential Equations (FourierIdent)
Authors:
Mengyi Tang,
Hao Liu,
Wenjing Liao,
Sung Ha Kang
Abstract:
We investigate the benefits and challenges of utilizing the frequency information in differential equation identification. Solving differential equations and Fourier analysis are closely related, yet there is limited work in exploring this connection in the identification of differential equations. Given a single realization of the differential equation perturbed by noise, we aim to identify the u…
▽ More
We investigate the benefits and challenges of utilizing the frequency information in differential equation identification. Solving differential equations and Fourier analysis are closely related, yet there is limited work in exploring this connection in the identification of differential equations. Given a single realization of the differential equation perturbed by noise, we aim to identify the underlying differential equation governed by a linear combination of linear and nonlinear differential and polynomial terms in the frequency domain. This is challenging due to large magnitudes and sensitivity to noise. We introduce a Fourier feature denoising, and define the meaningful data region and the core regions of features to reduce the effect of noise in the frequency domain. We use Subspace Pursuit on the core region of the time derivative feature, and introduce a group trimming step to refine the support. We further introduce a new energy based on the core regions of features for coefficient identification. Utilizing the core regions of features serves two critical purposes: eliminating the low-response regions dominated by noise, and enhancing the accuracy in coefficient identification. The proposed method is tested on various differential equations with linear, nonlinear, and high-order derivative feature terms. Our results demonstrate the advantages of the proposed method, particularly on complex and highly corrupted datasets.
△ Less
Submitted 28 November, 2023;
originally announced November 2023.
-
Central extensions of higher groups: Green-Schwarz mechanism and 2-connections
Authors:
Monica Jinwoo Kang,
Sungkyung Kang
Abstract:
We study the smooth $2$-group structure arising in the presence of quantum field theory with one-form symmetry. We acquire $2$-group structures obtained by a central extension of the zero-form symmetry by the one-form symmetry. We determine that the existence of a $2$-group structure is guaranteed by Chern--Simons levels. We further verify how we will be able to provide a fix to the current $2$-gr…
▽ More
We study the smooth $2$-group structure arising in the presence of quantum field theory with one-form symmetry. We acquire $2$-group structures obtained by a central extension of the zero-form symmetry by the one-form symmetry. We determine that the existence of a $2$-group structure is guaranteed by Chern--Simons levels. We further verify how we will be able to provide a fix to the current $2$-group problems by using the bibundle model. We outline the principal $2$-connection theory with respect to such $2$-group and compare it with the ansatz obtained from the Green--Schwarz mechanism. We further propose the existence of smooth $\infty$-group symmetries in quantum field theory.
△ Less
Submitted 27 November, 2023; v1 submitted 24 November, 2023;
originally announced November 2023.
-
A Performance Study of Horizontally Explicit Vertically Implicit (HEVI) Time-Integrators for Non-Hydrostatic Atmospheric Models
Authors:
Francis X. Giraldo,
Felipe Augusto Ventura de Braganca Alves,
James F. Kelly,
Soonpil Kang,
P. Alex Reiencke
Abstract:
We conduct a thorough study of different forms of horizontally explicit and vertically implicit (HEVI) time-integration strategies for the compressible Euler equations on spherical domains typical of nonhydrostatic global atmospheric applications. We compare the computational time and complexity of two nonlinear variants (NHEVI-GMRES and NHEVI-LU) and a linear variant (LHEVI). We report on the per…
▽ More
We conduct a thorough study of different forms of horizontally explicit and vertically implicit (HEVI) time-integration strategies for the compressible Euler equations on spherical domains typical of nonhydrostatic global atmospheric applications. We compare the computational time and complexity of two nonlinear variants (NHEVI-GMRES and NHEVI-LU) and a linear variant (LHEVI). We report on the performance of these three variants for a number of additive Runge-Kutta Methods ranging in order of accuracy from second through fifth, and confirm the expected order of accuracy of the HEVI methods for each time-integrator. To gauge the maximum usable time-step of each HEVI method, we run simulations of a nonhydrostatic baroclinic instability for 100 days and then use this time-step to compare the time-to-solution of each method. The results show that NHEVI-LU is 2x faster than NHEVI-GMRES, and LHEVI is 5x faster than NHEVI-LU, for the idealized cases tested. The baroclinic instability and inertia-gravity wave simulations indicate that the optimal choice of time-integrator is LHEVI with either second or third order schemes, as both schemes yield similar time to solution and relative L2 error at their maximum usable time-steps. In the future, we will report on whether these results hold for more complex problems using, e.g., real atmospheric data and/or a higher model top typical of space weather applications.
△ Less
Submitted 19 November, 2023;
originally announced November 2023.
-
Doubled Disks and Satellite Surfaces
Authors:
Gary Guth,
Kyle Hayden,
Sungkyung Kang,
JungHwan Park
Abstract:
Conjecturally, a knot is slice if and only if its positive Whitehead double is slice. We consider an analogue of this conjecture for slice disks in the four-ball: two slice disks of a knot are smoothly isotopic if and only if their positive Whitehead doubles are smoothly isotopic. We provide evidence for this conjecture, using a range of techniques. More generally, we consider when isotopy obstruc…
▽ More
Conjecturally, a knot is slice if and only if its positive Whitehead double is slice. We consider an analogue of this conjecture for slice disks in the four-ball: two slice disks of a knot are smoothly isotopic if and only if their positive Whitehead doubles are smoothly isotopic. We provide evidence for this conjecture, using a range of techniques. More generally, we consider when isotopy obstructions persist under satellite operations. In particular, we show that obstructions coming from knot Floer homology, Seiberg-Witten theory, and Khovanov homology often behave well under satellite operations.
We apply these strategies to give a systematic method for constructing vast numbers of exotic disks in the four-ball, including the first infinite family of pairwise exotic slice disks. These same techniques are then upgraded to produce exotic disks that remain exotic after any prescribed number of internal stabilizations. Finally, we show that the branched double covers of certain stably-exotic disks become diffeomorphic after a single stabilization with $S^2 \times S^2$, hence stabilizing them yields exotic surfaces that have diffeomorphic branched covers.
△ Less
Submitted 30 October, 2023;
originally announced October 2023.
-
Enhancing Low-Order Discontinuous Galerkin Methods with Neural Ordinary Differential Equations for Compressible Navier--Stokes Equations
Authors:
Shinhoo Kang,
Emil M. Constantinescu
Abstract:
The growing computing power over the years has enabled simulations to become more complex and accurate. While immensely valuable for scientific discovery and problem-solving, however, high-fidelity simulations come with significant computational demands. As a result, it is common to run a low-fidelity model with a subgrid-scale model to reduce the computational cost, but selecting the appropriate…
▽ More
The growing computing power over the years has enabled simulations to become more complex and accurate. While immensely valuable for scientific discovery and problem-solving, however, high-fidelity simulations come with significant computational demands. As a result, it is common to run a low-fidelity model with a subgrid-scale model to reduce the computational cost, but selecting the appropriate subgrid-scale models and tuning them are challenging. We propose a novel method for learning the subgrid-scale model effects when simulating partial differential equations augmented by neural ordinary differential operators in the context of discontinuous Galerkin (DG) spatial discretization. Our approach learns the missing scales of the low-order DG solver at a continuous level and hence improves the accuracy of the low-order DG approximations as well as accelerates the filtered high-order DG simulations with a certain degree of precision. We demonstrate the performance of our approach through multidimensional Taylor-Green vortex examples at different Reynolds numbers and times, which cover laminar, transitional, and turbulent regimes. The proposed method not only reconstructs the subgrid-scale from the low-order (1st-order) approximation but also speeds up the filtered high-order DG (6th-order) simulation by two orders of magnitude.
△ Less
Submitted 30 January, 2024; v1 submitted 29 October, 2023;
originally announced October 2023.
-
Diffusion laws select boundary conditions
Authors:
Jaywan Chung,
Seungmin Kang,
Ho-Youn Kim,
Yong-Jung Kim
Abstract:
The choice of boundary condition makes an essential difference in the solution structure of diffusion equations. The Dirichlet and Neumann boundary conditions and their combination have been the most used, but their legitimacy has been disputed. We show that the diffusion laws may select boundary conditions by themselves, and through this, we clarify the meaning of boundary conditions. To do that…
▽ More
The choice of boundary condition makes an essential difference in the solution structure of diffusion equations. The Dirichlet and Neumann boundary conditions and their combination have been the most used, but their legitimacy has been disputed. We show that the diffusion laws may select boundary conditions by themselves, and through this, we clarify the meaning of boundary conditions. To do that we extend the domain with a boundary into the whole space by giving a small diffusivity $\eps>0$ outside the domain. Then, we show that the boundary condition turns out to be Neumann or Dirichlet as $\eps\to0$ depending on the choice of a heterogeneous diffusion law. These boundary conditions are interpreted in terms of a microscopic-scale random walk model.
△ Less
Submitted 1 August, 2023;
originally announced August 2023.
-
Fick's las selects the Neumann boundary condition
Authors:
Danielle Hilhorst,
Seung-Min Kang,
Ho-Youn Kim,
Yong-Jung Kim
Abstract:
We study the appearance of a boundary condition along an interface between two regions, one with constant diffusivity $1$ and the other with diffusivity $\eps>0$, when $\eps\to0$. In particular, we take Fick's diffusion law in a context of reaction-diffusion equation with bistable nonlinearity and show that the limit of the reaction-diffusion equation satisfies the homogeneous Neumann boundary con…
▽ More
We study the appearance of a boundary condition along an interface between two regions, one with constant diffusivity $1$ and the other with diffusivity $\eps>0$, when $\eps\to0$. In particular, we take Fick's diffusion law in a context of reaction-diffusion equation with bistable nonlinearity and show that the limit of the reaction-diffusion equation satisfies the homogeneous Neumann boundary condition along the interface. This problem is developed as an application of heterogeneous diffusion laws to study the geometry effect of domain.
△ Less
Submitted 10 August, 2023; v1 submitted 1 August, 2023;
originally announced August 2023.
-
Alder-type partition inequality at the general level
Authors:
Haein Cho,
Soon-Yi Kang,
Byungchan Kim
Abstract:
A Known Alder-type partition inequality of level $a$, which involves the second Rogers-Ramanujan identity when the level $a$ is 2, states that the number of partitions of $n$ into parts differing by at least $d$ with the smallest part being at least $a$ is greater than or equal to that of partitions of $n$ into parts congruent to $\pm a \pmod{d+3}$, excluding the part $d+3-a$. In this paper, we pr…
▽ More
A Known Alder-type partition inequality of level $a$, which involves the second Rogers-Ramanujan identity when the level $a$ is 2, states that the number of partitions of $n$ into parts differing by at least $d$ with the smallest part being at least $a$ is greater than or equal to that of partitions of $n$ into parts congruent to $\pm a \pmod{d+3}$, excluding the part $d+3-a$. In this paper, we prove that for all values of $d$ with a finite number of exceptions, an arbitrary level $a$ Alder-type partition inequality holds without requiring the exclusion of the part $d+3-a$ in the latter partition.
△ Less
Submitted 4 August, 2023; v1 submitted 26 July, 2023;
originally announced July 2023.
-
A general model-checking procedure for semiparametric accelerated failure time models
Authors:
Dongrak Choi,
Woojung Bae,
Jun Yan,
Sangwook Kang
Abstract:
We propose a set of goodness-of-fit tests for the semiparametric accelerated failure time (AFT) model, including an omnibus test, a link function test, and a functional form test. This set of tests is derived from a multi-parameter cumulative sum process shown to follow asymptotically a zero-mean Gaussian process. Its evaluation is based on the asymptotically equivalent perturbed version, which en…
▽ More
We propose a set of goodness-of-fit tests for the semiparametric accelerated failure time (AFT) model, including an omnibus test, a link function test, and a functional form test. This set of tests is derived from a multi-parameter cumulative sum process shown to follow asymptotically a zero-mean Gaussian process. Its evaluation is based on the asymptotically equivalent perturbed version, which enables both graphical and numerical evaluations of the assumed AFT model. Empirical p-values are obtained using the Kolmogorov-type supremum test, which provides a reliable approach for estimating the significance of both proposed un-standardized and standardized test statistics. The proposed procedure is illustrated using the induced smoothed rank-based estimator but is directly applicable to other popular estimators such as non-smooth rank-based estimator or least-squares estimator.Our proposed methods are rigorously evaluated using extensive simulation experiments that demonstrate their effectiveness in maintaining a Type I error rate and detecting departures from the assumed AFT model in practical sample sizes and censoring rates. Furthermore, the proposed approach is applied to the analysis of the Primary Biliary Cirrhosis data, a widely studied dataset in survival analysis, providing further evidence of the practical usefulness of the proposed methods in real-world scenarios. To make the proposed methods more accessible to researchers, we have implemented them in the R package afttest, which is publicly available on the Comprehensive R Archieve Network.
△ Less
Submitted 19 May, 2023;
originally announced May 2023.
-
Topologically and rationally slice knots
Authors:
Jennifer Hom,
Sungkyung Kang,
JungHwan Park
Abstract:
A knot in $S^3$ is topologically slice if it bounds a locally flat disk in $B^4$. A knot in $S^3$ is rationally slice if it bounds a smooth disk in a rational homology ball. We prove that the smooth concordance group of topologically and rationally slice knots admits a $\mathbb{Z}^\infty$ subgroup. All previously known examples of knots that are both topologically and rationally slice were of orde…
▽ More
A knot in $S^3$ is topologically slice if it bounds a locally flat disk in $B^4$. A knot in $S^3$ is rationally slice if it bounds a smooth disk in a rational homology ball. We prove that the smooth concordance group of topologically and rationally slice knots admits a $\mathbb{Z}^\infty$ subgroup. All previously known examples of knots that are both topologically and rationally slice were of order two. As a direct consequence, it follows that there are infinitely many topologically slice knots that are strongly rationally slice but not slice.
△ Less
Submitted 13 April, 2023;
originally announced April 2023.
-
Group projected Subspace Pursuit for Identification of variable coefficient differential equations (GP-IDENT)
Authors:
Yuchen He,
Sung-Ha Kang,
Wenjing Liao,
Hao Liu,
Yingjie Liu
Abstract:
We propose an effective and robust algorithm for identifying partial differential equations (PDEs) with space-time varying coefficients from a single trajectory of noisy observations. Identifying unknown differential equations from noisy observations is a difficult task, and it is even more challenging with space and time varying coefficients in the PDE. The proposed algorithm, GP-IDENT, has three…
▽ More
We propose an effective and robust algorithm for identifying partial differential equations (PDEs) with space-time varying coefficients from a single trajectory of noisy observations. Identifying unknown differential equations from noisy observations is a difficult task, and it is even more challenging with space and time varying coefficients in the PDE. The proposed algorithm, GP-IDENT, has three ingredients: (i) we use B-spline bases to express the unknown space and time varying coefficients, (ii) we propose Group Projected Subspace Pursuit (GPSP) to find a sequence of candidate PDEs with various levels of complexity, and (iii) we propose a new criterion for model selection using the Reduction in Residual (RR) to choose an optimal one among the pool of candidates. The new GPSP considers group projected subspaces which is more robust than existing methods in distinguishing correlated group features. We test GP-IDENT on a variety of PDEs and PDE systems, and compare it with the state-of-the-art parametric PDE identification algorithms under different settings to illustrate its outstanding performance. Our experiments show that GP-IDENT is effective in identifying the correct terms from a large dictionary and the model selection scheme is robust to noise.
△ Less
Submitted 11 April, 2023;
originally announced April 2023.
-
Unified $hp$-HDG Frameworks for Friedrichs' PDE systems
Authors:
Jau-Uei Chen,
Shinhoo Kang,
Tan Bui-Thanh,
John N. Shadid
Abstract:
This work proposes a unified $hp$-adaptivity framework for hybridized discontinuous Galerkin (HDG) method for a large class of partial differential equations (PDEs) of Friedrichs' type. In particular, we present unified $hp$-HDG formulations for abstract one-field and two-field structures and prove their well-posedness. In order to handle non-conforming interfaces we simply take advantage of HDG b…
▽ More
This work proposes a unified $hp$-adaptivity framework for hybridized discontinuous Galerkin (HDG) method for a large class of partial differential equations (PDEs) of Friedrichs' type. In particular, we present unified $hp$-HDG formulations for abstract one-field and two-field structures and prove their well-posedness. In order to handle non-conforming interfaces we simply take advantage of HDG built-in mortar structures. With split-type mortars and the approximation space of trace, a numerical flux can be derived via Godunov approach and be naturally employed without any additional treatment. As a consequence, the proposed formulations are parameter-free. We perform several numerical experiments for time-independent and linear PDEs including elliptic, hyperbolic, and mixed-type to verify the proposed unified $hp$-formulations and demonstrate the effectiveness of $hp$-adaptation. Two adaptivity criteria are considered: one is based on a simple and fast error indicator, while the other is rigorous but more expensive using an adjoint-based error estimate. The numerical results show that these two approaches are comparable in terms of convergence rate even for problems with strong gradients, discontinuities, or singularities.
△ Less
Submitted 7 April, 2023;
originally announced April 2023.
-
One stabilization is not enough for closed knotted surfaces
Authors:
Kyle Hayden,
Sungkyung Kang,
Anubhav Mukherjee
Abstract:
In this brief note, we show that there exist smooth 4-manifolds (with nonempty boundary) containing pairs of exotically knotted 2-spheres that remain exotic after one (either external or internal) stabilization. It follows that the ``one is enough'' theorem of Auckly-Kim-Melvin-Ruberman-Schwartz does not hold for closed surfaces whose homology classes are characteristic.
In this brief note, we show that there exist smooth 4-manifolds (with nonempty boundary) containing pairs of exotically knotted 2-spheres that remain exotic after one (either external or internal) stabilization. It follows that the ``one is enough'' theorem of Auckly-Kim-Melvin-Ruberman-Schwartz does not hold for closed surfaces whose homology classes are characteristic.
△ Less
Submitted 3 April, 2023;
originally announced April 2023.
-
Verification and Synthesis of Robust Control Barrier Functions: Multilevel Polynomial Optimization and Semidefinite Relaxation
Authors:
Shucheng Kang,
Yuxiao Chen,
Heng Yang,
Marco Pavone
Abstract:
We study the problem of verification and synthesis of robust control barrier functions (CBF) for control-affine polynomial systems with bounded additive uncertainty and convex polynomial constraints on the control. We first formulate robust CBF verification and synthesis as multilevel polynomial optimization problems (POP), where verification optimizes -- in three levels -- the uncertainty, contro…
▽ More
We study the problem of verification and synthesis of robust control barrier functions (CBF) for control-affine polynomial systems with bounded additive uncertainty and convex polynomial constraints on the control. We first formulate robust CBF verification and synthesis as multilevel polynomial optimization problems (POP), where verification optimizes -- in three levels -- the uncertainty, control, and state, while synthesis additionally optimizes the parameter of a chosen parametric CBF candidate. We then show that, by invoking the KKT conditions of the inner optimizations over uncertainty and control, the verification problem can be simplified as a single-level POP and the synthesis problem reduces to a min-max POP. This reduction leads to multilevel semidefinite relaxations. For the verification problem, we apply Lasserre's hierarchy of moment relaxations. For the synthesis problem, we draw connections to existing relaxation techniques for robust min-max POP, which first use sum-of-squares programming to find increasingly tight polynomial lower bounds to the unknown value function of the verification POP, and then call Lasserre's hierarchy again to maximize the lower bounds. Both semidefinite relaxations guarantee asymptotic global convergence to optimality. We provide an in-depth study of our framework on the controlled Van der Pol Oscillator, both with and without additive uncertainty.
△ Less
Submitted 21 July, 2023; v1 submitted 17 March, 2023;
originally announced March 2023.
-
Learning Subgrid-scale Models with Neural Ordinary Differential Equations
Authors:
Shinhoo Kang,
Emil M. Constantinescu
Abstract:
We propose a new approach to learning the subgrid-scale model when simulating partial differential equations (PDEs) solved by the method of lines and their representation in chaotic ordinary differential equations, based on neural ordinary differential equations (NODEs). Solving systems with fine temporal and spatial grid scales is an ongoing computational challenge, and closure models are general…
▽ More
We propose a new approach to learning the subgrid-scale model when simulating partial differential equations (PDEs) solved by the method of lines and their representation in chaotic ordinary differential equations, based on neural ordinary differential equations (NODEs). Solving systems with fine temporal and spatial grid scales is an ongoing computational challenge, and closure models are generally difficult to tune. Machine learning approaches have increased the accuracy and efficiency of computational fluid dynamics solvers. In this approach neural networks are used to learn the coarse- to fine-grid map, which can be viewed as subgrid-scale parameterization. We propose a strategy that uses the NODE and partial knowledge to learn the source dynamics at a continuous level. Our method inherits the advantages of NODEs and can be used to parameterize subgrid scales, approximate coupling operators, and improve the efficiency of low-order solvers. Numerical results with the two-scale Lorenz 96 ODE, the convection-diffusion PDE, and the viscous Burgers' PDE are used to illustrate this approach.
△ Less
Submitted 12 April, 2023; v1 submitted 19 December, 2022;
originally announced December 2022.
-
Phase transitions for frame potentials]{Phase transitions for the minimizers of the $p^{th}$ frame potentials in $\mathbb{R}^2$
Authors:
Radel Ben Av,
Xuemei Chen,
Assaf Goldberger,
Shujie Kang,
Kasso A. Okoudjou
Abstract:
Given $N$ points $X=\{x_k\}_{k=1}^N$ on the unit circle in $\mathbb{R}^2$ and a number $0\leq p \leq \infty$ we investigate the minimizers of the functional $\sum_{k, \ell =1}^N |\langle x_k, x_\ell\rangle|^p$. While it is known that each of these minimizers is a spanning set for $\mathbb{R}^2$, less is known about their number as a function of $p$ and $N$ especially for relatively small $p$. In t…
▽ More
Given $N$ points $X=\{x_k\}_{k=1}^N$ on the unit circle in $\mathbb{R}^2$ and a number $0\leq p \leq \infty$ we investigate the minimizers of the functional $\sum_{k, \ell =1}^N |\langle x_k, x_\ell\rangle|^p$. While it is known that each of these minimizers is a spanning set for $\mathbb{R}^2$, less is known about their number as a function of $p$ and $N$ especially for relatively small $p$. In this paper we show that there is unique minimum for this functional for all $p\leq \log 3/\log 2$ and all odd $N\geq 3$. In addition, we present some numerical results suggesting the emergence of a phase transition phenomenon for these minimizers. More specifically, for $N\geq 3$ odd, there exists a sequence of number of points $\log 3/\log 2=p_1< p_2< \cdots < p_N\leq 2$ so that a unique (up to some isometries) minimizer exists on each sub-intervals $(p_k, p_{k+1})$. %In addition we conjecture that $\lim_{k\to \infty}p_{2k+1}=2$.
△ Less
Submitted 8 December, 2022;
originally announced December 2022.
-
Young wall construction of level-1 highest weight crystals over $U_q(D_4^{(3)})$ and $U_q(G_2^{(1)})$
Authors:
Zhaobing Fan,
Shaolong Han,
Seok-Jin Kang,
Yong-Su Shin
Abstract:
With the help of path realization and affine energy function, we give a Young wall construction of level-1 highest weight crystals $B(λ)$ over $U_{q}(G_{2}^{(1)})$ and $U_{q}(D_{4}^{(3)})$. Our construction is based on four different shapes of colored blocks, $\mathbf O$-block, $\mathbf I$-block, $\mathbf L$-block and $\mathbf{LL}$-block, obtained by cutting the unit cube in three different ways.
With the help of path realization and affine energy function, we give a Young wall construction of level-1 highest weight crystals $B(λ)$ over $U_{q}(G_{2}^{(1)})$ and $U_{q}(D_{4}^{(3)})$. Our construction is based on four different shapes of colored blocks, $\mathbf O$-block, $\mathbf I$-block, $\mathbf L$-block and $\mathbf{LL}$-block, obtained by cutting the unit cube in three different ways.
△ Less
Submitted 25 February, 2023; v1 submitted 13 November, 2022;
originally announced November 2022.
-
WeakIdent: Weak formulation for Identifying Differential Equations using Narrow-fit and Trimming
Authors:
Mengyi Tang,
Wenjing Liao,
Rachel Kuske,
Sung Ha Kang
Abstract:
Data-driven identification of differential equations is an interesting but challenging problem, especially when the given data are corrupted by noise. When the governing differential equation is a linear combination of various differential terms, the identification problem can be formulated as solving a linear system, with the feature matrix consisting of linear and nonlinear terms multiplied by a…
▽ More
Data-driven identification of differential equations is an interesting but challenging problem, especially when the given data are corrupted by noise. When the governing differential equation is a linear combination of various differential terms, the identification problem can be formulated as solving a linear system, with the feature matrix consisting of linear and nonlinear terms multiplied by a coefficient vector. This product is equal to the time derivative term, and thus generates dynamical behaviors. The goal is to identify the correct terms that form the equation to capture the dynamics of the given data. We propose a general and robust framework to recover differential equations using a weak formulation, for both ordinary and partial differential equations (ODEs and PDEs). The weak formulation facilitates an efficient and robust way to handle noise. For a robust recovery against noise and the choice of hyper-parameters, we introduce two new mechanisms, narrow-fit and trimming, for the coefficient support and value recovery, respectively. For each sparsity level, Subspace Pursuit is utilized to find an initial set of support from the large dictionary. Then, we focus on highly dynamic regions (rows of the feature matrix), and error normalize the feature matrix in the narrow-fit step. The support is further updated via trimming of the terms that contribute the least. Finally, the support set of features with the smallest Cross-Validation error is chosen as the result. A comprehensive set of numerical experiments are presented for both systems of ODEs and PDEs with various noise levels. The proposed method gives a robust recovery of the coefficients, and a significant denoising effect which can handle up to $100\%$ noise-to-signal ratio for some equations. We compare the proposed method with several state-of-the-art algorithms for the recovery of differential equations.
△ Less
Submitted 6 November, 2022;
originally announced November 2022.
-
Crystal bases and canonical bases for quantum Borcherds-Bozec algebras
Authors:
Zhaobing Fan,
Shaolong Han,
Seok-Jin Kang,
Young Rock Kim
Abstract:
Let $U_{q}^{-}(\mathfrak g)$ be the negative half of a quantum Borcherds-Bozec algebra $U_{q}(\mathfrak g)$ and $V(λ)$ be the irreducible highest weight module with $λ\in P^{+}$. In this paper, we investigate the structures, properties and their close connections between crystal bases and canonical bases of $U_{q}^{-}(\mathfrak g)$ and $V(λ)$. We first re-construct crystal basis theory with modifi…
▽ More
Let $U_{q}^{-}(\mathfrak g)$ be the negative half of a quantum Borcherds-Bozec algebra $U_{q}(\mathfrak g)$ and $V(λ)$ be the irreducible highest weight module with $λ\in P^{+}$. In this paper, we investigate the structures, properties and their close connections between crystal bases and canonical bases of $U_{q}^{-}(\mathfrak g)$ and $V(λ)$. We first re-construct crystal basis theory with modified Kashiwara operators. While going through Kashiwara's grand-loop argument, we prove several important lemmas, which play crucial roles in the later developments of the paper. Next, based on the theory of canonical bases on quantum Bocherds-Bozec algebras, we introduce the notion of primitive canonical bases and prove that primitive canonical bases coincide with lower global bases.
△ Less
Submitted 31 March, 2024; v1 submitted 5 November, 2022;
originally announced November 2022.
-
One stabilization is not enough for contractible 4-manifolds
Authors:
Sungkyung Kang
Abstract:
We construct an example of a cork that remains exotic after taking a connected sum with $S^2 \times S^2$. Combined with a work of Akbulut-Ruberman, this implies the existence of an exotic pair of contractible 4-manifolds which remains absolutely exotic after taking a connected sum with $S^2 \times S^2$.
We construct an example of a cork that remains exotic after taking a connected sum with $S^2 \times S^2$. Combined with a work of Akbulut-Ruberman, this implies the existence of an exotic pair of contractible 4-manifolds which remains absolutely exotic after taking a connected sum with $S^2 \times S^2$.
△ Less
Submitted 19 October, 2022; v1 submitted 14 October, 2022;
originally announced October 2022.
-
Hecke equivariance of generalized Borcherds products of type $O(2,1)$
Authors:
Daeyeol Jeon,
Soon-Yi Kang,
Chang Heon Kim
Abstract:
Recently, a weak converse theorem for Borcherds' lifting operator of type $O(2,1)$ for $\G_0(N)$ is proved and the logarithmic derivative of a modular form for $\G_0(N)$ is explicitly described in terms of the values of Niebur-Poincaré series at its divisors in the complex upper half-plane. In this paper, we prove that the generalized Borcherds' lifting operator of type $O(2,1)$ is Hecke equivaria…
▽ More
Recently, a weak converse theorem for Borcherds' lifting operator of type $O(2,1)$ for $\G_0(N)$ is proved and the logarithmic derivative of a modular form for $\G_0(N)$ is explicitly described in terms of the values of Niebur-Poincaré series at its divisors in the complex upper half-plane. In this paper, we prove that the generalized Borcherds' lifting operator of type $O(2,1)$ is Hecke equivariant under the extension of Guerzhoy's multiplicative Hecke operator on the integral weight meromorphic modular forms and the Hecke operator on half-integral weight vector-valued harmonic weak Maass forms. Additionally, we show that the logarithmic differential operator is also Hecke equivariant under the multiplicative Hecke operator and the Hecke operator on integral weight meromorphic modular forms. As applications of Hecke equivariance of the two operators, we obtain relations for twisted traces of singular moduli modulo prime powers and congruences for twisted class numbers modulo primes, including those associated to genus $1$ modular curves.
△ Less
Submitted 7 April, 2023; v1 submitted 8 August, 2022;
originally announced August 2022.
-
The $(2,1)$-cable of the figure-eight knot is not smoothly slice
Authors:
Irving Dai,
Sungkyung Kang,
Abhishek Mallick,
JungHwan Park,
Matthew Stoffregen
Abstract:
We prove that the $(2,1)$-cable of the figure-eight knot is not smoothly slice by showing that its branched double cover bounds no equivariant homology ball.
We prove that the $(2,1)$-cable of the figure-eight knot is not smoothly slice by showing that its branched double cover bounds no equivariant homology ball.
△ Less
Submitted 28 July, 2022;
originally announced July 2022.
-
Torsion in the knot concordance group and cabling
Authors:
Sungkyung Kang,
JungHwan Park
Abstract:
We define a nontrivial mod 2 valued additive concordance invariant defined on the torsion subgroup of the knot concordance group using involutive knot Floer package. For knots not contained in its kernel, we prove that their iterated $(\text{odd},1)$-cables have infinite order in the concordance group and, among them, infinitely many are linearly independent. Furthermore, by taking $(2,1)$-cables…
▽ More
We define a nontrivial mod 2 valued additive concordance invariant defined on the torsion subgroup of the knot concordance group using involutive knot Floer package. For knots not contained in its kernel, we prove that their iterated $(\text{odd},1)$-cables have infinite order in the concordance group and, among them, infinitely many are linearly independent. Furthermore, by taking $(2,1)$-cables of the aforementioned knots, we present an infinite family of knots which are strongly rationally slice but not slice.
△ Less
Submitted 24 July, 2022;
originally announced July 2022.
-
Numerical Identification of Nonlocal Potential in Aggregation
Authors:
Yuchen He,
Sung Ha Kang,
Wenjing Liao,
Hao Liu,
Yingjie Liu
Abstract:
Aggregation equations are broadly used to model population dynamics with nonlocal interactions, characterized by a potential in the equation. This paper considers the inverse problem of identifying the potential from a single noisy spatial-temporal process. The identification is challenging in the presence of noise due to the instability of numerical differentiation. We propose a robust model-base…
▽ More
Aggregation equations are broadly used to model population dynamics with nonlocal interactions, characterized by a potential in the equation. This paper considers the inverse problem of identifying the potential from a single noisy spatial-temporal process. The identification is challenging in the presence of noise due to the instability of numerical differentiation. We propose a robust model-based technique to identify the potential by minimizing a regularized data fidelity term, and regularization is taken as the total variation and the squared Laplacian. A split Bregman method is used to solve the regularized optimization problem. Our method is robust to noise by utilizing a Successively Denoised Differentiation technique. We consider additional constraints such as compact support and symmetry constraints to enhance the performance further. We also apply this method to identify time-varying potentials and identify the interaction kernel in an agent-based system. Various numerical examples in one and two dimensions are included to verify the effectiveness and robustness of the proposed method.
△ Less
Submitted 7 July, 2022;
originally announced July 2022.
-
Inverse problem for a planar conductivity inclusion
Authors:
Doosung Choi,
Johan Helsing,
Sangwoo Kang,
Mikyoung Lim
Abstract:
This paper concerns the inverse problem of determining a planar conductivity inclusion. Our aim is to analytically recover from the generalized polarization tensors (GPTs), which can be obtained from exterior measurements, a homogeneous inclusion with arbitrary constant conductivity. The primary outcome of recovering a homogeneous inclusion is an inversion formula in terms of the GPTs for conforma…
▽ More
This paper concerns the inverse problem of determining a planar conductivity inclusion. Our aim is to analytically recover from the generalized polarization tensors (GPTs), which can be obtained from exterior measurements, a homogeneous inclusion with arbitrary constant conductivity. The primary outcome of recovering a homogeneous inclusion is an inversion formula in terms of the GPTs for conformal mapping coefficients associated with the inclusion. To prove the formula, we establish matrix factorizations for the GPTs.
△ Less
Submitted 19 January, 2023; v1 submitted 11 June, 2022;
originally announced June 2022.
-
Yang-Mills connections and sigma-models on quantum Heisenberg manifolds
Authors:
Stine Marie Berge,
Sooran Kang,
Franz Luef
Abstract:
We construct a spectral triple on a quantum Heisenberg manifold, which generalizes the results of Chakraborty and Shinha, and associate to it an energy functional on the set of projections, following the approach of Mathai-Rosenberg to non-linear sigma models. The spectral triples that we construct extend the We derive a lower bound for this energy functional that is linked on the topological char…
▽ More
We construct a spectral triple on a quantum Heisenberg manifold, which generalizes the results of Chakraborty and Shinha, and associate to it an energy functional on the set of projections, following the approach of Mathai-Rosenberg to non-linear sigma models. The spectral triples that we construct extend the We derive a lower bound for this energy functional that is linked on the topological charge of the projection which depends on the curvature of a compatible connection. A detailed study of this lower bound is given for the Kang projection in quantum Heisenberg manifolds. These results display an intriguing interplay between non-linear sigma models and Yang-Mills theory on quantum Heisenberg manifolds, unlike in the well-studied case of noncommutative tori.
△ Less
Submitted 6 April, 2022; v1 submitted 27 March, 2022;
originally announced March 2022.
-
A Lifted $\ell_1 $ Framework for Sparse Recovery
Authors:
Yaghoub Rahimi,
Sung Ha Kang,
Yifei Lou
Abstract:
Motivated by re-weighted $\ell_1$ approaches for sparse recovery, we propose a lifted $\ell_1$ (LL1) regularization which is a generalized form of several popular regularizations in the literature. By exploring such connections, we discover there are two types of lifting functions which can guarantee that the proposed approach is equivalent to the $\ell_0$ minimization. Computationally, we design…
▽ More
Motivated by re-weighted $\ell_1$ approaches for sparse recovery, we propose a lifted $\ell_1$ (LL1) regularization which is a generalized form of several popular regularizations in the literature. By exploring such connections, we discover there are two types of lifting functions which can guarantee that the proposed approach is equivalent to the $\ell_0$ minimization. Computationally, we design an efficient algorithm via the alternating direction method of multiplier (ADMM) and establish the convergence for an unconstrained formulation. Experimental results are presented to demonstrate how this generalization improves sparse recovery over the state-of-the-art.
△ Less
Submitted 12 May, 2022; v1 submitted 9 March, 2022;
originally announced March 2022.
-
Involutive knot Floer homology and bordered modules
Authors:
Sungkyung Kang
Abstract:
We prove that, up to local equivalences, a suitable truncation of the involutive knot Floer homology of a knot in $S^3$ and the involutive bordered Heegaard Floer theory of its complement determine each other. In particular, given two knots $K_1$ and $K_2$, we prove that the $\mathbb{F}_2[U,V]/(UV)$-coefficient involutive knot Floer homology of $K_1 \sharp -K_2$ is $ι_K$-locally trivial if…
▽ More
We prove that, up to local equivalences, a suitable truncation of the involutive knot Floer homology of a knot in $S^3$ and the involutive bordered Heegaard Floer theory of its complement determine each other. In particular, given two knots $K_1$ and $K_2$, we prove that the $\mathbb{F}_2[U,V]/(UV)$-coefficient involutive knot Floer homology of $K_1 \sharp -K_2$ is $ι_K$-locally trivial if $\widehat{CFD}(S^3 \backslash K_1)$ and $\widehat{CFD}(S^2 \backslash K_2)$ satisfy a certain condition which can be seen as the bordered counterpart of $ι_K$-local equivalence. We further establish an explicit algebraic formula that computes the hat-flavored truncation of the involutive knot Floer homology of a knot from the involutive bordered Floer homology of its complement. It follows that there exists an algebraic satellite operator defined on the local equivalence group of knot Floer chain complexes, which can be computed explicitly up to a suitable truncation.
△ Less
Submitted 12 April, 2022; v1 submitted 25 February, 2022;
originally announced February 2022.
-
Multirate Partitioned Runge-Kutta Methods for Coupled Navier-Stokes Equations
Authors:
Shinhoo Kang,
Alp Dener,
Aidan Hamilton,
Hong Zhang,
Emil M. Constantinescu,
Robert L. Jacob
Abstract:
Earth system models are complex integrated models of atmosphere, ocean, sea ice, and land surface. Coupling the components can be a significant challenge due to the difference in physics, temporal, and spatial scales. This study explores new coupling strategies for the fluid-fluid interaction problem based on multirate partitioned Runge-Kutta methods. We consider compressible Navier-Stokes equatio…
▽ More
Earth system models are complex integrated models of atmosphere, ocean, sea ice, and land surface. Coupling the components can be a significant challenge due to the difference in physics, temporal, and spatial scales. This study explores new coupling strategies for the fluid-fluid interaction problem based on multirate partitioned Runge-Kutta methods. We consider compressible Navier-Stokes equations with gravity coupled through a rigid-lid interface. Our large-scale numerical experiments reveal that multirate partitioned Runge-Kutta coupling schemes (1) can conserve total mass; (2) have second-order accuracy in time; and (3) provide favorable strong- and weak-scaling performance on modern computing architectures. We also show that the speedup factors of multirate partitioned Runge-Kutta methods match theoretical expectations over their base (single-rate) method.
△ Less
Submitted 10 April, 2023; v1 submitted 23 February, 2022;
originally announced February 2022.
-
Monostatic sampling methods in limited-aperture configuration
Authors:
Sangwoo Kang,
Mikyoung Lim
Abstract:
We present monostatic sampling methods for limited-aperture scattering problems in two dimensions. The direct sampling method (DSM) is well known to provide a robust, stable, and fast numerical scheme for imaging inhomogeneities from multistatic measurements even with only one or two incident fields. However, in practical applications, monostatic measurements in limited-aperture configuration are…
▽ More
We present monostatic sampling methods for limited-aperture scattering problems in two dimensions. The direct sampling method (DSM) is well known to provide a robust, stable, and fast numerical scheme for imaging inhomogeneities from multistatic measurements even with only one or two incident fields. However, in practical applications, monostatic measurements in limited-aperture configuration are frequently encountered. A monostatic sampling method (MSM) was studied in full-aperture configuration in recent literature. In this paper, we develop MSM in limited-aperture configuration and derive an asymptotic formula of the corresponding indicator function. Based on the asymptotic formula, we then analyze the imaging performance of the proposed method depending on the range of measurement directions and the geometric, material properties of inhomogeneities. Furthermore, we propose a modified numerical scheme with multi-frequency measurements that improve imaging performance, especially for small anomalies. Numerical simulations are presented to validate the analytical results.
△ Less
Submitted 13 April, 2022; v1 submitted 28 October, 2021;
originally announced October 2021.
-
Stabilization and satellite construction of doubly slice links
Authors:
Hongtaek Jung,
Sungkyung Kang,
Seungwon Kim
Abstract:
A 2-component oriented link in $S^3$ is called weakly doubly slice if it is a cross-section of an unknotted sphere in $S^4$, and strongly doubly slice if it is a cross-section of a 2-component trivial spherical link in $S^4$. We give the first example of 2-component boundary links which are weakly doubly slice but not strongly doubly slice. We also introduce a new invariant $g_{st}$ of homotopical…
▽ More
A 2-component oriented link in $S^3$ is called weakly doubly slice if it is a cross-section of an unknotted sphere in $S^4$, and strongly doubly slice if it is a cross-section of a 2-component trivial spherical link in $S^4$. We give the first example of 2-component boundary links which are weakly doubly slice but not strongly doubly slice. We also introduce a new invariant $g_{st}$ of homotopically trivial links that measures the failure of a link from being strongly doubly slice and that bounds the doubly slice genus $g_{ds}$ from below. Our examples have arbitrarily large doubly slice genus but satisfy $g_{st}=1$. We also prove that the Conway-Orson signature lower bound on $g_{ds}$ is actually a lower bound on $g_{st}$.
△ Less
Submitted 26 October, 2021;
originally announced October 2021.
-
Counting Objects by Diffused Index: geometry-free and training-free approach
Authors:
Mengyi Tang,
Maryam Yashtini,
Sung Ha Kang
Abstract:
Counting objects is a fundamental but challenging problem. In this paper, we propose diffusion-based, geometry-free, and learning-free methodologies to count the number of objects in images. The main idea is to represent each object by a unique index value regardless of its intensity or size, and to simply count the number of index values. First, we place different vectors, refer to as seed vector…
▽ More
Counting objects is a fundamental but challenging problem. In this paper, we propose diffusion-based, geometry-free, and learning-free methodologies to count the number of objects in images. The main idea is to represent each object by a unique index value regardless of its intensity or size, and to simply count the number of index values. First, we place different vectors, refer to as seed vectors, uniformly throughout the mask image. The mask image has boundary information of the objects to be counted. Secondly, the seeds are diffused using an edge-weighted harmonic variational optimization model within each object. We propose an efficient algorithm based on an operator splitting approach and alternating direction minimization method, and theoretical analysis of this algorithm is given. An optimal solution of the model is obtained when the distributed seeds are completely diffused such that there is a unique intensity within each object, which we refer to as an index. For computational efficiency, we stop the diffusion process before a full convergence, and propose to cluster these diffused index values. We refer to this approach as Counting Objects by Diffused Index (CODI). We explore scalar and multi-dimensional seed vectors. For Scalar seeds, we use Gaussian fitting in histogram to count, while for vector seeds, we exploit a high-dimensional clustering method for the final step of counting via clustering. The proposed method is flexible even if the boundary of the object is not clear nor fully enclosed. We present counting results in various applications such as biological cells, agriculture, concert crowd, and transportation. Some comparisons with existing methods are presented.
△ Less
Submitted 15 October, 2021;
originally announced October 2021.
-
On the nonorientable four-ball genus of torus knots
Authors:
Fraser Binns,
Sungkyung Kang,
Jonathan Simone,
Paula Truöl
Abstract:
The nonorientable four-ball genus of a knot $K$ in $S^3$ is the minimal first Betti number of nonorientable surfaces in $B^4$ bounded by $K$. By amalgamating ideas from involutive knot Floer homology and unoriented knot Floer homology, we give a new lower bound on the smooth nonorientable four-ball genus $γ_4$ of any knot. This bound is sharp for several families of torus knots, including…
▽ More
The nonorientable four-ball genus of a knot $K$ in $S^3$ is the minimal first Betti number of nonorientable surfaces in $B^4$ bounded by $K$. By amalgamating ideas from involutive knot Floer homology and unoriented knot Floer homology, we give a new lower bound on the smooth nonorientable four-ball genus $γ_4$ of any knot. This bound is sharp for several families of torus knots, including $T_{4n,(2n\pm 1)^2}$ for even $n\ge 2$, a family Longo showed were counterexamples to Batson's conjecture. We also prove that, whenever $p$ is an even positive integer and $\frac{p}{2}$ is not a perfect square, the torus knot $T_{p,q}$ does not bound a locally flat Möbius band for almost all integers $q$ relatively prime to $p$.
△ Less
Submitted 24 June, 2024; v1 submitted 19 September, 2021;
originally announced September 2021.
-
Entropy-Preserving and Entropy-Stable Relaxation IMEX and Multirate Time-Stepping Methods
Authors:
Shinhoo Kang,
Emil M. Constantinescu
Abstract:
We propose entropy-preserving and entropy-stable partitioned Runge--Kutta (RK) methods. In particular, we extend the explicit relaxation Runge--Kutta methods to IMEX--RK methods and a class of explicit second-order multirate methods for stiff problems arising from scale-separable or grid-induced stiffness in a system. The proposed approaches not only mitigate system stiffness but also fully suppor…
▽ More
We propose entropy-preserving and entropy-stable partitioned Runge--Kutta (RK) methods. In particular, we extend the explicit relaxation Runge--Kutta methods to IMEX--RK methods and a class of explicit second-order multirate methods for stiff problems arising from scale-separable or grid-induced stiffness in a system. The proposed approaches not only mitigate system stiffness but also fully support entropy-preserving and entropy-stability properties at a discrete level. The key idea of the relaxation approach is to adjust the step completion with a relaxation parameter so that the time-adjusted solution satisfies the entropy condition at a discrete level. The relaxation parameter is computed by solving a scalar nonlinear equation at each timestep in general; however, as for a quadratic entropy function, we theoretically derive the explicit form of the relaxation parameter and numerically confirm that the relaxation parameter works the Burgers equation. Several numerical results for ordinary differential equations and the Burgers equation are presented to demonstrate the entropy-conserving/stable behavior of these methods. We also compare the relaxation approach and the incremental direction technique for the Burgers equation with and without a limiter in the presence of shocks.
△ Less
Submitted 19 July, 2022; v1 submitted 19 August, 2021;
originally announced August 2021.
-
Global bases for quantum Borcherds-Bozec algebras
Authors:
Zhaobing Fan,
Seok-Jin Kang,
Young Rock Kim,
Bolun Tong
Abstract:
We provide a construction of global bases for quantum Borcherds-Bozec algebras and their integrable highest weight representations.
We provide a construction of global bases for quantum Borcherds-Bozec algebras and their integrable highest weight representations.
△ Less
Submitted 10 August, 2021;
originally announced August 2021.