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 739 records · Page 41
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

Size-controlled synthesis of anatase nanobrush structures with higher crystal density

Highlights: • Application of AAO as template for TiO{sub 2} nanowires synthesis. • Non-hydrolytic sol-gel method as an implementation of TiO{sub 2} into AAO template. • Vertical anchorage of TiO2 nanowires to the polymer support. • Photocatalytic oxidation of methyl orange correlate with geometric surface area of TiO{sub 2}. - Abstract: Titanium dioxide was synthesized in multistep process applying of AAO (Anode Aluminum Oxide) as a template. The TiO{sub 2} nanowires were obtained by deposition of titanium precursors inside the pores of AAO membranes (non-hydrolytic sol-gel technique), calcined at 875 K with subsequent removal of alumina matrix by leaching it with phosphoric acid. During the AAO synthesis, it was proved that the electrical work can serve as the unified factor influencing the amount of formed AAO. The diameters of obtained titania nanowires were consequence of used AAO membranes. The higher density of titania elementary cell was observed for all synthesized nanostructures. It was confirmed by XRD and SAED measurements. The TiO{sub 2} nanowires attached vertically to polymer disc (50 mm diameter) were used in photocatalytic decomposition of methyl orange. All nanomaterials showed ten times higher photocatalytic activity than P-25. High photoactivity of obtained material is related to external area of vertically aligned titania nanowires.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Quantitative assessment of fitting errors associated with streak camera noise in Thomson scattering data analysis

Thomson scattering measurements in high energy density experiments are often recorded using optical streak cameras. In the low-signal regime, noise introduced by the streak camera can become an important and sometimes the dominant source of measurement uncertainty. In this paper, we present a formal method of accounting for the presence of streak camera noise in our measurements. We present a phenomenological description of the noise generation mechanisms and present a statistical model that may be used to construct the covariance matrix associated with a given measurement. This model is benchmarked against simulations of streak camera images. We demonstrate how this covariance may then be used to weight fitting of the data and provide quantitative assessments of the uncertainty in the fitting parameters determined by the best fit to the data and build confidence in the ability to make statistically significant measurements in the low-signal regime, where spatial correlations in the noise become apparent. These methods will have general applicability to other measurements made using optical streak cameras.

47 OTHER INSTRUMENTATION↗
DOE OSTI2020

Spin-phonon relaxation from a universal ab initio density-matrix approach

Designing new quantum materials with long-lived electron spin states urgently requires a general theoretical formalism and computational technique to reliably predict intrinsic spin relaxation times. We present a new, accurate and universal first-principles methodology based on Lindbladian dynamics of density matrices to calculate spin-phonon relaxation time of solids with arbitrary spin mixing and crystal symmetry. This method describes contributions of Elliott-Yafet and D’yakonov-Perel’ mechanisms to spin relaxation for systems with and without inversion symmetry on an equal footing. We show that intrinsic spin and momentum relaxation times both decrease with increasing temperature; however, for the D’yakonov-Perel’ mechanism, spin relaxation time varies inversely with extrinsic scattering time. We predict large anisotropy of spin lifetime in transition metal dichalcogenides. The excellent agreement with experiments for a broad range of materials underscores the predictive capability of our method for properties critical to quantum information science.

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

Nanostructural Characterization of the Fe{sub 3}O{sub 4}/ZnO Magnetic Nanocomposite as an Application in Medicine

In the present work, we synthesize Fe{sub 3}O{sub 4}/ZnO nanocomposite by sol-gel method using iron nitrate, ethylene glycol, zinc acetate dehydrate, and ammonium carbonate as precursors. Amorphous matrix and crystalline grains were formed in composite as well. The nanostructural properties of the nanocomposite were characterized by X-ray diffraction (XRD), Fourier transform infrared (FTIR), X-ray photoelectron spectroscopy (XPS), atomic force microscopy (AFM), transmission electron microscopy (TEM), vibrating sample magnetometer (VSM), and colloidal stability was checked out by Zeta potential. The biological effects of Fe{sub 3}O{sub 4}/ZnO nanocomposite were investigated by microculture tetrazolium test (MTT) against MCF-7 breast cancer cells as usage in medicine. Crystallite size, lattice strain, and dislocation density of the nanoparticles were estimated with Williamson-Hall method. The obtained results show that the highest death rate of Michigan Cancer Foundation-7 (MCF-7) cancer cells was for Fe{sub 3}O{sub 4}/ZnO nanocomposite with the molar ratio is 1:10 in the lower concentration (49 μ g/ml). Therefore, these nanocomposites due to their small size, colloidal stability, magnetic properties, and nontoxicity can be a good candidate as option in medicine application.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2022

Cold spray manufacturing of oxide-dispersion strengthened (ODS) steels using gas-atomized and ball-milled 14YWT powders

The cold spray deposition process has been investigated for the manufacture of 14YWT oxide-dispersion strengthened (ODS) steel, a nanostructured ferritic alloy (NFA), using gas-atomized and ball-milled feedstock powders. Cold spraying of the gas-atomized powder resulted in a thick dense deposit, but annealing them above 900 °C induced grain growth and heterogeneous precipitation of nanoparticles, thus resulting in lower hardness (~ 170 HV) compared to 14YWT ODS steel manufactured by the conventional method. The ball-milled powder (Fe-14Cr-3W-0.4Ti-0.2Y-0.125O) was cryomilled in liquid N 2 to achieve small particle sizes followed by annealing in a H 2 /Ar environment at 1000 °C and 1100 °C to reduce the hardness from ~ 8 GPa to ~ 5 GPa. Both smaller particle sizes and lower hardness of annealed powders made them amenable to the cold spray deposition process, resulting in thick and dense deposits, while retaining the favorable fully solutionized microstructure of the ball-milled powders. The microstructure of the as-deposited ODS exhibited very fine grain sizes (50–250 nm) with a uniform dispersion of Y-Ti-O nanoprecipitates in the ferritic steel matrix. The hardness values of the cold spray deposits using the powders annealed at 1000 °C and 1100 °C was 770 HV and 496 HV, respectively. Post-heat treatment of the as-deposited material up to 1100 °C showed a high-density microstructure with uniformly dispersed nanoscale oxide particles. The hardness of the annealed deposits up to 1000 °C had a notably higher hardness than bulk 14YWT ODS steel manufactured by conventional methods. The study demonstrates that cold spray deposition of the microstructurally-tailored feedstock powders can be used to successfully manufacture ODS steel with microstructure and properties comparable to those produced by the conventional methods, while providing for an attractive option for the more rapid and cost-effective manufacturing of ODS steel cladding tubes.

14YWT↗
DOE OSTI2021

Impact of three-body forces on elastic nucleon-nucleus scattering observables

In a previous series of papers we investigated the domain of applicability of chiral potentials to the construction of a microscopic optical potential (OP) for elastic nucleon-nucleus scattering. The OP was derived at the first order of the spectator expansion of the Watson multiple scattering theory and its final expression was a folding integral between the nucleon-nucleon ( N N ) t matrix and the nuclear density of the target. In the calculations N N and three-nucleon ( 3 N ) chiral interactions were used for the target density and only the N N interaction for the N N t matrix. The purpose of this study is to achieve another step towards the calculation of a more consistent OP introducing the 3 N force also in the dynamic part of the OP. The full treatment of the 3 N interaction is beyond our present capabilities. Thus, here it is approximated with a density dependent N N interaction obtained after the averaging over the Fermi sphere. In practice, in our model the 3 N force acts as a medium correction of the bare N N interaction used to calculate the t matrix. Even if the 3 N force is treated in an approximate way, this method naturally extends our previous model of the OP and allows a direct comparison of our present and previous results. We consider as case studies the elastic scattering of nucleons off C 12 and O 16 . We present results for the differential cross section and the spin observables for different values of the projectile energy. From the comparison with the experimental data and with the results of our previous model we assess the importance of the 3 N interaction in the dynamic part of the OP. Conclusions: Our analysis indicates that the contribution of the 3 N force in the t matrix is small for the differential cross section and it is sizable for the spin observables, in particular, for the analyzing power. We find that the two-pion exchange term is the major contributor to the 3 N force. A chiral expansion order-by-order analysis of the scattering observables confirms the convergence of our results at the next-to-next-to-next-to-leading-order, as already established in our previous work.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2015

SU-E-I-45: Feasibility for Using Iodine Quantification to Assist Diagnosis in Dual Energy Contrast-Enhanced Digital Mammography

Purpose: The objective of this study is to develop quantitative calibration between image quality indexes and iodine concentration with dual-energy (DE) contrast-enhanced digital mammography (CEDM) techniques and further serve as the assistance for diagnosis. Methods: Custom-made acrylic phantom with dimensions of 24×30 cm{sup 2} simulated breast thickness from 2 to 6 cm was used in the calibration. The phantom contained matrix of four times four holes of 3 mm deep with a diameter of 15 mm for filling contrast agent with area density ranged from 0.1 to 10 mg/cm{sup 2}. All the image acquisitions were performed on a full-field digital mammography system (Senographe Essential, GE) with dual energy acquisitions. Mean pixel value (MPV), and contrast-to-noise ratio (CNR) were used for evaluating the relationship between image quality indexes and iodine concentration. Iodine map and CNR map could further be constructed with these calibration curves applied pixel by pixel utilized MATLAB software. Minimum iodine concentration could also be calculated with the visibility threshold of CNR=5 according the Rose model. Results: When evaluating the DE subtraction images, MPV increased linearly as the iodine concentration increased with all the phantom thickness surveyed (R{sup 2} between 0.989 and 0.992). Lesions with increased iodine uptake could thus be enhanced in the color-encoded iodine maps, and the mean iodine concentration could be obtained through the ROI measurements. As for investigating CNR performance, linear relationships were also shown between the iodine concentration and CNR (R{sup 2} between 0.983 and 0.990). Minimum iodine area density of 1.45, 1.73, 1.80, 1.73 and 1.72 mg/cm{sup 2} for phantom thickness of 2, 3, 4, 5, 6 cm were calculated based on Rose’s visualization criteria. Conclusion: Quantitative calibration between image quality indexes and iodine concentrations may further serving as the assistance for analyzing contrast enhancement for patient participating the dual energy CEDM procedures.

60 APPLIED LIFE SCIENCES↗
NASA NTRS1976

Dielectric coated wire antennas

An electrically thin dielectric insulating shell on an antenna composed of electrically thin circular cylindrical wires is examined. A moment method solution is obtained, and the insulating shell is modeled by equivalent volume polarization currents. These polarization currents are related in a simple manner to the surface charge density on the wire antenna. In this way the insulating shell causes no new unknowns to be introduced, and the size of the impedance matrix is the same as for the uninsulated wires. The insulation is accounted for entirely through a modification of the symmetric impedance matrix. This modification influences the current distribution, impedance, efficiency, field patterns, and scattering properties. The theory is compared with measurement for dielectric coated antennas in air.

Richmond, J. H.↗
NASA NTRS2019

Benchmarking Mixed Mode Failure in Progressive Damage and Failure Analysis Methods

The verification and validation of progressive-damage-analysis finite element methods are difficult but critical tasks to undertake during their development. Verification exercises assess whether a predictive analysis tool produces results that are consistent with the fundamental concepts and assumptions of the tool under evaluation. Ideally, closed-form analytical solutions can be derived for which method verification results can be compared. Problems selected for computational tool verification are often simple and isolate individual features of the tool. In the case of progressive damage finite element methods, verifications should be performed to evaluate the ability of the model to predict the initiation of damage and its growth through the finite element mesh under a variety of conditions. Mabson et al. proposed a test case of a unidirectional, fiber-reinforced plate with a center crack subjected to tensile loads to evaluate matrix crack propagation predictions. The problem was modeled using the Abaqus Hashin continuum damage mechanics (CDM) model for fiber-reinforced composites. Different combinations of matrix strengths and element sizes were used in the simulations, and the results were compared to a closed-form solution based on linear elastic fracture mechanics (LEFM). It was determined that the Abaqus CDM model could predict the LEFM solution of Mode I cracks only when the finite element mesh density met specific requirements based on the material properties. This paper presents closed-form LEFM solutions for a center notch mixed mode (CNMM) verification problem. Parametric finite element analyses were developed using progressive damage analysis methods of both the Discrete Damage Mechanics (DDM) and CDM classes. The progressive damage analysis methods applied in the analyses of the CNMM problem include CompDam and the Floating Node Method. Analyses were conducted with various mode mixities and element sizes to verify that the damage models were working as intended and to identify any limits of applicability.

Leone, Frank↗
DOE OSTI2015

Implementation of a generalized coherent elastic scattering formulation for thermal neutron scattering analysis - 14487

The scattering interactions of thermal neutrons in matter are classified into inelastic and elastic components. In crystalline materials, coherent elastic scattering represents a dominant form of interaction. Traditionally, the treatment of this interaction in standard computer codes (e.g., the LEAPR module of the NJOY code system) was performed within the cubic approximation and restricted the analyst to a few materials that were hardwired into the code. In general, these restrictions can hinder the production of ENDF style thermal neutron scattering libraries for new materials outside the current portfolio of treated materials. In this work, a general formulation based on the exact theory of coherent elastic scattering was developed. This formulation is capable of processing the crystalline structure of any material. The cubic approximation is relaxed by accounting for the partial forces from asymmetric crystal structures in the calculation. The atom sites approximation is relaxed by utilizing a site-wise Debye-Waller factor. Specifically, the Debye-Waller factor is no longer regarded as a uniform constant across the crystal but treated as a matrix tensor that uses the integral of the partial phonon density of states. Therefore, the calculation of the coherent elastic scattering cross section is executed by applying the matrix-form of the Debye-Waller factor on each reciprocal space point and summing them up. Furthermore, the developed algorithm supports compounds with any number of elements and of any ratio to each other. This general formulation was successfully implemented and contributed to the production of ENDF style thermal neutron scattering libraries for silicon dioxide (α-SiO{sub 2}) and silicon carbide (3C-SiC). The needed input for the analysis was generated using density functional theory (DFT) atomistic simulations and dynamical matrix calculations for both materials. (authors)

42 ENGINEERING↗
DOE OSTI2018

A Finite-Field Approach for GW Calculations beyond the Random Phase Approximation

In this work, we describe a finite-field approach to compute density response functions, which allows for efficient G0W0 and G0W0Γ0 calculations beyond the random phase approximation. The method is easily applicable to density functional calculations performed with hybrid functionals. Lastly, we present results for the electronic properties of molecules and solids and we discuss a general scheme to overcome slow convergence of quasiparticle energies obtained from G0W0Γ0 calculations, as a function of the basis set used to represent the dielectric matrix.

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

Ab initio ultrafast spin dynamics in solids

Spin relaxation and decoherence is at the heart of spintronics and spin-based quantum information science. Currently, theoretical approaches that can accurately predict spin relaxation of general solids including necessary scattering pathways and are capable of nanosecond to millisecond simulation time are urgently needed. We present a first-principles real-time density-matrix approach based on Lindblad dynamics to simulate ultrafast spin dynamics for general solid-state systems. Through the complete first-principles descriptions of pump, probe, and scattering processes including electron-phonon, electron-impurity, and electron-electron scatterings with self-consistent electronic spin-orbit couplings, our method can directly simulate the ultrafast pump-probe measurements for coupled spin and electron dynamics over nanoseconds at any temperatures and doping levels. We first apply this method to a prototypical system GaAs and obtain excellent agreement with experiments. We found that the relative contributions of different scattering mechanisms and phonon modes differ considerably between spin and carrier relaxation processes. In sharp contrast to previous work based on model Hamiltonians, we point out that the electron-electron scattering is negligible at room temperature but becomes dominant at low temperatures for spin relaxation in n-type GaAs. We further examine ultrafast dynamics in novel spin-valleytronic materials: monolayer and bilayer WSe 2 with realistic defects. We find that spin relaxation is highly sensitive to local symmetry and chemical bonds around defects. For the bilayer WSe 2 , we identify the scattering pathways in ultrafast dynamics and determine relevant dynamical properties, essential to its utilization of unique spin-valley-layer locking effects. In conclusion, our work provides a predictive computational platform for spin dynamics in solids, which has potential for designing new materials ideal for spintronics and quantum information technology.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2018

Ab initio electron-phonon interactions using atomic orbital wave functions

The interaction among electrons and lattice vibrations determines key physical properties of materials, including their electrical and heat transport, excited electron dynamics, phase transitions, and superconductivity. We present an ab initio method that employs atomic orbital (AO) wave functions to compute the electron-phonon (e-ph) interactions in materials and interpolate the e-ph coupling matrix elements to fine Brillouin zone grids. We illustrate the numerical implementation of such AO-based e-ph calculations, and benchmark them against direct density functional theory calculations and Wannier function (WF) interpolation. The prime advantages of AOs over WFs for e-ph calculations are outlined. Since AOs are fixed basis functions associated with the atoms, they circumvent the need to generate a material-specific localized basis set with a trial-and-error approach, as is needed in WFs. Therefore, AOs are ideal to compute e-ph interactions in chemically and structurally complex materials for which WFs are challenging to generate, and are also promising for high-throughput materials discovery. Although our findings focus on AOs, the formalism we present generalizes e-ph calculations to arbitrary localized basis sets, with WFs recovered as a special case.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS2005

Energy Finite Element Analysis for Computing the High Frequency Vibration of the Aluminum Testbed Cylinder and Correlating the Results to Test Data

The Energy Finite Element Analysis (EFEA) is a finite element based computational method for high frequency vibration and acoustic analysis. The EFEA solves with finite elements governing differential equations for energy variables. These equations are developed from wave equations. Recently, an EFEA method for computing high frequency vibration of structures either in vacuum or in contact with a dense fluid has been presented. The presence of fluid loading has been considered through added mass and radiation damping. The EFEA developments were validated by comparing EFEA results to solutions obtained by very dense conventional finite element models and solutions from classical techniques such as statistical energy analysis (SEA) and the modal decomposition method for bodies of revolution. EFEA results have also been compared favorably with test data for the vibration and the radiated noise generated by a large scale submersible vehicle. The primary variable in EFEA is defined as the time averaged over a period and space averaged over a wavelength energy density. A joint matrix computed from the power transmission coefficients is utilized for coupling the energy density variables across any discontinuities, such as change of plate thickness, plate/stiffener junctions etc. When considering the high frequency vibration of a periodically stiffened plate or cylinder, the flexural wavelength is smaller than the interval length between two periodic stiffeners, therefore the stiffener stiffness can not be smeared by computing an equivalent rigidity for the plate or cylinder. The periodic stiffeners must be regarded as coupling components between periodic units. In this paper, Periodic Structure (PS) theory is utilized for computing the coupling joint matrix and for accounting for the periodicity characteristics.

Vlahopoulos, Nickolas↗
NASA NTRS2013

T-Matrix Modeling of Linear Depolarization by Morphologically Complex Soot and Soot-Containing Aerosols

We use state-of-the-art public-domain Fortran codes based on the T-matrix method to calculate orientation and ensemble averaged scattering matrix elements for a variety of morphologically complex black carbon (BC) and BC-containing aerosol particles, with a special emphasis on the linear depolarization ratio (LDR). We explain theoretically the quasi-Rayleigh LDR peak at side-scattering angles typical of low-density soot fractals and conclude that the measurement of this feature enables one to evaluate the compactness state of BC clusters and trace the evolution of low-density fluffy fractals into densely packed aggregates. We show that small backscattering LDRs measured with groundbased, airborne, and spaceborne lidars for fresh smoke generally agree with the values predicted theoretically for fluffy BC fractals and densely packed near-spheroidal BC aggregates. To reproduce higher lidar LDRs observed for aged smoke, one needs alternative particle models such as shape mixtures of BC spheroids or cylinders.

atmospheric radiation↗
DOE OSTI2018

Compositions of particles comprising rare-earth oxides in a metal alloy matrix and related methods

A composition includes a metal alloy matrix comprising iron and a plurality of nanoparticles dispersed within the metal alloy matrix. Each nanoparticle of the plurality comprises an oxide of a rare-earth metal and at least one metal selected from the group consisting of tantalum, niobium, vanadium, and titanium. Some compositions include a metal alloy matrix comprising iron and a plurality of nanoparticles comprising at least two different oxides of rare-earth metals dispersed within the metal alloy matrix. Some methods include mixing an oxide of a rare-earth metal with a first metal and a second metal. Other methods include mixing a plurality of particles comprising at least one oxide of a rare-earth metal with a molten metal comprising iron. Each particle of the plurality may exhibit a density between about 6.9 g/cm.sup.3 and about 9.0 g/cm.sup.3.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Analytical derivatives of the individual state energies in ensemble density functional theory. II. Implementation on graphical processing units (GPUs)

Conical intersections control excited state reactivity, and thus, elucidating and predicting their geometric and energetic characteristics are crucial for understanding photochemistry. Locating these intersections requires accurate and efficient electronic structure methods. Unfortunately, the most accurate methods (e.g., multireference perturbation theories such as XMS-CASPT2) are computationally challenging for large molecules. The state-interaction state-averaged restricted ensemble referenced Kohn–Sham (SI-SA-REKS) method is a computationally efficient alternative. The application of SI-SA-REKS to photochemistry was previously hampered by a lack of analytical nuclear gradients and nonadiabatic coupling matrix elements. We have recently derived analytical energy derivatives for the SI-SA-REKS method and implemented the method effectively on graphical processing units. We demonstrate that our implementation gives the correct conical intersection topography and energetics for several examples. Furthermore, our implementation of SI-SA-REKS is computationally efficient, with observed sub-quadratic scaling as a function of molecular size. This demonstrates the promise of SI-SA-REKS for excited state dynamics of large molecular systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗