Diffusion in a multi-component Lattice Boltzmann Equation model

Diffusion phenomena in a multiple component lattice Boltzmann Equation (LBE) model are discussed in detail. The mass fluxes associated with different mechanical driving forces are obtained using a Chapman-Enskog analysis. This model is found to have correct diffusion behavior and the multiple diffusion coefficients are obtained analytically. The analytical results are further confirmed by numerical simulations in a few solvable limiting cases. The LBE model is established as a useful computational tool for the simulation of mass transfer in fluid systems with external forces.Comment: To appear in Aug 1 issue of PR

A Lattice Boltzmann method for simulations of liquid-vapor thermal flows

We present a novel lattice Boltzmann method that has a capability of simulating thermodynamic multiphase flows. This approach is fully thermodynamically consistent at the macroscopic level. Using this new method, a liquid-vapor boiling process, including liquid-vapor formation and coalescence together with a full coupling of temperature, is simulated for the first time.Comment: one gzipped tar file, 19 pages, 4 figure

Multi-component lattice-Boltzmann model with interparticle interaction

A previously proposed [X. Shan and H. Chen, Phys. Rev. E {\bf 47}, 1815, (1993)] lattice Boltzmann model for simulating fluids with multiple components and interparticle forces is described in detail. Macroscopic equations governing the motion of each component are derived by using Chapman-Enskog method. The mutual diffusivity in a binary mixture is calculated analytically and confirmed by numerical simulation. The diffusivity is generally a function of the concentrations of the two components but independent of the fluid velocity so that the diffusion is Galilean invariant. The analytically calculated shear kinematic viscosity of this model is also confirmed numerically.Comment: 18 pages, compressed and uuencoded postscript fil

Lattice Boltzmann method with self-consistent thermo-hydrodynamic equilibria

Lattice kinetic equations incorporating the effects of external/internal force fields via a shift of the local fields in the local equilibria, are placed within the framework of continuum kinetic theory. The mathematical treatment reveals that, in order to be consistent with the correct thermo-hydrodynamical description, temperature must also be shifted, besides momentum. New perspectives for the formulation of thermo-hydrodynamic lattice kinetic models of non-ideal fluids are then envisaged. It is also shown that on the lattice, the definition of the macroscopic temperature requires the inclusion of new terms directly related to discrete effects. The theoretical treatment is tested against a controlled case with a non ideal equation of state.Comment: 10 pages, 1 figur

Determining the Mass of Dark Matter Particles with Direct Detection Experiments

In this article I review two data analysis methods for determining the mass (and eventually the spin-independent cross section on nucleons) of Weakly Interacting Massive Particles with positive signals from direct Dark Matter detection experiments: a maximum likelihood analysis with only one experiment and a model-independent method requiring at least two experiments. Uncertainties and caveats of these methods will also be discussed.Comment: 24 pages, 10 figures, 1 reference added, typos fixed, published version, to appear in the NJP Focus Issue on "Dark Matter and Particle Physics

Families of Graphs With Chromatic Zeros Lying on Circles

We define an infinite set of families of graphs, which we call $p$-wheels and denote $(Wh)^{(p)}_n$, that generalize the wheel ($p=1$) and biwheel ($p=2$) graphs. The chromatic polynomial for $(Wh)^{(p)}_n$ is calculated, and remarkably simple properties of the chromatic zeros are found: (i) the real zeros occur at $q=0,1,...p+1$ for $n-p$ even and $q=0,1,...p+2$ for $n-p$ odd; and (ii) the complex zeros all lie, equally spaced, on the unit circle $|q-(p+1)|=1$ in the complex $q$ plane. In the $n \to \infty$ limit, the zeros on this circle merge to form a boundary curve separating two regions where the limiting function $W(\{(Wh)^{(p)}\},q)$ is analytic, viz., the exterior and interior of the above circle. Connections with statistical mechanics are noted.Comment: 8 pages, Late

Effective continuous model for surface states and thin films of three dimensional topological insulators

Two-dimensional effective continuous models are derived for the surface states and thin films of the three-dimensional topological insulator (3DTI). Starting from an effective model for 3DTI based on the first principles calculation [Zhang \emph{et al}, Nat. Phys. 5, 438 (2009)], we present solutions for both the surface states in a semi-infinite boundary condition and in the thin film with finite thickness. An effective continuous model was derived for surface states and the thin film 3DTI. The coupling between opposite topological surfaces and structure inversion asymmetry (SIA) give rise to gapped Dirac hyperbolas with Rashba-like splittings in energy spectrum. Besides, the SIA leads to asymmetric distributions of wavefunctions along the film growth direction, making some branches in the energy spectra much harder than others to be probed by light. These features agree well with the recent angle-resolved photoemission spectra of Bi$_{2}$Se $_{3}$ films grown on SiC substrate [Zhang et al, arXiv: 0911.3706]. More importantly, we use the effective model to fit the experimental data and determine the model parameters. The result indicates that the thin film Bi$_{2}$Se$_{3}$ lies in quantum spin Hall region based on the calculation of the Chern number and the $Z_{2}$ invariant. In addition, strong SIA always intends to destroy the quantum spin Hall state.Comment: 12 pages, 7 figures, references are update

The specialization problem and the completeness of unfolding

We discuss the problem of specializing a definite program with respect to sets of positive and negative examples, following Bostrom and Idestam-Almquist. This problem is very relevant in the field of inductive learning. First we show that there exist sets of examples that have no correct program, i.e., no program which implies all positive and no negative examples. Hence it only makes sense to talk about specialization problems for which a solution (a correct program) exists. To solve such problems, we first introduce UD1-specialization, based upon the transformation rule unfolding. We show UD1-specialization is incomplete - some solvable specialization problems do not have a UD1-specialization as solution - and generalize it to the stronger UD2-specialization. UD2 also turns out to be incomplete. An analysis of program specialization, using the subsumption theorem for SLD-resolution, shows the reason for this incompleteness. Based on that analysis, we then define UDS-specialization (a generalization of UD2-specialization), and prove that any specialization problem has a UDS-specialization as a solution. We also discuss the relationship between this specialization technique, and the generalization technique based on inverse resolution. Finally, we go into several more implementational matters, which outline an interesting topic for future research

Inducing ferromagnetism and Kondo effect in platinum by paramagnetic ionic gating

Electrically controllable magnetism, which requires the field-effect manipulation of both charge and spin degrees of freedom, has attracted growing interests since the emergence of spintronics. In this work, we report the reversible electrical switching of ferromagnetic (FM) states in platinum (Pt) thin films by introducing paramagnetic ionic liquid (PIL) as the gating media. The paramagnetic ionic gating controls the movement of ions with magnetic moments, which induces itinerant ferromagnetism on the surface of Pt films with large coercivity and perpendicular anisotropy mimicking the ideal two-dimensional Ising-type FM state. The electrical transport of the induced FM state shows Kondo effect at low temperature suggesting spatially separated coexistence of Kondo scattering beneath the FM interface. The tunable FM state indicates that paramagnetic ionic gating could serve as a versatile method to induce rich transport phenomena combining field effect and magnetism at PIL-gated interfaces.Comment: 17 pages, 4 figure