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 new class of (2+1)(2+1)-d topological superconductor with Z8\mathbb{Z}_8 topological classification

The classification of topological states of matter depends on spatial dimension and symmetry class. For non-interacting topological insulators and superconductors the topological classification is obtained systematically and nontrivial topological insulators are classified by either integer or $Z_2$. The classification of interacting topological states of matter is much more complicated and only special cases are understood. In this paper we study a new class of topological superconductors in $(2+1)$ dimensions which has time-reversal symmetry and a $\mathbb{Z}_2$ spin conservation symmetry. We demonstrate that the superconductors in this class is classified by $\mathbb{Z}_8$ when electron interaction is considered, while the classification is $\mathbb{Z}$ without interaction.Comment: 5 pages main text and 3 pages appendix. 1 figur

Top-N Recommendation on Graphs

Recommender systems play an increasingly important role in online applications to help users find what they need or prefer. Collaborative filtering algorithms that generate predictions by analyzing the user-item rating matrix perform poorly when the matrix is sparse. To alleviate this problem, this paper proposes a simple recommendation algorithm that fully exploits the similarity information among users and items and intrinsic structural information of the user-item matrix. The proposed method constructs a new representation which preserves affinity and structure information in the user-item rating matrix and then performs recommendation task. To capture proximity information about users and items, two graphs are constructed. Manifold learning idea is used to constrain the new representation to be smooth on these graphs, so as to enforce users and item proximities. Our model is formulated as a convex optimization problem, for which we need to solve the well-known Sylvester equation only. We carry out extensive empirical evaluations on six benchmark datasets to show the effectiveness of this approach.Comment: CIKM 201

Braiding statistics approach to symmetry-protected topological phases

We construct a 2D quantum spin model that realizes an Ising paramagnet with gapless edge modes protected by Ising symmetry. This model provides an example of a "symmetry-protected topological phase." We describe a simple physical construction that distinguishes this system from a conventional paramagnet: we couple the system to a Z_2 gauge field and then show that the \pi-flux excitations have different braiding statistics from that of a usual paramagnet. In addition, we show that these braiding statistics directly imply the existence of protected edge modes. Finally, we analyze a particular microscopic model for the edge and derive a field theoretic description of the low energy excitations. We believe that the braiding statistics approach outlined in this paper can be generalized to a large class of symmetry-protected topological phases.Comment: 17 pages, 12 figures, reorganized section V, added a referenc

Relativistic description of J/\psi dissociation in hot matter

The mass spectra and binding radii of heavy quark bound states are studied on the basis of the reduced Bethe-Salpeter equation. The critical values of screening masses for $c\bar{c}$ and $b\bar{b}$ bound states at a finite temperature are obtained and compared with the previous results given by non-relativistic models.Comment: 13 latex pages, 2 figure

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