Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Density matrix methods”

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,261 records · Page 70
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Matrix product operators, matrix product states, and ab initio density matrix renormalization group algorithms

Current descriptions of the ab initio density matrix renormalization group (DMRG) algorithm use two superficially different languages: an older language of the renormalization group and renormalized operators, and a more recent language of matrix product states and matrix product operators. The same algorithm can appear dramatically different when written in the two different vocabularies. Here, we carefully describe the translation between the two languages in several contexts. First, we describe how to efficiently implement the ab initio DMRG sweep using a matrix product operator based code, and the equivalence to the original renormalized operator implementation. Next we describe how to implement the general matrix product operator/matrix product state algebra within a pure renormalized operator-based DMRG code. Finally, we discuss two improvements of the ab initio DMRG sweep algorithm motivated by matrix product operator language: Hamiltonian compression, and a sum over operators representation that allows for perfect computational parallelism. The connections and correspondences described here serve to link the future developments with the past and are important in the efficient implementation of continuing advances in ab initio DMRG and related algorithms.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2015

A general framework of noise suppression in material decomposition for dual-energy CT

Purpose: As a general problem of dual-energy CT (DECT), noise amplification in material decomposition severely reduces the signal-to-noise ratio on the decomposed images compared to that on the original CT images. In this work, the authors propose a general framework of noise suppression in material decomposition for DECT. The method is based on an iterative algorithm recently developed in their group for image-domain decomposition of DECT, with an extension to include nonlinear decomposition models. The generalized framework of iterative DECT decomposition enables beam-hardening correction with simultaneous noise suppression, which improves the clinical benefits of DECT. Methods: The authors propose to suppress noise on the decomposed images of DECT using convex optimization, which is formulated in the form of least-squares estimation with smoothness regularization. Based on the design principles of a best linear unbiased estimator, the authors include the inverse of the estimated variance–covariance matrix of the decomposed images as the penalty weight in the least-squares term. Analytical formulas are derived to compute the variance–covariance matrix for decomposed images with general-form numerical or analytical decomposition. As a demonstration, the authors implement the proposed algorithm on phantom data using an empirical polynomial function of decomposition measured on a calibration scan. The polynomial coefficients are determined from the projection data acquired on a wedge phantom, and the signal decomposition is performed in the projection domain. Results: On the Catphan{sup ®}600 phantom, the proposed noise suppression method reduces the average noise standard deviation of basis material images by one to two orders of magnitude, with a superior performance on spatial resolution as shown in comparisons of line-pair images and modulation transfer function measurements. On the synthesized monoenergetic CT images, the noise standard deviation is reduced by a factor of 2–3. By using nonlinear decomposition on projections, the authors’ method effectively suppresses the streaking artifacts of beam hardening and obtains more uniform images than their previous approach based on a linear model. Similar performance of noise suppression is observed in the results of an anthropomorphic head phantom and a pediatric chest phantom generated by the proposed method. With beam-hardening correction enabled by their approach, the image spatial nonuniformity on the head phantom is reduced from around 10% on the original CT images to 4.9% on the synthesized monoenergetic CT image. On the pediatric chest phantom, their method suppresses image noise standard deviation by a factor of around 7.5, and compared with linear decomposition, it reduces the estimation error of electron densities from 33.3% to 8.6%. Conclusions: The authors propose a general framework of noise suppression in material decomposition for DECT. Phantom studies have shown the proposed method improves the image uniformity and the accuracy of electron density measurements by effective beam-hardening correction and reduces noise level without noticeable resolution loss.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

Examination of Hybrid Metal Coatings for Mitigation of Fission Product Release and Corrosion Protection of LWR SiC/SiC

There is a need to increase the safety margins of current and future light water reactors (LWRs) due to the unfortunate events at Fukushima Daiichi Nuclear Plant. Safety is crucial to restore public confidence in nuclear energy, acknowledged as an economical, high-­density energy solution to climate change. The development of accident-­tolerant fuel (ATF) concepts is crucial to this endeavor. The objective of ATF is to delay the consequences of accident progression, being inset in high temperature steam and maintaining high thermomechanical strength for radionuclide retention. The use of advanced SiCf-­SiC composite as a substitute for zircaloy-­based cladding is being considered. However, at normal operations, SiC is vulnerable to the reactor coolant and may corrode at an unacceptable rate. As a ceramic-­matrix composite material, it is likely to undergo microcracking operation, which may compromise the ability to contain gaseous fission products. A proposed solution to both issues is the application of mitigation coatings for use in normal operations. At Oak Ridge National Laboratory (ORNL), three coating technologies have been investigated with industry collaborators and vendors. These are electrochemical deposition, cathodic arc physical vapor deposition (PVD hereafter) and vacuum plasma spray (VPS). The objective of this document is to summarize these processing technologies, the resultant as-­processed microstructures and properties of the coatings. In all processes, substrate constraint resulted in substantial tensile stresses within the coating layer. Each technology must mitigate this tensile stress. Electrochemical coatings use chromium as the coolant facing material, and are deposited on a nickel or carbon “bond coat”. This is economical but suffers microcracking in the chromium layer. PVD-­based coatings use chromium and titanium in both metallic form and nitrides, and can be deposited defense-­in-­depth as multilayers. This vapor method eliminated tensile stress during processing and coatings were up to ~30 μm thick without microcracking. VPS produced coatings based on Zircaloy-2, which has a proven reactor-­compatibility. The tensile stresses appearred to be partially mitigated by annealing. Analysis showed that VPS coatings required further optimizations to prevent adverse reactions with the substrate and need a minimum thickness of ~50 μm. In addition, development of coatings are constrained by neutronic depletion analysis, which clearly indicated enrichment as an issue if the coating is too thick. Based on the present work, the cathodic arc PVD technology was considered ready for the extensive testing and evaluation on cladding materials due to their ability to mitigate the excessive tensile stresses and the reasonable coating quality achieved. The VPS Zircaloy-­2 coating technology required additional development toward mitigation of issues related to the substrate reaction and porosity. In the future, PVD and VPS will have be scaled upon successful development and demonstration. Electrochemical coatings, which are proven scalability, currently require development to mitigate issues related to the tensile stress after deposition.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2013

Reactions between atomic chlorine and pyridine in solid para-hydrogen: Infrared spectrum of the 1-chloropyridinyl (C{sub 5}H{sub 5}N-Cl) radical

With infrared absorption spectra we investigated the reaction between Cl atom and pyridine (C{sub 5}H{sub 5}N) in a para-hydrogen (p-H{sub 2}) matrix. Pyridine and Cl{sub 2} were co-deposited with p-H{sub 2} at 3.2 K; a planar C{sub 5}H{sub 5}N-Cl{sub 2} complex was identified from the observed infrared spectrum of the Cl{sub 2}/C{sub 5}H{sub 5}N/p-H{sub 2} matrix. Upon irradiation at 365 nm to generate Cl atom in situ and annealing at 5.1 K for 3 min to induce secondary reaction, the 1-chloropyridinyl radical (C{sub 5}H{sub 5}N-Cl) was identified as the major product of the reaction Cl + C{sub 5}H{sub 5}N in solid p-H{sub 2}; absorption lines at 3075.9, 1449.7, 1200.6, 1148.8, 1069.3, 1017.4, 742.9, and 688.7 cm{sup -1} were observed. The assignments are based on comparison of observed vibrational wavenumbers and relative IR intensities with those predicted using the B3PW91/6-311++G(2d, 2p) method. The observation of the preferential addition of Cl to the N-site of pyridine to form C{sub 5}H{sub 5}N-Cl radical but not 2-, 3-, or 4-chloropyridine (ClC{sub 5}H{sub 5}N) radicals is consistent with the reported theoretical prediction that formation of the former proceeds via a barrierless path.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRSDate not supplied

Use of Design of Experiments and Rule-Based Inference in Determining Neural Network Architectures for Loss of Control Detection

In this work, we describe methods for selecting the neural network architectures and input spaces to implement belief state inference on generic commercial transport aircraft. First, we highlight a case study on the planning, execution, and analysis of a set of experiments to determine the configurations of a conditional variational autoencoder (CVAE). We present a structured method that can be used in a number of aerospace applications, to optimize the structure and training parameters of the CVAE for belief state inference, using Design of Experiments (DOE) statistical methodologies. The motivation for this specific DOE was to identify the appropriate hyperparameters for measuring the CVAE reconstruction probability and latent space, such that the measurements can be used to infer qualitative state changes for the aircraft. We demonstrate that this process yields information about a trained neural network’s utility for this specific application, along with a quantifiable range of certainty. We execute 84 experiments using loss-of-control flight maneuver data from a NASA T-2 aircraft, demonstrating that this empirical process allows us to construct cheap and simple models with specific attributes amenable to belief state inference in aerospace applications. While theoretically, we could create a single CVAE with an input space the size of all measurable flight variables and environmental dynamics, it becomes intractable to use such a neural network in an in-situ intelligent multi-agent system. Using the recommendations from our case study, we introduce a technical approach for feasibly describing the belief space by (1) identifying significant statistical relationships among flight variables using rule induction, (2) using a set of rules that cover all features to define the input space of multiple CVAEs, and (3) forming a belief space based on the joint probability density of their collective latent spaces. This results in a series of relatively small matrix multiplications that can be performed in real time, as opposed to large matrix computations in a single CVAE. We demonstrate the application of this approach on the T-2 flight loss-of control experiments, using the architecture and hyperparameter recommendations from the case study. We compare the utilities of an individual CVAE trained on all flight variables and multiple CVAEs defined on subsets of flight variables for detecting qualitative changes in flight. We demonstrate that the use of multiple CVAEs with smaller input spaces permits the CVAE to capture more granular relationships in the latent space, permitting better state space characterization and loss-of-control detection.

Design of experiments↗
DOE OSTI2018

Leaching Behavior of Rare Earth Elements in Fort Union Lignite Coals of North America

Rare earth elements are crucial materials in an incredible array of consumer goods, energy system components and military defense applications. However, the global production and entire value chain for rare earth elements is dominated by China, with the U.S. currently 100% import reliant for these critical materials. Traditional mineral ores including those previously mined in the U.S., however, have several challenges. Chief among these is that the content of the most critical and valuable of the rare earths are deficient, making mining uneconomical. Further, the supply of these most critical rare earths is nearly 100% produced in China from a single resource that is only projected to last another 10 to 20 years. The U.S. currently considers the rare earths market an issue of national security. It is imperative that alternative domestic sources of rare earths be identified and methods developed to produce them. Recently, coal and coal byproducts have been identified as one of these promising alternative resources. This paper details the results of a study on characterization of North Dakota lignite and lignite-related feedstocks as an assessment of their feasibility for rare earth element recovery. The abundance, distribution and modes of occurrence of the rare earth elements in the samples collected were determined in this initial study to inform the selection of appropriate extraction and concentration methods to recover the rare earth elements. Materials investigated include the lignite coals, clay-rich sediments associated with the coal seams, and materials associated with a lignite beneficiation system and power plant. The results show that high rare earth element levels exist both in lignite coals and associated sediments. The form of the rare earth elements in the clay materials is primarily as ultra-fine mineral grains. In the lignite coals, approximately 80-95% of the rare earths content is organically associated, primarily as coordination complexes.

Laudal, Daniel A.↗
DOE OSTI2015

Transport properties of alkali metal doped fullerides

We have studied the intercage interactions between the adjacent C{sub 60} cages and expansion of lattice due to the intercalation of alkali atoms based on the spring model to estimate phonon frequencies from the dynamical matrix for the intermolecular alkali-C{sub 60} phonons. We considered a two-peak model for the phonon density of states to investigate the nature of electron pairing mechanism for superconducting state in fullerides. Coulomb repulsive parameter and the electron phonon coupling strength are obtained within the random phase approximation. Transition temperature, T{sub c}, is obtained in a situation when the free electrons in lowest molecular orbital are coupled with alkali-C{sub 60} phonons as 5 K, which is much lower as compared to reported T{sub c} (20 K). The superconducting pairing is mainly driven by the high frequency intramolecular phonons and their effects enhance it to 22 K. The importance of the present study, the pressure effect and normal state transport properties are calculated within the same model leading superconductivity.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2012

Oxygen diffusion in UO 2+x and (U,Pu)O 2±x

In the first part of this report we revisit an earlier study of oxygen diffusion in UO 2+x , in which we used density functional theory (DFT) calculations to parameterize a kinetic Monte Carlo (kMC) model. The results from these earlier kMC simulations are reproduced in Fig. 1 and they indicate fairly good agreement with available experiments. This work was later expanded to include a larger temperature range. However, since the publication of this study there have been a number of advancements in DFT methodology for UO 2 and UO 2+x providing increased accuracy. We have also gained better understanding of the oxygen clustering phenomena occurring in UO 2+x . For these two reasons, the DFT calculations of the migration barriers of single oxygen interstitials and di-interstitial clusters have been repeated using the LDA+U and GGA+U methodologies. The earlier study used regular GGA and, even though this method captures similar trends as the more advanced LDA+ U and GGA+ U techniques, it does not fulfill the quantitative requirements set by some applications. Additionally, we have identified a mechanism for the most stable quad-interstitial clusters to migrate and here we calculate the corresponding barriers within both the LDA+ U and GGA+ U methodologies. The new LDA+ U and GGA+ U data sets are analyzed in terms of available experiments. In the second part of this report we present initial results for the impact of Pu on oxygen diffusion in UO 2 . The first step in understanding this process is to calculate the binding energies of oxygen vacancies and interstitials to a Pu ion in the UO 2 matrix. Possible diffusion mechanisms are discussed for (U,Pu)O 2-x , (U,Pu)O 2 and (U,Pu)O 2+x .

36 MATERIALS SCIENCE↗
DOE OSTI2018

Improved Characterization and Modeling of Tight Oil Formations for CO 2 Enhanced Oil Recovery Potential and Storage Capacity Estimation

Tight oil formations such as those in the Bakken petroleum system are known to hold hundreds of billions of barrels of oil in place; however, the primary recovery factor for these plays is typically less than 10%. Tight oil formations, including the Bakken Formation, therefore, may be attractive candidates for enhanced oil recovery (EOR) using CO 2 . Multiphase fluid behavior and flow in fluid-rich shales can vary substantially depending on the size of pore throats, and properties such as fluid viscosity and density are much different in nanoscale pores than in macroscale pores. Thus it is critical to understand the nature and distribution of nano-, micro-, and macroscale pores and fracture networks. To address these issues, the Energy & Environmental Research Center (EERC) has been conducting a research program entitled “Improved Characterization and Modeling of Tight Oil Formations for CO 2 Enhanced Oil Recovery Potential and Storage Capacity Estimation.” The objectives of the project are 1) the use of advanced characterization methods to better understand and quantify the petrophysical and geomechanical factors that control CO 2 and oil mobility within tight oil formation samples, 2) the determination of CO 2 permeation and oil extraction rates in tight reservoir rocks and organic-rich shales of the Bakken, and 3) the integration of the laboratory-based CO 2 permeation and oil extraction data and the characterization data into geologic models and dynamic simulations to develop predictions of CO 2 storage resource and EOR in the Bakken tight oil formation. A combination of standard and advanced petrophysical characterization techniques were applied to characterize samples of Bakken Formation tight reservoir rock and shales from multiple wells. Techniques included advanced computer tomography (CT) imaging, scanning electron microscopy (SEM) techniques, whole-core and micro x-ray CT imaging, field emission (FE) SEM, and focused ion beam (FIB) SEM. Selected samples were also analyzed for geomechanical properties. X-ray CT imaging yielded information on the occurrence of fractures, bedding planes, fossils, and bioturbation in core, as well as data on bulk density and photoelectric factor logs, which were used to interpret porosity, organic content, and mineralogy. FESEM was used for characterization of nano- and microscale features, including nanoscale pore visualization and micropore and pore throat mineralogy. FIBSEM yielded micro- to nanoscale visualization of fracture networks, porosity and pore-size distribution, connected versus isolated porosity, and distribution of organics. Results from the characterization activities provide insight on nanoscale fracture properties, pore throat mineralogy and connectivity, rock matrix characteristics, mineralogy, and organic content. Laboratory experiments demonstrated that CO 2 can permeate the tight matrix of Bakken shale and nonshale reservoir samples and mobilize oil from those samples. Geologic models were created at scales ranging from the core plug to the reservoir, and dynamic simulations were conducted. The data from the characterization and laboratory-based activities were integrated into modeling research activities to determine the fundamental mechanisms controlling fluid transport in the Bakken, which support EOR scheme design and estimation of CO 2 storage potential in tight oil formations. Simulation results suggest a CO 2 storage resource estimate range of 169 million to 1.5 billion tonnes for the Bakken in North Dakota, possibly resulting in 1.8 billion to 16 billion barrels of incremental oil.

04 OIL SHALES AND TAR SANDS↗
DOE OSTI2021

Fast and robust all-electron density functional theory calculations in solids using orthogonalized enriched finite elements

Here, we present a computationally efficient approach to perform systematically convergent real-space all-electron Kohn-Sham density functional theory calculations for solids using an enriched finite element (FE) basis. The enriched FE basis is constructed by augmenting the classical FE basis with atom-centered numerical basis functions, comprising of atomic solutions to the Kohn-Sham problem. Notably, to improve the conditioning, we orthogonalize the enrichment functions with respect to the classical FE basis, without sacrificing the locality of the resultant basis. In addition to improved conditioning, this orthogonalization procedure also renders the overlap matrix block diagonal, greatly simplifying its inversion. Subsequently, we use a Chebyshev polynomial based filtering technique to efficiently compute the occupied eigenspace in each self-consistent field iteration. We demonstrate the accuracy and efficiency of the proposed approach on periodic unit cells and supercells. The benchmark studies show a staggering 130× speedup of the orthogonalized enriched FE basis over the classical FE basis. We also present a comparison of the orthogonalized enriched FE basis with the linearized augmented plane-wave + local orbitals basis, both in terms of accuracy and efficiency. Notably, we demonstrate that the orthogonalized enriched FE basis can handle large system sizes of ~10 000 electrons. Finally, we observe good parallel scalability of our implementation with 92% efficiency at 22× speedup for a system with 620 electrons

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2017

A Robust and Conductive Black Tin Oxide Nanostructure Makes Efficient Lithium-Ion Batteries Possible

SnO 2 -based lithium-ion batteries have low cost and high energy density, but their capacity fades rapidly during lithiation/delithiation due to phase aggregation and cracking. These problems can be mitigated by using highly conducting black SnO 2-x , which homogenizes the redox reactions and stabilizes fine, fracture-resistant Sn precipitates in the Li 2 O matrix. Such fine Sn precipitates and their ample contact with Li 2 O proliferate the reversible Sn → Li x Sn → Sn → SnO 2 /SnO 2-x cycle during charging/discharging. SnO 2-x electrode has a reversible capacity of 1340 mAh g -1 and retains 590 mAh g -1 after 100 cycles. The addition of highly conductive, well-dispersed reduced graphene oxide further stabilizes and improves its performance, allowing 950 mAh g -1 remaining after 100 cycles at 0.2 A g -1 with 700 mAh g -1 at 2.0 A g -1 . Conductivity-directed microstructure development may offer a new approach to form advanced electrodes.

25 ENERGY STORAGE↗
DOE OSTI2016

As cast precipitation microstructures in twin-roller melt-spun Cu{sub 90}Co{sub 10} alloys

As cast Cu{sub 90}Co{sub 10} ribbons rapidly solidified by twin-roller melt spinning, exhibit special microstructure features. This processing method provides scenarios where a different phase selection takes place; coherent Co precipitates form directly from solidification, with neither a spinodal-like composition oscillation nor a discontinuously precipitated laminar phase. Samples are processed at tangential wheel speeds of 10 m/s (V10), 15 m/s (V15), 20 m/s (V20) and 30 m/s (V30). Microstructures resulting from this single step process are characterized and the hysteresis properties and the magnetoresistance effects evaluated. Samples V30 have a quite uniform density of coherent precipitates, with a narrow size distribution around 4 nm. On contrary, non-uniform precipitate distributions are found in samples cooled at lower rates; zones with a high density of coherent Co-rich precipitates are found forming colonies. These colonies are consistent with the extended compositional fluctuations occurring during very early stages in the cooling process. Samples may exhibit wide (V10) and even bimodal (V15) size distributions. Only samples V30 behave close to the ideal superparamagnetism. Samples V20 present relatively large coercivity and relative remanence and behave as an interacting superparamagnet, while the hysteresis loops of ribbons cooled at lower rates exhibit a ferromagnetic contribution in addition to the superparamagnetic-like one. This ferromagnetic component arises from blocked precipitates, larger than the upper bound size for superparamagnetic behavior at 300 K (12 nm). Room temperature magnetoresistance values associated to granular scattering units decrease as the mean precipitate size increases, but they remain below 2%, which is lower than that measured in samples annealed after rapid solidification, indicating that in this latter case contributions from the spinodally segregated matrix take place in addition to that of Co granules. - Highlights: •Microstructures of twin-rolled Cu{sub 90}Cu{sub 10} alloys are for the first time characterized. •Coherent Co-rich precipitates distribute forming dense colonies inside the grains. •No coexisting multilayer-like structures (spinodal or laminar phases) are detected. •Magnetic hysteresis arises from inter-particle interaction in these dense colonies. •Magnetoresistance level in these purely granular systems is low (1.6% at 300 K).

36 MATERIALS SCIENCE↗
NASA NTRS1995

Hot-wire calibration in subsonic/transonic flow regimes

A different approach for calibrating hot-wires, which simplifies the calibration procedure and reduces the tunnel run-time by an order of magnitude was sought. In general, it is accepted that the directly measurable quantities in any flow are velocity, density, and total temperature. Very few facilities have the capability of varying the total temperature over an adequate range. However, if the overheat temperature parameter, a(sub w), is used to calibrate the hot-wire then the directly measurable quantity, voltage, will be a function of the flow variables and the overheat parameter i.e., E = f(u,p,a(sub w), T(sub w)) where a(sub w) will contain the needed total temperature information. In this report, various methods of evaluating sensitivities with different dependent and independent variables to calibrate a 3-Wire hot-wire probe using a constant temperature anemometer (CTA) in subsonic/transonic flow regimes is presented. The advantage of using a(sub w) as the independent variable instead of total temperature, t(sub o), or overheat temperature parameter, tau, is that while running a calibration test it is not necessary to know the recovery factor, the coefficients in a wire resistance to temperature relationship for a given probe. It was deduced that the method employing the relationship E = f (u,p,a(sub w)) should result in the most accurate calibration of hot wire probes. Any other method would require additional measurements. Also this method will allow calibration and determination of accurate temperature fluctuation information even in atmospheric wind tunnels where there is no ability to obtain any temperature sensitivity information at present. This technique greatly simplifies the calibration process for hot-wires, provides the required calibration information needed in obtaining temperature fluctuations, and reduces both the tunnel run-time and the test matrix required to calibrate hotwires. Some of the results using the above techniques are presented in an appendix.

Nagabushana, K. A.↗
DOE OSTI2026

Influence of Markovianity and self-consistency on time-resolved spectral functions of driven quantum systems

We present a systematic comparison of the real-time Dyson expansion (RTDE) with established nonequilibrium Green's function (GF) approaches for simulating driven, interacting quantum systems. Focusing on density matrix dynamics, time-off-diagonal GFs, and time-resolved photoemission spectra, we benchmark RTDE against fully self-consistent Kadanoff-Baym equation (KBE) calculations, the generalized Kadanoff-Baym ansatz, and exact diagonalization for small systems using second-order many-body perturbation theory. Using a driven two-band Hubbard model, we show that mean-field single-particle density matrix trajectories provide a reliable baseline for RTDE across a broad range of interaction strengths and excited-carrier populations. Further, RTDE accurately captures correlation effects in the GFs, including long-lived oscillations and revivals that are strongly suppressed by the overdamping inherent to self-consistent KBE schemes. As a consequence, RTDE resolves rich nonequilibrium spectral structure in time-resolved photoemission, such as interaction- and population-dependent quasiparticle splittings and band gap renormalization, which are largely washed out in self-consistent approaches yet are present in exact solutions. Furthermore, our results demonstrate that RTDE bridges the gap between mean-field propagation and full two-time KBE simulations, retaining favorable linear scaling while capturing essential dynamical correlations relevant for ultrafast spectroscopy.

Electronic structure↗
DOE OSTI2022

The Effect of Geometry, Spin, and Orbital Optimization in Achieving Accurate, Correlated Results for Iron–Sulfur Cubanes

Iron–sulfur clusters comprise an important functional motif in the catalytic centers of biological systems, capable of enabling important chemical transformations at ambient conditions. This remarkable capability derives from a notoriously complex electronic structure that is characterized by a high density of states that is sensitive to geometric changes. The spectral sensitivity to subtle geometric changes has received little attention from correlated, large active space calculations, owing partly to the exceptional computational complexity for treating these large and correlated systems accurately. Here, to provide insight into this aspect, we report the first Complete Active Space Self Consistent Field (CASSCF) calculations for different geometries of the [Fe(II/III) 4 S 4 (SMe) 4 ] -2 clusters using two complementary, correlated solvers: spin-pure Adaptive Sampling Configuration Interaction (ASCI) and Density Matrix Renormalization Group (DMRG). We find that the previously established picture of a double-exchange driven magnetic structure, with minute energy gaps (<1 mHa) between consecutive spin states, has a weak dependence on the underlying geometry. However, the spin gap between the singlet and the spin state 2S + 1 = 19, corresponding to a maximal number of Fe-d electrons being unpaired and of parallel spin, is strongly geometry dependent, changing by a factor of 3 upon slight deformations that are still within biologically relevant parameters. The CASSCF orbital optimization procedure, using active spaces as large as 86 electrons in 52 orbitals, was found to reduce this gap compared to typical mean-field orbital approaches. Our results show the need for performing large active space calculations to unveil the challenging electronic structure of these complex catalytic centers and should serve as accurate starting points for fully correlated treatments upon inclusion of dynamical correlation outside the active space.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2016

SU-F-R-30: Interscanner Variability of Radiomics Features in Computed Tomography (CT) Using a Standard ACR Phantom

Purpose: A simple approach to investigate Interscanner variability of Radiomics features in computed tomography (CT) using a standard ACR phantom. Methods: The standard ACR phantom was scanned on CT scanners from three different manufacturers. Scanning parameters of 120 KVp, 200 mA were used while slice thickness of 3.0 mm on two scanners and 3.27 mm on third scanner was used. Three spherical regions of interest (ROI) from water, medium density and high density inserts were contoured. Ninety four Radiomics features were extracted using an in-house program. These features include shape (11), intensity (22), GLCM (26), GLZSM (11), RLM (11), and NGTDM (5) and 8 fractal dimensions features. To evaluate the Interscanner variability across three scanners, a coefficient of variation (COV) is calculated for each feature group. Each group is further classified according to the COV- by calculating the percentage of features in each of the following categories: COV less than 2%, between 2 and 10% and greater than 10%. Results: For all feature groups, similar trend was observed for three different inserts. Shape features were the most robust for all scanners as expected. 70% of the shape features had COV <2%. For intensity feature group, 2% COV varied from 9 to 32% for three scanners. All features in four groups GLCM, GLZSM, RLM and NGTDM were found to have Interscanner variability ≥2%. The fractal dimensions dependence for medium and high density inserts were similar while it was different for water inserts. Conclusion: We concluded that even for similar scanning conditions, Interscanner variability across different scanners was significant. The texture features based on GLCM, GLZSM, RLM and NGTDM are highly scanner dependent. Since the inserts of the ACR Phantom are not heterogeneous in HU values suggests that matrix based 2nd order features are highly affected by variation in noise. Research partly funded by NIH/NCI R01CA190105-01.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2021

From hard spheres to hard-core spins

A system of hard spheres exhibits physics that is controlled only by their density. This comes about because the interaction energy is either infinite or zero, so all allowed configurations have exactly the same energy. The low-density phase is liquid, while the high-density phase is crystalline, an example of “order by disorder” as it is driven purely by entropic considerations. Here we study a family of hard spin models, which we call hard-core spin models, where we replace the translational degrees of freedom of hard spheres with the orientational degrees of freedom of lattice spins. Their hard-core interaction serves analogously to divide configurations of the many spin system into allowed and disallowed sectors. We present detailed results on the square lattice in d = 2 for a set of models with $\mathbb{Z}_n$ symmetry, which generalize Potts models, and their U(1) limits, for ferromagnetic and antiferromagnetic senses of the interaction, which we refer to as exclusion and inclusion models. As the exclusion and inclusion angles are varied, we find a Kosterlitz-Thouless phase transition between a disordered phase and an ordered phase with quasi-long-ranged order, which is the form order by disorder takes in these systems. These results follow from a set of height representations, an ergodic cluster algorithm, and transfer matrix calculations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

A bi-layer barrier design for 122-type iron-based superconducting wires and tapes

Iron-based superconducting wires and tapes hold great promise for high-field magnet applications. A promising design for 122-type wires and tapes based on the powder-in-tube method is using silver and copper double-layer sheaths. For this design a heat treatment temperature below ~ 779 °C is required to prevent Ag-Cu liquid formation. However, this may be below the optimal heat treatment temperature for the critical current density, and still cannot prevent Ag-Cu interdiffusion occurring in the solid state. In this work we propose adding a niobium or tantalum or vanadium (or their alloys) barrier layer between the Ag and Cu to solve the Ag-Cu interdiffusion issue, given that the group-VB metals (vanadium, niobium, tantalum) are relatively inert to both Ag and Cu. To investigate the effectiveness of this design, BaFe 1.84 Co 0.16 As 2 wires and tapes with Ag/Cu and Ag/Ta/Cu sheaths, as well as Ba 0.6 K 0.4 Fe 2 As 2 wires and tapes with Ag/Cu and Ag/Nb/Cu sheaths, were fabricated. It was found that both the Ta and Nb layers kept integral after wire drawing, but after a large flat-rolling reduction the Ta layer broke while the Nb layer kept integral. In the tapes with Ag/Cu sheaths (without the Ta or Nb layer) Cu diffused through the Ag layer and into the powder cores during 740 °C heat treatment, while in the tapes with Ag/Nb/Cu sheaths the Nb layer effectively blocked Ag-Cu interdiffusion even at 900 °C. Further, this work demonstrates that Ta is a suitable barrier material for 122-type wires, while Nb is suitable for both wires and tapes. In this design using Ag/Nb (or Ta)/Cu sheaths, we can regard the outer Cu as the conductor matrix while the Ag and Nb (or Ta) serve as two layers of barriers that suppress reactions between the components. Thus, we call this design a “bi-layer barrier” design for 122-type wires and tapes.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗