Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “independent particle approximation”

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 19 records

Comparison of interband related optical transitions and excitons in ZnGeN 2 and GaN

The optical dielectric function of ZnGeN 2 is calculated from the interband transitions using the energy bands calculated in the quasiparticle self-consistent (QS)GW^ method using two different levels of approximation: the independent particle approximation (IPA) and the Bethe-Salpeter Equation (BSE) approach. The first allows us to relate peaks in ε 2 to specific bands and k-points but does not include electron-hole interaction effects. The second includes electron-hole interaction or excitonic effects. The corresponding changes in the shape of ε 2 (ω) are found to be similar to those in GaN. The screened Coulomb interaction W^ is here calculated including electron-hole interactions in the polarization function and gives a band structure already going beyond the random phase approximation. The static dielectric constants including only electronic screening, commonly referred to as ε∞ , were calculated separately by extrapolating the wave vector dependent macroscopic dielectric function in the long-wavelength limit. Below the quasiparticle gap, we find three bound excitons optically active for different polarization. The convergence of these bound excitons with respect to the density of the k-mesh used in the BSE is studied and found to require a fine mesh. It is also found that these bound excitons originate from only the lowest conduction band and the top three valence bands. To incorporate the lattice screening, we include a scaling factor (ε∞/ε 0 ) 2 , which allows us to obtain exciton binding energies of the correct order of magnitude similar to those in GaN. Here, the excitons are related to each of the three fold split valence bands and the splittings of the latter are also studied as function of strain. Finally, a relation between the anisotropic effective masses and the valence band splitting is pointed out and explained.

36 MATERIALS SCIENCE↗

Large second harmonic generation in a penta-CdO 2 sheet exfoliated from its bulk phase

Two-dimensional (2D) materials composed solely of pentagonal motifs are of particular interest due to their unique geometries and novel properties. Especially, the broken centrosymmetry and in-plane mirror symmetry in penta-graphene-like materials result in both in-plane and out-of-plane second harmonic generation (SHG). Transition metal cations with a filled d 10 shell (d 10 -TM) can significantly enhance the SHG response as found in many previously studied bulk materials. Here, in this study, based on first-principles calculations combined with independent particle approximation, we show that large SHG can also exist in some 2D d 10 -TM oxides, such as penta-CdO 2 , which can be chemically exfoliated from its bulk phase with its dynamical, thermal, and mechanical stability intact. We further find that penta-CdO 2 possesses extraordinary in-plane and out-of-plane SHG responses with large static SHG susceptibilities of χ sheet 14 (0) = χ sheet 25 (0) = χ sheet 36 (0) = 8.86 pm 2 V -1 due to the phase-matching between fundamental and second-harmonic light. In addition, the low thermal expansion and a large optical band gap of 3.24 eV endow the penta-CdO 2 sheet with a high laser-induced damage threshold. Discussions are also made on the penta-ZnO 2 sheet. This study expands the family of 2D materials with outstanding SHG performance.

36 MATERIALS SCIENCE↗

Optical activity of solids from first principles

Within the framework of independent particle approximation, the optical activity tensor of solids is formulated as from different contributions: the magnetic dipole, electric quadrupole, and band dispersion terms. The first two terms have similar counterparts in the theory of finite systems, whereas the last term is unique for crystals. Further, the magnetic dipole and electric quadrupole transition moments are calculated with a sum-over-states formulation. We apply the formulation to calculate and analyze the optical rotation of elemental tellurium and the circular dichroism of the (6,4) carbon nanotube. Decomposed optical activity into different contributions is discussed. The calculated spectra agree well with experiments. As a showcase of achiral crystals, we calculate the optical activity of wurtzite GaN.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

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

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

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Comparison of Electron Capture Rates in the N = 50 Region using 1D Simulations of Core-collapse Supernovae

Abstract Recent studies have highlighted the sensitivity of core-collapse supernovae (CCSNe) models to electron-capture (EC) rates on neutron-rich nuclei near the N = 50 closed-shell region. In this work, we perform a large suite of one-dimensional CCSN simulations for 200 stellar progenitors using recently updated EC rates in this region. For comparison, we repeat the simulations using two previous implementations of EC rates: a microphysical library with parametrized N = 50 rates (LMP), and an older independent-particle approximation (IPA). We follow the simulations through shock revival up to several seconds post-bounce, and show that the EC rates produce a consistent imprint on CCSN properties, often surpassing the role of the progenitor itself. Notable impacts include the timescale of core collapse, the electron fraction and mass of the inner core at bounce, the accretion rate through the shock, the success or failure of revival, and the properties of the central compact remnant. We also compare the observable neutrino signal of the neutronization burst in a DUNE-like detector, and find consistent impacts on the counts and mean energies. Overall, the updated rates result in properties that are intermediate between LMP and IPA, and yet slightly more favorable to explosion than both.

79 ASTRONOMY AND ASTROPHYSICS↗

Inclusive inelastic scattering of heavy ions in the independent particle model

We consider the inclusive inelastic scattering of heavy ions using the Glauber (1959) model and the independent particle approximation. Inclusive inelastic distributions for projectile excitation of the target and total inelastic scattering, where all projectile and target excited states are summed, are discussed using closure. The total inelastic distribution, when integrated, is shown to be equivalent to the absorption cross section, found from applying the optical theorem to the elastic scattering amplitude in the coherent approximation. Calculations are presented for several heavy-ion pairs, using realistic nuclear densities in a large mass number approximation.

Cucinotta, Francis A.↗

Simplified Ab Initio Molecular Dynamics-Based Raman Spectral Simulations

We describe a simplified approach to simulating Raman spectra from ab initio molecular dynamics (AIMD) calculations. The protocol relies on on-the-fly calculations of approximate molecular polarizabilities using the well-known sum over orbitals (as opposed to states) method. This approach bypasses the more accurate but computationally expensive approach to calculating molecular polarizabilities along AIMD trajectories, i.e., solving the coupled perturbed Hartree–Fock/Kohn–Sham equations. We demonstrate the advantages and limitations of our method through a few case studies targeting molecular systems of interest to surface- and/or tip-enhanced Raman spectroscopy practitioners.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

A test of the applicability of independent scattering to high albedo planetary regoliths

We show that 'independent scattering' is a useful approximation for high albedo particles whose size and packing density are similar to typical particles in the lunar regolith. Laboratory measurements of the intensity and linear polarization of light scattered from a laboratory sample of glass spheres of know size and composition are compared to radiative transfer calculations of the same observable quantities. Mie scattering is integrated over the size distribution of the particles to determine the mean phase and polarization functions, or phase matrix, of the particle. Assuming that the particles scatter independently, the 'doubling method' is used to rigorously calculate multiple scattering in an optically thick layer of these anisotropically scattering particles. All of the major features 'predicted' in the calculated intensity (double peaks at small phase angles) and polarization (negative branch at small phase angles, large positive peak near 20 degrees phase, and small polarization elsewhere) are observed in the laboratory measurements, with good quantitative agreement indicated at phase angles less than 90 degrees. Even though the particles are supported by physical contact with each other in the sample, as are the particles in planetary regoliths, the independent scattering calculation yields a good approximation to both the intensity and polarization. The physical parameters input to the calculation are only the size distribution of the particles and their complex index of refraction (composition). Significant advantages of this approach are that the phase matrix is calculated from basic physical principles and both the intensity and polarization are determined simultaneously. This model may have broad applications to the interpretation of photometry, spectroscopy, and polarimetry of the icy regoliths of high albedo satellites. The intent of this effort is to perform a controlled experiment that tests the utility of the independent scattering assumption

Goguen, J. D.↗

Probing late-stage hadronic interactions at high baryon density via $K^{⁎0}$ production in the RHIC beam energy scan program

A precision measurement of the $K^{⁎0}$ meson yield is reported in Au+Au collisions at $\sqrt{s_{NN}}$ = 7.7, 11.5, 14.6, 19.6, and 27 GeV using the high-statistics data sample collected by the STAR experiment during the Beam Energy Scan II (BES-II) program at RHIC. The transeverse momentum (p T )-integrated yield ratios $\large{(}K^{⁎0} + \overline{K^{⁎0}}\large{)}/(K^+ + K^{-})$ in central collisions show a suppression relative to peripheral collisions at the (1.7–3.6) σ level, while a thermal model without final-stage rescattering overpredicts this ratio with a deviation of (6.9–8.2) σ. These results indicate a loss of the measured $K^{⁎0}$ signal in central collisions due to re-scattering of its hadronic decay products in the hadronic phase. The p T -integrated yield of charged kaons exhibits an approximate scaling with charged-particle multiplicity, independent of collision energy and system size. A similar trend is observed for the short-lived $K^{⁎0}$ resonance, although significant deviations emerge at lower energies. At BES energies, the $K^{⁎0}/K$ ratio shows stronger suppression than at the highest RHIC and LHC energies within a given multiplicity bin, particularly in central and mid-central collisions. This behavior is consistent with changes in the effective hadronic interaction cross section and is supported by transport model calculations, which indicate dominant meson–baryon interactions at lower energies and meson–meson interactions at higher energies.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Quenching of Single-Particle Strength in A=15 Nuclei

Here, absolute cross sections for the addition of s- and d-wave neutrons to 14 C and 14 N have been determined simultaneously via the (d,p) reaction at 10 MeV/u. The difference between the neutron and proton separation energies, ΔS, is around -20 MeV for the 14 C+n system and +8 MeV for 14 N+n. The population of the 1s 1/2 and 0d 5/2 orbitals for both systems is reduced by a factor of approximately 0.5 compared with the independent single-particle model, or about 0.6 when compared with the shell model. This finding strongly contrasts with results deduced from intermediate-energy knockout reactions between similar nuclei on targets of 9 Be and 12 C. The simultaneous technique used removes many systematic uncertainties.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Asymmetry parameters of the phase function for densely packed scattering grains

Spatial correlation among densely packed particles can substantially change their single-scattering properties, thus making questionable the applicability of the independent scattering approximation in calculations of light scattering by planetary regoliths. The same problem arises in geophysics in light scattering computations for snow, frosts, and bare soil. In this paper, we use a dense-medium light-scattering theory based on the introduction of the static structure factor to calculate asymmetry parameters of the phase function for densely packed particles with real refractive indices 1.31 and 1.66, approximating water ice and soil particles, respectively, and imaginary refractive indices 0, 0.01, and 0.3. For sparsely distributed, independently scattering grains, the calculated asymmetry parameters are always positive and always larger than those for densely packed particles. For densely packed grains, the asymmetry parameters may be negative but only for radius-to-wavelength ratios from about 0.1 to about 0.4. With decreasing particle size, the calculated asymmetry parameters tend to zero independently of the compaction state. In the geometrical optics regime, the asymmetry parameters for densely packed scatterers are positive and very close to those for independently scattering grains. These results may have important implications for remote sensing of the Earth and solid planetary surfaces. In particular, it is demonstrated that negative asymmetry parameters derived with some approximate multiple-scattering theories may be physically irrelevant and can be the result of using an inaccurate bidirectional reflection function combined with the ill-conditionally of the inverse scattering problem.

Mishchenko, Michael I.↗

An Overview of Methods for Deriving the Radiative Transfer Theory from the Maxwell Equations. II: Approach Based on the Dyson and Bethe-Salpeter Equations

In this paper, the vector radiative transfer equation is derived by means of the vector integral Foldy equations describing the electromagnetic scattering by a group of particles. By assuming that in a discrete random medium the positions of the particles are statistically independent and by applying the Twersky approximation to the order-of-scattering expansion of the total field, we derive the Dyson equation for the coherent field and the ladder approximated Bethe–Salpeter equation for the dyadic correlation function. Then, under the far-field assumption for sparsely distributed particles, the Dyson equation is reduced to the Foldy integral equation for the coherent field, while the iterated solution of the Bethe–Salpeter equation ultimately yields the vector radiative transfer equation.

Electromagnetic scattering↗

Intermediate scattering functions of a rigid body monoclonal antibody protein in solution studied by dissipative particle dynamic simulation

In the past decade, there was increased research interest in studying internal motions of flexible proteins in solution using Neutron Spin Echo (NSE) as NSE can simultaneously probe the dynamics at the length and time scales comparable to protein domain motions. However, the collective intermediate scattering function (ISF) measured by NSE has the contributions from translational, rotational, and internal motions, which are rather complicated to be separated. Widely used NSE theories to interpret experimental data usually assume that the translational and rotational motions of a rigid particle are decoupled and independent to each other. To evaluate the accuracy of this approximation for monoclonal antibody (mAb) proteins in solution, dissipative particle dynamic computer simulation is used here to simulate a rigid-body mAb for up to about 200 ns. The total ISF together with the ISFs due to only the translational and rotational motions as well as their corresponding effective diffusion coefficients is calculated. The aforementioned approximation introduces appreciable errors to the calculated effective diffusion coefficients and the ISFs. For the effective diffusion coefficient, the error introduced by this approximation can be as large as about 10% even though the overall agreement is considered reasonable. Thus, we need to be cautious when interpreting the data with a small signal change. In addition, the accuracy of the calculated ISFs due to the finite computer simulation time is also discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Cosmic ray anisotropies late in a solar flare event

The detailed relationship between the anisotropy characteristics observed during late times in the decay of a solar flare event and the interplanetary magnetic field parameters is investigated. The anisotropy always is from 45 deg east of the earth-sun line. This direction is approximately perpendicular to the nominal Archimedean spiral, independent of the particle energy. The amplitude of the anisotropy increases as the magnetic field azimuthal direction shows greater departure from the radial direction. These results are discussed in terms of current ideas about solar particle propagation in the interplanetary space.

Allum, F. R.↗

Diffuse reflectance of the ocean - The theory of its augmentation by chlorophyll a fluorescence at 685 nm

The radiative transfer equation is modified to include the effect of fluorescent substances and solved in the quasi-single scattering approximation for a homogeneous ocean containing fluorescent particles with wavelength independent quantum efficiency and a Gaussian shaped emission line. The results are applied to the in vivo fluorescence of chlorophyll a (in phytoplankton) in the ocean to determine if the observed quantum efficiencies are large enough to explain the enhancement of the ocean's diffuse reflectance near 685 nm in chlorophyll rich waters without resorting to anomalous dispersion. The computations indicate that the required efficiencies are sufficiently low to account completely for the enhanced reflectance. The validity of the theory is further demonstrated by deriving values for the upwelling irradiance attenuation coefficient at 685 nm which are in close agreement with the observations.

Gordon, H. R.↗

Singly-ionized helium in the driver gas of an interplanetary shock wave

The interplanetary shock wave observed on Jan. 29, 1977 by the HELIOS-1 plasma instruments shows an unusual feature: in the cold tenuous piston plasma following this shock, there appears a third peak in the energy per charge (E/q) spectra, in addition to the normal proton and alpha-particle peaks. The peak is located at E/q ? 4 and persists for about 14 hours, with slowly varying intensities. Independent simultaneous measurement of these particles' charge yields a value of approximately 1. These ions are thought to be He-4(+) ions travelling with the same speed as protons and alpha particles. The occurrence of He-4(+) indicates the possibility that during eruptive prominences or other solar transients, 'cold' chromospheric plasma might escape from the sun without undergoing the normal coronal heating process.

Schwenn, R.↗

Calculations of the variability of ice cloud radiative properties at selected solar wavelengths

This study shows that there is surprising little difference in values of reflectance, absorptance, and transmittance for many of the intermediate-size particle spectra. Particle size distributions with mode radii ranging from approximately 50 to 300 microns, irrespective of particle shape and nearly independent of the choice of size distribution representation, give relatively similar flux values. The very small particle sizes, however, have significantly larger values of reflectance and transmittance with corresponding smaller values of absorptance than do the larger particle sizes. The very large particle modes produce very small values of reflectance and transmittance along with very large values of absorptance. Such variations are particularly noticeable when plotted as a function of wavelength.

Welch, R. M.↗

A particle model for magnetotail neutral sheet equilibria

A particle model of the magnetotail neutral sheet self-consistent structure was developed by assuming that the thickness of the neutral sheet is much less than the thickness of the plasma sheet. This made it possible to approximate the neutral sheet as an x-independent structure. The role of nonlinear particle motions in determining the conductivity and the global current sheet structure is investigated.

Burkhart, G. R.↗