Search NASA⌕ Search

SEARCH · Search NASA

Results for “time Green 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 253 records · Page 14

Kirchhoff's Integral Representation and a Cavity Wake Potential

A method is proposed for the calculation of the short-range wake field potentials of an ultra-relativistic bunch passing near some irregularities in a beam pipe. The method is based on the space-time domain integration of Maxwell's equations using Kirchhoff's formulation. We demonstrate this method on two cases where we obtain the wake potentials for the energy loss of a bunch traversing an iris-collimator in a beam pipe and for a cavity. Likewise, formulas are derived for Green's functions that describe the transverse force action of wake fields. Simple formulas for the total energy loss of a bunch with a Gaussian charge density distribution are derived as well. The derived estimates are compared with computer results and predictions of other models.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Enhancement of thermoelectric performance in InAs nanotubes by tuning quantum confinement effect

By using the nonequilibrium Green's function method, we study the thermoelectric properties of InAs nanotubes. The results show that InAs nanotube with a certain internal diameter has much higher ZT value than nanowire due to the enhancement of quantum confinement effect leading to the increase of the power factor S{sup 2}G. The ZT value of InAs nanotube can reach 1.74, which is about three times greater than that of nanowires. Moreover, it is found that the ZT values of InAs nanotubes decrease rapidly with the increase of internal diameter, which results from the rapid increase of phonons thermal conductance due to the “red shift” of low-frequency optical phonon modes.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

Green’s function coupled cluster simulation of the near-valence ionizations of DNA-fragments

Accurate description of the ionization process in DNA is crucial to the understanding of the DNA damage under exposure to ionizing radiation and the exploration of the potential application of DNA strands in nanoelectronics. In this work, by employing our recently developed Green’s function coupled-cluster library on supercomputing facilities, we have studied the spectral functions of several guanine–cytosine (G–C) base pair structures ([G–C]n, n = 1–3) for the first time in a relatively broad near-valence regime ([−25.0, −5.0] eV) in the coupled-cluster with singles and doubles level. Our focus is to give a preliminary many-body coupled-cluster understanding and guideline of the vertical ionization energy (VIE), spectral profile, and ionization feature changes of these systems as the system size expands in this near-valence regime. The results show that, as the system size expands, even though the lowest VIEs keep decreasing, the changes of spectral function profile and the relative peak positions get unexpectedly smaller. Further analysis of the ionized states associated with the most intensive peak in the spectral functions reveals non-negligible |2h, 1p⟩’s in the ionized wave functions of the considered G–C base pair systems. The leading |2h, 1p⟩’s associated with the main ionizations from the cytosine part of the G–C base pairs feature a transition from the intra-base-pair cytosine π → π* excitation to the inter-base-pair electron excitation as the size of G–C base pairs expands, which also indicates the minimum quantum region in the many-body calculations of DNA systems.

Peng, Bo (ORCID:0000000242267294)↗

Roadmap on methods and software for electronic structure based simulations in chemistry and materials

This Roadmap article provides a succinct, comprehensive overview of the state of electronic structure methods and software for molecular and materials simulations. Seventeen distinct sections collect insights by 51 leading scientists in the field. Each contribution addresses the status of a particular area, as well as current challenges and anticipated future advances, with a particular eye towards software related aspects and providing key references for further reading. Foundational sections cover density functional theory and its implementation in real-world simulation frameworks, Green's function based many-body perturbation theory, wave-function based and stochastic electronic structure approaches, relativistic effects and semiempirical electronic structure theory approaches. Subsequent sections cover nuclear quantum effects, real-time propagation of the electronic structure, challenges for computational spectroscopy simulations, and exploration of complex potential energy surfaces. The final sections summarize practical aspects, including computational workflows for complex simulation tasks, the impact of current and future high-performance computing architectures, software engineering practices, education and training to maintain and broaden the community, as well as the status of and needs for electronic structure based modeling from the vantage point of industry environments. Overall, the field of electronic structure software and method development continues to unlock immense opportunities for future scientific discovery, based on the growing ability of computations to reveal complex phenomena, processes and properties that are determined by the make-up of matter at the atomic scale, with high precision.

36 MATERIALS SCIENCE↗

Reverberation mapping of the broad line region: application to a hydrodynamical line-driven disk wind solution

The latest analysis efforts in reverberation mapping are beginning to allow reconstruction of echo images (or velocity-delay maps) that encode information about the structure and kinematics of the broad line region (BLR) in active galactic nuclei (AGNs). Such maps can constrain sophisticated physical models for the BLR. The physical picture of the BLR is often theorized to be a photoionized wind launched from the AGN accretion disk. Previously we showed that the line-driven disk wind solution found in an earlier simulation by Proga and Kallman is virialized over a large distance from the disk. This finding implies that, according to this model, black hole masses can be reliably estimated through reverberation mapping techniques. However, predictions of echo images expected from line-driven disk winds are not available. Here, after presenting the necessary radiative transfer methodology, we carry out the first calculations of such predictions. We find that the echo images are quite similar to other virialized BLR models such as randomly orbiting clouds and thin Keplerian disks. We conduct a parameter survey exploring how echo images, line profiles, and transfer functions depend on both the inclination angle and the line opacity. We find that the line profiles are almost always single peaked, while transfer functions tend to have tails extending to large time delays. The outflow, despite being primarily equatorially directed, causes an appreciable blueshifted excess on both the echo image and line profile when seen from lower inclinations (i≲45{sup ∘}). This effect may be observable in low ionization lines such as Hβ.

79 ASTRONOMY AND ASTROPHYSICS↗

The bispectrum of matter perturbations from cosmic strings

We present the first calculation of the bispectrum of the matter perturbations induced by cosmic strings. The calculation is performed in two different ways: the first uses the unequal time correlators (UETCs) of the string network - computed using a Gaussian model previously employed for cosmic string power spectra. The second approach uses the wake model, where string density perturbations are concentrated in sheet-like structures whose surface density grows with time. The qualitative and quantitative agreement of the two gives confidence to the results. An essential ingredient in the UETC approach is the inclusion of compensation factors in the integration with the Green's function of the matter and radiation fluids, and we show that these compensation factors must be included in the wake model also. We also present a comparison of the UETCs computed in the Gaussian model, and those computed in the unconnected segment model (USM) used by the standard cosmic string perturbation package CMBACT. We compare numerical estimates for the bispectrum of cosmic strings to those produced by perturbations from an inflationary era, and discover that, despite the intrinsically non-Gaussian nature of string-induced perturbations, the matter bispectrum is unlikely to produce competitive constraints on a population of cosmic strings.

79 ASTRONOMY AND ASTROPHYSICS↗

Time Domain Moment Tensor Inversion in Python

MTtime is a software tool developed for time domain inversion of complete seismic waveform data to obtain the seismic moment tensor. It supports both deviatoric and full moment tensor inversions as well as 1-D and 3-D basis Green's functions.

Chiang, Andrea↗

Limitations of the clump-correlation theories of shear-induced turbulence suppression

The clump theory, primarily constructed by Dupree [Phys. Fluids 15, 334 (1972)] based on the moment approach and then generalized to the correlation theory [Y. Z. Zhang and S. M. Mahajan, Phys. Fluids B 5, 2000 (1993)], has long served as a basis for constructing theories of turbulence suppression by shear flow. In order to reveal the “intrinsic approximation” invoked in the clump-correlation theory, we examine a model based on two dimensional magnetized drift waves. After a rigorous derivation of the exact response function—a key to average the Green function of the system—we show that the Dupree, Zhang-Mahajan approach is recovered as the lowest order approximation in a small dimensionless parameter $Υ$ which is a triple product of the correlation time, wave number, and fluctuating drift velocity. The clump-correlation theory, thus, constitutes the Gaussian and lowest order non-Markovian process for a homogeneous stationary turbulence. Lastly we also provide, especially for the tokamak community, a readily usable formula to evaluate the effectiveness of shear-flow suppression; this formula pertains regardless of the specific model of correlation time.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Improved Earthquake Source Parameters with 3D Wavespeed Models in California and Nevada

Seismic tomography harnesses earthquake data to explore the inaccessible structure of the Earth. Adjoint waveform tomography (AWT), a method of seismic tomography, updates the tomographic model by optimizing the fit between observed earthquake data and synthetic waveforms. The synthetic data are calculated by solving the wave equation through a given 3D model. An important requirement to calculating synthetics is the source information (location, centroid time, depth, and moment tensor). Errors in source information affect the quality of the synthetics produced, which in turn can limit how structure can be inferred in the AWT workflow. Here, to test the effect of updating source information, we used MTTime (Chiang, 2020), a time-domain full-waveform moment tensor inversion code, to calculate the moment tensors and depths of 118 earthquakes that occurred in California and Nevada over a 20-yr period. We calculated 3D Green’s functions using a 3D seismic wavespeed model of California and Nevada (Doody et al., 2023b). We show that the inverted solutions provide better waveform fits than the Global Centroid Moment Tensor catalog and increase usable, well-correlated data by up to 7%. Therefore, we argue that recalculating source parameters should be considered in AWT workflows, particularly for smaller magnitude events (⁠M w > 5.0).

58 GEOSCIENCES↗

Gravitational radiation from radial infall of a particle into a Schwarzschild black hole: A numerical study of the spectra, quasinormal modes, and power-law tails

The computation of the gravitational radiation emitted by a particle falling into a Schwarzschild black hole is a classic problem that was already studied in the 1970s. Here we present a detailed numerical analysis of the case of radial infall starting at infinity with no initial velocity. We compute the radiated waveforms, spectra, and energies for multipoles up to l=6, improving significantly on the numerical accuracy of existing results. This is done by integrating the Zerilli equation in the frequency domain using the Green's function method. The resulting wave exhibits a ''ring-down'' phase whose dominant contribution is a superposition of the quasinormal modes of the black hole. The numerical accuracy allows us to recover the frequencies of these modes through a fit of that part of the wave. Comparing with direct computations of the quasinormal modes, we reach a {approx}10{sup -4} to {approx}10{sup -2} accuracy for the first two overtones of each multipole. Our numerical accuracy also allows us to display the power-law tail that the wave develops after the ring-down has been exponentially cut off. The amplitude of this contribution is {approx}10{sup 2} to {approx}10{sup 3} times smaller than the typical scale of the wave.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Experimental warming differentially affects vegetative and reproductive phenology of tundra plants

Rapid climate warming is altering Arctic and alpine tundra ecosystem structure and function, including shifts in plant phenology. While the advancement of green up and flowering are well-documented, it remains unclear whether all phenophases, particularly those later in the season, will shift in unison or respond divergently to warming. Here, we present the largest synthesis to our knowledge of experimental warming effects on tundra plant phenology from the International Tundra Experiment. We examine the effect of warming on a suite of season-wide plant phenophases. Results challenge the expectation that all phenophases will advance in unison to warming. Instead, we find that experimental warming caused: (1) larger phenological shifts in reproductive versus vegetative phenophases and (2) advanced reproductive phenophases and green up but delayed leaf senescence which translated to a lengthening of the growing season by approximately 3%. Patterns were consistent across sites, plant species and over time. The advancement of reproductive seasons and lengthening of growing seasons may have significant consequences for trophic interactions and ecosystem function across the tundra.

54 ENVIRONMENTAL SCIENCES↗

A Linear-Complexity Tensor Butterfly Algorithm for Compressing High-Dimensional Oscillatory Integral Operators

This paper presents a multilevel tensor compression algorithm called tensor butterfly algorithm for efficiently representing large-scale and high-dimensional oscillatory integral operators, including Green's functions for wave equations and integral transforms such as Radon transforms and Fourier transforms. The proposed algorithm leverages a tensor extension of the so-called complementary low-rank property of existing matrix butterfly algorithms. The algorithm partitions the discretized integral operator tensor into subtensors of multiple levels and factorizes each subtensor at the middle level as a Tucker-type interpolative decomposition, whose factor matrices are formed in a multilevel fashion. For a d-dimensional (d > 1) integral operator discretized into a 2d-mode tensor with n2d entries, the overall CPU time and memory requirement scale as O(nd), in stark contrast to the O(nd log n) complexity of existing matrix algorithms such as matrix butterfly algorithms and fast Fourier transforms (FFTs), where n is the number of points per direction. When comparing with other tensor algorithms such as quantized tensor train (QTT), the proposed algorithm also shows superior CPU and memory performance for tensor contraction. Remarkably, the tensor butterfly algorithm can efficiently model high-frequency Green's function interactions between two unit cubes, each spanning 512 wavelengths per direction, which represents problems of scale over 512× larger than that existing butterfly algorithms can handle, with the same amount of computation resources. On the other hand, for a problem representing 64 wavelengths per direction, which is the largest size existing algebraic matrix algorithms can handle, our tensor butterfly algorithm exhibits 200x speedups and 30× memory reduction compared with existing ones. Moreover, the tensor butterfly algorithm also permits O(nd)-complexity FFTs and Radon transforms up to d = 6 dimensions.

Kielstra, P Michael↗

Generalized and efficient algorithm for computing multipole energies and gradients based on Cartesian tensors

Accurate representation of intermolecular forces has been the central task of classical atomic simulations, known as molecular mechanics. Recent advancements in molecular mechanics models have put forward the explicit representation of permanent and/or induced electric multipole (EMP) moments. The formulas developed so far to calculate EMP interactions tend to have complicated expressions, especially in Cartesian coordinates, which can only be applied to a specific kernel potential function. For example, one needs to develop a new formula each time a new kernel function is encountered. The complication of these formalisms arises from an intriguing and yet obscured mathematical relation between the kernel functions and the gradient operators. Here, I uncover this relation via rigorous derivation and find that the formula to calculate EMP interactions is basically invariant to the potential kernel functions as long as they are of the form f(r), i.e., any Green’s function that depends on inter-particle distance. I provide an algorithm for efficient evaluation of EMP interaction energies, forces, and torques for any kernel f(r) up to any arbitrary rank of EMP moments in Cartesian coordinates. The working equations of this algorithm are essentially the same for any kernel f(r). Recently, a few recursive algorithms were proposed to calculate EMP interactions. Depending on the kernel functions, the algorithm here is about 4–16 times faster than these algorithms in terms of the required number of floating point operations and is much more memory efficient. I show that it is even faster than a theoretically ideal recursion scheme, i.e., one that requires 1 floating point multiplication and 1 addition per recursion step. This algorithm has a compact vector-based expression that is optimal for computer programming. The Cartesian nature of this algorithm makes it fit easily into modern molecular simulation packages as compared with spherical coordinate-based algorithms. A software library based on this algorithm has been implemented in C++11 and has been released.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Analytical Solutions of Transient Drift-Diffusion in P-N Junction Pixel Sensors

Radiation detection in applications ranging from high energy physics to medical imaging rely on solid state detectors, often hybrid pixel detectors with (1) reverse biased p-n junction pixel sensors and (2) readout ASICs, attached by flip chip-bonding. Transient signals characteristics are important in, e.g., matching ASIC and sensor design, modeling and optimizing detector parameters and describing timing and charge sharing properties. Currently analytical forms of transient signals are available for only a few limited cases (e.g., drift or diffusion) or for the steady state (which is not relevant for high energy radiation detection). Tools are available for (relatively slow) numerical evaluation of the transient charge transport. We present here the first analytical solutions of partial differential equations de scribing drift-diffusion-recombination charge transport in planar p-n junction sensors in a variety of conditions: (1) undepleted, (2) fully depleted, (3) taking into account the gradual velocity saturation, and (4) overdepleted. We deduce the Green’s functions which can be applied to any detection problem through simple convolution with the initial conditions. We compare the analytical solutions with Monte Carlo simulations and industry standard simulations (Synopsys Sentaurus), demonstrating good agreement. Using the analytical equations enables fast modeling of the influence of various detector parameters on tracking, imaging and timing performance, describing performance and enabling optimizations for different applications. Finally, we illustrate this model with applications in 3D+T (x,y,z,time) photon tracking and 4D+T (x,y,θ,φ,time) relativistic charged particle tracking.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Functional characterization of Brassica napus DNA topoisomerase Iα-1 and its effect on flowering time when expressed in Arabidopsis thaliana

Previous studies have shown that DNA topoisomerase Iα (AtTOP1α) has specific developmental functions during growth and development in Arabidopsis thaliana. However, little is known about the roles of DNA topoisomerases in the closely related and commercially important plant, rapeseed (Brassica napus). Here, the full-length BnTOP1α-1 coding sequence was cloned from the A2 subgenome of the Brassica napus inbred line L111. We determine that all BnTOP1α paralogs showed differing patterns of expression in different organs of L111, and that when expressed in tobacco leaves as a fusion protein with green fluorescent protein, BnTOP1α-1 localized to the nucleus. We further showed that ectopic expression of BnTOP1α-1 in the A. thaliana top1α-7 mutant fully complemented the early flowering phenotype of the mutant. Moreover, altered expression levels in top1α-7 seedlings of several key genes controlling flowering time were restored to wild type levels by ectopic expression of BnTOP1α-1. These results provide valuable insights into the roles of rapeseed DNA topoisomerases in flowering time, and provide a promising target for genetic manipulation of this commercially significant process in rapeseed. - Highlights: • BnTOP1α-1 was cloned from the A2 subgenome of Brassica napus inbred line L111. • BnTOP1α-1 rescued the early flowering phenotype of the Attop1α-7 mutant. • BnTOP1α-1 rescued the altered expression of flowering time genes in the Attop1α-mutant. • The functions of BnTOP1α-1 and AtTOP1α are likely conserved.

60 APPLIED LIFE SCIENCES↗

Rice husks as a sustainable source of nanostructured silicon for high performance Li-ion battery anodes

The recovery of useful materials from earth-abundant substances is of strategic importance for industrial processes. Despite the fact that Si is the second most abundant element in the Earth’s crust, processes to form Si nanomaterials is usually complex, costly and energy-intensive. Here we show that pure Si nanoparticles (SiNPs) can be derived directly from rice husks (RHs), an abundant agricultural byproduct produced at a rate of 1.2 3 108 tons/year, with a conversion yield as high as 5% by mass. And owing to their small size (10– 40 nm) and porous nature, these recovered SiNPs exhibits high performance as Li-ion battery anodes, with high reversible capacity (2,790 mA h g21 , seven times greater than graphite anodes) and long cycle life (86% capacity retention over 300 cycles). Using RHs as the raw material source, overall energy-efficient, green, and large scale synthesis of low-cost and functional Si nanomaterials is possible.

MATERIALS CHEMISTRY↗

Green's function-stochastic methods framework for probing nonlinear evolution problems: Burger's equation, the nonlinear Schroedinger's equation, and hydrodynamic organization of near-molecular-scale vorticity

Research Highlights: > Systematic approach for physically probing nonlinear and random evolution problems. > Evolution of vortex sheets corresponds to evolution of an Ornstein-Uhlenbeck process. > Organization of near-molecular scale vorticity mediated by hydrodynamic modes. > Framework allows calculation of vorticity evolution within random strain fields. - Abstract: A framework which combines Green's function (GF) methods and techniques from the theory of stochastic processes is proposed for tackling nonlinear evolution problems. The framework, established by a series of easy-to-derive equivalences between Green's function and stochastic representative solutions of linear drift-diffusion problems, provides a flexible structure within which nonlinear evolution problems can be analyzed and physically probed. As a preliminary test bed, two canonical, nonlinear evolution problems - Burgers' equation and the nonlinear Schroedinger's equation - are first treated. In the first case, the framework provides a rigorous, probabilistic derivation of the well known Cole-Hopf ansatz. Likewise, in the second, the machinery allows systematic recovery of a known soliton solution. The framework is then applied to a fairly extensive exploration of physical features underlying evolution of randomly stretched and advected Burger's vortex sheets. Here, the governing vorticity equation corresponds to the Fokker-Planck equation of an Ornstein-Uhlenbeck process, a correspondence that motivates an investigation of sub-sheet vorticity evolution and organization. Under the assumption that weak hydrodynamic fluctuations organize disordered, near-molecular-scale, sub-sheet vorticity, it is shown that these modes consist of two weakly damped counter-propagating cross-sheet acoustic modes, a diffusive cross-sheet shear mode, and a diffusive cross-sheet entropy mode. Once a consistent picture of in-sheet vorticity evolution is established, a number of analytical results, describing the motion and spread of single, multiple, and continuous sets of Burger's vortex sheets, evolving within deterministic and random strain rate fields, under both viscous and inviscid conditions, are obtained. In order to promote application to other nonlinear problems, a tutorial development of the framework is presented. Likewise, time-incremental solution approaches and construction of approximate, though otherwise difficult-to-obtain backward-time GF's (useful in solution of forward-time evolution problems) are discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗