An advanced meshless method for time fractional diffusion equation

Recently, because of the new developments in sustainable engineering and renewable energy, which are usually governed by a series of fractional partial differential equations (FPDEs), the numerical modelling and simulation for fractional calculus are attracting more and more attention from researchers. The current dominant numerical method for modeling FPDE is Finite Difference Method (FDM), which is based on a pre-defined grid leading to inherited issues or shortcomings including difficulty in simulation of problems with the complex problem domain and in using irregularly distributed nodes. Because of its distinguished advantages, the meshless method has good potential in simulation of FPDEs. This paper aims to develop an implicit meshless collocation technique for FPDE. The discrete system of FPDEs is obtained by using the meshless shape functions and the meshless collocation formulation. The stability and convergence of this meshless approach are investigated theoretically and numerically. The numerical examples with regular and irregular nodal distributions are used to validate and investigate accuracy and efficiency of the newly developed meshless formulation. It is concluded that the present meshless formulation is very effective for the modeling and simulation of fractional partial differential equations

Fermionic concurrence in the extended Hubbard dimer

In this paper, we introduce and study the fermionic concurrence in a two-site extended Hubbard model. Its behaviors both at the ground state and finite temperatures as function of Coulomb interaction $U$ (on-site) and $V$ (nearest-neighbor) are obtained analytically and numerically. We also investigate the change of the concurrence under a nonuniform field, including local potential and magnetic field, and find that the concurrence can be modulated by these fields.Comment: 5 pages, 7 figure

Efficient electronic entanglement concentration assisted with single mobile electron

We present an efficient entanglement concentration protocol (ECP) for mobile electrons with charge detection. This protocol is quite different from other ECPs for one can obtain a maximally entangled pair from a pair of less-entangled state and a single mobile electron with a certain probability. With the help of charge detection, it can be repeated to reach a higher success probability. It also does not need to know the coefficient of the original less-entangled states. All these advantages may make this protocol useful in current distributed quantum information processing.Comment: 6pages, 3figure

Neutrino masses, leptogenesis and dark matter in hybrid seesaw

We suggest a hybrid seesaw model where relatively ``light''right-handed neutrinos give no contribution to the neutrino mass matrix due to a special symmetry. This allows their Yukawa couplings to the standard model particles to be relatively strong, so that the standard model Higgs boson can decay dominantly to a left and a right-handed neutrino, leaving another stable right-handed neutrino as cold dark matter. In our model neutrino masses arise via the type-II seesaw mechanism, the Higgs triplet scalars being also responsible for the generation of the matter-antimatter asymmetry via the leptogenesis mechanism.Comment: 4 page

Fast Monte Carlo Simulation for Patient-specific CT/CBCT Imaging Dose Calculation

Recently, X-ray imaging dose from computed tomography (CT) or cone beam CT (CBCT) scans has become a serious concern. Patient-specific imaging dose calculation has been proposed for the purpose of dose management. While Monte Carlo (MC) dose calculation can be quite accurate for this purpose, it suffers from low computational efficiency. In response to this problem, we have successfully developed a MC dose calculation package, gCTD, on GPU architecture under the NVIDIA CUDA platform for fast and accurate estimation of the x-ray imaging dose received by a patient during a CT or CBCT scan. Techniques have been developed particularly for the GPU architecture to achieve high computational efficiency. Dose calculations using CBCT scanning geometry in a homogeneous water phantom and a heterogeneous Zubal head phantom have shown good agreement between gCTD and EGSnrc, indicating the accuracy of our code. In terms of improved efficiency, it is found that gCTD attains a speed-up of ~400 times in the homogeneous water phantom and ~76.6 times in the Zubal phantom compared to EGSnrc. As for absolute computation time, imaging dose calculation for the Zubal phantom can be accomplished in ~17 sec with the average relative standard deviation of 0.4%. Though our gCTD code has been developed and tested in the context of CBCT scans, with simple modification of geometry it can be used for assessing imaging dose in CT scans as well.Comment: 18 pages, 7 figures, and 1 tabl

A Cosmological Model with Dark Spinor Source

In this paper, we discuss the system of Friedman-Robertson-Walker metric coupling with massive nonlinear dark spinors in detail, where the thermodynamic movement of spinors is also taken into account. The results show that, the nonlinear potential of the spinor field can provide a tiny negative pressure, which resists the Universe to become singular. The solution is oscillating in time and closed in space, which approximately takes the following form g_{\mu\nu}=\bar R^2(1-\delta\cos t)^2\diag(1,-1,-\sin^2r ,-\sin^2r \sin^2\theta), with $\bar R= (1\sim 2)\times 10^{12}$ light year, and $\delta=0.96\sim 0.99$. The present time is about $t\sim 18^\circ$.Comment: 13 pages, no figure, to appear in IJMP

Scattering on two Aharonov-Bohm vortices with opposite fluxes

The scattering of an incident plane wave on two Aharonov-Bohm vortices with opposite fluxes is considered in detail. The presence of the vortices imposes non-trivial boundary conditions for the partial waves on a cut joining the two vortices. These conditions result in an infinite system of equations for scattering amplitudes between incoming and outgoing partial waves, which can be solved numerically. The main focus of the paper is the analytic determination of the scattering amplitude in two limits, the small flux limit and the limit of small vortex separation. In the latter limit the dominant contribution comes from the S-wave amplitude. Calculating it, however, still requires solving an infinite system of equations, which is achieved by the Riemann-Hilbert method. The results agree well with the numerical calculations