An official website of the United States government
In this work, forcefield flexibility parameters were constructed and validated for more than 100 metal-organic frameworks (MOFs). We used atom typing to identify bond types, angle types, and dihedral types associated with bond stretches, angle bends, dihedral torsions, and other flexibility interactions. Our work used Manz’s angle-bending and dihedral-torsion model potentials. For a crystal structure containing Natoms in its unit cell, the number of independent flexibility interactions is 3(Natoms – 1). Because the number of bonds, angles, and dihedrals is normally much larger than 3(Natoms – 1), these internal coordinates are redundant. To reduce (but not eliminate) this redundancy, our protocol prunes dihedral types in a way that preserves symmetry equivalency. Next, each dihedral type is classified as non-rotatable, hindered, rotatable, or linear. We introduce a smart selection method that identifies which particular torsion modes are important for each rotatable dihedral type. Then, we computed the force constants for all flexibility interactions together via LASSO regression (i.e., regularized linear least-squares fitting) of the training dataset. LASSO automatically identifies and removes unimportant forcefield interactions. For each MOF, the reference dataset was quantum-mechanically-computed in VASP via DFT with dispersion and included: (i) finite-displacement calculations along every independent atom translation mode, (ii) geometries randomly sampled via ab initio molecular dynamics (AIMD), (iii) the optimized ground-state geometry using experimental lattice parameters, and (iv) rigid torsion scans for each rotatable dihedral type. After training, the flexibility model was validated across geometries that were not part of the training dataset. For each MOF, we computed the goodness of fit (R-squared value) and the root-mean-squared error (RMSE) separately for the training and validation datasets. We compared flexibility models with and without bond-bond cross terms. Even without cross terms, the model yielded R-squared values of 0.910 (avg across all MOFs) ± 0.018 (st. dev.) for atom-in-material forces in the validation datasets. Our SAVESTEPS protocol should find widespread applications to parameterize flexible forcefields for material datasets. We performed molecular dynamics simulations using these flexibility parameters to compute heat capacities and thermal expansion coefficients for two MOFs.
more »
« less
This content will become publicly available on March 6, 2026
Dihedral–torsion model potentials that include angle-damping factors
This article introduces the new angle-damped dihedral torsion (ADDT), angle-damped linear dihedral (ADLD), angle-damped cosine only (ADCO), and constant amplitude dihedral torsion (CADT) model potentials.
more »
« less
- Award ID(s):
- 1555376
- PAR ID:
- 10576314
- Publisher / Repository:
- Royal Society of Chemistry
- Date Published:
- Journal Name:
- RSC Advances
- Volume:
- 15
- Issue:
- 10
- ISSN:
- 2046-2069
- Page Range / eLocation ID:
- 7257 to 7306
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
In the title compound, C10H12N2O4, the four substituents lie out of the phenyl plane by varying degrees. The methyl C atom lies 0.019 (3) A ˚ out of plane, while the methoxy O and C atoms lie 0.067 (2) and 0.042 (3) A ˚ out of plane, respectively, with the C—C—O—C torsion angle being 3.3 (2). The plane of the nitro group is twisted out of the phenyl plane, forming a dihedral angle of 12.03 (9) with it. The acetamide substituent is twisted considerably more out of the phenyl plane, forming a dihedral angle of 47.24 (6) with it. In the extended structure, the acetamide NH group donates a hydrogen bond to an acetamide carbonyl O atom, thereby forming chains propagating in the [010] direction.more » « less
-
Two bis-carbamoylmethylphosphine oxide compounds, namely {[(3-{[2-(diphenylphosphinoyl)ethanamido]methyl}benzyl)carbamoyl]methyl}diphenylphosphine oxide, C 36 H 34 N 2 O 4 P 2 , (I), and diethyl [({2-[2-(diethoxyphosphinoyl)ethanamido]ethyl}carbamoyl)methyl]phosphonate, C 14 H 30 N 2 O 8 P 2 , (II), were synthesized via nucleophilic acyl substitution reactions between an ester and a primary amine. Hydrogen-bonding interactions are present in both crystals, but these interactions are intramolecular in the case of compound (I) and intermolecular in compound (II). Intramolecular π–π stacking interactions are also present in the crystal of compound (I) with a centroid–centroid distance of 3.9479 (12) Å and a dihedral angle of 9.56 (12)°. Intermolecular C—H...π interactions [C...centroid distance of 3.622 (2) Å, C—H...centroid angle of 146°] give rise to supramolecular sheets that lie in the ab plane. Key geometric features for compound (I) involve a nearly planar, trans- amide group with a C—N—C—C torsion angle of 169.12 (17)°, and a torsion angle of −108.39 (15)° between the phosphine oxide phosphorus atom and the amide nitrogen atom. For compound (II), the electron density corresponding to the phosphoryl group was disordered, and was modeled as two parts with a 0.7387 (19):0.2613 (19) occupancy ratio. Compound (II) also boasts a trans -amide group that approaches planarity with a C—N—C—C torsion angle of −176.50 (16)°. The hydrogen bonds in this structure are intermolecular, with a D ... A distance of 2.883 (2) Å and a D —H... A angle of 175.0 (18)° between the amide hydrogen atom and the P=O oxygen atom. These non-covalent interactions create ribbons that run along the b -axis direction.more » « less
-
In the title compound, C16H16N2O3, the phenyl groups are twisted away from coplanarity with the ether linkage, forming a dihedral angle of 59.49 (4) with each other. The ether oxygen atom lies somewhat out of both phenyl planes, by 0.066 (2) and 0.097 (2) A ˚ . The acetamide substituents have quite different conformations with respect to the phenyl groups on either side of the molecule. On one side, the C—C—N—C torsion angle is 21.0 (2), while on the other side it is 76.4 (2). In the crystal, the acetamide N—H groups form intermolecular N—H O hydrogen bonds to acetamide O atom, with both NH groups donating to the same molecule. Thus, ladder-like chains exist in the [101] direction. One of the methyl groups has its H atoms disordered into two orientations, and the crystal chosen for data collection was found to be twinned.more » « less
-
The structure of zymonic acid (systematic name: 4-hydroxy-2-methyl-5-oxo-2,5-dihydrofuran-2-carboxylic acid), C 6 H 6 O 5 , which had previously eluded crystallographic determination, is presented here for the first time. It forms by intramolecular condensation of parapyruvic acid, which is the product of aldol condensation of pyruvic acid. A redetermination of the crystal structure of pyruvic acid (systematic name: 2-oxopropanoic acid), C 3 H 4 O 3 , at low temperature (90 K) and with increased precision, is also presented [for the previous structure, see: Harata et al. (1977). Acta Cryst. B 33 , 210–212]. In zymonic acid, the hydroxylactone ring is close to planar (r.m.s. deviation = 0.0108 Å) and the dihedral angle between the ring and the plane formed by the bonds of the methyl and carboxylic acid carbon atoms to the ring is 88.68 (7)°. The torsion angle of the carboxylic acid group relative to the ring is 12.04 (16)°. The pyruvic acid molecule is almost planar, having a dihedral angle between the carboxylic acid and methyl-ketone groups of 3.95 (6)°. Intermolecular interactions in both crystal structures are dominated by hydrogen bonding. The common R 2 2 (8) hydrogen-bonding motif links carboxylic acid groups on adjacent molecules in both structures. In zymonic acid, this results in dimers about a crystallographic twofold of space group C 2/ c , which forces the carboxylic acid group to be disordered exactly 50:50, which scrambles the carbonyl and hydroxyl groups and gives an apparent equalization of the C—O bond lengths [1.2568 (16) and 1.2602 (16) Å]. The other hydrogen bonds in zymonic acid (O—H...O and weak C—H...O), link molecules across a 2 1 -screw axis, and generate an R 2 2 (9) motif. These hydrogen-bonding interactions propagate to form extended pleated sheets in the ab plane. Stacking of these zigzag sheets along c involves only van der Waals contacts. In pyruvic acid, inversion-related molecules are linked into R 2 2 (8) dimers, with van der Waals interactions between dimers as the only other intermolecular contacts.more » « less
