• ### Symmetric structure for the endomorphism algebra of projective-injective module in parabolic category(1702.05834)

Sept. 9, 2018 math.RT
We show that for any singular dominant integral weight $\lambda$ of a complex semisimple Lie algebra $\mathfrak{g}$, the endomorphism algebra $B$ of any projective-injective module of the parabolic BGG category $\mathcal{O}_\lambda^{\mathfrak{p}}$ is a symmetric algebra (as conjectured by Khovanov) extending the results of Mazorchuk and Stroppel for the regular dominant integral weight. Moreover, the endomorphism algebra $B$ is equipped with a homogeneous (non-degenerate) symmetrizing form. In the appendix, there is a short proof due to K. Coulembier and V. Mazorchuk showing that the endomorphism algebra $B_\lambda^{\mathfrak{p}}$ of the basic projective-injective module of $\mathcal{O}_\lambda^{\mathfrak{p}}$ is a symmetric algebra.
• ### Modified affine Hecke algebras and quiver Hecke algebras of type $A$(1608.05453)

June 16, 2019 math.RT
We introduce some modified forms for the degenerate and non-degenerate affine Hecke algebras of type $A$. These are certain subalgebras living inside the inverse limit of cyclotomic Hecke algebras. We construct faithful representations and standard bases for these algebras and give some explicit description of their centers. We show that there are algebra isomorphisms between some generalized Ore localizations of these modified affine Hecke algebras and of the quiver Hecke algebras of type $A$. As an application, we show that the center conjecture for the cyclotomic quiver Hecke algebra of type $A$ holds if and only if the center conjecture for the cyclotomic Hecke algebra of type $A$ holds.
• ### Grassmannian, symmetric functions and cyclotomic nilHecke algebras(1804.01467)

April 11, 2018 math.RT
Let $\ell, n$ be positive integers such that $\ell\geq n$. Let $\mathbb{G}_{n,\ell-n}$ be the Grassmannian manifold which consists of the set of $n$-dimensional subspaces of $\mathbb{C}^{\ell}$. There is an $\mathbb{Z}$-graded algebra isomorphism between the cohomology $\rm{H}*(\mathbb{G}_{n,\ell-n},\mathbb{Z})$ of $\mathbb{G}_{n,\ell-n}$ and a natural $\mathbb{Z}$-form $B$ of the $\mathbb{Z}$-graded basic algebra of the type $A$ cyclotomic nilHecke algebra $\mathscr{H}_{\ell,n}^{(0)}=< \psi_1,\cdots,\psi_{n-1},y_1,\cdots,y_n>$. In this paper, we show that the isomorphism can be chosen such that the image of each (geometrically defined) Schubert class $(a_1,\cdots,a_{n})$ coincides with the basis element $b_{\mathbf{\lambda}}$ constructed in an earlier paper by the second author and Xinfeng Liang using purely algebraic method, where $0\leq a_1\leq a_2\leq\cdots\leq a_{n}\leq \ell-n$ with $a_i\in\mathbb{Z}$ for each $i$, $\mathbb{\lambda}$ is the $\ell$-multipartition of $n$ associated to $\(\ell+1-(a_{n}+n), \ell+1-(a_{n-1}+n-1),\cdots,\ell+1-(a_1+1))$. A similar isomorphism between the cohomology $\rm{H}^*(\mathbb{G}_{\ell-n,n},\mathbb{Z})$ of the Grassmannian $\mathbb{G}_{\ell-n,n}$ and the natural $\mathbb{Z}$-form $B$ of the $\mathbb{Z}$-graded basic algebra of $\mathscr{H}_{\ell,n}^{(0)}$ is also obtained. As applications, we obtain a second version of Giambelli formula for Schubert classes and show that each basis element $z_{\mathbf{\lambda}}$ (constructed in the earlier paper mentioned above) of the center $Z$ of $\mathscr{H}_{\ell,n}^{(0)}$ is equal to the evaluation of a Schur (symmetric) polynomial at $y_1,\cdots,y_n$.
• ### The Flash ADC system and PMT waveform reconstruction for the Daya Bay Experiment(1707.03699)

April 9, 2018 physics.ins-det
To better understand the energy response of the Antineutrino Detector (AD), the Daya Bay Reactor Neutrino Experiment installed a full Flash ADC readout system on one AD that allowed for simultaneous data taking with the current readout system. This paper presents the design, data acquisition, and simulation of the Flash ADC system, and focuses on the PMT waveform reconstruction algorithms. For liquid scintillator calorimetry, the most critical requirement to waveform reconstruction is linearity. Several common reconstruction methods were tested but the linearity performance was not satisfactory. A new method based on the deconvolution technique was developed with 1% residual non-linearity, which fulfills the requirement. The performance was validated with both data and Monte Carlo (MC) simulations, and 1% consistency between them has been achieved.
• ### Engineering giant magnetic anisotropy in single-molecule magnets by dimerizing heavy transition-metal atoms(1803.08213)

March 22, 2018 cond-mat.mes-hall
Search for single-molecule magnets with large magnetic anisotropy energy (MAE) is essential for the development of molecular spintronics devices used at room temperature. Through systematic first-principles calculations, we found that an Os-Os or Ir-Ir dimer embedded in the (5,5'-Br2-Salophen) molecule gives rise to large MAE of 41.6 or 51.4 meV which is large enough to hold the spin orientation at room temperature. Analysis of electronic structures reveals that the top Os and Ir atoms are most responsible for the spin moments and large MAEs of the molecules.
• ### Invariant Wide Bandgaps in Honeycomb Monolayer and Single-Walled Nanotubes of IIB-VI Semiconductors(1803.08217)

March 22, 2018 cond-mat.mtrl-sci
Search for low-dimensional materials with unique electronic properties is important for the development of electronic devices in nano scale. Through systematic first-principles calculations, we found that the band gaps of the two-dimensional honeycomb monolayers and one-dimensional single-walled nanotubes of IIB-VI semiconductors (ZnO, CdO, ZnS and CdS) are nearly chirality-independent and weakly diameter-dependent. Based on analysis of the electronic structures, it was found that the conduction band minimum is contributed by the spherically symmetric $s$ orbitals of cations and the valence band maximum is dominated by the in-plane $(d_{xy}-p_y)$ and $(d_{x^2-y^2}-p_x)$ hybridizations. These electronic states are robust against radius curvature, resulting in the invariant feature of the band gaps for the structures changing from honeycomb monolayer to single-walled nanotubes. The band gaps of these materials range from 2.3 eV to 4.7 eV, which is of potential applications in electronic devices and optoelectronic devices. Our studies show that searching for and designing specific electronic structures can facilitate the process of exploring novel nanomaterials for future applications.
• ### Large Perpendicular Magnetocrystalline Anisotropy at Fe/Pb(001) interface(1803.08219)

March 22, 2018 cond-mat.mtrl-sci
Search for ultrathin magnetic film with large perpendicular magnetocrystalline anisotropy (PMA) has been inspired for years by the continuous miniaturization of magnetic units in spintronics devices. The common magnetic materials used in research and applications are based on Fe because the pure Fe metal is the best yet simple magnetic material from nature. Through systematic first-principles calculations, we explored the possibility to produce large PMA with ultrathin Fe on non-noble and non-magnetic Pb(001) substrate. Interestingly, huge magnetocrystalline anisotropy energy (MAE) of 7.6 meV was found in Pb/Fe/Pb(001) sandwich structure with only half monolayer Fe. Analysis of electronic structures reveals that the magnetic proximity effect at the interface is responsible for this significant enhancement of MAE. The MAE further increases to 13.6 meV with triply repeated capping Pb and intermediate Fe layers. Furthermore, the MAE can be tuned conveniently by charge injection.
• ### Two low-order nonconforming finite element methods for the Stokes flow in 3D(1710.06682)

March 16, 2018 math.NA
In this paper, we propose two low order nonconforming finite element methods (FEMs) for the three-dimensional Stokes flow that generalize the non-conforming FEM of Kouhia and Stenberg (1995, Comput. Methods Appl. Mech. Engrg.). The finite element spaces proposed in this paper consist of two globally continuous components (one piecewise affine and one enriched component) and one component that is continuous at the midpoints of interior faces. We prove that the discrete Korn inequality and a discrete inf-sup condition hold uniformly in the meshsize and also for a non-empty Neumann boundary. Based on these two results, we show the well-posedness of the discrete problem. Two counterexamples prove that there is no direct generalization of the Kouhia-Stenberg FEM to three space dimensions: The finite element space with one non-conforming and two conforming piecewise affine components does not satisfy a discrete inf-sup condition with piecewise constant pressure approximations, while finite element functions with two non-conforming and one conforming component do not satisfy a discrete Korn inequality.
• ### Probing the collective dynamics of nuclear spin bath in a rare-earth ion doped crystal(1803.03075)

March 9, 2018 quant-ph
Probing collective spin dynamics is a current challenge in the field of magnetic resonance spectroscopy and has important applications in material analysis and quantum information protocols. Recently, the rare-earth ion doped crystals are an attractive candidate for making long-lived quantum memory. Further enhancement of its performance would benefit from the direct knowledge on the dynamics of nuclear-spin bath in the crystal. Here we detect the collective dynamics of nuclear-spin bath located around the rare-earth ions in a crystal using dynamical decoupling spectroscopy method. From the measured spectrum, we analyze the configuration of the spin bath and characterize the flip-flop time between two correlated nuclear spins in a long time scale ($\sim$1s). Furthermore, we experimentally demonstrate that the rare-earth ions can serve as a magnetic quantum sensor for external magnetic field. These results suggest that the rare-earth ion is a useful probe for complex spin dynamics in solids and enable quantum sensing in the low-frequency regime, revealing promising possibilities for applications in diverse fields.
• ### High accuracy methods for eigenvalues of elliptic operators by nonconforming elements(1802.01896)

Feb. 6, 2018 math.NA
In this paper, three high-accuracy methods for eigenvalues of second order elliptic operators are proposed by using the nonconforming Crouzeix-Raviart(CR for short) element and the nonconforming enriched Crouzeix-Raviart(ECR for short) element. They are based on a crucial full one order superconvergence of the first order mixed Raviart-Thomas(RT for short) element. The main ingredient of such a superconvergence analysis is to employ a discrete Helmholtz decomposition of the difference between the canonical interpolation and the finite element solution of the RT element. In particular, it allows for some vital cancellation between terms in one key sum of boundary terms. Consequently, a full one order superconvergence follows from a special relation between the CR element and the RT element, and the equivalence between the ECR element and the RT element for these two nonconforming elements. These superconvergence results improve those in literature from a half order to a full one order for the RT element, the CR element and the ECR element. Based on the aforementioned superconvergence of the RT element, asymptotic expansions of eigenvalues are established and employed to achieve high accuracy extrapolation methods for these two nonconforming elements. In contrast to a classic analysis in literature, the novelty herein is to use not only the canonical interpolations of these nonconforming elements but also that of the RT element to analyze such asymptotic expansions. Based on the superconvergence of these nonconforming elements, asymptotically exact a posteriori error estimators of eigenvalues are constructed and analyzed for them. Finally, two post-processing methods are proposed to improve accuracy of approximate eigenvalues by employing these a posteriori error estimators.Numerical tests are provided to justify and compare the performance of the aforementioned methods.
• ### A Raman-Heterodyne Study of the Hyperfine Interaction of the Optically-Excited State $^5$D$_0$ of Eu$^{3+}$:Y$_2$SiO$_5$(1711.01814)

Jan. 28, 2018 quant-ph
The spin coherence time of $^{151}$Eu$^{3+}$ which substitutes the yttrium at site 1 in Y$_2$SiO$_5$ crystal has been extended to 6 hours in a recent work [\textit{Nature} \textbf{517}, 177 (2015)]. To make this long-lived spin coherence useful for optical quantum memory applications, we experimentally characterize the hyperfine interaction of the optically-excited state $^5$D$_0$ using Raman-heterodyne-detected nuclear magnetic resonance. The effective spin Hamiltonians for excited and ground state are fitted based on the experimental spectra obtained in 200 magnetic fields with various orientations. To show the correctness of the fitted parameters and potential application in quantum memory protocols, we also characterize the ground-state hyperfine interaction and predict the critical magnetic field which produces the 6-hour-long coherence time. The complete energy level structure for both the $^7$F$_0$ ground state and $^5$D$_0$ excited state at the critical magnetic field are obtained. These results enable the design of quantum memory protocols and the optimization of optical pumping strategy for realization of photonic quantum memory with hour-long lifetime.
• ### On the structure of cyclotomic nilHecke algebras(1709.08760)

Dec. 28, 2017 math.RT
In this paper we study the structure of the cyclotomic nilHecke algebras $\HH_{\ell,n}^{(0)}$, where $\ell,n\in\N$. We construct a monomial basis for $\HH_{\ell,n}^{(0)}$ which verifies a conjecture of Mathas. We show that the graded basic algebra of $\HH_{\ell,n}^{(0)}$ is commutative and hence isomorphic to the center $Z$ of $\HH_{\ell,n}^{(0)}$. We further prove that $\HH_{\ell,n}^{(0)}$ is isomorphic to the full matrix algebra over $Z$ and construct an explicit basis for the center $Z$. We also construct a complete set of pairwise orthogonal primitive idempotents of $\HH_{\ell,n}^{(0)}$. Finally, we present a new homogeneous symmetrizing form $\Tr$ on $\HH_{\ell,n}^{(0)}$ by explicitly specifying its values on a given homogeneous basis of $\HH_{\ell,n}^{(0)}$ and show that it coincides with Shan--Varagnolo--Vasserot's symmetrizing form $\Tr^{\text{SVV}}$ on $\HH_{\ell,n}^{(0)}$.
• ### Multigrid Methods for Hellan-Herrmann-Johnson Mixed Method of Kirchhoff Plate Bending Problems(1508.04289)

Dec. 24, 2017 math.NA
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.
• ### Record High Magnetic Anisotropy in Chemically Engineered Iridium Dimer(1709.01878)

Exploring giant magnetic anisotropy in small magnetic nanostructures is of both fundamental interest and technological merit for information storage. To prevent spin flipping at room temperature due to thermal fluctuation, large magnetic anisotropy energy (MAE) over 50 meV in magnetic nanostructure is desired for practical applications. We chose one of the smallest magnetic nanostructures-Ir2 dimer, to investigate its magnetic properties and explore possible approach to engineer the magnetic anisotropy. Through systematic first-principles calculations, we found that the Ir2 dimer already possesses giant MAE of 77 meV. We proposed an effective way to enhance the MAE of the Ir2 dimer to 223~294 meV by simply attaching a halogen atom at one end of the Ir-Ir bond. The underlying mechanism for the record high MAE is attributed to the modification of the energy diagram of the Ir2 dimer by the additional halogen-Ir bonding, which alters the spin-orbit coupling Hamiltonian and hence the magnetic anisotropy. Our strategy can be generalized to design other magnetic molecules or clusters with giant magnetic anisotropy.
• ### On the centers of cyclotomic quiver Hecke algebras(1707.02534)

July 13, 2017 math.RT
Let $n\in\mathbb{N}$ and $K$ be any field. For any symmetric generalized Cartan matrix $A$, any $\beta$ in the positive root lattice with height $n$ and any integral dominant weight $\Lambda$, one can associate a quiver Hecke algebras $R_{\beta}(K)$ and its cyclotomic quotient $R_{\beta}^{\Lambda}(K)$ over $K$. It has been conjectured that the natural map from $R_{\beta}(K)$ to $R_{\beta}^{\Lambda}(K)$ maps the center of $R_{\beta}(K)$ surjectively onto the center of $R_{\beta}^{\Lambda}(K)$. A similar conjecture claims that the center of the affine Hecke algebra of type $A$ maps surjectively onto the center of its cyclotomic quotient---the cyclotomic Hecke algebra $H_n^{\Lambda}$ of type $G(\ell,1,n)$ over $K$. In this paper, we prove these two conjectures affirmatively. As a consequence, we show that the center of $H_n^{\Lambda}$ is stable under base change and it has dimension equal to the number of $\ell$-partitions of $n$. Finally, as a byproduct, we also verify a conjecture of Shan, Varagnolo and Vasserot on the grading structure of the center of $R_{\beta}^{\Lambda}(K)$.
• ### Experimental Demonstration of an Electride as a 2D Material(1706.02774)

June 8, 2017 cond-mat.mtrl-sci
Because of their loosely bound electrons, electrides offer physical properties useful in chemical synthesis and electronics. For these applications and others, nano-sized electrides offer advantages, but to-date no electride has been synthesized as a nanomaterial. We demonstrate experimentally that Ca$_2$N, a layered electride in which layers of atoms are separated by layers of a 2D electron gas (2DEG), can be exfoliated into two-dimensional (2D) nanosheets using liquid exfoliation. The 2D flakes are stable in a nitrogen atmosphere or in select organic solvents for at least one month. Electron microscopy and elemental analysis reveal that the 2D flakes retain the crystal structure and stoichiometry of the parent 3D Ca$_2$N. In addition, the 2D flakes exhibit metallic character and an optical response that agrees with DFT calculations. Together these findings suggest that the 2DEG is preserved in the 2D material. With this work, we bring electrides into the nano-regime and experimentally demonstrate a 2D electride, Ca$_2$N.
• ### Fayers' conjecture and the socles of cyclotomic Weyl modules(1602.06631)

May 29, 2017 math.QA, math.RT, math.GR
Gordon James proved that the socle of a Weyl module of a classical Schur algebra is a sum of simple modules labelled by $p$-restricted partitions. We prove an analogue of this result in the very general setting of "Schur pairs". As an application we show that the socle of a Weyl module of a cyclotomic $q$-Schur algebra is a sum of simple modules labelled by Kleshchev multipartitions and we use this result to prove a conjecture of Fayers that leads to an efficient LLT algorithm for the higher level cyclotomic Hecke algebras of type $A$. Finally, we prove a cyclotomic analogue of the Carter-Lusztig theorem.
• ### Residual-Based A Posteriori Error Estimates for Symmetric Conforming Mixed Finite Elements for Linear Elasticity Problems(1705.04106)

May 11, 2017 math.NA
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.
• ### Nodal Finite Element de Rham Complexes(1611.02558)

March 21, 2017 math.NA
We construct 2D and 3D finite element de Rham sequences of arbitrary polynomial degree with extra smoothness. Some of these elements have nodal degrees of freedom (DoFs), and can be considered as generalisations of scalar Hermite and Lagrange elements. Using the nodal values, the number of global degrees of freedom is reduced compared with the classical N\'{e}d\'{e}lec and Brezzi-Douglas-Marini (BDM) finite elements, and the basis functions are more canonical and easier to write. Our finite elements for ${H}(\mathrm{div})$ with regularity $r=2$ coincide with the nonstandard elements given by Stenberg \cite{Stenberg2010}. We show how regularity decreases in the finite element complexes, so that they branch into known complexes. The standard de Rham complexes of Whitney forms and their higher order version can be considered as the family with lowest regularity. The construction of the new families is motivated by finite element systems (FES).
• ### High Order Hierarchical Divergence-free Constrained Transport $H(div)$ Finite Element Method for Magnetic Induction Equation(1702.01180)

Feb. 3, 2017 math.NA
In this paper, we will use the interior functions of an hierarchical basis for high order $BDM_p$ elements to enforce the divergence-free condition of a magnetic field $B$ approximated by the H(div) $BDM_p$ basis. The resulting constrained finite element method can be used to solve magnetic induction equation in MHD equations. The proposed procedure is based on the fact that the scalar $(p-1)$-th order polynomial space on each element can be decomposed as an orthogonal sum of the subspace defined by the divergence of the interior functions of the $p$-th order $BDM_p$ basis and the constant function. Therefore, the interior functions can be used to remove element-wise all higher order terms except the constant in the divergence error of the finite element solution of $B$-field. The constant terms from each element can be then easily corrected using a first order H(div) basis globally. Numerical results for a 3-D magnetic induction equation show the effectiveness of the proposed method in enforcing divergence-free condition of the magnetic field.
• ### A Butterfly-Based Direct Integral Equation Solver Using Hierarchical LU Factorization for Analyzing Scattering from Electrically Large Conducting Objects(1610.00042)

Jan. 15, 2017 math.NA
A butterfly-based direct combined-field integral equation (CFIE) solver for analyzing scattering from electrically large, perfect electrically conducting objects is presented. The proposed solver leverages the butterfly scheme to compress blocks of the hierarchical LU-factorized discretized CFIE operator and uses randomized butterfly reconstruction schemes to expedite the factorization. The memory requirements and computational cost of the direct butterfly-CFIE solver scale as $O(N\mathrm{log}^2N)$ and $O(N^{1.5}\mathrm{log}N)$, respectively. These scaling estimates permit significant memory and CPU savings when compared to those realized by low-rank (LR) decomposition-based solvers. The efficacy and accuracy of the proposed solver are demonstrated through its application to the analysis of scattering from canonical and realistic objects involving up to 14 million unknowns.
• ### On involutions in Weyl groups(1609.08494)

Nov. 10, 2016 math.CO, math.GR
Let $(W,S)$ be a Coxeter system and $\ast$ be an automorphism of $W$ with order $\leq 2$ such that $s^{\ast}\in S$ for any $s\in S$. Let $I_{\ast}$ be the set of twisted involutions relative to $\ast$ in $W$. In this paper we consider the case when $\ast=\text{id}$ and study the braid $I_\ast$-transformations between the reduced $I_\ast$-expressions of involutions. If $W$ is the Weyl group of type $B_n$ or $D_n$, we explicitly describe a finite set of basic braid $I_\ast$-transformations for all $n$ simultaneously, and show that any two reduced $I_\ast$-expressions for a given involution can be transformed into each other through a series of basic braid $I_\ast$-transformations. In both cases, these basic braid $I_\ast$-transformations consist of the usual basic braid transformations plus some natural "right end transformations" and plus exactly one extra transformation. The main result generalizes our previous work for the Weyl group of type $A_{n}$.
• ### Stabilized mixed finite element methods for linear elasticity on simplicial grids in $\mathbb{R}^{n}$(1512.03998)

Oct. 27, 2016 math.NA
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.
• ### A Penalized Crouzeix-Raviart Element Method for Second Order Elliptic Eigenvalue Problems(1608.03987)

Aug. 13, 2016 math.NA
In this paper we propose a penalized Crouzeix-Raviart element method for eigenvalue problems of second order elliptic operators. The key idea is to add a penalty term to tune the local approximation property and the global continuity property of the discrete eigenfunctions. The feature of this method is that by adjusting the penalty parameter, the resulted discrete eigenvalues can be in a state of "chaos", and consequently a large portion of them can be reliable and approximate the exact ones with high accuracy. Furthermore, we design an algorithm to select such a quasi-optimal penalty parameter. Finally, we provide numerical tests to demonstrate the performance of the proposed method.
• ### Fast Auxiliary Space Preconditioner for Linear Elasticity in Mixed Form(1604.02568)

May 21, 2016 math.NA
A block diagonal preconditioner with the minimal residual method and a block triangular preconditioner with the generalized minimal residual method are developed for Hu-Zhang mixed finite element methods of linear elasticity. They are based on a new stability result of the saddle point system in mesh-dependent norms. The mesh-dependent norm for the stress corresponds to the mass matrix which is easy to invert while the displacement it is spectral equivalent to Schur complement. A fast auxiliary space preconditioner based on the $H^1$ conforming linear element of the linear elasticity problem is then designed for solving the Schur complement. For both diagonal and triangular preconditioners, it is proved that the conditioning numbers of the preconditioned systems are bounded above by a constant independent of both the crucial Lam\'e constant and the mesh-size. Numerical examples are presented to support theoretical results. As byproducts, a new stabilized low order mixed finite element method is proposed and analyzed and superconvergence results of Hu-Zhang element are obtained.