skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


This content will become publicly available on November 12, 2026

Title: Gas-Phase-like Inertial Rotation of a Platonic Solid in Reticular Amphidynamic Crystals of CUB-5: Molecular Dynamics and Experiment
Not AvailableAmphidynamic crystals are a type of condensed matter that blends two extremes of the dynamic spectrum: rigid components forming a static lattice and rapidly moving parts. Among them, ordered rotor arrays within metal-organic frameworks (MOFs) constitute a promising platform to explore unchartered territories, such as gas phase-like dynamics in the crystalline state. Through quantum mechanical (QM) calculations and molecular dynamics (MD) simulations we verified that nearly barrierless cubane rotators in CUB-5 display rotational dynamics that transitions from continuous or inertial at high tempera-ture, to chaotic behavior, and ultimately to discrete jumps, as the temperature decreases from room temperature down to cry-ogenic conditions. 1H NMR spin-lattice (T1) relaxation measurements corroborate our theoretical predictions, with experi-mental rotational activation energy of 0.17 kcal/mol and an attempt frequency of 1.03×1012 s-1 that compare well with calcu-lated values of 0.15 kcal/mol and 0.38×1012 s-1, respectively.  more » « less
Award ID(s):
2203519
PAR ID:
10648292
Author(s) / Creator(s):
; ; ; ; ; ; ;
Publisher / Repository:
American Chemical Society
Date Published:
Journal Name:
Journal of the American Chemical Society
Volume:
147
Issue:
45
ISSN:
0002-7863
Page Range / eLocation ID:
41809 to 41818
Subject(s) / Keyword(s):
Metal organic framework, rotor arrays, inertial rotation, molecular dynamics simulations
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Scott J. Miller (Ed.)
    Ground state destabilization is a promising strategy to modulate rotational barriers in amphidynamic crystals. DFT studies of polar phenylenes installed as rotators in pillared-paddle wheel metal-organic frameworks were performed to investigate the effects of ground state destabilization on their rotational dynamics. We found that as the steric size of phenylene substituents increases the ground state destabilization effect is also increased. Specifically, a significant destabilization of the ground state energy occurred as the size of the substituents increased, with values ranging from 2 kcal/mol to 11.7 kcal/mol. An evalua-tion of the effects of substituents on dipole-dipole interaction energies and rotational barriers suggest that it should be possi-ble to engineer amphidynamic crystals where the dipole-dipole interaction energy becomes comparable to the rotational barri-ers. Notably, dipole-dipole interaction energies reached values ranging from 0.6 kcal/mol to 2.4 kcal/mol. We propose that careful selection of polar substituents with different size may help create temperature-responsive materials with switchable collective polarization. 
    more » « less
  2. Using a combination of two-dimensional infrared (2D IR) and variable temperature Fourier transform infrared (FTIR) spectroscopies the rapid structural isomerization of a five-coordinate ruthenium complex is investigated. In methylene chloride, three exchanging isomers were observed: (1) square pyramidal equatorial, ( 1 ); (2) trigonal bipyramidal, ( 0 ); and (3) square pyramidal apical, ( 2 ). Exchange between 1 and 0 was found to be an endergonic process (Δ H = 0.84 (0.08) kcal mol −1 , Δ S = 0.6 (0.4) eu) with an isomerization time constant of 4.3 (1.5) picoseconds (ps, 10 −12 s). Exchange between 0 and 2 however was found to be exergonic (Δ H = −2.18 (0.06) kcal mol −1 , Δ S = −5.3 (0.3) eu) and rate limiting with an isomerization time constant of 6.3 (1.6) ps. The trigonal bipyramidal complex was found to be an intermediate, with an activation barrier of 2.2 (0.2) kcal mol −1 and 2.4 (0.2) kcal mol −1 relative to the equatorial and apical square pyramidal isomers respectively. This study provides direct validation of the mechanism of Berry pseudorotation – the pairwise exchange of ligands in a five-coordinate complex – a process that was first described over fifty years ago. This study also clearly demonstrates that the rate of pseudorotation approaches the frequency of molecular vibrations. 
    more » « less
  3. With dual goals of efficient and accurate modeling of solvation thermodynamics in molten salt liquids, we employ ab initio molecular dynamics (AIMD) simulations, deep neural network interatomic potentials (NNIP), and quasichemical theory (QCT) to calculate the excess chemical potentials for the solute ions Na + and Cl − in the molten NaCl liquid. NNIP-based molecular dynamics simulations accelerate the calculations by 3 orders of magnitude and reduce the uncertainty to 1 kcal mol −1 . Using the Density Functional Theory (DFT) level of theory, the predicted excess chemical potential for the solute ion pair is −178.5 ± 1.1 kcal mol −1 . A quantum correction of 13.7 ± 1.9 kcal mol −1 is estimated via higher-level quantum chemistry calculations, leading to a final predicted ion pair excess chemical potential of −164.8 ± 2.2 kcal mol −1 . The result is in good agreement with a value of −163.5 kcal mol −1 obtained from thermo-chemical tables. This study validates the application of QCT and NNIP simulations to the molten salt liquids, allowing for significant insights into the solvation thermodynamics crucial for numerous molten salt applications. 
    more » « less
  4. Abstract Reaction of {LiC6H2−2,4,6‐Cyp3⋅Et2O}2(Cyp=cyclopentyl) (1) of the new dispersion energy donor (DED) ligand, 2,4,6‐triscyclopentylphenyl with SnCl2afforded a mixture of the distannene {Sn(C6H2−2,4,6‐Cyp3)2}2(2), and the cyclotristannane {Sn(C6H2−2,4,6‐Cyp3)2}3(3).2is favored in solution at higher temperature (345 K or above) whereas3is preferred near 298 K. Van't Hoff analysis revealed the3to2conversion has a ΔH=33.36 kcal mol−1and ΔS=0.102 kcal mol−1 K−1, which gives a ΔG300 K=+2.86 kcal mol−1, showing that the conversion of3to2is an endergonic process. Computational studies show that DED stabilization in3is −28.5 kcal mol−1per {Sn(C6H2−2,4,6‐Cyp3)2unit, which exceeds the DED energy in2of −16.3 kcal mol−1per unit. The data clearly show that dispersion interactions are the main arbiter of the3to2equilibrium. Both2and3possess large dispersion stabilization energies which suppress monomer dissociation (supported by EDA results). 
    more » « less
  5. The reactions of thioformaldehyde (H 2 CS) with OH radicals and assisted by a single water molecule have been investigated using high level ab initio quantum chemistry calculations. The H 2 CS + ˙OH reaction can in principle proceed through: (1) abstraction, and (2) addition pathways. The barrier height for the addition reaction in the absence of a catalyst was found to be −0.8 kcal mol −1 , relative to the separated reactants, which has a ∼1.0 kcal mol −1 lower barrier than the abstraction channel. The H 2 CS + ˙OH reaction assisted by a single water molecule reduces the barrier heights significantly for both the addition and abstraction channels, to −5.5 and −6.7 kcal mol −1 respectively, compared to the un-catalyzed H 2 CS + ˙OH reaction. These values suggest that water lowers the barriers by ∼6.0 kcal mol −1 for both reaction paths. The rate constants for the H 2 CS⋯H 2 O + ˙OH and OH⋯H 2 O + H 2 CS bimolecular reaction channels were calculated using Canonical Variational Transition state theory (CVT) in conjunction with the Small Curvature Tunneling (SCT) method over the atmospherically relevant temperatures between 200 and 400 K. Rate constants for the H 2 CS + ˙OH reaction paths for comparison with the H 2 CS + ˙OH + H 2 O reaction in the same temperature range were also computed. The results suggest that the rate of the H 2 CS + ˙OH + H 2 O reaction is slower than that of the H 2 CS + ˙OH reaction by ∼1–4 orders of magnitude in the temperatures between 200 and 400 K. For example, at 300 K, the rates of the H 2 CS + ˙OH + H 2 O and H 2 CS + ˙OH reactions were found to be 2.2 × 10 −8 s −1 and 6.4 × 10 −6 s −1 , respectively, calculated using [OH] = 1.0 × 10 6 molecules cm −3 , and [H 2 O] = 8.2 × 10 17 molecules cm −3 (300 K, RH 100%) atmospheric conditions. Electronic structure calculations on the H 2 C(OH)S˙ product in the presence of 3 O 2 were also performed. The results show that H 2 CS is removed from the atmosphere primarily by reacting with ˙OH and O 2 to form thioformic acid, HO 2 , formaldehyde, and SO 2 as the main end products. 
    more » « less