PQS (software)

PQS
Developer(s)	Parallel Quantum Solutions
Stable release	PQS ab initio v. 4.0
Operating system	Linux, Microsoft Windows, Mac OS
License	Commercial
Website	www.pqs-chem.com

Last updated June 13, 2023

PQS is a general purpose quantum chemistry program. Its roots go back to the first ab initio gradient program developed in Professor Peter Pulay's group but now it is developed and distributed commercially by Parallel Quantum Solutions. There is a reduction in cost for academic users and a site license. Its strong points are geometry optimization, NMR chemical shift calculations, and large MP2 calculations, and high parallel efficiency on computing clusters. It includes many other capabilities including Density functional theory, the semiempirical methods, MINDO/3, MNDO, AM1 and PM3, Molecular mechanics using the SYBYL 5.0 Force Field, the quantum mechanics/molecular mechanics mixed method using the ONIOM method, natural bond orbital (NBO) analysis and COSMO solvation models. Recently, a highly efficient parallel CCSD(T) code for closed shell systems has been developed. This code includes many other post Hartree–Fock methods: MP2, MP3, MP4, CISD, CEPA, QCISD and so on.

History

The origin of PQS program was developed by Meyer^{[ who? ]} and Pulay in the late 1960s. They both were at the Max-Planck Institute for Physics and Astrophysics in Munich when they began to write a new ab initio program.^[1] The main purpose was to establish new ab initio techniques. Pulay and Meyer had slightly different interests.^[1] Pulay was interested in implementing gradient geometry optimization, analytical energy derivatives (force), and force constant calculations via the numerical differentiation of analytical forces, while Meyer enthused about the coupled-electron pair approximation (CEPA), spin density calculation, and extremely accurate correlation methods like pseudonatural orbital-configuration interaction (PNO-CI).^[1] At that time, analytical gradients were limited to closed-shell Hartree-Fock wavefunctions. However, they were able to do it for unrestricted (UHF) and restricted open-shell (ROHF) methods in 1970. The first version of the code was completed in 1969 at Max-Planck institute and University of Stuttgart. Then, Meyer named it “MOLPRO” and used Gaussian lobe basis sets.^[1] In the 1970s, the current version of MOLPRO added a number of advanced methods such as multiconfiguration self-consistent field (MC-SCF) and internally contracted multireference configuration interaction (MR-CI). Simultaneously, in the 1980s, MOLPRO was extended and mostly rewritten by Hans-Joachim Werner, Peter Knowles and Meyer's coworkers.^[1]

Meanwhile, in 1976, Pulay was visiting Boggs at the University of Texas, Austin and Schaefer at the University of California. They wrote a new program called TEXAS based on the original MOLPRO and replaced Gaussian lobe functions with the standard Gaussian functions.^[1] TEXAS emphasized large molecules, SCF convergence, geometry optimization techniques, and vibrational spectroscopy-related calculations. From 1982 onward, the program was further developed at the University of Arkansas.^[1]

The primary significant expansion was the usage of a few new electron correlation methods by Saebo and a first- order MC-SCF program by Hamilton. A critical option was the implementation of the first practical gauge-invariant atomic orbital (GIAO) NMR program by Wolinski, who additionally included a highly efficient integral package.^[1] Bofill executed an unhindered natural orbital-complete active space (UNO-CAS) program including analytical gradients; this is a minimal-cost alternative to MC-SCF and works just as well in most cases. TEXAS was initially parallelized in 1995–1996 on a cluster of 10 IBM RS6000 workstations.^[1]

In 1996, Baker joined Pulay and, around the same time, Intel brought out the Pentium Pro, a processor for PC that was competitive with low-end workstations and less costly by around an order of magnitude. Understanding the capability of this improvement for computational chemistry, PQS was formed and a SBIR grant application was submitted in July 1997 for the commercial development of PC clusters for parallel ab initio calculations.^[1] In the meantime, the Pulay group, financed by a National Science Foundation grant, set about constructing a Linux cluster using 300 MHz Pentium II processors. By a fortunate circumstance, a few capable and PC proficient graduate students were available, specifically Magyarfalvi and Shirel.^[1] The PC cluster was a complete success, and significantly outperformed the IBM Workstation cluster that was the group's computational mainstay at a small amount of its expense.^[1]

The PQS programming was demonstrated on the TEXAS code and parts of it, chiefly the NMR code, were authorized to PQS from the University of Arkansas.^[1] Much of the code was widely changed to comply with the twin points of (a) having all major functionality fully parallel; and (b) having the capacity to routinely perform calculations on extensive systems.^[1] They aimed primarily for a modest level of parallelism (from 8 to 32 CPUs), as this is the most widely recognized size for an individual or group resource. Indeed, even on very large clusters it is normal for any given user to be allocated only a percentage of the available processors.^[1]

Features

The basis capabilities in high-level correlated energies for PQS ab initio v. 4.0 include MP3, MP4, CID, CISD, CEPA-0, CEPA-2, QCISD, QCISD(T), CCD, CCSD and CCSD(T) wavefunctions; enforced geometry optimization (used, among other things, to simulate the results of Atomic Force Microscopy (AFM) experiments); full-accuracy, canonical UMP2 energies and analytical polarizabilities and hyperpolarizabilities for HF and DFT wavefunctions.^[2]

An efficient vectorized Gaussian integral package allowing high angular momentum basis functions and general contractions.
Abelian point group symmetry throughout; utilizes full point group symmetry (up to Ih) for geometry optimization step and Hessian (2nd derivative) CPHF.
Closed-shell (RHF) and open-shell (UHF) SCF energies and gradients, including several initial wavefunction guess options. Improved SCF convergence for open-shell systems.
Closed-shell (RHF) and open-shell (UHF) density functional energies and gradients including all popular exchange-correlation functionals: VWN, B88, OPTX, LYP, P86, PW91, PBE, B97, HCTH, B3LYP, make up your own functional etc.
Fast and accurate pure DFT energies and gradients for large basis sets using the Fourier Transform Coulomb (FTC) method.
Productive, adaptable geometry optimization for every one of these methods including Eigenvector Following (EF) algorithm for minimization and saddle-point search, Pulay's GDIIS algorithm for minimization, use of Cartesian, Z-matrix and delocalized internal coordinates. Includes new coordinates for efficient optimization of molecular clusters and adsorption/reaction on model surfaces.^[3]
Full range of geometrical constraints including fixed distances, planar bends, torsions and out-of-plane bends between any atoms in the molecule and frozen (fixed) atoms. Atoms involved in constraints do not need to be formally bonded and - unlike with a Z matrix - desired constraints do not need to be satisfied in the starting geometry.
Expository second subsidiaries for every one of these systems, including the calculation of vibrational frequencies, IR intensities and thermodynamic analysis.^[3]
Efficient NMR Chemical Shifts for closed-shell HF and DFT wavefunctions.
A full range of effective core potentials (ECPs), both relativistic and non-relativistic, with energies, gradients, analytical second derivatives and NMR.
Closed-shell MP2 energies and analytical gradients and dual-basis MP2 energies; numerical MP2 second derivatives.
Potential scan, including scan + optimization of all other degrees of freedom.
Reaction Path (IRC) following using either Z-matrix, Cartesian or mass-weighted Cartesian coordinates.
Conductor-like screening solvation model (COSMO) including energies, analytical gradients, numerical second derivatives and NMR.
Population analysis, including bond orders and atomic valencies (with free valencies for open-shell systems); CHELP and Cioslowski charges.
Weinhold's Natural Bond Order (NBO) analysis, including natural population and steric analysis.
Properties module with charge, spin-density and electric field gradient at the nucleus.
Polarizabilities and dipole and polarizability derivatives; Raman intensities.
Full Semiempirical package, both open (unrestricted) and closed-shell energies and gradients, including MINDO/3, MNDO, AM1 and PM3. For the latter, all main group elements through the fourth row (except the noble gases) as well as Zinc and Cadmium, have been parametrized.
Molecular Mechanics using the Sybyl 5.2 and UFF Force Fields.
QM/MM using the ONIOM method.
Molecular dynamics using the simple Verlet algorithm.
Pople-style input for quick input generation and compatibility with other programs.
Graphical input generation and display
All major ab initio functionality is fully parallel (except MP2 gradients which is serial only - parallel version under development).
Calculate molecular structure and vibrational spectra for transition state, infrared (IR), Raman and Vibrational circular dichroism (VCD).^[4]

Related Research Articles

Computational chemistry is a branch of chemistry that uses computer simulation to assist in solving chemical problems. It uses methods of theoretical chemistry, incorporated into computer programs, to calculate the structures and properties of molecules, groups of molecules, and solids. It is essential because, apart from relatively recent results concerning the hydrogen molecular ion, the quantum many-body problem cannot be solved analytically, much less in closed form. While computational results normally complement the information obtained by chemical experiments, it can in some cases predict hitherto unobserved chemical phenomena. It is widely used in the design of new drugs and materials.

MOLPRO is a software package used for accurate ab initio quantum chemistry calculations. It is developed by Peter Knowles at Cardiff University and Hans-Joachim Werner at Universität Stuttgart in collaboration with other authors.

Gaussian is a general purpose computational chemistry software package initially released in 1970 by John Pople and his research group at Carnegie Mellon University as Gaussian 70. It has been continuously updated since then. The name originates from Pople's use of Gaussian orbitals to speed up molecular electronic structure calculations as opposed to using Slater-type orbitals, a choice made to improve performance on the limited computing capacities of then-current computer hardware for Hartree–Fock calculations. The current version of the program is Gaussian 16. Originally available through the Quantum Chemistry Program Exchange, it was later licensed out of Carnegie Mellon University, and since 1987 has been developed and licensed by Gaussian, Inc.

Q-Chem is a general-purpose electronic structure package featuring a variety of established and new methods implemented using innovative algorithms that enable fast calculations of large systems on various computer architectures, from laptops and regular lab workstations to midsize clusters and HPCC, using density functional and wave-function based approaches. It offers an integrated graphical interface and input generator; a large selection of functionals and correlation methods, including methods for electronically excited states and open-shell systems; solvation models; and wave-function analysis tools. In addition to serving the computational chemistry community, Q-Chem also provides a versatile code development platform.

In computational chemistry, post–Hartree–Fock (post-HF) methods are the set of methods developed to improve on the Hartree–Fock (HF), or self-consistent field (SCF) method. They add electron correlation which is a more accurate way of including the repulsions between electrons than in the Hartree–Fock method where repulsions are only averaged.

<span class="mw-page-title-main">MOLCAS</span> Computational chemistry software

MOLCAS is an ab initio computational chemistry program, developed as a joint project by a number of international institutes. MOLCAS is developed by scientists to be used by scientists. It is not primarily a commercial product and it is not sold in order to produce a fortune for its owner.

TURBOMOLE is an ab initio computational chemistry program that implements various quantum chemistry methods. It was initially developed by the group of Prof. Reinhart Ahlrichs at the University of Karlsruhe. In 2007, TURBOMOLE GmbH, founded by R. Ahlrichs, F. Furche, C. Hättig, W. Klopper, M. Sierka, and F. Weigend, took over the responsibility for the coordination of the scientific development of TURBOMOLE program, for which the company holds all copy and intellectual property rights. In 2018 David P. Tew joined the TURBOMOLE GmbH. Since 1987, this program is one of the useful tools as it involves in many fields of research including heterogeneous and homogeneous catalysis, organic and inorganic chemistry, spectroscopy as well as biochemistry. This can be illustrated by citation records of Ahlrich's 1989 publication which is more than 6700 times as of 18 July 2020. In the year 2014, the second Turbomole article has been published. The number of citations from both papers indicates that the Turbomole's user base is expanding.

CRYSTAL is a quantum chemistry ab initio program, designed primarily for calculations on crystals, slabs and polymers using translational symmetry, but it can also be used for single molecules. It is written by V.R. Saunders, R. Dovesi, C. Roetti, R. Orlando, C.M. Zicovich-Wilson, N.M. Harrison, K. Doll, B. Civalleri, I.J. Bush, Ph. D’Arco, and M. Llunell from Theoretical Chemistry Group at the University of Torino and the Computational Materials Science Group at the Daresbury Laboratory near Warrington in Cheshire, England. The current version is CRYSTAL17. Earlier versions were CRYSTAL88, CRYSTAL92, CRYSTAL95, CRYSTAL98, CRYSTAL03, CRYSTAL06, CRYSTAL09 and CRYSTAL14.

Jaguar is a computer software package used for ab initio quantum chemistry calculations for both gas and solution phases. It is commercial software marketed by the company Schrödinger. The program was originated in research groups of Richard Friesner and William Goddard and was initially called PS-GVB.

The COLUMBUS PROGRAMS are a computational chemistry software suite for calculating ab initio molecular electronic structures, designed as a collection of individual programs communicating through files. The programs focus on extended multi-reference calculations of atomic and molecular ground and excited states. In addition to standard classes of reference wave functions such as CAS and RAS, calculations can be performed with selected configurations. Some features employ the atomic orbital integrals and gradient routines from the Dalton as well as MOLCAS program suites. COLUMBUS is distributed open-source under the LGPL license.

General Atomic and Molecular Electronic Structure System (GAMESS-UK) is a computer software program for computational chemistry. The original code split in 1981 into GAMESS-UK and GAMESS (US) variants, which now differ significantly. Many of the early developments in the UK version arose from the earlier UK based ATMOL program, which, unlike GAMESS, lacked analytical gradients for geometry optimisation.

<span class="mw-page-title-main">Spartan (chemistry software)</span>

Spartan is a molecular modelling and computational chemistry application from Wavefunction. It contains code for molecular mechanics, semi-empirical methods, ab initio models, density functional models, post-Hartree–Fock models, and thermochemical recipes including G3(MP2) and T1. Quantum chemistry calculations in Spartan are powered by Q-Chem.

Ab initio quantum chemistry methods are computational chemistry methods based on quantum chemistry. The term ab initio was first used in quantum chemistry by Robert Parr and coworkers, including David Craig in a semiempirical study on the excited states of benzene. The background is described by Parr. Ab initio means "from first principles" or "from the beginning", implying that the only inputs into an ab initio calculation are physical constants. Ab initio quantum chemistry methods attempt to solve the electronic Schrödinger equation given the positions of the nuclei and the number of electrons in order to yield useful information such as electron densities, energies and other properties of the system. The ability to run these calculations has enabled theoretical chemists to solve a range of problems and their importance is highlighted by the awarding of the Nobel prize to John Pople and Walter Kohn.

The fragment molecular orbital method (FMO) is a computational method that can compute very large molecular systems with thousands of atoms using ab initio quantum-chemical wave functions.

Quantum chemistry composite methods are computational chemistry methods that aim for high accuracy by combining the results of several calculations. They combine methods with a high level of theory and a small basis set with methods that employ lower levels of theory with larger basis sets. They are commonly used to calculate thermodynamic quantities such as enthalpies of formation, atomization energies, ionization energies and electron affinities. They aim for chemical accuracy which is usually defined as within 1 kcal/mol of the experimental value. The first systematic model chemistry of this type with broad applicability was called Gaussian-1 (G1) introduced by John Pople. This was quickly replaced by the Gaussian-2 (G2) which has been used extensively. The Gaussian-3 (G3) was introduced later.

In the field of computational chemistry, energy minimization is the process of finding an arrangement in space of a collection of atoms where, according to some computational model of chemical bonding, the net inter-atomic force on each atom is acceptably close to zero and the position on the potential energy surface (PES) is a stationary point. The collection of atoms might be a single molecule, an ion, a condensed phase, a transition state or even a collection of any of these. The computational model of chemical bonding might, for example, be quantum mechanics.

CP2K is a freely available (GPL) quantum chemistry and solid state physics program package, written in Fortran 2008, to perform atomistic simulations of solid state, liquid, molecular, periodic, material, crystal, and biological systems. It provides a general framework for different methods: density functional theory (DFT) using a mixed Gaussian and plane waves approach (GPW) via LDA, GGA, MP2, or RPA levels of theory, classical pair and many-body potentials, semi-empirical and tight-binding Hamiltonians, as well as Quantum Mechanics/Molecular Mechanics (QM/MM) hybrid schemes relying on the Gaussian Expansion of the Electrostatic Potential (GEEP). The Gaussian and Augmented Plane Waves method (GAPW) as an extension of the GPW method allows for all-electron calculations. CP2K can do simulations of molecular dynamics, metadynamics, Monte Carlo, Ehrenfest dynamics, vibrational analysis, core level spectroscopy, energy minimization, and transition state optimization using NEB or dimer method.

The ONIOM method is a computational approach developed by Morokuma and co-workers. ONIOM is a hybrid method that enables different ab initio, semi-empirical, or molecular mechanics methods to be applied to different parts of a molecule/system in combination to produce reliable geometry and energy at reduced computational cost.

<span class="mw-page-title-main">Newton-X</span> Molecular dynamics simulation software

Newton-X is a general program for molecular dynamics simulations beyond the Born-Oppenheimer approximation. It has been primarily used for simulations of ultrafast processes in photoexcited molecules. It has also been used for simulation of band envelops of absorption and emission spectra.

Python-based Simulations of Chemistry Framework (PySCF) is an ab initio computational chemistry program natively implemented in Python program language. The package aims to provide a simple, light-weight and efficient platform for quantum chemistry code developing and calculation. It provides various functions to do the Hartree–Fock, MP2, density functional theory, MCSCF, coupled cluster theory at non-relativistic level and 4-component relativistic Hartree–Fock theory. Although most functions are written in Python, the computation critical modules are intensively optimized in C. As a result, the package works as efficient as other C/Fortran-based quantum chemistry program. PySCF is developed by Qiming Sun. PySCF2.0 is the latest version of the program.

References

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 Baker, Jon; Wolinski, Krzysztof; Malagoli, Massimo; Kinghorn, Don; Wolinski, Pawel; Magyarfalvi, Gábor; Saebo, Svein; Janowski, Tomasz; Pulay, Peter (2008). "Quantum chemistry in parallel with PQS". Journal of Computational Chemistry. 30 (2): 317–335. doi:10.1002/jcc.21052. PMID 18615419. S2CID 10795179.
↑ "PQS ab initio v. 4.0". Parallel Quantum Solutions. Turnkey Computational Chemistry.
1 2 "PQS Ab Initio Program Package" (PDF). Parallel Quantum Solution.
↑ Leszczynski, Jerzy; Kaczmarek-Kędziera, Anna; G. Papadopoulos, Manthos; Reis, Heribert; J. Sadlej, Andrzej; K. Shukla, Manoj (2012-01-13). Handbook of Computational Chemistry (2012 ed.). Springer. ISBN 9789400707115.

External links

Parallel Quantum Solutions

This page is based on this Wikipedia article
Text is available under the CC BY-SA 4.0 license; additional terms may apply.
Images, videos and audio are available under their respective licenses.

[pqs1-1] 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 Baker, Jon; Wolinski, Krzysztof; Malagoli, Massimo; Kinghorn, Don; Wolinski, Pawel; Magyarfalvi, Gábor; Saebo, Svein; Janowski, Tomasz; Pulay, Peter (2008). "Quantum chemistry in parallel with PQS". Journal of Computational Chemistry. 30 (2): 317–335. doi:10.1002/jcc.21052. PMID 18615419. S2CID 10795179.

[2] "PQS ab initio v. 4.0". Parallel Quantum Solutions. Turnkey Computational Chemistry.

[pqs2-3] 1 2 "PQS Ab Initio Program Package" (PDF). Parallel Quantum Solution.

[4] Leszczynski, Jerzy; Kaczmarek-Kędziera, Anna; G. Papadopoulos, Manthos; Reis, Heribert; J. Sadlej, Andrzej; K. Shukla, Manoj (2012-01-13). Handbook of Computational Chemistry (2012 ed.). Springer. ISBN 9789400707115.

[1]

[2]

[3]

[4]