Multigrid Methods for Hellan-Herrmann-Johnson Mixed Method of Kirchhoff Plate Bending Problems

A V-cycle multigrid method for the Hellan-Herrmann-Johnson (HHJ) discretization of the Kirchhoff plate bending problems is developed in this paper. It is shown that the contraction number of the V-cycle multigrid HHJ mixed method is bounded away from one uniformly with respect to the mesh size. The uniform convergence is achieved for the V-cycle multigrid method with only one smoothing step and without full elliptic regularity. The key is a stable decomposition of the kernel space which is derived from an exact sequence of the HHJ mixed method, and the strengthened Cauchy Schwarz inequality. Some numerical experiments are provided to confirm the proposed V-cycle multigrid method. The exact sequences of the HHJ mixed method and the corresponding commutative diagram is of some interest independent of the current context.Comment: 23 page

Stabilized mixed finite element methods for linear elasticity on simplicial grids in Rn\mathbb{R}^{n}

In this paper, we design two classes of stabilized mixed finite element methods for linear elasticity on simplicial grids. In the first class of elements, we use $\boldsymbol{H}(\mathbf{div}, \Omega; \mathbb{S})$-$P_k$ and $\boldsymbol{L}^2(\Omega; \mathbb{R}^n)$-$P_{k-1}$ to approximate the stress and displacement spaces, respectively, for $1\leq k\leq n$, and employ a stabilization technique in terms of the jump of the discrete displacement over the faces of the triangulation under consideration; in the second class of elements, we use $\boldsymbol{H}_0^1(\Omega; \mathbb{R}^n)$-$P_{k}$ to approximate the displacement space for $1\leq k\leq n$, and adopt the stabilization technique suggested by Brezzi, Fortin, and Marini. We establish the discrete inf-sup conditions, and consequently present the a priori error analysis for them. The main ingredient for the analysis is two special interpolation operators, which can be constructed using a crucial $\boldsymbol{H}(\mathbf{div})$ bubble function space of polynomials on each element. The feature of these methods is the low number of global degrees of freedom in the lowest order case. We present some numerical results to demonstrate the theoretical estimates.Comment: 16 pages, 1 figur

Quantum Game with Restricted Matrix Strategies

We study a quantum game played by two players with restricted multiple strategies. It is found that in this restricted quantum game Nash equilibrium does not always exist when the initial state is entangled. At the same time, we find that when Nash equilibrium exists the pay off function is usually different from that in the classical counterpart except in some special cases. This presents an explicit example where quantum game and classical game may differ. When designing a quantum game with limited strategies, the allowed strategy should be carefully chosen according to the type of initial state.Comment: 5 pages and 3 figure

Residual-Based A Posteriori Error Estimates for Symmetric Conforming Mixed Finite Elements for Linear Elasticity Problems

A posteriori error estimators for the symmetric mixed finite element methods for linear elasticity problems of Dirichlet and mixed boundary conditions are proposed. Stability and efficiency of the estimators are proved. Finally, we provide numerical examples to verify the theoretical results

Two-Loop integrals for CP-even heavy quarkonium production and decays: Elliptic Sectors

By employing the differential equations, we compute analytically the elliptic sectors of two-loop master integrals appearing in the NNLO QCD corrections to CP-even heavy quarkonium exclusive production and decays, which turns out to be the last and toughest part in the relevant calculation. The integrals are found can be expressed as Goncharov polylogarithms and iterative integrals over elliptic functions. The master integrals may be applied to some other NNLO QCD calculations about heavy quarkonium exclusive production, like $\gamma^*\gamma\rightarrow Q\bar{Q}$, $e^+e^-\rightarrow \gamma+ Q\bar{Q}$,~and~$H/Z^0\rightarrow \gamma+ Q\bar{Q}$, heavy quarkonium exclusive decays, and also the CP-even heavy quarkonium inclusive production and decays.Comment: 23 pages, 3 figures, more discussions and references adde

Advantages of the multinucleon transfer reactions based on 238U target for producing neutron-rich isotopes around N = 126

The mechanism of multinucleon transfer (MNT) reactions for producing neutron-rich heavy nuclei around N = 126 is investigated within two different theoretical frameworks: dinuclear system (DNS) model and isospin-dependent quantum molecular dynamics (IQMD) model. The effects of mass asymmetry relaxation, N=Z equilibration, and shell closures on production cross sections of neutron-rich heavy nuclei are investigated. For the first time, the advantages for producing neutron-rich heavy nuclei around N = 126 is found in MNT reactions based on 238U target. We propose the reactions with 238U target for producing unknown neutron-rich heavy nuclei around N = 126 in the future.Comment: 6 pages, 6 figure