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 145 records · Page 8

ManiNetCluster: a novel manifold learning approach to reveal the functional links between gene networks

BACKGROUND:The coordination of genomic functions is a critical and complex process across biological systems such as phenotypes or states (e.g., time, disease, organism, environmental perturbation). Understanding how the complexity of genomic function relates to these states remains a challenge. To address this, we have developed a novel computational method, ManiNetCluster, which simultaneously aligns and clusters gene networks (e.g., co-expression) to systematically reveal the links of genomic function between different conditions. Specifically, ManiNetCluster employs manifold learning to uncover and match local and non-linear structures among networks, and identifies cross-network functional links. RESULTS:We demonstrated that ManiNetCluster better aligns the orthologous genes from their developmental expression profiles across model organisms than state-of-the-art methods (p-value <2.2×10-16). This indicates the potential non-linear interactions of evolutionarily conserved genes across species in development. Furthermore, we applied ManiNetCluster to time series transcriptome data measured in the green alga Chlamydomonas reinhardtii to discover the genomic functions linking various metabolic processes between the light and dark periods of a diurnally cycling culture. We identified a number of genes putatively regulating processes across each lighting regime. CONCLUSIONS:ManiNetCluster provides a novel computational tool to uncover the genes linking various functions from different networks, providing new insight on how gene functions coordinate across different conditions. ManiNetCluster is publicly available as an R package at https://github.com/daifengwanglab/ManiNetCluster.

Comparative network analysis↗

Bayesian Inference for the Seismic Moment Tensor Using Regional Waveforms and Teleseismic- P Polarities with a Data-Derived Distribution of Velocity Models and Source Locations

The largest source of uncertainty in any source inversion is the velocity model used in the transfer function that relates observed ground motion to the seismic moment tensor. However, standard inverse procedure often does not quantify uncertainty in the seismic moment tensor due to error in the Green’s functions from uncertain event location and Earth structure. Here, we incorporate this uncertainty into an estimation of the seismic moment tensor using a data-derived distribution of velocity models based on complementary geophysical data sets, including thickness constraints, velocity profiles, gravity data, surface-wave group velocities, and regional body-wave travel times. The data-derived distribution of velocity models is then used as a prior distribution of Green’s functions for use in Bayesian inference of an unknown seismic moment tensor using regional and teleseismic-P waveforms. The use of multiple data sets is important for gaining resolution to different components of the moment tensor. The combined likelihood is estimated using data-specific error models and the posterior of the seismic moment tensor is estimated and interpreted in terms of the most probable source type.

58 GEOSCIENCES↗

Waveform inversion of acoustic waves for explosion yield estimation

We present a new waveform inversion technique to estimate the energy of near-surface explosions using atmospheric acoustic waves. Conventional methods often employ air blast models based on a homogeneous atmosphere, where the acoustic wave propagation effects (e.g., refraction and diffraction) are not taken into account, and therefore, their accuracy decreases with increasing source-receiver distance. In this study, three-dimensional acoustic simulations are performed with a finite difference method in realistic atmospheres and topography, and the modeled acoustic Green's functions are incorporated into the waveform inversion for the acoustic source time functions. The strength of the acoustic source is related to explosion yield based on a standard air blast model. The technique was applied to local explosions (<10 km) and provided reasonable yield estimates (<~30% error) in the presence of realistic topography and atmospheric structure. In conclusion, the presented method can be extended to explosions recorded at far distance provided proper meteorological specifications.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

Regional Moment Tensor Inversion Using Rotational Observations

There are benefits from the addition of rotational motions to translational displacements for moment tensor (MT) inversions. The rotational radiation pattern is orthogonal to the shear radiation pattern, thus incorporating rotations is equivalent to gaining another observation point on the focal sphere. We demonstrated this by simulating curl and displacement wavefields for a regional distance station. Thus, one 6-C station (3-Component translational + 3-Component rotational) gathers the same information on radiation pattern as two 3-C stations at 90° azimuth from one another along the focal plane axis, which is sometimes difficult to obtain when restricted to surface sensors. We added rotational Green's functions to a regional MT inversion scheme (long-period, time-domain, and linear inversion) by computing spatial gradients from f-κ reflectivity synthetics. For rotational data, we used Array Derived Rotations (ADR) from Piñon Flats Observatory Array and Golay array deployed during IRIS Community Wavefield Demonstration Experiment. The hope is to ultimately use compact and field deployable broadband rotational seismometers instead of ADRs. Rotational motions were predicted from well-constrained deviatoric MT solutions of nine earthquakes recorded by the arrays and other seismic networks. Additionally, we formed three station sparse datasets with (3-C and 6-C) and without rotational ground motions (3-C only) for eight earthquakes testing the benefits of including rotational waveforms in MT inversions when they are equally weighted with translational waveforms. Adding rotational motions improved the double-couple components and reduced compensated-linear-vector-dipole and isotropic components in full-MT solutions using sparse datasets with poor coverage.

58 GEOSCIENCES↗

A line source in Minkowski for the de Sitter spacetime scalar Green's function: Massless minimally coupled case

Motivated by the desire to understand the causal structure of physical signals produced in curved spacetimes – particularly around black holes – we show how, for certain classes of geometries, one might obtain its retarded or advanced minimally coupled massless scalar Green's function by using the corresponding Green's functions in the higher dimensional Minkowski spacetime where it is embedded. Analogous statements hold for certain classes of curved Riemannian spaces, with positive definite metrics, which may be embedded in higher dimensional Euclidean spaces. The general formula is applied to (d ≥ 2)-dimensional de Sitter spacetime, and the scalar Green's function is demonstrated to be sourced by a line emanating infinitesimally close to the origin of the ambient (d + 1)-dimensional Minkowski spacetime and piercing orthogonally through the de Sitter hyperboloids of all finite sizes. This method does not require solving the de Sitter wave equation directly. Only the zero mode solution to an ordinary differential equation, the “wave equation” perpendicular to the hyperboloid – followed by a one-dimensional integral – needs to be evaluated. A topological obstruction to the general construction is also discussed by utilizing it to derive a generalized Green's function of the Laplacian on the (d ≥ 2)-dimensional sphere.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Time Evolution of Correlation Functions in Quantum Many-Body Systems

We give rigorous analytical results on the temporal behavior of two-point correlation functions—also known as dynamical response functions or Green’s functions—in closed many-body quantum systems. We show that in a large class of translation-invariant models the correlation functions factorize at late times $\langle$A ( t ) B$\rangle$ β → $\langle$ A$\rangle$ β $\langle$ B$\rangle$ β , thus proving that dissipation emerges out of the unitary dynamics of the system. We also show that for systems with a generic spectrum the fluctuations around this late-time value are bounded by the purity of the thermal ensemble, which generally decays exponentially with system size. For autocorrelation functions we provide an upper bound on the timescale at which they reach the factorized late time value. Remarkably, this bound is only a function of local expectation values and does not increase with system size. We give numerical examples that show that this bound is a good estimate in nonintegrable models, and argue that the timescale that appears can be understood in terms of an emergent fluctuation-dissipation theorem. Our study extends to further classes of two point functions such as the symmetrized ones and the Kubo function that appears in linear response theory, for which we give analogous results.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Electromagnetic and neutral-weak response functions of 4 He and 12 C

A major goal of nuclear theory is to understand the strong interaction in nuclei as it manifests itself in terms of two- and many-body forces among the nuclear constituents, the protons and neutrons, and the interactions of these constituents with external electroweak probes via one- and many-body currents. The objective of the present work is to calculate the quasielastic electroweak response functions in light nuclei within the realistic dynamical framework outlined above. These response functions determine the inclusive cross section as function of the lepton momentum and energy transfers. Their ab initio calculation is a very challenging quantum many-body problem, since it requires summation over the entire excitation spectrum of the nucleus and inclusion in the electroweak currents of one- and many-body terms. Green's functions Monte Carlo methods allow one to circumvent both difficulties by computing the response in imaginary time (the so-called Euclidean response) and hence summing implicitly over the bound and continuum states of the nucleus, and by implementing specific algorithms designed to deal with the complicated spin-isospin structure of nuclear many-body operators. Theoretical predictions for 4 He and 12 C, confirmed by experiment in the electromagnetic case, show that two-body currents generate excess transverse strength from threshold to the quasielastic to the dip region and beyond. In conclusion, these results challenge the conventional picture of quasielastic inclusive scattering as being largely dominated by single-nucleon knockout processes.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Many-body effects in transport and energy transfer

Many-body physics is common in defining responses of nanoscale to external perturbations such as bias, optical driving, temperature gradient. When applied to open systems these perturbations lead to matter and energy transfer between the system and its surroundings (baths) and/or cause intra-system dynamics. While many different theoretical techniques have been developed to describe system responses, not all of them are convenient (or even appropriate) when applied at nanoscale. Traditional theoretical techniques are either based on kinetic schemes (e.g., Redfield-Lindblad quantum master equation) or utilize quasiparticle description (e.g., single-particle scattering theory, input-output formalism, standard nonequilibrium Green’s functions technique). The former fails to account for quantum mechanics of the system-bath interface and may fail qualitatively when used in open interacting systems. The latter is inconvenient when many-body states description of the system is preferable, which is usual, e.g., in molecular optical spectroscopy. This project was devoted to development of new nonequilibrium Green’s function (NEGF) techniques and application of standard and many-body flavors of the NEGF to description of charge and energy transport, optical response, nonadiabatic molecular dynamics, and thermodynamics of nanoscale devices in general, and, in particular, of single-molecule junctions. In particular, Hubbard NEGF (new many-body flavor of the NEGF) was introduced as a new technique especially convenient in treatment of optoelectronic experiments in single-molecule junctions and useful in discussion of bias induced nuclear dynamics in junctions. The research also led to further technique developments and applications beyond the timing of this Award.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

Instanton liquid properties from lattice QCD

Here, we examined the instanton contribution to the QCD configurations generated from lattice QCD for N F = 0, N F = 2 + 1 and N F = 2 + 1 + 1 dynamical quark flavors from two different and complementary approaches. First via the use of Gradient flow, we computed instanton liquid properties using an algorithm to localize instantons in the gauge field configurations and studied their evolution with flow time. Then, the analysis of the running at low momenta of gluon Green's functions serves as an independent confirmation of the instanton density which can also be derived without the use of the Gradient flow.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

NESSi : The N on- E quilibrium S ystems S imulation package

The nonequilibrium dynamics of correlated many-particle systems is of interest in connection with pump–probe experiments on molecular systems and solids, as well as theoretical investigations of transport properties and relaxation processes. Nonequilibrium Green’s functions are a powerful tool to study interaction effects in quantum many-particle systems out of equilibrium, and to extract physically relevant information for the interpretation of experiments. Here, we present the open-source software package NESSi (The Non-Equilibrium Systems Simulation package) which allows to perform many-body dynamics simulations based on Green’s functions on the L-shaped Kadanoff–Baym contour. NESSi contains the library libcntr which implements tools for basic operations on these nonequilibrium Green’s functions, for constructing Feynman diagrams, and for the solution of integral and integro-differential equations involving contour Green’s functions. The library employs a discretization of the Kadanoff–Baym contour into time points and a high-order implementation of integration routines. The total integrated error scales up to $\mathcal{O}(N^{-7})$, which is important since the numerical effort increases at least cubically with the simulation time. A distributed-memory parallelization over reciprocal space allows large-scale simulations of lattice systems. We provide a collection of example programs ranging from dynamics in simple two-level systems to problems relevant in contemporary condensed matter physics, including Hubbard clusters and Hubbard or Holstein lattice models. The libcntr library is the basis of a follow-up software package for nonequilibrium dynamical mean-field theory calculations based on strong-coupling perturbative impurity solvers.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Adjoint method and runaway electron avalanche

The adjoint method for the study of runaway electron dynamics in momentum space Liu et al (2016 Phys. Plasmas 23 010702) is rederived using the Green's function method, for both the runaway probability function (RPF) and the expected loss time (ELT). The RPF and ELT obtained using the adjoint method are presented, both with and without the synchrotron radiation reaction force. In conclusion, the adjoint method is then applied to study the runaway electron avalanche. Both the critical electric field and the growth rate for the avalanche are calculated using this fast and novel approach.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Summability of Connected Correlation Functions of Coupled Lattice Fields

We consider two nonindependent random fields ψ and ϕ defined on a countable set Z. For instance, Z=ℤ{sup d} or Z=ℤ{sup d}×I, where I denotes a finite set of possible “internal degrees of freedom” such as spin. We prove that, if the cumulants of ψ and ϕ enjoy a certain decay property, then all joint cumulants between ψ and ϕ are ℓ{sub 2}-summable in the precise sense described in the text. The decay assumption for the cumulants of ψ and ϕ is a restricted ℓ{sub 1} summability condition called ℓ{sub 1}-clustering property. One immediate application of the results is given by a stochastic process ψ{sub t}(x) whose state is ℓ{sub 1}-clustering at any time t: then the above estimates can be applied with ψ=ψ{sub t} and ϕ=ψ{sub 0} and we obtain uniform in t estimates for the summability of time-correlations of the field. The above clustering assumption is obviously satisfied by any ℓ{sub 1}-clustering stationary state of the process, and our original motivation for the control of the summability of time-correlations comes from a quest for a rigorous control of the Green–Kubo correlation function in such a system. A key role in the proof is played by the properties of non-Gaussian Wick polynomials and their connection to cumulants.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Approximating and incorporating model uncertainty in an inversion for seismic source functions: Preliminary results

We present preliminary work on propagating model uncertainty into the estimation of the time domain source time functions of the seismic source. Our method is based on an estimated model covariance function, which we estimate from the data. The model covariance function is then used to construct a suite of surrogate Greens functions which we use in a Monte Carlo type inversion scheme. The result is a probability density function of the six independent source time functions, each of which corresponds to an individual component of the seismic moment tensor. We compare the results of our method with those obtained using a computationally expensive finite difference Monte Carlo method and find that our new method produces results that are deficient in low frequencies. The advantage of our new method, which we term the Karhunen-Loeve Monte Carlo (KLMC) method, is that is several orders of magnitude faster than our current method, which uses a finite difference scheme to produce the suite of forward models.

42 ENGINEERING↗

Asymptotic neutron scattering laws for anomalously diffusing quantum particles

The paper deals with a model-free approach to the analysis of quasielastic neutron scattering intensities from anomalously diffusing quantum particles. All quantities are inferred from the asymptotic form of their time-dependent mean square displacements which grow ∝t{sup α}, with 0 ≤ α < 2. Confined diffusion (α = 0) is here explicitly included. We discuss in particular the intermediate scattering function for long times and the Fourier spectrum of the velocity autocorrelation function for small frequencies. Quantum effects enter in both cases through the general symmetry properties of quantum time correlation functions. It is shown that the fractional diffusion constant can be expressed by a Green-Kubo type relation involving the real part of the velocity autocorrelation function. The theory is exact in the diffusive regime and at moderate momentum transfers.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Light-Enhanced Spin Fluctuations and d -Wave Superconductivity at a Phase Boundary

Time-domain techniques have shown the potential of photomanipulating existing orders and inducing new states of matter in strongly correlated materials. Using time-resolved exact diagonalization, we perform in this paper numerical studies of pump dynamics in a Mott-Peierls system with competing charge and spin density waves. A light-enhanced d-wave superconductivity is observed when the system resides near a quantum phase boundary. By examining the evolution of spin, charge, and superconducting susceptibilities, we show that a subdominant state in equilibrium can be stabilized by photomanipulating the charge order to allow superconductivity to appear and dominate. Finally, this work provides an interpretation of light-induced superconductivity from the perspective of order competition and offers a promising approach for designing novel emergent states out of equilibrium.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

A Green’s function fast multipole method for computation of micromechanical fields in heterogeneous materials

Computation of micromechanical fields in heterogeneous materials is usually performed using either the finite element method or the Green’s function method based on FFTs. The finite element method allows for accurate discretization and for non-periodic boundary conditions but is computationally expensive. On the other hand, the FFT-based method is computationally efficient but requires discretization on a regular grid of hexahedral voxels. In this paper, a Green’s function method allowing for accurate discretization using tetrahedral elements and for non-periodic boundary conditions is proposed. The convolution is computed using the fast multipole method, which provides good accuracy even for low-order expansion due to the fast decay of interactions between elements. The proposed Green’s function fast multipole method is verified by comparison with analytical and FFT-based solutions. Furthermore, the computational time is analyzed and compared to the FFT-based method for non-periodic convolution. Finally, effective properties of an elastic polycrystalline microstructure containing thin intergranular cracks are computed and analyzed.

36 MATERIALS SCIENCE↗

Discrete Green’s functions and spectral graph theory for computationally efficient thermal modeling

Here, this work concerns solutions of the heat equation with the spectral graph method, for which the temperature is defined at discrete points in the domain and the spatial relationship among the points is described by a graph. The heat equation on the graph is solved using matrix techniques involving the eigenvectors and eigenvalues of the Laplacian matrix. The spectral graph approach precludes the computationally intensive meshing and numerous time-integration steps of the finite element method. In the present work, the spectral graph method is extended to include heat loss at the boundaries with a generalized boundary condition, and physics-based edge weights are introduced which simplify the calibration process. From this approach a discrete Green’s function is defined which allows for solutions under a variety of heating conditions including: space-varying initial conditions; time-and-space varying internal heating; and, time-and-space-varying heating at boundaries of type 1 (Dirichlet), type 2 (Neumann) and type 3 (Robin). Results are provided for benchmark heat transfer problems in one spatial dimension and in three spatial dimensions, and verification is provided by comparison with exact analytical solutions and finite difference solutions. The spectral graph method converges within 0.4% error of the analytical solution. The practical utility of the approach is demonstrated by thermal simulation of a multilayer additive manufacturing process. The spectral graph results are compared to experimentally-obtained temperature data for two metal parts, with error less than 5% of the experimental measurements, with computation time less than one minute on a desktop computer.

36 MATERIALS SCIENCE↗

Exact calculation of the time convolutionless master equation generator: Application to the nonequilibrium resonant level model

The generalized quantum master equation provides a powerful tool to describe the dynamics in quantum impurity models driven away from equilibrium. Two complementary approaches, one based on Nakajima–Zwanzig–Mori time-convolution (TC) and the other on the Tokuyama–Mori time-convolutionless (TCL) formulations provide a starting point to describe the time-evolution of the reduced density matrix. A key in both approaches is to obtain the so called “memory kernel” or “generator,” going beyond second or fourth order perturbation techniques. While numerically converged techniques are available for the TC memory kernel, the canonical approach to obtain the TCL generator is based on inverting a super-operator in the full Hilbert space, which is difficult to perform and thus, nearly all applications of the TCL approach rely on a perturbative scheme of some sort. Here, the TCL generator is expressed using a reduced system propagator which can be obtained from system observables alone and requires the calculation of super-operators and their inverse in the reduced Hilbert space rather than the full one. This makes the formulation amenable to quantum impurity solvers or to diagrammatic techniques, such as the nonequilibrium Green’s function. We implement the TCL approach for the resonant level model driven away from equilibrium and compare the time scales for the decay of the generator with that of the memory kernel in the TC approach. Furthermore, the effects of temperature, source-drain bias, and gate potential on the TCL/TC generators are discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗