Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “sum rules”

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 415 records · Page 23

Noncentrosymmetric rare-earth copper gallium chalcogenides RE{sub 3}CuGaCh{sub 7} (RE=La–Nd; Ch=S, Se): An unexpected combination

The quaternary rare-earth chalcogenides RE{sub 3}CuGaS{sub 7} and RE{sub 3}CuGaSe{sub 7} (RE=La–Nd) have been prepared by reactions of the elements at 1050 °C and 900 °C, respectively. They crystallize in the noncentrosymmetric La{sub 3}CuSiS{sub 7}-type structure (hexagonal, space group P6{sub 3}, Z=2) in which the a-parameter is largely controlled by the RE component (a=10.0–10.3 Å for the sulfides and 10.3–10.6 Å for the selenides) whereas the c-parameter is essentially fixed by the choice of Ga and chalcogen atoms within tetrahedral units (c=6.1 Å for the sulfides and 6.4 Å for the selenides). They extend the series RE{sub 3}MGaCh{sub 7}, previously known for divalent metal atoms (M=Mn–Ni), differing in that the Cu atoms in RE{sub 3}CuGaCh{sub 7} occupy trigonal planar sites instead of octahedral sites. Among quaternary chalcogenides RE{sub 3}MM′Ch{sub 7}, the combination of monovalent (M=Cu) and trivalent (M′=Ga) metals is unusual because it appears to violate the condition of charge balance satisfied by most La{sub 3}CuSiS{sub 7}-type compounds. The possibility of divalent Cu atoms was ruled out by bond valence sum analysis, magnetic measurements, and X-ray photoelectron spectroscopy. The electron deficiency in RE{sub 3}CuGaCh{sub 7} is accommodated through S-based holes at the top of the valence band, as shown by band structure calculations on La{sub 3}CuGaS{sub 7}. An optical band gap of about 2.0 eV was found for La{sub 3}CuGaSe{sub 7}. - Graphical abstract: The chalcogenides RE{sub 3}CuGaCh{sub 7} contain monovalent Cu in trigonal planes and trivalent Ga in tetrahedra; they are electron-deficient representatives of La{sub 3}CuSiS{sub 7}-type compounds, which normally satisfy charge balance. - Highlights: • Quaternary chalcogenides RE{sub 3}CuGaCh{sub 7} (RE=La–Nd; Ch=S, Se) were prepared. • Bond valence sums, magnetism, and XPS data give evidence for monovalent Cu. • Crystal structures reveal high anisotropy of Cu displacement. • Electron deficiency is accommodated by S-based holes in valence band.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A quasi-vortex-lattice method in thin wing theory

A quasi-continuous method is developed for solving thin-wing problems. For the purpose of satisfying the wing boundary conditions, the spanwise vortex distribution is assumed to be stepwise-constant, while the chordwise vortex integral is reduced to a finite sum through a modified trapezoidal rule and the theory of Chebyshev polynomials. Wing-edge and Cauchy singularities are acounted for. The total aerodynamic characteristics are obtained by an appropriate quadrature integration. The two-dimensional results for airfoils without flap deflection reproduce the exact solutions in lift and pitching moment coefficients, the leading edge suction, and the pressure difference at a finite number of points. For a flapped airfoil, the present results are more accurate than those given by the vortex-lattice method. The three-dimensional results also show an improvement over the results of the vortex-lattice method. Extension to nonplanar applications is discussed.

Lan, C. E.↗

Various Ambiguities in Re-constructing Laser Pulse Parameters

We think that mode lock laser pulses are generated by the summation process that take place between the monochromatic EM filed frequencies as if they interact with each other as shown in equation 1. In reality, the pulse generation is a collaborative interaction process between EM fields and various material medium. When we carry out the actual mode lock analysis, we do take into account of interpaly between all the temporal dynamics of the cavity gain medium, cavity round trip time and the response time of the intra cavity element (saturable absorber, Kerr medium, etc.). that really enforces the locking of the phase of the cavity spontaneous emissions. On a conceptual level, this simplistic representation of the mode locking by Eq.1 ignores all these critical physical processes. When we try to analyze a pulsed field, again we start by representing it very much like this equation, even though we can only detect the square modulus of this complex field and loose a lot of phase related information to the detectors quantum whims and their time constants. The key parameters for a light pulse are as follows. Foremost is the (i) carrier frequency, which cannot be described or imagined without its state of undulation expressed as its (ii) phase. Next is our imagined time finite (iii) carrier envelope that provides the temporal boundary of the field amplitude strength of the undulating E-field. The final parameter is the (iv) state of polarization or the unique plane along which the strength of the E-field gradient undulates. None of these filed characteristics are made self-evident to us by the fields themselves. We do not see light. Light does not see light. Light beams pass through each other without altering each others energy distribution unless there are interacting material molecules (dipoles) within the physical volume of superposition of the beams. In contrast, we can sense the material particles. Material particles sense each other and they cannot pass through each other without interacting with (scattering from) each other. Thus the interpretation of the superposition phenomenon of multiple fields on detectors should not be lumped under the mysterious "wave-particle duality" philosophy. The phenomenon of superposition can be understood better when we focus on the actual process experienced by the detecting dipoles when allowed by QM rules, they respond to and sum all the induced stimulations due to all the superposed fields followed by the proportionate energy absorption giving rise to the fringes we observe. We will present various experimental results to illustrate our arguments. Our position is that such detector behavior driven interpretations rather than the generally implied field-field interaction driven explanations, will help us better understand the ultimate nature of light and hence invent better and newer devices and instruments.

Roychoudhuri, Chandrasekhar↗

Derivation of Effective Properties Based on Porous Scale Simulations Using Filtering Techniques

This study presents a method for derivation of effective properties at the interface and in-depth of porous materials. The method defines a Representative Elementary Volume (REV) and applies filtering techniques to computer effective properties such as porosity and flow quantities, such as velocity and pressure. The script, developed to process the data was tested on the VTK type files that contain the mesh information and the flow solution. The method allows to choose between two types of filters, such as cellular and top-hat and define the size of the REV and number of samples along the domain. Extraction of the REV from the domain is performed to exact boundaries requested for the user. This is done using a triangulation technique and cutting through the cells to comply to the requested boundaries of the volume. The method can be applied to both structured and unstructured meshes. Filtering the material porosity and flow quantities involves integration of the numerical data. The algorithm provides three integration methods, such as Riemann sum, Monte Carlo and Quadrature rule to perform the integration. The Monte-Carlo technique permits the use of either uniform or linearly spaced distribution of points. The Quadrature rule is currently applicable to tetrahedral element types. The Monte Carlo and Quadrature rule methods require interpolation of the flow quantities at the sample points. For interpolation, two methods were tested and are readily available, Gaussian interpolation and re-sampling. It has been shown that re-sampling method has better consistency and acceptable accuracy in interpolation of the data. The algorithm was written in Python language and uses a number of modules. The major module besides numpy is PyVista. It is used to process the computational domain, clip the REV and interpolate the data. Quadrature rule integration was performed using a quadpy module. ParaView software was used externally to convert the flow solution to the VTK (or more specifically VTU) format. Integration of ParaView in the same environment with PyVista encountered problems and could not be implemented in this work. The developed algorithm is expected to be applicable to unstructured meshes and more complex porous structures as soon as the data can be passed in VTK type format. With the report is provided Python script for filtering the solution and a Matlab script for simple generation and processing of 2-D and 3-D porous channel geometries. The two scripts don't communicate.

Alexsander Zibitsker↗

Optimum frame synchronization.

This paper considers the optimum method for locating a sync word periodically imbedded in binary data and received over the additive white Gaussian noise channel. It is shown that the optimum rule is to select the location that maximizes the sum of the correlation and a correction term. Simulations are reported that show approximately a 3-dB improvement at interesting signal-to-noise ratios compared to a pure correlation rule. Extensions are given to the phase-shift keyed sync case where the detector output has a binary ambiguity and to the case of Gaussian data.

Massey, J. L.↗

Dislocation Patterning in Deforming Crystals: Theory, Computational Predictions and Validation (Final Technical Report)

This project was awarded for an initial period of three years, followed by a no-cost extension for one year and a funded one-year renewal. During the first four years, the project focused on investigating the role of dislocation reactions and dislocation correlation in dislocation patterning in FCC metals. During the fifth year, the scope of research was expanded to investigate the effects of composition inhomogeneity on the mesoscale plastic response of Body-Centered Cubic (BCC) alloys. In addition to its impact on metal hardening during deformation, dislocation patterning provides the microstructure information required to understand phenomena like fracture and recrystallization in metals. The Continuum Dislocation Dynamics (CDD) framework was the methodology used during the first four years, followed by the use of Discrete Dislocation Dynamics (DDD) during the fifth year. CDD is a density-based formalism of dislocation dynamics in which the plastic deformation of crystals is predicted together with the mesoscale dislocation patterns by tracking the space-time evolution of dislocations as driven by the applied stress, short-range and long-range interactions of dislocations, and cross slip. CDD is a crystal mechanics approach in which the plastic constitutive part is replaced with the equations of dislocation dynamics, which is driven by the internal stress via a mobility laws, while the evolution of the density gives the eigenstrain required to update the internal stress itself. The governing equations are thus those of crystal mechanics cast as an eigenstrain problem and those of dislocation transport and reactions. Our investigations during the first four years were driven by the hypothesis that that dislocation patterning is triggered by spatiotemporal dislocation density and internal stress fluctuations, and that such fluctuations influence the collective dislocation dynamic through their effects on the short-range reaction rates and the collective dislocation mobility. This hypothesis was tested by modeling the influence of the dislocation reactions and dislocation correlations on collective dislocation dynamics. The main research components during the first four years were: 1) Reformulation of the CDD framework to integrate the dislocation correlations and dislocation reactions. 2) Simulation of the dislocation correlations and quantifying their contribution to the long-range stress of complex dislocation systems. 3) Computational implementation of the updated CDD model for FCC single crystals and development of an efficient CDD code. 4) Investigation of dislocation patterning in FCC crystals based on the updated CDD model. Some key findings from these investigations were: 1) The patterning of dislocations is initiated by cross slip and internal stress fluctuations, and the refinement of the pattern results from junction formation. 2) Patterning is more prominent under high stress. 3) The correlation stress was found to be a significant part of the mean field stress in continuum representations of dislocation dynamics. During the last year of the project, we have established a method for performing dislocation dynamics in inhomogeneous alloys and demonstrated the impact of composition inhomogeneity on the characteristics of initial yielding and dislocation character in BCC alloys. We also tested this method using an irradiated ferritic alloy in which the composition and dislocation loops effects are both present. In doing so, we have considered two aspects of the role of inhomogeneity, the creation of internal coherency stress due to the dependence of the lattice parameter on the local composition, and the dependence of the dislocation mobility on the local composition. These two aspects resulted in an unexpected behavior of the dislocation system, and, in turn, the yielding behavior of BCC alloys. Key findings include: 1) The composition inhomogeneity in BCC alloys alters the well-known role of screw dislocations by forcing a level of waviness on the average dislocation character, thus destroying the screw character dominance in the yielding observed in pure or homogeneous BCC metals. 2) In the case of irradiated alloys, while composition inhomogeneity by itself and dislocation loops by themselves result in some hardening, the superposition of the two mechanisms does not result in the sum of the two contributions via any known rule. The latter result contradicts the classical works related hardening via multiple mechanisms, which state that various hardening mechanisms can be added linearly or using some Pythagorean additive rule. We were able to rationalize this unexpected result by a closer loop at the origin of the hardening mechanisms themselves, which both involved long-range elastic interactions with dislocations. We reached the conclusion that the hardening resulting from the total stress associated with these two mechanisms (current work) differs from the sum of the effects of the individual stress fields (classical literature). We have also found a major flaw in the classical theory of spinodal hardening. The latter theory never considered the impact of composition undulation on the solute hardening. We have built a new general theory accounting for the effect of solute friction together with the composition undulations on the dislocation configuration at Critical Resolved Shear Stress (CRSS), which is yielding a new definition of the spinodal hardening and new different results. This work is the first to recognize the lack of the role of solute in spinodal hardening theory of alloys.

36 MATERIALS SCIENCE↗

On Some New Properties of Quaternion Functions

Quaternions discovered by W. R. Hamilton made a great contribution to the progress in noncommutative algebra and vector analysis. However, the analysis of quaternion functions has not been duly developed. The matter is that the notion of a derivative of quaternion functions of a quaternion variable has not been known until recently. The author has succeeded in improving the situation. The present work contains an account of the results obtained by him in this direction. The notion of an ℍ-derivative is introduced for quaternion functions of a quaternion variable. The existence of an ℍ-derivative of elementary functions is established retaining the well-known formulas for the corresponding functions from complex (real) analysis. The rules on the ℍ-differentiation of a sum, a product, and an inverse function are formulated and proved. Necessary and sufficient conditions for the existence of an ℍ-derivative are established. The notions of ℂ{sup 2}-differentiation and ℂ{sup 2}-holomorphy are introduced for quaternion functions of a quaternion variable. Three equivalent conditions are found, each of them being a necessary and sufficient one for ℂ{sup 2}-differentiation. Representations by an integral and a power series are given for ℂ{sup 2}-holomorphic functions. It is proved that the harmonicity of functions f(z), z · f(z), and f(z) · z is the necessary and sufficient condition for a function f to be Fueter-regular.

97 MATHEMATICS AND COMPUTING↗

Exploiting Kronecker structure in exponential integrators: Fast approximation of the action of φ $-$functions of matrices via quadrature

Here, in this article, we propose an algorithm for approximating the action of φ $-$ functions of matrices against vectors, which is a key operation in exponential time integrators. In particular, we consider matrices with Kronecker sum structure, which arise from problems admitting a tensor product representation. The method is based on quadrature approximations of the integral form of the φ $-$ functions combined with a scaling and modified squaring method. Owing to the Kronecker sum representation, only actions of 1D matrix exponentials are needed at each quadrature node and assembly of the full matrix can be avoided. Additionally, we derive a priori bounds for the quadrature error, which show that, as expected by classical theory, the rate of convergence of our method is supergeometric. Guided by our analysis, we construct a fast and robust method for estimating the optimal scaling factor and number of quadrature nodes that minimizes the total cost for a prescribed error tolerance. We investigate the performance of our algorithm by solving several linear and semilinear time-dependent problems in 2D and 3D. The results show that our method is accurate and orders of magnitude faster than the current state-of-the-art.

97 MATHEMATICS AND COMPUTING↗

Cornering (quasi) degenerate neutrinos with cosmology

In light of the improved sensitivities of cosmological observations, we examine the status of quasi-degenerate neutrino mass scenarios. Within the simplest extension of the standard cosmological model with massive neutrinos, we find that quasi-degenerate neutrinos are severely constrained by present cosmological data and neutrino oscillation experiments. We find that Planck 2018 observations of cosmic microwave background (CMB) anisotropies disfavour quasi-degenerate neutrino masses at 2.4 Gaussian >σ’s, while adding baryon acoustic oscillations (BAO) data brings the rejection to 5.9 σ’s. The highest statistical significance with which one would be able to rule out quasi-degeneracy would arise if the sum of neutrino masses is Σm v = 60 meV (the minimum allowed by neutrino oscillation experiments); indeed a sensitivity of 15 meV, as expected from a combination of future cosmological probes, would further improve the rejection level up to 17σ. We discuss the robustness of these projections with respect to assumptions on the underlying cosmological model, and also compare them with bounds from β decay endpoint and neutrinoless double beta decay studies.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

New Dimensions in the Theory of Excited States and X-ray Spectra (Final Report)

This Final Technical Report briefly summarizes the achievements during the lifetime of our DOE BES grant DE-FG02-97ER45623. The long-term goal of this project has been the development of quantitative theories of the interaction between radiation and matter, with a focus on x-ray spectroscopies. X-ray spectra have long been among the most important probes of atomic-scale structure and properties of matter, ranging from atoms and molecular systems to condensed matter and exotic states. These spectroscopies are widely used in investigations at the major DOE synchrotron x-ray facilities and related centers world-wide. In addition to fundamental theory, a major goal of our project has been the development of computational software that implements the theory for calculations of x-ray spectra of various materials throughout the periodic table. Due to the complex nature of x-ray spectra, quantitative theory is essential for its interpretation. The theory is challenging since it involves excited state electronic structure and many-body correlation effects that go beyond independent particle approximations like DFT or Hartree-Fock. Moreover, the experimental investigations typically involve a broad range of energy, time, and temperature scales, from the UV-Vis to hard x-ray energies of order 10 4 eV, and temperatures T from ambient up to the warm-dense-matter regime where the Fermi energy kBTF is of order a few eV, i.e., temperatures of order 105 Kelvin. This broad range of experimental conditions has fostered many novel theoretical approaches and computational techniques, many of which we have developed systematically over the duration of the grant. In contrast to the traditional wave-function approach of quantum theory and electronic structure methods, our theoretical approach is based on modern Green's function techniques. This approach is better suited for aperiodic structures, excited states, and broad spectral ranges, since it avoids the computational bottlenecks of sum-over-states approaches, as in the Fermi golden rule. This theoretical framework has been incorporated into efficient, user-friendly x-ray spectroscopy software which is now used routinely worldwide to simulate and analyze spectra. These theoretical tools provide an essential complement to synchrotron and next-generation light sources, which are used to investigate complex materials with ever increasing precision. Moreover, the synergism between theory, computation and experiment contributed by our research enhances scientific understanding and creates opportunities for innovations in materials and energy science and in many fields. As documented in this Report, this research grant has been remarkably successful in achieving these goals. In particular, this grant has supported the development of the x-ray spectroscopy software suite known as FEFF (named for an effective scattering amplitude f eff in the theory). The FEFF codes have become one of the premier tools for quantitative simulations of x-ray spectra as documented by many thousands of citations in the Web of Science and Google-Scholar.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Cardiac chamber volumes by echocardiography using a new mathematical method: A promising technique for zero-G use

Accurate, reproducible, and non-invasive means for ventricular volume determination are needed for evaluating cardiovascular function zero-gravity. Current echocardiographic methods, particularly for the right ventricle, suffer from a large standard error. A new mathematical approach, recently described by Watanabe et al., was tested on 1 normal formalin-fixed human hearts suspended in a mineral oil bath. Volumes are estimated from multiple two-dimensional echocardiographic views recorded from a single point at sequential angles. The product of sectional cavity area and center of mass for each view summed over the range of angles (using a trapezoidal rule) gives volume. Multiple (8-14) short axis right ventricle and left ventricle views at 5.0 deg intervals were videotaped. The images were digitized by two independent observers (leading-edge to leading-edge technique) and analyzed using a graphics tablet and microcomputer. Actual volumes were determined by filling the chambers with water. These data were compared to the mean of the two echo measurements.

Buckey, J. C.↗

Combined influences of gravitoinertial force level and visual field pitch on visually perceived eye level

Psychophysical measurements of the level at which observers set a small visual target so as to appear at eye level (VPEL) were made on 13 subjects in 1.0 g and 1.5 g environments in the Graybiel Laboratory rotating room while they viewed a pitched visual field or while in total darkness. The gravitoinertial force was parallel to the z-axis of the head and body during the measurements. The visual field consisted of two 58 degrees high, luminous, pitched-from-vertical, bilaterally symmetric, parallel lines, viewed in otherwise total darkness. The lines were horizontally separated by 53 degrees and presented at each of 7 angles of pitch ranging from 30 degrees with the top of the visual field turned away from the subject (top backward) to 30 degrees with the top turned toward the subject (top forward). At 1.5 g, VPEL changed linearly with the pitch of the 2-line stimulus and was depressed with top backward pitch and elevated with top forward pitch as had been reported previously at 1.0 g (1,2); however, the slopes of the VPEL-vs-pitch functions at 1.0 g and 1.5 g were indistinguishable. As reported previously also (3,4), the VPEL in darkness was considerably lower at 1.5 g than at 1.0 g; however, although the y-intercept of the VPEL-vs-pitch function in the presence of the 2-line visual field (visual field erect) was also lower at 1.5 g than at 1.0 g as it was in darkness, the G-related difference was significantly attenuated by the presence of the visual field. The quantitative characteristics of the results are consistent with a model in which VPEL is treated as a consequence of an algebraic weighted average or a vector sum of visual and nonvisual influences although the two combining rules lead to fits that are equally good.

Non-NASA Center↗

Unitarity cuts of the worldsheet

We compute the imaginary parts of genus-one string scattering amplitudes. Following Witten’s i\varepsilon i ε prescription for the integration contour on the moduli space of worldsheets, we give a general algorithm for computing unitarity cuts of the annulus, Möbius strip, and torus topologies exactly in \alpha' α ′ . With the help of tropical analysis, we show how the intricate pattern of thresholds (normal and anomalous) opening up arises from the worldsheet computation. The result is a manifestly-convergent representation of the imaginary parts of amplitudes, which has the analytic form expected from Cutkosky rules in field theory, but bypasses the need for performing laborious sums over the intermediate states. We use this representation to study various physical aspects of string amplitudes, including their behavior in the (s,t) ( s , t ) plane, exponential suppression, decay widths of massive strings, total cross section, and low-energy expansions. We find that planar annulus amplitudes exhibit a version of low-spin dominance: at any finite energy, only a finite number of low partial-wave spins give an appreciable contribution to the imaginary part.

97 MATHEMATICS AND COMPUTING↗

Electron collisions with methyl-substituted ethylenes: Cross section measurements and calculations for 2-methyl–2-butene and 2,3-dimethyl–2-butene

We report electron-scattering cross sections determined for 2-methyl–2-butene [(H{sub 3}C)HC = C(CH{sub 3}){sub 2}] and 2,3-dimethyl–2-butene [(H{sub 3}C){sub 2}C = C(CH{sub 3}){sub 2}] molecules. Absolute grand-total cross sections (TCSs) were measured for incident electron energies in the 0.5–300 eV range, using a linear electron-transmission technique. The experimental TCS energy dependences for the both targets appear to be very similar with respect to the shape. In each TCS curve, three features are discernible: the resonant-like structure located around 2.6–2.7 eV, the broad distinct enhancement peaking near 8.5 eV, and a weak hump in the vicinity of 24 eV. Theoretical integral elastic (ECS) and ionization (ICS) cross sections were computed up to 3 keV by means of the additivity rule (AR) approximation and the binary-encounter-Bethe method, respectively. Their sums, (ECS+ICS), are in a reasonable agreement with the respective measured TCSs. To examine the effect of methylation of hydrogen sides in the ethylene [H{sub 2}C = CH{sub 2}] molecule on the TCS, we compared the TCS energy curves for the sequence of methylated ethylenes: propene [H{sub 2}C = CH(CH{sub 3})], 2-methylpropene [H{sub 2}C = C(CH{sub 3}){sub 2}], 2-methyl–2-butene [(H{sub 3}C)HC = C(CH{sub 3}){sub 2}], and 2,3-dimethyl–2-butene [(H{sub 3}C){sub 2}C = C(CH{sub 3}){sub 2}], measured in the same laboratory. Moreover, the isomeric effect is also discussed for the C{sub 5}H{sub 10} and C{sub 6}H{sub 12} compounds.

2-METHYLPROPENE↗

Multiple symbol partially coherent detection of MPSK

It is shown that by using the known (or estimated) value of carrier tracking loop signal to noise ratio (SNR) in the decision metric, it is possible to improve the error probability performance of a partially coherent multiple phase-shift-keying (MPSK) system relative to that corresponding to the commonly used ideal coherent decision rule. Using a maximum-likeihood approach, an optimum decision metric is derived and shown to take the form of a weighted sum of the ideal coherent decision metric (i.e., correlation) and the noncoherent decision metric which is optimum for differential detection of MPSK. The performance of a receiver based on this optimum decision rule is derived and shown to provide continued improvement with increasing length of observation interval (data symbol sequence length). Unfortunately, increasing the observation length does not eliminate the error floor associated with the finite loop SNR. Nevertheless, in the limit of infinite observation length, the average error probability performance approaches the algebraic sum of the error floor and the performance of ideal coherent detection, i.e., at any error probability above the error floor, there is no degradation due to the partial coherence. It is shown that this limiting behavior is virtually achievable with practical size observation lengths. Furthermore, the performance is quite insensitive to mismatch between the estimate of loop SNR (e.g., obtained from measurement) fed to the decision metric and its true value. These results may be of use in low-cost Earth-orbiting or deep-space missions employing coded modulations.

Simon, M. K.↗

Approximation of periodic Green's operator in real space using numerical integration and its use in fast Fourier transform-based micromechanical models

In this work, we propose an expression for the periodic first derivative of Green's function in real space. The proposed expression allows an alternative way of computing the periodic Green's operator based on periodically summing the free-space Green's operator in terms of an appropriate quadrature rule. We provide computational examples, which show the accuracy of the proposed approach, together with reduced spurious oscillations in the solution fields.

42 ENGINEERING↗

Sparse low rank approximation of potential energy surfaces with applications in estimation of anharmonic zero point energies and frequencies

We introduce a method that exploits sparse representation of potential energy surfaces (PES) on a polynomial basis set selected by compressed sensing . The method is useful for studies involving large numbers of PES evaluations, such as the search for local minima, transition states, or integration. We apply this method for estimating zero point energies and frequencies of molecules using a three step approach. In the first step, we interpret the PES as a sparse tensor on polynomial basis and determine its entries by a compressed sensing based algorithm using only a few PES evaluations. Then, we establish a rank reduction strategy to compress this tensor in a suitable low-rank canonical tensor format using standard tensor compression tools. This allows representing a high dimensional PES as a small sum of products of one dimensional functions. Finally, a low dimensional Gauss–Hermite quadrature rule is used to integrate the product of sparse canonical low-rank representation of PES and Green’s function in the second-order diagrammatic vibrational many-body Green’s function theory (XVH2) for estimation of zero-point energies and frequencies. Numerical tests on molecules considered in this work indicate a more efficient scaling of computational cost with molecular size as compared to other methods.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Compressed sparse tensor based quadrature for vibrational quantum mechanics integrals

A new method for fast evaluation of high dimensional integrals arising in quantum mechanics is proposed. Here, the method is based on sparse approximation of a high dimensional function followed by a low-rank compression. In the first step, we interpret the high dimensional integrand as a tensor in a suitable tensor product space and determine its entries by a compressed sensing based algorithm using only a few function evaluations. Secondly, we implement a rank reduction strategy to compress this tensor in a suitable low-rank tensor format using standard tensor compression tools. This allows representing a high dimensional integrand function as a small sum of products of low dimensional functions. Finally, a low dimensional Gauss–Hermite quadrature rule is used to integrate this low-rank representation, thus alleviating the curse of dimensionality. Finally, numerical tests on synthetic functions, as well as on energy correction integrals for water and formaldehyde molecules demonstrate the efficiency of this method using very few function evaluations as compared to other integration strategies.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗