Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “response theory”

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 91 records · Page 5

Theory of optical responses in clean multi-band superconductors

Electromagnetic responses in superconductors provide valuable information on the pairing symmetry as well as physical quantities such as the superfluid density. However, at the superconducting gap energy scale, optical excitations of the Bogoliugov quasiparticles are forbidden in conventional Bardeen-Cooper-Schrieffer superconductors when momentum is conserved. Accordingly, far-infrared optical responses have been understood in the framework of a dirty-limit theory by Mattis and Bardeen for over 60 years. Here we show, by investigating the selection rules imposed by particle-hole symmetry and unitary symmetries, that intrinsic momentum-conserving optical excitations can occur in clean multi-band superconductors when one of the following three conditions is satisfied: (i) inversion symmetry breaking, (ii) symmetry protection of the Bogoliubov Fermi surfaces, or (iii) simply finite spin-orbit coupling with unbroken time reversal and inversion symmetries. This result indicates that clean-limit optical responses are common beyond the straightforward case of broken inversion symmetry. We apply our theory to optical responses in FeSe, a clean multi-band superconductor with inversion symmetry and significant spin-orbit coupling. This result paves the way for studying clean-limit superconductors through optical measurements.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

A Theory for the Response of Tropical Moist Convection to Mechanical Orographic Forcing

Abstract Spatial patterns of tropical rainfall are strongly influenced by mountains. Although theories for precipitation induced by convectively stable upslope ascent exist for the midlatitudes, these do not represent the interaction of moist convection with orographic forcing. Here, we present a theory for convective precipitation produced by the mechanical interaction of a tropical ridge with a basic-state horizontal wind. Deviations from this basic state are represented as the sum of a “dry” perturbation, due to the stationary orographic gravity wave, and a “moist” perturbation that carries the convective response. The moist component dynamics are subject to the weak temperature gradient approximation; they are forced by the dry mode’s influence on lower-tropospheric moisture and temperature. Analytical solutions provide estimates of the precipitation distribution, including peak precipitation, upstream extent, and rain shadow extent. The theory can be used with several degrees of complexity depending on the technique used to compute the dry mode, which can be drawn from linear mountain wave theory or full numerical simulations. To evaluate the theory, we use a set of convection-permitting simulations with a flow-perpendicular ridge in a long channel. The theory makes a good prediction for the cross-slope precipitation profile, indicating that the organization of convective rain by orography can be quantitatively understood by considering the effect of stationary orographic gravity waves on a lower-tropospheric convective quasi-equilibrium state.

54 ENVIRONMENTAL SCIENCES↗

Testing unified theories for ozone response in C 4 species

Abstract There is tremendous interspecific variability in O 3 sensitivity among C 3 species, but variation among C 4 species has been less clearly documented. It is also unclear whether stomatal conductance and leaf structure such as leaf mass per area (LMA) determine the variation in sensitivity to O 3 across species. In this study, we investigated leaf morphological, chemical, and photosynthetic responses of 22 genotypes of four C 4 bioenergy species (switchgrass, sorghum, maize, and miscanthus) to elevated O 3 in side‐by‐side field experiments using free‐air O 3 concentration enrichment (FACE). The C 4 species varied largely in leaf morphology, physiology, and nutrient composition. Elevated O 3 did not alter leaf morphology, nutrient content, stomatal conductance, chlorophyll fluorescence, and respiration in most genotypes but reduced net CO 2 assimilation in maize and photosynthetic capacity in sorghum and maize. Species with lower LMA and higher stomatal conductance tended to show greater losses in photosynthetic rate and capacity in elevated O 3 compared with species with higher LMA and lower stomatal conductance. Stomatal conductance was the strongest determinant of leaf photosynthetic rate and capacity. The response of both area‐ and mass‐based leaf photosynthetic rate and capacity to elevated O 3 were not affected by LMA directly but negatively influenced by LMA indirectly through stomatal conductance. These results demonstrate that there is significant variation in O 3 sensitivity among C 4 species with maize and sorghum showing greater sensitivity of photosynthesis to O 3 than switchgrass and miscanthus. Interspecific variation in O 3 sensitivity was determined by direct effects of stomatal conductance and indirect effects of LMA. This is the first study to provide a test of unifying theories explaining variation in O 3 sensitivity in C 4 bioenergy grasses. These findings advance understanding of O 3 tolerance in C 4 grasses and could aid in optimal placement of diverse C 4 bioenergy feedstock across a polluted landscape.

54 ENVIRONMENTAL SCIENCES↗

An improved theory of the response of DIII-D H-mode discharges to static resonant magnetic perturbations and its implications for the suppression of edge localized modes

The plasma response to an externally generated, static, n = 2, resonant magnetic perturbation (RMP) in the pedestal region of DIII-D H-mode discharge #158115 is investigated using a toroidal generalization of the asymptotic matching model presented by Fitzpatrick [Phys. Plasmas 27, 042506 (2020)]. Just as in a recent paper [Q. M. Hu et al., Phys. Plasmas 26, 120702 (2019)], it is hypothesized that the density pump-out phenomenon is due to locked magnetic island chains induced at the bottom of the pedestal, whereas the suppression of edge localized modes is associated with a locked magnetic island chain induced at the top of the pedestal. Neutral penetration inside the last closed magnetic flux-surface is found to have a significant influence on locked magnetic island chain formation at the bottom of the pedestal. It is found that locked island formation at the top of the pedestal is only possible when q95 lies in certain narrow windows. Finally, it is inferred that, in order to successfully induce a locked island chain at the top of the pedestal, an external RMP field-coil system must generate a magnetic field that is simultaneously strongly amplified by the plasma (via the ideal kink response) and has a large resonant component at a rational surface that lies close to the top of the pedestal.

Fitzpatrick, R. (ORCID:0000000162379309)↗

Testing unified theories for ozone response in C4 species

Twenty-two genotypes of C4 species grown under ambient and elevated O3 concentration were studied at the SoyFACE (40°02’N, 88°14’W) in 2019. This dataset contains leaf morphology, photosynthesis and nutrient contents measured at three time points. The results of CO2 response curves are also included.

C4↗

Optimal stomatal theory predicts CO 2 responses of stomatal conductance in both gymnosperm and angiosperm trees

Optimal stomatal theory predicts that stomata operate to maximise photosynthesis (A net ) and minimise transpirational water loss to achieve optimal intrinsic water-use efficiency (iWUE). We tested whether this theory can predict stomatal responses to elevated atmospheric CO 2 (eCO 2 ), and whether it can capture differences in responsiveness among woody plant functional types (PFTs). We conducted a meta-analysis of tree studies of the effect of eCO 2 on iWUE and its components A net and stomatal conductance (g s ). We compared three PFTs, using the unified stomatal optimisation (USO) model to account for confounding effects of leaf–air vapour pressure difference (D). We expected smaller g s , but greater A net , responses to eCO 2 in gymnosperms compared with angiosperm PFTs. We found that iWUE increased in proportion to increasing eCO 2 in all PFTs, and that increases in A net had stronger effects than reductions in g s . The USO model correctly captured stomatal behaviour with eCO 2 across most datasets. The chief difference among PFTs was a lower stomatal slope parameter (g 1 ) for the gymnosperm, compared with angiosperm, species. Land surface models can use the USO model to describe stomatal behaviour under changing atmospheric CO 2 conditions.

54 ENVIRONMENTAL SCIENCES↗

Strong stretching theory for pH-responsive polyelectrolyte brushes in large salt concentrations

Here, we develop a theory for describing the thermodynamics, configuration, and electrostatics of strongly-stretched, pH-responsive polyelectrolyte (PE) brushes in the presence of large salt concentrations. The aim of the paper, therefore, is to study the properties of a PE brush in a salt concentration regime (namely, large concentrations of several molars) that has been hitherto unexplored theoretically in the context of PE brushes but can be routinely encountered in molecular scale simulations of the problem. The brushes are modelled using our recently developed augmented Strong Stretching Theory (SST), while the effect of the presence of the large salt concentration is accounted for by including the contributions of three different types of non-Poisson–Boltzmann (non-PB) effects to the free energy description of the PE brush induced electric double layer (EDL). These non-PB effects are ionic non-mean-field ion–ion correlations, solvent polarization, and the finite size effect of the ions and water dipoles. We study the individual influences of these different effects and show that the ion–ion correlations and solvent polarization effect reduce the brush height which consequently enhances the monomer density and leads to an electrostatic potential distribution of the brush induced EDL that has a larger magnitude at near-wall locations and becomes zero at shorter distances from the wall. The finite size effect, on the other hand, increases the brush height and therefore, weakens the monomer density and leads to a smaller near-wall magnitude of the EDL potential that becomes zero at larger distances from the wall. In time, we consider the impact of all the three non-PB effects simultaneously and show that the ion–ion correlations and solvent polarization effect dominate the size effects and dictate the overall brush configuration and the EDL electrostatics. We also point out that the influence of all the three non-PB effects becomes the largest for a larger salt concentration and a smaller bulk pH. Finally, we compare our theoretical predictions with those obtained from our recently developed all-atom MD simulation model and obtain an excellent match.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Interpreting angle-dependent magnetoresistance in layered materials: Application to cuprates

The evolution of the low-temperature electronic structure of the cuprate metals from the overdoped to the underdoped side has recently been addressed through angle-dependant magnetoresistance (ADMR) experiments in La 1.6-x Nd 0.4 Sr x CuO 4 . Here the results show a striking difference between hole dopings p = 0.24 and p = 0.21, which lie on either side of a putative quantum critical point at intermediate p. Motivated by this, we here study the theory of ADMR in correlated layered materials, paying special attention to the role of angle-dependent quasiparticle weights Z k . Such a Z k is expected to characterize a number of popular models of the cuprate materials, particularly when underdoped. Further, in the limit of weak interlayer hopping the quasiparticle weight will affect the c-axis transport measured in ADMR experiments. We show that proper inclusion of the quasiparticle weight does not support an interpretation of the data in terms of a (π,π) spin density wave ordered state, in agreement with the lack of direct evidence for such order. We show that a simple model of Fermi surface reconfiguring across a van Hove point captures many of the striking differences seen between p = 0.21 and p = 0.24. We comment on why such a model may be appropriate for interpreting the ADMR data, despite having a large Fermi surface at p = 0.21, seemingly in contradiction with other evidence for a small Fermi surface at that doping level.

36 MATERIALS SCIENCE↗

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↗

Can we use antipredator behavior theory to predict wildlife responses to high-speed vehicles?

Animals seem to rely on antipredator behavior to avoid vehicle collisions. There is an extensive body of antipredator behavior theory that have been used to predict the distance/time animals should escape from predators. These models have also been used to guide empirical research on escape behavior from vehicles. However, little is known as to whether antipredator behavior models are appropriate to apply to an approaching high-speed vehicle scenario. We addressed this gap by (a) providing an overview of the main hypotheses and predictions of different antipredator behavior models via a literature review, (b) exploring whether these models can generate quantitative predictions on escape distance when parameterized with empirical data from the literature, and (c) evaluating their sensitivity to vehicle approach speed using a simulation approach wherein we assessed model performance based on changes in effect size with variations in the slope of the flight initiation distance (FID) vs. approach speed relationship. The slope of the FID vs. approach speed relationship was then related back to three different behavioral rules animals may rely on to avoid approaching threats: the spatial, temporal, or delayed margin of safety. We used literature on birds for goals (b) and (c). Our review considered the following eight models: the economic escape model, Blumstein’s economic escape model, the optimal escape model, the perceptual limit hypothesis, the visual cue model, the flush early and avoid the rush (FEAR) hypothesis, the looming stimulus hypothesis, and the Bayesian model of escape behavior. We were able to generate quantitative predictions about escape distance with the last five models. However, we were only able to assess sensitivity to vehicle approach speed for the last three models. The FEAR hypothesis is most sensitive to high-speed vehicles when the species follows the spatial (FID remains constant as speed increases) and the temporal margin of safety (FID increases with an increase in speed) rules of escape. The looming stimulus effect hypothesis reached small to intermediate levels of sensitivity to high-speed vehicles when a species follows the delayed margin of safety (FID decreases with an increase in speed). The Bayesian optimal escape model reached intermediate levels of sensitivity to approach speed across all escape rules (spatial, temporal, delayed margins of safety) but only for larger (> 1 kg) species, but was not sensitive to speed for smaller species. Overall, no single antipredator behavior model could characterize all different types of escape responses relative to vehicle approach speed but some models showed some levels of sensitivity for certain rules of escape behavior. We derive some applied applications of our findings by suggesting the estimation of critical vehicle approach speeds for managing populations that are especially susceptible to road mortality. Overall, we recommend that new escape behavior models specifically tailored to high-speeds vehicles should be developed to better predict quantitatively the responses of animals to an increase in the frequency of cars, airplanes, drones, etc. they will face in the next decade.

54 ENVIRONMENTAL SCIENCES↗

Enhancing the efficiency of time-dependent density functional theory calculations of dynamic response properties

X-ray Thomson scattering (XRTS) constitutes an essential technique for diagnosing material properties under extreme conditions, such as high pressures and intense laser heating. Time-dependent density functional theory (TDDFT) is one of the most accurate available ab initio methods for modeling XRTS spectra, as well as a host of other dynamic material properties. However, strong thermal excitations, along with the need to account for variations in temperature and density as well as the finite size of the detector significantly increase the computational cost of TDDFT simulations compared to ambient conditions. In this work, we present a broadly applicable method for optimizing and enhancing the efficiency of TDDFT calculations. Our approach is based on a one-to-one mapping between the dynamic structure factor and the imaginary time density–density correlation function, which naturally emerges in Feynman’s path integral formulation of quantum many-body theory. Specifically, we combine rigorous convergence tests in the imaginary time domain with a constraints-based attenuation of narrow-band fluctuations to improve the efficiency of TDDFT modeling without the introduction of any significant bias. As a result, we can report a speed-up by up to an order of magnitude, thus substantially reducing the burden of computational cost required for XRTS analysis.

Moldabekov, Zhandos A. [Helmholtz-Zentrum Dresden-↗

A Gaussian field approach to the planar electric double layer structures in electrolyte solutions

Here, in this work, the planar, electric, double-layer structures of non-polarizable electrodes in electrolyte solutions are studied with Gaussian field theory. A response function with two Yukawa functions is used to capture the electrostatic response of the electrolyte solution, from which the modified response function in the planar symmetry is derived analytically. The modified response function is further used to evaluate the induced charge density and the electrostatic potential near an electrode. The Gaussian field theory, combined with a two-Yukawa response function, can reproduce the oscillatory decay behavior of the electric potentials in concentrated electrolyte solutions. When the exact sum rules for the bulk electrolyte solutions and the electric double layers are used as constraints to determine the parameters of the response function, the Gaussian field theory could at least partly capture the nonlinear response effect of the surface charge density. Comparison with results for a planar electrode with fixed surface charge densities from molecular simulations demonstrates the validity of Gaussian field theory.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Lectures on statistical mechanics

Presented here is a transcription of the lecture notes from Professor Allan N. Kaufman’s graduate statistical mechanics course Physics 212A and 212B at the University of California Berkeley from the 1972–1973 academic year. 212A addressed equilibrium statistical mechanics with topics: fundamentals (micro-canonical and sub-canonical ensembles, adiabatic law and action conservation, fluctuations, pressure, and virial theorem), classical fluids and other systems (equation of state, deviations from ideality, virial coefficients and van der Waals potential, canonical ensemble and partition function, quasistatic evolution, grand-canonical ensemble and partition function, chemical potential, simple model of a phase transition, quantum virial expansion, numerical simulation of equations of state, and phase transition), chemical equilibrium (systems with multiple species and chemical reactions, law of mass action, Saha equation, chemical equilibrium including ionization and excited states), and long-range interactions (including Coulomb, dipole, and gravitational interactions, Debye–Hückel theory, and shielding). 212B addressed nonequilibrium statistical mechanics with topics: fundamentals (definitions: realizations, moments, characteristic function, and discrete variables), Brownian motion (Langevin equation, fluctuation–dissipation theorem, spatial diffusion, Boltzmann’s H-theorem), Liouville and Klimontovich equations, Landau equation (derivation, elaboration, and H-theorem, and irreversibility), Markov processes and Fokker–Planck equation (derivations of the Fokker–Planck equation and a master equation), linear response and transport theory (linear Boltzmann equation, linear response theory of Kubo and Mori, relation of entropy production to electrical conductivity, transport relations and coefficients, normal mode solutions of the transport equations, sketch of a generalized Langevin equation method for transport theory), and an introduction to nonequilibrium quantum statistical mechanics.

plasma dynamics↗

On-the-fly response function generation method for composite coarse mesh

The hybrid stochastic deterministic transport code COMET, based on the incident response expansion theory, is used to model reactor cores with high fidelity and formidable computational speed. COMET models a reactor core using a library of incident flux response expansion coefficients that are pre computed for all the unique lattice cells (e.g., fuel assemblies, reflector blocks, etc.) in the core. In order to further improve its computational efficiency in pre-calculating the response library a new response function generation method is developed to compute the response functions for the composite coarse meshes made of a smaller set of unique lattices on the fly within the COMET's deterministic transport core sweep. The efficiency is achieved by eliminating a number of unique lattices that can be made up from the reduced set of unique meshes on the fly. The numerical process consists of the following steps. First, the boundary condition on composite coarse mesh boundaries is projected onto the expansion basis to compute the incident flux moments on external surfaces of all the basic (reduced set of unique) coarse meshes. Secondly, the deterministic sweeping solver in COMET is used to converge on the outgoing/incoming flux expansion moments crossing interfaces between the basic coarse meshes. Thirdly, the response functions for the composite coarse meshes are constructed as a superposition on the fly. The new response function generation method was tested on 88 composite coarse meshes consisting of CANDU fuel bundles and moderator blocks. It was found that response functions generated by the new method agree very well with those generated by direct Monte Carlo calculations. The average and maximum relative differences in the surface-to-surface response coefficients computed by the two methods are 0.10% and 0.20%, respectively. Similarly, the average and maximum relative differences in the response fission densities are 0.13% and 0.43%, respectively. These discrepancies are within one standard deviation of the stochastic uncertainties. The new method is five times faster than the original direct Monte Carlo method. The size of the response function library for the new method is five times smaller than that for the original method, leading to significantly less requirement for the computer hard drive space and memory. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

First-principles residual resistivity using a locally self-consistent multiple scattering method

The locally self-consistent multiple scattering (LSMS) method can perform efficient first-principles calculations of systems with a large number of atoms. Here, in this paper, we combine the Kubo-Greenwood equation with LSMS, enabling us to calculate the first-principles residual resistivity of large systems. This has been implemented in the open-source code LSMS . We apply this method to selected pure elements and binary random alloys. The results compare well with experiment, and with values obtained from a first-principles effective medium technique (the Korringa-Kohn-Rostoker coherent potential approximation). We discuss future applications of this method to complex systems where other methods are not applicable.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗