SEARCH · Search NASA

Results for “CHI”

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.

SourceResearch connectionFrom yearThrough year

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

At least 109 records · Page 6

All sources NASA NTRS DOE OSTI

DOE OSTI2024

Theoretical modeling of hepatitis C acute infection in liver-humanized mice support pre-clinical assessment of candidate viruses for controlled-human-infection studies

Designing and carrying out a controlled human infection (CHI) model for hepatitis C virus (HCV) is critical for vaccine development. However, key considerations for a CHI model protocol include understanding of the earliest viral-host kinetic events during the acute phase and susceptibility of the viral isolate under consideration for use in the CHI model to antiviral treatment before any infections in human volunteers can take place. Humanized mouse models lack adaptive immune responses but provide a unique opportunity to obtain quantitative understanding of early HCV kinetics and develop mathematical models to further understand viral and innate immune response dynamics during acute HCV infection. We show that the models reproduce the measured HCV kinetics in humanized mice, which are consistent with early acute HCV-host dynamics in immunocompetent chimpanzees. Our findings suggest that humanized mice are well-suited to support development of a CHI model. In-silico and in-vivo modeling estimates provide a starting point to characterize candidate viruses for testing in CHI model studies.

Agent-based modeling↗

DOE OSTI2022

Light $$Z^\prime $$ and Dirac fermion dark matter in the $$B-L$$ model

Abstract We consider a $$U(1)_{B-L}$$ U ( 1 ) B - L model with a $$Z^\prime $$ Z ′ portal Dirac fermion dark matter (DM) $$\chi $$ χ of low mass which couples very weakly to the $$B - L$$ B - L gauge boson $$Z^\prime $$ Z ′ . An arbitrary $$B-L$$ B - L charge $$Q\ne \pm 1, \pm 3$$ Q ≠ ± 1 , ± 3 of the DM $$\chi $$ χ ensures its stability. Motivated by the sensitivity reach of forthcoming “Lifetime Frontier” experiments, we focus on the $$Z^\prime $$ Z ′ mass, $$m_{Z^\prime }$$ m Z ′ , in the sub-GeV to few GeV range. To evaluate the DM relic abundance, we examine both the freeze-out and freeze-in DM scenarios. For the freeze-out scenario, we show that the observed DM abundance is reproduced near the $$Z^\prime $$ Z ′ resonance, $$m_\chi \simeq m_{Z^\prime }/2$$ m χ ≃ m Z ′ / 2 , where $$m_\chi $$ m χ is the DM mass. For the freeze-in scenario, we focus on $$m_\chi \ll m_{Z^\prime }$$ m χ ≪ m Z ′ . We show that for a fixed value of $$m_{Z^\prime }$$ m Z ′ , $$g_{BL}$$ g BL values roughly scale as 1/ Q to reproduce the observed DM abundance. For various Q values in the range between $$10^{-6}$$ 10 - 6 and $$10^2$$ 10 2 , we show that the gauge coupling values $$g_{BL}$$ g BL needed to reproduce the observed DM abundance lie in the search reach of future planned and/or proposed experiments such as FASER, Belle-II, LDMX, and SHiP. In the freeze-in case, the Q values to realize observable $$g_{BL}$$ g BL values are found to be much smaller than that in the freeze-out case.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

NASA NTRS1989

A finite difference treatment of Stokes-type flows: Preliminary report

The equations Laplacian operator omega = 0, (1.1a) and omega = Laplacian operator Chi, (1.1b) describe, in suitable units, 2-D Stokes flow of an incompressible fluid occupying a domain D in which omega is the vorticity and Chi is the stream function. The flow is uniquely determined by specifying the velocity on the boundary B of D, a condition which leads to specifying the stream function Chi and its normal derivative Chi sub n on B. A mathematically similar problem arises in describing the equilibrium of a flat plate in structural mechanics where a related 1-D problem by finite difference or finite element methods is to introduce effective methods for imposing the boundary conditions through which (1.1a) is coupled to (1.1b). These models thus provide a simple starting point for examining the general treatment of boundary conditions for more general time dependent Navier-Stokes incompressible flows. For the purpose of discussion it is assumed that D is a square domain. A standard finite difference method to solve (1.1) is to introduce a uniform grid and then use standard five point finite difference operators to express each equation in (1.1). At any point on the boundary B a value of Chi is specified by the boundary conditions but a value of omega at the same boundary mesh point will also be required to complete the computation. Methods are discussed which overcome the difficulty in solving these problems.

Rose, M. E.↗

NASA NTRS1996

Variance Analysis of Unevenly Spaced Time Series Data

We have investigated the effect of uneven data spacing on the computation of delta (sub chi)(gamma). Evenly spaced simulated data sets were generated for noise processes ranging from white phase modulation (PM) to random walk frequency modulation (FM). Delta(sub chi)(gamma) was then calculated for each noise type. Data were subsequently removed from each simulated data set using typical two-way satellite time and frequency transfer (TWSTFT) data patterns to create two unevenly spaced sets with average intervals of 2.8 and 3.6 days. Delta(sub chi)(gamma) was then calculated for each sparse data set using two different approaches. First the missing data points were replaced by linear interpolation and delta (sub chi)(gamma) calculated from this now full data set. The second approach ignored the fact that the data were unevenly spaced and calculated delta(sub chi)(gamma) as if the data were equally spaced with average spacing of 2.8 or 3.6 days. Both approaches have advantages and disadvantages, and techniques are presented for correcting errors caused by uneven data spacing in typical TWSTFT data sets.

Hackman, Christine↗

NASA NTRS1999

A Spherical Earth Solution for TOA Lightning Location Retrieval

With the advent of high-speed digital computers, the method of chi square minimization is a highly practical means for analyzing a wide variety of (otherwise intractable) nonlinear inversion problems in applied mathematical physics. Little thought or effort is required to apply the chi square method to obtain quick and reasonable estimates of a solution, and the method offers a means to assess retrieval errors. Because the method is simple and practical it is sometimes hastily applied to problems that can be solved by formal analytic or quasi-analytic means. Presently, Global Atmospherics Inc. (GAI) finds the minimum of a chi square function to analyze time-of-arrival (TOA) and magnetic bearing data derived from the National Lightning Detection Network (NLDN); ellipsoidal Earth geometry is assumed. An analytic solution to this problem has not yet been attained, but the consideration and solving of less general problem statements might eventually lead to a final solution. In the present study, the problem of retrieving lightning ground-strike location on a spherical Earth surface using a network of 4 or more time-of-arrival (TOA) sensors is considered. It is shown that this problem has an analytic solution and therefore does not require the use of nonlinear estimation theory (such as the chi square method mentioned above). The mathematical robustness of the analytic solution is tested using computer-generated lightning sources and simulated TOA measurement errors. A quasi-analytic extension of the spherical Earth solution for an oblate spheroidal Earth geometry is considered in a related study. The incorporation of magnetic bearing information into these analytic solutions would lead to a general and elegant analytic retrieval scheme that would most likely replace the chi square estimation theory currently employed by Global Atmospherics Inc. (GAI).

Koshak, William J,↗

DOE OSTI2023

Discovering and Designing a Chimeric Hyperthermophilic Chitinase for Crystalline Chitin Degradation

Chitin is one of the most abundant renewable biopolymers on earth. However, it is highly crystalline and recalcitrant to degrade. Here, we report a hyperthermophilic chitinase (ActChi) to directly hydrolyze crystalline chitin at its optimal temperature of 80 °C. It contains a malectin domain, a fibronectin type-III (Fn3) domain, and a catalytic domain (CD chi ). Both Fn3 and malectin have the function of chitin binding domain (ChBD) to increase the activity. Fn3 also significantly increases thermostability, but malectin decreases it. To enhance both activity and thermostability, here we introduced a heterogeneous and hyperthermophilic ChBD at the N-terminus of CD chi to obtain ChBD-CD chi . The activity of this hybrid enzyme is 201 U/μmol for crystalline chitin, which has increased 400% compared with that of ActChi. In addition, ChBD-CD chi can continuously degrade crystalline chitin for more than 4 days at 70 °C to increase the overall hydrolysis rate. The strategy is a good example of green sustainable degradation for crystalline biopolymer in nature.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

The role of electric dominance for particle injection in relativistic reconnection

ABSTRACT Magnetic reconnection in relativistic plasmas – where the magnetization $\sigma \gg 1$ – is regarded as an efficient particle accelerator, capable of explaining the most dramatic astrophysical flares. We employ two-dimensional (2D) particle-in-cell simulations of relativistic pair-plasma reconnection with vanishing guide field and outflow boundaries to quantify the impact of the energy gain occurring in regions of electric dominance ($E\gt B$) for the early stages of particle acceleration (i.e. the ‘injection’ stage). Given an injection threshold energy $\epsilon ^\ast =\sigma /4$ for the particles that eventually attain energy $\epsilon _{\rm T}$ by time T, we calculate the mean fractional contribution $\zeta (\epsilon ^\ast ,\epsilon _{\rm T})$ by $E\gt B$ fields to particle energization at the time when the threshold $\epsilon ^\ast$ is reached. We find that $\zeta$ monotonically increases with $\sigma$ and $\epsilon _{\rm T}$; for $\sigma \gtrsim 50$ and $\epsilon _{\rm T}/\sigma \gtrsim 8$, we find that $\gtrsim 80~{{\ \rm per\ cent}}$ of the energy gain obtained before reaching $\epsilon ^\ast =\sigma /4$ occurs in $E\gt B$ regions. We find that $\zeta$ is independent of simulation box size $L_x$, as long as $\epsilon _{\rm T}$ is normalized to the maximum particle energy, which scales as $\epsilon _{\rm max}\propto L_{\rm x}^{1/2}$ in 2D. The distribution of energy gains $\epsilon _{\chi }$ acquired in $E\gt B$ regions can be modelled as $\mathrm{ d}N/\mathrm{ d}\epsilon _{\chi }\propto \epsilon _{\chi }^{-0.35}\exp [-(\epsilon _{\chi }/0.06\, \sigma)^{0.5}]$. Our results help assess the role of electric dominance in relativistic reconnection with vanishing guide fields, which is realized in the magnetospheres of black holes and neutron stars.

Gupta, Sanya (ORCID:0000000151944384)↗

DOE OSTI2024

Nucleon Spin Sum Rules and Spin Polarizabilities at low $Q^2$

We report on recently published experimental studies on spin sum rules, namely the generalized Gerasimov-Drell-Hearn, Bjorken, Burkhardt-Cottingham, Schwinger, and generalized spin polarizability sum rules. The data were taken at Jefferson Lab in Halls A and B by experiments E97-110 and EG4, respectively. They covered the very low $Q^2$ domain, down to $Q^2$$\simeq$0.02 GeV$^2$, where Chiral Effective Field Theory ($\chi$EFT) predictions should be valid. While some of the obervables agree with the state-of-the-art $\chi$EFT predictions, others are in tensions, including the Longitudinal-Transverse interference polarizability $\delta_{\rm LT}^n(Q^2)$ for which $\chi$EFT prediction was expected to be robust. This suggests that $\chi$EFT does not yet consistently describes nucleon spin observables, even in the very low $Q^2$ domain covered by the experiments.

Deur, Alexandre↗

DOE OSTI2026

Diquark scalar production of a vectorlike quark pair at the LHC

We study the discovery potential of LHC experiments for resonantly produced vectorlike quarks ($\chi$) when the $s$-channel resonance is an ultraheavy scalar diquark ($S_{uu}$) with a mass in the 7–8.5 TeV range. Focusing on the process $pp \rightarrow S_{uu} \rightarrow \chi \chi \rightarrow (W^+b)(W^+b)$ we target the fully hadronic decay mode of both $W^+$ bosons, resulting in a six-jet final state. Signal–background separation is performed using Machine Learning algorithms trained to construct a multidimensional classifier. Our results show that ATLAS or CMS searches in this channel, with an integrated luminosity of 3000 fb$^{-1}$, could discover or exclude a scalar diquark with a mass near 8 TeV even for a relatively small Yukawa coupling to up quarks, $y_{uu} \simeq 0.2$. We also present preliminary studies of the four-jet final state from $pp \rightarrow S_{uu} \rightarrow u \chi \rightarrow u (W^+b)$.

Duminica, Ioana [Bucharest, IFIN-HH; Bucharest U.]↗

DOE OSTI2022

Rheology and Gelation of Hyaluronic Acid/Chitosan Coacervates

Hyaluronic acid (HA) and chitosan (CHI) are biopolyelectrolytes which are interesting for both the medical and polymer physics communities due to their biocompatibility and semi-flexibility, respectively. In this work, we demonstrate by rheology experiments that the linear viscoelasticity of HA/CHI coacervates depends strongly on the molecular weight of the polymers. Moduli for coacervates were found significantly higher than those of individual HA and CHI physical gels. A remarkable 1.5-fold increase in moduli was noted when catechol-conjugated HA and CHI were used instead. This was attributed to the conversion of coacervates to chemical gels by oxidation of 3,4-dihydroxyphenylalanine (DOPA) groups in HA and CHI to di-DOPA crosslinks. These rheological results put HA/CHI coacervates in the category of strong candidates as injectable tissue scaffolds or medical adhesives.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Tuning Solvation Dynamics of Electrolytes at Their Eutectic Point Through Halide Identity

Deep eutectic solvents (DESs) are regarded as highly promising solvent systems for redox flow batteries. DESs, composed of choline halides (ChX, X = F − , Cl − , Br − , I − ) and ethylene glycol (EG), exhibit distinct physicochemical properties at their eutectic points, including halide-dependent phase behavior, viscosity, polarity, conductivity, and solvation dynamics. In this study, we investigate the effects of the halide identity on the solvation properties of ChX:EG mixtures at varying mol % of ChX salt content. The solvatochromic polarity based on ET(30) measurements indicates higher polarity for larger halides (I − > Br − ) than for smaller halides (Cl − > F − ), which exhibit larger compensating solvation shells. The ionic conductivity follows the trend of the solvent fluidity (the inverse of the viscosity), namely ChCl > ChBr > ChI > ChF, influenced by the ion mobility and solvodynamic radii. Measurements of the liquidus temperatures (T L ) reveal that the system with ChCl exhibits the deepest eutectic point (at ~20 mol % ChCl), while ChBr and ChI have shallower minima at ~10 mol % ChBr and ~3 mol % ChI, respectively. ChF does not display a eutectic transition but instead appears to readily supercool at salt concentrations above 30 mol % ChF. Consistent with the phase transition measurements, femtosecond transient absorption spectroscopy shows that in the ChCl system, the solvation dynamics become faster with an increasing salt concentration up to ~16.67 mol %, after which the dynamics slow down with further increases in the salt content. The ChF-based system exhibits similar behavior, though with slower dynamics. In contrast, the solvation dynamics of the systems containing ChBr and ChI monotonously slow down with an increasing salt concentration, in agreement with the phase transition measurements, which show that the eutectic points occur at low salt concentrations. These measurements suggest that the solvent composition and, in particular, the identity of the halide anion play a significant role in the solvation behavior of these ethylene-glycol-based DESs, offering a foundation for tuning the DES properties for specific applications.

Biochemistry & Molecular Biology↗

NASA NTRS1987

Statistical tests for closure of plate motion circuits

Two methods, one based on a chi-square test and the second on an F-ratio test, of testing for the plate motion circuit closures are described and evaluated. The chi-square test is used to evaluate goodness of fit, and it is assumed that the assigned errors are accurate estimates of the true errors in the data. The F-ratio test is used to compare variances of distributions, and it is assumed that the relative values of assigned error are accurate. The two methods are applied to the data of Minster and Jordan (1978) on the motion of the three plates that meet at the Galapagos Triple Junction, and the motion of the three plates that meet at the Indian Ocean Triple Junction. It is noted that the F-ratio plate circuit closure test is more useful than the chi-square test for identifying systematic misfits in data because the chi-square test overestimates the errors of plate motion data.

Gordon, Richard G.↗

NASA NTRS1992

A model for inferring canopy and underlying soil temperatures from multi-directional measurements

A model is presented of thermal emission from a canopy/soil surface, where the soil and the leaves are at different temperatures, Tg and Tc, respectively. The temperature Tm corresponding to a radiometer reading is given by B sub lambda(Tm) = chi-B sub lambda(Tg) + (1-chi)B sub lambda(Tc), where B sub lambda denotes the Planck blackbody function at wavelength lambda, chi specifies the fraction of the field of view occupied by the soil at a given view direction, and an emissivity of 1.0 is assumed for the plants and the soil. It is observed that at large view zenith angles, only the plants are effectively seen, and therefore Tc can be determined from the observations a large zenith angles, to the extent that such observations are practical. Water stress can produce an increase of chi and thus tends to produce an exaggerated increase in the observed temperature compared to the actual increase in canopy temperature. These effects are analyzed for a simulated soybean canopy.

Otterman, J.↗

NASA NTRS1994

Comparison of theories for gravity wave induced fluctuations in airglow emissions

A comparison is undertaken of theories for the gravity wave induced fluctuations in the intensity of airglow emissions and the associated temperature of the source region. The comparison is made in terms of Krassovsky's ratio eta(sub E) for a vertically extended emission region (eta(sub E) is the ratio of the vertically integrated normalized intensity perturbation to the vertically integrated intensity-weighted temperature perturbation). It is shown that the formulas for eta(sub E) in the works by Tarasick and Hines (1990) and Schubert et al. (1991) are in agreement for the case of an inviscid atmosphere. The calculation of eta(sub E) using the theory of Tarasick and Hines (1990) requires determination of their function chi; we show that chi is simply related to the 'single-level' Krassovsky's ratio eta of Schubert et al. (1991). The general relationship between chi and eta is applied to a simple chemical-dynamical model of the O2 atmospheric airglow and the altitude dependence of these quantities is evaluated for nonsteady state chemistry. Though the Tarasick and Hines (1990) formula for eta(sub E) does not explicitly depend on the scale heights of the minor constituents involved in airglow chemistry, eta(sub E) implicitly depends upon these scale heights through its dependences on chemical production and loss contained in chi. We demonstrate this dependence of eta(sub E) for the OH nightglow on atomic oxygen scale height by direct numerical evaluation of eta(sub E) in this case the dependence originates in the chemical production of perturbed ozone.

Walterscheid, R. L.↗

NASA NTRS2005

Improvement of Raman lidar algorithm for quantifying aerosol extinction

Aerosols are particles of different composition and origin and influence the formation of clouds which are important in atmospheric radiative balance. At the present there is high uncertainty on the effect of aerosols on climate and this is mainly due to the fact that aerosol presence in the atmosphere can be highly variable in space and time. Monitoring of the aerosols in the atmosphere is necessary to better understanding many of these uncertainties. A lidar (an instrument that uses light to detect the extent of atmospheric aerosol loading) can be particularly useful to monitor aerosols in the atmosphere since it is capable to record the scattered intensity as a function of altitude from molecules and aerosols. One lidar method (the Raman lidar) makes use of the different wavelength changes that occur when light interacts with the varying chemistry and structure of atmospheric aerosols. One quantity that is indicative of aerosol presence is the aerosol extinction which quantifies the amount of attenuation (removal of photons), due to scattering, that light undergoes when propagating in the atmosphere. It can be directly measured with a Raman lidar using the wavelength dependence of the received signal. In order to calculate aerosol extinction from Raman scattering data it is necessary to evaluate the rate of change (derivative) of a Raman signal with respect to altitude. Since derivatives are defined for continuous functions, they cannot be performed directly on the experimental data which are not continuous. The most popular technique to find the functional behavior of experimental data is the least-square fit. This procedure allows finding a polynomial function which better approximate the experimental data. The typical approach in the lidar community is to make an a priori assumption about the functional behavior of the data in order to calculate the derivative. It has been shown in previous work that the use of the chi-square technique to determine the most likely functional behavior of the data prior to actually calculating the derivative eliminates the need for making a priori assumptions. We note that the a priori choice of a model itself can lead to larger uncertainties as compared to the method that is validated here. In this manuscript, the chi-square technique that determines the most likely functional behavior is validated through numerical simulation and by application to a large body of Raman lidar measurements. In general, we show that the chi-square approach to evaluate aerosol extinction yields lower extinction uncertainty than the traditional technique. We also use the technique to study the feasibility of developing a general characterization of the extinction uncertainty that could permit the uncertainty in Raman lidar aerosol extinction measurements to be estimated accurately without the use of the chi-square technique.

Russo, Felicita↗

NASA NTRS2008

Experimental and Coupled-channels Investigation of the Radiative Properties of the N2 c4 (sup 1)Sigma+(sub u) - X (sup 1)Sigma+(sub g) Band System

The emission properties of the N2 c(sup prime)(sub 4) (sup 1)Sigma+(sub u) - Chi (sup 1)Sigma+(sub g) band system have been investigated in a joint experimental and coupled-channels theoretical study. Relative intensities of the c(sup prime)(sub 4) (sup 1)Sigma+(sub u)(0) - Chi (sup 1)Sigma+(sub g)(v(sub i)) transitions, measured via electron-impact-induced emission spectroscopy, are combined with a coupled-channel Schroedinger equation (CSE) model of the N2 molecule, enabling determination of the diabatic electronic transition moment for the c(sup prime)(sub 4) (sup 1)Sigma+(sub u) - Chi (sup 1)Sigma+(sub g) system as a function of internuclear distance. The CSE probabilities are further verified by comparison with a high-resolution experimental spectrum. Spontaneous transition probabilities of the c(sup prime)(sub 4) (sup 1)Sigma+(sub u) - Chi (sup 1)Sigma+(sub g) modeling atmospheric emission, can now be calculated reliably.

internuclear distance↗

NASA NTRS2011

Titan's Aerosol and Stratospheric Ice Opacities Between 18 and 500 Micrometers: Vertical and Spectral Characteristics from Cassini CIRS

Vertical distributions and spectral characteristics of Titan's photochemical aerosol and stratospheric ices are determined between 20 and 560 per centimeter (500-18 micrometers) from the Cassini Composite Infrared Spectrometer (CIRS). Results are obtained for latitudes of 15 N, 15 S, and 58 S, where accurate temperature profiles can be independently determined. In addition, estimates of aerosol and ice abundances at 62 N relative to those at 15 S are derived. Aerosol abundances are comparable at the two latitudes, but stratospheric ices are approximately 3 times more abundant at 62 N than at 15 S. Generally, nitrile ice clouds (probably HCN and HC3N), as inferred from a composite emission feature at approximately 160 per centimeter, appear to be located over a narrow altitude range in the stratosphere centered at approximately 90 km. Although most abundant at high northern latitudes, these nitrile ice clouds extend down through low latitudes and into mid southern latitudes, at least as far as 58 S. There is some evidence of a second ice cloud layer at approximately 60 km altitude at 58 S associated with an emission feature at approximately 80 per centimeter. We speculate that the identify of this cloud may be due to C2H6 ice, which in the vapor phase is the most abundant hydrocarbon (next to CH4) in the stratosphere of Titan. Unlike the highly restricted range of altitudes (50-100 km) associated with organic condensate clouds, Titan's photochemical aerosol appears to be well-mixed from the surface to the top of the stratosphere near an altitude of 300 km, and the spectral shape does not appear to change between 15 N and 58 S latitude. The ratio of aerosol-to-gas scale heights range from 1.3-2.4 at about 160 km to 1.1-1.4 at 300 km, although there is considerable variability with latitude, The aerosol exhibits a very broad emission feature peaking at approximately 140 per centimeter. Due to its extreme breadth and low wavenumber, we speculate that this feature may be caused by low-energy vibrations of two-dimensional lattice structures of large molecules. Examples of such molecules include polycyclic aromatic hydrocarbons (PAHs) and nitrogenated aromatics. Finally, volume extinction coefficients N chi EPSILON derived from 15 S CIRS data at a wavelength of lambda = 62.5 micrometers are compared with those derived from the 10 S Huygens Descent Imager/Spectral Radiometer (DISR) data at 1.583 micrometers. This comparison yields volume extinction coefficient ratios N chi EPSILON (1.583 micrometers)/N chi EPSILON (62.5 micrometers) of roughly 70 and 20, respectively, for Titan's aerosol and stratospheric ices, The inferred particle cross-section ratios chi EPSILON(1.583 micrometers)/chi EPSILON (62.5 micrometers) appear to be consistent with sub-micron size aerosol particles, and effective radii of only a few microns for stratospheric ice cloud particles.

Anderson, Carrie M.↗

NASA NTRS2017

Data Reduction Pipeline for the CHARIS Integral-Field Spectrograph I: Detector Readout Calibration and Data Cube Extraction

We present the data reduction pipeline for CHARIS, a high-contrast integral-field spectrograph for the Subaru Telescope. The pipeline constructs a ramp from the raw reads using the measured nonlinear pixel response and reconstructs the data cube using one of three extraction algorithms: aperture photometry, optimal extraction, or chi-squared fitting. We measure and apply both a detector flatfield and a lenslet flatfield and reconstruct the wavelength- and position-dependent lenslet point-spread function (PSF) from images taken with a tunable laser. We use these measured PSFs to implement a chi-squared-based extraction of the data cube, with typical residuals of approximately 5 percent due to imperfect models of the under-sampled lenslet PSFs. The full two-dimensional residual of the chi-squared extraction allows us to model and remove correlated read noise, dramatically improving CHARIS's performance. The chi-squared extraction produces a data cube that has been deconvolved with the line-spread function and never performs any interpolations of either the data or the individual lenslet spectra. The extracted data cube also includes uncertainties for each spatial and spectral measurement. CHARIS's software is parallelized, written in Python and Cython, and freely available on github with a separate documentation page. Astrometric and spectrophotometric calibrations of the data cubes and PSF subtraction will be treated in a forthcoming paper.

Brandt, Timothy D.↗