Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Function approximation”

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 1,657 records · Page 92

Electronic structure, magnetic correlations, and superconducting pairing in the reduced Ruddlesden-Popper bilayer La 3 Ni 2 O 6 under pressure: Different role of d 3 z 2 - r 2 orbital compared with La 3 Ni 2 O 7

The recent discovery of superconductivity in bilayer La 3 Ni 2 O 7 (327-LNO) under pressure stimulated much interest in layered nickelates. However, superconductivity was not found in another bilayer nickelate system, La 3 Ni 2 O 6 (326-LNO), even under pressure. To understand the similarities and differences between 326-LNO and 327-LNO, using density functional theory and the random phase approximation (RPA), we systematically investigate 326-LNO under pressure. The large crystal-field splitting between the e g orbitals caused by the missing apical oxygen moves the d 3z 2 -r 2 orbital farther away from the Fermi level, implying that the d 3z 2 -r 2 orbital plays a less important role in 326-LNO than in 327-LNO. This also results in a smaller bandwidth for the d x 2 -y 2 orbital and a reduced energy gap for the bonding-antibonding splitting of the d 3z 2 -r 2 orbital in 326-LNO, as compared to 327-LNO. Moreover, the in-plane hybridization between the d x 2 -y 2 and d 3z 2 -r 2 orbitals is found to be small in 326-LNO, while it is much stronger in 327-LNO. Furthermore, the low-spin ferromagnetic state is found to be the likely ground state in 326-LNO under high pressure. The weak interlayer coupling suggests that s ± -wave pairing is unlikely in 326-LNO. The robust in-plane ferromagnetic coupling also suggests that d-wave superconductivity, which is usually caused by antiferromagnetic fluctuations of the d x 2 -y 2 orbital, is also unlikely in 326-LNO. These conclusions are supported by our many-body RPA calculations of the pairing behavior. Additionally, contrasting with the cuprates, for the bilayer cuprate HgBa 2 CaCu 2 O 6 , we find a strong self-doping effect of the d x 2 -y 2 orbital under pressure, with the charge of Cu being reduced by approximately 0.13 electrons from 0 GPa to 25 GPa. In contrast, we do not observe such a change in the electronic density in 326-LNO under pressure, establishing another important difference between the nickelates and the cuprates.

36 MATERIALS SCIENCE↗

Robust s ± pairing in CaK ( Fe 1 – x Ni x ) 4 As 4 ( x = 0 and 0.05) from the response to electron irradiation

Here, controlled pointlike disorder introduced by 2.5-MeV electron irradiation was used to probe the superconducting state of single crystals of CaK(Fe 1–x Ni x ) 4 As 4 superconductor at x = 0 and 0.05 doping levels. Both compositions show an increase of the residual resistivity and a decrease of the superconducting transition temperature, T c , at the rate of dT c /dρ(T c ) ≈ 0.19 K/(μΩcm) for x = 0 and 0.38 K/(μΩcm) for x = 0.05, respectively. In the Ni-doped compound (x = 0.05), the coexisting spin-vortex crystal (SVC) magnetic phase is suppressed at the rate of dT N /dρ(T N ) ≈ 0.16 K/(μΩcm). The low-temperature variation of London penetration depth is well approximated by the power-law function, Δλ(T) = AT n , with n ≈ 2.5 for x = 0 and n ≈ 1.9 for x = 0.05 in the pristine state. Detailed analysis of λ(T) and T c evolution with disorder is consistent with two effective nodeless energy gaps in the density of states due to robust s± pairing. Overall the behavior of CaK(Fe 1–x Ni x ) 4 As 4 at x = 0 is similar to a slightly overdoped Ba 1–y K y Fe 2 As 2 at y ≈ 0.5, and at x = 0.05 to an underdoped composition at y ≈ 0.2.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Pion and kaon valence quark distribution functions from Dyson-Schwinger equations

Using realistic quark propagators and meson Bethe-Salpeter amplitudes based on the Dyson-Schwinger equations, we calculate the pion and kaon's valence parton distribution functions (PDF) through the modified impulse approximation. The PDFs we obtained at hadronic scale have the purely valence characteristic and exhibit both dynamical chiral symmetry breaking and SU(3) flavor symmetry breaking effects. A new calculation technique is introduced to determine the valence PDFs with precision. Through next-to-leading order (NLO) evolution, our result is compared with pion and kaon valence PDF data at experimental scale. Good agreement is found in the case of pion, while deviation emerges for kaon. We point out this situation can potentially be resolved by incorporating gluon contributions into the mesons if the pion hosts more gluons than kaon nonperturbatively.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Logarithmic Corrections to Kerr Thermodynamics

Recent work has shown that loop corrections from massless particles generate $\frac{3}{2}$ ⁢log⁡𝑇 Hawking corrections to black hole entropy which dominate the thermodynamics of cold near-extreme charged black holes. Here we adapt this analysis to near-extreme Kerr black holes. Like AdS 2 ×𝑆 2 , the near-horizon extreme Kerr (NHEK) metric has a family of normalizable zero modes corresponding to reparametrizations of boundary time. The path integral over these zero modes leads to an infrared divergence in the one-loop approximation to the Euclidean NHEK partition function. We regulate this divergence by retaining the leading finite temperature correction in the NHEK scaling limit. This “not-NHEK” geometry lifts the eigenvalues of the zero modes, rendering the path integral infrared finite. The quantum-corrected near-extremal entropy exhibits $\frac{3}{2}$ ⁢log⁡𝑇 Hawking behavior characteristic of the Schwarzian model and predicts a lifting of the ground state degeneracy for the extremal Kerr black hole.

General relativity↗

Neural refinement of sample weights

Monte Carlo simulations are an essential tool in particle physics data analysis. Events are typically generated alongside weights that redistribute the cross section of the simulated process across the phase space. These weights can be negative, and several post hoc methods have been developed to eliminate or mitigate the negative values. All of these methods share the common strategy of approximating the average weight as a function of phase space. We introduce an alternative approach, which, instead of reweighting to the average, refines the initial weights with a scaling transformation, utilizing a phase space-dependent factor. Since this new refinement method does not need to model the full weight distribution, it can be more accurate. High-dimensional and unbinned phase space is processed using neural networks for the refinement method. In addition to the refinement method, we introduce a new resampling protocol, which can be used in conjunction with any weight transformation to not only preserve the average weight but also the statistical uncertainties of the initial distribution. Using both realistic and synthetic examples, we show that the new neural refinement method is able to match or exceed the accuracy of similar weight transformations and that the new resampling protocol is simpler in implementation than previous methods while exhibiting equivalent statistical properties.

Artificial neural networks↗

Towards Highly Scalable Ab Initio Molecular Dynamics (AIMD) Simulations on the Intel Knights Landing Manycore Processor

The Ab Initio Molecular Dynamics (AIMD) method allows scientists to treat the dynamics of molecular and condensed phase systems while retaining a first-principles-based description of their interactions. This extremely important method has tremendous computational requirements, because the electronic Schr¨odinger equation, approximated using Kohn-Sham Density Functional Theory (DFT), is solved at every time step. With the advent of manycore architectures, application developers have a significant amount of processing power within each compute node that can only be exploited through massive parallelism. A compute intensive application such as AIMD forms a good candidate to leverage this processing power. In this paper, we focus on adding thread level parallelism to the plane wave DFT methodology implemented in NWChem. Through a careful optimization of tall-skinny matrix products, which are at the heart of the Lagrange multiplier and nonlocal pseudopotential kernels, as well as 3D FFTs, our OpenMP implementation delivers excellent strong scaling on the latest Intel Knights Landing (KNL) processor. We assess the efficiency of our Lagrange multiplier kernels by building a Roofline model of the platform, and verify that our implementation is close to the roofline for various problem sizes. Finally, we present strong scaling results on the complete AIMD simulation for a 64 water molecules test case, that scales up to all 68 cores of the Knights Landing processor.

Environmental Molecular Sciences Laboratory↗

Decentralized Filtering Adaptive Neural Network Control for Uncertain Switched Interconnected Nonlinear Systems

This article presents a novel decentralized filtering adaptive neural network control framework for uncertain switched interconnected nonlinear systems. Each subsystem has its own decentralized controller based on the established decentralized state predictor. For each subsystem, the nonlinear uncertainties are approximated by a Gaussian radial basis function (GRBF) neural network incorporated with a piecewise constant adaptive law, where the adaptive law will update adaptive parameters from the error dynamics between the host system and the decentralized state predictor by discarding the unknowns, whereas a decentralized filtering control law is derived to cancel both local and mismatched uncertainties from other subsystems, as well as achieve the local objective tracking of the host system. The achievement of global objective depends on the achievement of local objective for each subsystem. The matched uncertainties are canceled directly by adopting their opposite in the control signal, whereas a dynamic inversion of the system is required to eliminate the effect of the mismatched uncertainties on the output. By exploiting the average dwell time principle, the error bounds between the real system and the virtual reference system, which defines the best performance that can be achieved by the closed-loop system, are derived. A numerical example is given to illustrate the effectiveness of the decentralized filtering adaptive neural network control architecture by comparing against the model reference adaptive control (MRAC).

Average dwell time, decentralized, filtering adapt↗

Multi-Stage Stochastic Programming to Joint Economic Dispatch for Energy and Reserve With Uncertain Renewable Energy

To address the uncertain renewable energy in the day-ahead optimal dispatch of energy and reserve, a multi-stage stochastic programming model is established in this paper to minimize the expected total costs. The uncertainties over the multiple stages are characterized by a scenario tree and the optimal dispatch scheme is cast as a decision tree which guarantees the flexibility to decide the reasonable outputs of generation and the adequate reserves accounting for different realizations of renewable energy. Most importantly, to deal with the “Curse of Dimensionality” of stochastic programming, stochastic dual dynamic programming (SDDP) is employed, which decomposes the original problem into several sub-problems according to the stages. Specifically, the SDDP algorithm performs forward pass and backward pass repeatedly until the convergence criterion is satisfied. At each iteration, the original problem is approximated by creating a linear piecewise function. Besides, an improved convergence criterion is adopted to narrow the optimization gaps. The results on the IEEE 118-bus system and real-life provincial power grid show the effectiveness of the proposed model and method.

economic dispatch↗

Root angle in maize influences nitrogen capture and is regulated by calcineurin B-like protein (CBL)-interacting serine/threonine-protein kinase 15 ( ZmCIPK15 )

Crops with reduced nutrient and water requirements are urgently needed in global agriculture. Root growth angle plays an important role in nutrient and water acquisition. A maize diversity panel of 481 genotypes was screened for variation in root angle employing a high-throughput field phenotyping platform. Genome-wide association mapping identified several single nucleotide polymorphisms (SNPs) associated with root angle, including one located in the root expressed CBL-interacting serine/threonine-protein kinase 15 (ZmCIPK15) gene (LOC100285495). Reverse genetic studies validated the functional importance of ZmCIPK15, causing a approximately 10° change in root angle in specific nodal positions. A steeper root growth angle improved nitrogen capture in silico and in the field. OpenSimRoot simulations predicted at 40 days of growth that this change in angle would improve nitrogen uptake by 11% and plant biomass by 4% in low nitrogen conditions. In field studies under suboptimal N availability, the cipk15 mutant with steeper growth angles had 18% greater shoot biomass and 29% greater shoot nitrogen accumulation compared to the wild type after 70 days of growth. We propose that a steeper root growth angle modulated by ZmCIPK15 will facilitate efforts to develop new crop varieties with optimal root architecture for improved performance under edaphic stress.

59 BASIC BIOLOGICAL SCIENCES↗

Distribution of ion current density on a rotating spherical cap substrate during ion-assisted deposition

The uniformity of ion density is critical for applications relying on the ion assisted deposition technique for the fabrication of the high quality thin films. The authors propose and describe here a method allowing one to calculate the ion density distribution on spherical substrate holders under stationary and rotating conditions for different positions of the ion source. The ion beam shape was approximated by a cos{sup n} function, and the ion current density was represented by a function inversely proportional to the distance from the ion source in accordance with our experimental results. As an example, a calculation of the current density distribution on the spherical cap substrate was performed for a broad beam ion source operated with an anode current of 3 A. The authors propose an approach for process optimization with respect to the ion source position and its inclination, in terms of uniformity and absolute value of the ion current density.

42 ENGINEERING↗

Design and performance of Soft Gamma-ray Detector onboard the Hitomi (ASTRO-H) satellite

Hitomi (ASTRO-H) was the sixth Japanese x-ray satellite that carried instruments with exquisite energy resolution of <7 eV and broad energy coverage of 0.3 to 600 keV. The Soft Gamma-ray Detector (SGD) was the Hitomi instrument that observed the highest energy band (60 to 600 keV). The SGD design achieves a low background level by combining active shields and Compton cameras where Compton kinematics is utilized to reject backgrounds coming from outside of the field of view. A compact and highly efficient Compton camera is realized using a combination of silicon and cadmium telluride semiconductor sensors with a good energy resolution. Compton kinematics also carries information for gamma-ray polarization, making the SGD an excellent polarimeter. Following several years of development, the satellite was successfully launched on February 17, 2016. After proper functionality of the SGD components were verified, the nominal observation mode was initiated on March 24, 2016. The SGD observed the Crab Nebula for approximately two hours before the spacecraft ceased to function on March 26, 2016. We present concepts of the SGD design followed by detailed description of the instrument and its performance measured on ground and in orbit.

79 ASTRONOMY AND ASTROPHYSICS↗

Ab initio GGA+U investigations of the structural, electronic, and magnetic properties of Cd 1-x Mn x Te alloy

The structural, electronic and magnetic properties of Cd 1-x Mn x Te alloy in the zinc-blende (ZB) phase within the framework of spin-polarized density functional theory using the generalized gradient approximation (GGA) are studied. In particular, we employ GGA+U SIC method that gives a better description of systems with strong Coulomb correlations of the d-electrons. First, we discuss implementation of GGA+U method to the ground state calculations of pure CdTe and ferromagnetic MnTe. The alloy is modelled at some selected compositions by ordered structures described in terms of periodically repeated supercells for the compositions x = 0.25, 0.5, and 0.75. The compositional dependence of the lattice constant, electronic band structure and partial densities of states of ferromagnetic Cd 1-x Mn x Te alloy are studied. We estimate the spin-exchange splitting energies produced by the Mn 3d-states and determine the exchange constants for the conduction and valence bands of the alloy. The energy positions of the occupied and unoccupied Mn 3d-bands in the electronic structure of ferromagnetic Cd 1-x Mn x Te alloy are also presented. Our calculations using GGA+U SIC approach agree well with the available experimental data and other calculations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Computer-aided pulmonary image analysis in small animal models

Purpose: To develop an automated pulmonary image analysis framework for infectious lung diseases in small animal models. Methods: The authors describe a novel pathological lung and airway segmentation method for small animals. The proposed framework includes identification of abnormal imaging patterns pertaining to infectious lung diseases. First, the authors’ system estimates an expected lung volume by utilizing a regression function between total lung capacity and approximated rib cage volume. A significant difference between the expected lung volume and the initial lung segmentation indicates the presence of severe pathology, and invokes a machine learning based abnormal imaging pattern detection system next. The final stage of the proposed framework is the automatic extraction of airway tree for which new affinity relationships within the fuzzy connectedness image segmentation framework are proposed by combining Hessian and gray-scale morphological reconstruction filters. Results: 133 CT scans were collected from four different studies encompassing a wide spectrum of pulmonary abnormalities pertaining to two commonly used small animal models (ferret and rabbit). Sensitivity and specificity were greater than 90% for pathological lung segmentation (average dice similarity coefficient > 0.9). While qualitative visual assessments of airway tree extraction were performed by the participating expert radiologists, for quantitative evaluation the authors validated the proposed airway extraction method by using publicly available EXACT’09 data set. Conclusions: The authors developed a comprehensive computer-aided pulmonary image analysis framework for preclinical research applications. The proposed framework consists of automatic pathological lung segmentation and accurate airway tree extraction. The framework has high sensitivity and specificity; therefore, it can contribute advances in preclinical research in pulmonary diseases.

60 APPLIED LIFE SCIENCES↗

WE-AB-207B-11: Optimizing Tumor Control Probability in Radiation Therapy Treatment - Application to HDR Cervical Cancer

Purpose: The ultimate goal of radiotherapy treatment planning is to find a treatment that will yield a high tumor-control-probability(TCP) with an acceptable normal-tissue-complication probability(NTCP). Yet most treatment planning today is not based upon optimization of TCPs and NTCPs, but rather upon meeting physical dose and volume constraints defined by the planner. We design treatment plans that optimize TCP directly and contrast them with the clinical dose-based plans. PET image is incorporated to evaluate gain in TCP for dose escalation. Methods: We build a nonlinear mixed integer programming optimization model that maximizes TCP directly while satisfying the dose requirements on the targeted organ and healthy tissues. The solution strategy first fits the TCP function with a piecewise-linear approximation, then solves the problem that maximizes the piecewise linear approximation of TCP, and finally performs a local neighborhood search to improve the TCP value. To gauge the feasibility, characteristics, and potential benefit of PET-image guided dose escalation, initial validation consists of fifteen cervical cancer HDR patient cases. These patients have all received prior 45Gy of external radiation dose. For both escalated strategies, we consider 35Gy PTV-dose, and two variations (37Gy-boost to BTV vs 40Gy-boost) to PET-image-pockets. Results: TCP for standard clinical plans range from 59.4% - 63.6%. TCP for dose-based PET-guided escalated-dose-plan ranges from 63.8%–98.6% for all patients; whereas TCP-optimized plans achieves over 91% for all patients. There is marginal difference in TCP among those with 37Gy-boosted vs 40Gy-boosted. There is no increase in rectum and bladder dose among all plans. Conclusion: Optimizing TCP directly results in highly conformed treatment plans. The TCP-optimized plan is individualized based on the biological PET-image of the patients. The TCP-optimization framework is generalizable and has been applied successfully to other external-beam delivery modalities. A clinical trial is on-going to gauge the clinical significance. Partially supported by the National Science Foundation.

60 APPLIED LIFE SCIENCES↗

Theoretical study of oxygen sorption and diffusion in the volume and on the surface of a γ-TiAl alloy

The oxygen sorption on the low-index (001), (100), and (110) surfaces of a γ-TiAl alloy is studied by the pseudopotential method with the generalized gradient approximation for the exchange-correlation functional. The most preferred sites for oxygen sorption in the bulk and on the surface of the alloy are determined. The titanium-rich octahedral site is shown to be preferred for oxygen sorption in the bulk material. The effect of the oxygen concentration on the atomic and electronic structures of the stoichiometric TiAl(100) surface is studied. It is shown that, at the first stage of oxidation, oxygen prefers to form bonds with titanium. The energy barriers for oxygen diffusion on the stoichiometric (100) surface and in the bulk of the material are calculated. The energy barriers are shown to depend substantially on the local environments of oxygen and to increase during diffusion from titanium-rich sites. The most possible mechanism of oxygen diffusion from the (100) surface to the bulk of the material is oxygen migration through tetrahedral sites.

36 MATERIALS SCIENCE↗

Trinal Nature of the Excitation of Bulk Plasmons by a Fast Charged Particle at Grazing Incidence on the Interface between Vacuum and a Metal with Strong Spatial Dispersion

The excitation and propagation of bulk and surface (surface waves of transition radiation of plasmons at frequencies above the plasma frequency) plasma waves by incident electrons moving both in vacuum toward the surface of a metal and inside the metal whose boundary elastically and spectacularly reflects internal nonequilibrium electrons have been analyzed. In contrast to work [B. N. Libenson, J. Exp. Theor. Phys. 113, 553 (2011)], attention in this work is focused on the influence of surface effects on bulk-plasmon excitation by incident electrons. The probabilities and spectra of single characteristic energy loss of an intermediate- energy electron (50–500 eV) moving at an angle to the surface of the medium in three regions: in vacuum, in the medium, and again in vacuum after the electron leaves the medium are calculated. The kinetic approximation is used for the dielectric function, where the entire range of the plasmon spectrum is taken into account correctly for the problem under consideration. In the indicated energy range of incident electrons, surface effects, on the one hand, significantly reduces the probability of excitation of bulk plasma waves in the medium with strong spatial dispersion, in particular, as compared to the results obtained in [B.N. Libenson, J. Exp. Theor. Phys. 113, 553 (2011)], where surface effects were disregarded and the probability of bulkplasmon excitation by a 200-eV electron incident and emitted perpendicularly to the boundary is about one third of that in the unbounded medium. On the other hand, at grazing incidence from vacuum, the probability of transition radiation of bulk plasmons increases significantly and can lead to a change in the character of the angular dependence of the intensity of bulk plasma energy loss. Thus, the main result of this work is that a decrease in the glancing angle of the fast electron with respect to the vacuum–metal interface is accompanied both by an increase in the contribution from the transition mechanism to the probability of bulk-plasmon excitation in the vacuum region and by a decrease in the contribution from Cherenkov and bremsstrahlung mechanisms of excitation in the medium. The probability of bulk-plasmon excitation in the vacuum region exceeds the probability of excitation at the further motion of the electron in metallic aluminum at angles of incidence larger than 65°, 70°, and 75° at the energies E = 200, 350, and 500 eV, respectively.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗