ANI-1, A data set of 20 million calculated off-equilibrium conformations for organic molecules

Justin S. Smith, Olexandr Isayev, Adrian E. Roitberg · Scientific Data · 2017 · chemistry · 10.1038/sdata.2017.193 · record smith2017-ani1-dataset-organic-molecules
tag layer · index · concepts · authors · עברית
Where this record's text came from. 2026-09-07 schema-0.3 pass: every quote verified as an exact substring of the staged text and anchored at merge time; concept terms verified present; every edge endpoint an existing node.
4 quotations on this page contain a reference number where the paper printed a superscript. This journal style prints citations as raised numerals — parameters12,13 — and the extraction flattened them into the line, so the sentence reads parameters 12 , 13 . The digits are citation markers and not part of what the author wrote. The quotes keep them for the same reason they keep every other extraction fault: a quote is reproduced exactly or it proves nothing.
How to read this page. The article text below is reproduced exactly as it was staged; nothing in it has been rewritten. Coloured stretches are tagged. A solid underline means the tag records the paper's own words. A dotted underline means the tag text is the tagger's sentence about the paper — the quotation is the paper's, the wording of the tag is not. Every superscript label links to that tag's row in the table at the foot of the page. The machine-readable edition of the same thing is the TEI file beside this one, where the tags live in standOff. When we have changed our mind. Some entries in the table begin with a word in capitals and the reason for it. EVIDENCE CORRECTED — the finding still stands, but we had quoted the wrong sentence from the paper, and the quotation has been replaced with the right one. RETRACTED — the entry itself was wrong; it stays on the page, marked, with the reason, rather than being deleted. WITHDRAWN AND REPLACED — the entry was wrong, and a corrected version of it appears in the entries below. Each mark carries its own reason, in brackets, because the reasons differ from one entry to the next. An entry we have WITHDRAWN is removed from this page and kept in the record, with its reason. A line reading Audit and a date is a later re-reading of the entry against its whole sentence: DOUBTFUL means it was kept, and the reason says what is loose about it. Quotation marks inside the tagger's sentence. Where a tag is the tagger's own sentence, a phrase in quotation marks inside it is borrowed from the paper — the tagger's wording around it is not. The paper's own sentence is the evidence quote shown on the same row, so the borrowing can be checked against the thing borrowed from without leaving this page. Nothing inside quotation marks anywhere on this page is ours: if it is quoted, it is the paper's. Whose words a concept carries. A concept entry records a term the paper itself uses, so the term, the sense given for it, and its evidence are the paper's own words, quoted where the paper fixes them. This is the opposite of a sense entry in the tag table, which is the tagger's sentence about how the paper uses a term. Where a concept entry could not be grounded in the paper's words, the entry says so on its own row.
Other fields define these differently the same word, defined another way in another discipline
machine learning 14 fields · DFT 4 fields · Species 3 fields · coordinates 3 fields · training 3 fields · SMILES 2 fields · and 4 more →

Tags on this paper

Every tag the record carries, each linking to every other paper in the corpus that carries it. A tag appears here whether or not its words occur in the article — the spans below are only the ones a sentence in this paper licenses. Of the 44 concept tags on this paper, all 44 are grounded in a quoted sentence. All carry the tagger's gloss of what the term means in this paper, which is not something the article can license.
ANI-1 data set non-equilibrium DFT total energy calculations machine learning ML potentials neural network potential ANI-1 reference data ab initio methods Kohn-Sham density-functional theory DFT chemical space conformational space configurational space GDB-11 database SMILES RDKit MMFF94 force field Gaussian 09 ωB97x hybrid-meta GGA functional 6–31 G(d) basis set normal mode sampling normal mode coordinates force constants Bernoulli distribution single point energy HDF5 Figshare data repository python dict Species python list of strings high-energy conformers atomic clashes low energy equilibrium data sets python library fitting training QMx datasets MD datasets benchmark datasets QM7 dataset coordinates energies normal mode displacement joint chemical and conformational sampling three-dimensional structure preparation geometry convergence filtering single-point convergence filtering high-energy partitioning per-molecule HDF5 organization reader-library access GDB molecules molecular conformation number of heavy atoms simulation
Method, as described: select 1–8-heavy-atom C/N/O molecules from GDB-11, generate and pre-optimize three-dimensional structures with RDKit/MMFF94, optimize them with ωB97x/6–31G(d), sample displaced normal modes, compute single-point DFT energies, exclude nonconverged calculations, and package coordinates, energies, and species in HDF5

Legend

relation to other workclaim_relationshow this paper stands to another: agrees, contradicts, extends, qualifies
claim typeclaim_typethe shape of the claim: existence, mechanism, method, critique, prediction, null
concept edgeconcept_edgesan asserted relation between two concepts, with the sentence asserting it
conceptconceptsan idea the paper works with, quoted where the paper uses it
content tagcontent_tagsa finding, definition, quantity, statistic or variable named in the paper
contributioncontributionwhat the paper claims to add
limitationlimitationswhat the paper says its evidence cannot support
mechanismmechanismshow the reported effect is said to come about
methodmethodhow the paper produced its evidence
object of studyobject_of_studywhat the paper is about — the thing it investigates
phenomenonphenomenonthe effect or regularity under study
research questionresearch_questionthe question the paper sets itself
sensesensesa distinct use of a term, quoted where the paper fixes it
unrecordable relationunhoused_relationsa relation the paper asserts that this vocabulary cannot express, with the reason

The paper, tagged

ANI-1, A data set of 20 million calculated off-equilibrium conformations for organic molecules

Abstract

concepts↗One of the grand challenges in modern theoretical chemistry is designing and implementing approximations that expedite ab initio methods without loss of accuracy. concepts↗Machine learning (ML) methods are emerging as a powerful approach to constructing various forms of transferable atomistic potentials. They have been successfully applied in a variety of applications in chemistry, biology, catalysis, and solid-state physics. unhoused_relationsHowever, these models are heavily dependent on the quality and quantity of data used in their fitting. concepts↗concepts↗Fitting highly flexible ML potentials, such as neural networks, comes at a cost: a vast amount of reference data is required to properly train these models. research_questioncontent_tagscontent_tagscontent_tagsWe address this need by providing access to a large computational DFT database, which consists of more than 20 M off equilibrium conformations for 57,462 small organic molecules. We believe it will become a new standard benchmark for comparison of current and future methods in the ML potential community.

Subject terms: Method development, Combinatorial libraries, Computational chemistry, Density functional theory

Background & Summary

Accurate descriptions of atomic and intermolecular interactions are a cornerstone of reliable computer simulations in biophysics, chemistry, and materials science. For the past 50 years we have seen tremendous progress in the development of theoretical methods and software tools aiming to describe more complex systems and allow for longer time scales. concepts↗concepts↗concept_edgesKohn-Sham density-functional theory (KS-DFT or DFT for short) has become by far the most popular electronic structure method in computational physics and chemistry 1 . DFT has found applications in many systems in organic chemistry 2 , 3 , biology 4 , catalysis 3 , 5 and solid state chemistry 6 , 7 . It is also frequently combined with molecular dynamics (AIMD) and classical force fields (quantum mechanics-molecular mechanics (QM-MM)) to describe chemical reactions in extended systems.

content_tagscontent_tagsAlthough DFT calculations have become affordable on modern supercomputers, we face a dilemma: standard computational algorithms representing the N electrons system require O (N 2 ) storage and O (N 3 ) arithmetic operations. limitationsunhoused_relationsThis O (N 3 ) complexity has become a critical bottleneck which limits capabilities to study larger realistic physical systems, as well as longer time scales relevant to actual experiment. Consequently, a lot of progress has been made in the development of atomistic potentials using machine learning (ML) 8 , 9 . The low numerical complexity and high accuracy of machine learning algorithms makes them very attractive as a pragmatic substitute for ab-initio and DFT methods. Thanks to their remarkable ability to find complex relationships among data, in many cases these ‘machine learned’ models out-perform more physically sound approximations (like force fields) and methods while also reducing the computational time required for a given application 9–15 . limitationsconcepts↗concepts↗sensesconcept_edgescontent_tags↗These models are heavily dependent on the quality and quantity of data used in their fitting, also called training. Neural networks are highly efficient and effective at modeling reference training data, due to their flexible functional form. However, this flexibility comes at a cost: a vast amount of reference data is required to properly train these models.

The Chemical Space Project 16 computationally enumerated all possible organic molecules up to a certain size, resulting in the creation of the GDB databases. content_tagsTheir latest GDB-17 database 17 contains 166.4 billion molecules of up to 17 atoms of C, N, O, S, and halogens. All molecules follow the valency rules and are filtered for unstable substructures, non-synthesizable and strained topologies. concepts↗unhoused_relationscontent_tags↗GDB molecules are stored as SMILES [ www.opensmiles.org ] strings representing the composition and connectivity of a molecule.

concepts↗concepts↗unhoused_relationsThe GDB databases were fundamental in creating the QM7 dataset 18 , one of the first benchmark datasets for training atomistic ML potentials. The QM7 dataset consists of 7,165 energy minimized (equilibrium) molecules calculated with the PBE0 functional. All structures are a small subset of GDB-13 (older GDB database of nearly 1 billion organic molecules) composed of molecules with up to 7 heavy atoms C, N, O, and S. Later, QM7 was extended to include 13 additional properties, like frontier molecular orbital energies, dipole moments, polarizability, and excitation energies 19 . The first ML model trained on QM7 used kernel ridge regression with the Coulomb matrix representation, which predicted atomization energies with a mean absolute error (MAE) of 9.9 kcal×mol −1 . This error was quickly reduced to 3.3 kcal×mol −1 (ref. 20 ) and eventually was under 1 kcal×mol −1 (ref. 21 ).

QM9, is perhaps the most well-known benchmark dataset 17 , 22 . It consists of 133,885 equilibrium organic molecules containing up to nine heavy atoms (CONF) from the GDB-17 database. In addition to energy minima it reports corresponding harmonic frequencies, dipole moments, polarizabilities, along with energies, enthalpies, and free energies of atomization. All properties were calculated at the B3LYP/6–31 G(2df,p) level of quantum chemistry. A subset of 6,095 constitutional isomers in QM9 corresponding to a brutto formula C7H10O2 was also calculated at the more accurate G4MP2 level of theory. Various molecular representations and ML methods were benchmarked against the QM9 dataset 20 , 21 , 23 , 24 . See also a recent survey of methods 23 . Later, a Message Passing Neural Network (MPNN) 10 achieved chemical accuracy in 11 out of 13 target properties in the QM9 dataset. Finally, the hierarchical interacting particle neural network (HIP-NN) 15 model of Lubbers et. al . achieved state-of-the-art accuracy of just 0.26 kcal×mol −1 MAE on total energy prediction.

limitationsconcepts↗sensesclaim_relationsA common feature of all QMx datasets is that they only explore chemical degrees of freedom by providing information about energy minimized (equilibrium) molecular configurations. In these molecules, the forces of all atoms are equal to zero. Therefore, considerable efforts were undertaken to produce off-equilibrium datasets using ab initio molecular dynamics (AIMD) simulations. The C7O2H10–17 dataset includes energies from AIMD trajectories of 113 isomers of C7O2H10 (5 k frames each). All simulations used the DFT/PBE level of theory and were carried out at 500 K. Very recently Schutt et al. 21 and Chmiela et al. 25 released MD17 dataset, a collection of eight AIMD/PBE+vdW-TS simulations for small organic molecules. Each of these consist of an MD trajectory for a single molecule extending from ~100 K to 900 K frames. concepts↗concept_edgesIn contrast to the QMx datasets, these MD datasets explore conformational space while keeping composition fixed.

concepts↗concepts↗claim_relationsWe recently introduced a neural network potential (NNP) called ANI-1, the first NNP for organic molecules shown to transfer to molecular systems well outside of its training set. concepts↗As presented, the ANI-1 potential was trained on a data set, which spans both conformational and configurational space, built from small organic molecules of up to 8-heavy atoms. We show its applicability to much larger systems, up to 50 atoms, including well known drug molecules and a random selection of molecules from the GDB-11 (refs 26 , 27 ) database with 10-heavy atoms. content_tagscontent_tagsANI-1 shows exceptional predictive power on the 10-heavy atom test set, with RMSE versus DFT relative energies as low as 0.57 kcal×mol −1 when only considering molecular conformations that are within 30 kcal×mol −1 of the energy minimum for each molecule. More recently, Gastegger et. al. 28 , showed similar results for large organic systems that were fragmented into smaller molecules and DFT data was generated on the fly for training. This was done in an active-learning fashion where the goal is to train the potential to a specific system during an MD simulation. Shortly after, Huang and Von Lilienfeld 29 used a fragmentation scheme for training an ML model to predict energies of large rigid drug molecules. Both studies back up the argument that information about the physics of large systems can be learned from data sets of small molecules.

object_of_studyconcepts↗sensesIn this data descriptor, we report a large dataset of non-equilibrium DFT total energy calculations for organic molecules. contributionmechanisms↗content_tagsIn total, we provide access to the total energies of ~20 M molecular conformations for 57,462 molecules from the GDB database 26 , 27 , which samples both chemical and conformational degrees of freedom at the same time. content_tagscontent_tagsAs the accuracy of modern ML methods for molecules in equilibrium on the QM9 benchmark achieved 1 kcal×mol −1 , ANI-1 provides 100x more data and a much more challenging task to learn. claim_typesensesTherefore, we expect it will become a new standard benchmark of comparison for current and future methods in the machine learned potential community. More importantly, it is a sound foundation for the development of future general-purpose machine learned potentials, providing an exhaustive head start on data generation, which can be augmented with future data sets covering relevant regions of chemical space.

Methods

QM calculations

concepts↗concepts↗All electronic structure calculations are carried out with the ωB97x (ref. 30 ) density functional and the 6–31 G(d) basis set 31 in the Gaussian 09 (ref. 32 ) electronic structure package. concepts↗concepts↗concept_edgesωB97x is a hybrid-meta GGA functional 30 , which has been shown to be chemically accurate compared to high-level CCSD(T) calculations 33–37 .

Molecular geometry generation

concepts↗The GDB-11 database 26 , 27 provides an exhaustive search of stable and chemically viable molecules, supplied in the SMILES [ www.opensmiles.org ] string format, containing C, N, O, and F atoms with up to 11 of these ‘heavy’ atoms. Hydrogen atoms are added through the RDKit cheminformatics software package [ www.rdkit.org ] to make molecular structures that are charge neutral and have a singlet electronic ground state. method↗concepts↗The ANI-1 data set presented here is built from an exhaustive sampling of a subset of the GDB-11 database containing molecules with between 1 and 8 heavy atoms and limiting the atomic species to C, N, and O. content_tagsThis leaves a subset of 57,947 starting molecules. All molecules are neutral and with a singlet electronic ground state. The conformation generation process is carried out in five steps starting with these 57,947 molecules. The steps are listed below and qualitatively depicted in Fig. 1 .

Schematic representation of data generation. Scheme for generating non-equilibrium conformations of 57,462 molecules from the GDB-11 database. The goal with this scheme is to generate a ‘window’ of the potential surface around each optimized equilibrium structure.

mechanisms↗concepts↗Smiles strings from the GDB-11 subset described above are used to generate 3D conformations using RDKit. Also with RDKit, all structures are saturated with hydrogens such that each has charge 0 and multiplicity 1. concepts↗The 3D structures are then pre-optimized to a stationary point using the MMFF94 force field 38 as implemented in RDKit.

At the chosen DFT or ab-initio level of theory, geometries are optimized until energy minima convergence. Optimization is carried out using Gaussian 09’s default method and convergence criteria. Obtained geometries correspond to the first stationary point reached on the potential surface and correspond to some local minima or in a rare case to a saddle point. mechanisms↗unhoused_relationsIf convergence fails, the structure is not included in the data set. content_tagscontent_tagsAt this step, 485 (0.84% of total) molecules failed to converge during the structural optimization. The final data set is built from these 57,462 equilibrium geometries. concepts↗content_tags↗Finally, for each of the 57,462 structurally optimized molecules, a normal mode calculation is performed in the Gaussian 09 package to obtained normal mode coordinates and their associated force constants. This is accomplished using the UltraFine DFT grid option with the ωB97x density functional.

Normal mode sampling (NMS)

To carry out normal mode sampling on an energy minimized molecule of N a atoms, first a set of N f normal mode coordinates, Q = { q 1 , q 2 , q 3 , … q N f } , is computed at the desired ab-initio level of theory, where N f =3 N a −5 for linear molecules and N f =3 N a −6 for all others. The corresponding force constants, K = { K 1 , K 2 , K 3 , ⋯ , K N f } , are obtained alongside Q . Then a set of N f uniformly distributed pseudo-random numbers, c i , are generated such that ∑ i N f c i is in the range [0,1]. Next, a displacement, R i , for each normal mode coordinate is computed by setting a harmonic potential equal to the c i scaled average energy of the system of particles at some temperature, T. Solving for the displacement gives,

where k b is Boltzmann’s constant. mechanisms↗concepts↗sensescontent_tagsThe sign of R i is determined randomly from a Bernoulli distribution where P =0.5 to ensure that both sides of the harmonic potential are sampled equally. concepts↗content_tags↗Each R i is used to scale the normalized normal mode coordinates by q i R = R i q i . content_tags↗Next, a new conformation of the molecule is generated by displacing the structurally optimized coordinates by Q R , the superposition of all q i R . concepts↗content_tags↗Finally, a single point energy at the desired level of theory is calculated using the newly displaced coordinates as input.

N data points (new conformations) are generated, representing a window of the potential surface. content_tags↗N is calculated by S × K where S is an empirically chosen value (See Table 1 ) based on the number of heavy atoms in each molecule and K is the number of degrees of freedom of the molecule. concepts↗content_tags↗The total energy, atomic symbols, and cartesian coordinates of the structure are stored as described in the Data Format section.

List of information and parameters used to generate the ANI-1 data set.

Data Records

concepts↗concepts↗unhoused_relationscontent_tagsThe data set is provided in an HDF5 based file in a Figshare data repository (Data Citation 1). A GitHub repository containing a README file with technical usage details and examples of how to access the data set is supplied online ( https://github.com/isayev/ANI1_dataset ).

File format

mechanisms↗Data is stored per molecule as described in Fig. 2 . concepts↗Data for each X molecule is stored in a python dict type containing all conformer data. concepts↗content_tags↗The keys shown in Fig. 2 : coordinates, energies, and species give access to containers of the type shown, and containing data described by the key. concepts↗concepts↗sensesconcept_edgescontent_tags↗content_tags↗Species is a python list of strings containing the atomic symbol of each atom and its order corresponds correctly to dimension 1 of the coordinates numpy array. Appending ‘HE’ to the end of the coordinates and energies keys will yield high energy structures as described in the technical validation section.

Data Structure. Description of the containers stored in the dictionary returned by iterating through molecules stored in the HDF5 file format. The ‘coordinates’ key gives access to a 3D array containing each conformer of the molecule in cartesian coordinates, while the ‘energies’ key gives the 1D array of energies for the conformers. The first dimension of each ‘coordinates’ and ‘energies’ array maps correctly to the corresponding structure. The ‘species’ key contains the atomic symbol of the atoms and is ordered to correspond to the correct atoms in the second dimension of the array returned by the ‘coordinates’ key. Other keys in the returned dictionary are: ‘coordinatesHE’, ‘energiesHE’, and ‘smiles’ for the high energy coordinates, high energy energies and SMILES string, respectively.

Technical Validation

limitationsconcepts↗concept_edgesSince normal mode sampling is used to generate the non-equilibrium structures, high-energy conformers exist in the data set. concepts↗These high energy conformations occur where the harmonic approximation of normal modes fail in anharmonic regions of a potential, and are caused by atomic clashes or other highly unfavorable molecular conformations. content_tagsThe distribution shown in Fig. 3b visualizes the energies in the dataset, which contains structures with energies as high as 15 Ha. mechanisms↗content_tagscontent_tagsFor this reason, energies greater than 275 kcal×mol −1 higher than the lowest energy conformer were not included into the training set of the ANI-1 potential. content_tagscontent_tagsThis removed 2,630,435 (10.7% of the original total) structures yielding 22,057,374 structures. content_tagsRegions this high in energy are generally not considered in bio-chemical research. However, this data might be useful for some purposes. Therefore, we include both the high-energy and the low energy datasets as described in the data description section. Figure 3c shows the new distribution of energies which are never larger than 0 Ha in total energy minus the sum of atomic contributions to the total energy.

Dataset energy distribution. ( a ) The distribution of total energies divided by the number of electrons from normal mode sampling conducted on each sub set (04 through 08) of GDB-11. Each distribution is scaled to have equal area. ( b ) Distribution of atomization energies from the completed data set with the inset showing a long tail reaching greater than 12 Ha. ( c ) Distribution of atomization energies after truncating any energies over 275 kcal×mol −1 from each molecule’s minimum energy.

limitationsDuring the structural optimization phase, we do not distinguish between optimized structures that might land at a saddle point in the potential surface and those that land at some structural minima. Given the goal of sampling conformational space, the fact that some structures might land at off equilibrium geometries (saddle points) could in fact help in using this data to fit potential surfaces, as it will help to cover regions of conformational space not covered by equilibrium molecule normal mode sampling. limitationsHowever, if the optimization fails to converge to a stationary point, as 485 molecules did, then these structures were not included in the training set, as the validity of their configuration could not immediately be confirmed. However, given the vast number of structures in the data set, it is likely any interaction found in these 485 molecules can be found elsewhere in the data set.

A similar process of not including information for unconverged calculations is used in the generation of the total energies. limitationsFor certain highly elongated bonds the molecular orbital optimization process, the self-consistent field procedure used in obtaining the total energy of the conformation, can fail to converge to a solution when two orbitals are too close in energy. mechanisms↗unhoused_relationsFor this reason, if a structure’s single point energy calculation failed to converge, then this data is not included in the data set.

phenomenonconcepts↗claim_relationsThe primary concept of including non-equilibrium data is to sample regions of chemical space that would be sparsely covered in equilibrium only data sets. Figure 3a provides validation of energy sampling by showing the distribution of total energies divided by the total number of electrons for each molecule in the GDB subsets from 4 to 8-heavy atoms. concepts↗sensesFigure 3b,c show the distribution of total energies minus the sum of all individual atomic energies (tabulated in Supplementary Information Table 1 ) for the full and ‘low energy’ (less than 275 kcal×mol −1 from the minimum energy) data sets, respectively.

Further validation of the of non-equilibrium sampling is to show the data set covers a large domain of the chemical degrees of freedom in conformational space. Figure 4 contains five panels representing the distribution of atomic distances in the resulting non-equilibrium data set (blue line) compared with a data set of equilibrium only conformations (red) of the same molecule. concepts↗concepts↗unhoused_relationssensescontent_tagsAs expected, the normal mode sampling method used to generate non-equilibrium conformations visits areas of conformational space not covered by equilibrium only data. A similar plot, Supplementary Information Fig. 1 shows distance distributions for the remaining atomic pairs. Figure 5 shows distributions involving the angles in the data sets, and tell a similar story in terms of coverage in conformational space for three body interactions. concepts↗concept_edgescontent_tagsThe blue background density plot shows that the ANI-1 data set covers far more angular space than the equilibrium data sets (red and orange). The remaining plots are included in the Supplementary Information Figs 2–4 .

Distance distribution for dataset. Distribution of atomic distances in the subset of the data set constructed from the molecules containing between 4 and 8 heavy atoms (GDB-04 to 08) of C, N, and O. The y-axis is the base 10 logarithm of the count of distances in each bin, normalized over the full domain so that the two sets can be compared. The x-axis represents the atomic distance (r) divided by the single bond equilibrium distance ( r 0 ) for the smallest possible molecule containing a single bond of the type shown, as calculated using the ωB97x density functional with the 6–31 g(d) basis set. The red histogram shows the full distribution of distances for a data set containing only equilibrium distances. The blue line shows the distribution of our non-equilibrium data set, with distances randomly sub sampled at a rate of 1%. As the figure shows, even 1% of the non-equilibrium data set covers vast areas of atomic distance space where the equilibrium data set fails to sample.

Angular distribution. Figure shows distributions involving the angles in the data sets, and tells a similar story in terms of coverage in conformational space for three body interactions. The blue background density plot shows that the ANI-1 data set better covers angle space than the equilibrium data sets (red and orange). The remaining figures for the angular distributions are included in the Supplementary Information .

Usage Notes

mechanisms↗concepts↗To ensure that all readers have easy access to the ANI-1 data set, we have developed a python library with an easy to use interface for extracting the data. Examples uses of this library are included in the ‘readers’ folder.

Additional information

How to cite this article: Smith, J. S. et al. ANI-1, A data set of 20 million calculated off-equilibrium conformations for organic molecules. Sci. Data 4:170193 doi: 10.1038/sdata.2017.193 (2017).

Publisher’s note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Material

Acknowledgments

J.S.S. acknowledges the University of Florida for funding through the Graduate School Fellowship (GSF). A.E.R. thanks NIH award GM110077. O.I. acknowledges support from DOD-ONR (N00014-16-1-2311) and Eshelman Institute for Innovation award. Part of this research was performed while O.I. was visiting the Institute for Pure and Applied Mathematics (IPAM), which is supported by the National Science Foundation (NSF). The authors acknowledge Extreme Science and Engineering Discovery Environment (XSEDE) award DMR110088, which is supported by National Science Foundation grant number ACI-1053575. We gratefully acknowledge the support of the U.S. Department of Energy through the LANL/LDRD Program for this work. We also acknowledge Nicholas Lubbers and Roman Zubatyuk for stimulating discussions and technical help with data organization.

Footnotes

Data Citations

References

Associated Data

Data Citations

Supplementary Materials

Every tag on this page

claim_relations 3 · claim_type 1 · concept_edges 7 · concepts 44 · content_tags 35 · contribution 1 · limitations 7 · mechanisms 8 · method 1 · object_of_study 1 · phenomenon 1 · research_question 1 · senses 8 · unhoused_relations 8

Relations on this page carry a quotation for the claim. Where a relation was separately evidenced — and so tested against the vocabulary for its type — the row says so; 7 of 7 relations here carry that mark. The rest are the tagger’s reading of the sentence quoted beside them.

This record holds no controls. A layer is empty when the tagger recorded nothing in it — which is not the same as the paper having nothing to record.

facettagfieldconf.whose words
object_of_studythe ANI-1 dataset of off-equilibrium molecular conformations, DFT energies, and associated molecular descriptorsobject_of_studythe paper's words
research_questionhow to construct and validate a large reference dataset spanning chemical and conformational space for training transferable machine-learned atomistic potentialsresearch_questionthe paper's words
contributiona reusable HDF5 dataset of roughly 20 million off-equilibrium conformations for 57,462 organic molecules, generated with normal-mode sampling and DFT calculations and accompanied by a reader librarycontributionthe paper's words
claim_typeexistenceclaim_typethe paper's words
phenomenoncoverage of off-equilibrium molecular energy landscapes across both chemical composition and molecular conformationphenomenonthe paper's words
methodsimulation
also in: simulation
methodthe paper's words
mechanismsRandom signed displacements scale normal-mode coordinates, whose superposition produces a new molecular conformation for a single-point energy calculation.mechanisms[0]0.95the tagger's words
mechanismsThe data combine many molecular compositions with non-equilibrium conformations so chemical and conformational degrees of freedom are sampled together.mechanisms[1]0.95the tagger's words
mechanismsRDKit generates hydrogen-saturated neutral singlet structures from SMILES and MMFF94 pre-optimizes them to a stationary point.mechanisms[2]0.95the tagger's words
mechanismsOnly structures whose DFT geometry optimization converges are retained for the final sampling set.mechanisms[3]0.99the tagger's words
mechanismsConformations whose self-consistent-field energy calculation fails to converge are omitted.mechanisms[4]0.99the tagger's words
mechanismsA 275 kcal/mol relative-energy cutoff separates the lower-energy ANI-1 training subset while the complete high-energy data remain available.mechanisms[5]0.95the tagger's words
mechanismsEach molecule is exposed as a Python dictionary whose coordinates, energies, and species keys access conformer-level containers.mechanisms[6]0.95the tagger's words
mechanismsA Python reader library supplies an interface and examples for extracting ANI-1 data.mechanisms[7]0.99the tagger's words
limitationsCubic arithmetic complexity is a bottleneck for conventional electronic-structure calculations on larger systems and longer experimental time scales.limitations[0]0.95the tagger's words
limitationsFlexible machine-learning potentials require a vast amount of high-quality reference data.limitations[1]0.95the tagger's words
limitationsEarlier QMx benchmark datasets cover equilibrium molecular configurations and therefore do not span off-equilibrium conformational degrees of freedom.limitations[2]0.95the tagger's words
limitationsThe structural-optimization stage does not distinguish local minima from saddle points.limitations[3]0.99the tagger's words
limitationsThe 485 starting molecules whose structural optimization failed were excluded because their configurations could not immediately be validated.limitations[4]0.95the tagger's words
limitationsSome highly elongated bonds cause self-consistent-field optimization failure, so those single-point energies are absent.limitations[5]0.95the tagger's words
limitationsNormal-mode sampling produces very high-energy conformers where the harmonic approximation fails in anharmonic regions.limitations[6]0.95the tagger's words
conceptsANI-1 data setconcepts[0]0.99the paper's words
conceptsnon-equilibrium DFT total energy calculationsconcepts[1]0.99the paper's words
conceptsmachine learningconcepts[2]0.99the paper's words
conceptsML potentials
also in: ML potentials
concepts[3]0.99the paper's words
conceptsneural network potentialconcepts[4]0.99the paper's words
conceptsANI-1
also in: ANI-1
concepts[5]0.99the paper's words
conceptsreference dataconcepts[6]0.95the paper's words
conceptsab initio methodsconcepts[7]0.95the paper's words
conceptsKohn-Sham density-functional theoryconcepts[8]0.99the paper's words
conceptsDFT
also in: DFT
concepts[9]0.99the paper's words
conceptschemical spaceconcepts[10]0.99the paper's words
conceptsconformational spaceconcepts[11]0.99the paper's words
conceptsconfigurational spaceconcepts[12]0.95the paper's words
conceptsGDB-11 databaseconcepts[13]0.99the paper's words
conceptsSMILES
also in: SMILES
concepts[14]0.99the paper's words
conceptsRDKit
also in: RDKit
concepts[15]0.99the paper's words
conceptsMMFF94 force fieldconcepts[16]0.99the paper's words
conceptsGaussian 09
also in: Gaussian 09
concepts[17]0.99the paper's words
conceptsωB97x
also in: ωB97x
concepts[18]0.99the paper's words
conceptshybrid-meta GGA functionalconcepts[19]0.99the paper's words
concepts6–31 G(d) basis setconcepts[20]0.99the paper's words
conceptsnormal mode samplingconcepts[21]0.99the paper's words
conceptsnormal mode coordinatesconcepts[22]0.99the paper's words
conceptsforce constantsconcepts[23]0.95the paper's words
conceptsBernoulli distributionconcepts[24]0.99the paper's words
conceptssingle point energyconcepts[25]0.99the paper's words
conceptsHDF5
also in: HDF5
concepts[26]0.99the paper's words
conceptsFigshare data repositoryconcepts[27]0.99the paper's words
conceptspython dict
also in: python dict
concepts[28]0.95the paper's words
conceptsSpecies
also in: Species
concepts[29]0.99the paper's words
conceptspython list of stringsconcepts[30]0.99the paper's words
conceptshigh-energy conformersconcepts[31]0.99the paper's words
conceptsatomic clashesconcepts[32]0.99the paper's words
conceptslow energy
also in: low energy
concepts[33]0.99the paper's words
conceptsequilibrium data setsconcepts[34]0.99the paper's words
conceptspython libraryconcepts[35]0.99the paper's words
conceptsfitting
also in: fitting
concepts[36]0.99the paper's words
conceptstraining
also in: training
concepts[37]0.99the paper's words
conceptsQMx datasets
also in: QMx datasets
concepts[38]0.95the paper's words
conceptsMD datasets
also in: MD datasets
concepts[39]0.95the paper's words
conceptsbenchmark datasetsconcepts[40]0.95the paper's words
conceptsQM7 dataset
also in: QM7 dataset
concepts[41]0.95the paper's words
conceptscoordinates
also in: coordinates
concepts[42]0.99the paper's words
conceptsenergies
also in: energies
concepts[43]0.99the paper's words
unhoused_relationsquality and quantity of data → fitted machine-learning models
Fitted models depend heavily on the quality and quantity of their training data. — the relation stated: The models depend heavily on training-data quality and quantity.
unhoused_relations[0]0.99the tagger's words
unhoused_relationsO (N 3 ) complexity → capabilities to study larger systems and longer time scales
Cubic computational complexity limits study of larger systems and longer relevant time scales. — the relation stated: Cubic complexity limits the tractable system size and time scale.
unhoused_relations[1]0.99the tagger's words
unhoused_relationsGDB molecules → SMILES strings
GDB molecules are stored as SMILES strings representing composition and connectivity. — the relation stated: GDB molecules are serialized as SMILES strings.
unhoused_relations[2]0.99the tagger's words
unhoused_relationsnormal mode sampling → areas of conformational space not covered by equilibrium only data
Normal-mode sampling visits conformational regions absent from equilibrium-only data. — the relation stated: Normal-mode sampling reaches conformational regions not covered by equilibrium-only data.
unhoused_relations[3]0.95the tagger's words
unhoused_relationsfailed structural optimization → inclusion in the data set
A structure is excluded when its geometry optimization does not converge. — the relation stated: Failed convergence excludes the structure from the dataset.
unhoused_relations[4]0.99the tagger's words
unhoused_relationssingle point energy calculation failed to converge → inclusion in the data set
A conformation is excluded when its single-point energy does not converge. — the relation stated: Failed single-point convergence excludes the data from the dataset.
unhoused_relations[5]0.99the tagger's words
unhoused_relationsHDF5 based file → Figshare data repository
The released HDF5 dataset file is hosted in Figshare. — the relation stated: The HDF5 file is distributed through Figshare.
unhoused_relations[6]0.99the tagger's words
unhoused_relationsQM7 dataset → benchmark datasets
QM7 is described as one of the first benchmark datasets for training atomistic machine-learning potentials. — the relation stated: QM7 is one of the first benchmark datasets for atomistic ML potentials.
unhoused_relations[7]0.99the tagger's words
sensesANI-1 as an off-equilibrium DFT reference datasetsenses[0]0.99the tagger's words
sensesANI-1 as a benchmark foundation for machine-learned potentialssenses[1]0.9the tagger's words
senseschemical-space coverage as molecular-composition diversitysenses[2]0.95the tagger's words
sensesconformational-space coverage as variation away from energy-minimized geometrysenses[3]0.95the tagger's words
sensesnormal mode sampling as randomized displacement of optimized geometriessenses[4]0.95the tagger's words
senseslow-energy ANI-1 as the under-275-kcal/mol partitionsenses[5]0.99the tagger's words
sensesSpecies as the ordered atomic-symbol list aligned to coordinatessenses[6]0.99the tagger's words
sensesfitting as machine-learning trainingsenses[7]0.99the tagger's words
claim_relationsbuilds-on → QM7 and QM9 equilibrium quantum-chemistry benchmark datasets
note: ANI-1 extends benchmark resources beyond equilibrium configurations.
claim_relations[0]0.95the tagger's words
claim_relationsbuilds-on → the previously introduced ANI-1 neural-network potential
note: This data descriptor releases and documents the training-data foundation underlying the prior potential.
claim_relations[1]0.95the tagger's words
claim_relationscontrasts-with → equilibrium-only molecular benchmark datasets
note: The dataset's motivating contrast is explicit.
claim_relations[2]0.95the tagger's words
concept_edgesKohn-Sham density-functional theory —same-as→ DFT
the relation itself is separately evidenced · strength of the relation: 1 · note: The abbreviation equivalence is explicit and categorical. · the relation stated: DFT is the stated short name for Kohn-Sham density-functional theory.
concept_edges[0]0.99the tagger's words
concept_edgesfitting —same-as→ training
the relation itself is separately evidenced · strength of the relation: 1 · note: The source gives a categorical naming equivalence. · the relation stated: Fitting is also called training in this modeling context.
also in: fitting · training
concept_edges[1]0.99the tagger's words
concept_edgesωB97x —is-a→ hybrid-meta GGA functional
the relation itself is separately evidenced · strength of the relation: 1 · note: The source makes the class relation explicitly. · the relation stated: ωB97x is a hybrid-meta GGA functional.
concept_edges[2]0.99the tagger's words
concept_edgesSpecies —is-a→ python list of strings
the relation itself is separately evidenced · strength of the relation: 1 · note: This is an explicit categorical representation statement. · the relation stated: The Species value is a Python list of strings.
concept_edges[3]0.99the tagger's words
concept_edgesMD datasets —contrasted-with→ QMx datasets
the relation itself is separately evidenced · strength of the relation: 0.95 · note: In contrast to is an explicit contrast cue. · the relation stated: MD datasets are contrasted with QMx datasets by varying conformation while holding composition fixed.
concept_edges[4]0.99the tagger's words
concept_edgesANI-1 data set —contrasted-with→ equilibrium data sets
the relation itself is separately evidenced · strength of the relation: 0.95 · note: Far more than is a direct gradable comparison. · the relation stated: ANI-1 is contrasted with equilibrium datasets by its much broader angular-space coverage.
concept_edges[5]0.95the tagger's words
concept_edgesatomic clashes —causes→ high-energy conformers
the relation itself is separately evidenced · strength of the relation: 0.95 · note: Caused by is an explicit positive causal cue; the nonexclusive alternative keeps weight below one. · the relation stated: Atomic clashes cause high-energy sampled conformations.
concept_edges[6]0.99the tagger's words
content_tagsANI-1 provides more than 20 million off-equilibrium conformations for 57,462 small organic molecules.content_tags.findings[0]0.99the tagger's words
content_tagsANI-1 samples chemical and conformational degrees of freedom simultaneously.content_tags.findings[1]0.99the tagger's words
content_tagsANI-1 contains roughly one hundred times as much data as QM9 and presents a harder learning task.content_tags.findings[2]0.95the tagger's words
content_tagsOf 57,947 starting molecules, 485 failed structural optimization and the final dataset uses 57,462 equilibrium geometries.content_tags.findings[3]0.99the tagger's words
content_tagsApplying the high-energy exclusion removed 2,630,435 structures and left 22,057,374 in the lower-energy training set.content_tags.findings[4]0.99the tagger's words
content_tagsThe released data include both high- and low-energy partitions despite the lower-energy cutoff used for ANI-1 potential training.content_tags.findings[5]0.95the tagger's words
content_tagsNormal-mode sampling reaches conformational regions absent from an equilibrium-only dataset of the same molecules.content_tags.findings[6]0.95the tagger's words
content_tagsANI-1 covers substantially more angular space than equilibrium datasets.content_tags.findings[7]0.95the tagger's words
content_tagsThe prior ANI-1 neural-network potential transferred to larger systems and reached relative-energy RMSE as low as 0.57 kcal/mol on a 10-heavy-atom test set under the stated energy window.content_tags.findings[8]0.9the tagger's words
content_tagsThe dataset is distributed as HDF5 through Figshare with usage documentation and examples on GitHub.content_tags.findings[9]0.99the tagger's words
content_tagstraining a model on reference data
also in: fitting
content_tags.definitions[0]0.99the tagger's words
content_tagsenumerated molecules stored as SMILES strings representing molecular composition and connectivity
also in: GDB molecules
content_tags.definitions[1]0.95the tagger's words
content_tagsa Python list of atomic-symbol strings ordered to match the first coordinate-array dimension
also in: Species
content_tags.definitions[3]0.99the tagger's words
content_tagsmore than 20 M — off-equilibrium conformations in the computational DFT databasecontent_tags.quantities[0]0.99the tagger's words
content_tags57,462 — small organic molecules represented in the released ANI-1 datasetcontent_tags.quantities[1]0.99the tagger's words
content_tagsO (N 2 ) — storage scaling of standard N-electron computational algorithmscontent_tags.quantities[2]0.95the tagger's words
content_tagsO (N 3 ) — arithmetic scaling of standard N-electron computational algorithmscontent_tags.quantities[3]0.95the tagger's words
content_tags166.4 billion — molecules reported in GDB-17content_tags.quantities[4]0.95the tagger's words
content_tags57,947 — starting molecules after selecting 1–8-heavy-atom C/N/O structures from GDB-11content_tags.quantities[5]0.99the tagger's words
content_tags485 (0.84% of total) — molecules failing convergence during structural optimizationcontent_tags.quantities[6]0.99the tagger's words
content_tagsP =0.5 — Bernoulli probability used to sample both signs of each normal-mode displacementcontent_tags.quantities[7]0.99the tagger's words
content_tags275 kcal×mol −1 — relative-energy threshold for the low-energy ANI-1 training partitioncontent_tags.quantities[8]0.99the tagger's words
content_tags2,630,435 (10.7% of the original total) — high-energy structures removed from the ANI-1 potential training setcontent_tags.quantities[9]0.99the tagger's words
content_tags22,057,374 — structures remaining in the lower-energy training partitioncontent_tags.quantities[10]0.99the tagger's words
content_tags15 Ha — upper energy reported for structures in the complete sampled distributioncontent_tags.quantities[11]0.95the tagger's words
content_tags100x more data — ANI-1 size relative to the QM9 benchmark as characterized by the authorscontent_tags.quantities[12]0.95the tagger's words
content_tags0.57 kcal×mol −1 — lowest reported ANI-1-potential RMSE versus DFT relative energies on the restricted 10-heavy-atom test conformationscontent_tags.quantities[13]0.9the tagger's words
content_tagsmolecular conformationcontent_tags.variables[1]0.95the tagger's words
content_tagsnormal mode coordinatescontent_tags.variables[2]0.95the tagger's words
content_tagsforce constantscontent_tags.variables[3]0.9the tagger's words
content_tagssingle point energycontent_tags.variables[4]0.95the tagger's words
content_tagsnumber of heavy atomscontent_tags.variables[5]0.95the tagger's words
content_tagscoordinates
also in: coordinates
content_tags.variables[6]0.95the tagger's words
content_tagsenergies
also in: energies
content_tags.variables[7]0.95the tagger's words
content_tagsspecies
also in: species
content_tags.variables[8]0.99the tagger's words

Interpretive choices

Tags that could not be placed in the text:

What was deliberately not tagged, and why, is recorded per paper by the tagger. Where a paper has no such record, this page says so rather than implying that nothing was declined.

© 2026 Shir Sivroni · ORCID 0009-0000-2597-2824