On the formation of current sheets in response to the compression or expansion of a potential magnetic field

The compression or expansion of a magnetic field that is initially potential is considered. It was recently suggested by Janse & Low [2009, ApJ, 690, 1089] that, following the volumetric deformation, the relevant lowest energy state for the magnetic field is another potential magnetic field that in general contains tangential discontinuities (current sheets). Here we examine this scenario directly using a numerical relaxation method that exactly preserves the topology of the magnetic field. It is found that of the magnetic fields discussed by Janse & Low, only those containing magnetic null points develop current singularities during an ideal relaxation, while the magnetic fields without null points relax toward smooth force-free equilibria with finite non-zero current.Comment: Accepted for publication in Ap

Gunrock: A High-Performance Graph Processing Library on the GPU

For large-scale graph analytics on the GPU, the irregularity of data access and control flow, and the complexity of programming GPUs have been two significant challenges for developing a programmable high-performance graph library. "Gunrock", our graph-processing system designed specifically for the GPU, uses a high-level, bulk-synchronous, data-centric abstraction focused on operations on a vertex or edge frontier. Gunrock achieves a balance between performance and expressiveness by coupling high performance GPU computing primitives and optimization strategies with a high-level programming model that allows programmers to quickly develop new graph primitives with small code size and minimal GPU programming knowledge. We evaluate Gunrock on five key graph primitives and show that Gunrock has on average at least an order of magnitude speedup over Boost and PowerGraph, comparable performance to the fastest GPU hardwired primitives, and better performance than any other GPU high-level graph library.Comment: 14 pages, accepted by PPoPP'16 (removed the text repetition in the previous version v5

Ballistic-Ohmic quantum Hall plateau transition in graphene pn junction

Recent quantum Hall experiments conducted on disordered graphene pn junction provide evidence that the junction resistance could be described by a simple Ohmic sum of the n and p mediums' resistances. However in the ballistic limit, theory predicts the existence of chirality-dependent quantum Hall plateaus in a pn junction. We show that two distinctively separate processes are required for this ballistic-Ohmic plateau transition, namely (i) hole/electron Landau states equilibration and (ii) valley iso-spin dilution of the incident Landau edge state. These conclusions are obtained by a simple scattering theory argument, and confirmed numerically by performing ensembles of quantum magneto-transport calculations on a 0.1um-wide disordered graphene pn junction within the tight-binding model. The former process is achieved by pn interface roughness, where a pn interface disorder with a root-mean-square roughness of 10nm was found to suffice under typical experimental conditions. The latter process is mediated by extrinsic edge roughness for an armchair edge ribbon and by intrinsic localized intervalley scattering centers at the edge of the pn interface for a zigzag ribbon. In light of these results, we also examine why higher Ohmic type plateaus are less likely to be observable in experiments.Comment: 9 pages, 6 figure

Theoretically Efficient Parallel Graph Algorithms Can Be Fast and Scalable

There has been significant recent interest in parallel graph processing due to the need to quickly analyze the large graphs available today. Many graph codes have been designed for distributed memory or external memory. However, today even the largest publicly-available real-world graph (the Hyperlink Web graph with over 3.5 billion vertices and 128 billion edges) can fit in the memory of a single commodity multicore server. Nevertheless, most experimental work in the literature report results on much smaller graphs, and the ones for the Hyperlink graph use distributed or external memory. Therefore, it is natural to ask whether we can efficiently solve a broad class of graph problems on this graph in memory. This paper shows that theoretically-efficient parallel graph algorithms can scale to the largest publicly-available graphs using a single machine with a terabyte of RAM, processing them in minutes. We give implementations of theoretically-efficient parallel algorithms for 20 important graph problems. We also present the optimizations and techniques that we used in our implementations, which were crucial in enabling us to process these large graphs quickly. We show that the running times of our implementations outperform existing state-of-the-art implementations on the largest real-world graphs. For many of the problems that we consider, this is the first time they have been solved on graphs at this scale. We have made the implementations developed in this work publicly-available as the Graph-Based Benchmark Suite (GBBS).Comment: This is the full version of the paper appearing in the ACM Symposium on Parallelism in Algorithms and Architectures (SPAA), 201

A_4 Symmetry and Lepton Masses and Mixing

Stimulated by Ma's idea which explains the tribimaximal neutrino mixing by assuming an A_4 flavor symmetry, a lepton mass matrix model is investigated. A Frogatt-Nielsen type model is assumed, and the flavor structures of the masses and mixing are caused by the VEVs of SU(2)_L-singlet scalars \phi_i^u and \phi_i^d (i=1,2,3), which are assigned to {\bf 3} and ({\bf 1}, {\bf 1}',{\bf 1}'') of A_4, respectively.Comment: 13 pages including 1 table, errors in Sec.7 correcte

On electroweak baryogenesis in the littlest Higgs model with T parity

We study electroweak baryogenesis within the framework of the littlest Higgs model with T parity. This model has shown characteristics of a strong first-order electroweak phase transition, which is conducive to baryogenesis in the early Universe. In the T parity symmetric theory, there are two gauge sectors, viz., the T-even and the T-odd ones. We observe that the effect of the T-parity symmetric interactions between the T-odd and the T-even gauge bosons on gauge-higgs energy functional is quite small, so that these two sectors can be taken to be independent. The T-even gauge bosons behave like the Standard Model gauge bosons, whereas the T-odd ones are instrumental in stabilizing the Higgs mass. For the T-odd gauge bosons in the symmetric and asymmetric phases and for the T-even gauge bosons in the asymmetric phase, we obtain, using the formalism of Arnold and McLerran, very small values of the ratio, (Baryon number violation rate/Universe expansion rate). We observe that this result, in conjunction with the scenario of inverse phase transition in the present work and the value of the ratio obtained from the lattice result of sphaleron transition rate in the symmetric phase, can provide us with a plausible baryogenesis scenario.Comment: 13 pages, 2 figures, published version, references modifie

Protostellar Collapse with Various Metallicities

The thermal and chemical evolution of gravitationally collapsing protostellar clouds is investigated, focusing attention on their dependence on metallicity. Calculations are carried out for a range of metallicities spanning the local interstellar value to zero. During the time when clouds are transparent to continuous radiation, the temperatures are higher for those with lower metallicity, reflecting lower radiative ability. However, once the clouds become opaque, in the course of the adiabatic contraction of the transient cores, their evolutionary trajectories in the density-temperature plane converge to a unique curve that is determined by only physical constants. The trajectories coincide with each other thereafter. Consequently, the size of the stellar core at the formation is the same regardless of the gas composition of the parent cloud.Comment: 30 pages. The Astrophysical Journal, 533, in pres

Mushroom-Shaped Structures as Tracers of Buoyant Flow in the Galactic Disk

Recent HI emission observations of the Southern Galactic hemisphere have revealed a mushroom-like structure extending from z=-70 to -450 pc, composed of a stem and a cap. Similar structures occur in three-dimensional simulations of a dynamic galactic disk driven by isolated and clustered supernovae. Using these simulations, we show that hot gas in the Galactic disk that is not evacuated through chimneys expands into the cooler gas of the thick disk, forming mushroom-shaped structures. This new class of objects traces buoyant flow of hot gas into the thick disk.Comment: Accepted for publication in ApJ Letters. Latex manuscript, 3 figures (4 postsript files

Star Formation in Isolated Disk Galaxies. I. Models and Characteristics of Nonlinear Gravitational Collapse

We model gravitational collapse leading to star formation in a wide range of isolated disk galaxies using a three-dimensional, smoothed particle hydrodynamics code. The model galaxies include a dark matter halo and a disk of stars and isothermal gas. Absorbing sink particles are used to directly measure the mass of gravitationally collapsing gas. They reach masses characteristic of stellar clusters. In this paper, we describe our galaxy models and numerical methods, followed by an investigation of the gravitational instability in these galaxies. Gravitational collapse forms star clusters with correlated positions and ages, as observed, for example, in the Large Magellanic Cloud. Gravitational instability alone acting in unperturbed galaxies appears sufficient to produce flocculent spiral arms, though not more organized patterns. Unstable galaxies show collapse in thin layers in the galactic plane; associated dust will form thin dust lanes in those galaxies, in agreement with observations. (abridged)Comment: 49 pages, 22 figures, to appear in ApJ (July, 2005), version with high quality color images can be fond in http://research.amnh.org/~yuexing/astro-ph/0501022.pd

An Interstellar Conduction Front Within a Wolf-Rayet Ring Nebula Observed with the GHRS

With the High Resolution Spectrograph aboard the Hubble Space Telescope we obtained high signal-to-noise (S/N > 200 - 600 per 17 km/s resolution element) spectra of narrow absorption lines toward the Wolf-Rayet star HD 50896. The ring nebula S308 that surrounds this star is thought to be caused by a pressure-driven bubble bounded by circumstellar gas (most likely from a red supergiant or luminous blue variable progenitor) pushed aside by a strong stellar wind. Our observation has shown for the first time that blueshifted (approximately 70 km/s relative to the star) absorption components of C IV and N V arise in a conduction front between the hot interior of the bubble and the cold shell of swept-up material. These lines set limits on models of the conduction front. Nitrogen in the shell appears to be overabundant by a factor ~10. The P Cygni profiles of N V and C IV are variable, possibly due to a suspected binary companion to HD 50896.Comment: 32 pages, Latex, to appear in the Astrophysical Journal, April, 199