Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “EXPECTATION VALUE”

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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Investigation of Ga 2 O 3 as a new transparent conductive oxide for photovoltaics applications

This small innovation project intended to leverage recent activity and advances in ultrawide-bandgap oxide power electronics toward the investigation of Ga 2 O 3 as a transparent conductive oxide (TCO) for use in photovoltaics (PV). Ga 2 O 3 ’s theoretical advantage over incumbent TCOs is its very large bandgap of 4.8 eV, ensuring optical transparency of photons with λ ≥ 260 nm, effectively the full terrestrial solar spectrum. At the time of proposal writing, literature on Ga 2 O 3 as a TCO/optical material was relatively sparse — what little existed was focused on for UV sensors and/or emitters, with none related to PV. As such, this project was intended to help determine the potential of Ga 2 O 3 as a PV-oriented TCO by investigating its deposition using tools common to the PV and electronics industries — atomic layer deposition (ALD) and RF sputtering (RS), both of which are already used in PV manufacturing — and its resultant optical and electronic properties. Final project goals were to test its application to Si and III-V solar cells. Controllable film thickness and excellent uniformity was established for both methods deposition methods, with ALD providing higher precision for thinner films and RS more effective for growing thicker films. The as-deposited films were found to be amorphous in nature. Spectroscopic ellipsometry (SE) confirmed bandgaps of at least 4.8 eV and non-parasitic absorption across both AM0 and AM1.5G spectra; a refractive index approaching the expected value of 1.8 was also observed, with some degree of tunability based on process parameters. Using this initial optical data, a transfer matrix model, which interfaces with in-house EQE and LIV models, was developed to simulate Ga 2 O 3 optical effects on various solar cells, included potential in antireflection coatings (ARCs). Despite promising optical properties, the resultant resistivity / conductivity metrics did not meet expectations, regardless of deposition process. Although thick RS-deposited films, using both undoped and 1 at% Ge (n-type) doped sintered targets — demonstrated high net carrier concentrations (via C-V), and low specific resistance Ohmic contacts were demonstrated, transmission line measurements (TLM) showed very high resistivity. Further analysis indicated high vertical film conductivity, but very low lateral conductivity. Deeper characterization revealed a relatively high density of nano/microcrystalline inclusions that were ostensibly the source of the vertical conductivity, with the amorphous matrix serving as an effective insulator (likely due to high concentrations of electronic trap states). Despite extensive work to increase film polycrystallinity — demonstrated via high temperature annealing — sufficient lateral conductivity for TCO use was not achieved. The final phase of the project shifted focus toward investigation of ALD-deposited Ga 2 O 3 ’s potential as a passivant and/or passivating contact for both Si and III-V solar cells, with expectation of performance similar to Al 2 O 3 . However, initial rounds of testing using our baseline ALD process yielded minimal (but not quite zero) passivation of both GaInP and Si surfaces. Dielectric passivation is well-known to be highly process sensitive; given the unoptimized nature of the process used, this work was deemed inconclusive. The final outlook with respect to feasibility of Ga 2 O 3 as a PV-oriented TCO is, ultimately, still uncertain. The work performed in this project confirmed the optical properties and deposition methods, but the achieved electrical properties do not meet technological needs; literature reports in recent years are somewhat inconsistent and potentially untrustworthy, but generally appear to be in line with our results. At the very least it is clear that, should Ga 2 O 3 still be under consideration for PV TCO and/or selective contact applications, a significant amount of optimization and study is still needed.

14 SOLAR ENERGY↗
DOE OSTI2024

Dynamics of dark matter misalignment through the Higgs portal

A light singlet scalar field feebly coupled through the super-renormalizable Higgs portal provides a minimal and well-motivated realization of ultra-light bosonic dark matter. We study the cosmological production of dark matter in this model by elucidating the dynamics of two sources of scalar field misalignment generated during the radiation era. For large scalar masses (above $\mathcal{O}$ (10 −3 eV)), dark matter is produced through thermal misalignment, by which the scalar field is driven towards large field values as a result of the finite-temperature effective potential. The dominance of thermal misalignment in this mass range leads to a sharp relic abundance prediction which is, to a significant extent, insensitive to the initial conditions of the scalar field. On the other hand, for low mass scalars (below $\mathcal{O}$ (10 −5 eV)), dark matter is produced via VEV misalignment, which is caused by the induced scalar field vacuum expectation value triggered by the electroweak phase transition. We show that the relic abundance in this low mass range is sensitive to the scalar field initial conditions. In the intermediate mass range, the relic abundance is a consequence of a competition between thermal misalignment and VEV misalignment, which can potentially lead to novel forced resonance effects which cause a recurring enhancement and suppression in the late time oscillation amplitude as a function of the scalar mass. We compare our relic abundance predictions with constraints and projections from equivalence principle and inverse square law tests, stellar cooling, resonant molecular absorption, and observations of extra-galactic background light and diffuse X-ray backgrounds. New experimental ideas are needed to probe most of the cosmologically motivated regions of parameter space.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1987

IRAS observations of supernova remnants - A comparison between their infrared and X-ray cooling rates

A comparison is presented between the total IR and X-ray cooling rates of nine selected Galactic supernova remnants. The observed IR-to-X-ray cooling ratio (IRX ratio) values are larger than unity for most remnants, ranging from five for the adiabatic remnant Puppis A, to about 1000 for RCW 86. Most of the observed IR emission from the remnants can be attributed to thermal emission from dust collisionally heated by the shocked plasma. A comparison between the theoretical and observed IRX ratio shows that only two of the nine remnants have IRX ratios within a factor of about three of the expected value. Puppis A, Kepler, Tycho, and SN 1006 have IRX ratios that are significantly smaller than the theoretically predicted value, suggesting that the dust is significantly depleted in the ambient medium into which they are expanding.

Dwek, Eli↗
NASA NTRS1997

Advection in Accretion Disk Boundary Layers

Recent numerical simulations of accretion disk boundary layers have shown qualitatively the importance tance of advected energy in the inner region of the disk. In this short paper we present quantitative results of advective boundary layers in the optically thick regime. Numerical results are obtained for various systems, by means of a one-dimensional time-dependent numerical code. At high accretion mass rates, dot-M approx. = 10(exp -4) solar mass/yr, or low values of the viscosity parameter, alpha approx. = 0.001-0.01 (characteristic of FU Orionis systems and some symbiotic stars), the optical thickness in the inner part of the disk becomes very large (tau much greater than 1). The disk, unable to cool efficiently, becomes geometrically thick (H/r approx. = 0.5). The energy dissipated in the dynamical boundary layer is radiated outward to larger radii and advected into the star. The boundary-layer luminosity is only a fraction of its expected value; the rest of the energy is advected into the star. The fraction of the advected energy is zeta = L(sub adv)/L(sub acc) approx. = 0.1 in symbiotic stars (accretion onto a low-mass main-sequence star) and zeta approx. = 0.2 in FU Ori systems (accretion onto a pre-main sequence star).

Godon, Patrick↗
DOE OSTI2020

Spin-1 Kitaev-Heisenberg model on a honeycomb lattice

We study the Kitaev-Heisenberg model with spin-1 local degree of freedom on a honeycomb lattice numerically by the infinite density matrix renormalization group method on a cylinder geometry. By tuning the relative value of the Kitaev and Heisenberg exchange couplings, we obtain the phase diagram with two spin liquid phases and four symmetry-broken phases. Here, we identify that the spin liquid phases are gapless by calculating the central charge at the pure Kitaev points without the Heisenberg interactions. Comparing to its spin-1/2 counterpart, the position and number of gapless modes of the spin-1 case are quite different. Due to the approximate Z 2 local conservations, the expectation value of the Wilson loop operator measuring the flux of each plaquette stays near to 1, and the static spin-spin correlations remain short range in the entire spin liquid phases.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Atomistic Simulations of the Elastic Compression of Platinum Nanoparticles

Abstract The elastic behavior of nanoparticles depends strongly on particle shape, size, and crystallographic orientation. Many prior investigations have characterized the elastic modulus of nanoscale particles using experiments or simulations; however their reported values vary widely depending on the methods for measurement and calculation. To understand these discrepancies, we used classical molecular dynamics simulation to model the compression of platinum nanoparticles with two different polyhedral shapes and a range of sizes from 4 to 20 nm, loaded in two different crystal orientations. Multiple standard methods were used to calculate the elastic modulus from stress-vs-strain data for each nanoparticle. The magnitudes and particle-size dependence of the resulting moduli varied with calculation method and, even for larger nanoparticles where bulk-like behavior may be expected, the effective elastic modulus depended strongly on shape and orientation. Analysis of per-atom stress distributions indicated that the shape- and orientation-dependence arise due to stress triaxiality and inhomogeneity across the particle. When the effective elastic modulus was recalculated using a representative volume element in the center of a large nanoparticle, the elastic modulus had the expected value for each orientation and was shape independent. It is only for single-digit nanoparticles that meaningful differences emerged, where even the very center of the particle had a lower modulus due to the effect of the surface. These findings provide better understanding of the elastic properties of nanoparticles and disentangle geometric contributions (such as stress triaxiality and spatial inhomogeneity) from true changes in elastic properties of the nanoscale material.

36 MATERIALS SCIENCE↗
NASA NTRS2003

Interleaved Observation Execution and Rescheduling on Earth Observing Systems

Observation scheduling for Earth orbiting satellites solves the following problem: given a set of requests for images of the Earth, a set of instruments for acquiring those images distributed on a collecting of orbiting satellites, and a set of temporal and resource constraints, generate a set of assignments of instruments and viewing times to those requests that satisfy those constraints. Observation scheduling is often construed as a constrained optimization problem with the objective of maximizing the overall utility of the science data acquired. The utility of an image is typically based on the intrinsic importance of acquiring it (for example, its importance in meeting a mission or science campaign objective) as well as the expected value of the data given current viewing conditions (for example, if the image is occluded by clouds, its value is usually diminished). Currently, science observation scheduling for Earth Observing Systems is done on the ground, for periods covering a day or more. Schedules are uplinked to the satellites and are executed rigorously. An alternative to this scenario is to do some of the decision-making about what images are to be acquired on-board. The principal argument for this capability is that the desirability of making an observation can change dynamically, because of changes in meteorological conditions (e.g. cloud cover), unforeseen events such as fires, floods, or volcanic eruptions, or un-expected changes in satellite or ground station capability. Furthermore, since satellites can only communicate with the ground between 5% to 10% of the time, it may be infeasible to make the desired changes to the schedule on the ground, and uplink the revisions in time for the on-board system to execute them. Examples of scenarios that motivate an on-board capability for revising schedules include the following. First, if a desired visual scene is completely obscured by clouds, then there is little point in taking it. In this case, satellite resources, such as power and storage space can be better utilized taking another image that is higher quality. Second, if an unexpected but important event occurs (such as a fire, flood, or volcanic eruption), there may be good reason to take images of it, instead of expending satellite resources on some of the lower priority scheduled observations. Finally, if there is unexpected loss of capability, it may be impossible to carry out the schedule of planned observations. For example, if a ground station goes down temporarily, a satellite may not be able to free up enough storage space to continue with the remaining schedule of observations. This paper describes an approach for interleaving execution of observation schedules with dynamic schedule revision based on changes to the expected utility of the acquired images. We describe the problem in detail, formulate an algorithm for interleaving schedule revision and execution, and discuss refinements to the algorithm based on the need for search efficiency. We summarize with a brief discussion of the tests performed on the system.

Khatib, Lina↗
DOE OSTI2021

Regulation of the normalized rate of driven magnetic reconnection through shocked flux pileup

Magnetic reconnection is explored on the Terrestrial Reconnection Experiment (TREX) for asymmetric inflow conditions and in a configuration where the absolute rate of reconnection is set by an external drive. Magnetic pileup enhances the upstream magnetic field of the high-density inflow, leading to an increased upstream Alfvén speed and helping to lower the normalized reconnection rate to values expected from theoretical consideration. In addition, a shock interface between the far upstream supersonic plasma inflow and the region of magnetic flux pileup is observed, important to the overall force balance of the system, thereby demonstrating the role of shock formation for configurations including a supersonically driven inflow. Despite the specialized geometry where a strong reconnection drive is applied from only one side of the reconnection layer, previous numerical and theoretical results remain robust and are shown to accurately predict the normalized rate of reconnection for the range of system sizes considered. This experimental rate of reconnection is dependent on system size, reaching values as high as 0.8 at the smallest normalized system size applied.

Physics↗
DOE OSTI2023

Development of Multiresolution Capabilities for the Holistic Energy Resource Optimization Network (HERON) tool A progress update

INL researchers work on technoeconomic analyses for integrated energy systems (IES) using the Framework for Optimization of ResourCes and Economics (FORCE). Within FORCE, researchers use the Holistic Energy Resource Optimization Network (HERON) tool to conduct optimization of grid portfolios under uncertain market conditions. These optimizations determine optimal capacities for all IES components and strategies for resource dispatch which maximize some economic metric (e.g., net present value). Resource dispatch occurs on finer timescales (typically hours) and thus are asked to respond to a given time series (e.g. hourly load demand profiles for a grid, or pre-determined electricity prices). Volatile and complex bidding dynamics as well as poorly forecasted weather events within deregulated markets add uncertainty to the time series; FORCE can address this uncertainty by training a reduced order model on historical time series and generate unique synthetic time series which represent individual scenarios or realizations of the market. The IES configuration can be simulated under these different sampled realizations and a stochastic optimization is conducted which optimizes the expected value of the desired economic metric. The training of a synthetic time series generator is limited by the chosen time resolution; dynamics can occur on different time scales. Seasonal demand trends can dominate faster dynamical events (such as power outages from certain sectors or severe weather events) which might not get captured correctly by the trained model. In this report, we investigate different ways of addressing the training and generation of time series on multiple time scales using three main algorithms: wavelet decomposition, dynamic mode decomposition, and generative adversarial networks for time series. We demonstrate a time series analysis that yields information on not just the frequency space but also temporal space: where a fast Fourier transform can provide what frequencies dominate, the new algorithms can provide when the frequencies dominate as well. These analyses can help improve IES optimization by allowing researchers to couple simulations at different timescales when it is most needed - seasonal, day-ahead, and real time optimization - with greater computational efficiency. Future work will include implementation of a subset of the proposed algorithms into the FORCE toolset and application of these analyses into multiple timescale optimization.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

Proof that Combining the Forced-collision and DXTRAN Monte Carlo Variance-reduction Techniques is Fair

This report provides mathematical proof that the MCNP DXTRAN (also known in other codes as forced-flight) and forced-collision variance-reduction techniques do not bias the expected value of Monte Carlo simulation estimates when combined. Proof is also provided that the techniques are unbiased when used independently. To prove that the techniques are unbiased, this report derives the first History Score Moment Equation (HSME) for non-multiplying media employing only forced collisions, only DXTRAN, and then both the variance-reduction techniques combined, as defined next. The HSMEs are found by forming the History Score Probability Density Functions (HSPDFs) and then taking the first score moment. Following its derivation, each HSME with a variance-reduction technique employed is reduced to the HSME for an analog simulation. Because the HSME represents the expected contribution to estimators in the simulation, reducing the HSME with variance reduction to the analog HSME shows that the simulation is unbiased despite the variance-reduction technique considered in the HSME. Analysis regarding higher-score moments of each technique is the subject of prior work and is not addressed herein. Throughout this work, the phase space p is defined to be the particle position x, direction-of-flight unit vector $\hat{Ω}$, energy E, and statistical weight w, $$p ≡ (x; \hat{Ω}; E; w).$$ A reduced phase-space excluding the statistical weight of the particle, $$r ≡ (x; \hat{Ω}; E),$$ is also used.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Dynamic, symmetry-preserving, and hardware-adaptable circuits for quantum computing many-body states and correlators of the Anderson impurity model

We present a hardware-reconfigurable ansatz on N q -qubits for the variational preparation of many-body states of the Anderson impurity model (AIM) with N imp + N bath = N q /2 sites, which conserves total charge and spin z component within each variational search subspace. The many-body ground state of the AIM is determined as the minimum over all minima of O(N$^2_ q$) distinct charge-spin sectors. Hamiltonian expectation values are shown to require ω(N q ) < N meas. $\leqslant$ O(N imp N bath ) symmetry-preserving, parallelizable measurement circuits, each amenable to postselection. To obtain the one-particle impurity Green’s function we show how initial Krylov vectors can be computed via midcircuit measurement and how Lanczos iterations can be computed using the symmetry-preserving ansatz. For a single-impurity Anderson model with a number of bath sites increasing from one to seven, we show using numerical emulation that the ease of variational ground-state preparation is suggestive of linear scaling in circuit depth and subquartic scaling in optimizer complexity. We therefore expect that, combined with time-dependent methods for Green’s function computation, our ansatz provides a useful tool to account for electronic correlations on early fault-tolerant processors. Finally, with a view towards computing real materials properties of interest like magnetic susceptibilities and electron-hole propagators, we provide a straightforward method to compute many-body, time-dependent correlation functions using a combination of time evolution, midcircuit measurement-conditioned operations, and the Hadamard test.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Hybrid Analytics Solution to Improve Coal Power Plant Operations

This project focused on developing advanced methods for thermal performance monitoring of a coal-fueled power plant. The specific goal was to develop and demonstrate a new thermal performance monitoring approach using a hybrid model that integrates a physics-based heat balance model with a machine learning-based pattern recognition model. The hybrid model enables increased accuracy and scope of the thermal analysis and an improved ability to monitor and detect changes in plant operation. This new approach takes full advantage of the individual model capabilities and creates an important new set of capabilities not previously possible using the two types of models separately. Using the heat balance model, a rich set of derived parameters (virtual sensors) are calculated from the measured plant operating data at each time point. The combined measured and derived data values are used by machine learning algorithms to create pattern recognition models over the range of normal unit operation. To create the monitoring models, historical data from normal operation of the plant is first processed by the heat balance model to compute the derived parameter data. The result is a greatly expanded set of normal operating data that can be used as input to create the pattern recognition model. Once the models are calibrated for normal operation, the hybrid model is suitable for use in continuous online monitoring. During online monitoring, new plant operating data is processed first by the heat balance model and then by the pattern recognition model. Results from the pattern recognition model quantify the deviation of each measured or derived parameter from its expected value in normal operation. The hybrid models can detect abnormal changes in plant operating data with very high accuracy and sensitivity. When abnormal behavior is detected, alerts are generated automatically for evaluation by the plant monitoring staff. The new hybrid solution product was developed and verified in the performance of the project. The hybrid solution was tested first in a simulation environment that mimicked the plant data systems and infrastructure used by U.S. power generating plants and utilities. The hybrid solution was then deployed for real-time, online monitoring of an operating coal-fueled power plant at a field test site. Field testing demonstrated that all hybrid solution development objectives were accomplished. The project work was based on combining the capabilities of two existing software products to create the new hybrid solution product. One of these was the existing MapEx® heat balance product and the other was the existing SureSense® advanced pattern recognition product. Each of these separate products was assessed to be at a Technology Readiness Level (TRL) of 9 at the start of the effort. The hybrid solution product was assessed to be at a TRL of 2 at the start of the project based on early feasibility work by the project team. At completion of the field testing performed in the project, the hybrid solution product was assessed to be at a TRL of 7. The project team expects that the hybrid solution product will be deployed commercially and will achieve a TRL of 9 within one year after completion of the project.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2022

A -type antiferromagnetic order in the Zintl-phase insulator EuZn 2 P 2

Zintl phases, containing strongly covalently bonded frameworks with separate ionically bonded ions, have emerged as a critical materials family in which to couple magnetism and strong spin-orbit coupling to drive diverse topological phases of matter. Here we report the single-crystal synthesis, magnetic, thermodynamic, transport, and theoretical properties of the Zintl compound EuZn 2 P 2 that crystallizes in the anti-La 2 O 3 (CaAl 2 Si 2 ) P-3m1 structure, containing triangular layers of Eu 2+ ions. In-plane resistivity measurements reveal insulating behavior with an estimated activation energy of E g = 0.11eV. Specific heat and magnetization measurements indicate antiferromagnetic ordering at T N = 23K. Curie-Weiss analysis of in-plane and out of plane magnetic susceptibility from T = 150 to 300 K yields p eff = 8.61 for μ 0 H⊥c and p eff = 7.74 for μ 0 H//c, close to the expected values for the 4f 7 J = S = 7/2 Eu 2+ ion and indicative of weak anisotropy. Below T N , a significant anisotropy of χ ⊥ /χ // ≈ 2.3 develops, consistent with A-type magnetic order as observed in isostructural analogs and as predicted by the density functional theory calculations reported herein. The positive Weiss temperatures of θ W =19.2K for μ 0 H⊥c and θ W =41.9K for μ 0 H//c show a similar anisotropy and suggest competing ferromagnetic and antiferromagnetic interactions. Comparing Eu magnetic ordering temperatures across trigonal EuM 2 X 2 (M= divalent metal, X= pnictide) shows that EuZn 2 P 2 exhibits the highest ordering temperature, with variations in T N correlating with changes in expected dipolar interaction strengths within and between layers and independent of the magnitude of electrical conductivity. These results provide experimental validation of the crystochemical intuition that the cation Eu 2+ layers and the anionic (M 2 X 2 ) 2– framework can be treated as electronically distinct subunits, enabling further predictive materials design.

36 MATERIALS SCIENCE↗
NASA NTRS1987

Heat transfer with very high free-stream turbulence and heat transfer with streamwise vortices

Two experimental programs related to augmentation of heat transfer by complex flow characteristics are reviewed. The first program deals with very high turbulence (up to 63 percent) which was shown to result in Stanton numbers as much as five times the expected values. Results from a number of trials show that fixing the free stream velocity, x-Reynolds number, turbulence intensity and integral length scale does not fix the Stanton number. Two such cases were found in which the Stanton number of one was 40 percent larger than the other. Mean velocity and mean temperature profiles are presented, as well as profiles of turbulence intensity within the boundary layer. The second program deals with vortices originating at bluff bodies and traveling downstream embedded in the wall boundary layer. Velocity vector maps from the boundary layers and distributions of Stanton number on the wall are presented for three types of bodies: square, cylindrical and teardrop. The heat transfer and velocity maps do not show evidence of the expected horseshoe vortices but, instead, show a strong common flow up vortex pair. The fluid mechanic mechanism responsible for this secondary flow field has not yet been identified.

Moffat, Robert J.↗
DOE OSTI2022

Quantum algorithm for electronic band structures with local tight-binding orbitals

While the main thrust of quantum computing research in materials science is to accurately measure the classically intractable electron correlation effects due to Coulomb repulsion, designing optimal quantum algorithms for simpler problems with well-understood solutions is a useful tactic to advance our quantum “toolbox”. With this in mind, we consider the quantum calculation of a periodic system’s single-electron band structure over a path through reciprocal space. Previous efforts have used the Variational Quantum Eigensolver algorithm to solve the energy of each band, which involves numerically optimizing the parameters of a variational quantum circuit to minimize a cost function, constructed as the expectation value of a Hamiltonian operator. Traditionally, a unique Hamiltonian operator is constructed for each k-point, so that many cost functions, each with their own parameter space, must be optimized to generate a single band. Similarly, calculating higher bands than the first has traditionally involved modifying the cost function with additional overlap terms to ensure higher-energy eigenstates are orthogonal to those of lower bands. In this paper, we adopt a direct space approach, using a novel hybrid first/second-quantized qubit mapping which allows us to construct a single Hamiltonian, and a single cost-function, suitable for solving the entire electronic band structure. In contrast to previous approaches, the k-point and the band index are selected by additional parameters in our quantum circuit, rather than through modifications to the cost function. The result is a technically and conceptually simpler approach to band structure calculations on a quantum computer. Moreover, we expect that the tools developed herein will motivate new strategies for tackling highly-correlated materials beyond the grasp of classical computing.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1999

ISCCP Cloud Properties Associated with Standard Cloud Types Identified in Individual Surface Observations

Individual surface weather observations from land stations and ships are compared with individual cloud retrievals of the International Satellite Cloud Climatology Project (ISCCP), Stage C1, for an 8-year period (1983-1991) to relate cloud optical thicknesses and cloud-top pressures obtained from satellite data to the standard cloud types reported in visual observations from the surface. Each surface report is matched to the corresponding ISCCP-C1 report for the time of observation for the 280x280-km grid-box containing that observation. Classes of the surface reports are identified in which a particular cloud type was reported present, either alone or in combination with other clouds. For each class, cloud amounts from both surface and C1 data, base heights from surface data, and the frequency-distributions of cloud-top pressure (p(sub c) and optical thickness (tau) from C1 data are averaged over 15-degree latitude zones, for land and ocean separately, for 3-month seasons. The frequency distribution of p(sub c) and tau is plotted for each of the surface-defined cloud types occurring both alone and with other clouds. The average cloud-top pressures within a grid-box do not always correspond well with values expected for a reported cloud type, particularly for the higher clouds Ci, Ac, and Cb. In many cases this is because the satellites also detect clouds within the grid-box that are outside the field of view of the surface observer. The highest average cloud tops are found for the most extensive cloud type, Ns, averaging 7 km globally and reaching 9 km in the ITCZ. Ns also has the greatest average retrieved optical thickness, tau approximately equal 20. Cumulonimbus clouds may actually attain far greater heights and depths, but do not fill the grid-box. The tau-p(sub c) distributions show features that distinguish the high, middle, and low clouds reported by the surface observers. However, the distribution patterns for the individual low cloud types (Cu, Sc, St) occurring alone overlap to such an extent that it is not possible to distinguish these cloud types from each other on the basis of tau-p(sub c) values alone. Other cloud types whose tau-p(sub c) distributions are indistinguishable are Cb, Ns, and thick As. However, the tau-p(sub c) distribution patterns for the different low cloud types are nevertheless distinguishable when all occurrences of a low cloud type are included, indicating that the different low types differ in their probabilities of co-occurrence with middle and high clouds.

Hahn, Carole J.↗
NASA NTRS1993

Low redshift Lyman alpha absorption lines and the dark matter halos of disk galaxies

Recent observations using the Hubble Space Telescope of the z = 0.156 QSO 3C 273 have discovered a surprisingly large number of Ly-alpha absorption lines. In particular, Morris et al. found 9 certain and 7 possible Ly-alpha lines with equivalent widths above 25 mA. This is much larger (by a factor of 5-10) than the number expected from extrapolation of the high-redshift behavior of the Ly-alpha forest. Within the context of pressure-confined models for the Ly-alpha clouds, this behavior can be understood if the ionizing background declines sharply between z is approximately 2 and z is approximately 0. However, this requires that the ionizing photon flux drop as rapidly as the QSO volume emissivity; moreover, the absorbers must have a space density n(sub O) is approximately 2.6(N/10)h/((D/100 kpc)(sup 2)) Mpc(sup -3) where D is the present-day diameter of the absorbers. It is somewhat surprising that such necessarily fragile objects could have survived in such numbers to the present day. It is shown that it is plausible that the atomic hydrogen extents of spiral and irregular galaxies are large enough to produce the observed number of Ly-alpha absorption lines toward 3C 273, and that the neutral column densities and doppler b-values expected under these conditions fall in the range found by Morris et al. (1991).

Maloney, Philip↗
NASA NTRS1995

Multiparameter linear least-squares fitting to Poisson data one count at a time

A standard problem in gamma-ray astronomy data analysis is the decomposition of a set of observed counts, described by Poisson statistics, according to a given multicomponent linear model, with underlying physical count rates or fluxes which are to be estimated from the data. Despite its conceptual simplicity, the linear least-squares (LLSQ) method for solving this problem has generally been limited to situations in which the number n(sub i) of counts in each bin i is not too small, conventionally more than 5-30. It seems to be widely believed that the failure of the LLSQ method for small counts is due to the failure of the Poisson distribution to be even approximately normal for small numbers. The cause is more accurately the strong anticorrelation between the data and the wieghts w(sub i) in the weighted LLSQ method when square root of n(sub i) instead of square root of bar-n(sub i) is used to approximate the uncertainties, sigma(sub i), in the data, where bar-n(sub i) = E(n(sub i)), the expected value of N(sub i). We show in an appendix that, avoiding this approximation, the correct equations for the Poisson LLSQ (PLLSQ) problems are actually identical to those for the maximum likelihood estimate using the exact Poisson distribution. We apply the method to solve a problem in high-resolution gamma-ray spectroscopy for the JPL High-Resolution Gamma-Ray Spectrometer flown on HEAO 3. Systematic error in subtracting the strong, highly variable background encountered in the low-energy gamma-ray region can be significantly reduced by closely pairing source and background data in short segments. Significant results can be built up by weighted averaging of the net fluxes obtained from the subtraction of many individual source/background pairs. Extension of the approach to complex situations, with multiple cosmic sources and realistic background parameterizations, requires a means of efficiently fitting to data from single scans in the narrow (approximately = 1.2 keV, HEAO 3) energy channels of a Ge spectrometer, where the expected number of counts obtained per scan may be very low. Such an analysis system is discussed and compared to the method previously used.

Wheaton, Wm. A.↗