Search NASA⌕ Search

SEARCH · Search NASA

Results for “Gaussian process”

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 1,711 records · Page 95

Barrier inhomogeneities and electronic transport of Pt contacts to relatively highly doped n-type 4H-SiC

The barrier characteristics of Pt contacts to relatively highly doped (∼1 × 10{sup 18 }cm{sup −3}) 4H-SiC were investigated using current-voltage (I-V) and capacitance-voltage (C-V) measurements in the temperature range of 160–573 K. The barrier height and ideally factor estimated from the I-V characteristics based on the thermionic emission model are abnormally temperature-dependent, which can be explained by assuming the presence of a double Gaussian distribution (GD) of inhomogeneous barrier heights. However, in the low temperature region (160–323 K), the obtained mean barrier height according to GD is lower than the actual mean value from C-V measurement. The values of barrier height determined from the thermionic field emission model are well consistent with those from the C-V measurements, which suggest that the current transport process could be modified by electron tunneling at low temperatures.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

SU-E-T-253: Open-Source Automatic Software for Quantifying Biological Assays of Radiation Effects

Purpose: Clonogenic cell survival is a common assay for quantifying the effect of drugs and radiation. Manual counting of surviving colonies can take 30–90seconds per plate, a major limitation for large studies. Currently available automatic counting tools are not easily modified for radiation oncology research. Our goal is to provide an open-source toolkit for precise, accurate and fast analysis of biological assays in radiation oncology. Methods: As an example analysis, we used HeLa cells incubated with gadolinium nanoparticles prior to irradiation. After treatment, the cells are grown for 14days to allow for colony formation. To analyze the colony growth, we capture images of each dish for archiving and automatic computer-based analysis. A FujifilmX20 camera is placed at the top of a box setup, 20cm above the sample, which is backlit by a LED lamp placed at the bottom of the box. We use a Gaussian filter (width=1.3mm) and color threshold (19–255). The minimum size for a colony to be counted is 1mm. For this example, 20 dishes with a large range of colonies were analyzed. Each dish was counted 3 times manually by 3 different users and then compared to our counter. Results: Automatic counting of cell colonies takes an average of 7seconds, enabling the analysis process to be accelerated 4–12 times. The average precision of the automatic counter was 1.7%. The Student t-test demonstrated the non-significant differences between the two counting methods (p=0.64). The ICC demonstrated the reliability of each method with ICC>0.999 (automatic) and ICC=0.95 (manual). Conclusion: We developed an open-source automatic toolkit for the analysis of biological assays in radiation oncology and demonstrated the accuracy, precision and effort savings for clonogenic cell survival quantification. This toolkit is currently being used in two laboratories for routine experimental analysis and will be made freely available on our departmental website.

60 APPLIED LIFE SCIENCES↗

Lorentz drift compensation in high harmonic generation in the soft and hard X-ray regions of the spectrum

Here, we present a semi-classical study of the effects of the Lorentz force on electrons during high harmonic generation in the soft and hard X-ray regions driven by near- and mid-infrared lasers with wavelengths from 0.8 to 20 μm, and at intensities below 10 15 W/cm 2 . The transverse extent of the longitudinal Lorentz drift is compared for both Gaussian focus and waveguide geometries. Both geometries exhibit a longitudinal electric field component that cancels the magnetic Lorentz drift in some regions of the focus, once each full optical cycle. We show that the Lorentz force contributes a super-Gaussian scaling which acts in addition to the dominant high harmonic flux scaling of λ -(5-6) due to quantum diffusion. We predict that the high harmonic yield will be reduced for driving wavelengths > 6 μm, and that the presence of dynamic spatial mode asymmetries results in the generation of both even and odd harmonic orders. Remarkably, we show that under realistic conditions, the recollision process can be controlled and does not shut off completely even for wavelengths >10 μm and recollision energies greater than 15 keV.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Structural Characterization of Crystalline Ice Nanoclusters

Water ice nanoclusters are useful analogs for studying a variety of processes that occur within icy grains in the extraterrestrial environment. The surface of ice nanoclusters prepared in the laboratory is similar to the surface of interstellar ice grains. In cold molecular clouds, the silicate cores of interstellar grains are typically approx. 100 nm in diameter and have a coating of impure amorphous water ice. Depositional, thermal and radiolytic processes leave the surface and subsurface molecules in a disordered state. In this state, structural defects become mobile and reactions of trapped gases and small molecules can occur. The large surface area of nanocluster deposits relative to their bulk allows for routine observation of such surface-mediated processes. Furthermore, the disordered surface and subsurface layers in nanocluster deposits mimic the structure of amorphous ice rinds found on interstellar dust grains. Transmission Electron Microscopy (TEM has been used tn characterize the crystallinity, growth mechanism, and size distribution of nanoclusters formed from a mixture of water vapor with an inert carrier gas that has been rapidly cooled to 77K. E M imaging reveals a Gaussian size distribution around a modal diameter that increases from approx. 15 to 30 nm as the percentage of water vapor within the mixture increases from 0.5 to 2.007, respectively . TEM bright and dark field imaging also reveals the crystalline nature of the clusters. h4any of the clusters show a mosaic structure in which crystalline domains originate at the center Other images show mirror planes that are separated by approx. 10 nm. Electron diffraction patterns of these clusters show that the clusters are composed of cubic ice with only a small hexagonal component. Further, the crystalline domain size is approximately the same as the modal diameter suggesting that the clusters are single crystals.

Blake, David↗

Mechanism of the G-protein mimetic nanobody binding to a muscarinic G-protein-coupled receptor

Protein–protein binding is key in cellular signaling processes. Molecular dynamics (MD) simulations of protein–protein binding, however, are challenging due to limited timescales. In particular, binding of the medically important G-protein-coupled receptors (GPCRs) with intracellular signaling proteins has not been simulated with MD to date. Here, we report a successful simulation of the binding of a G-protein mimetic nanobody to the M 2 muscarinic GPCR using the robust Gaussian accelerated MD (GaMD) method. Through long-timescale GaMD simulations over 4,500 ns, the nanobody was observed to bind the receptor intracellular G-protein-coupling site, with a minimum rmsd of 2.48 Å in the nanobody core domain compared with the X-ray structure. Binding of the nanobody allosterically closed the orthosteric ligand-binding pocket, being consistent with the recent experimental finding. In the absence of nanobody binding, the receptor orthosteric pocket sampled open and fully open conformations. The GaMD simulations revealed two low-energy intermediate states during nanobody binding to the M 2 receptor. The flexible receptor intracellular loops contribute remarkable electrostatic, polar, and hydrophobic residue interactions in recognition and binding of the nanobody. As a result, these simulations provided important insights into the mechanism of GPCR–nanobody binding and demonstrated the applicability of GaMD in modeling dynamic protein–protein interactions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Discovery of a quasar with double-peaked broad balmer emission lines

ABSTRACT Most massive galaxies contain a supermassive black hole (SMBH) at their centre. When galaxies merge, their SMBHs sink to the centre of the new galaxy, where they are thought to eventually merge. During this process, an SMBH binary is formed. The presence of two sets of broad emission lines in the optical spectrum of an active galactic nucleus (AGN) has been interpreted as evidence for two broad-line regions (BLRs), one surrounding each SMBH in a binary. We modelled the broad Balmer emission lines in the SDSS spectra of 373 extreme variability AGNs using one broad and several narrow Gaussian components. We report on the discovery of SDSS J021647.53 − 011341.5 (hereafter J0216) as a double-peaked broad emission line source. Among the 373 AGNs, there were five sources that are known as double-peaked emission line sources. Three of these have been reported as candidate SMBH binaries in previous studies. We present all six objects and their double-peaked broad Balmer emission lines, and discuss the implications for a tidal disruption event (TDE) interpretation of the extreme variability, assuming the double-peaked sources are SMBH binaries.

Terwel, Jacco H. (ORCID:0000000198343439)↗

A discrepancy within primate spatial vision and its bearing on the definition of edge detection processes in machine vision

The visual perception of form information is considered to be based on the functioning of simple and complex neurons in the primate striate cortex. However, a review of the physiological data on these brain cells cannot be harmonized with either the perceptual spatial frequency performance of primates or the performance which is necessary for form perception in humans. This discrepancy together with recent interest in cortical-like and perceptual-like processing in image coding and machine vision prompted a series of image processing experiments intended to provide some definition of the selection of image operators. The experiments were aimed at determining operators which could be used to detect edges in a computational manner consistent with the visual perception of structure in images. Fundamental issues were the selection of size (peak spatial frequency) and circular versus oriented operators (or some combination). In a previous study, circular difference-of-Gaussian (DOG) operators, with peak spatial frequency responses at about 11 and 33 cyc/deg were found to capture the primary structural information in images. Here larger scale circular DOG operators were explored and led to severe loss of image structure and introduced spatial dislocations (due to blur) in structure which is not consistent with visual perception. Orientation sensitive operators (akin to one class of simple cortical neurons) introduced ambiguities of edge extent regardless of the scale of the operator. For machine vision schemes which are functionally similar to natural vision form perception, two circularly symmetric very high spatial frequency channels appear to be necessary and sufficient for a wide range of natural images. Such a machine vision scheme is most similar to the physiological performance of the primate lateral geniculate nucleus rather than the striate cortex.

Jobson, Daniel J.↗

Preliminary Results of a Microgravity Investigation to Measure Net Charge on Granular Materials

Accurate characterization of the electrostatic charge on granular materials has typically been limited to materials with diameters on the order of 10 microns and below due to high settling velocities of larger particles. High settling velocities limit both the time and the acceptable uncertainty with which a measurement can be made. A prototype device has been developed at NASA Glenn Research Center (GRC) to measure coulombic charge on individual particles of granular materials that are 50 to 500 microns in diameter. This device, a novel extension of Millikan's classic oil drop experiment, utilizes the NASA GRC 2.2 second drop tower to extend the range of electrostatic charge measurements to accommodate moderate size granular materials. A dielectric material with a nominal grain diameter between 1.06 and 250 microns was tribocharged using a dry gas jet, suspended in a 5x10x10 cm enclosure during a 2.2 second period of microgravity and exposed to a known electric field. The response was recorded on video and post processed to allow tracking of individual particles. By determining the particle trajectory and velocity, estimates of the coulombic charge were made. Over 30 drops were performed using this technique and the analysis showed that first order approximations of coulombic charge could successfully be obtained, with the mean charge of 3.4E-14 coulombs measured for F-75 Ottawa quartz sand. Additionally, the measured charge showed a near-Gaussian distribution, with a standard deviation of 2.14E -14 coulombs.

Green, Robert D.↗

Bayesian optimization of laser wakefield acceleration via spectral pulse shaping

In this paper, we investigate the effect of spectral pulse shaping of the laser driver on the performance of channel-guided, laser–plasma accelerators. The study was carried out with the assistance of Bayesian optimization using particle-in-cell simulations. We used a realistic plasma profile based on a novel optical-field-ionized channel technique with ionization injection and low, on-axis plasma densities to maximize the energy gain of the electron bunch trailing the laser. Spectral shaping allows us to modify the temporal profile of the laser driver while keeping the laser energy constant, affecting the acceleration and injection processes. In addition, we consider how modifying the plasma channel parameters may affect the target outputs. Given the complexity and breadth of the parameter space in question, we used numerical optimization to identify high-performing configurations. In particular, we found laser profiles with additional spectral content that, when used with optimal plasma channel parameters, result in charge content an order of magnitude higher than the baseline Gaussian case while also increasing the mean energy of the electron bunch.

Physics - Plasma physics↗

Variational Dirichlet Blur Kernel Estimation

Blind image deconvolution involves two key objectives: 1) latent image and 2) blur estimation. For latent image estimation, we propose a fast deconvolution algorithm, which uses an image prior of nondimensional Gaussianity measure to enforce sparsity and an undetermined boundary condition methodology to reduce boundary artifacts. For blur estimation, a linear inverse problem with normalization and nonnegative constraints must be solved. However, the normalization constraint is ignored in many blind image deblurring methods, mainly because it makes the problem less tractable. In this paper, we show that the normalization constraint can be very naturally incorporated into the estimation process by using a Dirichlet distribution to approximate the posterior distribution of the blur. Making use of variational Dirichlet approximation, we provide a blur posterior approximation that considers the uncertainty of the estimate and removes noise in the estimated kernel. Furthermore, experiments with synthetic and real data demonstrate that the proposed method is very competitive to the state-of-the-art blind image restoration methods.

97 MATHEMATICS AND COMPUTING↗

A NON-PARAMETRIC APPROACH TO CONSTRAIN THE TRANSFER FUNCTION IN REVERBERATION MAPPING

Broad emission lines of active galactic nuclei stem from a spatially extended region (broad-line region, BLR) that is composed of discrete clouds and photoionized by the central ionizing continuum. The temporal behaviors of these emission lines are blurred echoes of continuum variations (i.e., reverberation mapping, RM) and directly reflect the structures and kinematic information of BLRs through the so-called transfer function (also known as the velocity-delay map). Based on the previous works of Rybicki and Press and Zu et al., we develop an extended, non-parametric approach to determine the transfer function for RM data, in which the transfer function is expressed as a sum of a family of relatively displaced Gaussian response functions. Therefore, arbitrary shapes of transfer functions associated with complicated BLR geometry can be seamlessly included, enabling us to relax the presumption of a specified transfer function frequently adopted in previous studies and to let it be determined by observation data. We formulate our approach in a previously well-established framework that incorporates the statistical modeling of continuum variations as a damped random walk process and takes into account long-term secular variations which are irrelevant to RM signals. The application to RM data shows the fidelity of our approach.

79 ASTRONOMY AND ASTROPHYSICS↗

Incorporation of shuttle CCT parameters in computer simulation models

Computer simulations of shuttle missions have become increasingly important during recent years. The complexity of mission planning for satellite launch and repair operations which usually involve EVA has led to the need for accurate visibility and access studies. The PLAID modeling package used in the Man-Systems Division at Johnson currently has the necessary capabilities for such studies. In addition, the modeling package is used for spatial location and orientation of shuttle components for film overlay studies such as the current investigation of the hydrogen leaks found in the shuttle flight. However, there are a number of differences between the simulation studies and actual mission viewing. These include image blur caused by the finite resolution of the CCT monitors in the shuttle and signal noise from the video tubes of the cameras. During the course of this investigation the shuttle CCT camera and monitor parameters are incorporated into the existing PLAID framework. These parameters are specific for certain camera/lens combinations and the SNR characteristics of these combinations are included in the noise models. The monitor resolution is incorporated using a Gaussian spread function such as that found in the screen phosphors in the shuttle monitors. Another difference between the traditional PLAID generated images and actual mission viewing lies in the lack of shadows and reflections of light from surfaces. Ray tracing of the scene explicitly includes the lighting and material characteristics of surfaces. The results of some preliminary studies using ray tracing techniques for the image generation process combined with the camera and monitor effects are also reported.

Huntsberger, Terry↗

Multilevel Monte Carlo for two phase flow and Buckley–Leverett transport in random heterogeneous porous media

Monte Carlo (MC) is a well known method for quantifying uncertainty arising for example in subsurface flow problems. Although robust and easy to implement, MC suffers from slow convergence. Extending MC by means of multigrid techniques yields the multilevel Monte Carlo (MLMC) method. MLMC has proven to greatly accelerate MC for several applications including stochastic ordinary differential equations in finance, elliptic stochastic partial differential equations and also hyperbolic problems. In this study, MLMC is combined with a streamline-based solver to assess uncertain two phase flow and Buckley–Leverett transport in random heterogeneous porous media. The performance of MLMC is compared to MC for a two dimensional reservoir with a multi-point Gaussian logarithmic permeability field. The influence of the variance and the correlation length of the logarithmic permeability on the MLMC performance is studied.

97 MATHEMATICS AND COMPUTING↗

AN AZIMUTHAL ASYMMETRY IN THE LkHα 330 DISK

Theory predicts that giant planets and low mass stellar companions shape circumstellar disks by opening annular gaps in the gas and dust spatial distribution. For more than a decade it has been debated whether this is the dominant process that leads to the formation of transitional disks. In this paper, we present millimeter-wave interferometric observations of the transitional disk around the young intermediate mass star LkHα 330. These observations reveal a lopsided ring in the 1.3 mm dust thermal emission characterized by a radius of about 100 AU and an azimuthal intensity variation of a factor of two. By comparing the observations with a Gaussian parametric model, we find that the observed asymmetry is consistent with a circular arc, that extends azimuthally by about 90° and emits about 1/3 of the total continuum flux at 1.3 mm. Hydrodynamic simulations show that this structure is similar to the azimuthal asymmetries in the disk surface density that might be produced by the dynamical interaction with unseen low mass companions orbiting within 70 AU from the central star. We argue that such asymmetries might lead to azimuthal variations in the millimeter-wave dust opacity and in the dust temperature, which will also affect the millimeter-wave continuum emission. Alternative explanations for the observed asymmetry that do not require the presence of companions cannot be ruled out with the existing data. Further observations of both the dust and molecular gas emission are required to derive firm conclusions on the origin of the asymmetry observed in the LkHα 330 disk.

79 ASTRONOMY AND ASTROPHYSICS↗

Presenting the Rain-Sea Interaction Facility

The new Rain-Sea Interaction Facility (RSIF) was established at GSFC/WFF and the first finds are presented. The unique feature of this laboratory is the ability to systematically study microwave scattering from a water surface roughened by artificial rain, for which the droplets are at terminal velocity. The fundamental instruments and systems (e.g., the rain simulator, scatterometers, and surface elevation probes) were installed and evaluated during these first experiments - so the majority of the data were obtained with the rain simulator at 1 m above the water tank. From these initial experiments, three new models were proposed: the square-root function for NCS vs. R, the log Gaussian model for ring-wave elevation frequency spectrum, and the Erland probability density distribution for back scattered power. Rain rate is the main input for these models, although the coefficients may be dependent upon other factors (drop-size distribution, fall velocity, radar configuration, etc.). The facility is functional and we foresee collaborative studies with investigators who are engaged in measuring and modeling rain-sea interaction processes.

Bliven, Larry F.↗

Combined iterative reconstruction and image-domain decomposition for dual energy CT using total-variation regularization

Purpose: Dual-energy CT (DECT) is being increasingly used for its capability of material decomposition and energy-selective imaging. A generic problem of DECT, however, is that the decomposition process is unstable in the sense that the relative magnitude of decomposed signals is reduced due to signal cancellation while the image noise is accumulating from the two CT images of independent scans. Direct image decomposition, therefore, leads to severe degradation of signal-to-noise ratio on the resultant images. Existing noise suppression techniques are typically implemented in DECT with the procedures of reconstruction and decomposition performed independently, which do not explore the statistical properties of decomposed images during the reconstruction for noise reduction. In this work, the authors propose an iterative approach that combines the reconstruction and the signal decomposition procedures to minimize the DECT image noise without noticeable loss of resolution. Methods: The proposed algorithm is formulated as an optimization problem, which balances the data fidelity and total variation of decomposed images in one framework, and the decomposition step is carried out iteratively together with reconstruction. The noise in the CT images from the proposed algorithm becomes well correlated even though the noise of the raw projections is independent on the two CT scans. Due to this feature, the proposed algorithm avoids noise accumulation during the decomposition process. The authors evaluate the method performance on noise suppression and spatial resolution using phantom studies and compare the algorithm with conventional denoising approaches as well as combined iterative reconstruction methods with different forms of regularization. Results: On the Catphan©600 phantom, the proposed method outperforms the existing denoising methods on preserving spatial resolution at the same level of noise suppression, i.e., a reduction of noise standard deviation by one order of magnitude. This improvement is mainly attributed to the high noise correlation in the CT images reconstructed by the proposed algorithm. Iterative reconstruction using different regularization, including quadratic orq-generalized Gaussian Markov random field regularization, achieves similar noise suppression from high noise correlation. However, the proposed TV regularization obtains a better edge preserving performance. Studies of electron density measurement also show that our method reduces the average estimation error from 9.5% to 7.1%. On the anthropomorphic head phantom, the proposed method suppresses the noise standard deviation of the decomposed images by a factor of ∼14 without blurring the fine structures in the sinus area. Conclusions: The authors propose a practical method for DECT imaging reconstruction, which combines the image reconstruction and material decomposition into one optimization framework. Compared to the existing approaches, our method achieves a superior performance on DECT imaging with respect to decomposition accuracy, noise reduction, and spatial resolution.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Thermal analysis for laser selective removal of metallic single-walled carbon nanotubes

Single-walled carbon nanotubes (SWNTs) have been envisioned as one of the best candidates for future semiconductors due to their excellent electrical properties and ample applications. However, SWNTs grow as mixture of both metallic and semiconducting tubes and this heterogeneity hampers their practical applications. Laser radiation shows promises to remove metallic SWNTs (m-SWNTs) in air under an appropriate condition. We established a scaling law, validated by finite element simulations, for the temperature rise of m-SWNTs under a pulsed laser with a Gaussian spot. It is shown that the maximum normalized m-SWNT temperature rise only depends on two non-dimensional parameters: the normalized pulse duration time and the normalized interfacial thermal resistance. In addition, the maximum temperature rise is inversely proportional to the square of spot size and proportional to the incident laser power. These results are very helpful to understand the underlying physics associated with the removal process and provides easily interpretable guidelines for further optimizations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

CXSFIT Code Application to Process Charge-Exchange Recombination Spectroscopy Data at the T-10 Tokamak

The applicability of the CXSFIT code to process experimental data from Charge-eXchange Recombination Spectroscopy (CXRS) diagnostics at the T-10 tokamak is studied with a view to its further use for processing experimental data at the ITER facility. The design and operating principle of the CXRS diagnostics are described. The main methods for processing the CXRS spectra of the 5291-Å line of C{sup 5+} ions at the T-10 tokamak (with and without subtraction of parasitic emission from the edge plasma) are analyzed. The method of averaging the CXRS spectra over several shots, which is used at the T-10 tokamak to increase the signal-to-noise ratio, is described. The approximation of the spectrum by a set of Gaussian components is used to identify the active CXRS line in the measured spectrum. Using the CXSFIT code, the ion temperature in ohmic discharges and discharges with auxiliary electron cyclotron resonance heating (ECRH) at the T-10 tokamak is calculated from the CXRS spectra of the 5291-Å line. The time behavior of the ion temperature profile in different ohmic heating modes is studied. The temperature profile dependence on the ECRH power is measured, and the dynamics of ECR removal of carbon nuclei from the T-10 plasma is described. Experimental data from the CXRS diagnostics at T-10 substantially contribute to the implementation of physical programs of studies on heat and particle transport in tokamak plasmas and investigation of geodesic acoustic mode properties.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗