{"title":"Strain tunable pudding-mold-type band structure and thermoelectric properties of SnP3 monolayer","authors":"Shasha Wei, Cong Wang, S. Fan, G. Gao","doi":"10.1063/5.0003241","DOIUrl":"https://doi.org/10.1063/5.0003241","url":null,"abstract":"Recent studies indicated the interesting metal-to-semiconductor transition when layered bulk GeP3 and SnP3 are restricted to the monolayer or bilayer, and SnP3 monolayer has been predicted to possess high carrier mobility and promising thermoelectric performance. Here, we investigate the biaxial strain effect on the electronic and thermoelectric properties of SnP3 monolayer. Our first-principles calculations combined with Boltzmann transport theory indicate that SnP3 monolayer has the pudding-mold-type valence band structure, giving rise to a large p-type Seebeck coefficient and a high p-type power factor. The compressive biaxial strain can decrease the energy gap and result in the metallicity. In contrast, the tensile biaxial strain increases the energy gap, and increases the n-type Seebeck coefficient and decreases the n-type electrical conductivity. Although the lattice thermal conductivity becomes larger at a tensile biaxial strain due to the increased maximum frequency of the acoustic phonon modes and the increased phonon group velocity, it is still low, only e.g. 3.1 W/(mK) at room temperature with the 6% tensile biaxial strain. Therefore, SnP3 monolayer is a good thermoelectric material with low lattice thermal conductivity even at the 6% tensile strain, and the tensile strain is beneficial to the increase of the n-type Seebeck coefficient.","PeriodicalId":8424,"journal":{"name":"arXiv: Computational Physics","volume":"160 2","pages":""},"PeriodicalIF":0.0,"publicationDate":"2019-12-27","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"91440220","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Changli Ma, He Cheng, Taisen Zuo, Guisheng Jiao, Zehua Han
{"title":"NeuDATool: An open source neutron data analysis tools, supporting GPU hardware acceleration, and across-computer cluster nodes parallel","authors":"Changli Ma, He Cheng, Taisen Zuo, Guisheng Jiao, Zehua Han","doi":"10.1063/1674-0068/CJCP2005077","DOIUrl":"https://doi.org/10.1063/1674-0068/CJCP2005077","url":null,"abstract":"Empirical potential structure refinement (EPSR) is a neutron scattering data analysis algorithm and a software package. It was developed by the British spallation neutron source (ISIS) Disordered Materials Group in 1980s, and aims to construct the most-probable atomic structures of disordered liquids. It has been extensively used during the past decades, and has generated reliable results. However, it is programmed in Fortran and implements a shared-memory architecture with OpenMP. With the extensive construction of supercomputer clusters and the widespread use of graphics processing unit (GPU) acceleration technology, it is now necessary to rebuild the EPSR with these techniques in the effort to improve its calculation speed. In this study, an open source framework NeuDATool is proposed. It is programmed in the object-oriented language C++, can be paralleled across nodes within a computer cluster, and supports GPU acceleration. The performance of NeuDATool has been tested with water and amorphous silica neutron scattering data. The test shows that the software could reconstruct the correct microstructure of the samples, and the calculation speed with GPU acceleration could increase by more than 400 times compared with CPU serial algorithm at a simulation box consists about 100 thousand atoms. NeuDATool provides another choice for scientists who are familiar with C++ programming and want to define specific models and algorithms for their analyses.","PeriodicalId":8424,"journal":{"name":"arXiv: Computational Physics","volume":"14 1","pages":""},"PeriodicalIF":0.0,"publicationDate":"2019-12-17","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"86700909","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Rigidly rotating gravitationally bound systems of point particles, compared to polytropes","authors":"Yngve Hopstad, J. Myrheim","doi":"10.1142/S0129183120500904","DOIUrl":"https://doi.org/10.1142/S0129183120500904","url":null,"abstract":"In order to simulate rigidly rotating polytropes we have simulated systems of $N$ point particles, with $N$ up to 1800. Two particles at a distance $r$ interact by an attractive potential $-1/r$ and a repulsive potential $1/r^2$. The repulsion simulates the pressure in a polytropic gas of polytropic index $3/2$. We take the total angular momentum $L$ to be conserved, but not the total energy $E$. The particles are stationary in the rotating coordinate system. The rotational energy is $L^2/(2I)$ where $I$ is the moment of inertia. Configurations where the energy $E$ has a local minimum are stable. In the continuum limit $Ntoinfty$ the particles become more and more tightly packed in a finite volume, with the interparticle distances decreasing as $N^{-1/3}$. We argue that $N^{-1/3}$ is a good parameter for describing the continuum limit. We argue further that the continuum limit is the polytropic gas of index $3/2$. For example, the density profile of the nonrotating gas approaches that computed from the Lane--Emden equation describing the nonrotating polytropic gas. In the case of maximum rotation the instability occurs by the loss of particles from the equator, which becomes a sharp edge, as predicted by Jeans in his study of rotating polytropes. We describe the minimum energy nonrotating configurations for a number of small values of $N$.","PeriodicalId":8424,"journal":{"name":"arXiv: Computational Physics","volume":"26 1","pages":""},"PeriodicalIF":0.0,"publicationDate":"2019-11-04","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"88278370","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Dissipative modes, Purcell factors, and directional beta factors in gold bowtie nanoantenna structures","authors":"Chelsea Carlson, S. Hughes","doi":"10.1103/physrevb.102.155301","DOIUrl":"https://doi.org/10.1103/physrevb.102.155301","url":null,"abstract":"We present a detailed quasinormal mode analysis of gold bowtie nanoantennas, and highlight the unusual role of the substrate and the onset of multi-mode behaviour. In particular, we show and explain why the directional radiatiave beta factor is completely dominated by emission into the substrate, and explain how the beta factors and quenching depend on the underlying mode properties. We also quantitatively explain the generalized Purcell factors and explore the role of gap size and substrate in detail. These rich modal features are essential to understand for future applications such as sensing, lasing, and quantum information processing, for example in the design of efficient single photon emitters.","PeriodicalId":8424,"journal":{"name":"arXiv: Computational Physics","volume":"296 1","pages":""},"PeriodicalIF":0.0,"publicationDate":"2019-10-22","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"90801236","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"A relativistic particle pusher for ultra-strong electromagnetic fields","authors":"Jérôme Pétri","doi":"10.1017/S0022377820000719","DOIUrl":"https://doi.org/10.1017/S0022377820000719","url":null,"abstract":"Abridged. Kinetic plasma simulations are nowadays commonly used to study a wealth of non-linear behaviours and properties in laboratory and space plasmas. In particular, in high-energy physics and astrophysics, the plasma usually evolves in ultra-strong electromagnetic fields produced by intense laser beams for the former or by rotating compact objects such as neutron stars and black holes for the latter. In these ultra-strong electromagnetic fields, the gyro-period is several orders of magnitude smaller than the timescale on which we desire to investigate the plasma evolution. Some approximations are required like for instance artificially decreasing the electromagnetic field strength which is certainly not satisfactory. The main flaw of this downscaling is that it cannot reproduce particle acceleration to ultra-relativistic speeds with Lorentz factor above $gamma approx 10^3-10^4$. In this paper, we design a new algorithm able to catch particle motion and acceleration to Lorentz factor up to $10^{15}$ or even higher by using Lorentz boosts to special frames where the electric and magnetic field are parallel. Assuming that these fields are locally uniform in space and constant in time, we solve analytically the equation of motion in a tiny region smaller than the length scale of the spatial and temporal gradient of the field.","PeriodicalId":8424,"journal":{"name":"arXiv: Computational Physics","volume":"1 1","pages":""},"PeriodicalIF":0.0,"publicationDate":"2019-10-10","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"88760791","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Maolin Bo, Hanze Li, Zhongkai Huang, Lei Li, Chuang Yao
{"title":"Bond relaxation and electronic properties of two-dimensional Sb/MoSe2 and Sb/MoTe2 van der Waals heterostructures","authors":"Maolin Bo, Hanze Li, Zhongkai Huang, Lei Li, Chuang Yao","doi":"10.1063/1.5130533","DOIUrl":"https://doi.org/10.1063/1.5130533","url":null,"abstract":"Van der Waals heterostructures have recently garnered interest for application in high-performance photovoltaic materials. Consequently, understanding the basic electronic characteristics of these heterostructures is important for their utilisation in optoelectronic devices. The electronic structures and bond relaxation of two-dimensional (2D) Sb/transition metal disulfides (TMDs, MoSe2, and MoTe2) van der Waals heterostructures were systematically studied using the bond-charge (BC) correlation and hybrid density functional theory. We found that the Sb/MoSe2 and Sb/MoTe2 heterostructures had indirect band gaps of 0.701 and 0.808 eV, respectively; further, these heterostructures effectively modulated the band gaps of MoSe2 (1.463 eV) and MoTe2 (1.173 eV). The BC correlation revealed four bonding and electronic contributions (electron-holes, antibonding, nonbonding, and bonding states) of the heterostructures. Our results provide an in-depth understanding of the Sb/TMD van der Waals heterojunction, which should be utilised to design 2D metal/semiconductor-based devices.","PeriodicalId":8424,"journal":{"name":"arXiv: Computational Physics","volume":"119 1","pages":""},"PeriodicalIF":0.0,"publicationDate":"2019-10-04","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"86801512","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Warm dense matter simulation via electron temperature dependent deep potential molecular dynamics","authors":"Yuzhi Zhang, Chang Gao, Qianrui Liu, Linfeng Zhang, Han Wang, Mohan Chen","doi":"10.1063/5.0023265","DOIUrl":"https://doi.org/10.1063/5.0023265","url":null,"abstract":"Simulating warm dense matter that undergoes a wide range of temperatures and densities is challenging. Predictive theoretical models, such as quantum-mechanics-based first-principles molecular dynamics (FPMD), require a huge amount of computational resources. Herein, we propose a deep learning based scheme, called electron temperature dependent deep potential molecular dynamics (TDDPMD), for efficiently simulating warm dense matter with the accuracy of FPMD. The TDDPMD simulation is several orders of magnitudes faster than FPMD, and, unlike FPMD, its efficiency is not affected by the electron temperature. We apply the TDDPMD scheme to beryllium (Be) in a wide range of temperatures (0.4 to 2500 eV) and densities (3.50 to 8.25 g/cm$^3$). Our results demonstrate that the TDDPMD method not only accurately reproduces the structural properties of Be along the principal Hugoniot curve at the FPMD level, but also yields even more reliable diffusion coefficients than typical FPMD simulations due to its ability to simulate larger systems with longer time.","PeriodicalId":8424,"journal":{"name":"arXiv: Computational Physics","volume":"30 1","pages":""},"PeriodicalIF":0.0,"publicationDate":"2019-08-31","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"76895030","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"A full Stokes subgrid model for simulation of grounding line migration in ice sheets using Elmer/ICE(v8.3)","authors":"Gong Cheng, Per Lötstedt, L. von Sydow","doi":"10.5194/gmd-2019-244","DOIUrl":"https://doi.org/10.5194/gmd-2019-244","url":null,"abstract":"Abstract. The full Stokes equations are solved by a finite element method for simulation of large ice sheets and glaciers. The simulation is particularly sensitive to the discretization of the grounding line which separates the ice resting on the bedrock and the ice floating on water and is moving in time. The boundary conditions at the ice base are enforced by Nitsche's method and a subgrid treatment of the elements in the discretization close to the grounding line. Simulations with the method in two dimensions for an advancing and a retreating grounding line illustrate the performance of the method. It is implemented in the two dimensional version of the open source code Elmer/ICE.\u0000","PeriodicalId":8424,"journal":{"name":"arXiv: Computational Physics","volume":"18 1","pages":""},"PeriodicalIF":0.0,"publicationDate":"2019-08-28","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"73328261","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"New stable, explicit, first order method to solve the heat conduction equation","authors":"E. Kovács","doi":"10.32973/JCAM.2020.001","DOIUrl":"https://doi.org/10.32973/JCAM.2020.001","url":null,"abstract":"We introduce a novel explicit and stable numerical algorithm to solve the spatially discretized heat or diffusion equation. We compare the performance of the new method with analytical and numerical solutions. We show that the method is first order in time and can give approximate results for extremely large systems faster than the commonly used explicit or implicit methods.","PeriodicalId":8424,"journal":{"name":"arXiv: Computational Physics","volume":"402 1","pages":""},"PeriodicalIF":0.0,"publicationDate":"2019-08-26","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"77158134","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"An implicit P-multigrid flux reconstruction method for simulation of locally preconditioned unsteady Navier–Stokes equations at low Mach numbers","authors":"Lai Wang, Meilin Yu","doi":"10.13016/M2ZXLR-YRBL","DOIUrl":"https://doi.org/10.13016/M2ZXLR-YRBL","url":null,"abstract":"The authors gratefully acknowledge the support of the Office of Naval \u0000Research through the award N00014-16-1-2735, and the faculty startup support from the department of mechanical engineering at the University of \u0000Maryland, Baltimore County (UMBC). The hardware used in the computational studies is part of the UMBC High Performance Computing Facility \u0000(HPCF). The facility is supported by the U.S. National Science Foundation \u0000through the MRI program (grant nos. CNS-0821258, CNS-1228778, and \u0000OAC-1726023) and the SCREMS program (grant no. DMS-0821311), with \u0000additional substantial support from UMBC.","PeriodicalId":8424,"journal":{"name":"arXiv: Computational Physics","volume":"69 1","pages":""},"PeriodicalIF":0.0,"publicationDate":"2019-08-11","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"84256305","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}