Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Non-covalent bond interaction”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

34 records · Page 2
All sourcesNASA NTRSDOE OSTI
NASA NTRS2012

Is Water a Universal Solvent for Life?

There are strong reasons to believe that the laws, principles and constraints of physics and chemistry are universal. It is much less clear how this universality translates into our understanding of the origins of life. Conventionally, discussions of this topic focus on chemistry that must be sufficiently rich to seed life. Although this is clearly a prerequisite for the emergence of living systems, I propose to focus instead on self-organization of matter into functional structures capable of reproduction, evolution and responding to environmental changes. In biology, most essential functions are largely mediated by noncovalent interactions (interactions that do not involve making or breaking chemical bonds). Forming chemical bonds is only a small part of what living systems do. There are specific implications of this point of view for universality. I will concentrate on one of these implications. Strength of non-covalent interactions must be properly tuned. If they were too weak, the system would exhibit undesired, uncontrolled response to natural fluctuations of physical and chemical parameters. If they were too strong kinetics of biological processes would be slow and energetics costly. This balance, however, is not a natural property of complex chemical systems. Instead, it has to be achieved with the aid of an appropriate solvent for life. In particular, potential solvents for life must be characterized by a high dielectric constant to ensure solubility of polar species and sufficient flexibility of biological structures stabilized by electrostatic interactions. Among these solvents, water exhibits a remarkable trait that it also promotes solvophobic (hydrophobic) interactions between non-polar species, typically manifested by a tendency of these species to aggregate and minimize their contacts with the aqueous solvent. Hydrophobic interactions are responsible, at least in part, for many self-organization phenomena in biological systems, such as the formation of cellular boundary structures or protein folding. Strengths of electrostatic and hydrophobic interactions are similar and can be balanced over a wide range of temperatures, which considerably increases the repertoire of interactions that can be used to modulate biological functions. Some properties of water, e.g. its chemical activity against polymerization reactions, are considered as unfavorable to life. In actuality, this might be a favorable trait because life requires a balance between constructive and destructive processes. For example, molecules synthesized in response to specific conditions must be degraded once these conditions change. Otherwise regulation of biological processes would be virtually impossible. Water might not be the only liquid with favorable properties for supporting life. It has been proposed that formamide, which might be present elsewhere in the universe in sufficient quantities to warrant interest, could be a potential alternative to water for the origin of life. However, this will remain highly hypothetical until it is demonstrated in further studies on its physical, chemical and biological properties it is capable of mediating self-organization of matter and providing proper balance between different types of non-covalent interactions.

Pohorill, Andrew↗
DOE OSTI2021

Single crystals of mechanically entwined helical covalent polymers

Double helical conformation of polymer chains is widely observed in biomacromolecules and plays an essential role in exerting their biological functions, such as molecular recognition and information storage. It has remained challenging, however, to prepare synthetic helical polymers, and those that exist have mainly been limited to single-stranded polymers or short oligomeric double helices. In this work, we report the synthesis of covalent helical polymers, with a high molecular weight, from the achiral monomer hexahydroxytriphenylene through to spiroborate formation. Polymerization and crystallization occurred simultaneously under solvothermal conditions to form single crystals of the resulting helical covalent polymers. Characterization by single-crystal X-ray diffraction showed that each crystal consisted of pairs of mechanically entwined polymers. No strong non-covalent interactions were observed between the two helical polymers that formed a pair; instead, each strand interacted with neighbouring pairs through hydrogen bonding. Each individual crystal was made up of helical polymers of the same handedness, but the crystallization process produced a racemic conglomerate, with equal amounts of right-handed and left-handed crystals.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Construction of nano slow-release systems for antibacterial active substances and its applications: A comprehensive review

At present, nano-carrier materials with antibacterial activity are of great significance. Due to the widespread resistance of many pathogenic microorganisms, it has seriously threatened human health. The natural antimicrobial substances extracted from fruits and vegetables can significantly improve their stability combined with nano-carrier materials. The resistance of pathogenic microorganisms will be substantially reduced, greatly enhancing the effect of active antimicrobial substances. Nanotechnology has excellent research prospects in the food industry, antibacterial preservation, food additives, food packaging, and other fields. This paper introduces nano-carrier materials and preparation techniques for loading and encapsulating active antibacterial substances in detail by constructing a nano-release system for active antibacterial substances. The antibacterial effect can be achieved by protecting them from adverse external conditions and destroying the membrane of pathogenic microorganisms. The mechanism of the slow release of the bacteriostatic active substance is also described. The mechanism of carrier loading and release is mainly through non-covalent forces between the bacteriostatic active substance and the carrier material, such as hydrogen bonding, π-π stacking, van der Waals forces, electrostatic interactions, etc., as well as the loading and adsorption of the bacteriostatic active substance by the chemical assembly. Finally, its wide application in food and medicine is introduced. It is hoped to provide a theoretical basis and technical support for the efficient utilization and product development of bacteriostatic active substances.

Cao, Jiayong↗
DOE OSTI2023

ORNL_AISD-Ex: Quantum chemical prediction of UV/Vis absorption spectra for over 10 million organic molecules

We performed calculations of electronic excitation energies and associated oscillator strengths based on the time-dependent density-functional tight-binding (TD-DFTB) method [1]. The SMILES (Simplified molecular-input line-entry system) strings of the molecules from the AISD HOMO-LUMO database [2] were converted to a 3D atomistic structure and stored in a PDB file after preliminary geometry optimization using the Merck Molecular Force Field (MMFF94) in RDKit [3,4]. The primary information stored in the PDB file archive consists of Cartesian coordinates for each atom of the molecule in their 3D location in space, along with summary information about the structure, sequence, and experiment. We then performed molecular geometry optimization using the density-functional tight-binding (DFTB) method [5] in the electronic ground state, followed by single-point excited states calculations, as described below. We note that, since RDKit employs a random choice for the generation of molecular conformers, the molecular geometries obtained in this dataset could be different from the ones that were generated when the AISD HOMO-LUMO dataset was generated. The computed excitation energies and associated oscillator strengths can be converted to predict UV/Vis absorption spectra, where excitation energies correspond to absorption peak positions, and oscillator strengths are a good measure of the probability of absorption of visible or UV light in transitions between electronic ground and excited states. The conversion of SMILES strings to 3D Cartesian coordinates of fully DFTB-optimized molecules was successful for 10,502,904 out of 10,502,917 molecules. For these molecules, both geometry optimizations and excited states calculations were successful. The DFTB calculations did not complete for 13 molecules of the original AISD HOMO-LUMO dataset. We still provide information about the geometry of these molecules. The molecules are diverse for chemical compositions (which span 5 non-hydrogen elements: oxygen, carbon, nitrogen, fluorine, sulfur) and molecular size (the smallest molecule contains 5 non-hydrogen atoms, and the largest molecule contains 71 non-hydrogen atoms). The DFTB method [5] is an approximation to density functional theory (DFT), utilizing a minimal basis set in conjunction with a two-center approximation to the electronic Hamiltonian and overlap matrix elements. The DFTB total energy is the sum of an electronic and a repulsive energy contribution, and their calculation requires optimized electronic parameters and diatomic repulsive potential energy functions. All DFTB calculations were performed using the DFTB+ code [6] (version 21.2) and the wrapper for DFTB+ in the Atomic Simulation Environment (ASE) (version 3.22.1) [7], which performed an internal conversion of Cartesian coordinates from PDB to the .gen file format. For the geometry optimizations on the electronic ground state potential energy surface of the molecules, we have chosen the third-order DFTB (DFTB3) method [5c] and employed the matching 3ob set of electronic parameters and repulsive potentials [8]. The empirical γ-damping for hydrogen bond correction, and Grimme's D3 empirical dispersion correction with Becke-Johnson damping (D3(BJ)) [9] dispersion correction was included to improve the description of non-covalent interactions. For excited states single-point energy calculations, we employed the TD-DFTB method in conjunction with the DFTB2 method [5b] and the matching mio [5b,10] and halorg [11] parameter sets. We opted to request the simultaneous calculation of 50 excited states for singlet transition to investigate sufficient number of excited states, based on linear response theory using the Casida equation [Ref: T. A. Niehaus, S. Suhai, F. Della Sala, P Lugli, M. Elstner, G. Seifert, and Th. Frauenheim. Tight-binding approach to time-dependent density-functional response theory. Phys. Rev. B, 63:085108, 2001] and the ARPACK diagonalizer [R. B. Lehoucq, D. C. Sorensen, and C. Yang. Arpack users guide: Solution of large-scale eigenvalue problems by implicitly restarted arnoldi methods, 1997. 46, 51]. The dataset contains 1001 tar.gz files. Tar files are named as “ornl_aisd_ex_1.tar.gz†through “ornl_aisd_ex_1000.tar.gzâ€. Additionally, the 13 failed molecules are in “ornl_aisd_ex_unprocessed.tar.gzâ€. Except for the tar files listed below, each tar file contains 10,500 molecules. Tar files numbered 34, 121, 128, 352, 360, 429, 495, 509, 518, 627, 676, 668, and 862 contain 10,499 molecules each. The last tar file numbered 1000 contains 13,417 molecules. The total size of the uncompressed dataset is over 283 Gigabytes. The code for calculating the electronic excitation energies and statistical analysis of the dataset is provided at the following GitLab repository: https://github.com/ORNL/Analysis-of-Large-Scale-Molecular-Datasets-with-Python Calculating the UV spectrum of a molecule requires performing 3 main operations: 1. Converting the smiles string representation of a molecule into a geometric structure where each atom is assigned XYZ coordinates. The geometric structure is written to the file smiles.pdb. 2. Using smiles.pdb to compute the relaxed geometry of the molecule, which corresponds with the position of the atoms at the position of equilibrium at the ground state. This generates the files band.out, detailed.out, and geo_end.gen. 3. Using geo_end.gen to calculate the UV spectrum of the molecule which is written into the file EXC.DAT. Every molecule in the dataset has its own directory. The files contained in each molecule directory are as follows: 1. geo_end.gen 2. detailed.out 3. band.out 4. EXC.DAT 5. smiles.pdb REFERENCES [1] Niehaus, T. A.; Suhai, S.; Della Salla, F.; Lugli, P.; Elstner, M.; Seifert, G.; Frauenheim, Th. Tight-binding approach to time-dependent density-functional response theory. Phys. Rev. B, 2001, 63, 085108/1-9. [2] Blanchard, A.; Gounley, J.; Metha, K.; Yoo, P.; Irle, S. AISD HOMO-LUMO. DOI: 10.13139/ORNLNCCS/1869409 [3] RDKit: Cheminformatics and Machine Learning Software. 2013, [http://www.rdkit.org] [4] Tosco, P.; Stiefl, N. and Landrum, G. Bringing the MMFF force field to the RDKit: implementation and validation. J Cheminform. 2014, 6, 1–4. [5] a) Porezag, D.; Frauenheim, T.; Kohler, T.; Seifert, G.; Kaschner, Construction of tight-binding-like potentials on the basis of density-functional theory: Application to carbon, R. Phys. Rev. B 1995, 51, 12947-12957; b) Elstner, M.; Porezag, D.; Jungnickel, G.; Elsner, J.; Haugk, M.; Frauenheim, Th.; Suhai, S.; Seifert, G.; Phys. Rev. B 1998, 58, 7260-7268; c) Gaus, M.; Cui, Q.; Elstner, M. DFTB3: Extension of the Self-Consistent-Charge Density-Functional Tight-Binding Method (SCC-DFTB), J. Chem. Theory Comput. 2011, 7, 931-948; d) Cui, Q.; Elstner, M. Density functional tight binding: values of semi-empirical methods in an ab initio era, Phys. Chem. Chem. Phys. 2014, 16, 14368-14377. [6] Hourahine, B. et al. DFTB+, a software package for efficient approximate density functional theory based atomistic simulations, J. Chem. Phys. 2020, 152, 124101/1-19. [7] Larsen, A. H. et al. The atomic simulation environment—a Python library for working with atoms. J. Phys.: Cond. Matter 2017, 29, 273002. [8] Kubillus, M.; Kubar, T.; Gaus, M.; Rezac, J.; Elstner, M. Parameterization of the DFTB3 Method for Br, Ca, Cl, F, I, K, and Na in Organic and Biological Systems, J. Chem. Theory Comput. 2015, 11, 332-342. [9] Brandenburg, J. G.; Grimme, S. Accurate Modeling of Organic Molecular Crystals by Dispersion-Corrected Density Functional Tight Binding (DFTB), J. Phys. Chem. Lett. 2014, 5, 1785−1789. [10] a) Niehaus, T. A.; Elstner, M.; Frauenheim, Th.; Suhai, S. Application of an approximate density-functional method to sulfur containing compounds. J. Mol. Struct.: THEOCHEM 2001, 541, 185-94; b) Elstner, M.; Hobza, P.; Frauenheim, Th.; Suhai, S.; Kaxiras, E. Hydrogen bonding and stacking interactions of nucleic acid base pairs: A density-functional-theory based treatment. J. Chem. Phys. 2001, 114, 5149-55. [11] Kubar, T.; Bodrog, Z.; Gaus, M.; Köhler, C.; Aradi, B.; Frauenheim, Th.; Elstner, M. Parametrization of the SCC-DFTB Method for Halogens. J. Chem. Theory Comput. 2013, 9, 2939-49.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Anion Binding by Macrocyclic Receptors: Computational Landscape of 1:1 and 2:1 Stoichiometries

Macrocyclic receptors play a crucial role in supramolecular chemistry, enabling the selective binding of guest molecules through non-covalent interactions. This study investigates the anion binding properties of three shape-persistent macrocycles, triazolophane, cyanostar, and tricarb, each featuring preorganized cavities with polarized CH hydrogen bond donors. In the 1:1 binding stoichiometries, the macrocycles preferentially bind anions that fit inside their two-dimensional cavities, with this preference being more pronounced in solution. Specifically, the triazolophane favors smaller anions like chloride and bromide, while the larger cavities offered by cyanostar and tricarb macrocycles preferentially bind the larger iodide. For large, polyatomic anions (SCN − , BF$^{−}_{4}$, ClO$^{−}_{4}$, and PF$^{−}_{6}$), the macrocycles self-assemble to form stable 2:1 sandwich complexes, exhibiting stronger binding. This shift in selectivity highlights the versatility of the macrocycles, making them promising candidates for anion sensing and regulation with various complexation stoichiometries.

Che, Minwei [Indiana University, Bloomington, IN (↗
DOE OSTI2022

Supramolecular assemblies of polybenzimidazole and aromatic polycarboxylic acids with superior mechanical and H 2 /CO 2 separation properties

Aromatic polycarboxylic acids complex with imidazoles, forming supramolecular assemblies. Herein, we demonstrate that polybenzimidazole (PBI) can be augmented by aromatic polycarboxylic acids, including phthalic acid (PA), trimesic acid (TMA), and pyromellitic acid (PMA), and the resulting supramolecular assemblies exhibit attractive H 2 /CO 2 separation performance for precombustion carbon capture. The acid doping decreases the free volume and gas permeability but increases H 2 /CO 2 selectivity, which can be correlated with the molar ratio of the protons to PBI repeating units in the assemblies. Increasing the temperature decreases the H 2 /CO 2 selectivity because of weakened interactions. When challenged with model gas mixtures, a supramolecular assembly based on TMA shows H 2 /CO 2 separation properties superior to state-of-the-art materials and above Robeson's upper bound. Furthermore, this study elucidates the fruitful harnessing of dynamic non-covalent bonds to design robust and highly selective membranes for industrial gas separations.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Improved accuracy and transferability of molecular-orbital-based machine learning: Organics, transition-metal complexes, non-covalent interactions, and transition states

We report molecular-orbital-based machine learning (MOB-ML) provides a general framework for the prediction of accurate correlation energies at the cost of obtaining molecular orbitals. The application of Nesbet’s theorem makes it possible to recast a typical extrapolation task, training on correlation energies for small molecules and predicting correlation energies for large molecules, into an interpolation task based on the properties of orbital pairs. We demonstrate the importance of preserving physical constraints, including invariance conditions and size consistency, when generating the input for the machine learning model. Numerical improvements are demonstrated for different datasets covering total and relative energies for thermally accessible organic and transition-metal containing molecules, non-covalent interactions, and transition-state energies. MOB-ML requires training data from only 1% of the QM7b-T dataset (i.e., only 70 organic molecules with seven and fewer heavy atoms) to predict the total energy of the remaining 99% of this dataset with sub-kcal/mol accuracy. This MOB-ML model is significantly more accurate than other methods when transferred to a dataset comprising of 13 heavy atom molecules, exhibiting no loss of accuracy on a size intensive (i.e., per-electron) basis. It is shown that MOB-ML also works well for extrapolating to transition-state structures, predicting the barrier region for malonaldehyde intramolecular proton-transfer to within 0.35 kcal/mol when only trained on reactant/product-like structures. Finally, the use of the Gaussian process variance enables an active learning strategy for extending the MOB-ML model to new regions of chemical space with minimal effort. We demonstrate this active learning strategy by extending a QM7b-T model to describe non-covalent interactions in the protein backbone–backbone interaction dataset to an accuracy of 0.28 kcal/mol.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Nickel(ii) complexes based on dithiolate–polyamine binary ligand systems: crystal structures, hirshfeld surface analysis, theoretical study, and catalytic activity study on photocatalytic hydrogen generation

To ascertain the influence of binary ligand systems [1,1-dicyanoethylene-2,2-dithiolate (i-mnt -2 ) and polyamine {tetraen = tris(2-aminoethyl)amine, tren = diethylene triamine and opda = o-phenylenediamine}] on the coordination modes of the Ni(II) metal center and resulting supramolecular architectures, a series of nickel(II) thiolate complexes [Ni(tetraen)(i-mnt)](DMSO) (1), [Ni 2 (tren) 2 (i-mnt) 2 ] (2), and [Ni 2 (i-mnt) 2 (opda) 2 ] n (3) have been synthesized in high yield in one step in water and structurally characterized by single crystal X-ray crystallography and spectroscopic techniques. X-ray diffraction studies disclose the diverse i-mnt -2 coordination to the Ni +2 center in the presence of active polyamine ligands, forming a slightly distorted octahedral geometry (NiN 4 S 2 ) in 1, square planar (NiS 4 ) and distorted octahedral geometries (NiN 6 ) in the bimetallic co-crystallized aggregate of cationic [Ni(tren) 2 ] +2 and anionic [Ni(i-mnt) 2 ] -2 in 2, and a one dimensional (1D) polymeric chain along the [100] axis in 3, having consecutive square planar (NiS 4 ) and octahedral (NiN 6 ) coordination kernels. The N–H···O, N–H···S, N–H···N, N–H···S, N–H···N, and N–H···O type hydrogen bonds stabilize the supramolecular assemblies in 1, 2, and 3 respectively imparting interesting graph-set-motifs. The molecular Hirshfeld surface analyses (HS) and 2D fingerprint plots were utilized for decoding all types of non-covalent contacts in the crystal networks. Atomic HS analysis of the Ni +2 centers reveals significant Ni–N metal–ligand interactions compared to Ni–S interactions. We have also studied the unorthodox interactions observed in the solid state structures of 1–3 by QTAIM and NBO analyses. Moreover, all the complexes proved to be highly active water reduction co-catalysts (WRC) in a photo-catalytic hydrogen evolution process involving iridium photosensitizers, wherein 2 and 3 having a square planar arrangement around the nickel center(s) – were found to be the most active ones, achieving 1000 and 1119 turnover numbers (TON), respectively.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Multi-scale computational screening and mechanistic insights of cyclic amines as solvents for improved lignocellulosic biomass processing

A computational screening workflow for the efficient deconstruction of cellulose, lignin and hemicellulose fractions of lignocellulosic biomass using cyclic amines as solvents. Lignocellulosic biomass is a promising feedstock for production of affordable fuels and chemicals from renewable resources. Effective solubilization and subsequent deconstruction of its cellulose, hemicellulose, and lignin fractions is essential for the viability of future biorefineries. This study used quantum chemistry-based equilibrium thermodynamics methods to evaluate the potential of 650 cyclic amines to solubilize cellulose, hemicellulose, and lignin. The activity coefficients of solvent - biopolymer interactions were predicted using the COSMO-RS (COnductor-like Screening MOdel for Real Solvents) method and used to identify cyclic amines that can efficiently dissolve and extract selective fractions of biopolymers during biomass pretreatment. Among the 650 cyclic amines, 1-piperazineethanmaine was predicted to be an effective solvent for extracting all three polymers and was experimentally shown to achieve the highest lignin removal (97.1%). Non-covalent interaction, reduced density gradient and quantum chemical calculations were performed to elucidate the dissolution mechanism of lignin, cellulose and hemicellulose and gain further molecular level insights into the interactions between the cyclic amines and biomass polymers that promote efficient solubilization and extraction. These analyses indicated that 1-piperazineethanmaine and 1-methylimidazole make noncovalent van der Waals, electrostatic interactions and hydrogen bonding with lignin, leading to enhanced lignin removal, while the strong intramolecular hydrogen bonding interactions in cellulose and hemicellulose result in weaker solvent-biopolymer interactions. Overall, the computational approach provided an efficient method for identifying cyclic amines tailored for optimal biomass pretreatment and resulted in the identification of a potential new class of solvents for effective biomass pretreatment.

Kumar, Nikhil↗
DOE OSTI2022

On the mechanism of intermolecular nitrogen-atom transfer from a lattice-isolated diruthenium nitride intermediate

Catalyst confinement within microporous media provides the opportunity to site isolate reactive intermediates, enforce intermolecular functionalization chemistry by co-localizing reactive intermediates and substrates in molecular-scale interstices, and harness non-covalent host–guest interactions to achieve selectivities that are complementary to those accessible in solution. As part of an ongoing program to develop synthetically useful nitrogen-atom transfer (NAT) catalysts, we have demonstrated intermolecular benzylic amination of toluene at a Ru 2 nitride intermediate confined within the interstices of a Ru2-based metal–organic framework (MOF), Ru 3 (btc) 2 X 3 (btc = 1,3,5-benzenetricarboxylate, i.e., Ru-HKUST-1 for X = Cl). Nitride confinement within the extended MOF lattice enabled intermolecular C–H functionalization of benzylic C–H bonds in preference to nitride dimerization, which was encountered with soluble molecular analogues. Detailed study of the kinetic isotope effects (KIEs, i.e., k H /k D ) of C–H amination, assayed both as intramolecular effects using partially labeled toluene and as intermolecular effects using a mixture of per-labeled and unlabeled toluene, provided evidence for restricted substrate mobility on the time scale of interstitial NAT. Analysis of these KIEs as a function of material mesoporosity provided approximate experimental values for functionalization in the absence of mass transport barriers. Here, we disclose a combined experimental and computational investigation of the mechanism of NAT from a Ru 2 nitride to the C–H bond of toluene. Computed kinetic isotope effects for a H-atom abstraction (HAA)/radical rebound (RR) mechanism are in good agreement with experimental data obtained for C–H amination at the rapid diffusion limit. These results provide the first detailed analysis of the mechanism of intermolecular NAT to a C–H bond, bolster the use of KIEs as a probe of confinement effects on NAT within MOF lattices, and provide mechanistic insights unavailable by experiment because rate-determining mass transport obscured the underlying chemical kinetics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

The approximate second order coupled-cluster method based on a size-consistent Brillouin–Wigner partitioning

We present a variant of the approximate second order coupled-cluster method (CC2) with a two-parameter size-consistent Brillouin–Wigner (BW-s) partitioning instead of a Møller–Plesset (MP) partitioning for the unperturbed Hamiltonian, which we refer to as BWs-CC2. The computational complexity of this model scales identically to CC2 with molecular size. Conventional CC2 and its regularized BWs-CC2 variants, as well as conventional MP2 and two of its regularized BW-s2 variants, were assessed on a 535 element database spanning thermochemistry, non-covalent interactions, barrier heights, and isomerization energies. To ensure a well-defined model chemistry, the assessment was performed using internally stable spin-polarized Hartree–Fock (HF) orbitals in the finite aug-cc-pVQZ basis without counterpoise corrections. As a result of using stable orbitals, contrary to conventional wisdom, we find that CC2 substantially outperforms MP2 on molecules with significantly spin contaminated reference orbitals without a significant increase in error on systems with a spin-pure reference, showing the value of its single substitutions. While no single choice of regularization parameters can be optimal for all datasets, we find that BWs-CC2 generally outperforms both CC2 and BW-s2 with a single judicious parameter choice. Additional tests on dipole moments and bond lengths of diatomics provide further support for the utility of this choice. Furthermore, the main outliers and poorest performing cases are associated with large amounts of spin-contamination in the HF reference, which is indicative of systems with either strong correlation or extensive artificial symmetry breaking. Overall, these findings argue that the perception of the quality of the CC2 ground state should be reevaluated and that it can be further improved upon by the soundly based BWs-CC2 variant with the recommended parameter choice.

Correlation energy↗
DOE OSTI2021

DFT Mechanism Studies: Biomimetic 1,4-NADH Chemoselective, Co-factor Regeneration with [Cp*Rh(bpy)H] + , in Tandem with the Biocatalysis Pathways of a Core Model of the (HLADH)-Zn(II) Mediated Enzyme, in the Enantioselective Reduction of Achiral Ketones to Chiral S-Alcohols

In this study, Quantum Chemical (QC) calculations, utilizing Density Functional Theory (DFT), were performed to investigate the mechanistic aspects of the chemoselective catalyzed reaction of [Cp*Rh(bpy)H] + with the biomimetic NAD + analogues, N-benzylnicotinamide triflate, 1, and β-nicotinamide ribose-5'-methyl phosphate, 2, in the conversion to their 1,4-NADH analogues, 1,4-dihydro-N-benzylnicotinamide, 4, and β-1,4-dihydronicotinamide-5'-ribose methyl phosphate, 5. This reaction was in tandem with the 1,4-NADH dependent HLADH-Zn(II)- catalyzed reduction of achiral ketones to chiral S-alcohols. The [Cp*Rh(bpy)H] + complex, and not its equilibrium tautomer, [η 4 -Cp*HRh(bpy)] + , was found to control the hydride transfer during the biomimetic NAD + /1,4-NADH conversion, through the non-covalent interactions of the biomimetic co-factors with [Cp*Rh(bpy)H] + . The thermodynamics and kinetics for the chiral reduction of the Zn(II) bound ketones, 2-pentanone and 4-phenyl-2-butanone, with co-factor, 4, catalyzed by Zn(SCH 3 ) 2 (Imidazole), a core model of the Zn(II)-based catalytic center of HLADH, was also investigated by the evaluation of two possible reaction pathways: (1) formation of a ZnH from the C4-H hydride transfer of co-factor, 4, followed by reaction of the postulated ZnH with the bound 2-pentanone or 4-phenyl-2-butanone substrate, and (2), the direct C4-H transfer to the bound achiral ketone substrates, to provide the dominant chiral alcohols, S-2-pentanol or S-4-phenyl-2-butanol. The latter pathway was found most viable, and DFT calculations also revealed an essential η 2 -coordination of the 5,6 double bond of co-factor, 4, to the HLADH-Zn(II) metal ion center, upon imidazole decomplexation, providing an asymmetric differentiation of S-η 2 -5,6-1,4-NADH-Zn(II) binding. A proposed new paradigm for the Zn(II)'s non-innocent role in the HLADH-Zn(II) biocatalysis reduction mechanism, for enantioselective hydride transfer to a Zn(II) bound ketone, providing S-alcohols.

1,4 NADH co-factors↗
DOE OSTI2021

NENCI-2021. I. A large benchmark database of non-equilibrium non-covalent interactions emphasizing close intermolecular contacts

In this work, we present NENCI-2021, a benchmark database of ~8000 Non-Equilibirum Non-Covalent Interaction energies for a large and diverse selection of intermolecular complexes of biological and chemical relevance. To meet the growing demand for large and high-quality quantum mechanical data in the chemical sciences, NENCI-2021 starts with the 101 molecular dimers in the widely used S66 and S101 databases and extends the scope of these works by (i) including 40 cation–π and anion–π complexes, a fundamentally important class of non-covalent interactions that are found throughout nature and pose a substantial challenge to theory, and (ii) systematically sampling all 141 intermolecular potential energy surfaces (PESs) by simultaneously varying the intermolecular distance and intermolecular angle in each dimer. Designed with an emphasis on close contacts, the complexes in NENCI-2021 were generated by sampling seven intermolecular distances along each PES (ranging from 0.7× to 1.1× the equilibrium separation) and nine intermolecular angles per distance (five for each ion–π complex), yielding an extensive database of 7763 benchmark intermolecular interaction energies (E int ) obtained at the coupled-cluster with singles, doubles, and perturbative triples/complete basis set [CCSD(T)/CBS] level of theory. The E int values in NENCI-2021 span a total of 225.3 kcal/mol, ranging from -38.5 to +186.8 kcal/mol, with a mean (median) E int value of -1.06 kcal/mol (-2.39 kcal/mol). In addition, a wide range of intermolecular atom-pair distances are also present in NENCI-2021, where close intermolecular contacts involving atoms that are located within the so-called van der Waals envelope are prevalent—these interactions, in particular, pose an enormous challenge for molecular modeling and are observed in many important chemical and biological systems. A detailed symmetry-adapted perturbation theory (SAPT)- based energy decomposition analysis also confirms the diverse and comprehensive nature of the intermolecular binding motifs present in NENCI-2021, which now includes a significant number of primarily induction-bound dimers (e.g., cation–π complexes). NENCI-2021 thus spans all regions of the SAPT ternary diagram, thereby warranting a new four-category classification scheme that includes complexes primarily bound by electrostatics (3499), induction (700), dispersion (1372), or mixtures thereof (2192). A critical error analysis performed on a representative set of intermolecular complexes in NENCI-2021 demonstrates that the E int values provided herein have an average error of ±0.1 kcal/mol, even for complexes with strongly repulsive E int values, and maximum errors of ±0.2–0.3 kcal/mol (i.e., ~±1.0 kJ/mol) for the most challenging cases. For these reasons, we expect that NENCI-2021 will play an important role in the testing, training, and development of next-generation classical and polarizable force fields, density functional theory approximations, wavefunction theory methods, and machine learning based intra- and inter-molecular potentials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Uncovering the mechanism of selective stabilization of high-energy diastereoisomers via inclusion

Abstract Supramolecular systems may be used to stabilize otherwise unstable isomers to find alternative synthetic pathways. It has been reported that cucurbit[8]uril can stabilize trans -I and trans -II Cu $$^{{\textrm{II}}}$$ II cyclam, whereas trans -III is the only non-substituted trans Cu $$^{{\textrm{II}}}$$ II cyclam diastereoisomer found outside of the host molecule experimentally. Quantum chemistry methods can provide valuable insight into the intermolecular interactions involved in these inclusion complexes. All five possible trans diastereoisomers of Cu $$^{{\textrm{II}}}$$ II cyclam were studied within the host molecule to calculate the interaction energy and free energy of association for each complex. The relative free energies of the five free cyclams confirm that trans -I and trans -II are the most energetically accessible diastereoisomers from the initial trans -III starting point. Energy decomposition analysis was used to identify the attractive and repulsive interactions between cyclam and cucurbit[8]uril and showed that trans -II encounters repulsive forces almost three times greater than trans -I, which may explain the 7:3 ratio of trans -I to trans -II within cucurbit[8]uril that occurs experimentally. Optimized complex geometries with trans -III, IV, and V show that the cyclams protrude out of cucurbit[8]uril, whereas trans -I and trans -II become more encapsulated and elongate the host, suggesting that the position of the cyclam is extremely important when forming non-covalent interactions. Our results agree with the experimental findings and provide greater insight into why the most stable isolated cyclam diastereoisomer, trans -III, does not form a complex.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2026

A Permanently Porous Chalcogen-Bonded Organic Framework

The nature of connectivity between constituent atomic or molecular building blocks is fundamental in shaping the properties and functionality of materials. The extrapolation of emergent interatomic interactions to enable functional materials has driven transformative technological advancements. However, the bonding interactions used in material design have been largely static since the emergence of dynamic covalent chemistry ~30 years ago. Here we demonstrate that non-covalent chalcogen bonding (Ch-bonding) is a distinct mode of interatomic connectivity for constructing functional materials by design. This is established by leveraging self-complementary assembly of 1,2,5-telluradiazole moieties to construct a honeycomb-type permanently porous Ch-bonded organic framework, assembled and stabilized solely through non-covalent Te...N contacts. Empirical and computational studies of electronic structure, structural healing and lattice dynamics highlight the pi-type electronic communication, controlled assembly and modulated lattice dynamics in Trip3Tez-I arising directly from the unique nature of the Te...N Ch-bonding that holds substantial implications for next generation crystalline semiconductors. In addition to introducing a distinct class of permanently porous frameworks, this work establishes Ch-bonding as a programmable molecular tool for constructing functional materials with distinct properties.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2011

The Importance of Water for High Fidelity Information Processing and for Life

Is water an absolute prerequisite for life? Life depends on a variety of non-covalent interactions among molecules, the nature of which is determined as much by the solvent in which they occur as by the molecules themselves. Catalysis and information processing, two essential functions of life, require non-covalent molecular recognition with very high specificity. For example, to correctly reproduce a string consisting of 600,000 units of information (e.g ., 600 kilobases, equivalent to the genome of the smallest free living terrestrial organisms) with a 90% success rate requires specificity > 107 : 1 for the target molecule vs. incorrect alternatives. Such specificity requires (i) that the correct molecular association is energetically stabilized by at least 40 kJ/mol relative to alternatives, and (ii) that the system is able to sample among possible states (alternative molecular associations) rapidly enough to allow the system to fall under thermodynamic control and express the energetic stabilization. We argue that electrostatic interactions are required to confer the necessary energetic stabilization vs. a large library of molecular alternatives, and that a solvent with polarity and dielectric properties comparable to water is required for the system to sample among possible states and express thermodynamic control. Electrostatic associations can be made in non-polar solvents, but the resulting complexes are too stable to be "unmade" with sufficient frequency to confer thermodynamic control on the system. An electrostatic molecular complex representing 3 units of information (e.g., 3 base pairs) with specificity > 107 per unit has a stability in non-polar solvent comparable to that of a carbon-carbon bond at room temperature. These considerations suggest that water, or a solvent with properties very like water, is necessary to support high-fidelity information processing, and can therefore be considered a critical prerequisite for life.

Hoehler, Tori M.↗