Showing papers in "arXiv: Computational Physics in 2014"

The design and verification of Mumax3

Abstract: We report on the design, verification and performance of mumax3, an open-source GPU-accelerated micromagnetic simulation program. This software solves the time- and space dependent magnetization evolution in nano- to micro scale magnets using a finite-difference discretization. Its high performance and low memory requirements allow for large-scale simulations to be performed in limited time and on inexpensive hardware. We verified each part of the software by comparing results to analytical values where available and to micromagnetic standard problems. mumax3 also offers specific extensions like MFM image generation, moving simulation window, edge charge removal and material grains.

LanHEP - a package for automatic generation of Feynman rules from the Lagrangian. Updated version 3.2

Abstract: We present a new version 3.2 of the LanHEP software package. New features include UFO output, color sextet particles and new substutution techniques which allow to define new routines.

Evaluation of Constant Potential Method in Simulating Electric Double-Layer Capacitors

Abstract: A major challenge in the molecular simulation of electric double layer capacitors (EDLCs) is the choice of an appropriate model for the electrode. Typically, in such simulations the electrode surface is modeled using a uniform fixed charge on each of the electrode atoms, which ignores the electrode response to local charge fluctuations induced by charge fluctuations in the electrolyte. In this work, we evaluate and compare this Fixed Charge Method (FCM) with the more realistic Constant Potential Method (CPM), [Reed, et al., J. Chem. Phys., 126, 084704 (2007)], in which the electrode charges fluctuate in order to maintain constant electric potential in each electrode. For this comparison, we utilize a simplified LiClO$_4$-acetonitrile/graphite EDLC. At low potential difference ($\Delta\Psi\le 2V$), the two methods yield essentially identical results for ion and solvent density profiles; however, significant differences appear at higher $\Delta\Psi$. At $\Delta\Psi\ge 4V$, the CPM ion density profiles show significant enhancement (over FCM) of "partially electrode solvated" Li$^+$ ions very close to the electrode surface. The ability of the CPM electrode to respond to local charge fluctuations in the electrolyte is seen to significantly lower the energy (and barrier) for the approach of Li$^+$ ions to the electrode surface.

Manifestation of unexpected semiconducting properties in few-layer orthorhombic arsenene

Abstract: In this express, we demonstrate few-layer orthorhombic arsenene is an ideal semiconductor. Due to the layer stacking, multilayer arsenenes always behave as intrinsic direct bandgap semiconductors with gap values of around 1 eV. In addition, these bandgaps can be further tuned in its nanoribbons. Based on the so-called acoustic phonon limited approach, the carrier mobilities are predicted to approach as high as several thousand square centimeters per volt-second and simultaneously exhibit high directional anisotropy. All these make few-layer arsenene promising for device applications in semiconducting industry.

Tiling Phosphorene

Abstract: We present a scheme to categorize the structure of different layered phosphorene allotropes by mapping their non-planar atomic structure onto a two-color 2D triangular tiling pattern. In the buckled structure of a phosphorene monolayer, we assign atoms in "top" positions to dark tiles and atoms in "bottom" positions to light tiles. Optimum $sp^3$ bonding is maintained throughout the structure when each triangular tile is surrounded by the same number $N$ of like-colored tiles, with $0{\le}N{\le}2$. Our ab initio density functional calculations indicate that both the relative stability and electronic properties depend primarily on the structural index $N$. The proposed mapping approach may also be applied to phosphorene structures with non-hexagonal rings and 2D quasicrystals with no translational symmetry, which we predict to be nearly as stable as the hexagonal network.

Diffusion-Based Coarse Graining in Hybrid Continuum-Discrete Solvers: Theoretical Formulation and A Priori Tests

Abstract: Coarse graining is an important ingredient in many multi-scale continuum-discrete solvers such as CFD--DEM (computational fluid dynamics--discrete element method) solvers for dense particle-laden flows. Although CFD--DEM solvers have become a mature technique that is widely used in multiphase flow research and industrial flow simulations, a flexible and easy-to-implement coarse graining algorithm that can work with CFD solvers of arbitrary meshes is still lacking. In this work, we proposed a new coarse graining algorithm for continuum--discrete solvers for dense particle-laden flows based on solving a transient diffusion equation. Via theoretical analysis we demonstrated that the proposed method is equivalent to the statistical kernel method with a Gaussian kernel, but the current method is much more straightforward to implement in CFD--DEM solvers. \textit{A priori} numerical tests were performed to obtain the solid volume fraction fields based on given particle distributions, the results obtained by using the proposed algorithm were compared with those from other coarse graining methods in the literature (e.g., the particle centroid method, the divided particle volume method, and the two-grid formulation). The numerical tests demonstrated that the proposed coarse graining procedure based on solving diffusion equations is theoretically sound, easy to implement and parallelize in general CFD solvers, and has improved mesh-convergence characteristics compared with existing coarse graining methods. The diffusion-based coarse graining method has been implemented into a CFD--DEM solver, the results of which are presented in a separate work (R. Sun and H. Xiao, Diffusion-based coarse graining in hybrid continuum-discrete solvers: Application in CFD-DEM solvers for particle laden flows).

Anomalous doping effect in black phosphorene from first-principles calculations

Abstract: Using first-principles density functional theory calculations, we investigate the geometries, electronic structures, and thermodynamic stabilities of substitutionally doped phosphorene sheets with group III, IV, V, and VI elements. We find that the electronic properties of phosphorene are drastically modified by the number of valence electrons in dopant atoms. The dopants with even number of valence electrons enable the doped phosphorenes to have a metallic feature, while the dopants with odd number of valence electrons keep the semiconducting feature with a larger band gap than the undoped phosphorene. This even-odd behavior is attributed to the peculiar bonding characteristics of phosphorene and the strong hybridization of sp orbitals between dopants and phosphorene. The calculated formation energies of various substitutional dopants in phosphorene show that such doped systems can be thermodynamically stable.

Unbiased Reduced Density Matrices and Electronic Properties from Full Configuration Interaction Quantum Monte Carlo

Abstract: Properties that are necessarily formulated within pure (symmetric) expectation values are difficult to calculate for projector quantum Monte Carlo approaches, but are critical in order to compute many of the important observable properties of electronic systems. Here, we investigate an approach for the sampling of unbiased reduced density matrices within the Full Configuration Interaction Quantum Monte Carlo dynamic, which requires only small computational overheads. This is achieved via an independent replica population of walkers in the dynamic, sampled alongside the original population. The resulting reduced density matrices are free from systematic error (beyond those present via constraints on the dynamic itself), and can be used to compute a variety of expectation values and properties, with rapid convergence to an exact limit. A quasi-variational energy estimate derived from these density matrices is proposed as an accurate alternative to the projected estimator for multiconfigurational wavefunctions, while its variational property could potentially lend itself to accurate extrapolation approaches in larger systems.

A comparison between parallelization approaches in molecular dynamics simulations on GPUs

Abstract: We test the performances of two different approaches to the computation of forces for molecular dynamics simulations on Graphics Processing Units. A "vertex-based" approach, where a computing thread is started per particle, is compared to a newly proposed "edge-based" approach, where a thread is started per each potentially non-zero interaction. We find that the former is more efficient for systems with many simple interactions per particle, while the latter is more efficient if the system has more complicated interactions or fewer of them. By comparing computation times on more and less recent GPU technology, we predict that, if the current trend of increasing the number of processing cores - as opposed to their computing power - remains, the "edge-based" approach will gradually become the most efficient choice in an increasing number of cases.

Calculating vibrational spectra with sum of product basis functions without storing full-dimensional vectors or matrices

Abstract: We propose an iterative method for computing vibrational spectra that significantly reduces the memory cost of calculations. It uses a direct product primitive basis, but does not require storing vectors with as many components as there are product basis functions. Wavefunctions are represented in a basis each of whose functions is a sum of products (SOP) and the factorizable structure of the Hamiltonian is exploited. If the factors of the SOP basis functions are properly chosen, wavefunctions are linear combinations of a small number of SOP basis functions. The SOP basis functions are generated using a shifted block power method. The factors are refined with a rank reduction algorithm to cap the number of terms in a SOP basis function. The ideas are tested on a 20-D model Hamiltonian and a realistic CH$_3$CN (12 dimensional) potential. For the 20-D problem, to use a standard direct product iterative approach one would need to store vectors with about $10^{20}$ components and would hence require about $8 \times 10^{11}$ GB. With the approach of this paper only 1 GB of memory is necessary. Results for CH$_3$CN agree well with those of a previous calculation on the same potential.

Modeling radio communication blackout and blackout mitigation in hypersonic vehicles

Abstract: A procedure for the modeling and analysis of radio communication blackout of hypersonic vehicles is presented. The weakly ionized plasma generated around the surface of a hypersonic reentry vehicle is simulated using full Navier-Stokes equations in multi-species single fluid form. A seven species air chemistry model is used to compute the individual species densities in air including ionization - plasma densities are compared with experiment. The electromagnetic wave's interaction with the plasma layer is modeled using multi-fluid equations for fluid transport and full Maxwell's equations for the electromagnetic fields. The multi-fluid solver is verified for a whistler wave propagating through a slab. First principles radio communication blackout over a hypersonic vehicle is demonstrated along with a simple blackout mitigation scheme using a magnetic window.

Crystalline 'Genes' in Metallic Liquids

Abstract: The underlying structural order that transcends the liquid, glass and crystalline states is identified using an efficient genetic algorithm (GA). GA identifies the most common energetically favorable packing motif in crystalline structures close to the alloy's Al-10 at.% Sm composition. These motifs are in turn compared to the observed packing motifs in the actual liquid structures using a cluster-alignment method which reveals the average topology. Conventional descriptions of the short-range order, such as Voronoi tessellation, are too rigid in their analysis of the configurational poly-types when describing the chemical and topological ordering during transition from undercooled metallic liquids to crystalline phases or glass. Our approach here brings new insight into describing mesoscopic order-disorder transitions in condensed matter physics.

Calculation of Energy Band Diagram of a Photoelectrochemical Water Splitting Cell

Abstract: A physical model is presented for a semiconductor electrode of a photoelectrochemical (PEC) cell, accounting for the potential drop in the Helmholtz layer. Hence both band edge pinning and unpinning are naturally included in our description. The model is based on the continuity equations for charge carriers and direct charge transfer from the energy bands to the electrolyte. A quantitative calculation of the position of the energy bands and the variation of the quasi-Fermi levels in the semiconductor with respect to the water reduction and oxidation potentials is presented. Calculated current-voltage curves are compared with established analytical models and measurement. Our model calculations are suitable to enhance understanding and improve properties of semiconductors for photoelectrochemical water splitting.

Concurrent Cuba

Abstract: The parallel version of the multidimensional numerical integration package Cuba is presented and achievable speed-ups discussed

Modeling electron dynamics coupled to continuum states in finite volumes

Abstract: Absorbing boundaries are frequently employed in real-time propagation of the Schr\"odinger equation to remove spurious reflections and efficiently emulate outgoing boundary conditions. These conditions are a fundamental ingredient for an implicit description of observables involving infinitely extended continuum states. In the literature, several boundary absorbers have been proposed. They mostly fall into three main families: mask function absorbers, complex absorbing potentials, and exterior complex-scaled potentials. To date none of the proposed absorbers is perfect, and all present a certain degree of reflections. Characterization of such reflections is thus a critical task with strong implications for time-dependent simulations of atoms and molecules. We introduce a method to evaluate the reflection properties of a given absorber and present a comparison of selected samples for each family of absorbers. Further, we discuss the connections between members of each family and show how the same reflection curves can be obtained with very different absorption schemes.

Compressible Generalized Hybrid Monte Carlo

Abstract: One of the most demanding calculations is to generate random samples from a specified probability distribution (usually with an unknown normalizing prefactor) in a high-dimensional configuration space. One often has to resort to using a Markov chain Monte Carlo method, which converges only in the limit to the prescribed distribution. Such methods typically inch through configuration space step by step, with acceptance of a step based on a Metropolis(-Hastings) criterion. An acceptance rate of 100% is possible in principle by embedding configuration space in a higher-dimensional phase space and using ordinary differential equations. In practice, numerical integrators must be used, lowering the acceptance rate. This is the essence of hybrid Monte Carlo methods. Presented is a general framework for constructing such methods under relaxed conditions: the only geometric property needed is (weakened) reversibility; volume preservation is not needed. The possibilities are illustrated by deriving a couple of explicit hybrid Monte Carlo methods, one based on barrier-lowering variable-metric dynamics and another based on isokinetic dynamics.

Photoionization of few electron systems with a hybrid Coupled Channels approach

Abstract: We present the hybrid anti-symmetrized coupled channels method for the calculation of fully differential photo-electron spectra of multi-electron atoms and small molecules interacting with strong laser fields. The method unites quantum chemical few-body electronic structure with strong-field dynamics by solving the time dependent Schrodinger equation in a fully anti-symmetrized basis composed of multi-electron states from quantum chemistry and a one-electron numerical basis. Photoelectron spectra are obtained via the time dependent surface flux (tSURFF) method. Performance and accuracy of the approach are demonstrated for spectra from the helium and berryllium atoms and the hydrogen molecule in linearly polarized laser fields at wavelength from 21 nm to 400 nm. At long wavelengths, helium and the hydrogen molecule at equilibrium inter-nuclear distance can be approximated as single channel systems whereas beryllium needs a multi-channel description.

An Auxiliary-Field Quantum Monte Carlo Study of the Chromium Dimer

Abstract: The chromium dimer (Cr2) presents an outstanding challenge for many-body electronic structure methods. Its complicated nature of binding, with a formal sextuple bond and an unusual potential energy curve, is emblematic of the competing tendencies and delicate balance found in many strongly correlated materials. We present a near-exact calculation of the potential energy curve (PEC) and ground state properties of Cr2, using the auxiliary-field quantum Monte Carlo (AFQMC) method. Unconstrained, exact AFQMC calculations are first carried out for a medium-sized but realistic basis set. Elimination of the remaining finite-basis errors and extrapolation to the complete basis set (CBS) limit is then achieved with a combination of phaseless and exact AFQMC calculations. Final results for the PEC and spectroscopic constants are in excellent agreement with experiment.

From Molecular Dynamics to Brownian Dynamics

Abstract: Three coarse-grained molecular dynamics (MD) models are investigated with the aim of developing and analyzing multiscale methods which use MD simulations in parts of the computational domain and (less detailed) Brownian dynamics (BD) simulations in the remainder of the domain. The first MD model is formulated in one spatial dimension. It is based on elastic collisions of heavy molecules (e.g. proteins) with light point particles (e.g. water molecules). Two three-dimensional MD models are then investigated. The obtained results are applied to a simplified model of protein binding to receptors on the cellular membrane. It is shown that modern BD simulators of intracellular processes can be used in the bulk and accurately coupled with a (more detailed) MD model of protein binding which is used close to the membrane.

SIESTA-PEXSI: Massively parallel method for efficient and accurate \textit{ab initio} materials simulation without matrix diagonalization

Abstract: We describe a scheme for efficient large-scale electronic-structure calculations based on the combination of the pole expansion and selected inversion (PEXSI) technique with the SIESTA method, which uses numerical atomic orbitals within the Kohn-Sham density functional theory (KSDFT) framework The PEXSI technique can efficiently utilize the sparsity pattern of the Hamiltonian and overlap matrices generated in SIESTA, and for large systems has a much lower computational complexity than that associated with the matrix diagonalization procedure The PEXSI technique can be used to evaluate the electron density, free energy, atomic forces, density of states and local density of states without computing any eigenvalue or eigenvector of the Kohn-Sham Hamiltonian It can achieve accuracy fully comparable to that obtained from a matrix diagonalization procedure for general systems, including metallic systems at low temperature The PEXSI method is also highly scalable With the recently developed massively parallel PEXSI technique, we can make efficient use of more than $10,000$ processors on high performance machines We demonstrate the performance and accuracy of the SIESTA-PEXSI method using several examples of large scale electronic structure calculations, including 1D, 2D and bulk problems with insulating, semi-metallic, and metallic character

Alya: Towards Exascale for Engineering Simulation Codes

Abstract: Alya is the BSC in-house HPC-based multi-physics simulation code. It is designed from scratch to run efficiently in parallel supercomputers, solving coupled problems. The target domain is engineering, with all its particular features: complex geome- tries and unstructured meshes, coupled multi-physics with exotic coupling schemes and Physical models, ill-posed problems, flexibility needs for rapidly including new models, etc. Since its conception in 2004, Alya has shown scaling behaviour in an increasing number of cores. In this paper, we present its performance up to 100.000 cores in Blue Waters, the NCSA supercomputer. The selected tests are representative of the engineering world, all the problematic features included: incompressible flow in a hu- man respiratory system, low Mach combustion problem in a kiln furnace and coupled electro-mechanical problem in a heart. We show scalability plots for all cases, discussing all the aspects of such kind of simulations, including solvers convergence.

PIBM: Particulate immersed boundary method for fluid-particle interaction problems

Abstract: It is well known that the number of particles should be scaled up to enable industrial scale simulation. The calculations are more computationally intensive when the motion of the surrounding fluid is considered. Besides the advances in computer hardware and numerical algorithms, the coupling scheme also plays an important role on the computational efficiency. In this study, a particle immersed boundary method (PIBM) for simulating the fluid-particle multiphase flow was presented and assessed in both two- and three-dimensional applications. The idea behind PIBM derives from the conventional momentum exchange-based immersed boundary method (IBM) by treating each Lagrangian point as a solid particle. This treatment enables LBM to be coupled with fine particles residing within a particular grid cell. Compared with the conventional IBM, dozens of times speedup in two-dimensional simulation and hundreds of times in three-dimensional simulation can be expected under the same particle and mesh number. Numerical simulations of particle sedimentation in the Newtonian flows were conducted based on a combined lattice Boltzmann method - particle immersed boundary method - discrete element method scheme, showing that the PIBM can capture the feature of particulate flows in fluid and is indeed a promising scheme for the solution of the fluid-particle interaction problems.

Molecular Dynamics in a Grand Ensemble: Bergmann-Lebowitz model and Adaptive Resolution Simulation

Abstract: This article deals with the molecular dynamics simulation of open systems that can exchange energy and matter with a reservoir; the physics of the reservoir and its interactions with the system are described by the model introduced by Bergmann and Lebowitz.Despite its conceptual appeal, the model did not gain popularity in the field of molecular simulation and, as a consequence, did not play a role in the development of open system molecular simulation techniques, even though it can provide the conceptual legitimation of simulation techniques that mimic open systems. We shall demonstrate that the model can serve as a tool to devise both numerical procedures and conceptual definitions of physical quantities that cannot be defined in a straightforward way by systems with a fixed number of molecules. In particular, we discuss the utility of the Bergmann-Lebowitz (BL) model for the calculation of equilibrium time correlation functions within the Grand Canonical Adaptive Resolution method (GC-AdResS) and report numerical results for the case of liquid water.

A three-dimensional spin-diffusion model for micromagnetics

Abstract: We implement a finite-element scheme that solves the Landau-Lifshitz-Gilbert equation coupled to a diffusion equation accounting for spin-polarized currents. The latter solves for the spin accumulation not only in magnetic materials but also in nonmagnetic conductors. The presented method incorporates the model by Slonczewski for the description of spin torque in magnetic multilayers as well as the model of Zhang and Li for the description of current driven domain-wall motion. Furthermore it is able to do both resolve the time evolution of the spin accumulation or treat it in an adiabatic fashion by the choice of sufficiently large time steps.

xTRAM: Estimating equilibrium expectations from time-correlated simulation data at multiple thermodynamic states

Abstract: Computing the equilibrium properties of complex systems, such as free energy differences, is often hampered by rare events in the dynamics. Enhanced sampling methods may be used in order to speed up sampling by, for example, using high temperatures, as in parallel tempering, or simulating with a biasing potential such as in the case of umbrella sampling. The equilibrium properties of the thermodynamic state of interest (e.g., lowest temperature or unbiased potential) can be computed using reweighting estimators such as the weighted histogram analysis method or the multistate Bennett acceptance ratio (MBAR). weighted histogram analysis method and MBAR produce unbiased estimates, the simulation samples from the global equilibria at their respective thermodynamic state--a requirement that can be prohibitively expensive for some simulations such as a large parallel tempering ensemble of an explicitly solvated biomolecule. Here, we introduce the transition-based reweighting analysis method (TRAM)--a class of estimators that exploit ideas from Markov modeling and only require the simulation data to be in local equilibrium within subsets of the configuration space. We formulate the expanded TRAM (xTRAM) estimator that is shown to be asymptotically unbiased and a generalization of MBAR. Using four exemplary systems of varying complexity, we demonstrate the improved convergence (ranging from a twofold improvement to several orders of magnitude) of xTRAM in comparison to a direct counting estimator and MBAR, with respect to the invested simulation effort. Lastly, we introduce a random-swapping simulation protocol that can be used with xTRAM, gaining orders-of-magnitude advantages over simulation protocols that require the constraint of sampling from a global equilibrium.

NCON: A tensor network contractor for MATLAB

Abstract: A fundamental process in the implementation of any numerical tensor network algorithm is that of contracting a tensor network. In this process, a network made up of multiple tensors connected by summed indices is reduced to a single tensor or a number by evaluating the index sums. This article presents a MATLAB function ncon(), or "Network CONtractor", which accepts as its input a tensor network and a contraction sequence describing how this network may be reduced to a single tensor or number. As its output it returns that single tensor or number. The function ncon() may be obtained by downloading the source of this preprint.

Higher-order finite-difference formulation of periodic Orbital-free Density Functional Theory

Abstract: We present a real-space formulation and higher-order finite-difference implementation of periodic Orbital-free Density Functional Theory (OF-DFT). Specifically, utilizing a local reformulation of the electrostatic and kernel terms, we develop a generalized framework for performing OF-DFT simulations with different variants of the electronic kinetic energy. In particular, we propose a self-consistent field (SCF) type fixed-point method for calculations involving linear-response kinetic energy functionals. In this framework, evaluation of both the electronic ground-state as well as forces on the nuclei are amenable to computations that scale linearly with the number of atoms. We develop a parallel implementation of this formulation using the finite-difference discretization. We demonstrate that higher-order finite-differences can achieve relatively large convergence rates with respect to mesh-size in both the energies and forces. Additionally, we establish that the fixed-point iteration converges rapidly, and that it can be further accelerated using extrapolation techniques like Anderson's mixing. We validate the accuracy of the results by comparing the energies and forces with plane-wave methods for selected examples, including the vacancy formation energy in Aluminum. Overall, the suitability of the proposed formulation for scalable high performance computing makes it an attractive choice for large-scale OF-DFT calculations consisting of thousands of atoms.

On the Analysis of the Discretized Kohn-Sham Density Functional Theory

Abstract: In this paper, we study a few theoretical issues in the discretized Kohn-Sham (KS) density functional theory (DFT). The equivalence between either a local or global minimizer of the KS total energy minimization problem and the solution to the KS equation is established under certain assumptions. The nonzero charge densities of a strong local minimizer are shown to be bounded below by a positive constant uniformly. We analyze the self-consistent field (SCF) iteration by formulating the KS equation as a fixed point map with respect to the potential. The Jacobian of these fixed point maps is derived explicitly. Both global and local convergence of the simple mixing scheme can be established if the gap between the occupied states and unoccupied states is sufficiently large. This assumption can be relaxed if the charge density is computed using the Fermi-Dirac distribution and it is not required if there is no exchange correlation functional in the total energy functional. Although our assumption on the gap is very stringent and is almost never satisfied in reality, our analysis is still valuable for a better understanding of the KS minimization problem, the KS equation and the SCF iteration.

Ultra--fast carriers relaxation in bulk silicon following photo--excitation with a short and polarized laser pulse

Abstract: A novel approach based on the merging of the out--of--equilibrium Green's function method with the ab-initio, Density--Functional--Theory is used to describe the ultra--fast carriers relaxation in Silicon. The results are compared with recent two photon photo--emission measurements. We show that the interpretation of the carrier relaxation in terms of L -> X inter--valley scattering is not correct. The ultra--fast dynamics measured experimentally is, instead, due to the scattering between degenerate $L$ states that is activated by the non symmetric population of the conduction bands induced by the laser field. This ultra--fast relaxation is, then, entirely due to the specific experimental setup and it can be interpreted by introducing a novel definition of the quasi--particle lifetimes in an out--of--equilibrium context.

Continuous surface force based lattice Boltzmann equation method for simulating thermocapillary flow

Abstract: In this paper, we extend a lattice Boltzmann equation (LBE) with continuous surface fore (CSF) to simulate thermocapillary flows. The model is designed on our previous CSF LBE for athermal two phase flow, in which the interfacial tension forces and the Marangoni stresses as the results of the interface interactions between different phases are described by a conception of CSF. In this model, the sharp interfaces between different phases are separated by a narrow transition layers, and the kinetics and morphology evolution of phase separation would be characterized by an order parameter visa Cahn-Hilliard equation which is solved in the frame work of LBE. The scalar convection-diffusion equation for temperature field is also solved by thermal LBE. The models are validated by thermal two layered Poiseuille flow, and a two superimposed planar fluids at negligibly small Reynolds and Marangoni numbers for the thermocapillary driven convection, which have analytical solutions for the velocity and temperature. Then thermocapillary migration of two dimensional deformable droplet are simulated. Numerical results show that the predictions of present LBE agreed with the analytical solution/other numerical results.