Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “EXPONENTIAL FUNCTION”

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 991 records · Page 55

Generalized Unitary Coupled Cluster Wave functions for Quantum Computation

We introduce a unitary coupled-cluster (UCC) ansatz termed k-UpCCGSD that is based on a family of sparse generalized doubles operators, which provides an affordable and systematically improvable unitary coupled-cluster wave function suitable for implementation on a near-term quantum computer. k-UpCCGSD employs k products of the exponential of pair coupled-cluster double excitation operators (pCCD), together with generalized single excitation operators. We compare its performance in both efficiency of implementation and accuracy with that of the generalized UCC ansatz employing the full generalized single and double excitation operators (UCCGSD), as well as with the standard ansatz employing only single and double excitations (UCCSD). k-UpCCGSD is found to show the best scaling for quantum computing applications, requiring a circuit depth of O(kN) , compared with O(N 3 ) for UCCGSD, and O((N-η) 2 η) for UCCSD, where N is the number of spin orbitals and η is the number of electrons. We analyzed the accuracy of these three ansätze by making classical benchmark calculations on the ground state and the first excited state of H 4 (STO-3G, 6-31G), H 2 O (STO-3G), and N 2 (STO-3G), making additional comparisons to conventional coupled cluster methods. The results for ground states show that k-UpCCGSD offers a good trade-off between accuracy and cost, achieving chemical accuracy for lower cost of implementation on quantum computers than both UCCGSD and UCCSD. UCCGSD is also found to be more accurate than UCCSD but at a greater cost for implementation. Excited states are calculated with an orthogonally constrained variational quantum eigensolver approach. This is seen to generally yield less accurate energies than for the corresponding ground states. Lastly, we demonstrate that using a specialized multideterminantal reference state constructed from classical linear response calculations allows these excited state energetics to be improved.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Uplift and Seismicity Driven by Magmatic Inflation at Sierra Negra Volcano, Galápagos Islands

Abstract Although episodes of surface uplift and elevated seismicity precede many volcanic eruptions, their temporal evolution is often complex, and apparently in contradiction to simple trends predicted by mechanical deformation models. Here, we use continuous global positioning system and seismic data recorded at Sierra Negra volcano, Galápagos Islands, to show how the edifice responded to stress changes driven by magma accumulation in a shallow sill. The rate of uplift varied during the 13 years and 6.5 m of inflation before the 2018 eruption. The number of earthquakes per unit of uplift increased exponentially with total uplift as the differential stress increased. Accordingly, the temporal seismicity rate varied in time as a function of both the total uplift and the uplift rate. The Gutenberg‐Richter b ‐value decreased as a function of total uplift. In the final six months before the eruption, a sequence of large ( M > 4) earthquakes regulated the state of stress on the fault, each being followed by 2–3 days of postseismic quiescence, and retarding the increase in seismicity rate. These earthquakes did not affect the overall uplift rate. Subsidence of 8.5 m accompanied the 2‐month eruption. On resumption of uplift, the number of earthquakes per unit of uplift was very low, and the b ‐value high, reflecting the relaxed stress state of the fault system. These observations show that crustal deformation becomes increasingly brittle at higher stress states, and supports theoretical models based on elastic‐brittle mechanics. They suggest that joint interpretation of deformation and seismicity is key for forecasting future eruptions in similar volcanic settings.

Bell, Andrew Forbes↗

Understanding the role of ultra-thin polymeric interlayers in improving efficiency of polymer light emitting diodes

Insertion of ultra-thin polymeric interlayers (ILs) between the poly(3,4-ethylenedioxythiophene):polystyrene sulphonate hole injection and poly(9,9-dioctylfluorene-alt-benzothiadiazole) (F8BT) light emission layers of polymer light emitting diodes (PLEDs) can significantly increase their efficiency. In this paper, we investigate experimentally a broad range of probable causes of this enhancement with an eye to determining which IL parameters have the most significant effects. The importance of hole injection and electron blocking was studied through varying the IL material (and consequently its electronic energy levels) for both PLED and hole-only diode structures. The role of IL conductivity was examined by introducing a varying level of charge-transfer doping through blending the IL materials with a strong electron-accepting small molecule in concentrations from 1% to 7% by weight. Depositing ILs with thicknesses below the exciton diffusion length of ∼15 nm allowed the role of the IL as a physical barrier to exciton quenching to be probed. IL containing PLEDs was also fabricated with Lumation Green Series 1300 (LG 1300) light emission layers. On the other hand, the PLEDs were modeled using a 3D multi-particle Kinetic Monte Carlo simulation coupled with an optical model describing how light is extracted from the PLED. The model describes charge carrier transport and interactions between electrons, holes, singlets, and triplets, with the current density, luminance, and recombination zone (RZ) locations calculated for each PLED. The model shows F8BT PLEDs have a narrow charge RZ adjacent to the anode, while LG 1300 PLEDs have a wide charge RZ that is evenly distributed across the light emitting layer. Varying the light emitting layer from F8BT to Lumation Green Series 1300, we therefore experimentally examine the dependence of the IL function, specifically in regard to anode-side exciton quenching, on the location of the RZ. We found an exponential dependence of F8BT PLED luminance on the difference, δ, in the highest occupied to lowest unoccupied molecular orbital energy gap between the light emitting polymer and a semiconducting polymeric IL, with δ consequently the most important parameter determining efficiency. Understanding the exponential effect that wider energy gap IL materials have on exciton quenching may allow δ to be used to better guide PLED structure design.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Derivation of a Markov state model of the dynamics of a protein-like chain immersed in an implicit solvent

A Markov state model of the dynamics of a protein-like chain immersed in an implicit hard sphere solvent is derived from first principles for a system of monomers that interact via discontinuous potentials designed to account for local structure and bonding in a coarse-grained sense. The model is based on the assumption that the implicit solvent interacts on a fast time scale with the monomers of the chain compared to the time scale for structural rearrangements of the chain and provides sufficient friction so that the motion of monomers is governed by the Smoluchowski equation. A microscopic theory for the dynamics of the system is developed that reduces to a Markovian model of the kinetics under well-defined conditions. Microscopic expressions for the rate constants that appear in the Markov state model are analyzed and expressed in terms of a temperature-dependent linear combination of escape rates that themselves are independent of temperature. Excellent agreement is demonstrated between the theoretical predictions of the escape rates and those obtained through simulation of a stochastic model of the dynamics of bond formation. Finally, the Markov model is studied by analyzing the eigenvalues and eigenvectors of the matrix of transition rates, and the equilibration process for a simple helix-forming system from an ensemble of initially extended configurations to mainly folded configurations is investigated as a function of temperature for a number of different chain lengths. For short chains, the relaxation is primarily single-exponential and becomes independent of temperature in the low-temperature regime. The profile is more complicated for longer chains, where multi-exponential relaxation behavior is seen at intermediate temperatures followed by a low temperature regime in which the folding becomes rapid and single exponential. It is demonstrated that the behavior of the equilibration profile as the temperature is lowered can be understood in terms of the number of relaxation modes or “folding pathways” that contribute to the evolution of the state populations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Evolution of the age-included nearest pair distribution in disperse multiphase flows

The age of the nearest particle pair is introduced as the difference between the current time and the most recent time when the nearest particle pair was formed. The evolution equation for the age-included nearest pair distribution function is derived. With the assumption of random destruction of the nearest particle pairs, the evolution equation predicts the exponential probability distribution of the ages of the nearest particle pairs. Particle-resolved numerical simulations with moving particles are performed to verify this prediction. The equation is then used to derive the evolution equation for the particle–fluid–particle (PFP) stress, which is known to be related to hyperbolicity of the two-fluid equations. It is found that the relaxation time of the age probability distribution is also the relaxation time for the PFP stress. Guided by the closure terms in the PFP stress evolution equation, we study kinematics of the nearest particle pairs in the particle-resolved simulations for flows caused by sedimentation of the particles with initially isotropic and homogeneous particle distributions. At the steady states, the particle Reynolds numbers are around 20. Anisotropy and inhomogeneity of particle distributions are seen to develop in these flows. The mean distances to the nearest particles and evolution of the distribution of the Voronoi cell volumes are studied. We also found the PFP stress is closely related to the changes in these inter-particle scale quantities.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

A MOLECULAR STAR FORMATION LAW IN THE ATOMIC-GAS-DOMINATED REGIME IN NEARBY GALAXIES

We use the IRAM HERACLES survey to study CO emission from 33 nearby spiral galaxies down to very low intensities. Using 21 cm line atomic hydrogen (H I) data, mostly from THINGS, we predict the local mean CO velocity based on the mean H I velocity. By re-normalizing the CO velocity axis so that zero corresponds to the local mean H I velocity we are able to stack spectra coherently over large regions. This enables us to measure CO intensities with high significance as low as I{sub CO} {approx} 0.3 K km s{sup -1} ({Sigma}{sub H{sub 2}}{approx}1 M{sub sun} pc{sup -2}), an improvement of about one order of magnitude over previous studies. We detect CO out to galactocentric radii r{sub gal} {approx} r{sub 25} and find the CO radial profile to follow a remarkably uniform exponential decline with a scale length of {approx}0.2 r{sub 25}. Here we focus on stacking as a function of radius, comparing our sensitive CO profiles to matched profiles of H I, H{alpha}, far-UV (FUV), and Infrared (IR) emission at 24 {mu}m and 70 {mu}m. We observe a tight, roughly linear relationship between CO and IR intensity that does not show any notable break between regions that are dominated by molecular gas ({Sigma}{sub H{sub 2}}>{Sigma}{sub H{sub i}}) and those dominated by atomic gas ({Sigma}{sub H{sub 2}}<{Sigma}{sub H{sub i}}). We use combinations of FUV+24 {mu}m and H{alpha}+24 {mu}m to estimate the recent star formation rate (SFR) surface density, {Sigma}{sub SFR}, and find approximately linear relations between {Sigma}{sub SFR} and {Sigma}{sub H{sub 2}}. We interpret this as evidence of stars forming in molecular gas with little dependence on the local total gas surface density. While galaxies display small internal variations in the SFR-to-H{sub 2} ratio, we do observe systematic galaxy-to-galaxy variations. These galaxy-to-galaxy variations dominate the scatter in relationships between CO and SFR tracers measured at large scales. The variations have the sense that less massive galaxies exhibit larger ratios of SFR-to-CO than massive galaxies. Unlike the SFR-to-CO ratio, the balance between atomic and molecular gas depends strongly on the total gas surface density and galactocentric radius. It must also depend on additional parameters. Our results reinforce and extend to lower surface densities, a picture in which star formation in galaxies can be separated into two processes: the assembly of star-forming molecular clouds and the formation of stars from H{sub 2}. The interplay between these processes yields a total gas-SFR relation with a changing slope, which has previously been observed and identified as a star formation threshold.

79 ASTRONOMY AND ASTROPHYSICS↗

Infrared properties of z = 7 galaxies from cosmological simulations

Three-dimensional panchromatic dust radiative transfer calculations are performed on a set of 198 galaxies of stellar masses in the range of 5 × 10{sup 8} to 3 × 10{sup 10} M {sub ☉} from a cosmological hydrodynamic simulation (resolved at 29 h {sup –1} pc) at z ∼ 7. In a companion paper, the stellar mass and UV luminosity functions and the UV-optical and FUV-NUV colors are shown to be in good agreement with observations if a Small Magellanic Cloud type dust extinction curve is adopted. Here, we make useful predictions, self-consistently, of the infrared properties of these z ∼ 7 simulated galaxies that can be confronted with upcoming ALMA data. Our findings are as follows. (1) The effective radius in the rest-frame MIPS70 μm band is in the range of 80-400 pc proper for z = 7 galaxies with L {sub FIR} = 10{sup 11.3}-10{sup 12} L {sub ☉}. (2) The median of the peak wavelength of the far-infrared (FIR) spectral energy distribution is in the range of 45-60 μm, depending on the dust-to-metal ratio. (3) For star formation rate in the range of 3-100 M {sub ☉} yr{sup –1}, the median FIR to bolometric luminosity ratio is 60%-90%. (4) The FIR luminosity function displays a power law in the high end with a slope of –3.1 ± 0.4 instead of the usual exponential decline.

79 ASTRONOMY AND ASTROPHYSICS↗

Equilibration of the planar modes of ultracold two-dimensional ion crystals in a Penning trap

In this work, planar thermal equilibration is studied using direct numerical simulations of ultracold two-dimensional ion crystals in a Penning trap with a rotating wall. The large magnetic field of the trap splits the modes that describe in-plane motion of the ions into two branches: high-frequency cyclotron modes dominated by kinetic energy and low-frequency E × B modes dominated by potential energy associated with thermal position displacements. Using an eigenmode analysis we extract the equilibration rate between these two branches as a function of the ratio of the frequencies that characterize the two branches and observe this equilibration rate to be exponentially suppressed as the ratio increases. Under experimental conditions relevant for current work at NIST, the predicted equilibration time is orders of magnitude longer than any relevant experimental timescales. We also study the coupling rate dependence on the initial energy of the planar modes and the number of ions. In addition, we show how increasing the rotating wall strength improves crystal stability. These details of in-plane mode dynamics help set the stage for developing strategies to efficiently cool the in-plane modes and improve the performance of single plane ion crystals for quantum information processing.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Echoes in x-ray speckles track nanometer-scale plastic events in colloidal gels under shear

We report x-ray photon correlation spectroscopy experiments on a concentrated nanocolloidal gel subject to in situ oscillatory shear strain. The strain causes periodic echoes in the speckle pattern that lead to peaks in the intensity autocorrelation function. Above a threshold strain that is near the first yield point of the gel, the peak amplitude decays exponentially with the number of shear cycles, signaling irreversible particle rearrangements. The wave-vector dependence of the decay rate reveals a power-law distribution in the size of regions undergoing shear-induced rearrangement. The gel also displays strain softening well below the threshold, indicating a range of strains at which the rheology is nonlinear but the microscopic deformations are reversible.

Rogers, Michael C.↗

Quantum simulation of real-space dynamics

Quantum simulation is a prominent application of quantum computers. While there is extensive previous work on simulating finite-dimensional systems, less is known about quantum algorithms for real-space dynamics. We conduct a systematic study of such algorithms. In particular, we show that the dynamics of a d-dimensional Schrödinger equation with η particles can be simulated with gate complexity O ~ (ηdFpoly(log(g'/ϵ))), where ϵ is the discretization error, g' controls the higher-order derivatives of the wave function, and F measures the time-integrated strength of the potential. Compared to the best previous results, this exponentially improves the dependence on ϵ and g' from poly(g'/ϵ) to poly(log(g'/ϵ)) and polynomially improves the dependence on T and d, while maintaining best known performance with respect to η. For the case of Coulomb interactions, we give an algorithm using η 3 (d + η)Tpoly(log(ηdTg'/(Δϵ)))/Δ one- and two-qubit gates, and another using η 3 (4d) d/2 Tpoly(log(ηdTg'/(Δϵ)))/Δ one- and two-qubit gates and QRAM operations, where T is the evolution time and the parameter Δ regulates the unbounded Coulomb interaction. We give applications to several computational problems, including faster real-space simulation of quantum chemistry, rigorous analysis of discretization error for simulation of a uniform electron gas, and a quadratic improvement to a quantum algorithm for escaping saddle points in nonconvex optimization.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Radar scattering laws and wavelength dependence of the lunar surface

Data from Apollo lunar bistatic radar experiments have been processed to give probability density functions for surface slopes. These show best agreement with a Hagfors scattering law, though data having both gaussian and exponential characteristics also exist. Surface roughness estimates range from 4 deg in maria to at least 8 deg in highlands, values which are appropriate to 25 m horizontal scales and which are areal averages over tens of square kilometers. Roughness varies with wavelength, most strongly in maria.

Simpson, R. A.↗

Design of vibration isolation systems using multiobjective optimization techniques

The design of vibration isolation systems is considered using multicriteria optimization techniques. The integrated values of the square of the force transmitted to the main mass and the square of the relative displacement between the main mass and the base are taken as the performance indices. The design of a three degrees-of-freedom isolation system with an exponentially decaying type of base disturbance is considered for illustration. Numerical results are obtained using the global criterion, utility function, bounded objective, lexicographic, goal programming, goal attainment and game theory methods. It is found that the game theory approach is superior in finding a better optimum solution with proper balance of the various objective functions.

Rao, S. S.↗

Bimodal star formation - Constraints from the solar neighborhood

The chemical evolution resulting from a simple model of bimodal star formulation is investigated, using constraints from the solar neighborhood to set the parameters of the initial mass function and star formation rate. The two modes are an exclusively massive star mode, which forms stars at an exponentially declining rate, and a mode which contains stars of all masses and has a constant star formation rate. Satisfactory agreement with the age-metallicity relation for the thin disk and with the metallicity structure of the thin-disk and spheroid stars is possible only for a small range of parameter values. The preferred model offers a resolution to several of the long-standing problems of galactic chemical evolution, including explanations of the age-metallicity relation, the gas consumption time scale, and the stellar cumulative metallicity distributions.

Wyse, Rosemary F. G.↗

Statistical properties of Charney-Hasegawa-Mima zonal flows

A theoretical interpretation of numerically generated probability density functions (PDFs) of intermittent plasma transport events in unforced zonal flows is provided within the Charney-Hasegawa-Mima (CHM) model. The governing equation is solved numerically with various prescribed density gradients that are designed to produce different configurations of parallel and anti-parallel streams. Long-lasting vortices form whose flow is governed by the zonal streams. It is found that the numerically generated PDFs can be matched with analytical predictions of PDFs based on the instanton method by removing the autocorrelations from the time series. In many instances, the statistics generated by the CHM dynamics relaxes to Gaussian distributions for both the electrostatic and vorticity perturbations, whereas in areas with strong nonlinear interactions it is found that the PDFs are exponentially distributed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Modulation of lens cell adhesion molecules by particle beams

Cell adhesion molecules (CAMs) are proteins which anchor cells to each other and to the extracellular matrix (ECM), but whose functions also include signal transduction, differentiation, and apoptosis. We are testing a hypothesis that particle radiations modulate CAM expression and this contributes to radiation-induced lens opacification. We observed dose-dependent changes in the expression of beta 1-integrin and ICAM-1 in exponentially-growing and confluent cells of a differentiating human lens epithelial cell model after exposure to particle beams. Human lens epithelial (HLE) cells, less than 10 passages after their initial culture from fetal tissue, were grown on bovine corneal endothelial cell-derived ECM in medium containing 15% fetal bovine serum and supplemented with 5 ng/ml basic fibroblast growth factor (FGF-2). Multiple cell populations at three different stages of differentiation were prepared for experiment: cells in exponential growth, and cells at 5 and 10 days post-confluence. The differentiation status of cells was characterized morphologically by digital image analysis, and biochemically by Western blotting using lens epithelial and fiber cell-specific markers. Cultures were irradiated with single doses (4, 8 or 12 Gy) of 55 MeV protons and, along with unirradiated control samples, were fixed using -20 degrees C methanol at 6 hours after exposure. Replicate experiments and similar experiments with helium ions are in progress. The intracellular localization of beta 1-integrin and ICAM-1 was detected by immunofluorescence using monoclonal antibodies specific for each CAM. Cells known to express each CAM were also processed as positive controls. Both exponentially-growing and confluent, differentiating cells demonstrated a dramatic proton-dose-dependent modulation (upregulation for exponential cells, downregulation for confluent cells) and a change in the intracellular distribution of the beta 1-integrin, compared to unirradiated controls. In contrast, there was a dose-dependent increase in ICAM-1 immunofluorescence in confluent, but not exponentially-growing cells. These results suggest that proton irradiation downregulates beta 1-integrin and upregulates ICAM-1, potentially contributing to cell death or to aberrant differentiation via modulation of anchorage and/or signal transduction functions. Quantification of the expression levels of the CAMs by Western analysis is in progress.

NASA Discipline Radiation Health↗

Controlling Combinatorial Complexity in Software and Malware Behavior Computation

Virtually all software is out of intellectual control in that no one knows its full behavior. Software Behavior Computation (SBC) is a new technology for understanding everything software does. SBC applies the mathematics of denotational semantics implemented by function composition in Functional Trace Tables (FTTs) to compute the behavior of programs, expressed as disjoint cases of conditional concurrent assignments. In some circumstances, combinatorial explosions in the number of cases can occur when calculating the behavior of sequences of multiple branching structures. This paper describes computational methods that avoid combinatorial explosions. The predicates that control branching structures such as ifthenelses can be organized into three categories: 1) Independent, resulting in no behavior case explosion, 2) Coordinated, resulting in two behavior cases, or 3) Goaloriented, with potential exponential growth in the number of cases. Traditional FTT-based behavior computation can be augmented by two additional computational methods, namely, Single-Value Function Abstractions (SVFAs) and, introduced in this paper, Relational Trace Tables (RTTs). These methods can be applied to the three predicate categories to avoid combinatorial growth in behavior cases while maintaining mathematical correctness.

Hyperion system↗

Vortex flux dynamics and harmonic ac magnetic response of Ba(Fe 0.94 Ni 0.06 ) 2 As 2 bulk superconductor

Vortex dynamics and nonlinear ac response are studied in a Ba(Fe 0.94 Ni 0.06 ) 2 As 2 ( T c = 18.5 K) bulk superconductor in magnetic fields up to 12 T via ac susceptibility measurements of the first ten harmonics. A comprehensive study of the ac magnetic susceptibility and its first ten harmonics finds shifts to higher temperatures with increasing ac measurement frequencies (10 to 10,000 Hz) for a wide range of ac (1, 5, and 10 Oe) and dc fields (0 to 12 T). The characteristic measurement time constant t1 is extracted from the exponential fit of the data and linked to vortex relaxation. The Anderson-Kim Arrhenius law is applied to determine flux activation energy E a /k as a function dc magnetic field. The de-pinning, or irreversibility lines, were determined by a variety of methods and extensively mapped. The ac response shows surprisingly weak higher harmonic components, suggesting weak nonlinear behavior. Lastly, our data does not support the Fisher model; we do not see an abrupt vortex glass to vortex liquid transition and the resistivity does not drop to zero, although it appears to approach zero exponentially.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Electronic couplings for molecular charge transfer: Benchmarking CDFT, FODFT, and FODFTB against high-level ab initio calculations

We introduce a database (HAB11) of electronic coupling matrix elements (H{sub ab}) for electron transfer in 11 π-conjugated organic homo-dimer cations. High-level ab inito calculations at the multireference configuration interaction MRCI+Q level of theory, n-electron valence state perturbation theory NEVPT2, and (spin-component scaled) approximate coupled cluster model (SCS)-CC2 are reported for this database to assess the performance of three DFT methods of decreasing computational cost, including constrained density functional theory (CDFT), fragment-orbital DFT (FODFT), and self-consistent charge density functional tight-binding (FODFTB). We find that the CDFT approach in combination with a modified PBE functional containing 50% Hartree-Fock exchange gives best results for absolute H{sub ab} values (mean relative unsigned error = 5.3%) and exponential distance decay constants β (4.3%). CDFT in combination with pure PBE overestimates couplings by 38.7% due to a too diffuse excess charge distribution, whereas the economic FODFT and highly cost-effective FODFTB methods underestimate couplings by 37.6% and 42.4%, respectively, due to neglect of interaction between donor and acceptor. The errors are systematic, however, and can be significantly reduced by applying a uniform scaling factor for each method. Applications to dimers outside the database, specifically rotated thiophene dimers and larger acenes up to pentacene, suggests that the same scaling procedure significantly improves the FODFT and FODFTB results for larger π-conjugated systems relevant to organic semiconductors and DNA.

74 ATOMIC AND MOLECULAR PHYSICS↗