SEARCH · Engineering Papers
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.
We develop expressions for electron density defined through the linear response for general density functional approximations, demonstrating results for orbital functionals and for many-body perturbation theory, and explore the connections to developments in DFT.
We develop the relationship between the strain derivative of the mean-inner potential and surface contributions to flexoelectricity, identifying the true surface-specific component of the flexoelectric response of finite samples. Density functional theory calculations on a range of experimentally observed, low energy SrTiO 3 , MgO, and Si surfaces demonstrate that the mean-inner potential and its contributions to flexoelectricity are sensitive to small differences in surface structure, chemistry, and adsorbates. We also introduce a method to estimate mean-inner potential contributions to flexoelectricity using electron scattering factors and use this approximation to predict total flexoelectric responses for a variety of insulators. Strategies to experimentally disentangle bulk and surface flexoelectric terms are also discussed.
The Bethe-Salpeter equation (BSE) based on GW quasiparticle levels is a successful approach for calculating the optical gaps and spectra of solids and also for predicting the neutral excitations of small molecules. We here present an all-electron implementation of the GW+BSE formalism for molecules, using numeric atom-centered orbital (NAO) basis sets. We present benchmarks for low-lying excitation energies for a set of small organic molecules, denoted in the literature as “Thiel’s set.” Literature reference data based on Gaussian-type orbitals are reproduced to about one millielectron-volt precision for the molecular benchmark set, when using the same GW quasiparticle energies and basis sets as the input to the BSE calculations. For valence correlation consistent NAO basis sets, as well as for standard NAO basis sets for ground state density-functional theory with extended augmentation functions, we demonstrate excellent convergence of the predicted low-lying excitations to the complete basis set limit. A simple and affordable augmented NAO basis set denoted “tier2+aug2” is recommended as a particularly efficient formulation for production calculations. We finally demonstrate that the same convergence properties also apply to linear-response time-dependent density functional theory within the NAO formalism.
Treating water as a linearly responding dielectric continuum on molecular length scales allows very simple estimates of the solvation structure and thermodynamics for charged and polar solutes. While this approach can successfully account for basic length and energy scales of ion solvation, computer simulations indicate not only its quantitative inaccuracies but also its inability to capture some basic and important aspects of microscopic polarization response. Here, we consider one such shortcoming, a failure to distinguish the solvation thermodynamics of cations from that of otherwise-identical anions, and we pursue a simple, physically inspired modification of the dielectric continuum model to address it. The adaptation is motivated by analyzing the orientational response of an isolated water molecule whose dipole is rigidly constrained. Its free energy suggests a Hamiltonian for dipole fluctuations that accounts implicitly for the influence of higher-order multipole moments while respecting constraints of molecular geometry. Finally, we propose a field theory with the suggested form, whose nonlinear response breaks the charge symmetry of ion solvation. An approximate variational solution of this theory, with a single adjustable parameter, yields solvation free energies that agree closely with simulation results over a considerable range of solute size and charge.
Human activities involving subsurface reservoirs—resource extraction, carbon and nuclear waste storage—alter thermal, mechanical, and chemical steady-state conditions in these systems. Because these systems exist at lithostatic pressures, even minor chemical changes can cause chemically assisted deformation. Therefore, understanding how chemical effects control geomechanical properties is critical to optimizing engineering activities. The grand challenge in predicting the effect of chemical processes on mechanical properties lays in the fact that these phenomena take place at molecular scales, while they manifest all the way to reservoir scales. To address this fundamental challenge, we investigated chemical effects on deformation in model and real systems spanning molecular- to centimeter scales. We used theory, experiment, molecular dynamics simulation, and statistical analysis to (1) identify the effect of simple reactions, such as hydrolysis, on molecular structures in interfacial regions of stressed geomaterials; (2) quantify chemical effects on the bulk mechanical properties, fracture and displacement for granular rocks and single crystals; (3) develop initial understanding of universal scaling for individual displacement events in layered geomaterials; and (4) develop analytic approximations for the single-chain mechanical response utilizing asymptotically correct statistical thermodynamic theory. Taken together, these findings advance the challenging field of chemo-mechanics.
Fabric, expressed by means of an evolving deviatoric fabric tensor F, plays a very important role in the anisotropic mechanical response of granular materials. The Anisotropic Critical State Theory (ACST) addresses fabric anisotropy by rendering dilatancy a function of F, in addition to other state variables. In this paper, 3D DEM is used to guide the specific grain-level definition of F, the formulation of its continuum evolution equation and its effect on anisotropic dilatancy within ACST. DEM provides stress-ratio and shear strain variations as input for ACST analytical calculations of evolving fabric tensor and dilatancy, which are then favourably compared with totally independent direct measurements of these quantities by DEM. Furthermore, dilatancy is shown to be strongly affected by the contact normal-based fabric tensor Fc, whose evolution is best described by a continuum equation within ACST that includes dilatancy and a quantity related to particle orientation-based fabric tensor Fp. The aforementioned favourable comparison of the results for fabric tensor and dilatancy obtained independently by ACST and DEM, confirms the validity of the core framework of ACST irrespective of any constitutive model that addresses the deviatoric stress-strain relations.
Abstract not provided.
A classic challenge in chemical sensing is selectivity. Metal–organic frameworks (MOFs) are an exciting class of materials because they can be tuned for selective chemical adsorption. Adsorption events trigger work-function shifts, which can be detected with a chemical-sensitive field-effect transistor (power ≈microwatts). Here, several case studies were used towards generalizing the sensing mechanism, ultimately towards our metal-centric hypothesis. HKUST-1 was used as a proof-of-principle humidity sensor. The response is thickness independent, meaning the response is surface localized. ZIF-8 is demonstrated to be an NO 2 -sensing material, and the response is dominated by adsorption at metal sites. Finally, MFM-300(In) shows how standard hard–soft acid–base theory can be used to qualitatively predict sensor responses. This paper sets the groundwork for using the tunability of metal–organic frameworks for chemical sensing with distributed, scalable devices.
The average ionization state is a critical parameter in plasma models for charged particle transport, equations of state, and optical response. The dynamical or nonadiabatic Born effective charge (NBEC), calculated via first principles time-dependent density-functional theory, provides exact ionic partitioning of bulk electron response for both metallic and insulating materials. The NBEC can be transformed into a “group conductivity,” i.e., the electron conductivity ascribed to a subset of ions. We show that for disordered metallic systems, such as warm dense matter (WDM) and hot dense plasma, the static limit of the NBEC is different from the average ionization states, but that the ionization state can be extracted from the group conductivity even in mixed systems. Here we demonstrate this approach using a set of archetypical examples, including cold and warm aluminium, low- and high-density WDM carbon, and a WDM carbon-beryllium-hydrogen mixture.
We present a real-space method for computing the random phase approximation (RPA) correlation energy within Kohn–Sham density functional theory, leveraging the low-rank nature of the frequency-dependent density response operator. In particular, we employ a cubic-scaling formalism based on density functional perturbation theory that circumvents the calculation of the response function matrix, instead relying on the ability to compute its product with a vector through the solution of the associated Sternheimer linear systems. We develop a large-scale parallel implementation of this formalism using the subspace iteration method in conjunction with the spectral quadrature method while employing the Kronecker product-based method for the application of the Coulomb operator and the conjugate orthogonal conjugate gradient method for the solution of the linear systems. We demonstrate convergence with respect to key parameters and verify the method’s accuracy by comparing with plane-wave results. We show that the framework achieves good strong scaling to many thousands of processors, reducing the time to solution for a lithium hydride system with 128 electrons to around 150 s on 4608 processors.
The extradimensional axion solves the strong C P problem while largely circumventing the quality problem that plagues its four-dimensional counterparts. Such high quality can be clearly understood in terms of the generalized global symmetries of the higher-dimensional theory. We emphasize that an electric one-form symmetry is entirely responsible for protecting the potential of axions arising from 5D gauge theories and use this to systematically characterize the extradimensional axion quality problem. We identify three ways of breaking this one-form symmetry to generate an axion potential: adding electrically charged matter, gauging a magnetic higher-form symmetry, and turning on an Adler-Bell-Jackiw anomaly. In the latter case, we identify new ways of generating an axion potential via extradimensional magnetic monopoles. We also discuss how the axion is modified if the protective electric one-form symmetry is itself gauged. Finally, we relate these effects to gravitational expectations for the quality problem via generalized weak gravity conjectures. The clarity that generalized symmetries bring to the extradimensional axion quality problem highlights their relevance to particle phenomenology. Published by the American Physical Society 2025
This research project in theoretical condensed-matter physics aims to study the collective charge and spin excitations of itinerant quasi-two-dimensional electronic systems such as doped semiconductor nanostructures and carbon-based materials, in the presence of impurities and spin-orbit coupling. In previous work under DOE support, we had discovered an interesting interplay between many-body effects and Rashba and Dresselhaus spin-orbit interactions: contrary to many situation in spintronics where spin-orbit interactions lead to loss of spin coherence, collective spin modes are robust against such coherence losses, and behave as macroscopic quantum objects with a collective magnetic moment precessing about an effective spin-orbit magnetic field. Experimental data for GaAs and CdMnTe quantum wells, obtained by our collaborators at the University Pierre et Marie Curie in Paris, confirmed our predictions and suggested that the observed effect is universal. However, a full theoretical explanation of the effect is still missing, and a wealth of new experimental data remains to be analyzed. We therefore plan to address the following three interrelated objectives: 1. Theory of collective modes in chiral systems. We will explore the various ways in which spin-orbit interactions affect collective spin-density and spin-flip excitations in semiconductor quantum wells. The goal will be to develop a universal formalism for chiral collective modes, with a particular emphasis on analytic results for dispersion relations, which will be compared with experimental results from inelastic light scattering. We will calculate the linewidths of spin plasmons and spin-flip waves due to impurity scattering and spin Coulomb drag from first principles. 2.Many-body effects in noncollinear magnetism. We will describe electronic many-body effects using time-dependent density-functional theory (TDDFT) in the linear-response regime. Quantum wells with spin-orbit coupling are noncollinear itinerant magnetic systems at the crossover between two and three dimensions, which poses subtle challenges for standard exchange-correlation (xc) functionals such as the local-density approximation (LDA). We will develop new orbital-dependent xc functionals for noncollinear magnetic systems by generalizing the so-called STLS approach. 3. Beyond single semiconductor quantum wells. We will consider quantum well bilayer or multilayer systems with spin-orbit coupling, where the coupling between the quasi-2D electron systems leads to acoustic and optical modes (in phase and out of phase). How are these modes affected if the two layers have the same or different chirality? Plasmon modes have been widely studied in graphene and other topological materials, and we will study collective spin excitations in graphene and other carbon-based materials in the presence of spin-orbit coupling and magnetic fields. The potential impact of the proposed research will be both on a fundamental and a practical level. A universal theory of chiral collective modes will advance our understanding of the interplay between many-body effects and spin-orbit interactions in a wide variety of materials and systems. A generalized STLS approach for noncollinear magnetism will be of wide interest beyond the systems studied here. This work may also point towards new ways of manipulating and controlling spin waves in low-dimensional itinerant electron systems: this could be used to encode, transport and process information, which would be of interest to the field of magnon spintronics.
An electrolyte transport theory connects its transport properties, evolution of spatiotemporal fields (e.g., concentration), and corresponding macroscopic current and voltage responses. Given this interconnection, the transport properties are typically inferred by analyzing the macroscopic response through the lens of a chosen electrolyte transport theory. Unfortunately, the same macroscopic measurements can be analyzed by using different theories to arrive at seemingly dissimilar properties that are inconsistent with each other. Here, we offer a resolution to this dilemma by analyzing multiple spatiotemporal (i.e., operando) fields for a given electrolyte. We show that predictive transport theory can analyze operando fields to estimate underlying transport properties and subsequently predict another operando field. A theory passing this test identifies meaningful transport properties such that the electrolyte behavior can be accurately predicted over a wide range of excitations, and is critical to property-based electrolyte screening and electrolyte discovery efforts.
Solar-induced chlorophyll fluorescence (SIF) has shown great potential in estimating gross primary production (GPP). However, their quantitative relationship is not invariant, which undermines the reliability of empirical SIF-based GPP estimation at fine spatiotemporal scales, especially under extreme conditions. In this study, we developed a parsimonious mechanistic model for SIF-based GPP estimation in evergreen needle forests (ENF) by employing the Mechanistic Light Response framework and Eco-Evolutionary theory to describe the light and dark reactions during photosynthesis, respectively. Specifically, we found that considering the seasonal variation in a key parameter of the MLR framework, the maximum photochemical efficiency of photosystem II (Φ PSIImax ), can avoid the GPP overestimation in winter and early spring due to the relatively low environmental sensitivity of SIF. Compared to the estimates from other benchmark models, our GPP estimates were closer to the 1: 1 line and had higher accuracy (average R 2 = 0.86, RMSE=1.99 μmol m -2 s -1 ) across sites. Furthermore, the changes in the relationship between SIF and J (refers to the electron transport rate) contribute a lot to the dynamic SIF–GPP relationship in this study, while the J–GPP relationship is less variant when the temperature drops. Further, the seasonal variation in the SIF–J relationship, especially the reduction in its slope at low temperatures, is found largely explained by the Φ PSIImax . These results indicate the importance of the uncertainty caused by the variation in the SIF–J relationship for SIF-based GPP estimation, and the consideration of changes in Φ PSIImax under extreme conditions (such as severe winter in this study) is crucial for the improvement of GPP estimation via SIF.
The circularly polarized photogalvanic effect (CPGE) is studied in chiral Weyl semimetals with short-range quenched disorder. Without disorder, the topological properties of chiral Weyl semimetals lead to quantization of the CPGE, which is a second-order optical response. Furthermore, using a combination of diagrammatic perturbation theory in the continuum and exact numerical calculations via the kernel polynomial method on a lattice model, we show that disorder perturbatively destabilizes the quantization of the CPGE.
Burst pressure is one of the critical strength parameters used in the design and operation of pressure vessels because it represents the maximum pressure that a vessel can withstand before failing. Historically, the Barlow formula was used as a design base for estimating burst pressure. However, it does not consider the plastic flow response for ductile steels and is applicable only to thin-walled cylinders (i.e., the diameter to thickness ratio D/t ≥ 20). A new multiaxial plastic yield theory was developed to consider the plastic flow response, and the associated theoretical (i.e., Zhu–Leis) solution of burst pressure was obtained and has gained extensive applications in the pipeline industry because it was validated by different full-scale burst test datasets for large-diameter, thin-walled pipelines in a variety of steel grades from Grade B to X120. The Zhu–Leis flow theory of plasticity was recently extended to thick-walled pressure vessels, and the associated exact flow solution of burst pressure was obtained and is applicable to both thin and thick-walled cylindrical shells. Many full-scale burst tests are available for thin-walled line pipes in the pipeline industry, but limited pressure burst tests exist for thick-walled vessels. To validate the newly developed exact solutions of burst pressure for thick-walled cylinders, this paper conducts a series of burst pressure tests on small-diameter, thick-walled pipes. In particular, six burst tests are carried out for three thick-walled pipes in Grade B carbon steel. These pipes have a nominal diameter of 2.375 inches (60.33 mm) and three nominal wall thicknesses of 0.154, 0.218, and 0.344 inches (3.91, 5.54, and 8.74 mm), leading to D/t = 15.4, 10.9, and 6.9, respectively. With the burst test data, comparisons show that the Zhu–Leis flow solution of burst pressure matches well the burst test data for thick-walled pipes. Thus, these burst tests validate the accuracy of the Zhu–Leis flow solution of burst pressure for thick-walled cylindrical vessels.
The electronic structure, lattice dynamics, and superconducting properties of SrHg and BaHg intermetallic compounds are investigated using plane wave pseudo-potential and linear response methods based on density functional theory. The ground state properties such as lattice parameters are in good agreement with previously reported results. The positive phonon frequencies reveal the stability of both compounds in B{sub 2} (CsCl-type) structure. The metallic nature of SrHg and BaHg is confirmed by the electronic band structure as well as phonon dispersion curves. The superconducting behavior described in terms of Eliashberg spectral function (α{sup 2}F(ω)) and electron-phonon coupling constant (λ). The calculated values of superconducting transition temperatures (T{sub c}) show good agreement with available experimental data.
Scientific and engineering advances are primarily driven by multi-tier conceptual constructs and conditional theoretical frameworks. The theories allow predictions of hypothetical system responses, given a set of approximate conditions (ranges of applicability) imposed on latent parameters that cannot be measured directly. Learning to estimate the latent variables (Latent Learning) helps to pinpoint the anticipated range-edge anomalies and improves the confidence in interpretation, interpolation and extrapolation of limited experimental data. Due to high dimensionality and extreme non-linearity of the materials science problems, very large datasets are typically required for conventional data-driven model development. The vital experimental data collection, particularly on microstructural phases, is very challenging, which makes it difficult to compile a high-quality database. Incorporation of the domain knowledge into the computational graph structure, initialization and optimization processes presents a viable mechanism for developing accurate models, with limited datasets. Furthermore, this study successfully utilized the approach to build the Deep Freeze Graph (DeepFreG) by mapping known causality relationships and by digitizing empirical domain knowledge for Latent Learning (LL), with specific applications in materials science.