PySCF: the Python‐based simulations of chemistry framework Sun, Qiming; Berkelbach, Timothy C.; Blunt, Nick S. ...
Wiley interdisciplinary reviews. Computational molecular science,
January/February 2018, Letnik:
8, Številka:
1
Journal Article
Recenzirano
Odprti dostop
Python‐based simulations of chemistry framework (PySCF) is a general‐purpose electronic structure platform designed from the ground up to emphasize code simplicity, so as to facilitate new method ...development and enable flexible computational workflows. The package provides a wide range of tools to support simulations of finite‐size systems, extended systems with periodic boundary conditions, low‐dimensional periodic systems, and custom Hamiltonians, using mean‐field and post‐mean‐field methods with standard Gaussian basis functions. To ensure ease of extensibility, PySCF uses the Python language to implement almost all of its features, while computationally critical paths are implemented with heavily optimized C routines. Using this combined Python/C implementation, the package is as efficient as the best existing C or Fortran‐based quantum chemistry programs. In this paper, we document the capabilities and design philosophy of the current version of the PySCF package. WIREs Comput Mol Sci 2018, 8:e1340. doi: 10.1002/wcms.1340
This article is categorized under:
Structure and Mechanism > Computational Materials Science
Electronic Structure Theory > Ab Initio Electronic Structure Methods
Software > Quantum Chemistry
The PySCF package provides a Python programming environment to study the electronic structure of molecules and solids.
In this work, we demonstrate the process for porting the cloud resolving model (CRM) used in the Energy Exascale Earth System Model Multi-Scale Modeling Framework (E3SM-MMF) from its original Fortran ...code base to C++ code using a portability library. This porting process is performed using the Yet Another Kernel Library (YAKL), a simplified C++ portability library that specializes in Fortran porting. In particular, we detail our step-by-step approach for porting the System for Atmospheric Modeling (SAM), the CRM used in E3SM-MMF, using a hybrid Fortran/C++ framework that allows for systematic reproduction and correctness testing of gradually ported YAKL C++ code. Additionally, analysis is done on the performance of the ported code using OLCF’s Summit supercomputer.
The iEBE-VISHNU code package performs event-by-event simulations for relativistic heavy-ion collisions using a hybrid approach based on (2+1)-dimensional viscous hydrodynamics coupled to a hadronic ...cascade model. We present the detailed model implementation, accompanied by some numerical code tests for the package. iEBE-VISHNU forms the core of a general theoretical framework for model-data comparisons through large scale Monte-Carlo simulations. A numerical interface between the hydrodynamically evolving medium and thermal photon radiation is also discussed. This interface is more generally designed for calculations of all kinds of rare probes that are coupled to the temperature and flow velocity evolution of the bulk medium, such as jet energy loss and heavy quark diffusion.
Program title: iEBE-VISHNU
Catalogue identifier: AEYA_v1_0
Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AEYA_v1_0.html
Program obtainable from: CPC Program Library, Queen’s University, Belfast, N. Ireland
Licensing provisions: Standard CPC licence, http://cpc.cs.qub.ac.uk/licence/licence.html
No. of lines in distributed program, including test data, etc.: 5257939
No. of bytes in distributed program, including test data, etc.: 262822421
Distribution format: tar.gz
Programming language: Fortran, C++, python, bash, SQLite.
Computer: Laptop, desktop, cluster.
Operating system: Tested on GNU/Linux Ubuntu 12.04 x64, Red Hat Linux 6, Mac OS X 10.8+.
RAM: 2G bytes
Classification: 17.11, 17.16, 17.20.
External routines: GNU Scientific Library (GSL), HDF5 (Fortran and C++ enabled), Numpy
Nature of problem: Relativistic heavy-ion collisions are tiny in size (V≈10−42m3) and live in a flash (∼5×10−23s). It is impossible to use external probes to study the properties of the quark-gluon plasma (QGP), a novel state of matter created during the collisions. Experiments can only measure the momentum information of stable hadrons, who are the remnants of the collisions. In order to extract the thermal and transport properties of the QGP one needs to rely on Monte-Carlo event-by-event model simulations, which reverse-engineer the experimental measurements to the early time dynamics of the relativistic heavy-ion collisions.
Solution method: Relativistic heavy-ion collisions contain multiple stages of evolution. The physics that governs each stage is implemented into individual code components. A general driver script glues all the modular packages as a whole to perform large-scale Monte-Carlo simulations. The final results are stored into SQLite database, which supports standard querying for massive data analysis. By tuning transport coefficients of the QGP as free parameters, e.g. the specific shear viscosity η/s, we can constrain various transport properties of the QGP through model-data comparisons.
Additional comments: !!! The distribution file for this program is over 260 Mbytes and therefore is not delivered directly when download or Email is requested. Instead a html file giving details of how the program can be obtained is sent. !!!
Running time: The following running time is tested on a laptop computer with a 2.4 GHz Intel Core i5 CPU, 4 GB memory. All the C++ and Fortran codes are compiled with the GNU Compiler Collection (GCC) 4.9.2 and -O3 optimization (Table 1).
The R package lcmm provides a series of functions to estimate statistical models based on linear mixed model theory. It includes the estimation of mixed models and latent class mixed models for ...Gaussian longitudinal outcomes (hlme), curvilinear and ordinal univariate longitudinal outcomes (lcmm) and curvilinear multivariate outcomes (multlcmm), as well as joint latent class mixed models (Jointlcmm) for a (Gaussian or curvilinear) longitudinal outcome and a time-to-event outcome that can be possibly left-truncated right-censored and defined in a competing setting. Maximum likelihood esimators are obtained using a modified Marquardt algorithm with strict convergence criteria based on the parameters and likelihood stability, and on the negativity of the second derivatives. The package also provides various post-fit functions including goodness-of-fit analyses, classification, plots, predicted trajectories, individual dynamic prediction of the event and predictive accuracy assessment. This paper constitutes a companion paper to the package by introducing each family of models, the estimation technique, some implementation details and giving examples through a dataset on cognitive aging
We discuss the systematic evaluation of three-loop vacuum integrals with arbitrary masses. Using integration by parts, the general integral of this type can be reduced algebraically to a few basis ...integrals. We define a set of modified finite basis integrals that are particularly convenient for expressing renormalized quantities. The basis integrals can be computed numerically by solving coupled first-order differential equations, using as boundary conditions the analytically known special cases that depend on only one mass scale. We provide the results necessary to carry this out and introduce an implementation in the form of a public software package called 3vil (three-loop Vacuum Integral Library), which efficiently computes the numerical values of the basis integrals for any specified masses. Note that 3vil is written in c, and can be linked from c, c++, or fortran code.
We present a fast and user friendly exoplanet transit light-curve modelling package pytransit, implementing optimized versions of the Gimenez and Mandel & Agol transit models. The package offers an ...object-oriented python interface to access the two models implemented natively in fortran with OpenMP parallelization. A partial OpenCL version of the quadratic Mandel-Agol model is also included for GPU-accelerated computations. The aim of pytransit is to facilitate the analysis of photometric time series of exoplanet transits consisting of hundreds of thousands of data points, and of multipassband transit light curves from spectrophotometric observations, as a part of a researcher's programming toolkit for building complex, problem-specific analyses.
•A GPU-based discrete-finite element method is developed by using the CUDA FORTRAN.•The effectiveness and scalability of the method are validated.•The developed method is applied to tire-sand ...interaction simulations.•A speedup of more than 15 is achieved in tire-sand interaction simulation.•Simulation results of tire-sand interaction cases agree well with experimental results.
Recently, the discrete-finite element method (DEM-FEM) has proved to be an advanced technique for solid-particle interaction simulations. However, the low computational efficiency limits its applications to practical engineering problems. To achieve this end, this paper proposes a Graphics Processing Unit (GPU) based DEM-FEM, including contact detection, force calculation and information update, in the context of Compute Unified Device Architecture (CUDA) FORTRAN environment. Three numerical examples are performed to validate the efficiency, effectiveness and scalability of the developed method. On this basis, the GPU-based DEM-FEM is extended to a computing platform, and programmed into our in-house code PDFP-OVS for a typical solid-particle interaction problem, i.e. the running performance of a pneumatic tire on granular sand. Numerical result shows that a speedup of more than 15 can be achieved. The simulation results are found to be in good agreement with the experiment results in terms of the gross tractive effort, the drawbar pull and the running resistance, which validates the capacity of the platform in the travel performance analysis of the tire on granular sand.
We present a new Fortran library to evaluate all harmonic polylogarithms up to weight four numerically for any complex argument. The algorithm is based on a reduction of harmonic polylogarithms up to ...weight four to a minimal set of basis functions that are computed numerically using series expansions allowing for fast and reliable numerical results.
Program title: Chaplin
Catalogue identifier: AETC_v1_0
Program summary URL: http://cpc.cs.qub.ac.uk/summaries/AETC_v1_0.html
Program obtainable from: CPC Program Library, Queen’s University, Belfast, N. Ireland
Licensing provisions: Standard CPC licence, http://cpc.cs.qub.ac.uk/licence/licence.html
No. of lines in distributed program, including test data, etc.: 101192
No. of bytes in distributed program, including test data, etc.: 729234
Distribution format: tar.gz
Programming language: Fortran 77.
Computer: Computing systems on which Fortran 77 compilers are available.
Operating system: Operating systems on which Fortran 77 compilers are available.
Classification: 11.1.
Nature of problem:
Numerical evaluation of harmonic polylogarithms.
Solution method:
Inside the unit circle: series expansion. Outside the unit circle: inversion relations.
Restrictions:
Only harmonic polylogarithms up to weight four are supported.
Unusual features:
Allows the evaluation of HPL’s numerically for any point in the complex plane.
Running time:
Depending on the weight vector and argument of the HPL, between 0.2 and 400 μs.
The “paleo calendar effect” is a common expression for
the impact that changes in the length of months or seasons over time,
related to changes in the eccentricity of Earth's orbit and precession, ...have
on the analysis or summarization of climate-model output. This effect can
have significant implications for paleoclimate analyses. In particular,
using a “fixed-length” definition of months (i.e., defined by a fixed
number of days), as opposed to a “fixed-angular” definition (i.e., defined
by a fixed number of degrees of the Earth's orbit), leads to comparisons of
data from different positions along the Earth's orbit when comparing paleo
with modern simulations. This effect can impart characteristic spatial
patterns or signals in comparisons of time-slice simulations that otherwise
might be interpreted in terms of specific paleoclimatic mechanisms, and we
provide examples for 6, 97, 116, and 127 ka. The calendar effect is
exacerbated in transient climate simulations in which, in addition to spatial
or map-pattern effects, it can influence the apparent timing of extrema in
individual time series and the characterization of phase relationships among
series. We outline an approach for adjusting paleo simulations that have
been summarized using a modern fixed-length definition of months and that
can also be used for summarizing and comparing data archived as daily data.
We describe the implementation of this approach in a set of Fortran 90
programs and modules (PaleoCalAdjust v1.0).