Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Linear response”

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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Fingerprinting fractons with pump-probe spectroscopy

We demonstrate how pump-probe techniques enable specific spectroscopic diagnostics of fracton phases of matter by studying lineon-planon braiding in the paradigmatic X-cube model. Our discussion builds on previous works probing anyonic exchange statistics in conventional spin liquids, but the extension to fracton phases reveals qualitatively different phenomena due to the restricted mobility of fractionalized excitations. A key feature is that nearby planons can form an emergent bound state by accessing different planes via alternative pairing configurations. We show that this bound state qualitatively modifies the long-time linear and nonlinear responses, leading to an asymptotic linear-in-𝑡 behavior for the pump-probe signal 𝜒 𝑍𝑍𝑋 . By contrast, swapping the pump and probe polarizations produces a nonlinear response 𝜒 𝑋𝑋𝑍 that is asymptotically 𝑡 independent. This asymmetry reflects the fact that the two species of fractionalized excitations live in different spatial dimensions. Thus the pump-probe signals studied here are sensitive to (i) nontrivial braiding statistics in three dimensions, (ii) the existence of bound states among fractionalized excitations, and (iii) the one-dimensional mobility of lineons. Furthermore, our results therefore provide spectroscopic signatures that distinguish fracton phases from conventional topologically ordered spin liquids.

braiding↗
DOE OSTI2024

Human–mouse comparison of the multistage nature of radiation carcinogenesis in a mathematical model

Abstract Mouse models are vital for assessing risk from environmental carcinogens, including ionizing radiation, yet the interspecies difference in the dose response precludes direct application of experimental evidence to humans. Herein, we take a mathematical approach to delineate the mechanism underlying the human–mouse difference in radiation‐related cancer risk. We used a multistage carcinogenesis model assuming a mutational action of radiation to analyze previous data on cancer mortality in the Japanese atomic bomb survivors and in lifespan mouse experiments. Theoretically, the model predicted that exposure will chronologically shift the age‐related increase in cancer risk forward by a period corresponding to the time in which the spontaneous mutational process generates the same mutational burden as that the exposure generates. This model appropriately fitted both human and mouse data and suggested a linear dose response for the time shift. The effect per dose decreased with increasing age at exposure similarly between humans and mice on a per‐lifespan basis (0.72‐ and 0.71‐fold, respectively, for every tenth lifetime). The time shift per dose was larger by two orders of magnitude in humans (7.8 and 0.046 years per Gy for humans and mice, respectively, when exposed at ~35% of their lifetime). The difference was mostly explained by the two orders of magnitude difference in spontaneous somatic mutation rates between the species plus the species‐independent radiation‐induced mutation rate. Thus, the findings delineate the mechanism underlying the interspecies difference in radiation‐associated cancer mortality and may lead to the use of experimental evidence for risk prediction in humans.

Oncology↗
DOE OSTI2024

Denoising and Extension of Response Functions in the Time Domain

Response functions of quantum systems, such as electron Green’s functions, magnetic, or charge susceptibilities, describe the response of a system to an external perturbation. They are the central objects of interest in field theories and quantum computing and measured directly in experiment. Further, response functions are intrinsically causal. In equilibrium and steady-state systems, they correspond to a positive spectral function in the frequency domain. Since response functions define an inner product on a Hilbert space and thereby induce a positive definite function, the properties of this function can be used to reduce noise in measured data and, in equilibrium and steady state, to construct positive definite extensions for data known on finite time intervals, which are then guaranteed to correspond to positive spectra.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Verification Testing For Solid-Element Material Models 11-19 in DYNA3D/ParaDyn

This technical report documents the creation and implementation of verification tests for solid-element material models 11 through 19 available in DYNA3D/ParaDyn. The verification tests covered all aspects of each material model, except for the Weibull distribution functionality in material models 15 and 19. General test cases were created to verify the elastic and plastic behavior of the material models. Other additional tests were developed to examine the intricacies of each material model. Each test involved the use a kinematic load case and specification of material parameters necessary to activate corresponding features of the material model. The load cases prescribed the full time history of the kinematic motion for the solid elements, and these loads are independent of the material model or element formulation. When possible, closed form solutions were then derived for each verification test in a continuum setting. The DYNA3D simulations for each test were carried out over a unit time interval, t ϵ [0, 1], and the as implemented DYNA3D response was compared to the closed form solutions evaluated at discrete points in time. A relative error measure was determined for each test to justify the proper implementation of the material model. The relative errors comparing the DYNA3D solution to the analytical solution were, in general, on the order of machine precision except where noted. This signifies the proper implementation of solid-element material models 11-19. In the development of these verification tests, six bugs were found and fixed in the source code. Additionally, this work generated eighteen DYNA3D input decks and answer extraction scripts in the DYNA3D/ParaDyn Software Quality Assurance test suite, which are comprised of a total of 285 solid-element tests. Testing for each material model utilizes two input decks and answer extraction scripts, where one focuses on the linear elastic response and the other examines the inelastic and remaining functionalities of the material model. In total, this work added 285 individual verification test problems in the DYNA3D/ParaDyn test suite.

42 ENGINEERING↗
DOE OSTI2020

A pedagogical derivation of dynamical susceptibilities

Dynamical two-particle susceptibilites are important for a wide range of different experiments in condensed-matter physics and beyond. Nevertheless, most textbooks avoid describing how to derive such response functions, perhaps because they are viewed as too complex. In the literature, most derivations work with generalized susceptibilities, which are more general, but require an even higher layer of complexity. In this work, we show a more direct derivation in the context of model Hamiltonians which can be mapped directly onto an impurity model. We also present an alternative derivation for the irreducible vertex in the context of the Falicov-Kimball model.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Probing the edge states of Chern insulators using microwave impedance microscopy

Microwave impedance microscopy (MIM) has been utilized to directly visualize topological edge states in many quantum materials. While the microwave response for conventional metals and insulators can be accurately quantified using simple lumped-element circuits, whose applicability to more exotic quantum systems remain limited. In this work, we present a general theoretical framework of the MIM response of arbitrary quantum materials. Applying it to topological edge states in a Chern insulator predicts an enhanced MIM response at the crystal boundaries due to collective edge magnetoplasmon (EMP) excitations. The unique resonance frequency of these plasmonic modes allows one to disentangle the signatures of topological versus trivial edge states. To benchmark our analytical predictions, we experimentally probe the MIM response of quantum anomalous Hall edge states in a Cr-doped (Bi,Sb) 2 Te 3 topological insulator and perform numerical simulations using a classical formulation of the EMP modes based on this realistic tip-sample geometry, both of which yield results consistent with our theoretical picture. Here we also show how the technique of MIM can be used to quantitatively extract the topological invariant of a Chern insulator and shed light on the microscopic nature of dissipation along the crystal boundaries.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Constraining the Varied Response of Northern Hemisphere Winter Circulation Waviness to Climate Change

Abstract While a large latitudinal displacement of the westerly jet brings about disproportionate socioeconomic impacts over Northern Hemisphere midlatitude continents, it is not well understood as to whether the winter circulation will become wavier or less in response to climate change. Here, using observations and large ensembles of climate models, we show that changes in atmospheric waviness can be estimated from the optimal structures of the westerly jet for wavier circulation, which are obtained from an advection‐diffusion model. Thus, the changes in westerly jet structure in climate models under climate change provide a physical constraint on changes in atmospheric waviness, indicating that the North Atlantic wave activity will experience a robust decline in a warmer climate, while future North Pacific wave activity is obscured by model uncertainty rather than internal variability. These findings highlight the changes to jet stream structure as a constraint for regional circulation waviness in a changing climate.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Multiscale formulation of frictional contact mechanics at the pore scale

Direct numerical simulation (DNS) yields the highest fidelity predictions of mechanical deformation at the pore scale, but is prohibitively expensive for analyzing large or many samples. Discrete element methods (DEM) are an efficient alternative, but are limited to granular media and incapable of estimating or controlling prediction errors. We present a pore-level multiscale method (PLMM) that approximates DNS efficiently and with controllable accuracy. We focus on the linear elastic response of a consolidated geologic porous medium with arbitrary microstructure, heterogeneous mineralogy, containing cracks or defects. PLMM decomposes the solid phase into non-overlapping subdomains, on which local basis functions are constructed. The bases are then coupled with a global interface problem that accounts for slip or stick contact conditions between the subdomains. PLMM produces an initial, but accurate, approximation to DNS that can be iteratively improved. It is amenable to parallelism and allows for different mesh, models, and physics in each subdomain. An algebraic interpretation of PLMM as a preconditioner is also presented to allow non-intrusive implementation into existing solvers. Lastly, this work extends previous developments of PLMM in fluid dynamics to solid mechanics and enables future extensions towards modeling coupled flow and mechanics problems.

42 ENGINEERING↗
DOE OSTI2025

Determining spin-dependent light dark matter rates from neutron scattering

The scattering and absorption rates of light dark matter with electron spin-dependent interactions depend on the target’s spin response. We show how this response is encoded by the target’s dynamical magnetic susceptibility, which can be measured using neutron scattering. We directly use existing neutron scattering data to compute the dark matter scattering rate in a candidate target material, finding close agreement with the previous first-principles calculation at MeV dark matter masses. Complementary experiments and measurements can extend the reach of this technique to other dark matter models and masses, and identify promising target materials for future experiments.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Computational assessment of smooth and rough parameter dependence of statistics in chaotic dynamical systems

An assumption of smooth response to small parameter changes, of statistics or long-time averages of a chaotic system, is generally made in the field of sensitivity analysis, and the parametric derivatives of statistical quantities are critically used in science and engineering. In this paper, we propose a numerical procedure to assess the differentiability of statistics with respect to parameters in chaotic systems. Here, we numerically show that the existence of the derivative depends on the Lebesgue integrability of a certain density gradient function, which we define as the derivative of logarithmic SRB density along the unstable manifold. We develop a recursive formula for the density gradient that can be efficiently computed along trajectories, and demonstrate its use in determining the differentiability of statistics. Our numerical procedure is illustrated on low-dimensional chaotic systems whose statistics exhibit both smooth and rough regions in parameter space.

42 ENGINEERING↗
DOE OSTI2023

Vortical effects in chiral band structures

The chiral vortical effect is a chiral anomaly-induced transport phenomenon characterized by an axial current in a uniformly rotating chiral fluid. It is well understood for Weyl fermions in high energy physics, but its realization in condensed matter band structures, including those of Weyl semimetals, has been controversial. In this work, we develop the Kubo response theory for electrons in a general band structure subject to space- and time-dependent rotation or vorticity relative to the background lattice. For continuum Hamiltonians, we recover the chiral vortical effect in the static limit and the transport or uniform limit when the fluid, strictly, is not a Fermi liquid. In the transport limit of a Fermi liquid, we describe an effect that we dub the gyrotropic vortical effect. The latter is governed by Berry curvature of the occupied bands while the former contains an additional contribution from the magnetic moment of electrons on the Fermi surface. The two vortical effects can be understood as kinematic analogs of the well-known chiral and gyrotropic magnetic effects in chiral band structures. As a result, we address recent controversies in the field and conclude by describing device geometries that exploit Ohmic or Seebeck transport to drive the vortical effects.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Intrinsic Anomalous Hall Conductivity in a Nonuniform Electric Field

We study how the intrinsic anomalous Hall conductivity is modified in two-dimensional crystals with broken time-reversal symmetry due to weak inhomogeneity of the applied electric field. Focusing on a clean noninteracting two-band system without band crossings, we derive the general expression for the Hall conductivity at small finite wave vector q to order q 2 , which governs the Hall response to the second gradient of the electric field. Using the Kubo formula, we show that the answer can be expressed through the Berry curvature, Fubini-Study quantum metric, and the rank-3 symmetric tensor which is related to the quantum geometric connection and physically corresponds to the gauge-invariant part of the third cumulant of the position operator. We further compare our results with the predictions made within the semiclassical approach. By deriving the semiclassical equations of motion, we reproduce the result obtained from the Kubo formula in some limits. We also find, however, that the conventional semiclassical description in terms of the definite position and momentum of the electron is not fully consistent because of singular terms originating from the Heisenberg uncertainty principle. We thus present a clear example of a case when the semiclassical approach inherently suffers from the uncertainty principle, implying that it should be applied to systems in nonuniform fields with extra care.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Molecular dipole moment learning via rotationally equivariant derivative kernels in molecular-orbital-based machine learning

This study extends the accurate and transferable molecular-orbital-based machine learning (MOB-ML) approach to modeling the contribution of electron correlation to dipole moments at the cost of Hartree–Fock computations. A MOB pairwise decomposition of the correlation part of the dipole moment is applied, and these pair dipole moments could be further regressed as a universal function of MOs. The dipole MOB features consist of the energy MOB features and their responses to electric fields. An interpretable and rotationally equivariant derivative kernel for Gaussian process regression (GPR) is introduced to learn the dipole moment more efficiently. The proposed problem setup, feature design, and ML algorithm are shown to provide highly accurate models for both dipole moments and energies on water and 14 small molecules. To demonstrate the ability of MOB-ML to function as generalized density-matrix functionals for molecular dipole moments and energies of organic molecules, we further apply the proposed MOB-ML approach to train and test the molecules from the QM9 dataset. The application of local scalable GPR with Gaussian mixture model unsupervised clustering GPR scales up MOB-ML to a large-data regime while retaining the prediction accuracy. In addition, compared with the literature results, MOB-ML provides the best test mean absolute errors of 4.21 mD and 0.045 kcal/mol for dipole moment and energy models, respectively, when training on 110 000 QM9 molecules. The excellent transferability of the resulting QM9 models is also illustrated by the accurate predictions for four different series of peptides.

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

Computing nuclear response functions with time-dependent coupled-cluster theory

We compute nuclear response functions by solving the time-dependent 𝐴-body Schrödinger equation, recording the time-dependent transition moment and extracting spectral information via Fourier transforms. The solution of the time-dependent many-body problem accounts for correlations on top of the mean field by taking advantage of a time-dependent formulation of coupled-cluster theory. As a validation, we focus on electric dipole transitions in 4 He and 16 O and compare moments of the response function distribution to the results of an equivalent static framework, finding negligible discrepancies. We investigate how proton and neutron densities evolve in time, and we see the traditional picture of soft and giant dipole resonances as collective oscillations of protons and neutrons emerging from our calculations in 16 O and 24 O. Furthermore, this method also allows us to investigate the behavior of the nucleus in the presence of a strong electric field. In that regime, the behavior of the system becomes chaotic. Qualitatively, the spectral information obtained in this limit is in line with previous time-dependent mean-field results.

Ab initio calculations↗
DOE OSTI2023

Experimental and modeling studies of IPDI-based polyurea elastomers – The role of hard segment fraction

Segmented polyureas (PUa) are industrially important class of polymers widely used in coatings, sealant, and adhesive applications. Here, in this paper, we report synthesis, characterization, and modeling of Isophorone Diisocyanate-Diethyl-Toluene-Diamine-Polyether amine (IPDI-DETDA-PO PUa) with varied hard segment contents of 20, 30, and 40 weight percent. For each of the three materials, we study its structure and phase behavior using FTIR, DSC, and TEM, and clearly show the presence of microphase separation between the hard and soft nanodomains. We then measure the linear viscoelastic response of the PUa-s using DMA (frequency sweeps at multiple temperatures). The DMA data are shown to obey the time-temperature superposition. Finally, we develop a new micromechanical model describing the DMA results; the model describes a phase-separated PUa as two “Fractional-order Maxwell gels” branches, connected in parallel, with the first FMG branch representing the “percolated hard phase" and the second one modeling the “filled soft phase". In agreement with the earlier thermodynamic theories, the volume-fraction of the percolated hard phase is related to the hard segment weight-fraction (HSWF), defined as the combined mass of IPDI and DETDA normalized to the total mass of the polymer. The data and model are found to be in a good qualitative and quantitative agreement.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Differentiating densities on smooth manifolds

Lebesgue integration of derivatives of strongly-oscillatory functions is a recurring challenge in computational science and engineering. Integration by parts is an effective remedy for huge computational costs associated with Monte Carlo integration schemes. In case of Lebesgue integrals over a smooth manifold, however, integration by parts gives rise to a derivative of the density implied by charts describing the domain manifold. This paper focuses on the computation of that derivative, which we call the density gradient function, on general smooth manifolds. We analytically derive formulas for the density gradient and present examples of manifolds determined by popular differential equation-driven systems. Furthermore, we highlight the significance of the density gradient by demonstrating a numerical example of Monte Carlo integration involving oscillatory integrands.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

A parametric study of slow dynamic nonlinear elasticity with comparisons to models

Several phenomenological models that aspire to quantitative description of anomalous nonlinear mesoscopic elasticity are reviewed and compared with laboratory measurements. This class of nonlinearity, best known perhaps for slow dynamics and aging, is seen widely in imperfectly consolidated granular solids but is not well understood. Typical slow dynamic tests show that a modest conditioning oscillatory "pump" strain depresses material stiffness, which then recovers like the logarithm of time after conditioning ceases. Several phenomenological models based on physical arguments have been proposed that predict the material stiffness response to arbitrary pump strain histories during conditioning and recovery. Approximate closed form and numerical solutions to the models are presented that predict the quantitative influence of three key pump parameters: the pump's strain amplitude, the pump's strain rate, and the pump’s duration. Laboratory measurements on Berea sandstone, concrete and a confined single aluminum bead find that slow dynamic responses are linear in pump strain and independent of pump frequency. Measurements also show that, after pump-off, stiffness recovers over times far longer than the pump duration. These observations and others are compared to model predictions. One of the considered models, based on a picture of fast brittle damage and slow healing, successfully matches all these behaviors.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Laser-assisted embedding of all-glass optical fiber sensors into bulk ceramics for high-temperature applications

We develop a laser-assisted sensor embedding process to embed all-glass optical fiber sensors into bulk ceramics for high-temperature applications. A specially designed two-step microchannel was fabricated on an Al 2 O 3 substrate for sensor embedment using a picosecond (ps) laser. An optical fiber Intrinsic Fabry-Perot Interferometer (IFPI) sensor was embedded at the bottom of the microchannel and covered by Al 2 O 3 slurry which was subsequently sintered by a CO 2 laser. The sensor spectrum was in-situ monitored during the laser sintering process to ensure the survival of the sensor and optimize the laser sintering parameters. By testing in furnace through high temperature, the embedded optical fiber shows improved stability after CO 2 laser sealing, resulting in the linear temperature response of the embedded optical fiber IFPI sensor. To improve the embedded IFPI sensor for thermal strain measurement, a dummy fiber was co-embedded with the sensing fiber to improve the mechanical bonding between the sensing fiber and the ceramic substrate so that the thermal strain of the ceramic substrate can apply on the sensing fiber. The response sensitivity, measurement repeatability and high-temperature long-term stability of the embedded optical fiber IFPI sensor were evaluated in this work.

47 OTHER INSTRUMENTATION↗