Search NASA⌕ Search

SEARCH · Search NASA

Results for “Greens 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 811 records · Page 45

Lepton scattering from Ar 40 and Ti 48 in the quasielastic peak region

Neutron and proton spectral functions of 40 Ar, 40 Ca, and 48 Ti isotopes are computed using the ab initio self-consistent Green's function approach. The resulting radii and charge distributions are in good agreement with available experimental data. The spectral functions of Ar and Ti are then utilized to calculate inclusive (e,e') cross sections within a factorization scheme and are found to correctly reproduce the recent Jefferson Lab measurements. Based on these successful agreements, the weak charged and neutral current double-differential cross sections for neutrino- 40 Ar scattering are predicted in the quasielastic region. Finally, results obtained by replacing the (experimentally inaccessible) neutron spectral distribution of 40 Ar with the (experimentally accessible) proton distribution of 48 Ti are compared and the accuracy of this approximation is assessed.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Electromagnetic and scalar diffraction by a right-angled wedge with a uniform surface impedance

The diffraction of an electromagnetic wave by a perfectly-conducting right-angled wedge with one surface covered by a dielectric slab or absorber is considered. The effect of the coated surface is approximated by a uniform surface impedance. The solution of the normally incident electromagnetic problem is facilitated by introducing two scalar fields which satisfy a mixed boundary condition on one surface of the wedge and a Neumann of Dirichlet boundary condition on the other. A functional transformation is employed to simplify the boundary conditions so that eigenfunction expansions can be obtained for the resulting Green's functions. The eigenfunction expansions are transformed into the integral representations which then are evaluated asymptotically by the modified Pauli-Clemmow method of steepest descent. A far zone approximation is made to obtain the scattered field from which the diffraction coefficient is found for scalar plane, cylindrical or sperical wave incident on the edge. With the introduction of a ray-fixed coordinate system, the dyadic diffraction coefficient for plane or cylindrical EM waves normally indicent on the edge is reduced to the sum of two dyads which can be written alternatively as a 2 X 2 diagonal matrix.

Hwang, Y. M.↗

The control and estimation of large space structures

A shape control or estimation problem for a large space structure can be modeled by a partial differential equation which represents changes in shape with respect to space and time, together with spatially discrete forcing functions or observations which represent the placement of actuators or sensors at discrete points along the structure. The use of Green's functions to convert boundary value problems into integral equations provides a convenient treatment of this mixture of continuous and discrete mathematics. Control and estimation algorithms are developed for the one-dimensional static beam to illustrate this technique.

Weeks, C.↗

Resonance frequency shift of an acoustic chamber containing a rigid sphere

The resonance frequency shift of an acoustic rectangular chamber due to the presence of a rigid sphere has been measured for l = 1,2 modes as a function of sphere size and position. The frequency shift is the results of volume exclusion and wave scattering. An analytical Green's function calculation was used to explain the data, providing excellent agreement between the measured and the calculated values. Also reported are similar measurements for a thin disk and the ratio of second harmonic to fundamental pressure as a function of sphere position. The measurement shows that the sphere reduces the first harmonic content, with sharply peaked suppression minima at specific sphere positions.

Leung, E.↗

The effects of earth model uncertainty on the inversion of seismic data for seismic source functions

SUMMARY We use Monte Carlo simulations to explore the effects of earth model uncertainty on the estimation of the seismic source time functions that correspond to the six independent components of the point source seismic moment tensor. Specifically, we invert synthetic data using Green’s functions estimated from a suite of earth models that contain stochastic density and seismic wave-speed heterogeneities. We find that the primary effect of earth model uncertainty on the data is that the amplitude of the first-arriving seismic energy is reduced, and that this amplitude reduction is proportional to the magnitude of the stochastic heterogeneities. Also, we find that the amplitude of the estimated seismic source functions can be under- or overestimated, depending on the stochastic earth model used to create the data. This effect is totally unpredictable, meaning that uncertainty in the earth model can lead to unpredictable biases in the amplitude of the estimated seismic source functions.

58 GEOSCIENCES↗

$S$ Hmax orientation in the Alpine region from observations of stress-induced anisotropy of nonlinear elasticity

The orientation of $S$ Hmax is commonly estimated from in situ borehole breakouts and earthquake focal mechanisms. Borehole measurements are expensive, and therefore sparse, and earthquake measurements can only be made in regions with many well-characterized earthquakes. Here, we derive the stress-field orientation using stress-induced anisotropy in nonlinear elasticity. In this method, we measure the strain derivative of velocity as a function of azimuth. We use a natural pump-probe (NPP) approach which consists of measuring elastic wave speed using empirical Green’s functions (probe) at different points of the earth tidal strain cycle (pump). The approach is validated using a larger data set in the Northern Alpine Foreland region where the orientation of maximum horizontal compressive stress is known from borehole breakouts and drilling-induced fractures. The technique resolves NNW-SSW to N-S directed $S$ Hmax which is in good agreement with conventional methods and the recent crustal stress model. We confirm that the NPP method can be applied to dense large-scale seismic arrays. The technique is then applied to the Southern Alps to understand the contemporary stress pattern associated with the ongoing deformation due to counterclockwise rotation of the Adriatic plate with respect to the European plate. Our results explain why the two major faults in Northeastern Italy, the Giudicarie Fault and the Periadriatic Line (Pustertal–Gailtal Fault) are currently inactive, while the currently acting stress field allows faults in Slovenia to deform actively. We have demonstrated that the pump-probe method has the potential to fill in the measurement gap left by conventional approaches, both in terms of regional coverage and in depth.

58 GEOSCIENCES↗

Examination of the sensitivity of quasifree reactions to details of the bound-state overlap functions

It is often stated that heavy-ion nucleon knockout reactions are mostly sensitive to the tails of the bound-state wave functions. In contrast, (p,2p) and (p,pn) reactions are known to access information on the full overlap functions within the nucleus. We analyze the oxygen isotopic chain and explore the differences between single-particle wave functions generated with potential models, used in the experimental analysis of knockout reactions, and ab initio computations from self-consistent Green's function theory. Contrary to common belief, we find that not only the tail of the overlap functions, but also their internal part is assessed in both reaction mechanisms, which are crucial to yield accurately determined spectroscopic information.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

MMIC devices for active phased array antennas

The study of printed circuit discontinuities is necessary in order to design, for example, transitions between rectangular waveguides and printed circuits. New developments with respect to the analytical approaches to this problem are discussed. A summary of the progress in the experimental approach is presented. The accurate solution for the modes in various millimeter-wave waveguides is essential in the analysis of many integrated circuit components, such as filters and impedance transformers. Problems associated with the numerical computation of these modes in two frequently used waveguide forms, namely, the finline and microstrip, are presented. The spectral domain method of formulation, with a moment method solution, is considered. This approach can be readily extended to analyze an arbitrary configuration of dielectric and metallized regions in a shielded enclosure. Galerkin's method is used, where the testing and basic functions are the same. It is shown that the mode functions, or eigenfunctions, are more sensitive to errors than the phase constants, or eigenvalues. The approximate mode functions do not satisfy the orthogonality relationship well, resulting in difficulties when these modal solutions are used to form an approximate Green's function or are used in a mode matching analysis.

Mittra, R.↗

Phonon scattering limited performance of monolayer MoS{sub 2} and WSe{sub 2} n-MOSFET

In this paper we show the effect of electron-phonon scattering on the performance of monolayer (1L) MoS{sub 2} and WSe{sub 2} channel based n-MOSFETs. Electronic properties of the channel materials are evaluated using the local density approximation (LDA) in density functional theory (DFT). For phonon dispersion we employ the small displacement / frozen phonon calculations in DFT. Thereafter using the non-equilibrium Green’s function (NEGF) formalism, we study the effect of electron-phonon scattering and the contribution of various phonon modes on the performance of such devices. It is found that the performance of the WSe{sub 2} device is less impacted by phonon scattering, showing a ballisticity of 83% for 1L-WSe{sub 2} FET for channel length of 10 nm. Though 1L-MoS{sub 2} FET of similar dimension shows a lesser ballisticity of 75%. Also in the presence of scattering there exist a a 21–36% increase in the intrinsic delay time (τ) and a 10–18% reduction in peak transconductance (g{sub m}) for WSe{sub 2} and MoS{sub 2} devices respectively.

36 MATERIALS SCIENCE↗

Topological insulator-metal transition and molecular electronics device based on zigzag phagraphene nanoribbon

Here, we explore the electronic transport properties of a graphene allotrope composed of 5–6-7 carbon aromatic rings called phagraphene and compare with the results of the transition-voltage spectroscopy (TVS) and propose the behavior at low voltage characteristic of a topological insulator. Phagraphene properties were compared to those of graphene in a zigzag nanoribbon configuration, zigzag graphene vs zigzag phagraphene nanoribbon (zzGNR and zzPGNR). The molecular geometry and the electronic properties were calculated by density functional theory (DFT) without spin, and the electronic transport and TVS were obtained by means of DFT combined with non-equilibrium Green´s function when we couple the optimized geometry of zzGNR and zzPGNR to the leads (left and right), forming the molecular junction that will be subjected to the action of an external bias voltage (V e ) to generate the molecular device. The findings exhibit (i) a metal-insulator transition when V e is increased until V e = 1.4 V which corresponds to the nonlinear region (resonance), showing the field effect transistor behaviour for zzGNR junctions; and (ii) two nonlinear regions (two negative differential resistances), showing a resonant tunnel diode behaviour with two operation windows (V e = 0.5 V and V e = 1.7 V) for the zzPGNR junction. In addition, the zzPGNR junction exhibits topological insulator characteristics upon introducing topological defects such as pentagons and heptagons in the hexagonal lattice of graphene, and when V e = 1.7 V, there occurs a topological insulator-metal transition that can be seen in the behaviour of the density of states, transmittance, and frontier molecular orbitals with V e .

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Experimental and theoretical study of Rayleigh-Lamb waves in a plate containing a surface-breaking crack

Scattering of Rayleigh-Lamb waves by a normal surface-breaking crack in a plate has been studied both theoretically and experimentally. The two-dimensionality of the far field, generated by a ball impact source, is exploited to characterize the source function using a direct integration technique. The scattering of waves generated by this impact source by the crack is subsequently solved by employing a Green's function integral expression for the scattered field coupled with a finite element representation of the near field. It is shown that theoretical results of plate response, both in frequency and time, are similar to those obtained experimentally. Additionally, implication for practical applications are discussed.

Paffenholz, Joseph↗

Equation of motion coupled-cluster cumulant approach for intrinsic losses in x-ray spectra

We present an equation of motion coupled cluster approach for calculating and understanding intrinsic inelastic losses in core level x-ray absorption spectra (XAS). The method is based on a factorization of the transition amplitude in the time-domain, which leads to a convolution of an effective one-body spectrum and the core-hole spectral function. The spectral function characterizes these losses in terms of shake-up excitations and satellites, and is calculated using a cumulant representation of the core-hole Green’s function that includes non-linear corrections. The one-body spectrum also includes orthogonality corrections that enhance the XAS at the edge.

Rehr, John J.↗

Equation of state of boron nitride combining computation, modeling, and experiment

The equation of state (EOS) of materials at warm dense conditions poses significant challenges to both theory and experiment. We report a combined computational, modeling, and experimental investigation leveraging new theoretical and experimental capabilities to investigate warm-dense boron nitride (BN). The simulation methodologies include path integral Monte Carlo (PIMC), several density functional theory (DFT) molecular dynamics methods [plane-wave pseudopotential, Fermi operator expansion (FOE), and spectral quadrature (SQ)], activity expansion (actex), and all-electron Green's function Korringa-Kohn-Rostoker (mecca), and compute the pressure and internal energy of BN over a broad range of densities and temperatures. Our experiments were conducted at the Omega laser facility and the Hugoniot response of BN to unprecedented pressures (1200–2650 GPa). The EOSs computed using different methods cross validate one another in the warm-dense matter regime, and the experimental Hugoniot data are in good agreement with our theoretical predictions. By comparing the EOS results from different methods, we assess that the largest discrepancies between theoretical predictions are ≲ 4% in pressure and ≲ 3% in energy and occur at 10 6 K, slightly below the peak compression that corresponds to the K -shell ionization regime. At these conditions, we find remarkable consistency between the EOS from DFT calculations performed on different platforms and using different exchange-correlation functionals and those from PIMC using free-particle nodes. This provides strong evidence for the accuracy of both PIMC and DFT in the high-pressure, high-temperature regime. Moreover, the recently developed SQ and FOE methods produce EOS data that have significantly smaller statistical error bars than PIMC, and so represent significant advances for efficient computation at high temperatures. The shock Hugoniot predicted by PIMC, actex, and mecca shows a maximum compression ratio of 4.55 ± 0.05 for an initial density of 2.26 g/cm 3 , higher than the Thomas-Fermi predictions by about 5%. In addition, we construct tabular EOS models that are consistent with the first-principles simulations and the experimental data. Our findings clarify the ionic and electronic structure of BN over a broad range of temperatures and densities and quantify their roles in the EOS and properties of this material. The tabular models may be utilized for future simulations of laser-driven experiments that include BN as a candidate ablator material.

36 MATERIALS SCIENCE↗

Remote sensing of earth terrain

In remote sensing, the encountered geophysical media such as agricultural canopy, forest, snow, or ice are inhomogeneous and contain scatters in a random manner. Furthermore, weather conditions such as fog, mist, or snow cover can intervene the electromagnetic observation of the remotely sensed media. In the modelling of such media accounting for the weather effects, a multi-layer random medium model has been developed. The scattering effects of the random media are described by three-dimensional correlation functions with variances and correlation lengths corresponding to the fluctuation strengths and the physical geometry of the inhomogeneities, respectively. With proper consideration of the dyadic Green's function and its singularities, the strong fluctuation theory is used to calculate the effective permittivities which account for the modification of the wave speed and attenuation in the presence of the scatters. The distorted Born approximation is then applied to obtain the correlations of the scattered fields. From the correlation of the scattered field, calculated is the complete set of scattering coefficients for polarimetric radar observation or brightness temperature in passive radiometer applications. In the remote sensing of terrestrial ecosystems, the development of microwave remote sensing technology and the potential of SAR to measure vegetation structure and biomass have increased effort to conduct experimental and theoretical researches on the interactions between microwave and vegetation canopies. The overall objective is to develop inversion algorithms to retrieve biophysical parameters from radar data. In this perspective, theoretical models and experimental data are methodically interconnected in the following manner: Due to the complexity of the interactions involved, all theoretical models have limited domains of validity; the proposed solution is to use theoretical models, which is validated by experiments, to establish the region in which the radar response is most sensitive to the parameters of interest; theoretically simulated data will be used to generate simple invertible models over the region. For applications to the remote sensing of sea ice, the developed theoretical models need to be tested with experimental measurements. With measured ground truth such as ice thickness, temperature, salinity, and structure, input parameters to the theoretical models can be obtained to calculate the polarimetric scattering coefficients for radars or brightness temperature for radiometers and then compare theoretical results with experimental data. Validated models will play an important role in the interpretation and classification of ice in monitoring global ice cover from space borne remote sensors in the future. We present an inversion algorithm based on a recently developed inversion method referred to as the Renormalized Source-Type Integral Equation approach. The objective of this method is to overcome some of the limitations and difficulties of the iterative Born technique. It recasts the inversion, which is nonlinear in nature, in terms of the solution of a set of linear equations; however, the final inversion equation is still nonlinear. The derived inversion equation is an exact equation which sums up the iterative Neuman (or Born) series in a closed form and, thus, is a valid representation even in the case when the Born series diverges; hence, the name Renormalized Source-Type Integral Equation Approach.

Yueh, Herng-Aung↗

Improving the Accuracy of Quadrature Method Solutions of Fredholm Integral Equations That Arise from Nonlinear Two-Point Boundary Value Problems

In this paper we are concerned with high-accuracy quadrature method solutions of nonlinear Fredholm integral equations of the form y(x) = r(x) + definite integral of g(x, t)F(t,y(t))dt with limits between 0 and 1,0 less than or equal to x les than or equal to 1, where the kernel function g(x,t) is continuous, but its partial derivatives have finite jump discontinuities across x = t. Such integral equations arise, e.g., when one applied Green's function techniques to nonlinear two-point boundary value problems of the form y "(x) =f(x,y(x)), 0 less than or equal to x less than or equal to 1, with y(0) = y(sub 0) and y(l) = y(sub l), or other linear boundary conditions. A quadrature method that is especially suitable and that has been employed for such equations is one based on the trepezoidal rule that has a low accuracy. By analyzing the corresponding Euler-Maclaurin expansion, we derive suitable correction terms that we add to the trapezoidal rule, thus obtaining new numerical quadrature formulas of arbitrarily high accuracy that we also use in defining quadrature methods for the integral equations above. We prove an existence and uniqueness theorem for the quadrature method solutions, and show that their accuracy is the same as that of the underlying quadrature formula. The solution of the nonlinear systems resulting from the quadrature methods is achieved through successive approximations whose convergence is also proved. The results are demonstrated with numerical examples.

Sidi, Avram↗

Improving the Accuracy of Quadrature Method Solutions of Fredholm Integral Equations that Arise from Nonlinear Two-Point Boundary Value Problems

In this paper we are concerned with high-accuracy quadrature method solutions of nonlinear Fredholm integral equations of the form y(x) = r(x) + integral(0 to 1) g(x,t) F(t, y(t)) dt, 0 less than or equal to x less than or equal to 1, where the kernel function g(x,t) is continuous, but its partial derivatives have finite jump discontinuities across x = t. Such integrals equations arise, e.g., when one applies Green's function techniques to nonlinear two-point boundary value problems of the form U''(x) = f(x,y(x)), 0 less than or equal to x less than or equal to 1, with y(0) = y(sub 0) and g(l) = y(sub 1), or other linear boundary conditions. A quadrature method that is especially suitable and that has been employed for such equations is one based on the trapezoidal rule that has a low accuracy. By analyzing the corresponding Euler-Maclaurin expansion, we derive suitable correction terms that we add to the trapezoidal thus obtaining new numerical quadrature formulas of arbitrarily high accuracy that we also use in defining quadrature methods for the integral equations above. We prove an existence and uniqueness theorem for the quadrature method solutions, and show that their accuracy is the same as that of the underlying quadrature formula. The solution of the nonlinear systems resulting from the quadrature methods is achieved through successive approximations whose convergence is also proved. The results are demonstrated with numerical examples.

Sidi, Avram↗

Analysis of Photothermal Characterization of Layered Materials: Design of Optimal Experiments

In this paper numerical calculations are presented for the steady-periodic temperature in layered materials and functionally-graded materials to simulate photothermal methods for the measurement of thermal properties. No laboratory experiments were performed. The temperature is found from a new Green s function formulation which is particularly well-suited to machine calculation. The simulation method is verified by comparison with literature data for a layered material. The method is applied to a class of two-component functionally-graded materials and results for temperature and sensitivity coefficients are presented. An optimality criterion, based on the sensitivity coefficients, is used for choosing what experimental conditions will be needed for photothermal measurements to determine the spatial distribution of thermal properties. This method for optimal experiment design is completely general and may be applied to any photothermal technique and to any functionally-graded material.

Cole, Kevin D.↗

Electronic structure of chromium trihalides beyond density functional theory

In this work, we explore the electronic band structure of freestanding monolayers of chromium trihalides Cr X 3 , X = Cl, Br, I, within an advanced ab initio theoretical approach based on the use of Green's function functionals. We compare the local density approximation with the quasiparticle self-consistent GW (QS GW ) approximation and its self-consistent extension ( QS G W ^ ) by solving the particle-hole ladder Bethe-Salpeter equations to improve the effective interaction W . We show that, at all levels of theory, the valence band consistently changes shape in the sequence Cl → Br → I , and the valence band maximum shifts from the M point to the Γ point. By analyzing the dynamic and momentum-dependent self-energy, we show that QS G W ^ adds to the localization of the systems in comparison with QS GW , thereby leading to a narrower band and reduced amount of halogens in the valence band manifold. Further analysis shows that X = Cl is most strongly correlated, and X = I is least correlated (most bandlike) as the hybridization between Cr d and X p enhances in the direction Cl → Br → I . For CrBr 3 and CrI 3 , we observe remarkable differences between the QS GW and QS G W ^ valence band structures, while their eigenfunctions are very similar. We show that weak perturbations, like moderate strain, weak changes to the d - p hybridization, and adding small U , can flip the valence band structures between these two solutions in these materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗