MoCalc: a new graphical user interface for molecular calculations

A new computer program called MoCalc (Molecular Calculations) has been designed to help the computational chemistry practitioner in the task of performing and analyzing molecular calculations. MoCalc is a graphical user interface for the MO calculation programs Gamess and Mopac, and uses Rasmol and Babel for molecule display and file conversion, respectively. In its initial version, MoCalc can execute the following operations: (a) create and handle Gamess and Mopac input files; (b) import any kind of molecular geometry supported by Babel and paste it as Cartesian, internal, or Gaussian-type coordinates on the input file; (c) convert Gamess and Mopac output files to inputs of both programs; (d) edit and validate the keywords that control the Gamess and Mopac calculation procedure; (e) display the input (Mopac) and output (Gamess and Mopac) molecular geometries; (f) run single or multiple (batch) calculations, either interactively or in background; (g) automatically open the output files as soon as the calculation finishes; (h) extract results from the output files, such as energy, charges, dipole, population analysis, wave function, bond orders, and valence analysis, and display them in spreadsheets; (i) calculate reactivity indices derived from the frontier orbital theory and the root-mean-square (rms) deviation of input and output geometries. All the results generated by MoCalc can be promptly transferred to text editors and electronic spreadsheets, which facilitate a detailed subsequent analysis and the publication of the results. MoCalc can also perform graphical and numerical comparative analysis of the some results when more than one output file is loaded. The program was coded in Visual Basic and runs in Windows 95/98/NT4/ME/2000/XP environments.     

Computational chemistry on a PC

We describe a package of some IBM PC programs that may find application in computer-aided molecular design. PCGEOM constructs and visualizes molecular models from bond lengths, bond angles, and dihedral angles, from Cartesian coordinates, or from stored fragments. It may prepare output files to be used as input for other programs, like CNDOB (conventional CNDO /2) or PCMEP using the bond increment (BI ) method for the calculation of molecular electrostatic potentials. PCPROT is in preparation and will use Protein Data Bank coordinates to visualize and manipulate protein molecular models. Starting from these, it will calculate electrostatic potentials using the BI method and/or monopoles adjusted to reproduce ab initio values for amino acid residues. FSCF is based on a CNDO -type approximation and uses strictly localized molecular orbitals in order to partition large molecules into a central fragment, a polarizable region, and a fully transferable environment. The partition allows one to handle relatively large systems with up to 200 atoms. To illustrate applications, we present estimation of relative inhibitory potencies of a series of substituted triazines on chicken liver dihydrofolate reductase.     

Programmed chemical systems: multiple subprograms and multiple processing/expression of molecular information

Programmed chemical systems rest on the structural information stored in a molecular framework and on its reading and processing through non-covalent interactional algorithms to yield specific supramolecular entities. Beyond single-code self-assembly, which generates exclusively a single, specific superstructure, several codes may be implemented in the same overall program, thus opening the possibility to perform multiprogramming. Furthermore, the reading and processing of the same structural information through different interactional algorithms may lead to several different output entities, amounting to multiple expression of molecular information. Such features are revealed in the formation of double helicates, the assembly of metallosupramolecular architectures, and the differential reading of hydrogen bonding patterns in a molecular strand. They open novel perspectives within the framework of programmed chemical systems, concerning multiple processing capacity, and have intriguing implications from the biological point of view.     

ORBKIT: A modular python toolbox for cross‐platform postprocessing of quantum chemical wavefunction data

ORBKIT is a toolbox for postprocessing electronic structure calculations based on a highly modular and portable Python architecture. The program allows computing a multitude of electronic properties of molecular systems on arbitrary spatial grids from the basis set representation of its electronic wavefunction, as well as several grid‐independent properties. The required data can be extracted directly from the standard output of a large number of quantum chemistry programs. ORBKIT can be used as a standalone program to determine standard quantities, for example, the electron density, molecular orbitals, and derivatives thereof. The cornerstone of ORBKIT is its modular structure. The existing basic functions can be arranged in an individual way and can be easily extended by user‐written modules to determine any other derived quantity. ORBKIT offers multiple output formats that can be processed by common visualization tools (VMD, Molden, etc.). Additionally, ORBKIT possesses routines to order molecular orbitals computed at different nuclear configurations according to their electronic character and to interpolate the wavefunction between these configurations. The program is open‐source under GNU‐LGPLv3 license and freely available at https://github.com/orbkit/orbkit/ . This article provides an overview of ORBKIT with particular focus on its capabilities and applicability, and includes several example calculations. © 2016 Wiley Periodicals, Inc.     

GenLocDip: A Generalized Program to Calculate and Visualize Local Electric Dipole Moments

Local dipole moments (i.e., dipole moments of atomic or molecular subsystems) are essential for understanding various phenomena in nanoscience, such as solvent effects on the conductance of single molecules in break junctions or the interaction between the tip and the adsorbate in atomic force microscopy. We introduce Gen Loc Dip , a program for calculating and visualizing local dipole moments of molecular subsystems. Gen Loc Dip currently uses the Atoms‐In‐Molecules (AIM) partitioning scheme and is interfaced to various AIM programs. This enables postprocessing of a variety of electronic structure output formats including cube and wavefunction files, and, in general, output from any other code capable of writing the electron density on a three‐dimensional grid. It uses a modified version of Bader's and Laidig's approach for achieving origin‐independence of local dipoles by referring to internal reference points which can (but do not need to be) bond critical points (BCPs). Furthermore, the code allows the export of critical points and local dipole moments into a POVray readable input format. It is particularly designed for fragments of large systems, for which no BCPs have been calculated for computational efficiency reasons, because large interfragment distances prevent their identification, or because a local partitioning scheme different from AIM was used. The program requires only minimal user input and is written in the Fortran 90 programming language. To demonstrate the capabilities of the program, examples are given for covalently and non‐covalently bound systems, in particular molecular adsorbates. © 2016 Wiley Periodicals, Inc.     

Visual Displays that Directly Interface and Provide Read‐Outs of Molecular States via Molecular Graphics Processing Units

The monitoring of molecular systems usually requires sophisticated technologies to interpret nanoscale events into electronic‐decipherable signals. We demonstrate a new method for obtaining read‐outs of molecular states that uses graphics processing units made from molecular circuits. Because they are made from molecules, the units are able to directly interact with molecular systems. We developed deoxyribozyme‐based graphics processing units able to monitor nucleic acids and output alphanumerical read‐outs via a fluorescent display. Using this design we created a molecular 7‐segment display, a molecular calculator able to add and multiply small numbers, and a molecular automaton able to diagnose Ebola and Marburg virus sequences. These molecular graphics processing units provide insight for the construction of autonomous biosensing devices, and are essential components for the development of molecular computing platforms devoid of electronics.     

WebProp: Web interface for ab initio calculation of molecular one-electron properties

This note describes the features and implementation issues of WebProp, a web-based interface for evaluating ab initio quality one-electron properties. The interface code is written in HTML and Python, while the backend is handled using Python and our indigenously developed code INDPROP for property evaluation. A novel feature of this setup is that it provides a simple interface for computing first principle one-electron properties of small to medium sized molecules. To facilitate computation of otherwise expensive calculations on large molecular systems, we employ the Molecular Tailoring Approach (MTA) developed in our laboratory to obtain the density matrix (DM). This DM is then employed for computing the one-electron properties of these systems. The backend transparently handles jobs submitted by the user and runs them either on a single machine or over a grid of compute nodes. The results of the calculations, which include the summary and the files necessary for visualization of one-electron properties, are e-mailed to the user. The user can either directly use the data or visualize it using visualization tools such as UNIVIS-2000 or Drishti.     

ExcelAutomat 1.3: Fragment analysis based on the distortion/interaction-activation strain model

The distortion/interaction-activation strain model (D/I-ASM), a fragment analysis method, is applied to study the structure–reactivity relationship in reactions. The application of D/I-ASM involves the generation of input files for points along a reaction profile, submission of input files to a quantum software package, processing of parameters from the resulting output files and generation of graphical plots. The ExcelAutomat tool (Laloo et al., J. Comput. Aided Mol. Des. 2017, 31, 667) provides a framework and library in Visual Basic for Application programming language to process such files. New routines were written in ExcelAutomat 1.3 to facilitate processing of files for D/I-ASM. The worksheet “ASM” was included where initial parameters needed can be defined. The routines for D/I-ASM were tested successfully on bimolecular nucleophilic substitution, cycloaddition, and barrierless reactions. The automation of fragment analysis by ExcelAutomat 1.3 is compatible with Microsoft Excel and LibreOffice Calc. The extensible tool processes files from Gaussian and GAMESS-US packages. © 2018 Wiley Periodicals, Inc.     

Molecular photonic logic gates

The possibility of performing logical operations at the molecular level is being actively investigated at present with the aim of developing molecular logic gates, which can be used in information technologies. In this minireview, the design algorithm of molecular logic gates is considered and the requirements on molecular systems for use as logic gates are specified. Examples of molecular logic gates performing different logical operations are given. Attention is focused on all-photonic molecular logic gates, in which light is used as an input signal for transferring the system from one state to another and for reading the output signal by absorption or luminescence. In addition, optoelectronic devices with light as the input signal and electric current as the output signal are briefly discussed.     

DataPipeline: Automated importing and fitting of large amounts of biophysical data

Raw data from experiments across the biological sciences comes in a large variety of text formats. In small or medium sized laboratories researchers often use an assorted collection of software to interpret, fit, and visualize their data. The spreadsheet is commonly the core component of such a workflow. The limitations of such programs for large amounts of heterogeneous data can be frustrating. We report the construction of DataPipeline, a desktop and command‐line application that automates the tasks of importing, fitting, and plotting of text‐based data. The software is designed to simplify the process of importing text data from various sources using simple configuration files to describe raw file formats. Once imported, curve fitting can be performed using custom fitting models designed by the user inside the application. Fitted parameters can be grouped together as new datasets to be fitted to other models and experimental uncertainties propagated to give error estimates. This software will be useful for processing of data from high through‐put biological experiments or for rapid visualization of pilot data without the need for a chain of different programs to carry out each step. DataPipeline and source code is available under an open source license. The software can be freely downloaded at http://code.google.com/p/peat/downloads/list . © 2012 Wiley Periodicals, Inc.     

Computational spectroscopy: status of problem

The possibility of constrution the system of artificial intelligence for processing of molecular spectroscopy data is discussed. The development of the theory of spectra and molecular structures as well as possibility to carry out sufficiently precise calculations make it possible to construct algorithmic systems operating with molecular models of different complication. The successive study of such models leads to the complete solution of inverse spectral problems. The state of the problem is considered.     

NSPEC: a chemical speciation program for personal computers.

The details of a computer program, NSPEC, which simulates chemical speciation in aqueous systems, are presented. NSPEC has been designed explicitly for the IBM PC family of computers and is a departure from normal speciation programs in that it is controlled through interaction with a series of menus. The program is more user-friendly than other comparable speciation programs, can write the results files in formats suitable for importing into other programs, and, in many instances gives results faster than speciation programs running on microcomputers.     

Lessons learned from comparing molecular dynamics engines on the SAMPL5 dataset

We describe our efforts to prepare common starting structures and models for the SAMPL5 blind prediction challenge. We generated the starting input files and single configuration potential energies for the host-guest in the SAMPL5 blind prediction challenge for the GROMACS, AMBER, LAMMPS, DESMOND and CHARMM molecular simulation programs. All conversions were fully automated from the originally prepared AMBER input files using a combination of the ParmEd and InterMol conversion programs. We find that the energy calculations for all molecular dynamics engines for this molecular set agree to better than 0.1 % relative absolute energy for all energy components, and in most cases an order of magnitude better, when reasonable choices are made for different cutoff parameters. However, there are some surprising sources of statistically significant differences. Most importantly, different choices of Coulomb’s constant between programs are one of the largest sources of discrepancies in energies. We discuss the measures required to get good agreement in the energies for equivalent starting configurations between the simulation programs, and the energy differences that occur when simulations are run with program-specific default simulation parameter values. Finally, we discuss what was required to automate this conversion and comparison.     

Gabedit—A graphical user interface for computational chemistry softwares

Gabedit is a freeware graphical user interface, offering preprocessing and postprocessing adapted (to date) to nine computational chemistry software packages. It includes tools for editing, displaying, analyzing, converting, and animating molecular systems. A conformational search tool is implemented using a molecular mechanics or a semiempirical potential. Input files can be generated for the computational chemistry software supported by Gabedit. Some molecular properties of interest are processed directly from the output of the computational chemistry programs; others are calculated by Gabedit before display. Molecular orbitals, electron density, electrostatic potential, nuclear magnetic resonance shielding density, and any other volumetric data properties can be displayed. It can display electronic circular dichroism, UV–visible, infrared, and Raman‐computed spectra after a convolution. Gabedit can generate a Povray file for geometry, surfaces, contours, and color‐coded planes. Output can be exported to a selection of popular image and vector graphics file formats; the program can also generate a series of pictures for animation. Quantum mechanical electrostatic potentials can be calculated using the partial charges on atoms, or by solving the Poisson equation using the multigrid method. The atoms in molecule charges can also be calculated. Gabedit is platform independent. The code is distributed under free open source X11 style license and is available at http://gabedit.sourceforge.net/ . © 2010 Wiley Periodicals, Inc. J Comput Chem, 2010     

MEPSIM: A computational package for analysis and comparison of molecular electrostatic potentials

Summary MEPSIM is a computational system which allows an integrated computation, analysis, and comparison of molecular electrostatic potential (MEP) distributions. It includes several modules. Module MEPPLA supplies MEP values for the points of a grid defined on a plane which is specified by a set of three points. The results of this program can easily be converted into MEP maps using third-parties graphical software. Module MEPMIN allows to find automatically the MEP minima of a molecular system. It supplies the cartesian coordinates of these minima, their values, and all the geometrical relationships between them (distances, angles, and dihedral angles). Module MEPCOMP computes a similarity coefficient between the MEP distributions of two molecules and finds their relative position that maximizes the similarity. Module MEPCONF performs the same process as MEPCOMP, considering not only the relative position of both molecules but also a conformational degree of freedom of one of them. The most recently developed module, MEPPAR, is another modification of MEPCOMP in order to compute the MEP similarity between two molecules, but only taking into account a particular plane. The latter module is particularly useful to compare MEP distributions generated by systems of aromatic rings. MEPSIM can use several wavefunction computation approaches to obtain MEP distributions. MEPSIM has a menu type interface to simplify the following tasks: creation of input files from output files of external programs (GAUSSIAN and AMPAC/MOPAC), setting the parameters for the current computation, and submitting jobs to the batch queues of the computer. MEPSIM has been coded in FORTRAN and its current version runs on VMS/VAX computers.     

NADA: A versatile PC based program for neutron activation data analysis

We have developed a PC based program for neutron activation data analysis using the FORTRAN and C languages. The routines are based on creating files associated with conventional ORTEC hardware and output software. The main features of the program include radionuclide identification, and the use of semi-automatic integration or the peak fitting SAMPO routine. Other developments are hard and soft copy records for detailed sample identification and particular irradiation, decay and counting procedures. Flux variations, high deadtime corrections, counting geometries, spectral and nuclear interferences, as well as uranium fission interferences are also automatically accounted for. The data output includes concentration values in %, ppm, g or ppb units with associated errors, while detection limits for each individual sample are indicated. Further data output can easily be generated which can be imported to most spreadsheet programs for various statistical uses. A future implementation to the program will include batch-file processing and automated self-absorption calculations for geological samples.