No abstract uploaded!

We construct a smooth algebraic stack of tuples consisting of genus two nodal curves, line bundles, and twisted fields. It leads to a desingularization of the moduli of genus two stable maps to projective spaces. The construction is based on systematical application of the theory of stacks with twisted fields (STF), which has its prototype appeared in arXiv:1906.10527 and arXiv:1201.2427 and is fully developed in this article. The results of this article are the second step of a series of works toward the resolutions of the moduli of stable maps of higher genera.

We construct a smooth Artin stack parameterizing the stable weighted curves of genus one with twisted fields and prove that it is isomorphic to the blowup stack of the moduli of genus one weighted curves studied by Hu and Li. This leads to a blowup-free construction of Vakil-Zinger's desingularization of the moduli of genus one stable maps to projective spaces. This construction provides the cornerstone of the theory of stacks with twisted fields, which is thoroughly studied in arXiv:2005.03384 and leads to a blowup-free resolution of the stable map moduli of genus two.

Optimal mass transportation has been widely applied in various fields, such as data compression, generative adversarial networks, and image processing. In this paper, we adopt the projected gradient method, combined with the homotopy technique, to nd a minimal volume-measure-preserving solution for a 3-manifold optimal mass transportation problem. The proposed projected gradient method is shown to be sublinearly convergent at a rate of O(1/k). Several numerical experiments indicate that our algorithms are able to signicantly reduce transportation costs. Some applications of the optimal mass transportation maps - to deformations and canonical normalizations between brains and solid balls - are demonstrated to show the robustness of our proposed algorithms.

In the present paper, we prove the existence of solutions $(\lambda_1,\lambda_2,u,v)\in\R^2\times H^1(\R^3,\R^2)$ to systems of coupled Schr\"odinger equations $$ \begin{cases} -\Delta u+\lambda_1u=\mu_1 u^3+\beta uv^2\quad &\hbox{in}\;\R^3\\ -\Delta v+\lambda_2v=\mu_2 v^3+\beta u^2v\quad&\hbox{in}\;\R^3\\ u,v>0&\hbox{in}\;\R^3 \end{cases} $$ satisfying the normalization constraint $ \displaystyle\int_{\R^3}u^2=a^2\quad\hbox{and}\;\int_{\R^3}v^2=b^2, $ which appear in binary mixtures of Bose-Einstein condensates or in nonlinear optics. The parameters $\mu_1,\mu_2,\beta>0$ are prescribed as are the masses $a,b>0$. The system has been considered mostly in the case of fixed frequencies $\lambda_1,\lambda_2$. When the masses are prescribed, the standard approach to this problem is variational with $\lambda_1,\lambda_2$ appearing as Lagrange multipliers. Here we present a new approach based on the fixed point index in cones, bifurcation theory, and the continuation method. We obtain the existence of normalized solutions for any given $a,b>0$ for $\beta$ in a large range. We also have a result about the nonexistence of positive solutions which shows that our existence theorem is almost optimal. Especially, if $\mu_1=\mu_2$ we prove that normalized solutions exist for all $\beta>0$ and all $a,b>0$.

No abstract uploaded!

No abstract uploaded!

A two species interacting system motivated by the density functional theory for triblock copolymers contains long range interaction that affects the two species differently. In a two species periodic assembly of discs, the two species appear alternately on a lattice. The minimal two species periodic assembly is the one with the least energy per lattice cell area. There is a parameter $b$ in $[0,1]$ and the type of the lattice associated with the minimal assembly varies depending on $b$. There are several threshold defined by a number $B=0.1867...$ If $b \in [0, B)$, the minimal assembly is associated with a rectangular lattice whose ratio of the longer side and the shorter side is in $[\sqrt{3}, 1)$; if $b \in [B, 1-B]$, the minimal assembly is associated with a square lattice; if $b \in (1-B, 1]$, the minimal assembly is associated with a rhombic lattice with an acute angle in $[\frac{\pi}{3}, \frac{\pi}{2})$. Only when $b=1$, this rhombic lattice is the hexagonal lattice. None of the other values of $b$ yields the hexagonal lattice, a sharp contrast to the situation for one species interacting systems, where the hexagonal lattice is ubiquitously observed.

We show the intersection of a compact almost complex subvariety of dimension 4 and a compact almost complex submanifold of codimension 2 is a J-holomorphic curve. This is a generalization of positivity of intersections for J-holomorphic curves in almost complex 4-manifolds to higher dimensions. As an application, we discuss pseudoholomorphic sections of a complex line bundle. We introduce a method to produce J-holomorphic curves using the differential geometry of almost Hermitian manifolds. When our main result is applied to pseudoholomorphic maps, we prove the singularity subset of a pseudoholomorphic map between almost complex 4-manifolds is J-holomorphic. Building on this, we show degree one pseudoholomorphic maps between almost complex 4-manifolds are actually birational morphisms in pseudoholomorphic category.

No abstract uploaded!

We present a new Matched Interface and Boundary (MIB) regularization method for treating charge singularity in solvated biomolecules whose electrostatics are described by the Poisson–Boltzmann (PB) equation. In a regularization method, by decomposing the potential function into two or three components, the singular component can be analytically represented by the Green’s function, while other components possess a higher regularity. Our new regularization combines the efficiency of two-component schemes with the accuracy of the three-component schemes. Based on this regularization, a new MIB finite difference algorithm is developed for solving both linear and nonlinear PB equations, where the nonlinearity is handled by using the inexact-Newton’s method. Compared with the existing MIB PB solver based on a three-component regularization, the present algorithm is simpler to implement by circumventing the work to solve a boundary value Poisson equation inside the molecular interface and to compute related interface jump conditions numerically. Moreover, the new MIB algorithm becomes computationally less expensive, while maintains the same second order accuracy. This is numerically verified by calculating the electrostatic potential and solvation energy on the Kirkwood sphere on which the analytical solutions are available and on a series of proteins with various sizes.

In this paper, we introduce multiscale persistent functions for biomolecular structure characterization. The essential idea is to combine our multiscale rigidity functions with persistent homology analysis, so as to construct a series of multiscale persistent functions, particularly multiscale persistent entropies, for structure characterization. To clarify the fundamental idea of our method, the multiscale persistent entropy model is discussed in great detail. Mathematically, unlike the previous persistent entropy or topological entropy, a special resolution parameter is incorporated into our model. Various scales can be achieved by tuning its value. Physically, our multiscale persistent entropy can be used in conformation entropy evaluation. More specifically, it is found that our method incorporates in it a natural classification scheme. This is achieved through a density filtration of a multiscale rigidity function built from bond and/or dihedral angle distributions. To further validate our model, a systematical comparison with the traditional entropy evaluation model is done. It is found that our model is able to preserve the intrinsic topological features of biomolecular data much better than traditional approaches, particularly for resolutions in the mediate range. Moreover, our method can be successfully used in protein classification. For a test database with around nine hundred proteins, a clear separation between all-alpha and all-beta proteins can be achieved, using only the dihedral and pseudo-bond angle information. Finally, a special protein structure index (PSI) is proposed, for the first time, to describe the "regularity" of protein structures. Essentially, PSI can be used to describe the "regularity" information in any systems.

Principal component analysis (PCA) is arguably the most widely used dimension-reduction method for vector-type data. When applied to a sample of images, PCA requires vectorization of the image data, which in turn entails solving an eigenvalue problem for the sample covariance matrix. We propose herein a two-stage dimension reduction (2SDR) method for image reconstruction from high-dimensional noisy image data. The first stage treats the image as a matrix, which is a tensor of order 2, and uses multilinear principal component analysis (MPCA) for matrix rank reduction and image denoising. The second stage vectorizes the reduced-rank matrix and achieves further dimension and noise reduction. Simulation studies demonstrate excellent performance of 2SDR, for which we also develop an asymptotic theory that establishes consistency of its rank selection. Applications to cryo-EM (cryogenic electronic microscopy), which has revolutionized structural biology, organic and medical chemistry, cellular and molecular physiology in the past decade, are also provided and illustrated with benchmark cryo-EM datasets. Connections to other contemporaneous developments in image reconstruction and high-dimensional statistical inference are also discussed.

Ligand-receptor binding and unbinding are fundamental biomolecular processes and particularly essential to drug efficacy. Environmental water fluctuations, however, impact the corresponding thermodynamics and kinetics and thereby challenge theoretical descriptions. Here, we devise a holistic, implicit-solvent, multi-method approach to predict the (un)binding kinetics for a generic ligand-pocket model. We use the variational implicit-solvent model (VISM) to calculate the solute-solvent interfacial structures and the corresponding free energies, and combine the VISM with the string method to obtain the minimum energy paths and transition states between the various metastable (“dry” and “wet”) hydration states. The resulting dry-wet transition rates are then used in a spatially-dependent multi-state continuous-time Markov chain Brownian dynamics simulations, and the related Fokker–Planck equation calculations, of the ligand stochastic motion, providing the mean first-passage times for binding and unbinding. We find the hydration transitions to significantly slow down the binding process, in semi-quantitative agreement with existing explicit-water simulations, but significantly accelerate the unbinding process. Moreover, our methods allow the characterization of non-equilibrium hydration states of pocket and ligand during the ligand movement, for which we find substantial memory and hysteresis effects for binding versus unbinding. Our study thus provides a significant step forward towards efficient, physics-based interpretation and predictions of the complex kinetics in realistic ligand-receptor systems.

No abstract uploaded!

Let $(V, 0)$ be an isolated hypersurface singularity defined by the holomorphic function $f: (\mathbb{C}^n, 0)\rightarrow (\mathbb{C}, 0)$. In our previous work, we introduced a series of novel Lie algebras associated to $(V,0)$, i.e., $k$-th Yau algebra $L^k(V), k\geq 0$. It was defined to be the Lie algebra of derivations of the $k$-th moduli algebras $A^k(V)= \mathcal{O}_n/(f, m^k J(f)), k\geq 0$, where $m$ is the maximal ideal of $\mathcal{O}_n$. I.e., $L^k(V):=\text{Der}(A^k(V), A^k(V))$. The dimension of $L^k(V)$ was denoted by $\lambda^k(V)$. The number $\lambda^k(V)$, which was called $k$-th Yau number, is a subtle numerical analytic invariant of $(V, 0)$. Furthermore, we formulated two conjectures for these $k$-th Yau number invariants: a sharp upper estimate conjecture of $\lambda^k(V)$ for weighted homogeneous isolated hypersurface singularities (see Conjecture \ref{conj2}) and an inequality conjecture $\lambda^{(k+1)}(V) > \lambda^k(V), k\geq0$ (see Conjecture \ref{conj1}). In this article, we verify these two conjectures when $k$ is small for large class of singularities.

We investigate the asymptotic distributions of coordinates of regression M-estimates in the moderate p/n regime, where the number of covariates p grows proportionally with the sample size n. Under appropriate regularity conditions, we establish the coordinate-wise asymptotic normality of regression M-estimates assuming a fixed-design matrix. Our proof is based on the second-order Poincare inequality (Chatterjee, 2009) and leave-one-out analysis (El Karoui et al., 2011). Some relevant examples are indicated to show that our regularity conditions are satisfied by a broad class of design matrices. We also show a counterexample, namely the ANOVA-type design, to emphasize that the technical assumptions are not just artifacts of the proof. Finally, the numerical experiments confirm and complement our theoretical results.

We analyze the asymptotic behavior of certain twisted orbital integrals arising from the study of affine Deligne-Lusztig varieties. The main tools include the Base Change Fundamental Lemma and q-analogues of the Kostant partition functions. As an application we prove a conjecture of Miaofen Chen and Xinwen Zhu, relating the set of irreducible components of an affine Deligne-Lusztig variety modulo the action of the σ-centralizer group to the Mirkovic-Vilonen basis of a certain weight space of a representation of the Langlands dual group.

No abstract uploaded!

No abstract uploaded!

We consider the fundamental problem of solving quadratic systems of m equations in n variables. We propose a novel method, which starting with an initial guess computed by means of a spectral method, proceeds by minimizing a nonconvex functional as in the Wirtinger flow approach. There are several key distinguishing features, most notably, a distinct objective functional and novel update rules, which operate in an adaptive fashion and drop terms bearing too much influence on the search direction. These careful selection rules provide a tighter initial guess, better descent directions, and thus enhanced practical performance. On the theoretical side, we prove that for certain unstructured models of quadratic systems, our algorithms return the correct solution in linear time, i.e. in time proportional to reading the data and as soon as the ratio m/n between the number of equations and unknowns exceeds a fixed numerical constant. We extend the theory to deal with noisy systems and prove that our algorithms achieve a statistical accuracy, which is nearly un-improvable. We complement our theoretical study with numerical examples showing that solving random quadratic systems is both computationally and statistically not much harder than solving linear systems of the same size---hence the title of this paper. For instance, we demonstrate empirically that the computational cost of our algorithm is about four times that of solving a least-squares problem of the same size.

Community detection is a fundamental problem in social network analysis. It is desirable to have fast community detection algorithm that applicable to large-scale networks we see today and also have reliable community detection results. This paper proposes Spectral Clustering On Ratios-of-eigenvectors (SCORE) as such an algorithm. The algorithm is fast in computation, is applicable to real large networks, and is competitive in real data performance when applied to some representative data sets in this area. We study the algorithm with the Degree-Corrected Block Model (DCBM). Compare to most works in this area that focus on the much narrower Block Models, our analysis is more difficult and requires new techniques. The novelty of our algorithm lies in that it uses entry-wise ratios between leading eigenvectors of the adjacency to reduce a seemingly hard problem to a problem that is very much tractable.

Let $(X,\alpha)$ be a K\"ahler manifold of dimension $n$, and let $[\omega] \in H^{1,1}(X,\mathbb{R})$. We study the problem of specifying the Lagrangian phase of $\omega$ with respect to $\alpha$, which is described by the nonlinear elliptic equation \[ \sum_{i=1}^{n} \arctan(\la_i)= h(x) \] where $\la_i$ are the eigenvalues of $\omega$ with respect to $\alpha$. When $h(x)$ is a topological constant, this equation corresponds to the deformed Hermitian-Yang-Mills (dHYM) equation, and is related by Mirror Symmetry to the existence of special Lagrangian submanifolds of the mirror. We introduce a notion of subsolution for this equation, and prove a priori $C^{2,\beta}$ estimates when $|h|>(n-2)\frac{\pi}{2}$ and a subsolution exists. Using the method of continuity we show that the dHYM equation admits a smooth solution in the supercritical phase case, whenever a subsolution exists. Finally, we discover some stability-type cohomological obstructions to the existence of solutions to the dHYM equation and we conjecture that when these obstructions vanish the dHYM equation admits a solution. We confirm this conjecture for complex surfaces.

We study the deformed Hermitian-Yang-Mills (dHYM) equation, which is mirror to the special Lagrangian equation, from the variational point of view via an infinite dimensional GIT problem mirror to Thomas' GIT picture for special Lagrangians. This gives rise to infinite dimensional manifold $\mathcal{H}$ mirror to Solomon's space of positive Lagrangians. In the hypercritical phase case we prove the existence of smooth approximate geodesics, and weak geodesics with $C^{1,\alpha}$ regularity. This is accomplished by proving sharp with respect to scale estimates for the Lagrangian phase operator on collapsing manifolds with boundary. We apply these results to the infinite dimensional GIT problem for deformed Hermitian-Yang-Mills. We associate algebraic invariants to certain birational models of $X\times \Delta$, where $\Delta \subset \mathbb{C}$ is a disk. Using the existence of regular weak geodesics we prove that these invariants give rise to obstructions to the existence of solutions to the dHYM equation. Furthermore, we show that these invariants fit into a stability framework closely related to Bridgeland stability. Finally, we use a Fourier-Mukai transform on toric K\"ahler manifolds to describe degenerations of Lagrangian sections of SYZ torus fibrations of Landau-Ginzburg models $(Y,W)$. We speculate on the resulting algebraic invariants, and discuss the implications for relating Bridgeland stability to the existence of special Lagrangian sections of $(Y,W)$.

In this paper, we derive estimates for scalar curvature type equations with more singular right hand side. As an application, we prove Donaldson's conjecture on the equivalence between geodesic stability and existence of cscK when $Aut_0(M,J)\neq0$. Moreover, we show that when $Aut_0(M,J)\neq0$, the properness of $K$-energy with respect to a suitably defined distance implies the existence of cscK.