Computational Physics (physics.comp-ph)

  • PDF
    Layered Li(Ni,Mn,Co,)O$_2$ (NMC) presents an intriguing ternary alloy design space for the optimization of performance as a cathode material in Li-ion batteries. Here, we present a high fidelity computational search of the ternary phase diagram with an emphasis on high-Ni containing compositional phases. This is done through the use of density functional theory training data fed into a reduced order model Hamiltonian that accounts for effective electronic and spin interactions of neighboring transition metal atoms at various lengths in a background of fixed lithium and oxygen atoms. This model can then be solved to include finite temperature thermodynamics into a convex hull analysis. We also provide a method to propagate the uncertainty at every level of the analysis to the final prediction of thermodynamically favorable compositional phases thus providing a quantitative measure of confidence for each prediction made. Due to the complexity of the three component system, as well as the intrinsic error of density functional theory, we argue that this propagation of uncertainty, particularly the uncertainty due to exchange-correlation functional choice is necessary to have reliable and interpretable results. With our final result, we recover the prediction of already known phases such as LiNi$_{0.33}$Mn$_{0.33}$Co$_{0.33}$O$_2$ (111) and LiNi$_{0.8}$Mn$_{0.1}$Co$_{0.1}$O$_2$ (811) in exact proportion while finding other proportions very close to the experimentally claimed LiNi$_{0.6}$Mn$_{0.2}$Co$_{0.2}$O$_2$ (622) and LiNi$_{0.5}$Mn$_{0.3}$Co$_{0.2}$O$_2$ (532) phases, and overall predict a total of 37 phases with reasonable confidence and 69 more phases with a lower level of confidence. Through our analysis, we also can identify the phases with the highest average operational voltage at a given Co composition.
  • PDF
    A quantum interatomic scattering is implemented in the direct simulation Monte Carlo (DSMC) method applied to transport phenomena in rarefied gases. In contrast to the traditional DSMC method based on the classical scattering, the proposed implementation allows us to model flows of gases over the whole temperature range beginning from 1 K up any high temperature when no ionization happens. To illustrate the new numerical approach, two helium isotopes $^3$He and $^4$He were considered in two canonical problems, namely, heat transfer between two planar surfaces and planar Couette flow. To solve these problems, the ab initio potential for helium is used, but the proposed technique can be used with any intermolecular potential. The problems were solved over the temperature range from 1 K to 3000 K and for two values of the rarefaction parameter 1 and 10. The former corresponds to the transitional regime and the last describes the temperature jump and velocity slip regime. No influence of the quantum effects was detected within the numerical error of 0.1 % for the temperature 300 K and higher. However, the quantum approach requires less computational effort than the classical one in this temperature range. For temperatures lower than 300 K, the influence of the quantum effects exceed the numerical error and reaches 67 % at the temperature of 1 K.
  • PDF
    In this paper we numerically analyze the 1D self-localized solutions of the Kundu-Eckhaus equation (KEE) in nonlinear waveguides using the spectral renormalization method (SRM) and compare our findings with those solutions of the nonlinear Schrodinger equation (NLSE). We show that single, dual and N-soliton solutions exist for the case with zero optical potentials, i.e. V=0. We also show that these soliton solutions do not exist, at least for a range of parameters, for the photorefractive lattices with optical potentials in the form of V=Io.cos^2(x) for cubic nonlinearity. However, self-stable solutions of the KEE with saturable nonlinearity do exist for some range of parameters. We compare our findings for the KEE with those of the NLSE and discuss our results.
  • PDF
    In this paper we present scaling results of a FFT library, FFTK, and a pseudospectral code, Tarang, on grid resolutions up to $8192^3$ grid using 65536 cores of Blue Gene/P and 196608 cores of Cray XC40 supercomputers. We observe that communication dominates computation, more so on the Cray XC40. The computation time scales as $T_\mathrm{comp} \sim p^{-1}$, and the communication time as $T_\mathrm{comm} \sim n^{-\gamma_2}$ with $\gamma_2$ ranging from 0.7 to 0.9 for Blue Gene/P, and from 0.43 to 0.73 for Cray XC40. FFTK, and the fluid and convection solvers of Tarang exhibit weak as well as strong scaling nearly up to 196608 cores of Cray XC40. We perform a comparative study of the performance on the Blue Gene/P and Cray XC40 clusters.
  • PDF
    Without setting any initial Bravais lattice, both the crystallization and non-crystallization of Lennard-Jones (LJ) particles have been investigated by molecular dynamics (MD) simulation with a simple LJ potential. To identify the Bravais lattice of simulated system, the distribution functions of both the angles between one particle and its nearest neighbors and the distances between the particles have been introduced. The final identification can be made by comparing these two calculated distribution functions with those of ideal Bravais lattices and checking the particle arrangement of simulated system. The simulated systems are found to show the face-centered cubic (fcc) and ideal hexagonal close-packed (hcp) lattices. The non-crystallization system also has been obtained, and its microstructure is found to be similar to that of LJ liquid at a temperature near the crystallization temperature, and show no order of the second nearest neighbors in comparison with that of crystallization system. This paper has proposed a new way of investigating the microstructure of material and its evolution, and paved the way for MD simulation of large scale particle system consisting of more than one million particles.
  • PDF
    The efficient method for computing the sensitivities is the adjoint method. The cost of solving an adjoint equation is comparable to the cost of solving the governing equation. Once the adjoint solution is obtained, the sensitivities to any number of parameters can be obtained with little effort. There are two methods to develop the adjoint equations: continuous method and discrete method. In the continuous method, the control theory is applied to the forward governing equation and produces an analytical partial differential equation for solving the adjoint variable; in the discrete method, the control theory is applied to the discrete form of the forward governing equation and produces a linear system of equations for solving the adjoint variable. In this article, an adjoint sensitivity analysis framework is developed using both the continuous and discrete methods. These two methods are assessed with the faucet flow for steady-state problem and one transient test case based on the BFBT benchmark for transient problem. Adjoint sensitivities from both methods are verified by sensitivities given by the perturbation method. Adjoint sensitivities from both methods are physically reasonable and match each. The sensitivities obtained with discrete method is found to be more accurate than the sensitivities from the continuous method. The continuous method is computationally more efficient than the discrete method because of the analytical coefficient matrices and vectors. However, difficulties are observed in solving the continuous adjoint equation for cases where the adjoint equation contains sharp discontinuities in the source terms; in such cases, the continuous method is not as robust as the discrete adjoint method.
  • PDF
    Contact dynamics (CD) is a powerful method to solve the dynamics of large systems of colliding rigid bodies. CD can be computationally more efficient than classical penalty-based discrete element methods (DEM) for simulating contact between stiff materials such as rock, glass, or engineering metals. However, by idealizing bodies as perfectly rigid, contact forces computed by CD can be non-unique due to indeterminacy in the contact network, which is a common occurence in dense granular flows. We propose a CD method that is designed to identify only the unique set of contact forces that would be predicted by a soft particle method, such as DEM, in the limit of large stiffness. The method involves applying an elastic compatibility condition to the contact forces, which maintains no-penetration constraints but filters out force distributions that could not have arisen from stiff elastic contacts. The method can be used as a post-processing step that could be integrated into existing CD codes with minimal effort. We demonstrate its efficacy in a variety of indeterminate problems, including some involving multiple materials, non-spherical shapes, and nonlinear contact constitutive laws.

Recent comments

Dan Elton Mar 16 2018 04:36 UTC

Comments are appreciated. Message me here or on twitter @moreisdifferent

Code is open source and available at :


Chris Ferrie May 01 2017 12:52 UTC

I've posted a public referee report on this paper here:

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.


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

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 - (posting the link here as I assume not many people are subbed to cs.DL)