Search NASA⌕ Search

SEARCH · Search NASA

Results for “weighted relaxation”

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 19 records

Weighted relaxation for multigrid reduction in time

Current trends in computer architectures now mean that faster computation speed must come primarily from increased concurrency, not faster clock speeds, which are stagnating. Thus, this situation creates bottlenecks for serial algorithms, including the well-known bottleneck for sequential time-integration, where each individual time-value (i.e., time-step) is computed sequentially. One approach to alleviate this and achieve parallelism in time is with multigrid. Here, in this work, we consider multigrid-reduction-in-time (MGRIT), a multilevel method applied to the time dimension that computes multiple time-steps in parallel. Like all multigrid methods, MGRIT relies on the complementary relationship between relaxation on a fine-grid and a correction from the coarse grid to solve the problem. All current MGRIT implementations are based on unweighted-Jacobi relaxation; here we introduce the concept of weighted relaxation to MGRIT. We derive new convergence bounds for weighted relaxation, and use this analysis to guide the selection of relaxation weights. Numerical results then demonstrate that by choosing appropriate non-unitary relaxation weights, one can achieve faster convergence rates and lower iteration counts for MGRIT when compared with unweighted relaxation. In most cases, weighted relaxation yields a 10%–20% saving in iterations, which is significant when using large high-performance computers. For A-stable integration schemes, results also illustrate that under-relaxation can restore convergence in some cases where unweighted relaxation is not convergent.

97 MATHEMATICS AND COMPUTING↗

Aging and Lifetimes Milestone 7777, G.C. #4: Assess the ability of MRI to identify chemical and mechanical heterogeneity in polymers at various length scales

Polysiloxane foams utilized for stockpile applications experience changes over time induced by exposure to thermal, radiative, and mechanical stresses. The interplay of these stresses can lead to chemical and/or physical changes in the polymer, such as compression set and material degradation over time. Understanding the mechanisms of material aging and identifying material degradation prior to failure are critical in maintaining the integrity of the nuclear stockpile. Our prior work using solid-state NMR spectroscopy has shown that we can observe signatures of thermal, radiative, and mechanical aging of silicone elastomers using differences in NMR parameters such as relaxation (T1 and T2). However, most of our measurements have been on bulk samples and have only provided an average aging behavior of each sample. Here, we have used magnetic resonance imaging (MRI) to provide relaxation-weighted images that are a metric of structural heterogeneity (cross-link density) in polymers across multiple length scales. In addition, the larger sample size that can be accommodated in animal MRI scanners has the potential to make the technique adaptable to detecting heterogeneities whole parts.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

LLNL FY22 Aging and Lifetimes Exit Criteria: Milestone 8119, GC#01

Polysiloxane foams utilized for stockpile applications experience changes over time induced by exposure to thermal, radiative, and mechanical stresses. The interplay of these stresses can lead to chemical and/or physical changes in the polymer, such as compression set and material degradation over time. Understanding the mechanisms of material aging and identifying material degradation prior to failure are critical in maintaining the integrity of the nuclear stockpile. Our prior work using solid-state NMR spectroscopy has shown that we can observe signatures of thermal, radiative, and mechanical aging of silicone elastomers using differences in NMR parameters such as relaxation (T1 and T2). However, most of our measurements have been on bulk samples and have only provided an average aging behavior of each sample. Here, we have used magnetic resonance imaging (MRI) to provide relaxation-weighted images that are a metric of structural heterogeneity (crosslink density) in polymers across multiple length scales. In addition, the larger sample size that can be accommodated in animal MRI scanners has the potential to make the technique adaptable to detecting heterogeneities in whole parts.

36 MATERIALS SCIENCE↗

A Near-infrared Chemical Inventory of the Atmosphere of 55 Cancri e

We present high-resolution near-infrared spectra taken during eight transits of 55 Cancri e, a nearby low-density super-Earth with a short orbital period (<18 hr). While this exoplanet’s bulk density indicates a possible atmosphere, one has not been detected definitively. Our analysis relies on the Doppler cross-correlation technique, which takes advantage of the high spectral resolution and broad wavelength coverage of our data, to search for the thousands of absorption features from hydrogen-, carbon-, and nitrogen-rich molecular species in the planetary atmosphere. Although we are unable to detect an atmosphere around 55 Cancri e, we do place strong constraints on the levels of HCN, NH{sub 3}, and C{sub 2}H{sub 2} that may be present. In particular, at a mean molecular weight of 5 amu, we can rule out the presence of HCN in the atmosphere down to a volume mixing ratio (VMR) of 0.02%, NH{sub 3} down to a VMR of 0.08%, and C{sub 2}H{sub 2} down to a VMR of 1.0%. If the mean molecular weight is relaxed to 2 amu, we can rule out the presence of HCN, NH{sub 3}, and C{sub 2}H{sub 2} down to VMRs of 0.001%, 0.0025%, and 0.08%, respectively. Our results reduce the parameter space of possible atmospheres consistent with the analysis of Hubble Space Telescope/WFC3 observations by Tsiaras et al. and indicate that if 55 Cancri e harbors an atmosphere, it must have a high mean molecular weight or clouds.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

An analysis of the spatio-temporal resolution of the immersed boundary method with direct forcing

The immersed boundary method (IBM) with direct forcing is very popular in the simulation of rigid particulate flows. In the IBM, an interaction force is introduced at the interface between fluid and particle in order to approximate the no-slip boundary condition. The interaction force is calculated through dividing the velocity difference (or error) between fluid and particle at the interface by the time step. Here, a dynamic equation for the velocity difference is derived. Additionally, analyses on the dynamic equation provide a few new findings: (i) The interaction force is the solution of a least-squares error problem, with the direct implication that the Lagrangian marker distribution has no effect on the large scale flow structure once the distribution of Lagrangian markers become saturated along the interface (i.e., each marker remains properly correlated with all its neighbors); (ii) The Lagrangian volume-weight is a relaxation factor to control how fast the velocity error decays to the ideal value of zero; (iii) The optimal choice of the Lagrangian volume-weight is the largest value permissible by a stability condition. A comprehensive convergence analysis with regard to the spatial and temporal resolution is presented for the velocity error and also for the shear-stress and surface pressure. In three simple canonical problems, it is analytically and numerically shown that the IBM results converge to the theoretical solutions obtained with precise imposition of no-slip and no-penetration boundary conditions. It is observed that it is not necessary to match the Lagrangian marker volume-weight to that of the local Eulerian cell volume and in fact this matching leads to lower than optimal computational efficiency. However, it is found that extremely high Eulerian grid resolution and small time step have to be used to obtain high precision simulation results. Especially, the time step should be inversely proportional to the particle Reynolds number for low Reynolds number flows. For high frequency oscillation problems, the grid size needs to be reduced by a factor of the square root of the frequency, and the time step to be reduced by a factor of the frequency. The theoretical findings here can be used to alleviate the technical difficulties in simulating non-spherical particles by not requiring the Lagrangian marker distribution to match the Eulerian grids and also in the implementation of IBM on non-uniform Eulerian grids. The present work also provides simple practical guidance on the choice of temporal and spatial resolution so as to control the simulation error a priori.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Operator Relaxation and the Optimal Depth of Classical Shadows

Classical shadows are a powerful method for learning many properties of quantum states in a sample-efficient manner, by making use of randomized measurements. Here we study the sample complexity of learning the expectation value of Pauli operators via “shallow shadows,” a recently proposed version of classical shadows in which the randomization step is effected by a local unitary circuit of variable depth t. Here we show that the shadow norm (the quantity controlling the sample complexity) is expressed in terms of properties of the Heisenberg time evolution of operators under the randomizing (“twirling”) circuit—namely the evolution of the weight distribution characterizing the number of sites on which an operator acts nontrivially. For spatially contiguous Pauli operators of weight k, this entails a competition between two processes: operator spreading (whereby the support of an operator grows over time, increasing its weight) and operator relaxation (whereby the bulk of the operator develops an equilibrium density of identity operators, decreasing its weight). From this simple picture we derive (i) an upper bound on the shadow norm which, for depth t~log⁡(k), guarantees an exponential gain in sample complexity over the t=0 protocol in any spatial dimension, and (ii) quantitative results in one dimension within a mean-field approximation, including a universal subleading correction to the optimal depth, found to be in excellent agreement with infinite matrix product state numerical simulations. Our Letter connects fundamental ideas in quantum many-body dynamics to applications in quantum information science, and paves the way to highly optimized protocols for learning different properties of quantum states.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Electronic interactions in Dirac fluids visualized by nano-terahertz spacetime interference of electron-photon quasiparticles

Ultraclean graphene at charge neutrality hosts a quantum critical Dirac fluid of interacting electrons and holes. Interactions profoundly affect the charge dynamics of graphene, which is encoded in the properties of its electron-photon collective modes: surface plasmon polaritons (SPPs). Here, we show that polaritonic interference patterns are particularly well suited to unveil the interactions in Dirac fluids by tracking polaritonic interference in time at temporal scales commensurate with the electronic scattering. Spacetime SPP interference patterns recorded in terahertz (THz) frequency range provided unobstructed readouts of the group velocity and lifetime of polariton that can be directly mapped onto the electronic spectral weight and the relaxation rate. Our data uncovered prominent departures of the electron dynamics from the predictions of the conventional Fermi-liquid theory. The deviations are particularly strong when the densities of electrons and holes are approximately equal. The proposed spacetime imaging methodology can be broadly applied to probe the electrodynamics of quantum materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Reversed dynamics of bottlebrush polymers with stiff backbone and flexible side chains

The segmental dynamics of bottlebrush polymers with a stiff backbone and flexible side chains has been studied. The segmental relaxation time of side chains attached to a flexible backbone follows the same trend as linear polymers, an increase with the increasing molecular weight, but is slowed down compared to their linear counterparts. Theoretical work predicts a reversal of the molecular weight dependence of the relaxation time for stiff backbones. As a model for a stiff-g-flexible system, bottlebrushes with poly(norbornene) backbone and poly(propylene oxide) side chains, PNB-g-PPO, at a uniform grafting density have been synthesized and characterized with quasi-elastic neutron scattering. Indeed, the anticipated reversed dynamics was found. Increasing the side chain length decreases the segmental relaxation time. This indicates the importance of the characteristics of the grafting site beyond a simplified picture of an attached side chain. The mean square displacement shows a similar trend with longer side chains exhibiting a larger displacement.

Chemistry↗

Monitoring Depolymerization in Mesopores Using Dynamic Properties of Polymeric Melt Accessed via Dielectric Spectroscopy

Traditional design principles for heterogeneous catalysis guide the use of catalytic particles with mesosized (∼2–50 nm) pores to increase the number of surface-active sites by way of an increased surface area. However, the entry of long-chain polymers into such pores may be significantly limited by the size and entanglement of polymers in the melt state, thereby decreasing the number of accessible sites. Assessment of catalyst performance from traditional reactor-based studies averages over intrapore reaction events as well as reactions on the surface of a particle, resulting in an inability to distinguish between differences in site accessibility and activity. Techniques that assess the intrapore performance can inform the design of future heterogeneous catalysts for polymer upcycling. In this work, we demonstrate the use of broadband dielectric spectroscopy to monitor depolymerization of a polymer melt within mesopores via changes in the segmental relaxation time scale of amorphous polymer chains. In particular, we highlight the use of an anodic aluminum oxide (AAO) membrane as a readily available model for catalyst pores with a well-characterized pore morphology. The decrease in the segmental relaxation (α-relaxation) time of the melt with increasing chain scission emerges as a measure of the extent of polymer deconstruction inside mesopores. To demonstrate the utility of this technique, we demonstrate the decomposition of two commercial poly(propylene carbonate) polymers with different decomposition rates within mesopores. As the polymers depolymerize, their segmental relaxation time decreases as the molecular weight decreases (as predicted by the Fox–Flory equation). The BDS-measured change in segmental relaxation time mirrors the expected trend based on change in molecular weight measured by size exclusion chromatography.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Chain dynamics of alternating polymers P(C n EG 4 )

Alternating or sequence defined polymers attract the attention of an increasing number of researchers recently. Due to their different blocks, they are very customizable and material properties can be tuned. In this publication, we present dynamical studies with focus on polymer dynamics, investigated by rheology and fast field cycling (FFC) NMR. The molecular weight dependencies of the relaxation time and zero-shear viscosity could hint to entanglement effects; however, the spectral shape of the viscoelastic data resembles the polymer dynamics of unentangled melts. Furthermore, taking both techniques into account, i.e. rheology and FFC NMR leads to the conclusion that the relaxation spectra of the alternating polymer are similar to those of linear polymers.

36 MATERIALS SCIENCE↗

Segmental relaxation of sequence defined polymers

For this work, the dynamical behavior of sequence defined polymers, P(C n EG 4 ), was studied using dielectric spectroscopy showing one segmental relaxation in addition to a secondary relaxation. In case of segmental relaxation, the relaxation times strongly depend on the molecular weight at low temperatures, while at higher temperatures, unlike to linear homo polymers, this effect levels out. With increasing length of C-units, the segmental relaxation accelerates. This is also reflected in the glass transition temperature, extracted from dielectric spectroscopy. With increasing length of C-units the glass transition temperature decreases, compatible with the accelerated segmental relaxation.

36 MATERIALS SCIENCE↗

Effects of Asymmetric Molecular Architecture on Chain Stretching and Dynamics in Miktoarm Star Copolymers

We use broad-band dielectric spectroscopy (BDS) and small-angle X-ray scattering (SAXS) to investigate the impact of architectural asymmetry in the miktoarm star copolymer on the dielectric polymer relaxations. The miktoarm copolymers studied contain one or two polystyrene (PS) chains of constant molecular weight and two identical poly(cis-1,4-isoprene) (PI) chains with varied molecular weights. Using the chains in the PI block as dielectric probes, we find that the architecturally asymmetric miktoarm star copolymer systems (PSPI 2 ) feature distributions in chain relaxation times and dielectric relaxation strengths that are not dependent on molecular weight or morphology, in stark contrast to the effects of morphological confinement observed for symmetric diblock systems (PS-b-PI or PS 2 PI 2 ). Along with evidence from the SAXS measurements regarding phase separation, these results are attributed to influences from chain stretching within the framework of the Gaussian chain model for block copolymer systems. As such, the molecular architecture in block copolymers happens to be a versatile handle to control polymer chain dynamics and, ultimately, the macroscopic physicochemical properties in architecturally complex polymers.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Joule heating in bad and slow metals

Heat supplied to a metal is absorbed by the electrons and then transferred to the lattice. In conventional metals energy is released to the lattice by phonons emitted from the Lindhard continuum. However in a 'bad' metal, with short mean free path, the low energy Lindhard continuum is destroyed. To describe energy transfer to the lattice in these cases we obtain a general Kubo formula for the energy relaxation rate in terms of the electronic density spectral weight \text{lm} \, G^R_{nn}(\omega_{k},k) lm G n n R ( ω k , k ) evaluated on the phonon dispersion \omega_k ω k . We apply our Kubo formula to the high temperature Hubbard model, using recent data from quantum Monte Carlo and experiments in ultracold atoms to characterize \text{lm} \, G^R_{nn}(\omega_{k},k) lm G n n R ( ω k , k ) . We furthermore use recent data from electron energy-loss spectroscopy to estimate the energy relaxation rate of the cuprate strange metal to a high energy optical phonon. As a second, distinct, application of our formalism we consider 'slow' metals. These are defined to have Fermi velocity less than the sound velocity, so that particle-hole pairs are kinematically unable to emit phonons. We obtain an expression for the energy relaxation rate of a slow metal in terms of the optical conductivity.

36 MATERIALS SCIENCE↗

Dataset_for_Molecular_Motion_Below_the_Glass_Transition_A_Solid-State_NMR_Study_of_Siloxane_Polymer_Dynamics Study

This dataset contains solid-state 1H and 13C NMR relaxometry data, differential scanning calorimetry (DSC) data, and size exclusion chromatography (SEC/GPC) data supporting the study of sub-glass-transition (sub-Tg) molecular dynamics in a composition- and sequence-controlled series of diphenyl-substituted polysiloxanes (PDMS, 14Ph, 33Ph, 50Ph, 67Ph, and 100Ph; 0–100% diphenylsiloxane content by mole).All solid-state NMR data were acquired on a 200 MHz Bruker Avance III HD spectrometer using a static 7 mm HX probe or a 4 mm HX probe under 4 kHz magic-angle spinning. Raw Bruker TopSpin experiment folders are included for: (1) variable-temperature 1H lineshape measurements used to determine linewidth (FWHM) as a function of temperature across the glass transition; (2) 1H T1 (saturation recovery with solid-echo detection), probing nanosecond-scale dynamics near the 1H Larmor frequency; (3) 1H T1rho (direct spin-lock, 62.5 kHz), probing microsecond-scale segmental dynamics; (4) 13C-detected Lee–Goldburg cross-polarization 1H T1rho (LGCPH T1rho) for 33Ph and 50Ph, resolving aromatic and aliphatic proton environments; and (5) 13C T1 relaxation for 33Ph and 50Ph. Differential scanning calorimetry data (TA Instruments DSC 25, −150 to +120 °C, up to +300 °C for 100Ph, 10 °C/min) are included for all six compositions and support the glass-transition temperatures in Table 1 and Figure 1. Size exclusion chromatography data (Agilent 1200 Series, PL-Gel 300 mixed-C column, THF mobile phase, polystyrene calibration standards) are included for the three synthesized copolymers (33Ph, 50Ph, 67Ph) and support the number-average molecular weights in Table 1. Processed data include per-composition relaxation-time summaries (Excel), curve-fitting and Bloembergen-Purcell-Pound (BPP) model analysis notebooks (Jupyter/Python), and Igor Pro (.pxp) master files used to generate the manuscript's figures.

Bloembergen-Purcell-Pound theory↗

Improving Excited-State Potential Energy Surfaces via Optimal Orbital Shapes

In this work, we demonstrate that, rather than resorting to high-cost dynamic correlation methods, qualitative failures in excited-state potential energy surface predictions can often be remedied at no additional cost by ensuring that optimal molecular orbitals are used for each individual excited state. This approach also avoids the weighting choices required by state-averaging and dynamic weighting and obviates their need for expensive wave function response calculations when relaxing excited-state geometries. Although multistate approaches are of course preferred near conical intersections, other features of excited-state potential energy surfaces can benefit significantly from our single-state approach. In three different systems, including a double bond dissociation, a biologically relevant amino hydrogen dissociation, and an amino-to-ring intramolecular charge transfer, we show that state-specific orbitals offer qualitative improvements over the state-averaged status quo.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

On relaxations of the max k -cut problem formulations

Here, a tight continuous relaxation is a crucial factor in solving mixed integer formulations of many NP-hard combinatorial optimization problems. The (weighted) max k-cut problem is a fundamental combinatorial optimization problem with multiple notorious mixed integer optimization formulations. In this paper, we explore four existing mixed integer optimization formulations of the max k-cut problem. Specifically, we show that the continuous relaxation of a binary quadratic optimization formulation of the problem is: (i) stronger than the continuous relaxation of two mixed integer linear optimization formulations and (ii) at least as strong as the continuous relaxation of a mixed integer semidefinite optimization formulation. We also conduct a set of experiments on multiple sets of instances of the max k-cut problem using state-of-the-art solvers that empirically confirm the theoretical results in item (i). Furthermore, these numerical results illustrate the advances in the efficiency of global non-convex quadratic optimization solvers and more general mixed integer nonlinear optimization solvers. As a result, these solvers provide a promising option to solve combinatorial optimization problems. Our codes and data are available on GitHub.

97 MATHEMATICS AND COMPUTING↗

ICDARTS: Improving the Stability of Cyclic DARTS

Cyclic DARTS (CDARTS) is a Differentiable Architecture Search (DARTS)-based approach to neural architecture search (NAS) that uses a cyclic feedback mechanism to train search and evaluation networks concurrently. This training protocol aims to optimize the search process and evaluate the deep evaluation network comprised of discretized candidate operations. However, this approach introduces a loss function for the evaluation network dependent on the search network. The dissimilarity between the evaluation network’s loss function used during the search and retraining phases results in a search network that is a sub-optimal proxy for the final evaluation network accessed during retraining. We present a revised approach that removes the dependency of the evaluation network weights upon those of the search network. In addition, we introduce a modified process for relaxing the search network’s zero operations that allows these operations to be retained in the final evaluation networks.

Herron, Emily↗