Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “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 451 records · Page 25
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

A numerical Poisson solver with improved radial solutions for a self-consistent locally scaled self-interaction correction method

Abstract The universal applicability of density functional approximations is limited by self-interaction error made by these functionals. Recently, a novel one-electron self-interaction-correction (SIC) method that uses an iso-orbital indicator to apply the SIC at each point in space by scaling the exchange-correlation and Coulomb energy densities was proposed. The locally scaled SIC (LSIC) method is exact for the one-electron densities, and unlike the well-known Perdew–Zunger SIC (PZSIC) method recovers the uniform electron gas limit of the uncorrected density functional approximation, and reduces to PZSIC method as a special case when isoorbital indicator is set to the unity. Here, we present a numerical scheme that we have adopted to evaluate the Coulomb potential of the electron density scaled by the iso-orbital indicator required for the self-consistent LSIC calculations. After analyzing the behavior of the finite difference method (FDM) and the green function solution to the radial part of the Poisson equation, we adopt a hybrid approach that uses the FDM for the Coulomb potential due to the monopole and the GF for all higher-order terms. The performance of the resultant hybrid method is assessed using a variety of systems. The results show improved accuracy than earlier numerical schemes. We also find that, even with a generic set of radial grid parameters, accurate energy differences can be obtained using a numerical Coulomb solver in standard density functional studies.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

First-principles study of L -shell iron and chromium opacity at stellar interior temperatures

Recently developed free-energy density functional theory (DFT)-based methodology for optical property calculations of warm dense matter has been applied for studying L-shell opacity of iron and chromium at T = 182 eV. We use Mermin–Kohn–Sham density functional theory with a ground-state and a fully temperature-dependent generalized gradient approximation exchange-correlation (XC) functionals. It is demonstrated that the role of XC at such a high-T is negligible due to the total free-energy of interacting system being dominated by the noninteracting free-energy term in agreement with estimations for the homogeneous electron gas. Furthermore, our DFT predictions are compared to the radiative emissivity and opacity of dense plasmas model, to the real-space Greens function method, and to experimental measurements. Good agreement is found between all three theoretical methods, and in the bound–continuum region for Cr when compared to the experiment, while the discrepancy between direct DFT calculations and the experiment for Fe remains essentially the same as for plasma-physics models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Efficient separation and coprecipitation for simplified cathode recycling

Hydrometallurgical recycling of spent lithium-ion batteries is among the most promising recovery approaches. The current hydrometallurgical method for battery recycling faces challenges such as complex separation and precipitation processes, and environmental concerns from the use of caustic inorganic acid and hazardous hydrogen peroxide. Our proposed modified method, namely polyol-metallurgy, uses citric acid in ethylene glycol as dual-function green solution to overcome these obstacles. This bifunctional solution leaches valuable metal ions from cathode materials (e.g., LiCoO 2 ), and then acts as chelating agents to selectively precipitate cobalt through an esterification reaction, without the need for additional precipitation agents. The leaching efficiency of cobalt and lithium reaches 99.55% and 97.65%, respectively, and more than 96% of cobalt could be directly self-precipitated and recovered. Further, the unique characteristics of the dual-function solution also avoids impurities from Al foils and PVDF/carbon black films, enabling the simple separation process. Therefore, the process can be completed in one-pot system with efficient leaching, separation and coprecipitation, making it more efficient and feasible for operation than existing alternatives.

25 ENERGY STORAGE↗
DOE OSTI2018

Quantum oscillation in carrier transport in two-dimensional junctions

Two-dimensional (2D) junction devices have recently attracted considerable attention. In this report we show that most 2D junction structures, whether vertical or lateral, act as a lateral monolayer–bilayer–monolayer junction in their operation. In particular, a vertical structure cannot function as a vertical junction as having been widely believed in the literature. Due to a larger electrostatic screening, the bilayer region in the junction always has a smaller bandgap than its monolayer counterpart. As a result, a potential well, aside from the usual potential barrier, will form universally in the bilayer region to affect the hole or electron quantum transport in the form of transmission or reflection. Taking black phosphorus as an example, our calculations using a non-equilibrium Green function combined with density functional theory show a distinct oscillation in the transmission coefficient in a two-electrode prototypical device, and the results can be qualitatively understood using a simple quantum well model.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Geometries of edge and mixed dislocations in bcc Fe from first principles calculations

We use density functional theory (DFT) to compute the core structures of $a_0[100](010)$ edge, $a_0[100](011)$ edge, $a_0/2[\bar{1}\bar{1}1](1\bar{1}0)$ edge, and $a_0/2[111](1\bar{1}0)$ $71^{\circ}$ mixed dislocations in body-centered cubic (bcc) Fe. The calculations are performed using flexible boundary conditions (FBC), which effectively allow the dislocations to relax as isolated defects by coupling the DFT core to an infinite harmonic lattice through the lattice Green function (LGF). We use the LGFs of the dislocated geometries in contrast to most previous FBC-based dislocation calculations that use the LGF of the bulk crystal. The dislocation LGFs account for changes in the topology of the crystal in the core as well as local strain throughout the crystal lattice. A simple bulk-like approximation for the force constants in a dislocated geometry leads to dislocation LGFs that optimize the core structures of the $a_0[100](010)$ edge, $a_0[100](011)$ edge, and $a_0/2[111](1\bar{1}0)$ $71^{\circ}$ mixed dislocations. This approximation fails for the $a_0/2[\bar{1}\bar{1}1](1\bar{1}0)$ dislocation however, so in this case we derive the LGF from more accurate force constants computed using a Gaussian approximation potential. The standard deviations of the dislocation Nye tensor distributions quantify the widths of the dislocation cores. The relaxed cores are compact, and the local magnetic moments on the Fe atoms closely follow the volumetric strain distributions in the cores. We also compute the core structures of these dislocations using eight different classical interatomic potentials, and quantify symmetry differences between the cores using the Fourier coefficients of their Nye tensor distributions. Most of the core structures computed using the classical potentials agree well with the DFT results. Furthermore the DFT core geometries provide benchmarking for classical potential studies of work-hardening, as well as substitutional and interstitial sites for computing solute-dislocation interactions that serve as inputs for mesoscale models of solute strengthening and solute diffusion near dislocations.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Chemical structure–property relationships in nanocelluloses

Abstract Nanocelluloses, nanoscale cellulosic particles or fibrils, are an appealing class of nanomaterials with enormous potential for sustainable development. They have great promise to target some of the 17 Sustainable Development Goals outlined by the United Nations. Nanocelluloses are derived from the most abundant biopolymer in the world, cellulose, which have multiple hydroxyl groups on their surface, enabling a vast range of chemical modifications. The interplay between the chemical structure and physicochemical properties of nanocelluloses is crucial to engineering the next generation of green, functional nanomaterials. In this review paper, we provide a comprehensive account of the structure–property relationships in nanocelluloses, which may establish the basis for the design of precision bio‐based materials. Chemical modifications of nanocelluloses, including hydrolysis, oxidation, esterification, amidation, and polymer grafting, provide an array of chemical and physical properties that open opportunities in a diverse spectrum of applications. We review how chemical modifications affect the physicochemical properties of nanocelluloses, such as thermal stability, hydrophobicity, and dispersibility in nonpolar media. Understanding nanocellulose chemical structure–property relationships is integral to the development of sustainable material platforms based on the most renewable biopolymer on earth.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Notes on the complex Sachdev-Ye-Kitaev model

We describe numerous properties of the Sachdev-Ye-Kitaev model for complex fermions with N $\gg$ 1 flavors and a global U(1) charge. We provide a general definition of the charge in the ( G, Σ) formalism, and compute its universal relation to the infrared asymmetry of the Green function. The same relation is obtained by a renormalization theory. The conserved charge contributes a compact scalar field to the effective action, from which we derive the many-body density of states and extract the charge compressibility. We compute the latter via three distinct numerical methods and obtain consistent results. Finally, we present a two dimensional bulk picture with free Dirac fermions for the zero temperature entropy.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

Time-dependent observables in heavy ion collisions. Part II. In search of pressure isotropization in the φ 4 theory

To understand the dynamics of thermalization in heavy ion collisions in the perturbative framework it is essential to first find corrections to the free-streaming classical gluon fields of the McLerran-Venugopalan model. The corrections that lead to deviations from free streaming (and that dominate at late proper time) would provide evidence for the onset of isotropization (and, possibly, thermalization) of the produced medium. To find such corrections we calculate the late-time two-point Green function and the energy-momentum tensor due to a single 2 → 2 scattering process involving two classical fields. To make the calculation tractable we employ the scalar φ 4 theory instead of QCD. We compare our exact diagrammatic results for these quantities to those in kinetic theory and find disagreement between the two. The disagreement is in the dependence on the proper time τ and, for the case of the two-point function, is also in the dependence on the space-time rapidity η: the exact diagrammatic calculation is, in fact, consistent with the free streaming scenario. Kinetic theory predicts a build-up of longitudinal pressure, which, however, is not observed in the exact calculation. We conclude that we find no evidence for the beginning of the transition from the free-streaming classical fields to the kinetic theory description of the produced matter after a single 2 → 2 rescattering.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2019

Generalized asymptotics for gauge fields

An interesting question is to characterize the general class of allowed boundary conditions for gauge theories, including gravity, at spatial and null infinity. This has played a role in discussions of soft charges, where antipodal symmetry has typically been assumed. However, the existence of electric and gravitational line operators, arising from gauge-invariant dressed observables, for example associated to axial or Fefferman-Graham like gauges, indicates the existence of non-antipodally symmetric initial data. This note studies aspects of the solutions corresponding to such non-symmetric initial data. The explicit evolution can be found, via a Green function, and bounds can be given on the asymptotic behavior of such solutions, evading arguments for singular behavior. Likewise, objections to such solutions based on infinite symplectic form are also avoided, although these solutions may be superselected. Soft charge conservation laws, and their modification, are briefly examined for such solutions. This discussion strengthens (though is not necessary for) arguments that soft charges characterize gauge field degrees of freedom, but not necessarily the degrees of freedom associated to the matter sourcing the field.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2019

Three-dimensional nonlocal anisotropic elasticity: a generalized continuum theory of Ångström-mechanics

In this paper, based on Eringen’s theory of nonlocal anisotropic elasticity, the three-dimensional nonlocal anisotropic elasticity of generalized Helmholtz type is developed. The derivation of a new three-dimensional nonlocal anisotropic kernel, which is the Green function of the three-dimensional anisotropic Helmholtz equation, enables to capture anisotropic length scale effects by means of a length scale tensor, which is a symmetric tensor of rank two. The derived nonlocal kernel function possesses up to six internal characteristic lengths on the Ångström-scale. The presented theory of nonlocal elasticity possesses the appropriate property to be a generalized continuum theory of Ångström-mechanics, since the range of its validity and applicability is up to the Ångström-scale. The connection between the theory of nonlocal anisotropic elasticity and lattice theory is established. The tensor function of nonlocal elastic moduli as well as the nonlocal kernel function is given in terms of the Hessian matrix in the lattice approach. In the framework of the considered theory, the modeling of dislocations in anisotropic materials taking into consideration anisotropic dislocation core effects is presented. Important dislocation key formulas, namely the anisotropic Peach–Koehler stress formula, the Peach–Koehler force and the anisotropic Blin’s formula, are derived. A major tool used in deriving the expression of anisotropic Blin’s formula is Kirchner’s so-called FF-tensor, which is here generalized toward nonlocal anisotropic elasticity. The main feature and advantage of the derived fields, compared with the corresponding ones in classical anisotropic elasticity, is that they are free of singularities. Numerical applications to straight dislocations in bcc Fe are given, revealing the ability and advantage of the considered theory to describe adequately nonsingular anisotropic stress and self-energy fields capturing the effects of anisotropy on the Ångström-scale.

42 ENGINEERING↗
DOE OSTI2016

Mechanisms of NH3 and NO2 detection in carbon-nanotube-based sensors: An ab initio investigation

The mechanisms of NH3 and NO2 detection by single-walled carbon nanotube-based devices are investigated by ab initio calculations and the non-equilibrium Greens function (NEGF) methodology. While both NH3 and NO2 can physisorb to a pristine carbon nanotube, we show that their adsorption only results in small current changes through the device. For a carbon nanotube (CNT) attached to gold nanowire leads, the most sensitive detection site is at the CNT near the CNT-Au contact, where chemisorption occurs. The resulting change in electron transmission and low-bias current can lead to over 30% sensitivity. While both NH3 and NO2 can also chemisorb at the Au electrodes, their adsorption results in only a small change in the plurality of the conducting levels of the gold layers, and thus a small effect on current. In order to enhance the detection sensitivity, it is thus beneficial to mask the electrodes to prevent chemisorption. Furthermore, the length of the pure CNT segment does not strongly affect the relative sensitivity. Our results suggest that a short-CNT device with exposed contact regions and masked electrodes would have the greatest sensitivity.

Chemistry↗
DOE OSTI2020

Monte Carlo simulations for understanding the transport properties of metallic nanowires

We report the conductance calculation of metallic nanowires of different metals under the application of a tensile force. The elongation induced by the tensile force was simulated by means of the canonical Monte Carlo method combined with the embedded atom method to compute the energy of the system. We studied nanowires composed of Ag, Au, Pt and Cu. The conductance was calculated by means of the Landauer formula and the non-equilibrium Green function formalism. As previously verified experimentally and using other simulation techniques, the conductance is quantized. In addition to this, we have found that there is a strong correlation between the magnitude of the conductance and the width of the nanowire at its thinnest point. Conductance histograms were analyzed for all the different metals.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Photonics and spectroscopy in nanojunctions: a theoretical insight

The progress of experimental techniques at the nanoscale in the last decade made optical measurements in current-carrying nanojunctions a reality, thus indicating the emergence of a new field of research coined optoelectronics. Optical spectroscopy of open nonequilibrium systems is a natural meeting point for (at least) two research areas: nonlinear optical spectroscopy and quantum transport, each with its own theoretical toolbox. We review recent progress in the field comparing theoretical treatments of optical response in nanojunctions as is accepted in nonlinear spectroscopy and quantum transport communities. A unified theoretical description of spectroscopy in nanojunctions is presented. Here, we argue that theoretical approaches of the quantum transport community (and in particular, the Green function based considerations) yield a convenient tool for optoelectronics when the radiation field is treated classically, and that differences between the toolboxes may become critical when studying the quantum radiation field in junctions.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2019

Improving the efficiency of G 0 W 0 calculations with approximate spectral decompositions of dielectric matrices

Recently, it was shown that the calculation of quasiparticle energies using the G 0 W 0 approximation can be performed without computing explicitly any virtual electronic states, by expanding the Green function and screened Coulomb interaction in terms of the eigenstates of the static dielectric matrix. Avoiding the evaluation of virtual electronic states leads to improved efficiency and ease of convergence of G 0 W 0 calculations. Here, we propose a further improvement of the efficiency of these calculations, based on an approximation of density-density response functions of molecules and solids. The approximation relies on the calculation of a subset of eigenvectors of the dielectric matrix using the kinetic operator instead of the full Hamiltonian, and it does not lead to any substantial loss of accuracy for the quasiparticle energies. The computational savings introduced by this approximation depend on the system, and they become more substantial as the number of electrons increases.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Two-dimensional shaping of Solov'ev equilibrium with vacuum using external coils

In this work, we demonstrate a method for constructing the Solov'ev equilibrium with any given 2D shape surrounded by a vacuum region using external poloidal field coils, based on the method previously developed by T. Xu and R. Fitzpatrick [Nuclear Fusion, 59, 064002 (2019)]. The computational domain consists of two parts: the plasma region, where the solution is the same as the Solov'ev solution, and the vacuum region, where the magnetic field generated by external coils as well as plasma current is determined using the Green function method through a matching condition near the separatrix. However, the method is not limited to the Solov'ev equilibrium in particular. In conclusion, the accuracy, efficiency, and robustness of such a scheme suggest that this method may be applied to the 2D shaping of a tokamak plasma with a vacuum region using external coils in general.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

The effect of curvature induced broken potential vorticity conservation on drift wave turbulences

Abstract The correlation theory of turbulence suppression (Zhang and Mahajan 1993 Phys. Fluids B5 2000) by velocity shear was constructed by invoking the ansatz of potential vorticity conservation (PVC) that holds for relatively simple (slab) models of drift wave turbulence. It is, therefore, surprising that a detailed modern simulation of the H mode pedestal, using the gyrokinetic code GENE (Hatch et al 2018 Plasma Phys. Control. Fusion 60 084003), found ‘striking agreement’ with the predictions of the analytic model. To understand the reasons for this remarkable agreement, an extended theory that contains finite (magnetic) curvature, and which does not conserve potential vorticity, is developed and ‘solved’ by calculating an inhomogeneous Green function reflecting the fact that the new system has a potential vorticity source. It is, then, demonstrated that the effect of the broken PVC is insignificant for the normal operation parameters in tokamaks; the correction due to curvature is at the order of 2 L 0 / R , where L 0 is the scale length of the local gradient and R is the major radius. The excellent agreement between simulation and slab model 1993 theory is therefore not accidental; the latter can be applied with confidence to the tokamak pedestal.

Physics↗
DOE OSTI2022

Imaging the crust and uppermost mantle structure of Portugal (West Iberia) with seismic ambient noise

SUMMARY We present a new high-resolution 3-D shear wave velocity (Vs) model of the crust and uppermost mantle beneath Portugal, inferred from ambient seismic noise tomography. We use broad-band seismic data from a dense temporary deployment covering the entire Portuguese mainland between 2010 and 2012 in the scope of the WILAS project. Vertical component data are processed using phase correlation and phase weighted stack to obtain empirical Green functions (EGFs) for 2016 station pairs. Further, we use a random sampling and subset stacking strategy to measure robust Rayleigh-wave group velocities in the period range 7–30 s and associated uncertainties. The tomographic inversion is performed in two steps: First, we determine group-velocity lateral variations for each period. Next, we invert them at each grid point using a new trans-dimensional inversion scheme to obtain the 3-D shear wave velocity model. The final 3-D model extends from the upper crust (5 km) down to the uppermost mantle (60 km) and has a lateral resolution of ∼50 km. In the upper and middle crusts, the Vs anomaly pattern matches the tectonic units of the Variscan Massif and Alpine basins. The transition between the Lusitanian Basin and the Ossa Morena Zone is marked by a contrast between moderate- and high-velocity anomalies, in addition to two arched earthquake lineations. Some faults, namely, the Manteigas–Vilariça–Bragança fault and the Porto–Tomar–Ferreira do Alentejo fault, have a clear signature from the upper crust down to the uppermost mantle (60 km). Our 3-D shear wave velocity model offers new insights into the continuation of the main tectonic units at depth and contributes to better understanding the seismicity of Portugal.

Silveira, Graça (ORCID:0000000221102554)↗
DOE OSTI2020

Pairing instability on a Luttinger surface: A non-Fermi liquid to superconductor transition and its Sachdev-Ye-Kitaev dual

Superconductivity results from an instability of the Fermi surface—contour of poles of the single-particle propagator—to an infinitesimally small attraction between electrons. Here, we instead discuss the analogous problem on a model Luttinger surface, or contour of zeros of the Green function. At zero temperature ($β → \infty$) and a critical interaction strength ($u_{c∞}$) characterized by the residue of self-energy pole, we find that the pair susceptibility diverges leading to a superconducting instability. We evaluate the pair fluctuation partition function and find that the spectral density in the normal state has an interaction-driven, power-law $\frac{1}{\sqrt{ω}}$ type, van-Hove singularity (vHS), indicating non-Fermi liquid (NFL) physics. Crucially, in the strong coupling limit ($βu \gg 1$), the leading-order fluctuation free-energy terms in the normal state of this NFL-SC transition resembles the equivalent [$\textit{O}(1)$] terms of the Sachdev-Ye-Kitaev (SYK) model. This free-energy contribution takes a simple form $–βF = βu_{c∞} – γ \text{ln}(βu_{c∞}$), where $\gamma$ is a constant equal to $\frac{1}{2}$. Weak impurity scattering ($τ \gg β^{–1}$) leaves the low-energy spectral density unaffected but leads to an interaction-driven enhancement of superconductivity. Our results in this work shed light on the role played by order-parameter fluctuations in providing the key missing link between Mott physics and strongly coupled toy models exhibiting gravity duals.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗