Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “molecular functionalization”

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.

At least 1,099 records · Page 61
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Importance of imposing gauge invariance in time-dependent density functional theory calculations with meta-generalized gradient approximations

It has been known for more than a decade that the gauge variance of the kinetic energy density τ leads to additional terms in the magnetic orbital rotation Hessian used in linear-response time-dependent density functional theory (TDDFT), affecting excitation energies obtained with τ-dependent exchange–correlation functionals. While previous investigations found that a correction scheme based on the paramagnetic current density has a small effect on benchmark results, we report more pronounced effects here, in particular, for the popular M06-2X functional and for some other meta-generalized gradient approximations (mGGAs). In the first part of this communication, this is shown by a reassessment of a set of five Ni(II) complexes for which a previous benchmark study that did not impose gauge invariance has found surprisingly large errors for excitation energies obtained with M06-2X. These errors are more than halved by restoring gauge invariance. The variable importance of imposing gauge invariance for different mGGA-based functionals can be rationalized by the derivative of the mGGA exchange energy integrand with respect to τ. In the second part, a large set of valence excitations in small main-group molecules is analyzed. For M06-2X, several selected n → π* and π→π$^{*}_{⊥}$ excitations are heavily gauge-dependent with average changes of –0.17 and –0.28 eV, respectively, while π→π$^{*}_{∥}$ excitations are marginally affected (–0.04 eV). Similar patterns, but of the opposite signs, are found for SCAN0. Here, the results suggest that reevaluation of previous gauge variant TDDFT results based on M06-2X and other mGGA functionals is warranted.

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

A multi-center quadrature scheme for the molecular continuum

We report a common way to evaluate electronic integrals for polyatomic molecules is to use Becke’s partitioning scheme (Becke and Chem, 1988) in conjunction with overlapping grids centered at each atomic site. The Becke scheme was designed for integrands that fall off rapidly at large distances, such as those approximating bound electronic states. When applied to states in the electronic continuum, however, Becke scheme exhibits slow convergence and it is highly redundant. Here, we present a modified version of Becke scheme that is applicable to functions of the electronic continuum, such as those involved in molecular photoionization and electron–molecule scattering, and which ensures convergence and efficiency comparable to those realized in the calculation of bound states. In this modified scheme, the atomic weights already present in Becke’s partition are smoothly switched off within a range of few bond lengths from their respective nuclei, and complemented by an asymptotically unitary weight. The atomic integrals are evaluated on small spherical grids, centered on each atom, with size commensurate to the support of the corresponding atomic weight. The residual integral of the interstitial and long-range region is evaluated with a central master grid. The accuracy of the method is demonstrated by evaluating integrals involving integrands containing Gaussian Type Orbitals and Yukawa potentials, on the atomic sites, as well as spherical Bessel functions centered on the master grid. These functions are representative of those encountered in realistic electron-scattering and photoionization calculations in polyatomic molecules.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Catalytic N 2 -to-NH 3 (or -N 2 H 4 ) Conversion by Well-Defined Molecular Coordination Complexes

Nitrogen fixation, the six-electron/six-proton reduction of N 2 , to give NH 3 , is one of the most challenging and important chemical transformations. Notwithstanding the barriers associated with this reaction, significant progress has been made in developing molecular complexes that reduce N 2 into its bioavailable form, NH 3 . This progress is driven by the dual aims of better understanding biological nitrogenases and improving upon industrial nitrogen fixation. In this review, we highlight both mechanistic understanding of nitrogen fixation that has been developed, as well as advances in yields, efficiencies, and rates that make molecular alternatives to nitrogen fixation increasingly appealing. We begin with a historical discussion of N 2 functionalization chemistry that traverses a timeline of events leading up to the discovery of the first bona fide molecular catalyst system and follow with a comprehensive overview of d-block compounds that have been targeted as catalysts up to and including 2019. Furthermore, we end with a summary of lessons learned from this significant research effort and last offer a discussion of key remaining challenges in the field.

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

Machine-learning accelerated geometry optimization in molecular simulation

Geometry optimization is an important part of both computational materials and surface science because it is the path to finding ground state atomic structures and reaction pathways. These properties are used in the estimation of thermodynamic and kinetic properties of molecular and crystal structures. This process is slow at the quantum level of theory because it involves an iterative calculation of forces using quantum chemical codes such as density functional theory (DFT), which are computationally expensive and which limit the speed of the optimization algorithms. It would be highly advantageous to accelerate this process because then one could do either the same amount of work in less time or more work in the same time. Here, we provide a neural network (NN) ensemble based active learning method to accelerate the local geometry optimization for multiple configurations simultaneously. We illustrate the acceleration on several case studies including bare metal surfaces, surfaces with adsorbates, and nudged elastic band for two reactions. In all cases, the accelerated method requires fewer DFT calculations than the standard method. In addition, we provide an Atomic Simulation Environment (ASE)-optimizer Python package to make the usage of the NN ensemble active learning for geometry optimization easier.

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

Targeted seed EMS mutagenesis reveals a basic helix–loop–helix transcription factor underlying male sterility in sorghum

Abstract Forward genetic screens of mutant populations are fundamental for functional genomics studies. However, isolating independent mutant alleles to molecularly identify causal genes is challenging in species recalcitrant to genetic manipulation. Here, we demonstrate that classic seed ethyl methanesulfonate (EMS) mutagenesis coupled with genome sequencing can overcome this limitation in sorghum. We used this method to generate new mutant alleles of sorghum MALE STERILE 8 (MS8) and identified the causal locus for the ms8 phenotype as Sobic.004G270900, which encodes the sorghum ortholog of maize bhlh122, a basic helix–loop–helix (bHLH) transcription factor required for male fertility in maize. Bulked segregant analysis mapped ms8-1 to a region on chromosome 4 containing Sobic.004G270900. Seeds from heterozygous MS8/ms8-1 plants were mutagenized and screened for chimeric inflorescences containing sectors with white, sterile anthers resembling the ms8-1 homozygous phenotype. DNA sequencing of sterile and fertile sectors from a single chimeric inflorescence revealed two mutations in Sobic.004G270900 within the sterile sector, but not the fertile sector. Isolation of this loss-of-function allele (ms8-2) established Sobic.004G270900 as the causative locus for male sterility in the ms8 mutant. We generated additional alleles of MS8 in a different genetic background using CRISPR/Cas9-based gene editing, where deletions in Sobic.004G270900 also resulted in male sterility. Our work identified a gene underlying male sterility in sorghum and provides a novel and straightforward genetic tool for researchers who lack access to advanced transformation facilities to validate gene candidates. Unlike gene editing, no prior knowledge of candidate genes is required for targeted seed EMS mutagenesis to aid identification of causal loci.

Genetics & Heredity↗
DOE OSTI2025

Selective Adsorption of Thiol-Containing Molecules on Copper Sulfide Surfaces via Molecule–Surface Disulfide Bridges

Recent results in the fields of nanoenhanced agriculture and expanding interest in prebiotic chemistry have placed increased emphasis on understanding the chemically selective interaction of small molecules with the surfaces of metal sulfides. Here, we present an integrated experimental and computational study of the interaction of thiol-containing molecules with copper sulfide (covellite) surfaces in aqueous media. In situ Fourier-transform infrared (FTIR) measurements and ex situ X-ray photoelectron spectroscopy (XPS) measurements show that molecules bearing free thiol groups, including glutathione and cysteine, bind strongly to CuS (covellite) nanoparticles and to CuS (001) single crystals, while control studies show that similar molecules lacking the free thiol group exhibit much less binding. Additional experiments show that these thiol-containing molecules interact transiently with CuO nanoparticle surfaces but are readily removed by rinsing. The FTIR and XPS experiments demonstrate that adsorption of molecular thiols to CuS surfaces occurs in a chemically selective manner. Further experimental studies and density functional calculations show that the preferred mode of binding is through the surface S atoms, forming a Solid–S–S–Molecule disulfide linkage. While the role of disulfide linkages in controlling structure and function of proteins and other biomolecules is widely known, the formation of surface disulfide linkages as a motif for covalent molecular binding at surfaces has not been established previously.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Spatiotemporal mapping of mesoscopic liquid dynamics

The study of liquid dynamics at mesoscopic scales is still strewn with difficulty due to limitations in theory and experiment. Historically, significant attention has been given to the analysis of space-time correlation functions and their frequency-Fourier transforms at a few discrete wave numbers. The massive computing power afforded by modern high performance computing clusters and the advent of a wide-angle neutron spin-echo spectrometer, however, have unlocked a more intuitive and fruitful approach to this problem. Using molecular dynamics simulations, here we demonstrate the benefits of spatiotemporally mapping intermediate scattering functions on a dense grid of correlation times and wave numbers. Four model systems are investigated: a Lennard-Jones liquid, a coarse-grained bead-spring polymer, a molten sodium chloride, and a poly(ethylene oxide) melt. We show that the spatiotemporal mapping approach is particularly useful for elucidating the mesoscopic dynamics in these liquids, where several underlying mechanisms, such as molecular relaxations, hydrodynamic modes, and nonhydrodynamic excitations, are potentially at play. Therefore, compared to the traditional method, direct visualization of density space-time correlation functions on two-dimensional color maps permits appraisals of complicated dynamical behavior at mesoscales in a global manner. For example, the scaling relations between space and time for different types of molecular motions can be straightforwardly identified on these plots, without any model-dependent analysis. Additionally, we show how theoretical ideas regarding collective mesoscopic dynamics, such as the classical hydrodynamic theory, the convolution approximation, and a recently proposed phenomenological model, can be discussed in terms of the global features of spatiotemporal maps of intermediate scattering functions. The new perspective offered by the spatiotemporal mapping method should prove useful for the study of liquid dynamics in general.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Targeting Aspergillus allergen oryzin with a chemical probe at atomic precision

We report the molecular basis of Aspergillus fumigatus oryzin, allergen Asp f 13, or alkaline proteinase ALP1, containing the sequence motif His–Asp–Ser of the subtilisin family, structure, and function at atomic detail. Given the resolution of the data (1.06 Å), we use fragment molecular replacement with ideal polyalanine α-helices to determine the first crystal structure of oryzin. We probe the catalytic serine through formation of an irreversible bond to a small molecule compound, specifically labeling it, describing the amino acid residues performing the catalytic function. Defined by a self-processed pro-peptide, the active site architecture shapes up pocket-like subsites that bind to and unveil the S1′–S4′ substrate binding preferences. We use molecular modeling to dock a model of the pro-peptide in the S1–S4 region and to dock collagen along the active site cleft. Opposite to the face harboring the catalytic serine, the enzyme binds to a calcium ion in a binding site created by backbone flipping. We use thermal unfolding to show that this metal ion provides structural stability. With no known host inhibitor identified thus far, this structure may hasten the progress of developing new therapeutic agents for diseases caused by pathogenic fungi.

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

The molecular basis of Human FN3K mediated phosphorylation of glycated substrates

Abstract Glycation, a non-enzymatic post-translational modification occurring on proteins, can be actively reversed via site-specific phosphorylation of the fructose-lysine moiety by FN3K kinase, to impact the cellular function of the target protein. A regulatory axis between FN3K and glycated protein targets has been associated with conditions like diabetes and cancer. However, the molecular basis of this relationship has not been explored so far. Here, we determined a series of crystal structures of HsFN3K in the apo-state, and in complex with different nucleotide analogs together with a sugar substrate mimic to reveal the features important for its kinase activity and substrate recognition. Additionally, the dynamics in sugar substrate binding during the kinase catalytic cycle provide important mechanistic insights into HsFN3K function. Our structural work provides the molecular basis for rational small molecule design targeting FN3K.

Science & Technology - Other Topics↗
DOE OSTI2021

Thermal conductivity tensor of β -1,3,5,7-tetranitro-1,3,5,7-tetrazoctane ( β -HMX) as a function of pressure and temperature

Here, we have used reverse non-equilibrium molecular dynamics (RNEMD) simulations to determine the full thermal conductivity tensor for the monoclinic high explosive crystal β-1,3,5,7-tetranitro-1,3,5,7-tetrazoctane (β-HMX). In order to do so for the monoclinic crystal, four directions for heat propagation are used. Effects of the temperature and pressure are investigated between 200 and 500 K and 0 and 5 GPa, respectively, which approximately covers the range where the β polymorph is stable. Simulations are carried out with the Smith–Bharadwaj non-reactive empirical potential [Smith and Bharadwaj, J. Phys. Chem. B 103, 3570 (1999)], which is known to reproduce well the thermo-elastic properties of HMX. Our results indicate that the thermal conductivity, κ, is highly anisotropic, with 36% difference between the two extreme values at 300 K and 0 GPa. A simple function is used to interpolate κ in the pressure-temperature regime considered in this study, which can be used in continuum models. The results from RNEMD simulations compare well with available experimental results from the literature and allow the determination of κ for any direction and temperature and pressure within or around the fitting interval.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Defluorination Mechanisms and Real-Time Dynamics of Per- and Polyfluoroalkyl Substances on Electrified Surfaces

Per- and polyfluoroalkyl substances (PFAS) are persistent environmental contaminants found in groundwater sources and a wide variety of consumer products. In recent years, electrochemical approaches for the degradation of these harmful contaminants have garnered a significant amount of attention due to their efficiency and chemical-free modular nature. However, these electrochemical processes occur in open, highly non-equilibrium systems, and a detailed understanding of PFAS degradation mechanisms in these promising technologies is still in its infancy. To shed mechanistic insight into these complex processes, we present the first constant-electrode potential (CEP) quantum calculations of PFAS degradation on electrified surfaces. These advanced CEP calculations provide new mechanistic details about the intricate electronic processes that occur during PFAS degradation in the presence of an electrochemical bias, which cannot be gleaned from conventional density functional theory calculations. We complement our CEP calculations with large-scale ab initio molecular dynamics simulations in the presence of an electrochemical bias to provide time scales for PFAS degradation on electrified surfaces. Taken together, our CEP-based quantum calculations provide critical reaction mechanisms for PFAS degradation in open electrochemical systems, which can be used to prescreen candidate material surfaces and optimal electrochemical conditions for remediating PFAS and other environmental contaminants.

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

Metastable Materials Discovery in the Age of Large-Scale Computation

Computational materials discovery has been successful in predicting novel, technologically relevant materials. However, it has remained focused almost exclusively on finding ground-state structures. Now that the lower-hanging fruit has been found in many fields of application, materials exploration is moving toward metastable materials: higher energy phases that are stable at practical time scales. Because of the challenges associated with predicting which phases are realistic, this class of materials has remained relatively unexplored, despite numerous examples of metastable structures with unmatched properties (e.g., diamond). This article highlights recent advances in developing computational and theoretical methods for predicting useful and realizable metastable materials. Topics discussed cover (1) the latest strategies for identifying potential metastable phases, (2) methodologies for assessing which phases can be realized experimentally, and (3) current approaches to estimate the lifetime of metastable materials.

amorphous materials↗
DOE OSTI2020

Community structure – Ecosystem function relationships in the Congo Basin methane cycle depend on the physiological scale of function

Belowground ecosystem processes can be highly variable and difficult to predict using microbial community data. Here, we argue that this stems from at least three issues: (a) complex covariance structure of samples (with environmental conditions or spatial proximity) can make distinguishing biotic drivers a challenge; (b) communities can control ecosystem processes through multiple mechanisms, making the identification of these controls a challenge; and (c) ecosystem function assessments can be broad in physiological scale, encapsulating multiple processes with unique microbially mediated controls. In this work, we test these assertions using methane (CH 4 )-cycling processes in soil samples collected along a wetland-to-upland habitat gradient in the Congo Basin. We perform our measurements of function under controlled laboratory conditions and statistically control for environmental covariates to aid in identifying biotic drivers. We divide measurements of microbial communities into four attributes (abundance, activity, composition, and diversity) that represent different forms of community control. Lastly, our process measurements differ in physiological scale, including broader processes (gross methanogenesis and methanotrophy) that involve more mediating groups, to finer processes (hydrogenotrophic methanogenesis and high-affinity CH 4 oxidation) with fewer mediating groups. We observed that finer scale processes can be more readily predicted from microbial community structure than broader scale processes. In addition, the nature of those relationships differed, with broad processes limited by abundance while fine-scale processes were associated with diversity and composition. These findings demonstrate the importance of carefully defining the physiological scale of ecosystem function and performing community measurements that represent the range of possible controls on ecosystem processes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Theoretical Modeling of Electrochemical Proton-Coupled Electron Transfer

Proton-coupled electron transfer (PCET) plays an essential role in a wide range of electrocatalytic processes. Here, a vast array of theoretical and computational methods have been developed to study electrochemical PCET. These methods can be used to calculate redox potentials and pKa values for molecular electrocatalysts, proton-coupled redox potentials and bond dissociation free energies for PCET at metal and semiconductor interfaces, and reorganization energies associated with electrochemical PCET. Periodic density functional theory can also be used to compute PCET activation energies and perform molecular dynamics simulations of electrochemical interfaces. Various approaches for maintaining a constant electrode potential in electronic structure calculations and modeling complex interactions in the electric double layer (EDL) have been developed. Theoretical formulations for both homogeneous and heterogeneous electrochemical PCET spanning the adiabatic, nonadiabatic, and solvent-controlled regimes have been developed and provide analytical expressions for the rate constants and current densities as functions of applied potential. The quantum mechanical treatment of the proton and inclusion of excited vibronic states have been shown to be critical for describing experimental data, such as Tafel slopes and potential dependent kinetic isotope effects. The calculated rate constants can be used as input to microkinetic models and voltammogram simulations to elucidate complex electrocatalytic processes.

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

Machine learning-accelerated path integral molecular dynamics simulations of reactive organic electrolytes

Hydrogen bonded electrolytes that exhibit accelerated proton transport via sequential reactive hops have drawn interest for their promise in clean energy applications. Molecular dynamics simulations of these electrolytes offer the opportunity to uncover microscopic mechanistic details that could be used to design and tune the properties of candidate electrolyte technologies. However, accurately modeling the proton transfer reactions and transport properties that give rise to high charge conductivites in these electrolytes proves computationally challenging because of the need to perform lengthy condensed phase simulations, treating both the electronic and nuclear degrees of freedom quantum mechanically. In this paper, we demonstrate that such a modeling task can be efficiently achieved with the use of density functional theory (DFT)-trained machine learning potentials (MLP) to accelerate path integral molecular dynamics (PIMD) simulations. We highlight the practical utility of this approach by using it to benchmark how closely PIMD simulations employing different DFT exchange–correlation functionals reproduce the composition-dependent densities, diffusion coefficients, and electrical conductivities of mixtures consisting of imidazole and levulinic acid. Even with the speedup afforded by our MLPs, PIMD simulations remain quite expensive. Furthermore, in order to render PIMD more computationally tractable, we introduce and benchmark the accuracy of a ring polymer contraction approach that leverages a computationally efficient short-range MLP to accelerate our PIMD simulations by an additional factor of four.

Chemical bonding↗
DOE OSTI2026

Raman and Terahertz Spectroscopy of Low-Frequency Chiral Phonons in Amino Acids

Chiral phonons are mirror-symmetric vibrations with nonzero angular momenta that correspond to twisting and rotational motions of multiple atoms. In chiral crystals, these include low-energy terahertz (THz)-range vibrations of the molecular segments involving dozens of atoms with energies sensitive to molecular chirality. Here, we present spectral signatures of chiral phonons in circularly polarized Raman optical activity (ROA) spectra from enantiomers of amino acid crystals. Along with complementary THz circular dichroism (TCD) measurements, our ROA data reveal several vibrational bands in enantiomers of valine, alanine, tyrosine, and proline between 30 and 150 cm –1 (∼1–4.5 THz) that exhibit opposite intensities. Density functional theory calculations confirm their assignment to twisting and shearing molecular motions. The simultaneous registration of chiral phonon modes by ROA and TCD demonstrates the necessity of these complementary techniques to identify complex mirror-asymmetric vibrational modes and offers new insights into their interactions with circularly polarized light.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Engineering the formation of spin-defects from first principles

The full realization of spin qubits for quantum technologies relies on the ability to control and design the formation processes of spin defects in semiconductors and insulators. We present a computational protocol to investigate the synthesis of point-defects at the atomistic level, and we apply it to the study of a promising spin-qubit in silicon carbide, the divacancy (VV). Our strategy combines electronic structure calculations based on density functional theory and enhanced sampling techniques coupled with first principles molecular dynamics. We predict the optimal annealing temperatures for the formation of VVs at high temperature and show how to engineer the Fermi level of the material to optimize the defect’s yield for several polytypes of silicon carbide. Our results are in excellent agreement with available experimental data and provide novel atomistic insights into point defect formation and annihilation processes as a function of temperature.

36 MATERIALS SCIENCE↗