Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userInclusion and validation of electronic stopping in the open source LAMMPS code
58
0
0.0
(
0
)
Added by
Manoj Kumar Warrier Dr.
Publication date
2020
fields
Physics
and research's language is
English
Ask ChatGPT about the research
No Arabic abstract
Electronic stopping (ES) of energetic atoms is not taken care of by the interatomic potentials used in molecular dynamics (MD) simulations when simulating collision cascades. The Lindhard-Scharff (LS) formula for electronic stopping is therefore included as a drag term for energetic atoms in the open source large scale atomic molecular massively parallel simulator (LAMMPS) code. In order to validate the ES implementation, MD simulations of collision cascades at primary knock-on atom (PKA) energies of 5, 10 and 20 keV are carried out in W and Fe in 100 random directions. The total ES losses from the MD simulations show energy straggling due to the stochastic nature of the phenomena. Thelosses due to ES are compared with that predicted by theNorgett-Robinson-Torrens (NRT) model to validate our implementation. It is seenthat the root mean square deviation of ES losses from the MD implementation is around 10 % for both W and Fe compared to the NRT model. The velocity threshold above which electronic stopping is important is explored. The effect of ES on the number of defects in collision cascades is presented for Fe and W.
rate research
Read More
QMCPACK is an open source quantum Monte Carlo package for ab-initio electronic structure calculations. It supports calculations of metallic and insulating solids, molecules, atoms, and some model Hamiltonians. Implemented real space quantum Monte Carlo algorithms include variational, diffusion, and reptation Monte Carlo. QMCPACK uses Slater-Jastrow type trial wave functions in conjunction with a sophisticated optimizer capable of optimizing tens of thousands of parameters. The orbital space auxiliary field quantum Monte Carlo method is also implemented, enabling cross validation between different highly accurate methods. The code is specifically optimized for calculations with large numbers of electrons on the latest high performance computing architectures, including multicore central processing unit (CPU) and graphical processing unit (GPU) systems. We detail the programs capabilities, outline its structure, and give examples of its use in current research calculations. The package is available at http://www.qmcpack.org .
The surface code is highly practical, enabling arbitrarily reliable quantum computation given a 2-D nearest-neighbor coupled array of qubits with gate error rates below approximately 1%. We describe an open source library, Polyestimate, enabling a user with no knowledge of the surface code to specify realistic physical quantum gate error models and obtain logical error rate estimates. Functions allowing the user to specify simple depolarizing error rates for each gate have also been included. Every effort has been made to make this library user-friendly.
We implement two recently developed fast Coulomb solvers, HSMA3D [J. Chem. Phys. 149 (8) (2018) 084111] and HSMA2D [J. Chem. Phys. 152 (13) (2020) 134109], into a new user package HSMA for molecular dynamics simulation engine LAMMPS. The HSMA package is designed for efficient and accurate modeling of electrostatic interactions in 3D and 2D periodic systems with dielectric effects at the O(N) cost. The implementation is hybrid MPI and OpenMP parallelized and compatible with existing LAMMPS functionalities. The vectorization technique following AVX512 instructions is adopted for acceleration. To establish the validity of our implementation, we have presented extensive comparisons to the widely used particle-particle particle-mesh (PPPM) algorithm in LAMMPS and other dielectric solvers. With the proper choice of algorithm parameters and parallelization setup, the package enables calculations of electrostatic interactions that outperform the standard PPPM in speed for a wide range of particle numbers.
Routine applications of electronic structure theory to molecules and periodic systems need to compute the electron density from given Hamiltonian and, in case of non-orthogonal basis sets, overlap matrices. System sizes can range from few to thousands or, in some examples, millions of atoms. Different discretization schemes (basis sets) and different system geometries (finite non-periodic vs. infinite periodic boundary conditions) yield matrices with different structures. The ELectronic Structure Infrastructure (ELSI) project provides an open-source software interface to facilitate the implementation and optimal use of high-performance solver libraries covering cubic scaling eigensolvers, linear scaling density-matrix-based algorithms, and other reduced scaling methods in between. In this paper, we present recent improvements and developments inside ELSI, mainly covering (1) new solvers connected to the interface, (2) matrix layout and communication adapted for parallel calculations of periodic and/or spin-polarized systems, (3) routines for density matrix extrapolation in geometry optimization and molecular dynamics calculations, and (4) general utilities such as parallel matrix I/O and JSON output. The ELSI interface has been integrated into four electronic structure code projects (DFTB+, DGDFT, FHI-aims, SIESTA), allowing us to rigorously benchmark the performance of the solvers on an equal footing. Based on results of a systematic set of large-scale benchmarks performed with Kohn-Sham density-functional theory and density-functional tight-binding theory, we identify factors that strongly affect the efficiency of the solvers, and propose a decision layer that assists with the solver selection process. Finally, we describe a reverse communication interface encoding matrix-free iterative solver strategies that are amenable, e.g., for use with planewave basis sets.
In classical analyses of $gamma$-ray data from IACTs, such as H.E.S.S., aperture photometry, or photon counting, is applied in a (typically circular) region of interest (RoI) encompassing the source. A key element in the analysis is to estimate the amount of background in the RoI due to residual cosmic ray-induced air showers in the data. Various standard background estimation techniques have been developed in the last decades, most of them rely on a measurement of the background from source-free regions within the observed field of view. However, in particular in the Galactic plane, source analysis and background estimation are hampered by the large number of, sometimes overlapping, $gamma$-ray sources and large-scale diffuse $gamma$-ray emission. For complicated fields of view, a three-dimensional (3D) likelihood analysis shows the potential to be superior to classical analysis. In this analysis technique, a spectromorphological model, consisting of one or multiple source components and a background component, is fitted to the data, resulting in a complete spectral and spatial description of the field of view. For the application to IACT data, the major challenge of such an approach is the construction of a robust background model. In this work, we apply the 3D likelihood analysis to various test data recently made public by H.E.S.S., using the open analysis frameworks ctools and Gammapy. First, we show that, when using these tools in a classical analysis approach and comparing to the proprietary H.E.S.S. analysis framework, virtually identical high-level analysis results are obtained. We then describe the construction of a generic background model from data of H.E.S.S. observations, and demonstrate that a 3D likelihood analysis using this background model yields high-level analysis results that are highly compatible with those obtained from the classical analyses. (abridged)
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
