Search NASASearch

SEARCH · Search NASA

Results for “pair correlation function”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 145 records · Page 8

Statistical and linguistic features of DNA sequences

We present evidence supporting the idea that the DNA sequence in genes containing noncoding regions is correlated, and that the correlation is remarkably long range--indeed, base pairs thousands of base pairs distant are correlated. We do not find such a long-range correlation in the coding regions of the gene. We resolve the problem of the "non-stationary" feature of the sequence of base pairs by applying a new algorithm called Detrended Fluctuation Analysis (DFA). We address the claim of Voss that there is no difference in the statistical properties of coding and noncoding regions of DNA by systematically applying the DFA algorithm, as well as standard FFT analysis, to all eukaryotic DNA sequences (33 301 coding and 29 453 noncoding) in the entire GenBank database. We describe a simple model to account for the presence of long-range power-law correlations which is based upon a generalization of the classic Levy walk. Finally, we describe briefly some recent work showing that the noncoding sequences have certain statistical features in common with natural languages. Specifically, we adapt to DNA the Zipf approach to analyzing linguistic texts, and the Shannon approach to quantifying the "redundancy" of a linguistic text in terms of a measurable entropy function. We suggest that noncoding regions in plants and invertebrates may display a smaller entropy and larger redundancy than coding regions, further supporting the possibility that noncoding regions of DNA may carry biological information.

Non-NASA Center

Surface erythemal UV irradiance in the continental United States derived from ground-based and OMI observations: quality assessment, trend analysis and sampling issues

Surface full-sky erythemal dose rate (EDR) from the Ozone Monitoring Instrument (OMI) at both satellite overpass time and local noon time is evaluated against ground measurements at 31 sites from the US Department of Agriculture's (USDA) UV-B Monitoring and Research Program (UVMRP) over the period of 2005–2017. We find that both OMI overpass and solar noon time EDR are highly correlated with the measured counterparts (with a linear correlation coefficient of 0.90 and 0.88, respectively). Although the comparison statistics are improved with a longer time window (0.5–1.0 h) for pairing surface and OMI measurements, both OMI overpass and local noon time EDRs have 7 % overestimation that is larger than 6 % uncertainty in the ground measurements and show different levels of dependence on solar zenith angle (SZA) and to lesser extent on cloud optical depth. The ratio of EDR between local noon and OMI overpass time is often (95 % in frequency) larger than 1 with a mean of 1.18 in the OMI product; in contrast, the same ratio from surface observation is normally distributed with 22 % of the times less than 1 and a mean of 1.38. This contrast in part reflects the deficiency in the OMI surface UV algorithm that assumes constant atmospheric conditions between overpass and noon time. The probability density functions (PDFs) for both OMI and ground measurements of noontime EDR are in statistically significant agreement, showing dual peaks at ∼20 and ∼200 mW/sq. m, respectively; the latter is lower than 220 mW/ sq. m, the value at which the PDF of daily EDR from ground measurements peaks, and this difference indicates that the largest EDR value for a given day may not often occur at local noon. Lastly, statistically significant positive trends of EDR are found in the northeastern US in OMI data, but opposite trends are found within ground-based data (regardless of sampling for either noontime or daily averages). While positive trends are consistently found between OMI and surface data for EDR over the southern Great Plains (Texas and Oklahoma), their values are within the uncertainty of ground measurements. Overall, no scientifically sound trends can be found among OMI data for aerosol total and absorbing optical depth, cloud optical depth and total ozone to explain coherently the surface UV trends revealed either by OMI or ground-based estimates; these data also cannot reconcile trend differences between the two estimates (of EDR from OMI and surface observations). Future geostationary satellites with better spatiotemporal resolution data should help overcome spatiotemporal sampling issues inherent in OMI data products and therefore improve the estimates of surface UV flux and EDR from space.

Ozone Monitoring Instrument (OMI)

Convergent properties of vestibular-related brain stem neurons in the gerbil

Three classes of vestibular-related neurons were found in and near the prepositus and medial vestibular nuclei of alert or decerebrate gerbils, those responding to: horizontal translational motion, horizontal head rotation, or both. Their distribution ratios were 1:2:2, respectively. Many cells responsive to translational motion exhibited spatiotemporal characteristics with both response gain and phase varying as a function of the stimulus vector angle. Rotationally sensitive neurons were distributed as Type I, II, or III responses (sensitive to ipsilateral, contralateral, or both directions, respectively) in the ratios of 4:6:1. Four tested factors shaped the response dynamics of the sampled neurons: canal-otolith convergence, oculomotor-related activity, rotational Type (I or II), and the phase of the maximum response. Type I nonconvergent cells displayed increasing gains with increasing rotational stimulus frequency (0.1-2.0 Hz, 60 degrees /s), whereas Type II neurons with convergent inputs had response gains that markedly decreased with increasing translational stimulus frequency (0.25-2.0 Hz, +/-0.1 g). Type I convergent and Type II nonconvergent neurons exhibited essentially flat gains across the stimulus frequency range. Oculomotor-related activity was noted in 30% of the cells across all functional types, appearing as burst/pause discharge patterns related to the fast phase of nystagmus during head rotation. Oculomotor-related activity was correlated with enhanced dynamic range compared with the same category that had no oculomotor-related response. Finally, responses that were in-phase with head velocity during rotation exhibited greater gains with stimulus frequency increments than neurons with out-of-phase responses. In contrast, for translational motion, neurons out of phase with head acceleration exhibited low-pass characteristics, whereas in-phase neurons did not. Data from decerebrate preparations revealed that although similar response types could be detected, the sampled cells generally had lower background discharge rates, on average one-third lower response gains, and convergent properties that differed from those found in the alert animals. On the basis of the dynamic response of identified cell types, we propose a pair of models in which inhibitory input from vestibular-related neurons converges on oculomotor neurons with excitatory inputs from the vestibular nuclei. Simple signal convergence and combinations of different types of vestibular labyrinth information can enrich the dynamic characteristics of the rotational and translational vestibuloocular responses.

Non-NASA Center

Antigen S1, encoded by the MIC1 gene, is characterized as an epitope of human CD59, enabling measurement of mutagen-induced intragenic deletions in the AL cell system

S1 cell membrane antigen is encoded by the MIC1 gene on human chromosome 11. This antigen has been widely used as a marker for studies in gene mapping or in analysis of mutagen-induced gene deletions/mutations, which utilized the human-hamster hybrid cell-line, AL-J1, carrying human chromosome 11. Evidence is presented here which identifies S1 as an epitope of CD59, a cell membrane complement inhibiting protein. E7.1 monoclonal antibody, specific for the S1 determinant, was found to react strongly with membrane CD59 in Western blotting, and to bind to purified, urinary form of CD59 in ELISAs. Cell membrane expression of S1 on various cell lines always correlated with that of CD59 when examined by immunofluorescent staining. In addition, E7.1 antibody inhibited the complement regulatory function of CD59. Identification of S1 protein as CD59 has increased the scope of the AL cell system by enabling analysis of intragenic mutations, and multiplex PCR analysis of mutated cells is described, showing variable loss of CD59 exons.

Non-NASA Center

High-Altitude Emission from Pulsar Slot Gaps: The Crab Pulsar

We present results of a 3D model of optical to gamma-ray emission from the slot gap accelerator of a rotation-powered pulsar. Primary electrons accelerating to high-altitudes in the unscreened electric field of the slot gap reach radiation-reaction limited Lorentz factors of approx. 2 x 10(exp 7), while electron-positron pairs from lower-altitude cascades flow along field lines interior to the slot gap. The curvature, synchrotron and inverse Compton radiation of both primary electrons and pairs produce a broad spectrum of emission from infra-red to GeV energies. Both primaries and pairs undergo cyclotron resonant absorption of radio photons, allowing them to maintain significant pitch angles. Synchrotron radiation from pairs with a power-law energy spectrum from gamma = 10(exp 2) - 10(exp 5), dominate the spectrum up to approx. 10 MeV. Synchrotron and curvature radiation of primaries dominates from 10 MeV up to a few GeV. We examine the energy-dependent pulse profiles and phase-resolved spectra for parameters of the Crab pulsar as a function of magnetic inclination alpha and viewing angle zeta, comparing to broad-band data. In most cases, the pulse profiles are dominated by caustics on trailing field lines. We also explore the relation of the high-energy and the radio profiles, as well as the possibility of caustic formation in the radio cone emission. We find that the Crab pulsar profiles and spectrum can be reasonably well reproduced by a model with alpha = 45deg and zeta approx. 100deg or 80deg. This model predicts that the slot gap emission below 200 MeV will exhibit correlations in time and phase with the radio emission.

Harding, Alice K.

Modeling Dynamic Crush Behavior of Carbon Fiber Reinforced Polymer Composite Structures Using MAT213

Modeling crushing of carbon fiber reinforced polymer (CFRP) composites is challenging, and current simulation methodologies involve tuning of non-physical parameters. MAT213, a next-generation material model, has advanced functionality to better simulate dynamic impact loading. The objective of the present investigation is to evaluate the potential for using MAT213 to simulate dynamic crushing of CFRPs. Two sets of simulations were performed: one for calibration based on a set of coupon-level experiments and another for prediction of the response of structural elements. Simulations involving dynamic crushing of flat specimens were iteratively run to calibrate model parameters. The calibration demonstrated that MAT213 could produce a simulated force-displacement response within experimental scatter. The simulated failure morphology was also comparable to the experiments. After successful calibration, predictive simulations of dynamic crushing of C-channel shaped specimens were completed using a simulated crash sled test rig and two pairs of impactor mass/velocity conditions. The simulated force-displacement curve in the crash sled simulations for the lower-velocity condition fell within the experimental scatter, but the stable crush force was underpredicted by 27% in the higher-velocity simulations. Better correlation in the lower-velocity test condition likely results from the calibration condition being a similar velocity to the lower-velocity crash sled condition.

Dynamic Crush

The Effect of the Wall Contact and Post-Growth, Cool-Down on Defects in CdTe Crystals Grown By 'Contactless' PVT

To take a maximum advantage of materials processing in microgravity for understanding the effects of gravity, gravity-independent effects should be minimized. In crystal growth, the quality of the grown crystals may depend, among other factors, on their interaction with the walls of the processing container during and after growth, and on the rate of the crystal cool-down at the end of the process. To investigate the above phenomena, a series of CdTe crystal growth processes was carried out. The crystals were grown by physical vapor transport without contact with the side walls of the silica glass ampoules. To eliminate the effect of the seed quality, and to reduce the number of nuclei and related crystal grains, the Low Supersaturation Nucleation technique was applied. The source temperature was 930 C, the undercooling was a few degrees. The crystals, having the diameter of 25 mm, grew at the rate of a few mm per day. The post-growth cool-down to the room temperature was conducted at different rates, and lasted from a few minutes to four days. The crystals were characterized using chemical etching, low temperature luminescence, and Synchrotron White Beam X-ray Topography techniques. The dislocation (etch pit) density was measured and its distribution was analyzed by comparison with Poisson curves and with the Normalized Radial Distribution Correlation Function. In the regions where the crystal is in contact with silica, the materials show a considerable strain field which extends for a few mm or more from the silica-crystal interface. In the reference crystal grown with contact with the ampoule walls, and when the crystals are cooled at the highest rates, the etch pit/dislocation density is in the high 10(exp 5) per square centimeter region. Typical EPD values for lower cool-down rates are in the lower 10(exp 4) per square centimeter region. In some areas the actual dislocation density was about 10(exp 3) per square centimeter or even less. No apparent effect of the cool-down rate on polygonization was observed. Low temperature PL spectra show, that the dominant peak is (D(sup 0), h) and (A(sup 0), e) for samples with low and high dislocation densities, respectively. For low EPD crystals a peak at 1.45 eV with 21 meV phonon replicas was observed and attributed to donor-acceptor pair to neutral copper-acceptor transition. In high EPD crystals this PL structure was not observed, apparently due to the masking effect of the strong contribution from the dislocation band.

Palosz, W.

Using Digital Image Correlation for Material Testing at Cryogenic Temperatures

Recent efforts were made to parametrize an aluminum alloy at cryogenic temperatures to calibrate finite element models simulating material deformation. Digital image correlation (DIC) was used to measure material properties using front/back strains of specimens during tensile tests. Data acquired from DIC was compared to physical extensometer measurements. Two camera pairs (front and back) were mounted via aluminum extrusion to the exterior of an environmental chamber of a hydraulic test stand. Images were taken through single-pane optical glass windows of the environment chamber. Positive pressure bags encased the cameras, attaching to the perimeter of the glass to keep out moisture and prevent windows from frosting. Figure 1 shows the environmental setup and the camera setup during cryogenic testing. Some of the challenges during testing at cryogenic temperatures include the gaseous flow of liquid nitrogen in the optical path of the mounted cameras, frost buildup on the windows as seen in Figure 1b, and frost on the test coupon. The temperature gradient present during the introduction of liquid nitrogen into the chamber can increase the projection error of images captured of the specimen. Figure 2 shows the estimated projection error from DIC as a function of temperature during the cool down process. This presentation addresses some of the concerns and observations when using DIC for testing at cryogenic temperatures using methods of isolated camera system and techniques for maintaining the glass panel clear from frost.

Joseph Cochran

Proteomic Assessment of Fluid Shifts and Association with Visual Impairment and Intracranial Pressure in Twin Astronauts

BACKGROUND: Astronauts participating in long duration space missions are at an increased risk of physiological disruptions. The development of visual impairment and intracranial pressure (VIIP) syndrome is one of the leading health concerns for crew members on long-duration space missions; microgravity-induced fluid shifts and chronic elevated cabin CO2 may be contributing factors. By studying physiological and molecular changes in one identical twin during his 1-year ISS mission and his ground-based co-twin, this work extends a current NASA-funded investigation to assess space flight induced "Fluid Shifts" in association with the development of VIIP. This twin study uniquely integrates physiological and -omic signatures to further our understanding of the molecular mechanisms underlying space flight-induced VIIP. We are: (i) conducting longitudinal proteomic assessments of plasma to identify fluid regulation-related molecular pathways altered by long-term space flight; and (ii) integrating physiological and proteomic data with genomic data to understand the genomic mechanism by which these proteomic signatures are regulated. PURPOSE: We are exploring proteomic signatures and genomic mechanisms underlying space flight-induced VIIP symptoms with the future goal of developing early biomarkers to detect and monitor the progression of VIIP. This study is first to employ a male monozygous twin pair to systematically determine the impact of fluid distribution in microgravity, integrating a comprehensive set of structural and functional measures with proteomic, metabolomic and genomic data. This project has a broader impact on Earth-based clinical areas, such as traumatic brain injury-induced elevations of intracranial pressure, hydrocephalus, and glaucoma. HYPOTHESIS: We predict that the space-flown twin will experience a space flight-induced alteration in proteins and peptides related to fluid balance, fluid control and brain injury as compared to his pre-flight protein/peptide signatures. Conversely, the trajectory of these protein signatures will remain relatively constant in his ground based co-twin. METHODS: We are using proteomic and standard immunoelectrophoresis techniques to delineate the change in protein signatures throughout the course of a long duration space flight in relation to the development of VIIP. We are also applying a novel cell-based metaboloic organ system assay ("Organs on a Plate") to address how these circulating biomarkers affect physiological processes at the cellular and organ level which could result in VIIP symptoms. These molecular data will be correlated with physiological measures (eg. extra and intracellular fluid volume, vascular filling/flow patterns, MRI, and Optic Coherence Tomography. DISCUSSION: Pre- and in-flight data collection is in progress for the space-flown twin, and similar data have been obtained from the ground-based twin. Biosamples will be batch processed when received from ISS after the conclusion of the 1-year mission. Omic and Physiological measures from the twin astronauts will be compared to similar data being collected on twin subjects who participated in simulated microgravity study. bed rest study.

Rana, Brinda K.

A method of classification for multisource data in remote sensing based on interval-valued probabilities

An axiomatic approach to intervalued (IV) probabilities is presented, where the IV probability is defined by a pair of set-theoretic functions which satisfy some pre-specified axioms. On the basis of this approach representation of statistical evidence and combination of multiple bodies of evidence are emphasized. Although IV probabilities provide an innovative means for the representation and combination of evidential information, they make the decision process rather complicated. It entails more intelligent strategies for making decisions. The development of decision rules over IV probabilities is discussed from the viewpoint of statistical pattern recognition. The proposed method, so called evidential reasoning method, is applied to the ground-cover classification of a multisource data set consisting of Multispectral Scanner (MSS) data, Synthetic Aperture Radar (SAR) data, and digital terrain data such as elevation, slope, and aspect. By treating the data sources separately, the method is able to capture both parametric and nonparametric information and to combine them. Then the method is applied to two separate cases of classifying multiband data obtained by a single sensor. In each case a set of multiple sources is obtained by dividing the dimensionally huge data into smaller and more manageable pieces based on the global statistical correlation information. By a divide-and-combine process, the method is able to utilize more features than the conventional maximum likelihood method.

Kim, Hakil

Geomagnetic Secular Variation Prediction with Thermal Heterogeneous Boundary Conditions

It has long been conjectured that thermal heterogeneity at the core-mantle boundary (CMB) affects the geodynamo substantially. The observed two pairs of steady and strong magnetic flux lobes near the Polar Regions and the low secular variation in the Pacific over the past 400 years (and perhaps longer) are likely the consequences of this CMB thermal heterogeneity. There are several studies on the impact of the thermal heterogeneity with numerical geodynamo simulations. However, direct correlation between the numerical results and the observations is found very difficult, except qualitative comparisons of certain features in the radial component of the magnetic field at the CMB. This makes it difficult to assess accurately the impact of thermal heterogeneity on the geodynamo and the geomagnetic secular variation. We revisit this problem with our MoSST_DAS system in which geomagnetic data are assimilated with our geodynamo model to predict geomagnetic secular variations. In this study, we implement a heterogeneous heat flux across the CMB that is chosen based on the seismic tomography of the lowermost mantle. The amplitude of the heat flux (relative to the mean heat flux across the CMB) varies in the simulation. With these assimilation studies, we will examine the influences of the heterogeneity on the forecast accuracies, e.g. the accuracies as functions of the heterogeneity amplitude. With these, we could be able to assess the model errors to the true core state, and thus the thermal heterogeneity in geodynamo modeling.

Kuang, Weijia

Decorrelation Times of Photospheric Fields and Flows

We use autocorrelation to investigate evolution in flow fields inferred by applying Fourier Local Correlation Tracking (FLCT) to a sequence of high-resolution (0.3 "), high-cadence (approx = 2 min) line-of-sight magnetograms of NOAA active region (AR) 10930 recorded by the Narrowband Filter Imager (NFI) of the Solar Optical Telescope (SOT) aboard the Hinode satellite over 12 - 13 December 2006. To baseline the timescales of flow evolution, we also autocorrelated the magnetograms, at several spatial binnings, to characterize the lifetimes of active region magnetic structures versus spatial scale. Autocorrelation of flow maps can be used to optimize tracking parameters, to understand tracking algorithms f susceptibility to noise, and to estimate flow lifetimes. Tracking parameters varied include: time interval Delta t between magnetogram pairs tracked, spatial binning applied to the magnetograms, and windowing parameter sigma used in FLCT. Flow structures vary over a range of spatial and temporal scales (including unresolved scales), so tracked flows represent a local average of the flow over a particular range of space and time. We define flow lifetime to be the flow decorrelation time, tau . For Delta t > tau, tracking results represent the average velocity over one or more flow lifetimes. We analyze lifetimes of flow components, divergences, and curls as functions of magnetic field strength and spatial scale. We find a significant trend of increasing lifetimes of flow components, divergences, and curls with field strength, consistent with Lorentz forces partially governing flows in the active photosphere, as well as strong trends of increasing flow lifetime and decreasing magnitudes with increases in both spatial scale and Delta t.

Welsch, B. T.

Prediction of Aerodynamic Coefficient using Genetic Algorithm Optimized Neural Network for Sparse Data

Wind tunnels use scale models to characterize aerodynamic coefficients, Wind tunnel testing can be slow and costly due to high personnel overhead and intensive power utilization. Although manual curve fitting can be done, it is highly efficient to use a neural network to define the complex relationship between variables. Numerical simulation of complex vehicles on the wide range of conditions required for flight simulation requires static and dynamic data. Static data at low Mach numbers and angles of attack may be obtained with simpler Euler codes. Static data of stalled vehicles where zones of flow separation are usually present at higher angles of attack require Navier-Stokes simulations which are costly due to the large processing time required to attain convergence. Preliminary dynamic data may be obtained with simpler methods based on correlations and vortex methods; however, accurate prediction of the dynamic coefficients requires complex and costly numerical simulations. A reliable and fast method of predicting complex aerodynamic coefficients for flight simulation I'S presented using a neural network. The training data for the neural network are derived from numerical simulations and wind-tunnel experiments. The aerodynamic coefficients are modeled as functions of the flow characteristics and the control surfaces of the vehicle. The basic coefficients of lift, drag and pitching moment are expressed as functions of angles of attack and Mach number. The modeled and training aerodynamic coefficients show good agreement. This method shows excellent potential for rapid development of aerodynamic models for flight simulation. Genetic Algorithms (GA) are used to optimize a previously built Artificial Neural Network (ANN) that reliably predicts aerodynamic coefficients. Results indicate that the GA provided an efficient method of optimizing the ANN model to predict aerodynamic coefficients. The reliability of the ANN using the GA includes prediction of aerodynamic coefficients to an accuracy of 110% . In our problem, we would like to get an optimized neural network architecture and minimum data set. This has been accomplished within 500 training cycles of a neural network. After removing training pairs (outliers), the GA has produced much better results. The neural network constructed is a feed forward neural network with a back propagation learning mechanism. The main goal has been to free the network design process from constraints of human biases, and to discover better forms of neural network architectures. The automation of the network architecture search by genetic algorithms seems to have been the best way to achieve this goal.

Rajkumar, T.

Crystal Growth of ZnSe and Related Ternary Compound Semiconductors by Vapor Transport

The objective of the project is to determine the relative contributions of gravity-driven fluid flows to the compositional distribution, incorporation of impurities and defects, and deviation from stoichiometry observed in the crystals grown by vapor transport as results of buoyance-driven convection and growth interface fluctuations caused by irregular fluid-flows. ZnSe and related ternary compounds, such as ZnSeS and ZnSeTe, were grown by vapor transport technique with real time in-situ non-invasive monitoring techniques. The grown crystals were characterized extensively to correlate the grown crystal properties with the growth conditions. The following are the research progress in the past two years. In-situ monitoring of partial pressure by optical absorption technique and visual observation of the growing crystal were performed during vapor growth of ZnSe. Low-temperature photoluminescence (PL) spectra and glow discharge mass spectroscopy (GDMS) were measured on ZnSe starting materials provided by various vendors and on bulk crystals grown from these starting materials by physical vapor transport (PVT) to study the effects of purification and contamination during crystal growth process. Optical characterization was performed on wafers sliced from the grown crystals of ZnSe, ZnTe and ZnSe(1-x),Te(x), (0<x<0.4). Energy band gaps at room temperature were determined from optical transmission measurements and a best fit curve to the band gap vs. composition, x, data gives a bowing parameter of 1.45. Low-temperature photoluminescence (PL) spectra of ZnSe and ZnTe were dominated by near band edge emissions and no deep donor-acceptor pairs were observed. The PL spectrum exhibited a broad emission for the ZnSe(1-x),Te(x), samples, 0.09<x<0.39. The single broad PL emission spectra and the spectra measured as a function of temperature were interpreted as being associated with the exciton bound to Te clusters because of the high Te content in these samples. To validate numerical codes, in-situ monitoring during the PVT of HgI2 was designed. Optical absorption spectra of the vapor phase over HgI2 were measured for wavelengths between 200 and 600nm at sample temperatures between 349 and 610K. The Beer's Law constants for 15 wavelengths between 200 and 440 nm were determined. From these constants the vapor pressure of HgI2 was established as a function of temperature for the liquid and the solid Beta-phases To characterize the growth conditions during the PVT growth of In-doped ZnSe the optical absorbance of the vapor phase over the In-Se system were measured and were used to obtain the partial pressures of Se2(g) and In2Se(g).

Su, Ching-Hua

The structure and energetics of Cr(CO)6 and Cr(CO)5

The geometric structure of Cr(CO)6 is optimized at the modified coupled pair functional (MCPF), single and double excitation coupled-cluster (CCSD) and CCSD(T) levels of theory (including a perturbational estimate for connected triple excitations), and the force constants for the totally symmetric representation are determined. The geometry of Cr(CO)5 is partially optimized at the MCPF, CCSD, and CCSD(T) levels of theory. Comparison with experimental data shows that the CCSD(T) method gives the best results for the structures and force constants, and that remaining errors are probably due to deficiencies in the one-particle basis sets used for CO. The total binding energies of Cr(CO)6 and Cr(CO)5 are also determined at the MCPF, CCSD, and CCSD(T) levels of theory. The CCSD(T) method gives a much larger total binding energy than either the MCPF or CCSD methods. An analysis of the basis set superposition error (BSSE) at the MCPF level of treatment points out limitations in the one-particle basis used. Calculations using larger basis sets reduce the BSSE, but the total binding energy of Cr(CO)6 is still significantly smaller than the experimental value, although the first CO bond dissociation energy of Cr(CO)6 is well described. An investigation of 3s3p correlation reveals only a small effect. In the largest basis set, the total CO binding energy of Cr(CO)6 is estimated to be 140 kcal/mol at the CCSD(T) level of theory, or about 86 percent of the experimental value. The remaining discrepancy between the experimental and theoretical value is probably due to limitations in the one-particle basis, rather than limitations in the correlation treatment. In particular an additional d function and an f function on each C and O are needed to obtain quantitative results. This is underscored by the fact that even using a very large primitive set (1042 primitive functions contracted to 300 basis functions), the superposition error for the total binding energy of Cr(CO)6 is 22 kcal/mol at the MCPF level of treatment.

Barnes, Leslie A.

Positive selection moments identify potential functional residues in human olfactory receptors

Correlated mutation analysis and molecular models of olfactory receptors have provided evidence that residues in the transmembrane domains form a binding pocket for odor ligands. As an independent test of these results, we have calculated positive selection moments for the alpha-helical sixth transmembrane domain (TM6) of human olfactory receptors. The moments can be used to identify residues that have been preferentially affected by positive selection and are thus likely to interact with odor ligands. The results suggest that residue 622, which is commonly a serine or threonine, could form critical H-bonds. In some receptors a dual-serine subsite, formed by residues 622 and 625, could bind hydroxyl determinants on odor ligands. The potential importance of these residues is further supported by site-directed mutagenesis in the beta-adrenergic receptor. The findings should be of practical value for future physiological studies, binding assays, and site-directed mutagenesis.

Non-NASA Center

Photocatalytic Active Radiation Measurements and Use

Photocatalytic materials are being used to purify air, to kill microbes, and to keep surfaces clean. A wide variety of materials are being developed, many of which have different abilities to absorb various wavelengths of light. Material variability, combined with both spectral illumination intensity and spectral distribution variability, will produce a wide range of performance results. The proposed technology estimates photocatalytic active radiation (PcAR), a unit of radiation that normalizes the amount of light based on its spectral distribution and on the ability of the material to absorb that radiation. Photocatalytic reactions depend upon the number of electron-hole pairs generated at the photocatalytic surface. The number of electron-hole pairs produced depends on the number of photons per unit area per second striking the surface that can be absorbed and whose energy exceeds the bandgap of the photocatalytic material. A convenient parameter to describe the number of useful photons is the number of moles of photons striking the surface per unit area per second. The unit of micro-einsteins (or micromoles) of photons per m2 per sec is commonly used for photochemical and photoelectric-like phenomena. This type of parameter is used in photochemistry, such as in the conversion of light energy for photosynthesis. Photosynthetic response correlates with the number of photons rather than by energy because, in this photochemical process, each molecule is activated by the absorption of one photon. In photosynthesis, the number of photons absorbed in the 400 700 nm spectral range is estimated and is referred to as photosynthetic active radiation (PAR). PAR is defined in terms of the photosynthetic photon flux density measured in micro-einsteins of photons per m2 per sec. PcAR is an equivalent, similarly modeled parameter that has been defined for the photocatalytic processes. Two methods to measure the PcAR level are being proposed. In the first method, a calibrated spectrometer with a cosine receptor is used to measure the spectral irradiance. This measurement, in conjunction with the photocatalytic response as a function of wavelength, is used to estimate the PcAR. The photocatalytic response function is determined by measuring photocatalytic reactivity as a function of wavelength. In the second method, simple shaped photocatalytic response functions can be simulated with a broad-band detector with a cosine receptor appropriately filtered to represent the spectral response of the photocatalytic material. This second method can be less expensive than using a calibrated spectrometer.

Davis, Bruce A.

Direct Numerical Simulation of Transitional Multicomponent-Species Gaseous and Multicomponent-Liquid Drop-Laden Mixing

A model of multicomponent-liquid (MC-liquid) drop evaporation in a three-dimensional mixing layer is here exercised at larger Reynolds numbers than in a previous study, and transitional states are obtained. The gas phase is followed in an Eulerian frame and the multitude of drops is described in a Lagrangian frame. Complete coupling between phases is included with source terms in the gas conservation equations accounting for the drop/flow interaction in terms of drop drag, drop heating and species evaporation. The liquid composition, initially specified as a single-Gamma (SG) probability distribution function (PDF) depending on the molar mass is allowed to evolve into a linear combination of two SGPDFs, called the double-Gamma PDF (DGPDF). The compositions of liquid and vapor emanating from the drops are calculated through four moments of the DGPDFs, which are drop-specific and location-specific, respectively. The mixing layer is initially excited to promote the double pairing of its four initial spanwise vortices into an ultimate vortex in which small scales proliferate. Simulations are performed for four liquids of different compositions and the effect of the initial mass loading and initial free-stream gas temperature are explored. For reference, Simulations are also performed for gaseous multicomponent mixing layers for which the effect of Reynolds number is investigated. The results encompass examination of the global layer characteristics, flow visualizations and homogeneous-plane statistics at transition. Comparisons are performed with previous pre-transitional MC-liquid simulations and with transitional single-component (SC) liquid studies. It is found that MCC flows at transition, the classical energy cascade is of similar strength, but that the smallest scales contain orders of magnitude less energy than SC flows, which is confirmed by the larger viscous dissipation in the former case. Contrasting to pre-transitional MC flows, the vorticity and drop organization depend on the initial gas temperature, this being due to the drop/turbulence coupling. The vapor-composition mean molar mass and standard deviation distributions strongly correlate with the initial liquid-composition PDF; such a correlation only exists for the magnitude of the mean but not for that of the standard deviation. Unlike in pre-transitional situations, regions of large composition standard deviation no longer necessarily coincide with regions of large mean molar mass. The kinetic energy, rotational and composition characteristics, and dissipation are liquid specific and the variation among liquids is amplified with increasing free-stream gas temperature. Eulerian and Lagrangian statistics of gas-phase quantities show that the different. Observation framework may affect the perception of the flow characteristics. The gas composition, of which the first four moments are calculated, is shown to be close to, but distinct from a SGPDF. The PDF of the scalar dissipation rate is calculated for drop-laden layers and is shown to depart more significantly from the typically assumed Gaussian in gaseous flows than experimentally measured gaseous scalar dissipation rates, this being attributed to the increased heterogeneity due to drop/flow interactions.

tubulent multicompenent-species