Computational Physics (physics.comp-ph)

  • PDF
    While the equilibrium properties, states, and phase transitions of interacting systems are well described by statistical mechanics, the lack of suitable state parameters has hindered the understanding of non-equilibrium phenomena in divers settings, from glasses to driven systems to biology. Here we introduce a simple idea enabling the quantification of organization in non-equilibrium and equilibrium systems, even when the form of order is unknown. The length of a losslessly compressed data file is a direct measure of its information content. Here we use data compression to study several out-of-equilibrium systems, and show that it both identifies ordering and reveals critical behavior in dynamical phase transitions. Our technique should provide a quantitative measure of organization in systems ranging from condensed matter systems in and out of equilibrium, to cosmology, biology and possibly economic and social systems.
  • PDF
    We present a simple closed expression for determining the condition for Coherent Perfect Absorption derived through electromagnetic wave treatment and interface boundary conditions. Apart from providing physical insight this expression can be used to estimate the values of various parameters required for observation of coherent perfect absorption in a given medium characterized by a complex dielectric constant. The results of the theoretical expression are found to be in good agreement with those obtained through numerical simulations.
  • PDF
    In search of new prospects for thermoelectric materials, using ab-initio calculations and semi-classical Boltzmann theory, we have systematically investigated the electronic structure and transport properties of 18-valence electron count cobalt based half-Heusler alloys with prime focus on CoVSn, CoNbSn, CoTaSn, CoMoIn, and CoWIn. The effect of doping on transport properties has been studied under the rigid band approximation. The maximum power factor, S$^2\sigma$, for all systems is obtained on hole doping and is comparable to the existing thermoelectric material CoTiSb. The stability of all the systems is verified by phonon calculations. Based on our calculations, we suggest that CoVSn, CoNbSn, CoTaSn, CoMoIn and CoWIn could be potential candidates for high temperature thermoelectric materials.
  • PDF
    Three complementary methods have been implemented in the code Denovo that accelerate neutral particle transport calculations with methods that use leadership-class computers fully and effectively: a multigroup block (MG) Krylov solver, a Rayleigh Quotient Iteration (RQI) eigenvalue solver, and a multigrid in energy (MGE) preconditioner. The MG Krylov solver converges more quickly than Gauss Seidel and enables energy decomposition such that Denovo can scale to hundreds of thousands of cores. RQI should converge in fewer iterations than power iteration (PI) for large and challenging problems. RQI creates shifted systems that would not be tractable without the MG Krylov solver. It also creates ill-conditioned matrices. The MGE preconditioner reduces iteration count significantly when used with RQI and takes advantage of the new energy decomposition such that it can scale efficiently. Each individual method has been described before, but this is the first time they have been demonstrated to work together effectively. The combination of solvers enables the RQI eigenvalue solver to work better than the other available solvers for large reactors problems on leadership class machines. Using these methods together, RQI converged in fewer iterations and in less time than PI for a full pressurized water reactor core. These solvers also performed better than an Arnoldi eigenvalue solver for a reactor benchmark problem when energy decomposition is needed. The MG Krylov, MGE preconditioner, and RQI solver combination also scales well in energy. This solver set is a strong choice for very large and challenging problems.
  • PDF
    I present an accurate and efficient technique for numerical evaluation of singular 6-dimensional integrals over tetrahedon-product domains, with applications to calculation of Galerkin matrix elements for discretized volume-integral-equation (VIE) solvers using Schaubert-Wilton-Glisson (SWG) and other tetrahedral basis functions. My method extends the generalized Taylor-Duffy strategy---used to handle the singular \textittriangle-product integrals arising in discretized surface-integral-equation (SIE) formulations---to the tetrahedron-product case; it effects an exact transformation of a singular 6-dimensional integral to an nonsingular lower-dimensional integral that may be evaluated by simple numerical cubature The method is highly general and may---with the aid of automatic code generation facilitated by computer-algebra systems---be applied to a wide variety of singular integrals arising in various VIE formulations with various types of tetrahedral basis function, of which I present several examples. To demonstrate the accuracy and efficiency of my method, I apply it to the calculation of matrix elements for the volume electric-field integral equation (VEFIE) discretized with SWG basis functions, where the method yields 12-digit or higher accuracy with low computational cost---an improvement of many orders of magnitude compared to existing techniques.

Recent comments

Chris Ferrie May 01 2017 12:52 UTC

I've posted a public referee report on this paper here: https://nbviewer.jupyter.org/github/csferrie/openreviews/blob/master/arxiv.1703.10743/arxiv.1703.10743.ipynb

Salvatore Mandrà Nov 21 2014 13:50 UTC

The manuscript has been widely revised to focus the reader's attention on the proposed method and its application in presence of local disorder.

Best,

Salvatore, Gian Giacomo and Alán

Salvatore Mandrà Aug 01 2014 19:11 UTC

Thanks Dr. Hastings for your comment. It is true that the transverse field Ising model does not satisfy the requirements to apply our method with an exponential reduction. Indeed, the opposite would be quite impressive since the random Ising model is a NP-Hard problem and we ourselves would be prett

...(continued)
Matt Hastings Aug 01 2014 16:43 UTC

The "quite general" conditions seem not to include the transverse field Ising model, the subject of most of the intensive numerical work previously. Incidentally, the terms "Lanczos" and "Krylov subspace" might be helpful.

Steve Flammia Jan 10 2014 08:01 UTC

This paper is an important first step towards making numerical and computational studies more easily checkable and reproducible. Very interesting stuff.

Noon van der Silk Jan 10 2014 16:26 UTC

Interestingly, another thematically-similar paper popped up today - https://scirate3.herokuapp.com/1401.2134 (posting the link here as I assume not many people are subbed to cs.DL)