Search NASA⌕ Search

SEARCH · Search NASA

Results for “Complex reflection coefficients”

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.

157 records · Page 9
All sourcesNASA NTRSDOE OSTI
NASA NTRS2003

Plumes and Blooms: Observations, Analysis and Modeling for SIMBIOS

The goal of the Plumes and Blooms (PnB) project is to develop, validate and apply to imagery state-of-theart ocean color algorithms for quantifying sediment plumes and phytoplankton blooms for the Case I1 environment of the Santa Barbara Channel. We conduct monthly to twice-monthly transect observations across the Santa Barbara Channel to develop an algorithm development and product validation data set. The PnB field program started in the summer of 1996. At each of the 7 PnB stations, a complete verification bio-geo-optical data set is collected. Included are redundant measures of apparent optical properties (remote sensing reflectance and diffuse attenuation spectra), as well as in situ profiles of spectral absorption, beam attenuation and backscattering coefficients. Water samples are analyzed for component in vivo absorption spectra, fluorometric chlorophyll, phytoplankton pigment (by the SDSU CHORS laboratory), and inorganic nutrient concentrations (Table 1). A primary goal is to use the PnB field data set to objectively tune semi-analytical models of ocean color for this site and apply them using available satellite imagery (SeaWiFS and MODIS). In support of this goal, we have also been addressing SeaWiFS ocean color and AVHRR SST imagery (Otero and Siegel, 2003). We also are using the PnB data set to address time/space variability of water masses in the Santa Barbara Channel and its relationship to the 1997/1998 El Niiio. However, the comparison between PnB field observations and satellite estimates of primary products has been disappointing. We find that field estimates of water-leaving radiance, LwN(h), correspond poorly to satellite estimates for both SeaWiFS and MODIS local area coverage imagery. We believe this is due to poor atmospheric correction due to complex mixtures of aerosol types found in these near-coastal regions. Last, we remain active in outreach activities.

Siegel, D. A.↗
NASA NTRS2014

A Flexible Parameterization for Shortwave Optical Properties of Ice Crystals

A parameterization is presented that provides extinction cross section sigma (sub e), single-scattering albedo omega, and asymmetry parameter (g) of ice crystals for any combination of volume, projected area, aspect ratio, and crystal distortion at any wavelength in the shortwave. Similar to previous parameterizations, the scheme makes use of geometric optics approximations and the observation that optical properties of complex, aggregated ice crystals can be well approximated by those of single hexagonal crystals with varying size, aspect ratio, and distortion levels. In the standard geometric optics implementation used here, sigma (sub e) is always twice the particle projected area. It is shown that omega is largely determined by the newly defined absorption size parameter and the particle aspect ratio. These dependences are parameterized using a combination of exponential, lognormal, and polynomial functions. The variation of (g) with aspect ratio and crystal distortion is parameterized for one reference wavelength using a combination of several polynomials. The dependences of g on refractive index and omega are investigated and factors are determined to scale the parameterized (g) to provide values appropriate for other wavelengths. The parameterization scheme consists of only 88 coefficients. The scheme is tested for a large variety of hexagonal crystals in several wavelength bands from 0.2 to 4 micron, revealing absolute differences with reference calculations of omega and (g) that are both generally below 0.015. Over a large variety of cloud conditions, the resulting root-mean-squared differences with reference calculations of cloud reflectance, transmittance, and absorptance are 1.4%, 1.1%, and 3.4%, respectively. Some practical applications of the parameterization in atmospheric models are highlighted.

parameterization↗
DOE OSTI2015

Berkeley Lab WINDOW

WINDOW features include: - Microsoft Windows TM interface - algorithms for the calculation of total fenestration product U-values and Solar Heat Gain Coefficient consistent with ASHRAE SPC 142, ISO 15099, and the National Fenestration Rating Council - a Condensation Resistance Index in accordance with the NFRC 500 Standard - and integrated database of properties - imports data from other LBNL window analysis software: - Import THERM file into the Frame Library - Import records from IGDB and OPtics5 into the Glass Library for the optical properties of coated and uncoated glazings, laminates, and applied films. Program Capabilities WINDOW 7.2 offers the following features: The ability to analyze products made from any combination of glazing layers, gas layers, frames, spacers, and dividers under any environmental conditions and at any tilt; The ability to model complex glazing systems such as venetian blinds and roller shades. Directly accessible libraries of window system components, (glazing systems, glazing layers, gas fills, frame and divider elements), and environmental conditions; The choice of working in English (IP), or Systeme International (SI) units; The ability to specify the dimensions and thermal properties of each frame element (header, sills, jamb, mullion) in a window; A multi-band (wavelength-by-wavelength) spectral model; A Glass Library which can access spectral data files for many common glazing materials from the Optics5database; A night-sky radiative model; A link with the DOE-2.1E and Energy Plus building energy analysis program. Performance Indices and Other Results For a user-defined fenestration system and user-defined environmental conditions, WINDOW calculates: The U-value, solar heat gain coefficient, shading coefficient, and visible transmittance for the complete window system; The U-value, solar heat gain coefficient, shading coefficient, and visible transmittance for the glazing system (center-of-glass values); The U-values of the frame and divider elements and corresponding edge-of-glass areas (based on generic correlations); The total solar and visible transmittance and reflectances of the glazing system. Color properties, i.e. L*, a*, and b* color coordinates, dominant wavelength, and purity for transmitted and reflected (outdoor) solar radiation; The damage-weighted transmittance of the glazing system between 0.3 an 0.38 microns; The angular dependence of the solar and visible transmittances, solar and visible reflectances, solar absorptance, and solar heat gain coefficient of the glazing system; The percent relative humidity of the inside and outside air for which condensation will occur on the interior and exterior glazing surfaces respectively; The center-of-glass temperature distribution.

Curcija, DraganCharlie↗
DOE OSTI2018

Soft x-ray optical constants of sputtered chromium thin films with improved accuracy in the L and M absorption edge regions

In this paper, we determine with improved accuracy the complex index of refraction n = 1 - δ + iβ of sputtered chromium thin films for photon energies ranging from 25 eV to 813 eV. These data include the first absolute measurements of the absorption fine structure near the Cr-L edge. First, we verified by combining Rutherford Backscattering Spectrometry and grazing-incidence x-ray reflectometry that the sputtered thin films were pure Cr with a density consistent with tabulated values. Then, we demonstrated that the Cr surface oxide layer remains stable when the samples are exposed to air for up to 4 years. The Cr absorption coefficient β was determined from the transmittance of freestanding Cr thin films with various thicknesses, measured at the ALS synchrotron radiation source. A model is proposed to correct the transmittance data from the spectral contamination of the source. Finally, we used the new β values, combined with theoretical and tabulated data from the literature, in order to calculate the δ values by the Kramers-Kronig relation. The improvement in the accuracy of β values is demonstrated by the f-sum rule. Finally, an additional validation of the new Cr optical constants (δ, β) is performed by comparing the simulated and experimental reflectance of a Cr/B 4 C multilayer mirror near the Cr-L 2,3 edge.

36 MATERIALS SCIENCE↗
NASA NTRS2024

Exploring Standing and Reflected Slow-mode Waves in Flaring Coronal Loops: A Parametric Study Using 2.5D MHD Modeling

Abstract Recent observations of reflected propagating and standing slow-mode waves in hot flaring coronal loops have spurred our investigation into their underlying excitation and damping mechanisms. To understand these processes, we conduct 2.5D magnetohydrodynamic (MHD) simulations using an arcade active region model that includes a hot and dense loop. Our simulations allow for in-depth parametric investigations complementing and expanding our previous 3D MHD modeling results. We excite these waves using a large-amplitude, flow pulse applied at one footpoint of the loop in two distinct models as motivated by observations from the Solar Dynamics Observatory/Atmospheric Imaging Assembly (SDO/AIA). The first model (Model 1) incorporates classical compressive viscosity coefficient, while the second model (Model 2) adopts a 10-times enhanced viscosity coefficient. We obtain the following major results: (1) Our 2.5D MHD simulations reinforce previous conclusions derived from 1D and 3D MHD models that significantly enhanced viscosity is crucial for the rapid excitation of standing slow waves with damping times consistent with observations by Wang et al. (2015). (2)We uncover that nonlinearity in Model 1 delays the conversion of a reflected propagating wave into a standing wave. In contrast, Model 2 exhibits a much weak influence of nonlinearity on the excitation time of standing waves, thanks to the suppression of these effects by enhanced viscosity. (3) Our results reveal that the transverse temperature structure holds more influence on wave behavior than the density structure. In Model 1, increased loop temperature contrast significantly enhances wave trapping within the structure, mitigating the impact of temperature-dependent viscous damping. Conversely, in Model 2, the impact of temperature structure on wave behavior weakens in comparison to the effect of viscosity. (4) Model 1 displays evident nonlinear coupling to the fast and kink magnetoacoustic waves and pronounced wave leakage into the corona. Model 2 exhibits significantly weaker effects in this regard. Analyzing three observed wave events by SDO/AIA aligns with Model 2 predictions, providing further support for the substantial viscosity increase. Our 2.5D study unravels the complex interplay of wave-flow phenomena and nonlinear processes in coronal loops, extending our previous 1D modeling results to incorporate more realistic loop geometry. This provides insights into scenarios where 3D effects may be neglected, thereby enhancing our understanding of the intricate dynamics of the solar corona.

Tongjiang Wang↗
NASA NTRS2003

Development of Numerical Codes for Modeling Electromagnetic Behavior at High Frequencies Near Large Objects

A study into the problem of determining electromagnetic solutions at high frequencies for problems involving complex geometries, large sizes and multiple sources (e.g. antennas) has been initiated. Typical applications include the behavior of antennas (and radiators) installed on complex conducting structures (e.g. ships, aircrafts, etc..) with strong interactions between antennas, the radiation patterns, and electromagnetic signals is of great interest for electromagnetic compatibility control. This includes the overall performance evaluation and control of all on-board radiating systems, electromagnetic interference, and personnel radiation hazards. Electromagnetic computational capability exists at NASA LaRC, and many of the codes developed are based on the Moment Method (MM). However, the MM is computationally intensive, and this places a limit on the size of objects and structures that can be modeled. Here, two approaches are proposed: (i) a current-based hybrid scheme that combines the MM with Physical optics, and (ii) an Alternating Direction Implicit-Finite Difference Time Domain (ADI-FDTD) method. The essence of a hybrid technique is to split the overall scattering surface(s) into two regions: (a) a MM zone (MMZ) which can be used over any part of the given geometry, but is most essential over irregular and "non-smooth" geometries, and (b) a PO sub-region (POSR). Currents induced on the scattering and reflecting surfaces can then be computed in two ways depending on whether the region belonged to the MMZ or was part of the POSR. For the MMZ, the current calculations proceed in terms of basis functions with undetermined coefficients (as in the usual MM method), and the answer obtained by solving a system of linear equations. Over the POSR, conduction is obtained as a superposition of two contributions: (i) currents due to the incident magnetic field, and (ii) currents produced by the mutual induction from conduction within the MMZ. This effectively leads to a reduction in the size of linear equations from N to N - Npo with N being the total number of segments for the entire surface and Npo the number of segments over the POSR. The scheme would be appropriate for relatively large, flat surfaces, and at high frequencies. The ADI-FDTD scheme provides for both transient and steady state analyses. The restrictive Courant-Friedrich-Levy (CFL) condition on the time-step is removed, and so large time steps can be chosen even though the spatial grids are small. This report includes the problem definition, a detailed discussion of both the numerical techniques, and numerical implementations for simple surface geometries. Numerical solutions have been derived for a few simple situations.

Joshi, R. P.↗
NASA NTRS2021

A Three-Step Semi Analytical Algorithm (3SAA) for Estimating Inherent Optical Properties Over Oceanic, Coastal, and Inland Waters From Remote Sensing Reflectance

We present a three-step inverse model (3SAA) for estimating the inherent optical properties (IOPs) of surface waters from the remote sensing reflectance spectra, Rrs(). The derived IOPs include the total (a()), phytoplankton (aphy()), and colored detrital matter (acdm()), absorption coefficients, and the total (bb()) and particulate (bbp()) backscattering coefficients. The first step uses an improved neural network approach to estimate the diffuse attenuation coefficient of downwelling irradiance from Rrs. a() and bbp() are then estimated using the LS2 model (Loisel et al., 2018), which does not require spectral assumptions on IOPs and hence can assess a() and bb() at any wavelength at which Rrs() is measured. Then, an inverse optimization algorithm is combined with an optical water class (OWC) approach to assess aphy() and acdm() from anw().The proposed model is evaluated using an in situ dataset collected in open oceanic, coastal, and inland waters. Comparisons with other standard semi-analytical algorithms (QAA and GSM), as well as match-up exercises, have also been performed. The applicability of the algorithm on OLCI observations was assessed through the analysis of global IOPs spatial patterns derived from 3SAA and GSM. The good performance of 3SAA is manifested by median absolute percentage differences (MAPD) of 13%, 23%, 34% and 34% for bbp(443), anw(443), aphy(443) and acdm(443), respectively for oceanic waters. Due to the absence of spectral constraints on IOPs in the inversion of total IOPs, and the adoption of an OWC-based approach, the performance of 3SAA is only slightly degraded in bio-optical complex inland waters.

ocean color↗
DOE OSTI2016

A novel software and conceptual design of the hardware platform for intensity modulated radiation therapy

Purpose: To deliver high quality intensity modulated radiotherapy (IMRT) using a novel generalized sparse orthogonal collimators (SOCs), the authors introduce a novel direct aperture optimization (DAO) approach based on discrete rectangular representation. Methods: A total of seven patients—two glioblastoma multiforme, three head & neck (including one with three prescription doses), and two lung—were included. 20 noncoplanar beams were selected using a column generation and pricing optimization method. The SOC is a generalized conventional orthogonal collimators with N leaves in each collimator bank, where N = 1, 2, or 4. SOC degenerates to conventional jaws when N = 1. For SOC-based IMRT, rectangular aperture optimization (RAO) was performed to optimize the fluence maps using rectangular representation, producing fluence maps that can be directly converted into a set of deliverable rectangular apertures. In order to optimize the dose distribution and minimize the number of apertures used, the overall objective was formulated to incorporate an L2 penalty reflecting the difference between the prescription and the projected doses, and an L1 sparsity regularization term to encourage a low number of nonzero rectangular basis coefficients. The optimization problem was solved using the Chambolle–Pock algorithm, a first-order primal–dual algorithm. Performance of RAO was compared to conventional two-step IMRT optimization including fluence map optimization and direct stratification for multileaf collimator (MLC) segmentation (DMS) using the same number of segments. For the RAO plans, segment travel time for SOC delivery was evaluated for the N = 1, N = 2, and N = 4 SOC designs to characterize the improvement in delivery efficiency as a function of N. Results: Comparable PTV dose homogeneity and coverage were observed between the RAO and the DMS plans. The RAO plans were slightly superior to the DMS plans in sparing critical structures. On average, the maximum and mean critical organ doses were reduced by 1.94% and 1.44% of the prescription dose. The average number of delivery segments was 12.68 segments per beam for both the RAO and DMS plans. The N = 2 and N = 4 SOC designs were, on average, 1.56 and 1.80 times more efficient than the N = 1 SOC design to deliver. The mean aperture size produced by the RAO plans was 3.9 times larger than that of the DMS plans. Conclusions: The DAO and dose domain optimization approach enabled high quality IMRT plans using a low-complexity collimator setup. The dosimetric quality is comparable or slightly superior to conventional MLC-based IMRT plans using the same number of delivery segments. The SOC IMRT delivery efficiency can be significantly improved by increasing the leaf numbers, but the number is still significantly lower than the number of leaves in a typical MLC.

61 RADIATION PROTECTION AND DOSIMETRY↗
NASA NTRS2014

Improvement of Aerosol Optical Depth Retrieval over Hong Kong from a Geostationary Meteorological Satellite Using Critical Reflectance with Background Optical Depth Correction

Despite continuous efforts to retrieve aerosol optical depth (AOD) using a conventional 5-channelmeteorological imager in geostationary orbit, the accuracy in urban areas has been poorer than other areas primarily due to complex urban surface properties and mixed aerosol types from different emission sources. The two largest error sources in aerosol retrieval have been aerosol type selection and surface reflectance. In selecting the aerosol type from a single visible channel, the season-dependent aerosol optical properties were adopted from longterm measurements of Aerosol Robotic Network (AERONET) sun-photometers. With the aerosol optical properties obtained fromthe AERONET inversion data, look-up tableswere calculated by using a radiative transfer code: the Second Simulation of the Satellite Signal in the Solar Spectrum (6S). Surface reflectance was estimated using the clear sky composite method, awidely used technique for geostationary retrievals. Over East Asia, the AOD retrieved from the Meteorological Imager showed good agreement, although the values were affected by cloud contamination errors. However, the conventional retrieval of the AOD over Hong Kong was largely underestimated due to the lack of information on the aerosol type and surface properties. To detect spatial and temporal variation of aerosol type over the area, the critical reflectance method, a technique to retrieve single scattering albedo (SSA), was applied. Additionally, the background aerosol effect was corrected to improve the accuracy of the surface reflectance over Hong Kong. The AOD retrieved froma modified algorithmwas compared to the collocated data measured by AERONET in Hong Kong. The comparison showed that the new aerosol type selection using the critical reflectance and the corrected surface reflectance significantly improved the accuracy of AODs in Hong Kong areas,with a correlation coefficient increase from0.65 to 0.76 and a regression line change from τMI [basic algorithm] = 0.41τAERONET + 0.16 to τMI [new algorithm] = 0.70τAERONET + 0.01.

aerosols↗
DOE OSTI2016

Particle-image Velocimetry Sensitivity through Automatic Differentiation

A commonly encountered difficulty in modern experimental fluid dynamics experiments is the need for full-flow visualization and measurement; traditional techniques such as hot-wire probes and laser Doppler velocimetry cannot provide such information without difficult scanning procedures. Particle-image velocimetry (PIV) is a measurement technique which can overcome this spatial limitation to provide real-time three-dimensional fluid velocity vectors for an entire plane within the flow, as opposed to a single point that is provided by other measurement techniques. To accomplish this, a PIV system seeds the fluid flow with tracer particles; these must have nearly the same density as that of the fluid and a high drag coefficient to ensure that they follow the fluid flow without deviation. A laser is sent through optics to expand its beam into a sheet which illuminates the plane of interest in the fluid. The particles reflect more light than the fluid, and their image is recorded using a high speed camera. Due to limitations in bandwidth, video cannot be taken fast enough to be useful, and instead pairs of images are captured with a known separation in time. The illuminated particles will move between the two images taken by the camera, and their relative motion can be used to compute the fluid velocity at that point. When a single camera is used, only two-dimensional velocity vectors can be computed, but with the addition of another camera this can be improved to full three-dimensional vectors. The flow's time evolution can be measured by taking series of image pairs and processing each pair to get the velocity field at the time between each image pair. Only the two-dimensional version will be considered in this work. To compute the velocity vectors of the flow, the images are broken up into interrogation regions, usually in a uniform grid. In the most basic PIV algorithm, corresponding regions in the image pairs are cross-correlated to compute the most likely displacement from one image to the next. This is used since the cross correlation will achieve a maximum value when the two interrogation regions line up the best. The location of the peak in the correlation map and scaling information for the optical setup is used to compute the particle displacement from one image to the next. This process is similar to holding two semi-transparent images and moving one until they match as well as possible. Sub-pixel displacements between images can be obtained by fitting a Gaussian function to the correlation peak. Knowing the displacement vector for each set of interrogation regions in the image pair, it is a simple matter to divide by the time gap between the two images to compute the fluid velocity at the point in time and space. Since there are many interrogation regions in each image pair, the velocity field for a sheet is sampled according to the arrangement of the interrogation regions. The sensitivity of the final computed vector to the input images is difficult to compute, and has led most considerations of uncertainty in the PIV field to apply Monte Carlo methods and correlations to quantify the systematic and random errors which result from application of various PIV algorithms. In this work, automatic differentiation (AD) has been used to trace the dependence of velocity vectors computed using the correlation-peak PIV algorithm back to the input values given to an artificial image generator. These initial results form the foundation for continued work in the application of AD to algorithms in PIV. They demonstrate that AD methods can be successfully applied to PIV algorithms without undue complexity. The derivatives of each individual vector with regard to the coefficients of the velocity field used to create the artificial images can be computed and the statistics of them considered collectively, or correlated between each other to quantify the effectiveness of a PIV algorithm. With the initial implementation complete, adding extra algorithms is relatively simple and considerable data can be generated for analysis. It is hoped that application of these methods can lead to improved PIV algorithms which demonstrate better sensitivity to small variations in particle displacement. While these methods can be used as an analysis tool for refinement of new algorithms, computation of actual resultant vector sensitivities to input images' individual pixels is the next logical step of this work. Such information will be analyzed with experimental data to produce a better understanding of how output vectors' uncertainty relates to the actual image pairs, and lead to automatic bad vector removal. (authors)

42 ENGINEERING↗
DOE OSTI2019

Electron Cloud Studies at Fermilab

The presence of unwanted electrons in an accelerator vacuum chamber, known as E-cloud (E-cloud) can potentially cause operational problems in the Fermilab Main Injector (MI) and Recycler Ring (RR). E-cloud has caused instability in the RR in the past, and although it is not currently a problem, there is measurable E-cloud in the MI accelerator. There are planned beam intensity increases due to upgrades of the Fermilab accelerator complex, so E-cloud could become a problem. Some work has been done by others previously to understand how low SEY (Secondary Electron Yield coefficient) coatings might mitigate production of E-cloud, and to model the mechanism whereby E-cloud causes beam instability. Using previous studies as a base, this research took several approaches toward understanding the risk of E-cloud at Fermilab. The evolution of the SEY of the SS316L (stainless steel), TiN coated SS316L, and amorphous carbon coated SS316L were measured in-situ using a SEY meas urement station in the MI tunnel. The SEY of these materials change over time either due to bombardment of the E-cloud, or disruption of vacuum conditions. The SEY evolution was tracked over a several year period to find out how long it takes for the SEY of each material to reach its lowest level, and how much the SEY rises during deconditioning periods of poor vacuum. The SEY measurement results can be used to determine whether the SS316L will be a problem at upgrade intensities, and if so, whether or not TiN and A-Carbon coatings can mitigate E-cloud related problems sufficiently. Direct measurements of the E-cloud were done as well, and compared to simulation. The E-cloud bombardment rate was measured at different beam intensities and bunch lengths. It was possible to get detailed information on how the E-cloud varies over the acceleration cycle, where sensitivity to bunch length is reflected in the evolution of the E-cloud. The Retarding Field Analyzer (RFA) measuring th e E-cloud bombardment rate was near the instrument that is u! sed to measure the SEY of the beam pipe material. This proximity provided an accurate SEY value for simulations, so that the simulated E-cloud bombardment rate could be better compared to measurement. Bunch length measurements and computations generated accurate bunch length information also needed as input for simulations. After this careful control of the input parameters, the POSINST simulations of E-cloud were a good match to direct measurements. This gave confidence that predictions could be made concerning the E-cloud densities at upgrade intensities. These densities were compared against corresponding proton densities to predict the SEY required to avoid instabilities. That prediction and the information provided by the SEY measurements provide helpful information regarding the risks of E-cloud effects at future beam intensities at Fermilab.

43 PARTICLE ACCELERATORS↗
DOE OSTI2022

Methods for Determining the Infrared Complex Refractive Indices n (λ) and k (λ) from Organic Solid Powders: Comparison of a Sucrose Single Crystal and a Pressed Pellet

There exists a need to measure optical constants n and k for many materials, especially organic solids in the infrared spectral region, but at present, there are only two established methods for making such measurements: single-angle (SA) reflectance and infrared spectroscopic ellipsometry (IRSE). The former derives the n/k vectors via measurement of the change in the amplitude of light reflected from the sample and then uses the Kramers–Kronig transform, while the latter measures the change in polarization of the reflected light from the sample to derive n/k. Here we make a direct comparison between the n/k vectors derived from SA measurements of a sucrose single crystal and the n/k values obtained after the crystal had been crushed, ground to powder, and pressed into a pellet form. The pellet was measured not only by SA but also via the IRSE method, and we found that pressing specular-quality pellets was critical, especially for the SA data. The n/k values of this typical organic molecule all agree to within 0.04 root-mean-square error, confirming that a pressed pellet can be used instead of a crystal to derive the intrinsic n/k vectors. In conclusion, this result is important for the many organic materials that, unlike sugar, do not have centimeter-sized crystals readily available.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1994

The Balmer-dominated northeast limb of the Cygnus loop supernova remnant

We present a comprehensive investigation of the Balmar-dominated northeast limb of the Cygnus Loop supernova remnant. Data presented include H alpha (O III), and X-ray images, UV and visible spectrophotometry, and high-resolution spectroscophy. The two relatively bright Balmer-dominated filaments visible on the POSS prints are seen to be part of a very smooth and regular complex of filaments. These filaments mark the current location of the blast wave and are seen to bound the sharply limb-brightened X-ray emission, including the previously reported X-ray, 'halo.' The (O III)/h beta ratio throughout the region is approximately 0.1, except for regions in which the shock is undergoing a transition from nonradiative to incomplete radiative to incomplete radiative conditions. At these locations (O III) emission from the cooling region is quite strong, while collisionally excited Balmer-line emission can be weak because of photoionization of the preshock medium by UV from the nascent cooling region. As a result (O III)/H beta is greater than 100 in some locations. The nonradiative/radiative transition is best studied along the length of the northwestern of the two brightest filaments, where the shock velocity and swept-up column go from approximately 180 km/s and 10(exp 17)/sq cm at one end to approximately 140 km/s and 8 x 10(exp 17)/cm at the other. There are also a number of locations of such incomplete radiative emission where the shock has recently encountered denser regions with characteristic sizes of approximately 10(exp 18) cm. There is a considerable amount of evidence that the shock has decelerated from approximately 400 km/s to less than 200 km/s in the last 1000 yr. We interpret this as the result of the blast wave hitting the wall of a cavity which surround the supernova precursor and succeed in matching a wide range of data with a reflected shock model in which the density ofthe cavity wall is approximately 1.2/cu cm and the density in the interior of the cavity is about 0.09/cu cm. This analysis supports previous suggestions that the Cygnus Loop was a cavity explosion. We opredict that the X-ray shell will continue to brighten for the next 1000 years or so, during which time the shock front will become radiative throughout the northeast region. The similarity of the time scales for post-shock cooling, shock deceleration, and the 1ressure evolution of the region violates assumptions of steady flow, and greatly complicated interpretation of shock specptra in the region. This may be a common and important effect in many supernova remnants. High-resolutoin spectroscopy of the H alpha line profile shows that the broad component is surprisingly narrow (approximately 130 km/s). Consistency between this line width and the shock velocity inferred from spectrophotometry requires either rapid postshock equilibration of T(sub i) and T(sub e), or that the shock front has declerated from 180 km/s to 150 km/s in the last approximately 200 yr. While we cannot unambiguously distinquish between these models, the required deceleration is consistent with extrapolation of the shock decleration rate inferred over the last 1000 years. The width of the narrow component of the h alpha line profile is approximately 3 km/s, which indicates that the preshock medium has been heated to approximately 25,000 K by a shock precursor. The most promising types of precursor are heating either by fast neutrals which overtake the shock or by cosmic rays. In the latter case the observed properties of the precursor constrain the cosmic-ray diffusion coefficient to be between about 10(exp 23) sq cm/s. The characteristics of the precursor seem to be about right to explain the presence of weak (S II) and (N II) emission previously reported by Fesen & Itoh (1985).

Hester, J. Jeff↗