Vortex structures of rotating spin-orbit coupled Bose-Einstein condensates

We consider the quasi-2D two-component Bose-Einstein condensates with Rashba spin-orbit (SO) coupling in a rotating trap. An external Zeeman term favoring spin polarization along the radial direction is also considered, which has the same form as the non-canonical part of the mechanical angular momentum. The rotating condensate exhibits rich structures as varying the strengths of trapping potential and interaction. With a strong trapping potential, the condensate exhibits a half-quantum vortex-lattice configuration. Such a configuration is driven to the normal one by introducing the external radial Zeeman field. In the case of a weak trap potential, the condensate exhibits a multi-domain pattern of plane-wave states under the external radial Zeeman field.Comment: 8 pages, 7 figures, two figures are adde

Understanding and Improving the Wang-Landau Algorithm

We present a mathematical analysis of the Wang-Landau algorithm, prove its convergence, identify sources of errors and strategies for optimization. In particular, we found the histogram increases uniformly with small fluctuation after a stage of initial accumulation, and the statistical error is found to scale as $\sqrt{\ln f}$ with the modification factor $f$. This has implications for strategies for obtaining fast convergence.Comment: 4 pages, 2 figures, to appear in Phys. Rev.

Critical domain-wall dynamics of model B

With Monte Carlo methods, we simulate the critical domain-wall dynamics of model B, taking the two-dimensional Ising model as an example. In the macroscopic short-time regime, a dynamic scaling form is revealed. Due to the existence of the quasi-random walkers, the magnetization shows intrinsic dependence on the lattice size $L$. A new exponent which governs the $L$-dependence of the magnetization is measured to be $\sigma=0.243(8)$.Comment: 10pages, 4 figure

A non-local vector calculus,non-local volume-constrained problems,and non-local balance laws

A vector calculus for nonlocal operators is developed, including the definition of nonlocal divergence, gradient, and curl operators and the derivation of the corresponding adjoints operators. Nonlocal analogs of several theorems and identities of the vector calculus for differential operators are also presented. Relationships between the nonlocal operators and their differential counterparts are established, first in a distributional sense and then in a weak sense by considering weighted integrals of the nonlocal adjoint operators. The nonlocal calculus gives rise to volume-constrained problems that are analogous to elliptic boundary-value problems for differential operators; this is demonstrated via some examples. Another application is posing abstract nonlocal balance laws and deriving the corresponding nonlocal field equations

Constraining the Origin of Local Positrons with HAWC TeV Gamma-Ray Observations of Two Nearby Pulsar Wind Nebulae

The HAWC Gamma-Ray Observatory has reported the discovery of TeV gamma-ray emission extending several degrees around the positions of Geminga and B0656+14 pulsars. Assuming these gamma rays are produced by inverse Compton scattering off low-energy photons in electron halos around the pulsars, we determine the diffusion of electrons and positrons in the local interstellar medium. We will present the morphological and spectral studies of these two VHE gamma-ray sources and the derived positron spectrum at Earth.Comment: Presented at the 35th International Cosmic Ray Conference (ICRC2017), Bexco, Busan, Korea. See arXiv:1708.02572 for all HAWC contribution

Creep motion of a domain wall in the two-dimensional random-field Ising model with a driving field

With Monte Carlo simulations, we study the creep motion of a domain wall in the two-dimensional random-field Ising model with a driving field. We observe the nonlinear fieldvelocity relation, and determine the creep exponent {\mu}. To further investigate the universality class of the creep motion, we also measure the roughness exponent {\zeta} and energy barrier exponent {\psi} from the zero-field relaxation process. We find that all the exponents depend on the strength of disorder.Comment: 5 pages, 4 figure

Parallel Self-Consistent-Field Calculations via Chebyshev-Filtered Subspace Acceleration

Solving the Kohn-Sham eigenvalue problem constitutes the most computationally expensive part in self-consistent density functional theory (DFT) calculations. In a previous paper, we have proposed a nonlinear Chebyshev-filtered subspace iteration method, which avoids computing explicit eigenvectors except at the first SCF iteration. The method may be viewed as an approach to solve the original nonlinear Kohn-Sham equation by a nonlinear subspace iteration technique, without emphasizing the intermediate linearized Kohn-Sham eigenvalue problem. It reaches self-consistency within a similar number of SCF iterations as eigensolver-based approaches. However, replacing the standard diagonalization at each SCF iteration by a Chebyshev subspace filtering step results in a significant speedup over methods based on standard diagonalization. Here, we discuss an approach for implementing this method in multi-processor, parallel environment. Numerical results are presented to show that the method enables to perform a class of highly challenging DFT calculations that were not feasible before