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,279 records · Page 71
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

Effect of heat input on the microstructure, residual stresses and corrosion resistance of 304L austenitic stainless steel weldments

Austenitic stainless steels are widely used in high performance pressure vessels, nuclear, chemical, process and medical industry due to their very good corrosion resistance and superior mechanical properties. However, austenitic stainless steels are prone to sensitization when subjected to higher temperatures (673 K to 1173 K) during the manufacturing process (e.g. welding) and/or certain applications (e.g. pressure vessels). During sensitization, chromium in the matrix precipitates out as carbides and intermetallic compounds (sigma, chi and Laves phases) decreasing the corrosion resistance and mechanical properties. In the present investigation, 304L austenitic stainless steel was subjected to different heat inputs by shielded metal arc welding process using a standard 308L electrode. The microstructural developments were characterized by using optical microscopy and electron backscattered diffraction, while the residual stresses were measured by X-ray diffraction using the sin{sup 2}ψ method. It was observed that even at the highest heat input, shielded metal arc welding process does not result in significant precipitation of carbides or intermetallic phases. The ferrite content and grain size increased with increase in heat input. The grain size variation in the fusion zone/heat affected zone was not effectively captured by optical microscopy. This study shows that electron backscattered diffraction is necessary to bring out changes in the grain size quantitatively in the fusion zone/heat affected zone as it can consider twin boundaries as a part of grain in the calculation of grain size. The residual stresses were compressive in nature for the lowest heat input, while they were tensile at the highest heat input near the weld bead. The significant feature of the welded region and the base metal was the presence of a very strong texture. The texture in the heat affected zone was almost random. - Highlights: • Effect of heat input on microstructure, residual stresses and corrosion is studied. • HAZ and width of dendrite in the welded region increase with heat input. • Residual stresses are tensile near the welded region after the highest heat input. • Welded region has the highest pit density after highest heat input. • Dendrites and δ-ferrite were highly oriented in the welded region.

36 MATERIALS SCIENCE↗
NASA NTRS2012

Demonstration of Minimally Machined Honeycomb Silicon Carbide Mirrors

Honeycomb silicon carbide composite mirrors are made from a carbon fiber preform that is molded into a honeycomb shape using a rigid mold. The carbon fiber honeycomb is densified by using polymer infiltration pyrolysis, or through a reaction with liquid silicon. A chemical vapor deposit, or chemical vapor composite (CVC), process is used to deposit a polishable silicon or silicon carbide cladding on the honeycomb structure. Alternatively, the cladding may be replaced by a freestanding, replicated CVC SiC facesheet that is bonded to the honeycomb. The resulting carbon fiber-reinforced silicon carbide honeycomb structure is a ceramic matrix composite material with high stiffness and mechanical strength, high thermal conductivity, and low CTE (coefficient of thermal expansion). This innovation enables rapid, inexpensive manufacturing. The web thickness of the new material is less than 1 millimeter, and core geometries tailored. These parameters are based on precursor carbon-carbon honeycomb material made and patented by Ultracor. It is estimated at the time of this reporting that the HoneySiC(Trademark) will have a net production cost on the order of $38,000 per square meter. This includes an Ultracor raw material cost of about $97,000 per square meter, and a Trex silicon carbide deposition cost of $27,000 per square meter. Even at double this price, HoneySiC would beat NASA's goal of $100,000 per square meter. Cost savings are estimated to be 40 to 100 times that of current mirror technologies. The organic, rich prepreg material has a density of 56 kilograms per cubic meter. A charred carbon-carbon panel (volatile organics burnt off) has a density of 270 kilograms per cubic meter. Therefore, it is estimated that a HoneySiC panel would have a density of no more than 900 kilograms per cubic meter, which is about half that of beryllium and about onethird the density of bulk silicon carbide. It is also estimated that larger mirrors could be produced in a matter of weeks. Each cell is completely uniform, maintaining the shape of the inserted mandrel. Furthermore, the layup creates pressure that insures node bond strength. Each node is a composite laminate using only the inherent resin system to form the bond. This contrasts starkly with the other known method of producing composite honeycomb, in which individual corrugations are formed, cured, and then bonded together in a secondary process. By varying the size of the mandrels within the layup, varying degrees of density can be achieved. Typical sizes are 3/8 and 3/16 in. (approximately 10 and 5 millimeters). Cell sizes up to 1 in. (approximately 25 millimeters) have been manufactured. Similarly, the shape of the core can be altered for a flexible honeycomb structure.

Goodman, William↗
DOE OSTI2024

Nonadiabatic dynamics with classical trajectories: The problem of an initial coherent superposition of electronic states

Advances in coherent light sources and development of pump–probe techniques in recent decades have opened the way to study electronic motion in its natural time scale. When an ultrashort laser pulse interacts with a molecular target, a coherent superposition of electronic states is created and the triggered electron dynamics is coupled to the nuclear motion. A natural and computationally efficient choice to simulate this correlated dynamics is a trajectory-based method where the quantum-mechanical electronic evolution is coupled to a classical-like nuclear dynamics. These methods must approximate the initial correlated electron–nuclear state by associating an initial electronic wavefunction to each classical trajectory in the ensemble. Different possibilities exist that reproduce the initial populations of the exact molecular wavefunction when represented in a basis. We show that different choices yield different dynamics and explore the effect of this choice in Ehrenfest, surface hopping, and exact-factorization-based coupled-trajectory schemes in a one-dimensional two-electronic-state model system that can be solved numerically exactly. Furthermore, this work aims to clarify the problems that standard trajectory-based techniques might have when a coherent superposition of electronic states is created to initialize the dynamics, to discuss what properties and observables are affected by different choices of electronic initial conditions and to point out the importance of quantum-momentum-induced electronic transitions in coupled-trajectory schemes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2000

Assessment of Internal Fabric of Particulate Materials

Particle arrangement and distribution within a soil matrix has long been recognized as having significant influence on the mechanical behavior of cohesionless soils. It is well known that two soil specimens having the same grain type (e.g., quartz, feldspar, etc.), same grain size distribution and relative density (or void ratio) can display completely different mechanical behavior. Because of the different fabric configurations in the otherwise similar specimens, they are likely to have different mechanical properties such as stress-strain response, anisotropy, dilatancy, etc. Soil Fabric is defined as the arrangement of particles, particle groups and associated pore space. In the literature, fabric analysis techniques are mainly classified as destructive (e.g., specimen stabilization, thin-sectioning, and microscopy), and nondestructive techniques (e.g., magnetic resonance imaging, ultrasonic testing, x-ray radiography, and computed tomography). Quantifying the void ratio and its distribution is the main parameter used to describe the fabric of particulate materials. This paper presents a comprehensive literature review of fabric analysis techniques applied to particulate materials. In addition, the results of a comprehensive investigation to quantify void ratio of sand specimens will be presented and discussed. The sand used in the experiments in a natural, uniform rounded to sub-rounded silica sand known as F-75 banding sand with mean particle size of 0.22 mm. Uniform specimens and specimens subjected to different axial-strain levels tested under triaxial and biaxial conditions are examined to evaluate void ratio evolution and its distribution using destructive thin-sectioning and nondestructive Computed Tomography (CT) techniques. Details of a new innovative polygon generation technique called Voronoi tessellation used to quantify void ratio of microscopic images of sand grains will be presented and compared to classical Oda's method. Finally, frequency distribution of void ratio calculated using CT analysis will be presented and discussed.

Alshibi, Khalid A.↗
DOE OSTI2026

Testing of High S Matrix Glasses to Expand DFHLW Glass Compositional Ranges (Rev.1)

Gaps in glass composition-property data for direct-feed high-level waste (DFHLW) have recently been identified. One such gap is the region of high sulfur solubility since previous, pretreated, high-level wastes contained very little sulfur. Filling this data gap will significantly broaden the range of process flowsheet options including minimal washing and will allow for optimized waste loading in DFHLW glasses. This report summarizes the data collected during the characterization of the DFHLW High S Glass Matrix (HS24). A glass matrix of 50 glass compositions was developed to evenly cover the DFHLW composition region for high sulfur glass. Matrix glasses were designed to expand the composition region outside the current component concentration and property limits so as to reduce uncertainties at the limits. The 50 matrix glasses were fabricated and tested for properties important to the success of DFHLW vitrification including: compositions, canister centerline cooling (CCC) crystallinity and isothermal crystallinity, density, viscosity, electrical conductivity (EC), product consistency test (PCT) response, toxicity, and sulfate solubility. Melter materials corrosion testing is reported elsewhere. These glasses were intentionally designed to have high SO 3 solubilities (0.7 to 2.2 SO 3 wt%) in compositional regions that had not been previously explored. Forty-eight glasses showed the measured SO 3 content retained >80% of the target SO 3 and the densities of all the glasses ranged from 2.49 g·cm -3 to 2.74 g·cm -3 . While the model predicted nepheline formation in 5 glasses, one of the tested 50 CCC glasses formed nepheline, and 35 glasses formed Cr-containing phases such as spinels and eskolaite. Only five glasses were amorphous after CCC treatment where 44 glasses with detectable crystals contained =10 wt% crystals and only one glass had > 10 wt% crystals. None of the glasses exceeded the allowable T 2% for spinel crystal formation at 950 ºC (i.e., no glasses had >2 wt% spinel at 950 ºC) during isothermal crystal fraction tests where 10 glasses showed no crystalline phases at or below 950 ºC. All the glasses (except one which failed being slightly lower than the target) satisfied the SO 3 constraint while 98 glasses did not meet the viscosity constraints and 4 failed the EC constraints. Six quenched (Q) and six CCC glasses failed the Defense Waste Processing Facility (DWPF) Environmental Assessment (EA) glass PCT threshold and 3 Q and 4 CCC failed the PCT design constraint. One glass exceeded the WTP delisting limits for Cr via EPA Method 1311 (i.e., Toxicity Characteristic Leaching Procedure, TCLP). It should be emphasized that some of these glasses were specifically designed to approach or even exceed certain property constraints, as filling data gaps in these regions will provide the greatest benefit for future model development by improving accuracy and reducing uncertainties. These insights will ultimately support the development of more robust glass formulation strategies, enabling higher waste loading, reducing operational risks, and expanding the processing envelope.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
NASA NTRS2013

The Digital Astronaut Project Computational Bone Remodeling Model (Beta Version) Bone Summit Summary Report

Under the conditions of microgravity, astronauts lose bone mass at a rate of 1% to 2% a month, particularly in the lower extremities such as the proximal femur [1-3]. The most commonly used countermeasure against bone loss in microgravity has been prescribed exercise [4]. However, data has shown that existing exercise countermeasures are not as effective as desired for preventing bone loss in long duration, 4 to 6 months, spaceflight [1,3,5,6]. This spaceflight related bone loss may cause early onset of osteoporosis to place the astronauts at greater risk of fracture later in their lives. Consequently, NASA seeks to have improved understanding of the mechanisms of bone demineralization in microgravity in order to appropriately quantify this risk, and to establish appropriate countermeasures [7]. In this light, NASA's Digital Astronaut Project (DAP) is working with the NASA Bone Discipline Lead to implement well-validated computational models to help predict and assess bone loss during spaceflight, and enhance exercise countermeasure development. More specifically, computational modeling is proposed as a way to augment bone research and exercise countermeasure development to target weight-bearing skeletal sites that are most susceptible to bone loss in microgravity, and thus at higher risk for fracture. Given that hip fractures can be debilitating, the initial model development focused on the femoral neck. Future efforts will focus on including other key load bearing bone sites such as the greater trochanter, lower lumbar, proximal femur and calcaneus. The DAP has currently established an initial model (Beta Version) of bone loss due to skeletal unloading in femoral neck region. The model calculates changes in mineralized volume fraction of bone in this segment and relates it to changes in bone mineral density (vBMD) measured by Quantitative Computed Tomography (QCT). The model is governed by equations describing changes in bone volume fraction (BVF), and rates of changes in bone cell populations that remove and replace bone in packets within the bone region. The DAP bone model is unique in several respects. In particular in takes former models of volume fraction changes one step higher in fidelity and separates BVF into separate equations for mineralized and osteoid volume fractions governed by a mineralization rate. This more closely follows the physiology of the remodeling unit cycles where bone is first resorbed and then followed by the action of osteoblasts to lay down collagen matrix which eventually becomes mineralized. In another respect, the modules allow the functional description of the time rate of change of other parameters and variables in the model during a computational simulation. More detailed description of the model, preliminary validation results, current limitation and caveats, and planned advancements are provided in sections 2 through 5. The DAP bone model is being developed primarily as a research tool, and not as a clinical tool like QCT. Even if it transitions to a clinical tool, it is not intended to replace QCT or any other clinical tool. Moreover, the DAP bone model does not predict bone fracture. Its purpose is to provide valuable additional data via "forward prediction" simulations for during and after spaceflight missions to gain insight on, (1) mechanisms of bone demineralization in microgravity, and (2) the volumetric changes at the various bone sites in response to in-flight and post-flight exercise countermeasures. This data can then be used as input to the Keyak [8] (or equivalent) FE analysis method to gain insight on how bone strength may change during and after flight. This information can also be useful to help optimize exercise countermeasure protocols to minimize changes in bone strength during flight, and improve regain of bone strength post-flight. To achieve this goal, the bone model will be integrated with DAP's exercise countermeasure models to simulate the effect of exercise prescriptions on preserving bone. More specifically, the model will accept loading history due to muscle and joint force on bone and produce quantified remodeling within the bone region under influence of the applied stress. Furthermore, because they tend to respond differently, the bone remodeling model includes both trabecular bone and cortical bone.

Pennline, James↗
DOE OSTI2021

Benchmarking the multiconfigurational Hartree method by the exact wavefunction of two harmonically trapped bosons with contact interaction

Highlights: • Comparison of exact and numerical single-particle density matrices. • Benchmark for multiconfigurational Hartree. • Convergence for large number of orbitals. We consider two bosons in a one-dimensional harmonic trap, interacting by a contact potential, and compare the exact solution of this problem to a self-consistent numerical solution by using the multiconfigurational time-dependent Hartree (MCTDH) method. We thereby benchmark the predictions of the MCTDH method with a few-body problem that has an analytical solution for the most commonly experimentally realized interaction potential in ultracold quantum gases. It is found that exact ground state energy and first order correlations are accurately reproduced by MCTDH up to the intermediate dimensionless coupling strengths corresponding to typical background scattering lengths of magnetically trapped ultracold dilute Bose gases. For larger couplings, established for example by (a combination of) Feshbach resonances and optical trapping, the MCTDH approach overestimates the depth of the trap-induced correlation dip of first order correlations in position space, as well as overestimates the fragmentation, defined as the average relative occupation of orbitals other than the energetically lowest one. We anticipate that qualitatively similar features in the correlation function may arise for larger particle numbers, paving the way for a quantitative assessment of the accuracy of MCTDH by experiments with ultracold atoms.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2013

Quantifying the Impact of Single Bit Flips on Floating Point Arithmetic

In high-end computing, the collective surface area, smaller fabrication sizes, and increasing density of components have led to an increase in the number of observed bit flips. If mechanisms are not in place to detect them, such flips produce silent errors, i.e. the code returns a result that deviates from the desired solution by more than the allowed tolerance and the discrepancy cannot be distinguished from the standard numerical error associated with the algorithm. These phenomena are believed to occur more frequently in DRAM, but logic gates, arithmetic units, and other circuits are also susceptible to bit flips. Previous work has focused on algorithmic techniques for detecting and correcting bit flips in specific data structures, however, they suffer from lack of generality and often times cannot be implemented in heterogeneous computing environment. Our work takes a novel approach to this problem. We focus on quantifying the impact of a single bit flip on specific floating-point operations. We analyze the error induced by flipping specific bits in the most widely used IEEE floating-point representation in an architecture-agnostic manner, i.e., without requiring proprietary information such as bit flip rates and the vendor-specific circuit designs. We initially study dot products of vectors and demonstrate that not all bit flips create a large error and, more importantly, expected value of the relative magnitude of the error is very sensitive on the bit pattern of the binary representation of the exponent, which strongly depends on scaling. Our results are derived analytically and then verified experimentally with Monte Carlo sampling of random vectors. Furthermore, we consider the natural resilience properties of solvers based on the fixed point iteration and we demonstrate how the resilience of the Jacobi method for linear equations can be significantly improved by rescaling the associated matrix.

Elliott, James J.↗
DOE OSTI2016

Micro arc oxidation of aluminum for nuclear applications

Dry cask storage and transport of used (spent) nuclear fuel rely upon the basket and shell of the cask to transfer the decay heat from the fuel assemblies. Cask licensees have considered the use of light metal (aluminum) oxidation to enhance the emissivity of aluminum alloy or aluminum based Metal Matrix Composites (MMC) used in the cask basket. However, traditional anodization techniques suffer from issues that make fabrication and regulatory acceptance problematic. The Micro Arc Oxidation (MAO) technology overcomes the problems associated with traditional anodization in the production of consistent, uniform and reliable high emissivity coatings of aluminum alloys. Samples of AA 2219 and an aluminum-based MMC were treated using a proprietary MAO process and subjected to a number of tests to confirm suitability of the coated material. The treated material was inspected using macro photography and Scanning Electron Microscopy (SEM). A consistent, uniform black, porous surface with a typical surface roughness of Ra = 7-8 μm was observed. Emissivity of treated surfaces with thicknesses ranging from 24 microns to 40 microns was measured using an integrating sphere. In all cases, the emissivity coefficient at room temperature in the ranges of 400 to 700 nanometers and 1000 to 1650 nanometers was above 0.8, reaching 0.94 depending on the thickness. The MAO process yield was tested by coating 400 spent fuel basket components to ensure a minimum specified thickness. 97% of the tested points met or exceeded the minimum thickness required. The defect density was very low and 100% of the components met the quality assurance criteria. In-service surface performance was tested by exposing the treated surface to a number of high temperature environments (air, vacuum, nitrogen) and thermal cycling in air. In all cases, the surface remained unchanged. Test samples were exposed to 90 Celsius degrees borated and demineralized water for 48 hours with no change to the sample treated surface. Finally, tests sample were exposed to 60 C. degrees borated water for 168 hours with no change to the sample treated surface. MAO technology appears to be a viable method to provide consistent, durable, and high emissivity coatings, that can be used in the fabrication and licensing of dry storage/transportation casks.

36 MATERIALS SCIENCE↗
DOE OSTI2014

Laboratory infrared spectroscopy of gaseous negatively charged polyaromatic hydrocarbons

Based largely on infrared spectroscopic evidence, polycyclic aromatic hydrocarbon (PAH) molecules are now widely accepted to occur abundantly in the interstellar medium. Laboratory infrared spectra have been obtained for a large variety of neutral and cationic PAHs, but data for anionic PAHs are scarce. Nonetheless, in regions with relatively high electron densities and low UV photon fluxes, PAHs have been suggested to occur predominantly as negatively charged ions (anions), having substantial influence on cloud chemistry. While some matrix spectra have been reported for radical anion PAHs, no data is available for even-electron anions, which are more stable against electron detachment. Here we present the first laboratory infrared spectra of deprotonated PAHs ([PAH-H]{sup –}) in the wavelength ranges between 6 and 16 μm and around 3 μm. Wavelength-dependent infrared multiple-photon electron detachment is employed to obtain spectra for deprotonated naphthalene, anthracene, and pyrene in the gas phase. Spectra are compared with theoretical spectra computed at the density functional theory level. We show that the relative band intensities in different ranges of the IR spectrum deviate significantly from those of neutral and positively charged PAHs, and moreover from those of radical anion PAHs. These relative band intensities are, however, well reproduced by theory. An analysis of the frontier molecular orbitals of the even- and odd-electron anions reveals a high degree of charge localization in the deprotonated systems, qualitatively explaining the observed differences and suggesting unusually high electric dipole moments for this class of PAH molecules.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2011

Continuous distribution model for the investigation of complex molecular architectures near interfaces with scattering techniques

Biological membranes are composed of a thermally disordered lipid matrix and therefore require non-crystallographic scattering approaches for structural characterization with x-rays or neutrons. Here we develop a continuous distribution (CD) model to refine neutron or x-ray reflectivity data from complex architectures of organic molecules. The new model is a flexible implementation of the composition-space refinement of interfacial structures to constrain the resulting scattering length density profiles. We show this model increases the precision with which molecular components may be localized within a sample, with a minimal use of free model parameters. We validate the new model by parameterizing all-atom molecular dynamics (MD) simulations of bilayers and by evaluating the neutron reflectivity of a phospholipid bilayer physisorbed to a solid support. The determination of the structural arrangement of a sparsely-tethered bilayer lipid membrane (stBLM) comprised of a multi-component phospholipid bilayer anchored to a gold substrate by a thiolated oligo(ethylene oxide) linker is also demonstrated. From the model we extract the bilayer composition and density of tether points, information which was previously inaccessible for stBLM systems. The new modeling strategy has been implemented into the ga lowbar refl reflectivity data evaluation suite, available through the National Institute of Standards and Technology (NIST) Center for Neutron Research (NCNR).

60 APPLIED LIFE SCIENCES↗
DOE OSTI2013

Infrared absorption of 3-propenonyl (⋅CH{sub 2}CHCO) radical generated upon photolysis of acryloyl chloride [CH{sub 2}CHC(O)Cl] in solid para-H{sub 2}

Irradiation at 193 nm of a p-H{sub 2} matrix containing acryloyl chloride CH{sub 2}CHC(O)Cl at 3.2 K yielded infrared absorption lines at 3143.6 (ν{sub 1}), 3057.0 (ν{sub 2}), 3048.0 (ν{sub 3}), 2103.1 (ν{sub 4}), 1461.0 (ν{sub 5}), 1349.8 (ν{sub 6}), 1223.7 (ν{sub 11}+ν{sub 12} or 2ν{sub 12}), 1092.8 (ν{sub 8}), 918.1 (ν{sub 9}), 691.0 (ν{sub 10}), 624.3 (ν{sub 11}), and 597.1 (ν{sub 12}) cm{sup −1} that are assigned to the 3-propenonyl (⋅CH{sub 2}CHCO) radical. The assignments are based on the photolytic behavior and a comparison of observed vibrational wavenumbers and infrared intensities with those predicted with the B3PW91/aug-cc-pVDZ method. The observation is consistent with a major radical formation channel CH{sub 2}CHCO + Cl followed by escape of the Cl atom from the original p-H{sub 2} cage. The observation of 3-propenonyl (⋅CH{sub 2}CHCO) radical but not 3-propenalyl (s-cis- or s-trans-CH{sub 2}CHĊO) radical indicates that the former is the most stable isomer and that the barrier heights for conversion from s-cis- or s-trans-CH{sub 2}CHĊO to ⋅CH{sub 2}CHCO are small.

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

Downward quantum learning from element 118: Automated generation of Fermi–Löwdin orbitals for all atoms

A new algorithm based on a rigorous theorem and quantum data computationally mined from element 118 guarantees automated construction of initial Fermi–Löwdin-Orbital (FLO) starting points for all elements in the Periodic Table. It defines a means for constructing a small library of scalable FLOs for universal use in molecular and solid-state calculations. The method can be systematically improved for greater efficiency and for applications to excited states such as x-ray excitations and optically silent excitations. FLOs were introduced to recast the Perdew–Zunger self-interaction correction (PZSIC) into an explicit unitarily invariant form. The FLOs are generated from a set of N quasi-classical electron positions, referred to as Fermi-Orbital descriptors (FODs), and a set of N-orthonormal single-electron orbitals. FOD positions, when optimized, minimize the PZSIC total energy. However, creating sets of starting FODs that lead to a positive definite Fermi orbital overlap matrix has proven to be challenging for systems composed of open-shell atoms and ions. The proof herein guarantees the existence of a FLOSIC solution and further guarantees that if a solution for N electrons is found, it can be used to generate a minimum of N – 1 and a maximum of 2 N – 2 initial starting points for systems composed of a smaller number of electrons. As a result, applications to heavy and super-heavy atoms are presented. All starting solutions reported here were obtained from a solution for element 118, Oganesson.

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

Interface properties and built-in potential profile of a LaCr O 3 / SrTi O 3 superlattice determined by standing-wave excited photoemission spectroscopy

LaCrO 3 (LCO)/SrTiO 3 (STO) heterojunctions are intriguing due to a polar discontinuity along [001], exhibiting two distinct and controllable charged interface structures [(LaO) + /(TiO 2 ) 0 and (SrO) 0 /(CrO 2 ) - ] with induced polarization, and a resulting depth-dependent potential. In this study, we have used soft- and hard-x-ray standing-wave excited photoemission spectroscopy (SW-XPS) to quantitatively determine the elemental depth profile, interface properties, and depth distribution of the polarization-induced built-in potentials. We observe an alternating charged interface configuration: a positively charged (LaO) + /(TiO 2 ) 0 intermediate layer at the LCO top /STO bottom interface and a negatively charged (SrO) 0 /(CrO 2 ) - intermediate layer at the STO top /LCO bottom interface. Using core-level SW data, we have determined the depth distribution of species, including through the interfaces, and these results are in excellent agreement with scanning transmission electron microscopy and electron energy-loss spectroscopy mapping of local structure and composition. SW-XPS also enabled deconvolution of the LCO and STO contributions to the valence-band (VB) spectra. Using a two-step analytical approach involving first SW-induced core-level binding-energy shifts and then VB modeling, the variation in potential across the complete superlattice is determined in detail. As a result, this potential is in excellent agreement with density functional theory models, confirming this method as a generally useful tool for interface studies.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Thermal nature of confining strings

We investigate the quantum statistical properties of the confining string connecting a static fermion-antifermion pair in the massive Schwinger model. By analyzing the reduced density matrix of the subsystem located in between the fermion and antifermion, we demonstrate that as the interfermion separation approaches the string-breaking distance, the overlap between the microscopic density matrix and an effective thermal density matrix exhibits a pronounced, narrow peak, approaching unity at the onset of string breaking. This behavior reveals that the confining flux tube evolves toward a genuinely thermal state as the separation between the charges grows, even in the absence of an external heat bath. In other words, one cannot tell whether a reduced state of the subsystem arises from a surrounding heat bath or from entanglement with the rest of the system. The entanglement spectrum near the critical string-breaking distance exhibits a rapid transition from the dominance of a single state describing the confining electric string towards a strongly entangled state containing virtual fermion-antifermion pairs. Our findings establish a quantitative link between confinement, entanglement, and emergent thermality, and suggest that string breaking corresponds to a microscopic thermalization transition within the flux tube.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Second harmonic generation in topological insulators under quantizing magnetic fields

We theoretically investigate the second harmonic generation (SHG) of topological insulator surface states in a perpendicular magnetic field. Our theory is based on the microscopic expression of the second-order magneto-optical conductivity developed from the density matrix formalism, taking into account hexagonal warping effects on the surface states' band structure. Using numerically exact Landau-level energies and wave functions including hexagonal warping, we calculate the spectrum of SHG conductivities under normal incidence for different values of magnetic field and chemical potential. The imaginary parts of the SHG conductivities show prominent resonant peaks corresponding to one-photon and two-photon inter-Landau-level transitions. Treating the hexagonal warping term perturbatively, these transitions are clarified analytically within a perturbation theory from which approximate selection rules for the allowable optical transitions for SHG are determined. Furthermore, our results show extremely high SHG susceptibility that is easily tunable by magnetic field and doping level for topological surface states in the far-infrared regime, exceeding that of many conventional nonlinear materials. This work highlights the key role of hexagonal warping effects in generating second-order optical responses and provides insights into the nonlinear magneto-optical properties of the topological insulators.

Landau levels↗
DOE OSTI2016

Hydrodynamic description of an unmagnetized plasma with multiple ion species. I. General formulation

A generalization of the Braginskii ion fluid description [S. I. Braginskii, Sov. Phys. JETP 6, 358 (1958)] to the case of an unmagnetized collisional plasma with multiple ion species is presented. An asymptotic expansion in the ion Knudsen number is used to derive the individual ion species continuity, as well as the total ion mass density, momentum, and energy evolution equations accurate through the second order. Expressions for the individual ion species drift velocities with respect to the center of mass reference frame, as well as for the total ion heat flux and viscosity, which are required to close the fluid equations, are evaluated in terms of the first-order corrections to the lowest order Maxwellian ion velocity distribution functions. A variational formulation for evaluating such corrections and its relation to the plasma entropy are presented. Employing trial functions for the corrections, written in terms of expansions in generalized Laguerre polynomials, and maximizing the resulting functionals produces two systems of linear equations (for “vector” and “tensor” portions of the corrections) for the expansion coefficients. A general matrix formulation of the linear systems as well as expressions for the resulting transport fluxes are presented in forms convenient for numerical implementation. The general formulation is employed in the companion paper [A. N. Simakov and K. Molvig, Hydrodynamic description of an unmagnetized plasma with multiple ion species. II. Two and three ion species plasmas, submitted to Phys. Plasmas (2015)] to evaluate the individual ion drift velocities and the total ion heat flux and viscosity for specific cases of two and three ion species plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

SU-E-T-516: Investigation of a Novel Radiochromic Radiation Reporting System Utilizing the Reduction of Ferric Ion

Purpose To introduce and characterize a new “reverse-Fricke” radiation reporting system utilizing the reduction of ferric ions (Fe{sup 3+}) to ferrous ions (Fe{sup 2+}). Methods Two formulations of the radiochromic reporting system, referred to as A and B, were prepared for investigation. Formulation-A consisted of 14 mM 1,10-phenanthroline, 42 mM ethanol, and 57 mM ammonium ferric oxalate in water. Formulation-B consisted of 27 mM 1,10-phenanthroline, 42 mM ethanol, and 28 mM ammonium ferric oxalate in water. Solutions were prepared immediately prior to irradiation with a Cobalt-60 unit with radiation doses of 0, 1, 5, 10, 15, 20, and 25 Gy. The change in optical density over the visible range of 450–650 nm was measured using a spectrophotometer immediately after irradiation. The effective atomic numbers of the formulations were calculated using Mayneord’s formula. Results Ionizing radiation energy absorbed in the solutions causes the reduction of ferric ions (Fe{sup 3+}) into ferrous ions (Fe{sup 2+}), which then forms a 1:3 red colored complex with 1,10-phenanthroline ([(C{sub 1} {sub 2}H{sub 8}N{sup 2}){sub 3}Fe]{sup 2+}) that can be measured spectrophotometrically. The absorbance spectra of the resulting complex displayed a peak maximum at 512 nm with a greater change in absorbance for Formulation-B after receiving comparable radiation doses. The change in absorbance relative to dose exhibited a linear response up to 25 Gy for both Formulation-A (R{sup 2} = 0.98) and Formulation-B (R{sup 2} = 0.97). The novel formulations were also nearly water equivalent (Zeff = 7.42) with effective atomic numbers of 7.65 and 7.52 and mass densities within 0.2% of water. Conclusion Both formulations displayed visible Fe{sup 2+} complex formation with 1,10-phenanthroline after irradiation using a Cobalt-60 source. The higher sensitivity measured for Formulation-B is attributed to the increase in 1,10-phenanthroline concentration and the increase in the 1,10-phenanthroline to ammonium ferric oxalate ratio. Further investigation of this radiation reporting system in a 3D matrix material is encouraged. NSF GRFP Grant Award #LH-102SPS.

60 APPLIED LIFE SCIENCES↗