Search NASA⌕ Search

SEARCH · Search NASA

Results for “displacement functions”

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

Shear Softening in a Metallic Glass: First-Principles Local-Stress Analysis

Metallic glasses deform elastically under stress. However, the atomic-level origin of elastic properties of metallic glasses remain unclear. In this Letter using ab initio molecular dynamics simulations of the Cu 5 0Zr 50 metallic glass under shear strain, we show that the heterogeneous stress relaxation results in the increased charge transfer from Zr to Cu atoms, enhancing the softening of the shear modulus. Changes in compositional short-range order and atomic position shifts due to the nonaffine deformation are discussed. It is shown that the Zr subsystem exhibits a stiff behavior, whereas the displacements of Cu atoms from their initial positions, induced by the strain, provide the stress drop and softening.

36 MATERIALS SCIENCE↗

Out-of-plane polarization and topological magnetic vortices in multiferroic CrPSe 3

Two-dimensional multiferroic materials are ideal systems for exploring new coupling mechanisms between different ferroic orders and producing novel quantum phenomena with potential applications. Here we employ first-principles density functional theory calculations to discover intrinsic ferroelectric and antiferroelectric phases of Cr⁢P⁢Se 3 , which show ferromagnetic order and compete with the centrosymmetric phase with an antiferromagnetic order. Our analysis shows that the electrical dipoles of such type-I multiferroic phases come from the out-of-plane displacements of phosphorus ions due to the stereochemically active lone pairs. The coupling between polar and magnetic orders creates the opportunity for tuning the magnetic ground state by switching from the centrosymmetric to the ferroelectric phase using an out-of-plane electric field. In ferroelectric and antiferroelectric phases, the combination of easy-plane anisotropy and Dzyaloshinskii-Moriya interactions indicates that they can host topological magnetic vortices such as meron pairs.

36 MATERIALS SCIENCE↗

Investigating parameter trainability in the SNAP-Displacement protocol of a qudit system.

Our main focus is to investigate the sensitivity of training any of the SNAP parameters in the SNAP-Displacement protocol. We analyze conditions that could potentially lead to the Barren Plateau problem in a qudit system and draw comparisons with multi-qubit systems. The parameterized ansatz we consider consists of SNAP -Displacement blocks. We utilize techniques similar to those in [8] and [2] along with the concept of $t-$design. Through this analysis, we make the following key observations: (a) The trainability of a SNAP-parameter does not exhibit a preference for any particular direction within our cost function landscape, (b) using Haar measure properties, we establish new lemmas concerning the expectation of certain polynomial functions, and (c) utilizing these new lemmas, we identify a general condition that indicates an expected trainability advantage in a qudit system when compared to multi-qubit systems.

43 PARTICLE ACCELERATORS↗

Electromechanical memcapacitor model offering biologically plausible spiking.

In this article, we introduce a new nanoscale electromechanical device - a leaky memcapacitor - and show that it may be useful for the hardware implementation of spiking neurons. The leaky memcapacitor is a movableplate capacitor that becomes quite conductive when the plates come close to each other. The equivalent circuit of the leaky memcapacitor involves a memcapacitive and memristive system connected in parallel. In the leaky memcapacitor, resistance and capacitance depend on the same internal state variable, which is the displacement of the movable plate. We have performed a comprehensive analysis showing that several types of spiking observed in biological neurons can be implemented with the leaky memcapacitor. Significant attention is paid to the dynamic properties of the model. As in leaky memcapacitors the capacitive, leaking resistive, and reset functionalities are implemented naturally within the same device structure, their use will simplify the creation of spiking neural networks.

Zhang, Zixi↗

Dynamic Manipulation in Piezoresponse Force Microscopy: Creating Nonequilibrium Phases with Large Electromechanical Response

We discuss how domain walls and topological defects in ferroelectric materials have emerged as a powerful tool for functional electronic devices including memory and logic. Similarly, wall interactions and dynamics underpin a broad range of mesoscale phenomena ranging from giant electromechanical responses to memory effects. Exploring the functionalities of individual domain walls, their interactions, and controlled modifications of the domain structures is crucial for applications and fundamental physical studies. However, the dynamic nature of these features severely limits studies of their local physics since application of local biases or pressures in piezoresponse force microscopy induce wall displacement as a primary response. Here, we introduce an approach for the control and modification of domain structures based on automated experimentation, whereby real-space image-based feedback is used to control the tip bias during ferroelectric switching, allowing for modification routes conditioned on domain states under the tip. This automated experiment approach is demonstrated for the exploration of domain wall dynamics and creation of metastable phases with large electromechanical response.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Ring Pull Strain Analysis

Ring Pull Strain Analysis (RPSA) is a framework for analysis of Digital Image Correlation (DIC) data of gaugeless ring pull testing. It is meant for analysis of rings cut from tube cross sections and pulled in a displacement-controlled tensile test. This code takes load frame data, images take during testing, and csv files from the DIC analysis of those images and allows the user to analyze and plot data from the test. Some important features of the software include getting the strain values as a function of radial distance and angle, plotting of the strain around the ring, and stress-strain analysis of the load frame data.

Beck, Peter↗

Ring Pull Strain Analysis

Ring Pull Strain Analysis (RPSA) is a framework for analysis of Digital Image Correlation (DIC) data of gaugeless ring pull testing. It is meant for analysis of rings cut from tube cross sections and pulled in a displacement-controlled tensile test. This code takes load frame data, images take during testing, and csv files from the DIC analysis of those images and allows the user to analyze and plot data from the test. Some important features of the software include getting the strain values as a function of radial distance and angle, plotting of the strain around the ring, and stress-strain analysis of the load frame data. Additional capability is included with Ring Pull Coating Analysis (RPCA, also part of the software). Additional images are taken of samples with an outer coating and are used to correlate failure of the coating with strain values measured with DIC.

Beck, Peter↗

Active learning with multifidelity modeling for efficient rare event simulation

Here, while multifidelity modeling provides a cost-effective way to conduct uncertainty quantification with computationally expensive models, much greater efficiency can be achieved by adaptively deciding the number of required high-fidelity (HF) simulations, depending on the type and complexity of the problem and the desired accuracy in the results. We propose a framework for active learning with multifidelity modeling emphasizing the efficient estimation of rare events. Our framework works by fusing a low-fidelity (LF) prediction with an HF-inferred correction, filtering the corrected LF prediction to decide whether to call the high-fidelity model, and for enhanced subsequent accuracy, adapting the correction for the LF prediction after every HF model call. The framework does not make any assumptions as to the LF model type or its correlations with the HF model. In addition, for improved robustness when estimating smaller failure probabilities, we propose using dynamic active learning functions that decide when to call the HF model. We demonstrate our framework using several academic case studies (including some high-dimensional problems) and two finite element model case studies: estimating Navier-Stokes velocities using the Stokes approximation and estimating stresses in a transversely isotropic model subjected to displacements via a coarsely meshed isotropic model. Across these case studies, not only did the proposed framework estimate the failure probabilities accurately, but compared with either Monte Carlo or a standard variance reduction method, it also required only a small fraction of the calls to the HF model.

42 ENGINEERING↗

Switching of Hybrid Improper Ferroelectricity in Oxide Double Perovskites

In ABO 3 -type perovskite oxides with Pnma symmetry, rotation (Q R+ , a 0 a 0 c + ) and tilt (Q T , a – a – c 0 ) of BO 6 octahedra are the two primary order parameters. These order parameters establish an inherent trilinear coupling with anti-ferroelectric A-site displacement (Q AFE ) to form the low-symmetry phase. The symmetry is further lowered in double perovskite oxides (DPOs) due to A/A' cation ordering. It in turn makes these systems polar via hybrid improper ferroelectric mechanism, primarily driven by Q R+ and Q T . Naturally, it has been believed that functionalities such as polarization can also be switched by tuning these primary order parameters. However, mystery around finding switching mechanism still remains. Our study based on density functional theory calculations combined with finite-temperature molecular dynamics simulations shows that the polarization switching is a two-step process, driven by out-of-phase rotation (Q R– , a 0 a 0 c – when Q T = 0 or, a – a – b – when Q T ≠ 0). A series of polar DPOs such as KLnFeOsO 6 [Ln = Sm, Gd, Dy, Tm (lanthanides) and Y (rare earth)], all belonging to P2 1 symmetry, are considered in this investigation. The polarization switching P ($\overrightarrow{P}$) occurs at a very high temperature of ~1150 K through a phase transition, from a polar (P2 1 ) phase with $\overrightarrow{P}$(+) to $\overrightarrow{P}$(-) via a non-polar P4/n phase. The switching itself is metastable in nature. The switching (both polarization and spin state) is only observed for a very short period of time (~23 ps) that poses limitation on using such a mechanism in memory device realization. We demonstrate a concurrent heating–cooling procedure to overcome such shortcoming. In conclusion, simulations conducted at 600 K further imply that long lasting switching can be achieved, at least for 1.2 ns for 600 K, and ideally for an infinite time, if the material is heated just above the T c followed by rapid cooling to a temperature below T c .

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

High-throughput computation of electric polarization in solids via Berry flux diagonalization

Electric polarization in the absence of an externally applied electric field is a key property of polar materials, but the standard interpolation-based ab initio approach to compute polarization differences within the modern theory of polarization presents challenges for automated high-throughput calculations. Berry flux diagonalization [J. Bonini et al., Phys. Rev. B 102, 045141 (2020)] has been proposed as an efficient and reliable alternative, though it has yet to be widely deployed. Here, we assess Berry flux diagonalization using ab initio calculations of a large set of materials, introducing and validating heuristics that ensure branch alignment with a minimal number of intermediate interpolated structures. Our automated implementation of Berry flux diagonalization succeeds in cases where prior interpolation-based workflows fail due to band-gap closures or branch ambiguities. Benchmarking with ab initio calculations of 176 candidate ferroelectrics, we demonstrate the efficacy of the approach on a broad range of insulating materials and obtain accurate effective polarization values with fewer interpolated structures than prior automated interpolation-based workflows. Our real-space heuristics that can predict gauge stability a priori from ionic displacements enable a general automated framework for reliable polarization calculations and efficient high-throughput screening of chemically and structurally diverse polar insulators. These results establish Berry flux diagonalization as a robust and efficient method to compute the effective polarization of solids and to accelerate the data-driven discovery of functional polar materials.

Poteshman, Abigail N. [University of Chicago, IL (↗

A Lagrangian perturbation theory in the presence of massive neutrinos

We develop a Lagrangian Perturbation Theory (LPT) framework to study the clustering of cold dark matter (CDM) in cosmologies with massive neutrinos. We follow the trajectories of CDM particles with Lagrangian displacements fields up to third order in perturbation theory. Once the neutrinos become non-relativistic, their density fluctuations are modeled as being proportional to the CDM density fluctuations, with a scale-dependent proportionality factor. This yields a gravitational back-reaction that introduces additional scales to the linear growth function, which is accounted for in the higher order LPT kernels. Through non-linear mappings from Eulerian to Lagrangian frames, we ensure that our theory has a well behaved large scale behavior free of unwanted UV divergences, which are common when neutrino and CDM densities are not treated on an equal footing, and in resummation schemes that manifestly break Galilean invariance. We use our theory to construct correlation functions for both the underlying matter field, as well as for biased tracers using ConvolutionLPT. Redshift-space distortions effects are modeled using the Gaussian Streaming Model. When comparing our analytical results to simulated data from the Quijote 1 simulation suite, we find good accuracy down to r = 20 Mpc h -1 at redshift z = 0.5, for the real space and redshift space monopole particle correlation functions with no free parameters. The same accuracy is reached for the redshift space quadrupole if we additionally consider an effective field theory parameter that shifts the pairwise velocity dispersion. For modeling the correlation functions of tracers we adopt a simple Lagrangian biasing scheme with only density and curvature operators, which we find sufficient to reach down to r = 20 Mpc h -1 when comparing to simulated halos.

79 ASTRONOMY AND ASTROPHYSICS↗

A mixed, unified forward/inverse framework for earthquake problems: fault implementation and coseismic slip estimate

SUMMARY We introduce a new finite-element (FE) based computational framework to solve forward and inverse elastic deformation problems for earthquake faulting via the adjoint method. Based on two advanced computational libraries, FEniCS and hIPPYlib for the forward and inverse problems, respectively, this framework is flexible, transparent and easily extensible. We represent a fault discontinuity through a mixed FE elasticity formulation, which approximates the stress with higher order accuracy and exposes the prescribed slip explicitly in the variational form without using conventional split node and decomposition discrete approaches. This also allows the first order optimality condition, that is the vanishing of the gradient, to be expressed in continuous form, which leads to consistent discretizations of all field variables, including the slip. We show comparisons with the standard, pure displacement formulation and a model containing an in-plane mode II crack, whose slip is prescribed via the split node technique. We demonstrate the potential of this new computational framework by performing a linear coseismic slip inversion through adjoint-based optimization methods, without requiring computation of elastic Green’s functions. Specifically, we consider a penalized least squares formulation, which in a Bayesian setting—under the assumption of Gaussian noise and prior—reflects the negative log of the posterior distribution. The comparison of the inversion results with a standard, linear inverse theory approach based on Okada’s solutions shows analogous results. Preliminary uncertainties are estimated via eigenvalue analysis of the Hessian of the penalized least squares objective function. Our implementation is fully open-source and Jupyter notebooks to reproduce our results are provided. The extension to a fully Bayesian framework for detailed uncertainty quantification and non-linear inversions, including for heterogeneous media earthquake problems, will be analysed in a forthcoming paper.

58 GEOSCIENCES↗

Surfactant inhibition mechanisms of carbonate mineral dissolution in shale

Surfactants are common additives to hydraulic fracturing and enhanced oil recovery (EOR) fluids, and are under consideration for amendment to supercritical carbon dioxide for geological carbon sequestration (GCS). The effect of a common anionic surfactant, internal olefin sulfonate (IOS), on mineral dissolution from shale into brine was evaluated. When added to brine at concentrations exceeding the critical micelle concentration (94 mg/L), IOS inhibited carbonate mineral dissolution in an Eagle Ford shale, as well as dissolution of optical quality calcite (the dominant carbonate in the shale). Laser profilometry images provide spatial resolution across > 3 orders of magnitude, and indicate that IOS addition to brine both enhances the formation of new etch pits in calcite, and impedes their further growth. Time-of-flight secondary ion mass spectrometry surface profiles show for the first time that IOS preferentially adsorbs at calcite pit edges versus flat calcite surfaces (i.e., terraces). Surface pressure calculations, sulfur K-edge near edge X-ray absorption fine structure (NEXAFS) spectroscopy results, and density functional theory (DFT) calculations support this observation; the DFT results indicate that the sulfonate head group of the IOS molecule binds strongly to the calcite step site as compared to the terrace site. The S K-edge NEXAFS results indicate that IOS adsorbed more to etched calcite surfaces compared to smooth calcite surfaces. Overall, the results indicate that weak adsorption on flat calcite surfaces (i.e., terraces) disrupts water structure and enhances mass transfer of dissolution, while strong adsorption on calcite pit edges displaces adsorbed water and inhibits further etch pit growth. This work provides the first direct evidence of preferential adsorption of IOS to etched calcite surfaces and links it to macroscopic dissolution kinetics. Finally, this work has implications for surfactant-containing fluids used in hydraulic fracturing, EOR and potentially GCS for subsurface injection into carbonate rich reservoirs.

36 MATERIALS SCIENCE↗

Ab initio self-consistent many-body theory of polarons at all couplings

We present a theoretical framework to describe polarons from first principles within a many-body Green’s function formalism. Starting from a general electron-phonon Hamiltonian, we derive a self-consistent Dyson equation in which the phonon-mediated self-energy is composed by two distinct terms. One term is the Fan-Migdal self-energy and describes dynamic electron-phonon processes, the other term is a contribution to the self-energy originating from the static displacements of the atomic nuclei in the polaronic ground state. The lowest-order approximation to the present theory yields the standard many-body perturbation theory approach to electron-phonon interactions in the limit of large polarons, and the ab initio polaron equations introduced in the limit of small polarons. Here, a practical recipe to implement the present unifying formalism in first-principles calculations is outlined. We apply our method to the Fröhlich model, and obtain remarkably accurate polaron energies at all couplings, in line with Feynman’s polaron theory and diagrammatic Monte Carlo calculations. We also recover the well-known results of Fröhlich and Pekar at weak and strong coupling, respectively. The present approach enables predictive many-body calculations of polarons in real materials at all couplings.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Ultrafast heat transfer in single palladium nanocrystals seen with an X-ray free-electron laser

We report transient highly strained structural states in individual palladium (Pd) nanocrystals, electronically heated using an optical laser, which precede their uniform thermal expansion. Using an X-ray free-electron laser probe, the evolution of individual 111 Bragg peaks is measured as a function of delay time at various laser fluences. Above a laser fluence threshold at a sufficient pump-probe delay, the Bragg peak splits into multiple peaks, indicating heterogeneous strain, before returning to a single peak, corresponding to even heat distribution throughout the lattice expanded crystal. Our findings are supported by a lattice displacement and strain model of a single nanocrystal at different delay times, which agrees with the experimental data. Our observations have implications for understanding femtosecond laser interactions with metals and the potential photo-catalytic performance of Pd.

36 MATERIALS SCIENCE↗

Activated penetrant dynamics in glass forming liquids: size effects, decoupling, slaving, collective elasticity and correlation with matrix compressibility

In this work, we employ the microscopic self-consistent cooperative hopping theory of penetrant activated dynamics in glass forming viscous liquids and colloidal suspensions to address new questions over a wide range of high matrix packing fractions and penetrant-to-matrix particle size ratios. The focus is on the mean activated relaxation time of smaller tracers in a hard sphere fluid of larger particle matrices. This quantity also determines the penetrant diffusion constant and connects directly with the structural relaxation time probed in an incoherent dynamic structure factor measurement. The timescale of the non-activated fast dissipative process is also studied and is predicted to follow power laws with the contact value of the penetrant–matrix pair correlation function and the penetrant–matrix size ratio. For long time penetrant relaxation, in the relatively lower packing fraction metastable regime the local cage barriers are dominant and matrix collective elasticity effects unimportant. As packing fraction and/or penetrant size grows, much higher barriers emerge and the collective elasticity associated with the correlated matrix dynamic displacement that facilitates penetrant hopping becomes important. This results in a non-monotonic variation with packing fraction of the degree of decoupling between the matrix and penetrant alpha relaxation times. The conditions required for penetrant hopping to become slaved to the matrix alpha process are determined, which depend mainly on the penetrant to matrix particle size ratio. By analyzing the absolute and relative importance of the cage and elastic barriers we establish a mechanistic understanding of the origin of the predicted exponential growth of the penetrant hopping time with size ratio predicted at very high packing fractions. A dynamics-thermodynamics power law connection between the penetrant activation barrier and the matrix dimensionless compressibility is established as a prediction of theory, with different scaling exponents depending on whether matrix collective elasticity effects are important. Quantitative comparisons with simulations of the penetrant relaxation time, diffusion constant, and transient localization length of tracers in dense colloidal suspensions and cold viscous liquids reveal good agreements. Multiple new predictions are made that are testable via future experiments and simulations. Extension of the theoretical approach to more complex systems of high experimental interest (nonspherical molecules, semiflexible polymers, crosslinked networks) interacting via variable hard or soft repulsions and/or short range attractions is possible, including under external deformation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Signatures of linearized gravity in atom interferometers: A simplified computational framework

We develop a general framework for calculating the leading-order, general relativistic contributions to the gravitational phase shift in single-photon atom interferometers within the context of linearized gravity. We show that the atom gradiometer observable, which only depends on the atom interferometer propagation phase, can be written in terms of three distinct contributions: the Doppler phase shift, which accounts for the tidal displacement of atoms along the baseline, the Shapiro phase shift, which accounts for the delay in the arrival time of photons at atom-light interaction points, and the Einstein phase shift, which accounts for the gravitational redshift measured by the atoms. For specific atom gradiometer configurations, we derive the signal and response functions for two physically motivated scenarios: (i) transient gravitational waves in the transverse-traceless gauge and, for the first time, in the proper detector frame, and (ii) transient massive objects sourcing weak and slow-varying Newtonian potentials. Here, we find that the Doppler contribution of realistic Newtonian noise sources (e.g., a freight truck or a piece of space debris) at proposed atom gradiometer experiments, such as AION, MAGIS, and AEDGE, can exceed the shot noise level and thus affect physics searches if not properly subtracted. Note: 44 pages including appendices, 4 figures; v2: different formatting, fixed typos, 28 pages including appendices, 5 figures, agrees with published version

Atom interferometry↗

Radiation-Enhanced Anion Transport in Hematite

The influence of radiation-induced (1 MeV energy H+ to ~0.1 dpa at 450°C), non-equilibrium point defect populations on mass transport is studied with an integrated campaign of experimental and theoretical methods. Using epitaxial thin films of hematite with embedded 18O tracer layers and nanoscale atom probe tomography measurements, it is shown that anion self-diffusion is enhanced by at least 2 orders of magnitude under irradiation compared to thermal diffusion alone. Complementary scanning transmission electron microscopy of vacuum annealed specimens, reveals associated microstructural changes in the oxide films, including local phase transformation to Fe3O4 and the development of nanoscale voids from vacancy coalescence. Point defect formation and migration energies were computed from density functional theory and applied within the context of chemical rate theory to analyze contributions from both interstitial and vacancy mechanisms to self-diffusion in thermal and irradiation conditions. Comparisons are made between calculated, literature and newly measured self-diffusion values, revealing good agreement on the magnitude of radiation-enhanced anion diffusion. Further, the model suggests a transition from vacancy to interstitialcy mechanisms at low temperatures and high oxygen activity, consistent with the varied activation energies reported from prior studies.

iron oxides, isotopic tracers, atom probe tomograp↗