Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “affine approximation”

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 37 records · Page 2

Franck-Condon electron emission from polar semiconductor photocathodes

An analytical formulation of (optical-)phonon-mediated and momentum-resonant Franck-Condon emission of photoexcited electrons from polar semiconductors is shown to be very consistent with (i) the observed emission properties of a cesiated Ga⁢As⁢(001) photocathode at 808 nm [J. Phys. D: Appl. Phys. 54, 205301 (2021)] and (ii) the measured spectral emission properties of a Ga⁢N(0001) photocathode from just below its bandgap energy to 5 eV. The theoretical analysis in the parabolic band approximation predicts the form of both the quantum efficiency and mean transverse energy of photoemission as a function of the photocathode’s electron affinity and the electron temperature in the vicinity of its emission face. The good agreement between theory and experimental data also suggests that sub-10-nm rms surface roughness effects are not significant for polar semiconductor photocathodes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Stochastic DC optimal power flow with reserve saturation

We propose an optimization framework for stochastic optimal power flow with uncertain loads and renewable generator capacity. Our model follows previous work in assuming that generator outputs respond to load imbalances according to an affine control policy, but introduces a model of saturation of generator reserves by assuming that when a generator’s target level hits its limit, it abandons the affine policy and produces at that limit. This is a particularly interesting feature in models where wind power plants, which have uncertain upper generation limits, are scheduled to provide reserves to balance load fluctuations. Here, the resulting model is a nonsmooth nonconvex two-stage stochastic program, and we use a stochastic approximation method to find stationary solutions to a smooth approximation. Computational results on 6-bus and 118-bus test instances demonstrate that by considering the effects of saturation, our model can yield solutions with lower expected generation costs (at the same target line violation probability level) than those obtained from a model that enforces the affine policy to stay within generator limits with high probability.

17 WIND ENERGY↗

Enhanced Binding of SARS-CoV-2 Spike Protein to Receptor by Distal Polybasic Cleavage Sites

The receptor-binding domain (RBD) of the SARS-CoV-2 spike protein plays a crucial role in binding the human cell receptor ACE2 that is required for viral entry. Many studies have been conducted to target the structures of RBD–ACE2 binding and to design RBD-targeting vaccines and drugs. Nevertheless, mutations distal from the SARS-CoV-2 RBD also impact its transmissibility and antibody can target non-RBD regions, suggesting the incomplete role of the RBD region in the spike protein–ACE2 binding. Here, in order to elucidate distant binding mechanisms, we analyze complexes of ACE2 with the wild-type spike protein and with key mutants via large-scale all-atom explicit solvent molecular dynamics simulations. We find that though distributed approximately 10 nm away from the RBD, the SARS-CoV-2 polybasic cleavage sites enhance, via electrostatic interactions and hydration, the RBD–ACE2 binding affinity. A negatively charged tetrapeptide (GluGluLeuGlu) is then designed to neutralize the positively charged arginine on the polybasic cleavage sites. We find that the tetrapeptide GluGluLeuGlu binds to one of the three polybasic cleavage sites of the SARS-CoV-2 spike protein lessening by 34% the RBD–ACE2 binding strength. This significant binding energy reduction demonstrates the feasibility to neutralize RBD–ACE2 binding by targeting this specific polybasic cleavage site. Our work enhances understanding of the binding mechanism of SARS-CoV-2 to ACE2, which may aid the design of therapeutics for COVID-19 infection.

60 APPLIED LIFE SCIENCES↗

Exploring Camelina sativa lipid metabolism regulation by combining gene co-expression and DNA affinity purification analyses

Camelina (Camelina sativa) is an annual oilseed plant that is gaining momentum as a biofuel cover crop. Understanding gene regulatory networks is essential to deciphering plant metabolic pathways, including lipid metabolism. Furthermore, we take advantage of a growing collection of gene expression datasets to predict transcription factors (TFs) associated with the control of Camelina lipid metabolism. We identified approximately 350 TFs highly co-expressed with lipid-related genes (LRGs). These TFs are highly represented in the MYB, AP2/ERF, bZIP, and bHLH families, including a significant number of homologs of well-known Arabidopsis lipid and seed developmental regulators. After prioritizing the top 22 TFs for further validation, we identified DNA-binding sites and predicted target genes for 16 out of the 22 TFs tested using DNA affinity purification followed by sequencing (DAP-seq). Enrichment analyses of targets supported the co-expression prediction for most TF candidates, and the comparison to Arabidopsis revealed some common themes, but also aspects unique to Camelina. Within the top potential lipid regulators, we identified CsaMYB1, CsaABI3AVP1-2, CsaHB1, CsaNAC2, CsaMYB3, and CsaNAC1 as likely involved in the control of seed fatty acid elongation and CsaABI3AVP1-2 and CsabZIP1 as potential regulators of the synthesis and degradation of triacylglycerols (TAGs), respectively. Altogether, the integration of co-expression data and DNA-binding assays permitted us to generate a high-confidence and short list of Camelina TFs involved in the control of lipid metabolism during seed development.

59 BASIC BIOLOGICAL SCIENCES↗

Effect of activation temperature on quantum efficiency and lifetime of NEA truncated nanocone array GaAs photocathode

This study investigates the quantum efficiency (QE) and operational lifetime of a negative electron affinity GaAs truncated nanocone array (TNCA) photocathode benchmarked against a conventional flat GaAs photocathode under varying activation temperatures. The TNCA structure demonstrated a QE of up to 13.6% at 590 nm with room temperature (RT) activation—approximately 1.5 times higher than its flat counterpart. This enhancement is due to Mie resonance effects within the nanostructure, as confirmed by finite-difference time-domain simulations. Moreover, the TNCA photocathode exhibits significantly extended charge lifetime, with enhancement factors of ∼6.1 and ∼19.8 under RT and 50 °C activations, respectively. These gains are primarily attributed to increased effective surface area and optimized dipole layer formation at elevated temperatures. In addition, shorter excitation wavelengths further contribute to lifetime improvements. These findings underscore the TNCA GaAs photocathode’s potential as a high QE, long lifetime electron source for many large-scale electron accelerators.

Cs-NF3 activation↗

Formation of a Strong Heterogeneous Aluminum Lewis Acid on Silica

Abstract Al(OC(CF 3 ) 3 )(PhF) reacts with silanols present on partially dehydroxylated silica to form well‐defined ≡SiOAl(OC(CF 3 ) 3 ) 2 (O(Si≡) 2 ) ( 1 ). 27 Al NMR and DFT calculations with a small cluster model to approximate the silica surface show that the aluminum in 1 adopts a distorted trigonal bipyramidal coordination geometry by coordinating to a nearby siloxane bridge and a fluorine from the alkoxide. Fluoride ion affinity (FIA) calculations follow experimental trends and show that 1 is a stronger Lewis acid than B(C 6 F 5 ) 3 and Al(OC(CF 3 ) 3 )(PhF) but is weaker than Al(OC(CF 3 ) 3 ) and i Pr 3 Si + . Cp 2 Zr(CH 3 ) 2 reacts with 1 to form [Cp 2 ZrCH 3 ][≡SiOAl(OC(CF 3 ) 3 ) 2 (CH 3 )] ( 3 ) by methide abstraction. This reactivity pattern is similar to reactions of organometallics with the proposed strong Lewis acid sites present on Al 2 O 3 .

Samudrala, Kavyasripriya K.↗

Formation of a Strong Heterogeneous Aluminum Lewis Acid on Silica

Abstract Al(OC(CF 3 ) 3 )(PhF) reacts with silanols present on partially dehydroxylated silica to form well‐defined ≡SiOAl(OC(CF 3 ) 3 ) 2 (O(Si≡) 2 ) ( 1 ). 27 Al NMR and DFT calculations with a small cluster model to approximate the silica surface show that the aluminum in 1 adopts a distorted trigonal bipyramidal coordination geometry by coordinating to a nearby siloxane bridge and a fluorine from the alkoxide. Fluoride ion affinity (FIA) calculations follow experimental trends and show that 1 is a stronger Lewis acid than B(C 6 F 5 ) 3 and Al(OC(CF 3 ) 3 )(PhF) but is weaker than Al(OC(CF 3 ) 3 ) and i Pr 3 Si + . Cp 2 Zr(CH 3 ) 2 reacts with 1 to form [Cp 2 ZrCH 3 ][≡SiOAl(OC(CF 3 ) 3 ) 2 (CH 3 )] ( 3 ) by methide abstraction. This reactivity pattern is similar to reactions of organometallics with the proposed strong Lewis acid sites present on Al 2 O 3 .

Samudrala, Kavyasripriya K.↗

Interaction of Th with H 0/–/+ : Combined Experimental and Theoretical Thermodynamic Properties

In this work, high-level electronic structure calculations of the lowlying energy electronic states for ThH, ThH – , and ThH + are reported and compared to experimental measurements. The inclusion of spin–orbit coupling is critical to predict the ground-state ordering as inclusion of spin–orbit switches the coupled-cluster CCSD(T) ordering of the two lowest energy states for ThH and ThH + . At the multireference spin–orbit SO-CASPT2 level, the ground states of ThH, ThH – , and ThH + are predicted to be the 2 Δ 3/2 , 3 Φ 2 , and 3 Δ 1 states, respectively. The adiabatic electron affinity is calculated to be 0.820 eV, and the vertical detachment energy is calculated to be 0.832 eV in comparison to an experimental value of 0.87 ± 0.02 eV. The observed ThH – photoelectron spectrum has many transitions, which approximately correlate with excitations of Th + and/or Th. The adiabatic ionization energy of ThH including spin–orbit corrections is calculated to be 6.181 eV. The natural bond orbital results are consistent with a significant contribution of the Th + H – ionic configuration to the bonding in ThH. The bond dissociation energies for ThH, ThH – , and ThH + using the Feller–Peterson–Dixon approach were calculated to be similar for all three molecules and lie between 259 and 280 kJ/mol.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Sulfonyl PIM–1: A diverse separation membrane with dilation resistance

A sulfone modified variant (soPIM–1) of the first polymer of intrinsic microporosity has been studied through molecular simulations to analyze its applicability for adsorption–based separations of diverse nonpolar, quadrupolar, and dipolar adsorbate species. Single component gas phase adsorption isotherms of each adsorbate are provided. Here, the adsorptive properties of soPIM–1 are reported from the application of two methods: (a) traditional grand canonical Monte Carlo (GCMC) simulations and (b) a combination of Monte Carlo and molecular dynamics (GCMC/MD) techniques, which accounts for sorption–induced polymer dynamics. The calculated isosteric heats of adsorption are compared to those from the parent PIM–1 structure and reveals increased CO2 affinity with relatively no change in hydrocarbon affinity. Moreover, to quantify soPIM–1's resistance to dilation the evolution of the microporous structure across the adsorption isotherm loading pressures has been evaluated. Relatedly, nonequilibrium MD simulations have been utilized to determine that soPIM–1 is approximately twice as stiff as PIM–1.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A contorted nanographene shelter

Nanographenes have kindled considerable interest in the fields of materials science and supramolecular chemistry as a result of their unique self-assembling and optoelectronic properties. Encapsulating the contorted nanographenes inside artificial receptors, however, remains challenging. Herein, we report the design and synthesis of a trigonal prismatic hexacationic cage, which has a large cavity and adopts a relatively flexible conformation. It serves as a receptor, not only for planar coronene, but also for contorted nanographene derivatives with diameters of approximately 15 Å and thicknesses of 7 Å. A comprehensive investigation of the host-guest interactions in the solid, solution and gaseous states by experimentation and theoretical calculations reveals collectively an induced-fit binding mechanism with high binding affinities between the cage and the nanographenes. Notably, the photostability of the nanographenes is improved significantly by the ultrafast deactivation of their excited states within the cage. Encapsulating the contorted nanographenes inside the cage provides a noncovalent strategy for regulating their photoreactivity.

36 MATERIALS SCIENCE↗

The performance of ensemble-based free energy protocols in computing binding affinities to ROS1 kinase

Optimization of binding affinities for compounds to their target protein is a primary objective in drug discovery. Herein we report on a collaborative study that evaluates a set of compounds binding to ROS1 kinase. We use ESMACS (enhanced sampling of molecular dynamics with approximation of continuum solvent) and TIES (thermodynamic integration with enhanced sampling) protocols to rank the binding free energies. The predicted binding free energies from ESMACS simulations show good correlations with experimental data for subsets of the compounds. Consistent binding free energy differences are generated for TIES and ESMACS. Although an unexplained overestimation exists, we obtain excellent statistical rankings across the set of compounds from the TIES protocol, with a Pearson correlation coefficient of 0.90 between calculated and experimental activities.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A weighted shifted boundary method for immersed moving boundary simulations of Stokes' flow

The Shifted Boundary Method (SBM) belongs to the class of unfitted (or immersed, or embedded) finite element methods, and relies on reformulating the original boundary value problem over a surrogate (approximate) computational domain. The surrogate domain is constructed so as to avoid cut cells and the associated problematic implementation and numerical integration issues. Accuracy is maintained by modifying the original boundary conditions using Taylor expansions: hence the name of the method, that shifts the location and values of the boundary conditions. Here, in this article, we extend the SBM to the simulation of incompressible Stokes flow, by appropriately weighting its variational form with the elemental volume fraction of active fluid. This approach allows to drastically reduce spurious pressure oscillations in time, which are produced if the total volume of active fluid were to change abruptly over a time step. The proposed Weighted SBM (W-SBM) exactly preserves states of hydrostatic equilibrium, and induces small mass and momentum conservation errors, which converge as the grid is refined. This is in analogy to cutFEMs and related unfitted approaches, which rely on an affine representation of cut boundaries. We demonstrate the robustness and accuracy of the proposed method with an extensive suite of two-dimensional tests.

97 MATHEMATICS AND COMPUTING↗

RrA, an enzyme from Rhodospirillum rubrum , is a prototype of a new family of short‐chain L‐asparaginases

Abstract L‐Asparaginases (ASNases) catalyze the hydrolysis of L‐Asn to L‐Asp and ammonia. Members of the ASNase family are used as drugs in the treatment of leukemia, as well as in the food industry. The protomers of bacterial ASNases typically contain 300–400 amino acids (typical class 1 ASNases). In contrast, the chain of ASNase from Rhodospirillum rubrum , reported here and referred to as RrA, consists of only 172 amino acid residues. RrA is homologous to the N‐terminal domain of typical bacterial class 1 ASNases and exhibits millimolar affinity for L‐Asn. In this study, we demonstrate that RrA belongs to a unique family of cytoplasmic, short‐chain ASNases (scASNases). These proteins occupy a distinct region in the sequence space, separate from the regions typically assigned to class 1 ASNases. The scASNases are present in approximately 7% of eubacterial species, spanning diverse bacterial lineages. They seem to be significantly enriched in species that encode for more than one class 1 ASNase. Here, we report biochemical, biophysical, and structural properties of RrA, a member of scASNases family. Crystal structures of the wild‐type RrA, both with and without bound L‐Asp, as well as structures of several RrA mutants, reveal topologically unique tetramers. Moreover, the active site of one protomer is complemented by two residues (Tyr21 and Asn26) from another protomer. Upon closer inspection, these findings clearly outline scASNases as a stand‐alone subfamily of ASNases that can catalyze the hydrolysis of L‐Asn to L‐Asp despite the lack of the C‐terminal domain that is present in all ASNases described structurally to date.

59 BASIC BIOLOGICAL SCIENCES↗

Investigation and Analysis of Dolomite Dissolution in Variable Strength Systems Relevant to the WIPP

The Waste Isolation Pilot Plant (WIPP) was opened in 1999 as a solution to the long-term disposal of the nation's nuclear defense-generated transuranic (TRU) waste. The salt beds were recommended as the ideal location because they are free from flowing fresh water, easily mined, impermeable, and geologically stable. However, for performance assessment modeling, release scenarios must be considered. The most likely release scenario was determined to be human intrusion leading to direct and/or long-term brine release (US DOE 1995, US DOE 1996, Perkins et al., 1999). The focus of this research is the Rustler formation because it is located above the WIPP and is the most transmissive layer (Perkins et al., 1999). Objective: To determine the effect of chelating ligands ( e.g. EDTA) on dissolution of dolomite in high ionic strength systems. Do brines impact dolomite dissolution and how is it affected by EDTA binding with cations? By varying the ionic strengths, different brine conditions can be simulated similar to those in and around the Culebra member and evaluate trends of dissolved Ca and Mg in the aqueous phase from the quantified data. Conclusions: <2% of dolomite dissolved over one week at pH 8.5. Greater dissolution occurred at higher ionic strength and in the presence of NaCl with approximately 0.4 vs. 0.9% of Mg removed from dolomite in 0.1 and 1.0 M Na. Further, 0.3 vs. 0.9% of Mg was removed from dolomite in 1.0 M IS CaCl{sub 2} versus NaCl. EDTA increased dissolution of dolomite at low ionic strength due to its strong complexation affinity for Ca/Mg with slightly greater removal of Ca than Mg. Significant differences were not observed at high ionic strength due to the presence of large concentrations of competing cations (Na, Ca, Mg). Slightly higher concentrations of Ca were removed from dolomite than Mg after one week, suggesting either incongruent dissolution or release through ion exchange of Ca from trace mineral impurities. Further, when EDTA was present, it may have increased Ca removal due to its greater complexation for Ca than Mg. Environmental Implications: High ionic strength brines and low ionic strength solutions with organic ligands increase dissolution of dolomite. Due to the long half-lives of radionuclides in the TRU waste, it may be important to consider the long term impact of dolomite dissolution for potential release scenarios.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

A comparison of QTP functionals against coupled-cluster methods for EAs of small organic molecules

EA-EOM-CCSD electron affinities and LUMO energies of various Kohn–Sham density functional theory (DFT) methods are calculated for an a priori IP benchmark set of 64 small, closed-shell molecules. The purpose of these calculations was to investigate whether the QTP KS-DFT functionals can emulate EA-EOM-CC with only a mean-field approximation. We show that the accuracy of DFT—relative to CCSD—improves significantly when elements of correlated orbital theory are introduced into the parameterization to define the QTP family of functionals. In particular, QTP(02), which has only a single range separation parameter, provides results accurate to a MAD of <0.15 eV for the whole set of 64 molecules compared to EA-EOM-CCSD, far exceeding the results from the non-QTP family of density functionals.

Chemistry↗

Electron-Affinity Time-Dependent Density Functional Theory: Formalism and Applications to Core-Excited States

The lack of particle-hole attraction and orbital relaxation within time-dependent density functional theory (TDDFT) lead to extreme errors in the prediction of K-edge X-ray absorption spectra (XAS). We derive a linear-response formalism that uses optimized orbitals of the n - 1-electron system as the reference, building orbital relaxation and a proper hole into the initial density. Here, our approach is an exact generalization of the static-exchange approximation that ameliorates the particle-hole interaction error associated with the adiabatic approximation and reduces errors in TDDFT XAS by orders of magnitude. With a statistical performance of just 0.5 eV root-mean-square error and the same computational scaling as TDDFT under the core-valence separation approximation, we anticipate that this approach will be of great utility in XAS calculations of large systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Robust and Simple ADMM Penalty Parameter Selection

We present a new method for online selection of the penalty parameter for the alternating direction method of multipliers (ADMM) algorithm. ADMM is a widely used method for solving a range of optimization problems, including those that arise in signal and image processing. In its standard form, ADMM includes a scalar hyperparameter, known as the penalty parameter, which usually has to be tuned to achieve satisfactory empirical convergence. In this work, we develop a framework for analyzing the ADMM algorithm applied to a quadratic problem as an affine fixed point iteration. Using this framework, we develop a new method for automatically tuning the penalty parameter by detecting when it has become too large or small. We analyze this and several other methods with respect to their theoretical properties, i.e., robustness to problem transformations, and empirical performance on several optimization problems. Our proposed algorithm is based on a theoretical framework with clear, explicit assumptions and approximations, is theoretically covariant/invariant to problem transformations, is simple to implement, and exhibits competitive empirical performance.

42 ENGINEERING↗

Structural basis of the American mink ACE2 binding by Y453F trimeric spike glycoproteins of SARS‐CoV‐2

Abstract Severe Acute Respiratory Syndrome Coronavirus 2 (SARS‐CoV‐2) enters the host cell by binding to angiotensin‐converting enzyme 2 (ACE2). While evolutionarily conserved, ACE2 receptors differ across various species and differential interactions with Spike (S) glycoproteins of SARS‐CoV‐2 viruses impact species specificity. Reverse zoonoses led to SARS‐CoV‐2 outbreaks on multiple American mink (Mustela vison) farms during the pandemic and gave rise to mink‐associated S substitutions known for transmissibility between mink and zoonotic transmission to humans. In this study, we used bio‐layer interferometry (BLI) to discern the differences in binding affinity between multiple human and mink‐derived S glycoproteins of SARS‐CoV‐2 and their respective ACE2 receptors. Further, we conducted a structural analysis of a mink variant S glycoprotein and American mink ACE2 (mvACE2) using cryo‐electron microscopy (cryo‐EM), revealing four distinct conformations. We discovered a novel intermediary conformation where the mvACE2 receptor is bound to the receptor‐binding domain (RBD) of the S glycoprotein in a “down” position, approximately 34° lower than previously reported “up” RBD. Finally, we compared residue interactions in the S‐ACE2 complex interface of S glycoprotein conformations with varying RBD orientations. These findings provide valuable insights into the molecular mechanisms of SARS‐CoV‐2 entry.

Virology↗