Search NASA⌕ Search

SEARCH · Search NASA

Results for “quantum information and processing”

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,405 records · Page 78

Dimensionality and noise in energy selective x-ray imaging

Purpose: To develop and test a method to quantify the effect of dimensionality on the noise in energy selective x-ray imaging.Methods: The Cramèr-Rao lower bound (CRLB), a universal lower limit of the covariance of any unbiased estimator, is used to quantify the noise. It is shown that increasing dimensionality always increases, or at best leaves the same, the variance. An analytic formula for the increase in variance in an energy selective x-ray system is derived. The formula is used to gain insight into the dependence of the increase in variance on the properties of the additional basis functions, the measurement noise covariance, and the source spectrum. The formula is also used with computer simulations to quantify the dependence of the additional variance on these factors. Simulated images of an object with three materials are used to demonstrate the trade-off of increased information with dimensionality and noise. The images are computed from energy selective data with a maximum likelihood estimator.Results: The increase in variance depends most importantly on the dimension and on the properties of the additional basis functions. With the attenuation coefficients of cortical bone, soft tissue, and adipose tissue as the basis functions, the increase in variance of the bone component from two to three dimensions is 1.4 × 10{sup 3}. With the soft tissue component, it is 2.7 × 10{sup 4}. If the attenuation coefficient of a high atomic number contrast agent is used as the third basis function, there is only a slight increase in the variance from two to three basis functions, 1.03 and 7.4 for the bone and soft tissue components, respectively. The changes in spectrum shape with beam hardening also have a substantial effect. They increase the variance by a factor of approximately 200 for the bone component and 220 for the soft tissue component as the soft tissue object thickness increases from 1 to 30 cm. Decreasing the energy resolution of the detectors increases the variance of the bone component markedly with three dimension processing, approximately a factor of 25 as the resolution decreases from 100 to 3 bins. The increase with two dimension processing for adipose tissue is a factor of two and with the contrast agent as the third material for two or three dimensions is also a factor of two for both components. The simulated images show that a maximum likelihood estimator can be used to process energy selective x-ray data to produce images with noise close to the CRLB.Conclusions: The method presented can be used to compute the effects of the object attenuation coefficients and the x-ray system properties on the relationship of dimensionality and noise in energy selective x-ray imaging systems.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Reliability of unstable periodic orbit based control strategies in biological systems

Presence of recurrent and statistically significant unstable periodic orbits (UPOs) in time series obtained from biological systems is now routinely used as evidence for low dimensional chaos. Extracting accurate dynamical information from the detected UPO trajectories is vital for successful control strategies that either aim to stabilize the system near the fixed point or steer the system away from the periodic orbits. A hybrid UPO detection method from return maps that combines topological recurrence criterion, matrix fit algorithm, and stringent criterion for fixed point location gives accurate and statistically significant UPOs even in the presence of significant noise. Geometry of the return map, frequency of UPOs visiting the same trajectory, length of the data set, strength of the noise, and degree of nonstationarity affect the efficacy of the proposed method. Results suggest that establishing determinism from unambiguous UPO detection is often possible in short data sets with significant noise, but derived dynamical properties are rarely accurate and adequate for controlling the dynamics around these UPOs. A repeat chaos control experiment on epileptic hippocampal slices through more stringent control strategy and adaptive UPO tracking is reinterpreted in this context through simulation of similar control experiments on an analogous but stochastic computer model of epileptic brain slices. Reproduction of equivalent results suggests that far more stringent criteria are needed for linking apparent success of control in such experiments with possible determinism in the underlying dynamics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Theoretically informed Monte Carlo simulation of liquid crystals by sampling of alignment-tensor fields

A theoretically informed coarse-grained Monte Carlo method is proposed for studying liquid crystals. The free energy functional of the system is described in the framework of the Landau-de Gennes formalism. The alignment field and its gradients are approximated by finite differences, and the free energy is minimized through a stochastic sampling technique. The validity of the proposed method is established by comparing the results of the proposed approach to those of traditional free energy minimization techniques. Its usefulness is illustrated in the context of three systems, namely, a nematic liquid crystal confined in a slit channel, a nematic liquid crystal droplet, and a chiral liquid crystal in the bulk. It is found that for systems that exhibit multiple metastable morphologies, the proposed Monte Carlo method is generally able to identify lower free energy states that are often missed by traditional approaches. Importantly, the Monte Carlo method identifies such states from random initial configurations, thereby obviating the need for educated initial guesses that can be difficult to formulate.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Out-of-Distribution Generalization for Learning Quantum Channels with Low-Energy Coherent States

When experimentally learning the action of a continuous-variable quantum process by probing it with inputs, there will often be some restriction on the input states used. One experimentally simple way to probe a quantum channel is to use low-energy coherent states. Learning a quantum channel in this way presents difficulties, due to the fact that two channels may act similarly on low-energy inputs but very differently for high-energy inputs. They may also act similarly on coherent-state inputs but differently on nonclassical inputs. Extrapolating the behavior of a channel for more general input states from its action on the far more limited set of low-energy coherent states is a case of out-of-distribution generalization. To be sure that such generalization gives meaningful results, one needs to relate error bounds for the training set to bounds that are valid for all inputs. We show that for any pair of channels that act sufficiently similarly on low-energy coherent-state inputs, one can bound how different the input-output relations are for any (high-energy or highly nonclassical) input. This proves that out-of-distribution generalization is always possible for learning quantum channels using low-energy coherent states, as long as enough samples are used.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Bayesian 3D X-ray computed tomography image reconstruction with a scaled Gaussian mixture prior model

In order to improve quality of 3D X-ray tomography reconstruction for Non Destructive Testing (NDT), we investigate in this paper hierarchical Bayesian methods. In NDT, useful prior information on the volume like the limited number of materials or the presence of homogeneous area can be included in the iterative reconstruction algorithms. In hierarchical Bayesian methods, not only the volume is estimated thanks to the prior model of the volume but also the hyper parameters of this prior. This additional complexity in the reconstruction methods when applied to large volumes (from 512{sup 3} to 8192{sup 3} voxels) results in an increasing computational cost. To reduce it, the hierarchical Bayesian methods investigated in this paper lead to an algorithm acceleration by Variational Bayesian Approximation (VBA) [1] and hardware acceleration thanks to projection and back-projection operators paralleled on many core processors like GPU [2]. In this paper, we will consider a Student-t prior on the gradient of the image implemented in a hierarchical way [3, 4, 1]. Operators H (forward or projection) and H{sup t} (adjoint or back-projection) implanted in multi-GPU [2] have been used in this study. Different methods will be evalued on synthetic volume 'Shepp and Logan' in terms of quality and time of reconstruction. We used several simple regularizations of order 1 and order 2. Other prior models also exists [5]. Sometimes for a discrete image, we can do the segmentation and reconstruction at the same time, then the reconstruction can be done with less projections.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Comparative Density Functional Theory Study of Magnetic Exchange Couplings in Dinuclear Transition-Metal Complexes

Multicenter transition-metal complexes (MCTMs) with magnetically interacting ions have been proposed as components for information-processing devices and storage units. For any practical application of MCTMs as magnetic units, it is crucial to characterize their magnetic behavior, and in particular, the isotropic magnetic exchange coupling, J, between its magnetic centers. Due to the large size of typical MCTMs, density functional theory is the only practical electronic structure method for evaluating the J coupling. Here, we assess the accuracy of different density functional approximations for predicting the magnetic couplings of eight dinuclear transition-metal complexes, including five dimanganese, two dicopper, and one divanadium with known reliable experimental J couplings spanning from ferromagnetic to strong antiferromagnetic. The density functionals considered include global hybrid functionals which mix semilocal density functional approximations and exact exchange with a fixed admixing parameter, six local hybrid functionals where the admixing parameters are extended to be spatially dependent, the SCAN and r 2 SCAN meta-generalized gradient approximations (GGAs), and two widely used GGAs. We found that global hybrids tested in this work have a tendency to over-correct the error in magnetic coupling parameters from the Perdew–Burke–Ernzerhof (PBE) GGA as seen for manganese complexes. The performance of local hybrid density functionals shows no improvement in terms of bias and is scattered without a clear trend, suggesting that more efforts are needed for the extension from global to local hybrid density functionals for this particular property. The SCAN and r 2 SCAN meta-GGAs are found to perform as well as benchmark global hybrids on most tested complexes. We further analyze the charge density redistribution of meta-GGAs as well as global and local hybrid density functionals with respect to that of PBE, in connection to the self-interaction error or delocalization error.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Time-resolved interferometric measurements of plasma density evolution in laser-driven capacitor-coil targets

Laser-driven capacitor-coil targets provide a compact platform for generating strong magnetic fields and are widely used in magnetized high-energy-density (HED) plasma experiments. In addition to magnetic-field generation, these targets also produce plasma in the coil region, which can influence the subject physical processes and interact with secondary targets or external plasmas in their applications. However, direct, time-resolved measurements of the plasma density surrounding the coil remain limited. Here, we report interferometric measurements of the plasma density evolution in laser-driven capacitor-coil targets irradiated by the University of Osaka LFEX laser. Two-dimensional electron density maps reveal two distinct plasma sources loading the coil region: plasma generated in the coil itself and plasma produced by laser ablation of the target plates. These results provide quantitative information on plasma loading and evolution in capacitor-coil targets and are directly relevant to the design and modeling of magnetized HED plasma experiments.

Zhang, Yang [Princeton University, NJ (United Stat↗

Probing Linewidths and Biexciton Quantum Yields of Single Cesium Lead Halide Nanocrystals in Solution

Cesium lead halide (CsPbX 3 , X = Cl, Br, I) perovskite nanocrystals (PNCs) have recently become a promising material for optoelectronic applications due to their high emission quantum yields and facile band gap tunability via both halide composition and size. The spectroscopy of single PNCs enhances our understanding of the effect of confinement on excitations in PNCs in the absence of obfuscating ensemble averaging and can also inform synthetic efforts. However, single PNC studies have been hampered by poor PNC photostability under confocal excitation, precluding interrogation of all but the most stable PNCs, and leading to a lack of understanding of PNCs in the regime of high confinement. Here, we report the first comprehensive spectroscopic investigation of single PNC properties using solution-phase photon-correlation methods, including both highly confined and blue-emitting PNCs, previously inaccessible to single NC techniques. With minimally perturbative solution-phase photon-correlation Fourier spectroscopy (s-PCFS), we establish that the ensemble emission linewidth of PNCs of all sizes and compositions is predominantly determined by the intrinsic single PNC linewidth (homogeneous broadening). The single PNC linewidth, in turn, dramatically increases with increasing confinement, consistent with what has been found for II–VI semiconductor nanocrystals. With solution-phase photon antibunching measurements, we survey the biexciton-to-exciton quantum yield ratio (BX/X QY) in the absence of user-selection bias or photodegradation. Remarkably, the BX/X QY ratio depends both on the PNC size and halide composition, with values between ~2% for highly confined bromide PNCs and ~50% for intermediately confined iodide PNCs. Our results suggest a wide range of underlying Auger rates, likely due to transitory charge carrier separation in PNCs with relaxed confinement.

single nanocrystal spectroscopy↗

Purification of Zirconium Tetrachloride from UNF Cladding

The overarching objective of this project was to design and demonstrate (laboratory scale) a sublimation-based purification protocol to obtain pure zirconium chloride (ZrCl 4 ) from impure feeds derived from the direct chlorination of Zircaloy cladding materials. To achieve this objective the following goals were identified. (1) Design, construct and test an apparatus capable of subliming and purifying ZrCl 4 at room pressure and temperatures between 600 and 700K (sublimation point of ZrCl 4 606 K). (2) Investigate reactions that may be used to chemically alter the makeup of the impure ZrCl 4 feed materials to enhance their removal from ZrCl 4 in the sublimation-purification process. (3) Computationally model the reactions of metal chlorides in the gas phase under sublimation conditions (1 atm, 600 K) to determine if oligomers are formed and, if formed, how they would influence sublimation-based purification strategies. All experimental investigations performed as part of this project used stable (non-radioactive) compounds and materials as surrogates for the impure ZrCl 4 feedstocks that would be derived from the direct chlorination of Zircaloy claddings from spent nuclear fuel rods. No radiologically active materials were used at any time in this project. The impure ZrCl 4 feedstocks used in this study were derived from the direct chlorination of virgin Zircaloy claddings (Zirc-2, -4 and high Nb). The main goals of the project have been achieved during the three year period of this NEUP sponsored project. The impurities of greatest concern in the feedstocks were identified as the chloride salts of the alloying metals originally found in Zircaloy (Fe, Nb) and radiological contaminants (Cs, Sr, and Sb). A series of prototype sublimation apparatuses and heating arrangements were designed, constructed and tested. A general batch-sublimator design was identified and a 1.0 kg batch sublimator/oven system was fabricated from glass. Purification protocols developed on small scale runs were found to scale to 1.0 kg batches without problem. Purification factor ranges (PFs: wt% final / wt% initial ) for the impurity chlorides identified above have been measured: Fe (40 – 500); Nb (10 – 70); Cs (520 – 5000); Sr (250 – 1900); Sb (220 – 2300). Removal of iron is sensitive to the partial pressure of hydrogen (in nitrogen) flowing through the system. Lower PFs are obtained when the hydrogen partial pressure is lower. Lower flows are required to raise the iron PF. Removal of niobium is very sensitive to the endpoint determination of the forerun. When judged properly (color change yellow to white in sublimate) reasonable PFs (500) can be achieved. However, as discussed in the conclusions section, endpoint determination can be difficult and significantly lower PFs may be obtained when the endpoint is not correctly identified. Several methods of obtaining quantitative cesium, strontium and antimony PFs were investigated, the best being ICP-MS. The results reported here reflect very limited access to and ICP-MS instrument during the course of this project. It should be noted that, while significantly reduced amounts of antimony are observed in the product, antimony is also found in the exhaust trap after the condenser. Quantum chemical computational modeling studies of the reactions of metal chlorides in the gas phase provided important information about a chemical transport mechanism that could defeat sublimation-based purification strategies. Many metal chlorides are known to oligomerize in the gas phase to reduce their coordinative unsaturation. Dimerization of ZrCl 4 with iron or niobium chlorides would provide a mechanism by which both metal impurities could be carried over with the product during sublimation of the bulk material. Computational modeling conclusively showed that under sublimation conditions in the gas phase, dimerization reactions are not favored thermodynamically due to a significant entropic penalty for the reaction. Therefore, consistent with experimental observations, sublimation-based purification strategies remain viable. The experimental part of the report concludes with a short discussion of the issues that still remain in developing a technologically relevant Zircaloy recycling strategy. The final section of the report is a summary of the quantum chemical investigation of the reactions of metal chlorides in the gas phase under sublimation conditions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Unsteady fluid flow in a slightly curved pipe: A comparative study of a matched asymptotic expansions solution with a single analytical solution

The present work is motivated by the fact that blood flow in the aorta and the main arteries is governed by large finite values of the Womersley number α and for such values of α there is not any analytical solution in the literature. The existing numerical solutions, although accurate, give limited information about the factors that affect the flow, whereas an analytical approach has an advantage in that it can provide physical insight to the flow mechanism. Having this in mind, we seek analytical solution to the equations of the fluid flow driven by a sinusoidal pressure gradient in a slightly curved pipe of circular cross section when the Womersley number varies from small finite to infinite values. Initially the equations of motion are expanded in terms of the curvature ratio δ and the resulting linearized equations are solved analytically in two ways. In the first, we match the solution for the main core to that for the Stokes boundary layer. This solution is valid for very large values of α. In the second, we derive a straightforward single solution valid to the entire flow region and for 8 ≤ α < ∞, a range which includes the values of α that refer to the physiological flows. Each solution contains expressions for the axial velocity, the stream function, and the wall stresses and is compared to the analogous forms presented in other studies. The two solutions give identical results to each other regarding the axial flow but differ in the secondary flow and the circumferential wall stress, due to the approximations employed in the matched asymptotic expansion process. The results on the stream function from the second solution are in agreement with analogous results from other numerical solutions. The second solution predicts that the atherosclerotic plaques may develop in any location around the cross section of the aortic wall unlike to the prescribed locations predicted by the first solution. In addition, it gives circumferential wall stresses augmented by approximately 100% with respect to the matched asymptotic expansions, a factor that may contribute jointly with other pathological factors to the faster aging of the arterial system and the possible malfunction of the aorta.

42 ENGINEERING↗

Temporal breakdown and Borel resummation in the complex Langevin method

We reexamine the Parisi-Klauder conjecture for complex e{sup i{theta}/2}{phi}{sup 4} measures with a Wick rotation angle 0{<=}{theta}/2{<=}{pi}/2 interpolating between Euclidean signature and Lorentzian signature. Our main result is that the asymptotics for short stochastic times t encapsulates information also about the equilibrium aspects. The moments evaluated with the complex measure and with the real measure defined by the stochastic Langevin equation have the same t{yields}0 asymptotic expansion which is shown to be Borel summable. The Borel transform correctly reproduces the time dependent moments of the complex measure for all t, including their t{yields}{infinity} equilibrium values. On the other hand the results of a direct numerical simulation of the Langevin moments are found to disagree from the 'correct' result for t larger than a finite t{sub c}. The breakdown time t{sub c} increases powerlike for decreasing strength of the noise's imaginary part but cannot be excluded to be finite for purely real noise. To ascertain the discrepancy we also compute the real equilibrium distribution for complex noise explicitly and verify that its moments differ from those obtained with the complex measure. - Highlights: Black-Right-Pointing-Pointer The Parisi-Klauder conjecture is reexamined for complex e{sup i{theta}/2}{phi}{sup 4} measures. Black-Right-Pointing-Pointer The time dependent moments are evaluated by temporal Borel resummation. Black-Right-Pointing-Pointer The results disagree with the Langevin simulations beyond a critical time t{sub c}. Black-Right-Pointing-Pointer t{sub c} increases with decreasing strength of the noise's imaginary part. Black-Right-Pointing-Pointer The technical reason for the breakdown is identified.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Near-infrared spectral tomography integrated with digital breast tomosynthesis: Effects of tissue scattering on optical data acquisition design

Purpose: Design optimization and phantom validation of an integrated digital breast tomosynthesis (DBT) and near-infrared spectral tomography (NIRST) system targeting improvement in sensitivity and specificity of breast cancer detection is presented. Factors affecting instrumentation design include minimization of cost, complexity, and examination time while maintaining high fidelity NIRST measurements with sufficient information to recover accurate optical property maps. Methods: Reconstructed DBT slices from eight patients with abnormal mammograms provided anatomical information for the NIRST simulations. A limited frequency domain (FD) and extensive continuous wave (CW) NIRST system was modeled. The FD components provided tissue scattering estimations used in the reconstruction of the CW data. Scattering estimates were perturbed to study the effects on hemoglobin recovery. Breast mimicking agar phantoms with inclusions were imaged using the combined DBT/NIRST system for comparison with simulation results. Results: Patient simulations derived from DBT images show successful reconstruction of both normal and malignant lesions in the breast. They also demonstrate the importance of accurately quantifying tissue scattering. Specifically, 20% errors in optical scattering resulted in 22.6% or 35.1% error in quantification of total hemoglobin concentrations, depending on whether scattering was over- or underestimated, respectively. Limited frequency-domain optical signal sampling provided two regions scattering estimates (for fat and fibroglandular tissues) that led to hemoglobin concentrations that reduced the error in the tumor region by 31% relative to when a single estimate of optical scattering was used throughout the breast volume of interest. Acquiring frequency-domain data with six wavelengths instead of three did not significantly improve the hemoglobin concentration estimates. Simulation results were confirmed through experiments in two-region breast mimicking gelatin phantoms. Conclusions: Accurate characterization of scattering is necessary for quantification of hemoglobin. Based on this study, a system design is described to optimally combine breast tomosynthesis with NIRST.

60 APPLIED LIFE SCIENCES↗

A virtual trial framework for quantifying the detectability of masses in breast tomosynthesis projection data

Purpose: Digital breast tomosynthesis (DBT) is a promising breast cancer screening tool that has already begun making inroads into clinical practice. However, there is ongoing debate over how to quantitatively evaluate and optimize these systems, because different definitions of image quality can lead to different optimal design strategies. Powerful and accurate tools are desired to extend our understanding of DBT system optimization and validate published design principles. Methods: The authors developed a virtual trial framework for task-specific DBT assessment that uses digital phantoms, open-source x-ray transport codes, and a projection-space, spatial-domain observer model for quantitative system evaluation. The authors considered evaluation of reconstruction algorithms as a separate problem and focused on the information content in the raw, unfiltered projection images. Specifically, the authors investigated the effects of scan angle and number of angular projections on detectability of a small (3 mm diameter) signal embedded in randomly-varying anatomical backgrounds. Detectability was measured by the area under the receiver-operating characteristic curve (AUC). Experiments were repeated for three test cases where the detectability-limiting factor was anatomical variability, quantum noise, or electronic noise. The authors also juxtaposed the virtual trial framework with other published studies to illustrate its advantages and disadvantages. Results: The large number of variables in a virtual DBT study make it difficult to directly compare different authors' results, so each result must be interpreted within the context of the specific virtual trial framework. The following results apply to 25% density phantoms with 5.15 cm compressed thickness and 500 {mu}m{sup 3} voxels (larger 500 {mu}m{sup 2} detector pixels were used to avoid voxel-edge artifacts): 1. For raw, unfiltered projection images in the anatomical-variability-limited regime, AUC appeared to remain constant or increase slightly with scan angle. 2. In the same regime, when the authors fixed the scan angle, AUC increased asymptotically with the number of projections. The threshold number of projections for asymptotic AUC performance depended on the scan angle. In the quantum- and electronic-noise dominant regimes, AUC behaviors as a function of scan angle and number of projections sometimes differed from the anatomy-limited regime. For example, with a fixed scan angle, AUC generally decreased with the number of projections in the electronic-noise dominant regime. These results are intended to demonstrate the capabilities of the virtual trial framework, not to be used as optimization rules for DBT. Conclusions: The authors have demonstrated a novel simulation framework and tools for evaluating DBT systems in an objective, task-specific manner. This framework facilitates further investigation of image quality tradeoffs in DBT.

60 APPLIED LIFE SCIENCES↗

Electrical switching of a p-wave magnet

Magnetic states with zero magnetization but non-relativistic spin splitting are outstanding candidates for the next generation of spintronic devices. Their electronvolt (eV)-scale spin splitting, ultrafast spin dynamics and nearly vanishing stray fields make them particularly promising for several applications. A variety of such magnetic states with non-trivial spin textures have been identified recently, including even-parity d-wave, g-wave or i-wave altermagnets and odd-parity p-wave magnets. Achieving voltage-based control of the non-uniform spin polarization of these magnetic states is of great interest for realizing energy-efficient and compact devices for information storage and processing. Spin-spiral type II multiferroics are optimal candidates for such voltage-based control, as they exhibit an inversion-symmetry-breaking magnetic order that directly induces ferroelectric polarization, allowing for symmetry-protected cross-control between spin chirality and polar order. Here we combine photocurrent measurements, first-principles calculations and group-theory analysis to provide direct evidence that the spin polarization of the spin-spiral type II multiferroic NiI 2 exhibits odd-parity character connected to the spiral chirality. The symmetry-protected coupling between chirality and polar order enables electrical control of a primarily non-relativistic spin polarization. Our findings represent an observation of p-wave magnetism in a spin-spiral type II multiferroic, which may lead to the development of voltage-based switching of non-relativistic spin polarization in compensated magnets.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Modeling X-Ray Thomson Scattering Spectra of Warm Dense Matter

Warm dense matter (WDM) is sometimes described as a \solid density plasma" where neither weakly coupled plasma theory nor condensed matter theory are applicable. More specifically, WDM has densities of ~ 0:1 - 50 times solid density, temperatures of ~ 1 - 100 eV, is partially ionized, globally neutral, has strongly correlated particles (due to interactions), and partial electron degeneracy. In nature, WDM may be found in the interiors of planets such as Jupiter or in the envelopes of white dwarf stars, and in man-made processes such as inertial confinement fusion (ICF), WDM exists as a transient state of matter. To model these physical systems, the equation of state as well as opacity, conductivity, and diffusion coefficients must be known. Experimentally WDM is probed with an experimental technique known as X-Ray Thomson Scattering (XRTS). In such an experiment X-rays photons are scattered off of WDM and the resulting photons are measured at a given solid angle. The resulting spectrum gives information on the density and temperature of WDM, but a model of WDM and the scattering experiment is often required to interpret the results. As a result, these experiments also serve as a method of validation for potential models of WDM. The work presented here focuses on the model of WDM proposed by Starrett and Saumon and a model for scattering based off of Chihara’s work. The model proposed by Starrett and Saumon is an average-atom model that is coupled with a two component plasma model, hereafter referred to as an AA+TCP model. The AA+TCP model is an ab initio model with no free parameters, only requiring knowledge of the element under investigation, the density, and the temperature. The outputs from the model completely determine all the input parameters needed for Chihara’s structure function (chemical potential, elastic feature, screening potential, bound-state wavefunctions, and free wavefunctions), which in turn determines the double differential scattering cross-section of the experiment. This approach to determining the outcome of an XRTS experiment parallels that of Johnson, et al., but expands the work in many respects. For example, the AA+TCP model completely determines all the input parameters to the structure function; hence no ad hoc approximations are necessary to close the system. Also, the free-electron feature and the bound-free feature are modified to include collisions and occupation of free electron states, respectively.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Optical and microstructural studies of erbium-doped TiO 2 thin films on silicon, SrTiO 3 , and sapphire

Rare-earth ion doped oxide thin films integrated on silicon substrates provide a route toward scalable, chip-scale platforms for quantum coherent devices. Erbium-doped TiO 2 is an attractive candidate: the Er 3+ optical transition is compatible with C-band optical fiber communications, while TiO 2 is an insulating dielectric compatible with silicon process technology. Through structural and optical studies of Er-doped TiO 2 thin films grown via molecular beam deposition on silicon, SrTiO 3 , and sapphire substrates, we have explored the impact of polycrystallinity and microstructure on the optical properties of the Er emission. Comparing polycrystalline TiO 2 (rutile)/Si with single-crystalline TiO 2 (rutile)/r-sapphire and polycrystalline TiO 2 (anatase)/Si with single-crystalline TiO 2 (anatase)/SrTiO 3 , we observe that the inhomogeneous linewidth (Γ inh ) of the most prominent peak in the Er spectrum (the Y 1 –Z 1 transition, 1520 and 1533 nm in rutile and anatase TiO 2 ) is significantly narrower in the polycrystalline case. This implies a relative insensitivity to extended structural defects and grain boundaries in such films (as opposed to, e.g., point defects). We show that the growth of an undoped, underlying TiO 2 buffer on Si can reduce Γ inh by a factor of 4–5. Expectedly, Γ inh also reduces with decreasing Er concentrations: we observe a ∼2 order of magnitude reduction from ∼1000 ppm Er to ∼10 ppm Er. Γ inh then gets limited to a residual value of ∼5 GHz that is insensitive to further reduction in the Er concentration. Based upon the above results, we argue that the optical properties in these thin films are limited by the presence of high “grown-in” point defect concentrations.

Chemical elements↗

Evaluation of the microangiographic fluoroscope (MAF) using generalized system performance metrics

Purpose: The performance of a newly developed, high resolution, microangiographic fluoroscope (MAF) (35 {mu}m pixel pitch and 300 {mu}m thick CsI phosphor) was evaluated using a generalized linear system analysis and compared with that of a standard amorphous Si thin film transistor flat panel detector (FPD) (194 {mu}m pixel pitch and 600 {mu}m thick CsI phosphor). The linear system metrics such as modulation transfer function (MTF), noise power spectrum, and detection quantum efficiency (DQE) are commonly used to gauge the intrinsic detector performance in the detector plane. However, these linear system metrics do not provide information about the image receptor performance in a real system since they do not include the effects of other parameters such as focal spot distribution, scatter radiation, and geometric unsharpness, which may compromise detector performance characteristics. Use of generalized linear system metrics [generalized modulation transfer function (GMTF), generalized normalized noise power spectrum (GNNPS), and generalized detection quantum efficiency (GDQE)] that include these effects gives a more meaningful, complete, and appropriate evaluation of detector performance as part of the imaging system. Methods: A uniform head equivalent phantom was used to simulate realistic clinical parameters and x-ray spectra. The detector MTFs were measured using the slanted edge method and the focal spot MTFs were measured using a pinhole assembly. The scatter MTF was simulated and the scatter fraction was measured for a head-equivalent phantom. The generalized system metrics were calculated for different combinations of three choices of focal spots and three different magnifications with two different air-gaps. The performance of the MAF was also illustrated using stent images obtained with different focal spots under similar conditions. Results: Results for the generalized metrics provide a quantitative description of the performance of the imaging system for both detectors. This generalized analysis demonstrated that both detectors have similar imaging capabilities at lower spatial frequencies, but that the MAF has superior performance over the FPD at higher frequencies even when considering focal spot blurring and scatter. Conclusions: This generalized performance analysis demonstrates the significance of focal spot size, magnification, and scatter on the system performance metrics (GMTF, GNNPS, and GDQE). Although the ideal detector performance characteristics of the MAF are not fully realized due to these other system factors, it still retains an advantage in DQE at high spatial frequencies over the FPD. Similar studies based on the generalized linear system metrics can serve as an efficient tool to evaluate total system capabilities under different realistic conditions to enable optimal design for specific imaging tasks.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

A physics-informed and hierarchically regularized data-driven model for predicting fluid flow through porous media

This paper presents a new deep learning data-driven model for predicting structure dependent pore-fluid velocity fields in rock. The model is based on a Convolutional Auto-Encoder (CAE) artificial neural network capable of learning from image data generated by direct numerical simulations of fluid flow through pore-structures, such as by Lattice Boltzmann or molecular dynamics methods. The main novelty of the model in comparison to previous CAE-based data-driven approaches consists of three parts. The first is a methodology for decomposing the full-domain of the porous media into sub-regions, or “sub-domains”, in order to reduce the overall size of the CAE, batch process the sub-domains in parallel, and enable the CAE to learn local and generalizable nonlinear mappings of pore-fluid velocities. The second consists of embedding the finite difference solutions of the incompressible Navier-Stokes and continuity equations into convolutional layers prior to the CAE in order to provide the CAE with knowledge of fluid dynamics physics (PhyFlow). The third main novelty is that the training of the CAE is regularized with a hierarchical loss function that encourages the learning of fluid flow patterns (in a way similar to ranked modes in principal component analysis), ranking from most to least important. This is shown to increase the stability in learning, reduce over-fitting, and promote interpretability of the CAE neural network layers (HierCAE). The comprehensive new data-driven model, which we call the PhyFlow-HierCAE model, is shown to exhibit improved accuracy and generalizability of flow field predictions over conventional CAE models, attributable to the embedded physical knowledge and the hierarchical regularization, as well as realize orders of magnitude speed-ups in computation times as a surrogate for the direct numerical simulations. Examples of training and forward predictions on unseen pore-structures are provided and evaluated for data from Lattice Boltzmann and molecular dynamics simulations of pore-fluid flow. The model is shown to be a fast and accurate emulator (or “surrogate”) for predicting effective permeability of unseen pore-structures based on learning from relatively small direct numerical simulation datasets.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗