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 685 records · Page 38

Valley Stoner instability of the composite Fermi sea

We study two-component electrons in the lowest Landau level at total filling factor ν T = 1/2 with anisotropic mass tensors and principal axes rotated by π/2 as realized in aluminum arsenide (AlAs) quantum wells. Combining exact diagonalization and the density matrix renormalization group we demonstrate that the system undergoes a quantum phase transition from a gapless state in which both flavors are equally populated to another gapless state in which all the electrons spontaneously polarize into a single flavor beyond a critical mass anisotropy of m x /m y ~7. We propose that this phase transition is a form of itinerant Stoner transition between a two-component and a single-component composite Fermi sea states and describe a set of trial wave functions which successfully capture the quantum numbers and shell filling effects in finite size systems as well as providing a physical picture for the energetics of these states. Our estimates indicate that the composite Fermi sea of AlAs is the analog of an itinerant Stoner magnet with a finite spontaneous valley polarization. Here, we pinpoint experimental evidence indicating the presence of Stoner magnetism in the Jain states surrounding ν = 1/2.

2-dimensional systems↗

Annual report for DOE VTO

Carbon fiber (CF)/polymer composites are a transformative class of high-performance, lightweight material, where high aspect-ratio CFs reinforce a polymer matrix and exceed the strength of steel alloys at a fraction of the density. Despite the advantages of such a class of material, the broader implementation of CF composites in a range of automotive, aerospace, and energy applications is hindered by limitations of current manufacturing methods. These current techniques (e.g., hand lay-up, wet filament winding) are costly and impose severe limitations on fiber placement, orientation, and angle, and thus a composite’s ultimate properties. Today’s CF composites are expensive to manufacture, limited in form factor, and utilize costly and sub-optimal continuous filament CF. Advanced additive manufacturing (AM) processes, combined with computational design optimization and new approaches to resin development, offer alternative design and manufacturing paradigms that have the realistic potential to lift these constraints. Such integrated AM approaches could thus help to realize the full potential of CF composite materials. One relevant application of CF composite materials where manufacturing constraints limit the cost-benefit ratio is in the manufacture of high-performance composite pressure vessels for onboard compressed natural gas (CNG) storage. Current CNG storage vessels (Types 3–5) are made from load-bearing filament-wound carbon-fiber composite and are ~3.5 times as expensive as an all-metallic Type-1 vessel. This cost is invariably tied to the complex and labor-intensive nature of conventional filament winding processes and the large volumes of expensive high tensile-strength CF tow feedstock required in manufacture. Our proposed approach to CNG storage vessel manufacture is based on a combination of AM technologies for CF composite printing and design optimization tools that were pioneered at Lawrence Livermore National Laboratory (LLNL), with advances in resin/composite formulation enabled by chemical and nano-material modification. Through the successful development of this technology, LLNL seeks to demonstrate the capability for advanced CNG storage vessel manufacture at reduced cost with no reduction in performance versus the most advanced, extant Type-5 designs.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Ab Initio Many Body Quantum Embedding and Local Correlation in Crystalline Materials using Interpolative Separable Density Fitting

We present an efficient implementation of ab initio many-body quantum embedding and local correlation methods for infinite periodic systems through translational symmetry adapted interpolative separable density fitting, an approach which reduces the scaling of the calculations to only linear with the number of k-points. Employing this methodology, we compute correlated ground-state coupled cluster energies within density matrix embedding and local natural orbital correlation frameworks for both weakly and strongly correlated solids, using up to 1000 k-points. By extrapolating the local correlation domains and k-point sampling we further obtain estimates of the full coupled cluster with singles, doubles, and perturbative triples ground-state energies in the thermodynamic limit.

Chemical Physics (physics.chem-ph)↗

Considerations to achieve directionality for gamma ray lasers

A method of alignment of nuclei for a gamma ray laser and a means of achieving preferential emission of radiation along the crystal axis are studied. Atomic alignment was achieved by materials researchers who made composite structures composed of needle-like single crystals all with the same orientation and all pointing in the same direction contained in a matrix of cobalt or nickel. The proposed method of preferential emission of radiation along the aligned needles is to have a symmetric field gradient at the nucleus and a sequence of excited levels of spin and parity 2(+) and 0(+). The proposed scheme reduces the density of excited states required for lasing and reduces the linewidth due to inhomogenous broadening. Mossbauer absorption experiments intended to test these ideas are outlined.

Jha, S.↗

High Order Modulation Protograph Codes

Digital communication coding methods for designing protograph-based bit-interleaved code modulation that is general and applies to any modulation. The general coding framework can support not only multiple rates but also adaptive modulation. The method is a two stage lifting approach. In the first stage, an original protograph is lifted to a slightly larger intermediate protograph. The intermediate protograph is then lifted via a circulant matrix to the expected codeword length to form a protograph-based low-density parity-check code.

Nguyen, Thuy V.↗

Graphene oxide/metal nanocrystal multilaminates the atomic limit for safe, selective hydrogen storage

Various embodiments of the invention describe an environmentally stable, and exceptionally dense hydrogen storage (6.5 wt % and 0.105 kg H2/L in the total composite, 7.56 wt % in Mg) using atomically thin and gas-selective reduced graphene oxide sheets as encapsulants. Other approaches to protecting reactive materials involve energy intensive introduction of considerable amounts of inactive, protective matrix which compromises energy density. However, these multilaminates are able to deliver exceptionally dense hydrogen storage far-exceeding 2020 DOE target metrics for gravimetric capacity (5.5 wt %), and ultimate full-fleet volumetric targets (0.070 kg H 2 /L) for fuel cell electric vehicles. Methods of stabilizing reactive nanocrystalline metals in zero-valency also has wide-ranging applications for batteries, catalysis, encapsulants, and energetic materials.

Urban, Jeffrey J.↗

Design and performance testing of a novel far field Gamma system to assay radioactive waste in 400 and 200 litre drums - 15339

A Far Field Gamma ray Measurement system has been developed for measuring and assaying both 400 and 200 litre waste drums with a wide range of activity at nuclear power stations. These include radioactive waste drums with surface dose-rates of less than 2 mSv/h (200 mrem/h) and, at the higher range of activity, 400 litre drums with surface dose-rates of up to 50 mSv/h (5 rem/h). The assay system incorporates an electro-mechanically cooled high purity Germanium profile detector and advanced digital spectroscopic analysis electronics. A novel optimised conical-trapezoidal collimator has been incorporated to reduce the effect of the radioactive background in addition to the conventional detector lead shield with tin and copper graded lining to reduce the effect of lead X-rays. Two features have been incorporated to reduce detector dead time and maintain detector resolution when measuring high activity drums. The first is the inclusion of two automated tungsten filters of different thickness as an integral part of the detector collimator. The filters are used in conjunction with an automated rail system on which the detector platform is mounted. The detector can be positioned at different distances from the drum surface in order to reduce detector dead time. The deployment of the filters and the adjustment of the drum - detector position are automatically controlled based on a user defined dose-rate limit table. The data is provided by two Geiger-Muller dose-rate detectors, which are used to measure the drum surface dose-rate. This data is subsequently used to adjust the detector position and in cases of higher dose-rate deploy one or other of the tungsten filters. The drum rotation platform incorporates a load cell to determine drum weight and hence the drum density. This information is used to calculate an attenuation correction based on the assumption of uniform drum density and uniform distribution of activity within the drum. Where the waste drums have a more complicated and known regular internal structure, such as a small radioactive region surrounded by an annulus of shielding material, the analysis algorithms incorporated into the far field geometry and attenuation correction spectroscopic data analysis code are able to make a more accurate determination of drum activity than would be possible assuming a completely uniform drum. The detector energy and efficiency calibration is achieved using a point source with multiple gamma ray energy peaks. The response of the high purity Germanium detector is also modelled using the Monte Carlo Neutron Photon code and the model is benchmarked using the measured point source calibration data. The performance of the instrument has been determined for a range of uniform density matrix drums containing a set of Eu-152 line or rod sources located in re-entrant tubes positioned on an equal volume basis. When rotated, the test drums with volume distributed line sources simulate waste drums of different density each with a uniformly distributed radioactive source. The measurements have been validated using Monte Carlo simulations. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Wave function continuity and the diagonal Born-Oppenheimer correction at conical intersections

We demonstrate that though exact in principle, the expansion of the total molecular wave function as a sum over adiabatic Born-Oppenheimer (BO) vibronic states makes inclusion of the second-derivative nonadiabatic energy term near conical intersections practically problematic. In order to construct a well-behaved molecular wave function that has density at a conical intersection, the individual BO vibronic states in the summation must be discontinuous. When the second-derivative nonadiabatic terms are added to the Hamiltonian, singularities in the diagonal BO corrections (DBOCs) of the individual BO states arise from these discontinuities. In contrast to the well-known singularities in the first-derivative couplings at conical intersections, these singularities are non-integrable, resulting in undefined DBOC matrix elements. Though these singularities suggest that the exact molecular wave function may not have density at the conical intersection point, there is no physical basis for this constraint. Instead, the singularities are artifacts of the chosen basis of discontinuous functions. We also demonstrate that continuity of the total molecular wave function does not require continuity of the individual adiabatic nuclear wave functions. We classify nonadiabatic molecular dynamics methods according to the constraints placed on wave function continuity and analyze their formal properties. Based on our analysis, it is recommended that the DBOC be neglected when employing mixed quantum-classical methods and certain approximate quantum dynamical methods in the adiabatic representation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Thermofield Theory for Finite-Temperature Electronic Structure

Wave function methods have offered a robust, systematically improvable means to study ground-state properties in quantum many-body systems. Theories like coupled cluster and their derivatives provide highly accurate approximations to the energy landscape at a reasonable computational cost. Analogues of such methods to study thermal properties, though highly desirable, have been lacking because evaluating thermal properties involve a trace over the entire Hilbert space, which is a formidable task. Besides, excited-state theories are generally not as well studied as ground-state ones. In this mini-review, we present an overview of a finite-temperature wave function formalism based on thermofield dynamics to overcome these difficulties. Thermofield dynamics allows us to map the equilibrium thermal density matrix to a pure state, i.e., a single wave function, albeit in an expanded Hilbert space. Ensemble averages become expectation values over this so-called thermal state. Around this thermal state, we have developed a procedure to generalize ground-state wave function theories to finite temperatures. As explicit examples, we highlight formulations of mean-field, configuration interaction, and coupled cluster theories for thermal properties of Fermions in the grand-canonical ensemble. To assess the quality of these approximations, we also show benchmark studies for the one-dimensional Hubbard model, while comparing against exact results. We will see that the thermal methods perform similarly to their ground-state counterparts, while merely adding a prefactor to the asymptotic computational cost. Furthermore, they also inherit all the properties, good or bad, from the ground-state methods, signifying the robustness of our formalism and the scope for future development.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Hydrogen Permeability of Polymer Matrix Composites at Cryogenic Temperatures

This paper presents experimental methods and results of an ongoing study of the correlation between damage state and hydrogen gas permeability of laminated composite materials under mechanical strains and thermal loads. A specimen made from IM-7/977-2 composite material has been mechanically cycled at room temperature to induce microcrack damage. Crack density and tensile modulus were observed as functions of number of cycles. Damage development was found to occur most quickly in the off-axis plies near the outside of the laminate. Permeability measurements were made after 170,000 cycles and 430,000 cycles. Leak rate was found to depend on applied mechanical strain, crack density, and test temperature.

Grenoble, Ray W.↗

Crystal Plasticity Model of Reactor Pressure Vessel Embrittlement in GRIZZLY

The integrity of reactor pressure vessels (RPVs) is of utmost importance to ensure safe operation of nuclear reactors under extended lifetime. Microstructure-scale models at various length and time scales, coupled concurrently or through homogenization methods, can play a crucial role in understanding and quantifying irradiation-induced defect production, growth and their influence on mechanical behavior of RPV steels. A multi-scale approach, involving atomistic, meso- and engineering-scale models, is currently being pursued within the GRIZZLY project to understand and quantify irradiation-induced embrittlement of RPV steels. Within this framework, a dislocation-density based crystal plasticity model has been developed in GRIZZLY that captures the effect of irradiation-induced defects on the flow stress behavior and is presented in this report. The present formulation accounts for the interaction between self-interstitial loops and matrix dislocations. The model predictions have been validated with experiments and dislocation dynamics simulation.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Fully Recyclable CFRP Particles (CRADA 519)

The proposed resin-enhanced materials, which utilize fully recycled CFRP (carbon-fiber-reinforced polymer) particles, are enabled by PARC's proprietary Chemically Linked Particles Networks (13 granted patents and patent applications). PNNL and PARC have collaborated to develop manufacturing methods, such as vacuum-assisted wet compression molding and filament winding, for fabricating fiber-reinforced composites with the PARC matrix. PNNL has conducted various mechanical, physical, and thermal characterizations to compare their differences between the composites with baseline and PARC matrices. Those characterizations include density, void volume fraction, optical micrography, uniaxial tension, four-point bending/flexure, interlaminar shear, coefficient of thermal expansion, glass transition temperature, and curing kinetics.

36 MATERIALS SCIENCE↗

The Density Matrix of H20 - N2 In the Coordinate Representation: A Monte Carlo Calculation of the Far-Wing Line Shape

The far-wing line shape theory within the binary collision and quasistatic framework has been developed using the coordinate representation. Within this formalism, the main computational task is the evaluation of multidimensional integrals whose variables are the orientational angles needed to specify the initial and final positions of the system during transition processes. Using standard methods, one is able to evaluate the 7-dimensional integrations required for linear molecular systems, or the 7-dimensional integrations for more complicated asymmetric-top (or symmetric-top) molecular systems whose interaction potential contains cyclic coordinates. In order to obviate this latter restriction on the form of the interaction potential, a Monte Carlo method is used to evaluate the 9-dimensional integrations required for systems consisting of one asymmetric-top (or symmetric-top) and one linear molecule, such as H20-N2. Combined with techniques developed previously to deal with sophisticated potential models, one is able to implement realistic potentials for these systems and derive accurate, converged results for the far-wing line shapes and the corresponding absorption coefficients. Conversely, comparison of the far-wing absorption with experimental data can serve as a sensitive diagnostic tool in order to obtain detailed information on the short-range anisotropic dependence of interaction potentials.

Ma, Q.↗

Monte Carlo modeling of ultrasound probes for image guided radiotherapy

Purpose: To build Monte Carlo (MC) models of two ultrasound (US) probes and to quantify the effect of beam attenuation due to the US probes for radiation therapy delivered under real-time US image guidance. Methods: MC models of two Philips US probes, an X6-1 matrix-array transducer and a C5-2 curved-array transducer, were built based on their megavoltage (MV) CT images acquired in a Tomotherapy machine with a 3.5 MV beam in the EGSnrc, BEAMnrc, and DOSXYZnrc codes. Mass densities in the probes were assigned based on an electron density calibration phantom consisting of cylinders with mass densities between 0.2 and 8.0 g/cm{sup 3}. Beam attenuation due to the US probes in horizontal (for both probes) and vertical (for the X6-1 probe) orientation was measured in a solid water phantom for 6 and 15 MV (15 × 15) cm{sup 2} beams with a 2D ionization chamber array and radiographic films at 5 cm depth. The MC models of the US probes were validated by comparison of the measured dose distributions and dose distributions predicted by MC. Attenuation of depth dose in the (15 × 15) cm{sup 2} beams and small circular beams due to the presence of the probes was assessed by means of MC simulations. Results: The 3.5 MV CT number to mass density calibration curve was found to be linear with R{sup 2} > 0.99. The maximum mass densities in the X6-1 and C5-2 probes were found to be 4.8 and 5.2 g/cm{sup 3}, respectively. Dose profile differences between MC simulations and measurements of less than 3% for US probes in horizontal orientation were found, with the exception of the penumbra region. The largest 6% dose difference was observed in dose profiles of the X6-1 probe placed in vertical orientation, which was attributed to inadequate modeling of the probe cable. Gamma analysis of the simulated and measured doses showed that over 96% of measurement points passed the 3%/3 mm criteria for both probes placed in horizontal orientation and for the X6-1 probe in vertical orientation. The X6-1 probe in vertical orientation caused the highest attenuation of the 6 and 15 MV beams, which at 10 cm depth accounted for 33% and 43% decrease compared to the respective (15 × 15) cm{sup 2} open fields. The C5-2 probe in horizontal orientation, on the other hand, caused a dose increase of 10% and 53% for the 6 and 15 MV beams, respectively, in the buildup region at 0.5 cm depth. For the X6-1 probe in vertical orientation, the dose at 5 cm depth for the 3-cm diameter 6 MV and 5-cm diameter 15 MV beams was attenuated compared to the corresponding open fields to a greater degree by 65% and 43%, respectively. Conclusions: MC models of two US probes used for real-time image guidance during radiotherapy have been built. Due to the high beam attenuation of the US probes, the authors generally recommend avoiding delivery of treatment beams that intersect the probe. However, the presented MC models can be effectively integrated into US-guided radiotherapy treatment planning in cases for which beam avoidance is not practical due to anatomy geometry.

60 APPLIED LIFE SCIENCES↗

Thermal fluctuations and carrier localization induced by dynamic disorder in MAPbI 3 described by first-principles based tight-binding model

Halide perovskites are strongly influenced by large amplitude anharmonic lattice fluctuations at room temperature. We develop a tight binding model for dynamically disordered MAPbI$_3$ based on density functional theory (DFT) calculations to calculate electronic structure for finite temperature crystal structures at the length scale of thermal disorder and carrier localization. The model predicts individual Hamiltonian matrix elements and band structures with high accuracy, owing to the inclusion of additional matrix elements and descriptors for non-Coulombic interactions. We apply this model to electronic structure at length and time scales inaccessible to first principles methods, finding an increase in band gap, carrier mass, and the sub-picosecond fluctuations in these quantities with increasing temperature as well as the onset of carrier localization in large supercells induced by thermal disorder at 300 K. Finally, we identify the length scale $L^*= 5$ nm as the onset of localization in the electronic structure, associated with associated with decreasing band edge fluctuations, increasing carrier mass, and Rashba splitting approaching zero.

36 MATERIALS SCIENCE↗

WE-DE-BRA-03: Construction of An Ultrasound Guidance Platform for Image-Guided Radiotherapy with the Intent to Treat Transitional Cell Carcinoma

Purpose: Ultrasound (US) is a noninvasive, nonradiographic imaging technique with high spatial and temporal resolution that can be used for localizing soft-tissue structures and tumors in real-time during radiotherapy (inter- and intra-fraction). A detailed methodology integrating 3D-US within RT is presented. This method is easier to adopt into current treatment protocol than current US based systems and reduces user variability for image acquisition, thus eliminating transducer induced changes that limit CT planning system. Methods: We designed an in-house integrated US manipulator and platform to relate CT, 3D-US and linear accelerator coordinate systems. To validate the platform, an agar-based phantom with measured densities and speed-of-sound consistent with tissues surrounding the bladder, was rotated (0–45°) resulting in translations (up to 55mm) relative to the CT and US coordinate systems. After acquiring and integrating CT and US images into the treatment planning system, US-to-US and US-to-CT images were co-registered to re-align the phantom relative to the linear accelerator. Errors in the transformation matrix components were calculate to determine precision of this method under different patient positions. Results: Statistical errors from US-US registrations for different patient orientations ranged from 0.06–1.66mm for x, y, and z translational components, and 0.00–1.05° for rotational components. Statistical errors from US-CT registrations were 0.23–1.18mm for the x, y and z translational components, and 0.08–2.52° for the rotational components. Conclusion: Based on our result, this is consistent with currently used techniques for positioning prostate patients if couch re-positioning is less than a 5 degree rotation. We are now testing this on a dog patient to obtain both inter and intra-fractional positional errors. Additional design considerations include the future use of ultrasound-based functionality (photoacoustics, radioacoustics, Doppler) to monitor blood flow and hypoxia and/or in-vivo dosimetry for applications in other therapeutic techniques, such as hyperthermia, anti-angiogenesis, and particle therapy.

60 APPLIED LIFE SCIENCES↗

Pre-conceptual Development and characterization of an extruded graphite composite fuel for the TREAT Reactor

In an effort to explore fuel systems that are more robust under accident scenarios, the DOE-NE has identified the need to resume transient testing. The Transient Reactor Test (TREAT) facility has been identified as the preferred option for the resumption of transient testing of nuclear fuel in the United States. In parallel, NNSA’s Global Threat Reduction Initiative (GTRI) Convert program is exploring the needs to replace the existing highly enriched uranium (HEU) core with low enriched uranium (LEU) core. In order to construct a new LEU core, materials and fabrication processes similar to those used in the initial core fabrication must be identified, developed and characterized. In this research, graphite matrix fuel blocks were extruded and materials properties of were measured. Initially the extrusion process followed the historic route; however, the project was expanded to explore methods to increase the graphite content of the fuel blocks and explore modern resins. Materials properties relevant to fuel performance including density, heat capacity and thermal diffusivity were measured. The relationship between process defects and materials properties will be discussed.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Highly Uniform Resistive Switching in HfO 2 Films Embedded with Ordered Metal Nanoisland Arrays

Memristors enter a critical developmental stage where emerging large-scale integration methods face major challenges with severe switching instabilities in the oxide layer. In this work, the superior uniformity is achieved within HfO 2 films by embedding highly ordered metal nanoisland (NI) arrays. Embedded films exhibit a significant reduction in both SET and RESET while displaying enhanced uniformity in operating voltages and resistance states. This behavior is attributed to the concentration of electric fields along Pt and Ti NIs and their interactions with the surrounding oxide film matrix environment, which induce separate and distinct filamentary formation mechanisms that affect the stability. A method is reported to further optimize the uniformity of the SET voltage by translating the NI array position down the film-thickness dimension towards the bottom electrode. A comparison of the density and distribution of the oxygen vacancies responsible for the formation/dissolution of conducting filaments is made via combined electrostatic force microscopy and conductive atomic force microscopy (c-AFM) studies. Finally, complete observation of the morphological evolution of conducting filaments produced by Pt and Ti is enabled by 3D c-AFM nanotomography and cross-sectional scanning transmission electron microscopy–energy dispersive spectroscopy to provide direct correlations between NI-oxide interactions and overall switching performance.

36 MATERIALS SCIENCE↗