Search NASA⌕ Search

SEARCH · Search NASA

Results for “Analytic 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 127 records · Page 7

The four-loop remainder function and multi-Regge behavior at NNLLA in planar $ \mathcal{N} $ = 4 super-Yang-Mills theory

We present the four-loop remainder function for six-gluon scattering with maximal helicity violation in planar NN = 4 super-Yang-Mills theory, as an analytic function of three dual-conformal cross ratios. The function is constructed entirely from its analytic properties, without ever inspecting any multi-loop integrand. We employ the same approach used at three loops, writing an ansatz in terms of hexagon functions, and fixing coefficients in the ansatz using the multi-Regge limit and the operator product expansion in the near-collinear limit. We express the result in terms of multiple polylogarithms, and in terms of the coproduct for the associated Hopf algebra. From the remainder function, we extract the BFKL eigenvalue at next-to-next-to-leading logarithmic accuracy (NNLLA), and the impact factor at N3LLA. We plot the remainder function along various lines and on one surface, studying ratios of successive loop orders. As seen previously through three loops, these ratios are surprisingly constant over large regions in the space of cross ratios, and they are not far from the value expected at asymptotically large orders of perturbation theory.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

SU-E-T-517: Analytic Formalism to Compute in Real Time Dose Distributions Delivered by HDR Units

Purpose: Develop an analytical algorithm to compute the dose delivered by Ir-192 dwell positions with high accuracy using the 3-dimensional (3D) dose distribution of an HDR source. Using our analytical function, the dose delivered by an HDR unit as treatment progresses can be determined using the actual delivered temporal and positional data of each individual dwell. Consequently, true delivered dose can be computed when each catheter becomes active. We hypothesize that the knowledge of such analytical formulation will allow developing HDR systems with a real time treatment evaluation tool to avoid mistreatments. Methods: In our analytic formulation, the dose is computed by using the full anisotropic function data of the TG 43 formalism with 3D ellipsoidal function. The discrepancy between the planned dose and the delivered dose is computed using an analytic perturbation method over the initial dose distribution. This methodology speeds up the computation because only changes in dose discrepancies originated by spatial and temporal deviations are computed. A dose difference map at the point of interest is obtained from these functions and this difference can be shown during treatment in real time to examine the treatment accuracy. Results: We determine the analytical solution and a perturbation function for the 3 translational 3 rotational, and 1D temporal errors in source distributions. The analytic formulation is a sequence of simple equations that can be processed in any modern computer in few seconds. Because computations are based in an analytical solution, small deviations of the dose when sub-millimeter positional changes occur can be detected. Conclusions: We formulated an analytical method to compute 4D dose distributions and dose differences based on an analytical solution and perturbations to the original dose. This method is highly accurate and can be.

60 APPLIED LIFE SCIENCES↗

Charged current Drell-Yan production at N 3 LO

We present the production cross section for a lepton-neutrino pair at the Large Hadron Collider computed at next-to-next-to-next-to-leading order (N 3 LO) in QCD perturbation theory. We compute the partonic coefficient functions of a virtual W± boson at this order. We then use these analytic functions to study the progression of the perturbative series in different observables. In particular, we investigate the impact of the newly obtained corrections on the inclusive production cross section of W ± bosons, as well as on the ratios of the production cross sections for W + , W – and/or a virtual photon. Finally, we present N3LO predictions for the charge asymmetry at the LHC.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Ex Machina Determination of Structural Correlation Functions

Determining the structural properties of condensed-phase systems is a fundamental problem in theoretical statistical mechanics. In this study, we present a machine learning method that is able to predict structural correlation functions with significantly improved accuracy in comparison with traditional approaches. The usefulness of this ex machina (from the machine) approach is illustrated by predicting the radial distribution functions of two paradigmatic condensed-phase systems, a Lennard-Jones fluid and a hard-sphere fluid, and then comparing those results to the results obtained using both integral equation methods and empirically motivated analytical functions. We find that application of the developed ex machina method typically decreases the predictive error by more than an order of magnitude in comparison with traditional theoretical methods.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Characters, quasinormal modes, and Schwinger pairs in dS 2 with flux

An integral representation of the 1-loop partition function for charged scalars and spinors, minimally coupled to a uniform U(1) field on S 2 , is given in terms of SO(1, 2) Harish-Chandra group characters and evaluated exactly in terms of Hurwitz ζ-functions. Analytically continuing the U(1) field, we interpret the path integrals as quasicanonical partition functions in dS 2 with an electric field. The character itself is obtained as a trace over states living at the future boundary of de Sitter and has a quasinormal mode expansion. The imaginary part of the partition function captures Schwinger pair creation in the static patch at finite temperature. The thermal enhancement is most noticeable for scalar masses below Hubble and leads to non-monotonicity of the current as a function of the field. This parameter range, when dimensionally reducing from a charged or rotating Nariai spacetime, is excluded by Swampland-inspired bounds. Around the AdS 2 black hole, in contrast to dS 2 , there is a threshold to pair creation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Stability of uncertain systems

The asymptotic properties of feedback systems are discussed, containing uncertain parameters and subjected to stochastic perturbations. The approach is functional analytic in flavor and thereby avoids the use of Markov techniques and auxiliary Lyapunov functionals characteristic of the existing work in this area. The results are given for the probability distributions of the accessible signals in the system and are proved using the Prohorov theory of the convergence of measures. For general nonlinear systems, a result similar to the small loop-gain theorem of deterministic stability theory is given. Boundedness is a property of the induced distributions of the signals and not the usual notion of boundedness in norm. For the special class of feedback systems formed by the cascade of a white noise, a sector nonlinearity and convolution operator conditions are given to insure the total boundedness of the overall feedback system.

Blankenship, G. L.↗

Kinetic corrections from analytic non-Maxwellian distribution functions in magnetized plasmas

In magnetized plasma physics, almost all developed analytic theories assume a Maxwellian distribution function (MDF) and in some cases small deviations are described using the perturbation theory. The deviations with respect to the Maxwellian equilibrium, called kinetic effects, are required to be taken into account especially for fusion reactor plasmas. Generally, because the perturbation theory is not consistent with observed steady-state non-Maxwellians, these kinetic effects are numerically evaluated by very central processing unit (CPU)-expensive codes, avoiding the analytic complexity of velocity phase space integrals. We develop here a new method based on analytic non-Maxwellian distribution functions constructed from non-orthogonal basis sets in order to (i) use as few parameters as possible, (ii) increase the efficiency to model numerical and experimental non-Maxwellians, (iii) help to understand unsolved problems such as diagnostics discrepancies from the physical interpretation of the parameters, and (iv) obtain analytic corrections due to kinetic effects given by a small number of terms and removing the numerical error of the evaluation of velocity phase space integrals. This work does not attempt to derive new physical effects even if it could be possible to discover one from the better understandings of some unsolved problems, but here we focus on the analytic prediction of kinetic corrections from analytic non-Maxwellians. As applications, examples of analytic kinetic corrections are shown for the secondary electron emission, the Langmuir probe characteristic curve, and the entropy. This is done by using three analytic representations of the distribution function: the Kappa distribution function, the bi-modal or a new interpreted non-Maxwellian distribution function (INMDF). The existence of INMDFs is proved by new understandings of the experimental discrepancy of the measured electron temperature between two diagnostics in JET. As main results, it is shown that (i) the empirical formula for the secondary electron emission is not consistent with a MDF due to the presence of super-thermal particles, (ii) the super-thermal particles can replace a diffusion parameter in the Langmuir probe current formula, and (iii) the entropy can explicitly decrease in presence of sources only for the introduced INMDF without violating the second law of thermodynamics. Moreover, the first order entropy of an infinite number of super-thermal tails stays the same as the entropy of a MDF. The latter demystifies the Maxwell's demon by statistically describing non-isolated systems.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A method of revitalizing Sampson's theory of the Galilean satellites

A method is developed by which Sampson's theory of motion of Jupiter's Galilean satellites can be revitalized by use of algebraic manipulation software on a digital computer. The technique seeks (1) to remove algebraic errors existing in the current Sampson theory, (2) to introduce some neglected effects due to solar interactions and the 3-7 commensurability between the outer two satellites, (3) to allow for nonzero amplitude and phase of the libration, (4) to allow future revision of the arbitrary constants of integration, (5) to express the final results as analytic functions of variations in the numerous arbitrary constants of integration and arbitrary parameters, and (6) to provide analytic partial derivatives by means of which the numerical values of coefficients in the expressions for the coordinates can be adjusted. The level of precision desired is one arc second (Jovicentric) for the coordinates (2, 3.3, 5.2, and 9.1 km, respectively, for satellites I through IV).

Lieske, J. H.↗

Internal conversion and intersystem crossing dynamics based on coupled potential energy surfaces with full geometry-dependent spin–orbit and derivative couplings. Nonadiabatic photodissociation dynamics of NH 3 (A) leading to the NH(X 3 Σ – , a 1 Δ) + H 2 channel

Here we simulate the photodissociation of NH 3 originating from its first excited singlet state S 1 into the NH 2 + H (radical) and NH + H 2 (molecular) channels. The states considered are the ground singlet state S 0 , the first excited singlet state S 1 and the lowest-lying triplet state T 1 , which permit for the first time a uniform treatment of the internal conversion and intersystem crossing. The simulations are based on a diabatic potential energy matrix (DPEM) of S 0 , S 1 coupled by a conical intersection seam, as well as a potential energy surface (PES) for T 1 coupled by spin-orbit coupling (SOC) to the two singlet states. The DPEM and PES are fitted to ab initio electronic structure data (ESD) including energies, energy gradients, and derivative couplings. The DPEM also defines an adiabatic to diabatic state (AtD) transformation, which is used to transform the singular adiabatic SOC into a smooth function of the nuclear coordinates in the diabatic representation, allowing the diabatic SOC to be fit to an analytical functional form. ESD and SOC data obtained from these surfaces can serve as input for either quantum or semi-classical characterization of the nonadiabatic dynamics. Using the SHARC suite of programs, nonadiabatic simulations based on over 40 000 semi-classical trajectories assess the convergence of our results. The production of NH + H 2 is not direct, but is only achieved through a quasi-statistical dissociation mechanism after internal conversion to the ground electronic state. This leads to a much lower yield comparing with the main NH 2 + H channel. The NH(X 3 Σ_) radical produced through the intersystem crossing from S 0 to T 1 is rare (~0.2%) compared to NH(a 1 Δ) due to the process being spin forbidden.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Ginzburg--Landau functionals in the large-graph limit

Ginzburg–Landau (GL) functionals on graphs, which are relaxations of graph-cut functionals on graphs, have yielded a variety of insights in image segmentation and graph clustering. In this paper, we study large-graph limits of GL functionals by taking a functional-analytic view of graphs as nonlocal kernels. For a graph Wn with n nodes, the corresponding graph GL functional GL W n ϵ is an energy for functions on Wn. We minimize GL functionals on sequences of growing graphs that converge to functions called graphons. For such sequences of graphs, we show that the graph GL functional Γ-converges to a continuous and nonlocal functional that we call the graphon GL functional. We investigate the sharp-interface limits of the graph GL and graphon GL functionals, and we relate these limits to a nonlocal total-variation (TV) functional. We express the limiting GL functional in terms of Young measures and thereby obtain a probabilistic interpretation of the minimization problem in the large-graph limit. Finally, to develop intuition about graphon GL functionals, we determine the GL minimizer for several example families of graphons.

Zhang, Edith↗

A Simplified Model for an Evolving Protoplanetary Nebula

The dynamical evolution of the protoplanetary nebula is investigated using analytical solutions of the surface density transport equations. Constant and beta viscosity turbulence models are compared with a functional analytical model and the well-known alpha viscosity formulation. The beta viscosity model, heretofore used for steady-state disks, is shown to be a viable tool for separating dynamic and thermodynamic properties of an evolving disk.

Davis, Sanford S.↗

Enabling Computational Nanotechnology through JavaGenes in a Cycle Scavenging Environment

A genetic algorithm procedure is developed and implemented for fitting parameters for many-body inter-atomic force field functions for simulating nanotechnology atomistic applications using portable Java on cycle-scavenged heterogeneous workstations. Given a physics based analytic functional form for the force field, correlated parameters in a multi-dimensional environment are typically chosen to fit properties given either by experiments and/or by higher accuracy quantum mechanical simulations. The implementation automates this tedious procedure using an evolutionary computing algorithm operating on hundreds of cycle-scavenged computers. As a proof of concept, we demonstrate the procedure for evaluating the Stillinger-Weber (S-W) potential by (a) reproducing the published parameters for Si using S-W energies in the fitness function, and (b) evolving a "new" set of parameters using semi-empirical tightbinding energies in the fitness function. The "new" parameters are significantly better suited for Si cluster energies and forces as compared to even the published S-W potential.

Globus, Al↗

Educating local radial basis functions using the highest gradient of interest in three dimensional geometries: EDUCATING LOCALIZED RBF-SCHEMES

We present a novel methodology to effectively localize radial basis function approximation methods in three dimensions. The local scheme requires shape parameter-dependent functions that can be used to approximate gradients of scattered data and to solve partial differential operators. The optimum shape parameter is obtained from the highest gradient of interest, where a known analytical function, when boundary conditions are not present, or a shape parameter-free global approximation are used to educate the localized scheme. The later option is applicable to problems where the operator needs to be solved multiple times, like in time evolution or stochastic integration. Past shape parameter's optimizations, for two-dimensional domains, based on the condition number of the interpolant matrix, were unable to provide satisfactory approximations. The applicability of our method is illustrated in the context of an analytical expression interpolation and during a Ginzburg-Landau relaxation of a free energy functional. In general, the optimum shape parameter depends on geometry, node distribution, and density, whereas the approximation errors decrease as the node density and the local stencil size increase. The effective localization of radial basis functions motivates its use in moving boundary problems and accelerates solutions through sparse matrix solvers.

free energy functional↗

Numerical Optimization of Density Functional Tight Binding Models: Application to Molecules Containing Carbon, Hydrogen, Nitrogen, and Oxygen

New parametrizations for semiempirical density functional tight binding (DFTB) theory have been developed by the numerical optimization of adjustable parameters to minimize errors in the atomization energy and interatomic forces with respect to ab initio calculated data. Initial guesses for the radial dependences of the Slater- Koster bond integrals and overlap integrals were obtained from minimum basis density functional theory calculations. The radial dependences of the pair potentials and the bond and overlap integrals were represented by simple analytic functions. The adjustable parameters in these functions were optimized by simulated annealing and steepest descent algorithms to minimize the value of an objective function that quantifies the error between the DFTB model and ab initio calculated data. The accuracy and transferability of the resulting DFTB models for the C, H, N, and O system were assessed by comparing the predicted atomization energies and equilibrium molecular geometries of small molecules that were not included in the training data from DFTB to ab initio data. The DFTB models provide accurate predictions of the properties of hydrocarbons and more complex molecules containing C, H, N, and O.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Parameterization of electron attachment rate constants for impurities in LArTPC detectors

The ability of free electrons to drift long distances at high velocities in pure liquid argon under an applied electric field has been exploited for the past forty years to implement detectors with increasingly larger volumes for high energy physics research. The attachment of free electrons to impurities in the LAr is an important limit on the free instrumented volume of these extremely large detectors, and impurity concentrations as small as 100 ppt can reduce their resolution and efficiency. In this paper, we summarize the electron attachment rate constants as a function of the applied electric field, for common impurities in LArTPCs, obtained from data in the literature. Here, we further provide analytical functions to parameterize the data, which are useful to compare with new measurements, to model and analyze the performance of existing detectors, and to predict the performance of new detectors.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

A functional form for liquid scintillator pulse shapes

Digitization of detector signals enables analysis of the original waveform to extract timing, particle identification, and energy deposition information. Here we present the use of analytical functions based on sigmoids to model and fit such pulse shapes from liquid organic scintillators, though the method should also be applicable to other detector systems. Neutron and gamma interactions in NE213 detectors were digitized from the phototube anode and fit using a sigmoid-based function. The acuity of the fit in extracting timing information and performing neutron-gamma pulse-shape discrimination are presented and discussed.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Unsteady supersonic aerodynamic theory for interfering surfaces by the method of potential gradient

A generalized solution of the hyperbolic wave equation was further developed to relate the velocity components at a field point to the potential gradient distribution in the dependence domain. Singular integrals were evaluated in closed form, with numerical integration methods for more complex but analytic functions. Idealization of the lifting surfaces by trapezoidal elements with two sides parallel to the streamlines is computationally efficient. Streamwise integrals were performed analytically, and spanwise integrals were neccessary only on element leading and trailing sides. All integrands vanish on the Mach cone. Pressure distribution on a double delta wing and generalized aerodynamic coefficients for three AGARD planforms were calculated and compared with available results.

Jones, W. P.↗

Parameterization of photon beam dosimetry for a linear accelerator

Purpose: In radiation therapy, accurate data acquisition of photon beam dosimetric quantities is important for (1) beam modeling data input into a treatment planning system (TPS), (2) comparing measured and TPS modeled data, (3) the quality assurance process of a linear accelerator’s (Linac) beam characteristics, (4) the establishment of a standard data set for comparison with other data, etcetera. Parameterization of the photon beam dosimetry creates a data set that is portable and easy to implement for different applications such as those previously mentioned. The aim of this study is to develop methods to parameterize photon beam dosimetric quantities, including percentage depth doses (PDDs), profiles, and total scatter output factors (S{sub cp}). Methods: S{sub cp}, PDDs, and profiles for different field sizes, depths, and energies were measured for a Linac using a cylindrical 3D water scanning system. All data were smoothed for the analysis and profile data were also centered, symmetrized, and geometrically scaled. The S{sub cp} data were analyzed using an exponential function. The inverse square factor was removed from the PDD data before modeling and the data were subsequently analyzed using exponential functions. For profile modeling, one halfside of the profile was divided into three regions described by exponential, sigmoid, and Gaussian equations. All of the analytical functions are field size, energy, depth, and, in the case of profiles, scan direction specific. The model’s parameters were determined using the minimal amount of measured data necessary. The model’s accuracy was evaluated via the calculation of absolute differences between the measured (processed) and calculated data in low gradient regions and distance-to-agreement analysis in high gradient regions. Finally, the results of dosimetric quantities obtained by the fitted models for a different machine were also assessed. Results: All of the differences in the PDDs’ buildup and the profiles’ penumbra regions were less than 2 and 0.5 mm, respectively. The differences in the low gradient regions were 0.20% ± 0.20% (<1% for all) and 0.50% ± 0.35% (<1% for all) for PDDs and profiles, respectively. For S{sub cp} data, all of the absolute differences were less than 0.5%. Conclusions: This novel analytical model with minimum measurement requirements was proved to accurately calculate PDDs, profiles, and S{sub cp} for different field sizes, depths, and energies.

61 RADIATION PROTECTION AND DOSIMETRY↗