Search NASA⌕ Search

SEARCH · Search NASA

Results for “Material control and accounting”

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 397 records · Page 22

Phosphate adsorption kinetics and equilibria on natural iron and manganese oxide composites

We report that it is well known that phosphate retention in soils and sediments is strongly influenced by binding to secondary iron oxides, there have been relatively few studies examining its adsorption/desorption behavior on multicomponent particles of realistic natural complexity. In this study, natural Mn-rich limonite (LM), was used to prepare naturally complex Fe- and Mn-oxide composite materials to examine phosphate adsorption/desorption. To clarify the role of the Mn-oxides, results for the LM sample were compared to those for an acid treated version (LAT), in which the acid-extractable Mn-oxide fraction has been selectively eliminated while leaving the Fe-oxide fraction intact. The saturated adsorption capacity on LAT was almost double that on LM, suggesting that phosphate adsorption to the iron oxides is strongly occluded by the Mn-oxide fraction. This result is reinforced by the comparing the pH dependence and fits to adsorption isotherms, and by desorption experiments and STEM-EDS mapping showing that phosphate loading on Mn-oxides was limited. Hence, although the collective results confirm that phosphate uptake and strong binding is selectively controlled by the Fe-oxide fraction, our study reveals that the Mn-oxide fraction strongly interferes with this process. Therefore, phosphate uptake behavior on metal oxides cannot be predicted solely on the basis of the Fe-oxide fraction present, but instead must take into account the deleterious impacts of other intimately associated phases. For co-diagenetic Fe/Mn-oxide composites in particular, Mn-oxides appear to severely limit phosphate uptake on the Fe-oxide fraction, either by hindering access to binding sites on the Fe-oxide or by lowering their affinity for P.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Understanding the origin of softness in structurally tailored and engineered macromolecular (STEM) gels: A DPD study

Structurally tailored and engineered macromolecular (STEM) gels are polymer networks containing a primary network (the “STEM-0” gel) and side chains that are anchored to this network (forming the “STEM-1” gel). The ability to vary the features of the side chains provides a robust means of tailoring the macroscopic properties of the material. We use dissipative particle dynamics (DPD) simulations to determine the mechanical response of STEM-1 gels to uniaxial compression for various values of the side chain length, n sc . While the Young’s modulus of the material is significantly decreased with increasing n sc (at a fixed density of grafted side chains), above a certain saturation value, increasing n sc does not lead to any further softening of the sample. Using the simulations, we calculate the relevant stresses, which are related to the Young’s modulus and the number of entanglements in the network. We show that the backbone chains become more spread out upon addition of sidechains, leading to a decrease in the physical entanglements between backbone segments. This observation accounts for the decrease of stress in the backbone network and the softening of the networks through the addition of sidechains. For long chains, however, the number density of physical crosslinks between side chains is independent of side chain length; this observation explains the observed stress saturation in the STEM-1 gels. Our approach allows us to correlate the molecular architecture of the gels to the resultant macroscopic mechanical behavior and provide guidelines for fabricating STEM gels with well-defined mechanical properties, which allow the materials to be used for a variety of applications.

36 MATERIALS SCIENCE↗

Centrosymmetric or Noncentrosymmetric? Transition Metals Talking in K 2 TGe 3 S 8 (T = Co, Fe)

Two new quaternary sulfides K 2 TGe 3 S 8 (T=Co, Fe) have been synthesized by a high-temperature solid-state routine and flux growth method. The crystal growth process of K 2 TGe 3 S 8 (T=Co, Fe) was elucidated by in-situ powder X-ray diffraction and DSC thermal analysis. The mm-sized crystals of K 2 TGe 3 S 8 (T=Co, Fe) were grown. K 2 CoGe 3 S 8 crystallizes in a new structure type in centrosymmetric space group P1¯ (No. 2) with unit cell parameters of a = 7.016(1) Å, b = 7.770(1)Å, c = 14.342(1) Å, α = 93.80(1)°, β = 92.65(1)°, γ = 114.04(1)°. K 2 FeGe 3 S 8 crystallizes in K 2 FeGe 3 Se 8 structure type and the noncentrosymmetric space group P2 1 (No. 4) with unit cell parameters of a = 7.1089(5)Å, b = 11.8823(8)Å, c = 16.7588(11)Å, β = 96.604(2)°. There is a high structural similarity between K 2 CoGe 3 S 8 and K 2 FeGe 3 S 8 . The larger volume coupled with higher degrees of distortion of [FeS 4 ] tetrahedra compared to [CoS 4 ] tetrahedra accounts for the structure’s shift from centrosymmetric to noncentrosymmetric. The theory simulation confirms that [TS 4 ]T= Co or Fe tetrahedra play a crucial role in controlling the structure and properties of K 2 TGe 3 S 8 (T = Co, Fe). The measured optical bandgaps of K 2 CoGe 3 S 8 and K 2 FeGe 3 S 8 are 2.1(1) eV and 2.6(1) eV respectively. K 2 FeGe 3 S 8 shows antiferromagnetic ordering at 24K while no magnetic ordering was detected in K 2 CoGe 3 S 8 . In conclusion, the magnetic measurements also demonstrate the divalent nature of transition metals in K 2 TGe 3 S 8 (T = Co, Fe).

36 MATERIALS SCIENCE↗

On the elastic anisotropy of the entropy-stabilized oxide (Mg, Co, Ni, Cu, Zn)O compound

Herein, we study the elastic properties of the entropy-stabilized oxide (Mg, Co, Ni, Cu, Zn)O using experimental and first principles techniques. Our measurements of the indentation modulus on grains with a wide range of crystallographic orientations of the entropy-stabilized oxide revealed a high degree of elastic isotropy at ambient conditions. First principles calculations predict mild elastic anisotropy for the paramagnetic structure, which decreases when the system is considered to be non-magnetic. When the antiferromagnetic state of CoO, CuO, and NiO is accounted for in the calculations, a slight increase in elastic anisotropy is observed, suggesting a coupling between magnetic ordering and the orientation dependent elastic properties. Furthermore, an examination of the local structure reveals that the isotropy is favored through local ionic distortions of Cu and Zn—due to their tendencies to form tenorite and wurtzite phases. The relationships between the elastic properties of the multicomponent oxide and those of its constituent binary oxides are reviewed. These insights open up new avenues for controlling isotropy for technological applications through tuning composition and structure in the entropy-stabilized oxide or the high-entropy compounds in general.

Entropy↗

Pseudo-equilibrium theory for extrinsic doping control of the topological semimetal Cd 3 As 2

The standard approach for predicting defect equilibria from first principles assumes that the solid-state system is initially in a thermodynamic equilibrium with the external atomic reservoirs. This “growth step” is then often followed by a temperature quench in a “pseudo-equilibrium” in which some or all defect concentrations are frozen in until only the Fermi level E F remains to be equilibrated. However, this protocol does not account for the possibility of site exchanges which can create important defect redistributions as long as short-range defect migration is kinetically permissible. To model this redistribution, we developed an approach to solve for the non-equilibrium chemical potentials as a function of temperature while maintaining the overall defect stoichiometry. We then apply this approach to the Dirac semimetal Cd 3 As 2 to model extrinsic doping with group 1/11 and 14 elements. Undoped Cd 3 As 2 exhibits an undesirable mismatch between E F and the Dirac point. This unintentional electron doping originates from intrinsic defects and is difficult to overcome through adjustment of synthesis conditions alone. Employing our pseudo-equilibrium modeling, we identify extrinsic doping strategies for realizing doping-balanced Cd 3 As 2 at the relatively low temperatures accessible in thin-film growth of this material.

36 MATERIALS SCIENCE↗

Diagnostics: Chapter 8 of the special issue: on the path to tokamak burning plasma operation

This chapter presents the activity conducted by the ITPA topical group (TG) on Diagnostics over about the last 15 years. Following a general introduction of the ITER Diagnostics led by their measurement roles, the document is organized in several subchapters detailing the design support, research and development activity conducted by each of the specialist working groups (WGs) of the TG. Please note that the magnetic diagnostics were supported at the TG without a specific WG. Their status is included in the general introduction. In the following some highlights of the subchapter’s contents are provided. Recent advances in ITER first wall (FW) diagnostics for the measurements of plasma-metallic wall interaction in support of the ITER research plan are reported. An InfraRed imaging Video Bolometer for ITER has been developed and tested on several tokamaks to measure the radiated power loss. A laser-induced breakdown spectroscopy (LIBS) technique which utilizes a pulsed laser beam to ablate locally by forming a crater, will measure local tritium inventory in the FW material. Real-time Residual Gas Analyzers will measure the neutral gas composition in a divertor port and an equatorial port during plasma operation. Due to the full metallic FW environment, the plasma-wall interaction in ITER will face several challenges such as the compromised radiated power and divertor heat flux measurements by reflection. Ray tracing and analysis codes have been developed to eliminate and correct the effects of reflection in the measurements. The characteristics of the reflecting surfaces depending on the roughness and angle of the incidence have been measured by dedicated experiments, and the results were applied to the reflection elimination. For the measurement of the metallic impurity radiation induced by eroded metallic atoms, a vacuum ultraviolet spectrometer has been developed and tested. An extensive thermonuclear diagnostic suite will be required to support the operation of ITER and the planned experimental program for future burning plasma experiments. Due to the harsh environmental conditions, the implementation of diagnostic systems in ITER is a major challenge. These conditions include high levels of neutron and gamma fluxes, neutron heating, particle bombardment. Therefore, the selection and design of diagnostic systems must take into account a number of phenomena previously unseen in diagnostic design. For this reason, the measurement of neutrons and confined or lost fast ions, with particular emphasis on alpha particles, is critical to ITER. The diagnostics associated with these measurements will be important for future plasma-burning experiments at ITER. The high neutron emission and very large plasma size in ITER make neutron diagnostics the main diagnostic method used to measure plasma parameters such as fusion power, fusion power density, ion temperature, energy of fast ions and their spatial distributions in the plasma core. Active spectroscopy techniques are methods where a neutral particle beam is injected into the plasma and information on plasma parameters is extracted from the measurement of line emission resulting from the beam-plasma interaction, either by plasma ions or by beam atoms. Spatial localization is achieved by crossing the beamline and multiple observation lines. The ITER plasma will be a high temperature, moderately dense, fully ionized collisional plasma. The plasma facing surfaces are principally metallic being fashioned from beryllium or tungsten but many other elements, arising from either structural or from operational needs, may enter this plasma. The energy range of the emitted photons range from meV (infra-red) to multi keV (x-rays) and originate from all areas of the plasma volume. The primary role of passive emission diagnostics is to identify what is in the plasma from spectral signatures. Extracting quantitative information from these measurements such as impurity content, ion temperature, rotation, degree of detachment and radiated power depends on calibrated instruments, a physics model of the atomic and molecular processes and plasma transport and an analysis workflow that takes into account environmental effects such as reflections. The particular needs for ITER have prompted a multi-machine, many-year effort to address all these aspects and this chapter reviews the work on diagnostic design, experiments and new analysis techniques. An overview of the laser diagnostics to be implemented on ITER is also provided in this paper. This includes descriptions of the Thomson scattering in the core, edge and divertor regions, polarimetry and interferometry diagnostics used for measuring plasma density and also measurements of helium density in the divertor using Laser Induced Flourescence. Techniques which can allow improvements on current measurements are also addressed in particular expanding poloidal polarimetry measurements to measure field fluctuations and proposed use of dispersion interferometery which has a number of advantages over existing methods. This paper identifies particular areas where further research and testing on existing tokamaks is useful even at this advanced stage to inform the design of diagnostics for ITER. Outstanding areas of concern for the implementation of laser diagnostics, in particular with a view to reliable operation are identified. An overview of the latest developments of microwave diagnostic systems and techniques is given. The primary focus is the contributions for ITER—the next step burning plasma experiment—which is supplemented by describing recent progress of techniques applicable for fusion experiments beyond ITER. The contributions are intentionally kept concise, and are being supplemented by a rich list of references for further studies. Radiation induced effects are receiving continuous and well-deserved attention of the ITER diagnostic community and they are in many cases one of the primary design drivers of the ITER diagnostic systems. The paper summarizes recent progress in this area focusing primarily on the ITER diagnostics but in some cases provides also outlook for the possible solutions for even more demanding radiation environment of fusion reactors beyond ITER. Despite advancements in the area of modeling and simulation of various radiation induced effects, experimental testing in a nuclear environment as close as possible to the target one is still seen as unavoidable for proper qualification of particular diagnostic functional elements. Recent advancement within three diagnostic areas: optical diagnostics, magnetics and bolometers is covered. Encouraging results on qualification of silica glass vacuum window assemblies are presented. In the area of magnetic sensors, progress of irradiation tests performed on ITER in-vessel LTCC inductive sensors is presented with outlook for novel technological approaches to inductive sensors utilizing thick printing and photolithography technologies being highlighted. Summary of advancements in the area of steady state magnetic field sensors based on Hall effect is given. New results of neutron irradiation test of the ITER borosilicate glass inserts for vacuum electrical feedthroughs are summarized finding negligible swelling at target level of neutron fluence. Off-line irradiation tests of fiber optic current sensors for plasma current measurement demonstrated that both for gamma doses up to 5 MGy and a total neutron fluence up to 10 15 cm −2 , radiation induced changes are still compatible with required measurement accuracy on ITER. The ITER bolometers are given as an example how considering radiation effects may influence the diagnostic design. Finally, outlook for future main R&D directions is outlined. All optical and laser-based diagnostics in ITER will be using mirrors to guide plasma radiation toward detectors, cameras and sensors. In the hostile plasma, radiation and particle environment the optical characteristics of diagnostic mirrors will degrade directly affecting the entire performance of involved diagnostic systems. An assessment of factors affecting mirror performance is provided. Among the prime adverse factors are deposition of plasma impurities, sputtering of mirror surface and steam ingress in the vicinity of mirrors. Within the International Tokamak Physics Activity with active support by ITER central team and domestic agencies, the structured research and development (R&D) program on mitigation of risks for diagnostic mirrors is underway. Within this program the mirror material development, the passive mitigation of mirror degradation by using diagnostic ducts and shutters along with an active mirror recovery program comprising the in-situ mirror cleaning and calibration is underway. Recent developments in diagnostic mirror R&D are described in this Chapter along with an example of their implementation of R&D solutions in ITER Infrared Thermography diagnostic. An assessment of still open engineering and physics questions, considerations on mirror risks during an early phase of ITER operation are given along with an overview of diagnostic mirror evolution in the late ITER operation stage toward the demonstration fusion power plant. Several crucial areas of diagnostic R&D outlined in ITER Research Plan are addressed. The basic control groups in a fusion reactor can be broken-down in five categories: (1) plasma position, magnetic configuration, and plasma current control, (2) profile control and confinement optimization, (3) MHD control and suppression, (4) edge dissipation control, radiation and plasma exhaust control and (5) break-down optimization. These categories are coupled via the physics (a control action in one domain will affect the other domains) and via shared actuators (e.g. ECRH for impurity accumulation avoidance, current density distribution control and MHD suppression). Consequently, a supervisory control system should determine the priority of the various control tasks, their couplings, and the interfaces with the safety and interlock system. For the systematic development of the various controllers taking the complexity of the plasma and the control system into account, a model-based approach is required. A short historical overview is given of the developments in systems and control theory and control engineering with special emphasis on those developments that are most relevant for Nuclear Fusion research and operation. An overview is given of the state of the field of fusion plasma control for the control categories. It will be shown how synthetic diagnostics are being developed in ITER and how they are used in diagnostic design and design validation and how they can be in model-based controller synthesis using relatively simple models. In modern control methods, multiple diagnostics are used to constrain relatively simple models. The constrained models provide an estimate for the state. This opens the route to state controllers, such as model predictive control. A major challenge in nuclear fusion research is the coherent combination of data from heterogeneous diagnostics and modeling codes for machine control and safety as well as physics studies. Measured data from different diagnostics often provide information about the same subset of physical parameters. Additionally, information provided by some diagnostics might be needed for the analysis of other diagnostics. A joint analysis of complementary and redundant data allows, e.g. to improve the reliability of parameter estimation, to increase the spatial and temporal resolution of profiles, to obtain synergistic effects, to consider diagnostics interdependencies and to find and resolve data inconsistencies. Physics-based modeling and parameter relationships provide additional information improving the treatment of ill-posed inversion problems. A coherent combination of all kind of available information within a probabilistic framework allows for improved data analysis results. The concept of integrated data analysis (IDA) in the framework of Bayesian probability theory is outlined and contrasted with conventional data analysis. Components of the probabilistic approach are summarized and specific ingredients beneficial for data analysis at fusion devices are discussed.

ITER↗

Computational Fermi Level Engineering and Doping-Type Conversion of Ga2O3 Via Three-Step Processing

Ga2O3 is being actively explored for power electronics, deep-ultraviolet optoelectronics, and other applications due to its ultra-wide bandgap and low projected fabrication cost of large-size and high-quality crystals. N-type doping of Ga2O3 can be achieved and tuned, but p-type doping faces fundamental obstacles due to deep character of acceptor levels and polaron transport of resulting holes. However, successful engineering of Ga2O3 based devices requires critical control of doping density, Fermi level position, and free carrier concentration, providing opportunities for predictive process simulation. We use first-principles defect theory and defect equilibrium calculations to simulate a 3-step growth-annealing-quench protocol for hydrogen assisted Mg doping in Ga2O3, taking into account the hydrogen-oxygen-water gas phase equilibrium. We predict type conversion to a net p-type regime following O-rich annealing after growth under reducing conditions in the presence of H2. This process is similar to the Mg acceptor activation by H removal in GaN. We show that there is an optimal temperature that maximizes the net acceptor density during the equilibrium annealing step for a given Mg doping level. Quenching of non-equilibrium annealed samples then results in a Fermi level EF below mid-gap down to about EV +1.5 eV, creating a significant number of uncompensated neutral MgGa0 acceptors. The resulting free hole concentration in quenched samples is very low (109 – 1012 cm-3) due to deep energy level of these Mg acceptors, but this type converted Ga2O3 material can create a significant built-in field in a p-n junction with an adjoining n-type material. Additionally, the electron concentration is greatly suppressed in such acceptor-doped Ga2O3, which could enable the use as a current blocking layer to fabricate normally-off (enhancement-mode) vertical Ga2O3 based metal-oxide-semiconductor field effect transistors (MOSFETs), and as a guard ring for edge termination in MOSFETs and Schottky barrier diodes (SBDs) with increased breakdown voltage.

CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSICS,M↗

Phase field dislocation dynamics (PFDD) modeling of non-Schmid behavior in BCC metals informed by atomistic simulations

Body-Centered Cubic (BCC) metals exhibit anomalous mechanical properties such as twinning/anti-twinning and tension/compression asymmetry, which are attributed to asymmetric behavior of screw dislocations. Unlike Face-Centered Cubic (FCC) metals, the Critical Resolved Shear Stress (CRSS) of BCC metals deviates from the Schmid law. Here, we use a mesoscale modeling approach called Phase Field Dislocation Dynamics (PFDD) to understand the critical factors that control this non-Schmid behavior and reproduce atomistic predictions of the CRSS (Peierls stress). All inputs to the PFDD model are obtained from Molecular Statics (MS) simulations. Multiple pathways for modeling the non-Schmid behavior are investigated by incorporating the representative dislocation properties into different energy terms in the PFDD model. One way to understand non-Schmid behavior is to incorporate stress components projected on inclined planes into the external energy term within PFDD. Alternatively, we propose that non-Schmid behavior can also be accounted for by considering the variation of the {110} unstable stacking fault energy and the dislocation core width as a function of the applied tensorial stress field. The CRSS predicted using PFDD modeling is in excellent agreement with MS predictions.

36 MATERIALS SCIENCE↗

Influence of Transition Metal Ion Contaminants on the Performance of Amine-Based Solid Sorbents in Direct Air Capture

Amine-functionalized solid sorbents are a class of sorbent materials proposed for direct air capture (DAC) of CO 2 , yet their long-term performance is susceptible to degradation under realistic operating conditions. Many amines are not thermodynamically stable in air, and amine sorbents oxidize while in use during DAC temperature swing adsorption processes. In this study, we investigate the role of transition metal ion contaminants, specifically Cu 2+ , Fe 2+ , and Ni 2+ , on the oxidative degradation of poly(ethylenimine) (PEI)-impregnated SBA-15 sorbents. By introducing metal ions via different modes mimicking both synthesis-related impurities and impurities derived from environmental exposure, we systematically evaluate sorbent stability after exposure to dry air at an elevated temperature. Thermogravimetric CO 2 uptake measurements reveal that even trace levels of Cu and Fe (as low as ∼4 ppm) can lead to measurable sorbent deactivation after oxidative aging, despite negligible loss in the performance of the control samples. In situ infrared, UV–vis, and X-ray photoelectron spectroscopies indicate that these metals catalyze radical-driven oxidation pathways, altering the chemical structure of the sorbent and accelerating degradation. Our findings underscore the need to account for trace metal contamination during DAC sorbent synthesis and deployment and highlight the importance of environmental contamination pathways.

CO2 capture↗

Tuning 3-D Nanomaterial Architectures Using Atomic Layer Deposition to Direct Solution Synthesis

The ability to synthesize nanoarchitected materials with tunable geometries provides a means to control their functional properties, with applications in biological, environmental, and energy fields. To this end, various bottom-up and top-down synthesis processes have been developed. However, many of these processes require prepatterning or etching steps, making them challenging to scale-up to complex, nonplanar substrates. Furthermore, the ability to integrate nanomaterials into hierarchical arrays with precise control of feature spacing and orientation remains a challenge. One approach to overcome these patterning challenges is the use of surface modification layers to guide the resulting geometry of nanomaterial architectures grown from the substrate. A powerful strategy to accomplish this is what we will refer to as “surface-directed assembly,” where the resulting geometric parameters (feature size, shape, orientation) are predetermined by the initial surface layer. In particular, the use of Atomic Layer Deposition (ALD) to form a surface layer, followed by solution-based growth processes, has the ability to synthesize architected structures with tunable geometries on complex, nonplanar surfaces. Over the past decade, we have reported a series of studies where surface-directed assembly is used to synthesize ZnO nanowires (NWs) on top of a variety of substrates. In this case, a thin film of ZnO is deposited onto the substrate using ALD, which can guide the NW diameter, spacing, and angular orientation with respect to the substrate by controlling epitaxial relationships. Furthermore, we have shown that by depositing a submonolayer overcoat of a secondary material (e.g., amorphous TiO 2 ), nucleation sites are partially blocked, which can further tune the spacing between nanowires while minimizing changes to their other geometric properties. This approach can be used to generate multilevel hierarchical structures, such as hyperbranched NW arrays with tunable control of each level of hierarchy using ALD. Finally, we have demonstrated that the tunable control of geometric parameters can be scaled-up to curved, nonplanar substrates. This highlights the power of ALD to conformally and uniformly deposit the seed layers on complex substrates with subnanometer precision. To complement these seeded hydrothermal approaches, we expanded this strategy to include conversion chemistry of the initial ALD seed layers. For example, by replacing ZnO with Al 2 O 3 as the seed layer without changing the hydrothermal growth conditions, Al–Zn layered-double hydroxide nanosheets can be formed instead of nanowires. In another example of conversion chemistry, a solution anion-exchange process was used to incorporate sulfur into ALD metal oxide films. In both of these conversion processes, the properties of the initial ALD film enabled tuning of the resulting nanostructure geometry. In this Account, we describe the use of ALD to guide the growth of diverse nanomaterial systems, with tunable control over their geometry and composition. We further show how these approaches can be used to tune functional properties for a range of applications, including superomniphobic surfaces, antibiofouling coatings, and photocatalysis. In conclusion, we conclude with an outlook on how the combination of ALD and solution synthesis can enable future directions in scalable nanomanufacturing to overcome the limitations of traditional top-down and bottom-up approaches.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Modeling of THz Lasers Based on Intersubband Transitions in Semiconductor Quantum Wells

In semiconductor quantum well structures, the intersubband energy separation can be adjusted to the terahertz (THz) frequency range by changing the well width and material combinations. The electronic and optical properties of these nanostructures can also be controlled by an applied dc electric field. These unique features lead to a large frequency tunability of the quantum well devices. In the on-going project of modeling of the THz lasers, we investigate the possibility of using optical pumping to generate THz radiation based on intersubband transitions in semiconductor quantum wells. We choose the optical pumping because in the electric current injection it is difficult to realize population inversion in the THz frequency range due to the small intersubband separation (4-40 meV). We considered both small conduction band offset (GaAs/AlGaAs) and large band offset (InGaAs/AlAsSb) quantum well structures. For GaAs/AlGaAs quantum wells, mid-infrared C02 lasers are used as pumping sources. For InGaAs/AlAsSb quantum wells, the resonant intersubband transitions can be excited by the near-infrared diode lasers. For three- and four-subband quantum wells, we solve the pumpfield-induced nonequilibrium distribution function for each subband of the quantum well system from a set of rate equations that include both intrasubband and intersubband relaxation processes. Taking into account the coherent interactions between pump and THz (signal) waves, we calculate the optical gain for the THz field. The gain arising from population inversion and stimulated Raman processes is calculated in a unified manner. A graph shows the calculated THz gain spectra for three-subband GaAs/AlGaAs quantum wells. We see that the coherent pump and signal wave interactions contribute significantly to the gain. The pump intensity dependence of the THz gain is also studied. The calculated results are shown. Because of the optical Stark effect and pump-induced population redistribution, the maximum THz gain saturates at larger pump intensities.

Liu, Ansheng↗

Closed-Loop Control of Active Nematic Flows

Stabilizing and shaping autonomous flows of active fluids is a fundamental challenge and a prerequisite for applications. We embed a light-responsive microtubule-based nematic in a proportional-integral control loop that adjusts the applied light intensity in response to real-time measurements of the spatially averaged flow speed. The self-regulating hardware-software-wetware system maintains a target flow speed against external or internal perturbations, including protein aging and aggregation, sample-to-sample variability, and temperature variation. Varying the controller’s gains reveals antagonistic roles between feedback and intrinsic processes, leading to nontrivial dynamics observed in fluctuation spectra. In particular, oscillations emerge from the interplay between the controller, motor binding kinetics, and active hydrodynamic relaxation. Accounting for the underlying binding timescale, our coarse-grained model and nematohydrodynamics simulations corroborate these observations. This work provides insight into the coupled dynamics of controlled active matter, laying the foundation for spatiotemporal patterning of active stress to generate and stabilize new dynamical configurations.

Active nematics↗

Simulating the Formation of Lunar Floor-Fracture Craters Using Elastoviscoplastic Relaxation

Lunar floor-fracture craters formed during the height of mare basalt emplacement. Due to a general temporal and spatial relation with the maria, these craters, numbering some 200, may be diagnostic of the thermal structure of the crust during this time. As the name suggests, these craters exhibit brittle failure, generally limited to the central floor region. That, and a shallower depth than fresh lunar craters, has led to two main theories as to their formation: laccolith emplacement under the crater and viscous relaxation. The implications of each model for the state of the Moon's crust during this time are quite different, so the viability of each model must be checked. Laccolith emplacement has been treated elsewhere. However, previous attempts to study the relaxation of the craters have assumed only a uniform, Newtonian viscous response of the near surface to the topographic driving forces, and simply postulated that the fractures resulted from tensile stresses associated with floor uplift. Here, we use a more sophisticated rheological model that includes not only non-Newtonian viscous behavior (i.e., the viscosity is stress-dependent), but also incorporates elastic behavior and a plastic component to the rheology to directly simulate the formation of the floor fractures. The results of our simulations show that while elastoviscoplastic relaxation is potentially viable for larger floor-fracture craters, it is not viable for craters with diameters < or = 60 km, the size of the majority of floor-fracture craters. We employ the finite element method, a numerical technique well suited for boundary-value problems, via the commercially available MARC software package. To test the viability of topographic relaxation, our goal is to prepare the simulations as to maximize the amount of relaxation. We take advantage of the natural axisymmetry of craters, simulating one radial plane. Initial shapes are based on data for fresh craters from Pike. To simplify implementation, a fourth order polynomial is used for the basin, while a third order inverse function is used for the rim. This form closely approximates the long-wavelength behavior of complex craters, while ignoring higher-frequency topography, save the rim. This approximation is appropriate because crater relaxation is strongly controlled by long-wavelength topography. Loading is accomplished assuming a uniform gravity field (1.62 m/s-square) and a uniform density of 2900 kg/cubic m. The initial stress state is set to be hydrostatic, with an additional pressure term to account for any overlying topography. This additional pressure term is tapered exponentially with depth. While the simulations quickly settle on a preferred stress state, and while the final solution is fairly insensitive to the choice of the e-folding depth of the taper, selecting an e-folding depth close to the diameter of the crater sets the initial stress state near the preferred state. We assume a diuranally averaged surface temperature of -20 C, and allow temperature to increase with depth at a rate of 50 K//km. Assuming a thermal conductivity of 2 W/in/K, this gradient translates to a heat flow of 100mW/square m, an extremely high value for the Moon. Temperature, of course, will not increase without bound. To maximize relaxation, we allow our temperature profile to increase linearly until it reaches the solidus (assumed to be 1200C) at a depth of 24.4 km, at which point it is kept constant. The presence of melt will drop the bulk viscosity; however, we have no rheological control for partial melts. Therefore, we make no attempt to simulate this situation. Elastoviscoplastic rheological model. In general, geologic materials can behave in three main ways: elastically, viscously (via solid-state creep), and brittly (plasticity is a continuum approach to simulate this phenomenon). We combine these three deformation mechanisms in an extended Maxwell solid, where the total strain can be broken down into a simple summation of the elastic, creep, and plastic strains. In relaxation phenomena in general, the system takes advantage of any means possible to eliminate deviatoric stresses by relaxing away the topography. Previous analyses have only modeled the viscous response. Comparatively, the elastic response in our model can augment the relaxation, to a point. This effect decreases as the elastic response becomes stiffer; indeed, in the limit of infinite elastic Young's modulus (and with no plasticity), the solution converges on the purely viscous solution. Igneous rocks common to the lunar near-surface have Young's modulii in the range of 10-100 GPa. To maximize relaxation, we use a Young's modulus of 10 GPa. (There is negligible sensitivity to the other elastic modulus, the Poisson's ratio; we use 0.25.) For the viscous response, we use a flow law for steady-state creep in thoroughly dried Columbia diabase, because the high plagioclase (about 70 vol%) and orthopyroxene (about 17 vol%) content is similar to the composition of the lunar highland crust as described by remote sensing and sample studies: noritic anorthosite. This flow law is highly non-Newtonian, i.e., the viscosity is highly stress dependent. That, and the variability with temperature, stands in strong contrast to previous examinations of lunar floor-fracture crater relaxation. To model discrete, brittle faulting, we assume "Byerlee's rule," a standard geodynamical technique. We implement this "rule" with an-angle of internal friction of about 40 deg, and a higher-than-normal cohesion of about 3.2 MPa (to approximate the breaking of unfractured rock). The actual behavior of geologic materials is more complex than in our rheological model, so the uncertainties in the plasticity do not represent the state-of-the-art error. Additional information is contained in the original.

Dombard, A. J.↗

A correlation between grain boundary character and deformation twin nucleation mechanism in coarse-grained high-Mn austenitic steel

Abstract In polycrystalline materials, grain boundaries are known to be a critical microstructural component controlling material’s mechanical properties, and their characters such as misorientation and crystallographic boundary planes would also influence the dislocation dynamics. Nevertheless, many of generally used mechanistic models for deformation twin nucleation in fcc metal do not take considerable care of the role of grain boundary characters. Here, we experimentally reveal that deformation twin nucleation occurs at an annealing twin (Σ3{111}) boundary in a high-Mn austenitic steel when dislocation pile-up at Σ3{111} boundary produced a local stress exceeding the twining stress, while no obvious local stress concentration was required at relatively high-energy grain boundaries such as Σ21 or Σ31. A periodic contrast reversal associated with a sequential stacking faults emission from Σ3{111} boundary was observed by in-situ transmission electron microscopy (TEM) deformation experiments, proving the successive layer-by-layer stacking fault emission was the deformation twin nucleation mechanism, different from the previously reported observations in the high-Mn steels. Since this is also true for the observed high Σ-value boundaries in this study, our observation demonstrates the practical importance of taking grain boundary characters into account to understand the deformation twin nucleation mechanism besides well-known factors such as stacking fault energy and grain size.

36 MATERIALS SCIENCE↗

Meso-scale modelling of compressive fracture in concrete with irregularly shaped aggregates

This paper presents a meso-scale modelling framework to investigate the fracture process in concrete subjected to uniaxial and biaxial compression accounting for its mesostructural characteristics. 3D mesostructure of concrete consisting of coarse aggregates, mortar and interfacial transition zone between them was developed using an in-house code based on the Voronoi tessellation and splining method, which enables to generate the realistic-look aggregates with controllable structural features such as content, location, size and shape. Based on the generated 3D mesostructure, the concrete damage plasticity approach was employed to simulate the compressive fracture behaviour of concrete in terms of crack morphology and stress-strain response against the shape parameters of aggregate. Results indicate that the shape of aggregate has a negligible effect on compressive strength of concrete, which is highly associated with the random location and size distribution of aggregate. The aggregate irregularity has a significant influence on crack initiation and growth of concrete.

36 MATERIALS SCIENCE↗

Time-dependent THMC properties and microstructural evolution of damaged rocks in excavation damage zone

Modeling coupled thermo-hydro-mechanical-chemical (THMC) processes in host rocks near high-level nuclear waste (HLW) repositories at various time scales is an extremely challenging task. The current study integrates experimental, theoretical, and numerical methods in assessing the evolution of excavation damage zone (EDZ) over time and its implication on the long-term migration of hazardous species. Argillite and rock salt and are the focus of this study. The first part of the report presents a novel time-dependent directional microcrack damage theory for generic brittle rocks. It features detailed statistical description of the microcracks within a damaged solid, permitting a direct upscaling of microscale processes such as crack growth kinetics, crack closure/opening, sliding friction to explain the macroscopic creep, nonlinear elasticity, shear dilation, and loading-unloading hysteresis. This provides a basic platform for describing the anisotropic mechanical and transport properties of damaged rocks during excavation and subsequent THMC loadings. The model is validated and numerically implemented to Finite Element (FE) package ABAQUS through the user-defined material (UMAT) interface and have demonstrated great potential in resolving the time-dependent and anisotropic evolution of damage in the EDZ. The second part of the report presents a multi-scale experimental effort in characterizing the thermal, hydraulic, and mechanical properties of Mancos shale and Avery Island salt. For the Mancos shale, triaxial compression tests are performed at different confining pressures and temperatures to probe its thermomechanical properties relevant to HLW repositories. The obtained stress-strain data are interpreted using the proposed directional damage theory. Post-test specimens are subjected to gas permeability tests to reveal the correlation between permeability and the degree of microcracking. At microscale, temperature-controlled nanoindentation tests are performed and found a linear correlation between fracture toughness and elastic modulus from 25°C to 300°C. For Avery Island salt, we have designed and manufactured a novel relative-humidity controlled uniaxial creep device. Long-term creep tests at low stresses (< 5 MPa) are performed at different levels of relative humidity (RH). Besides confirming the much higher creep rates as one would expect through extrapolating the high-stress creep data, the results reveal that the steady-state creep rate of rock salt is strongly dependent on the ambient RH, an aspect that is often neglected in the literature. Both behaviors can be attributed to the pressure-solution creep mechanism which dominates at low stress and high RH levels. The third part of the report explores a set of numerical strategies in modeling the THMC behavior of porous geomaterials. A fully implicit, monolithic FE solution that can flexibly interface with different material models and coupling mechanisms for THMC problems is developed and verified through the ABAQUS user-defined element (UEL) interface. The scheme is used to study the THM response of a hypothetical HLW storage site with reference to an existing in-situ heater test. Strategies for integrating the UEL and the microcrack UMAT are suggested. Another numerical endeavor of this study is to implement a higher-order asymptotic homogenization method to account for the heterogeneous porous structures. The same method is then extended to perform microstructure-informed thermo-mechanical modeling of generalized continua. The above outcomes of this project provide a strong thrust towards enhancing the fundamental understanding and modeling capability of the long-term evolution of host rocks in EDZ, thus helping achieve the design goal of 1-million-year isolation of high-level nuclear wastes.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Modeling the Effects of Microwave Heating on Densification in Chemical Vapor Infiltration

Microwave heating has great potential to accelerate the synthesis of ceramic matrix composites (CMCs) by a process called Chemical Vapor Infiltration (CVI). In CVI, reactive gases ingress a porous preform and undergo chemical transformation to deposit solid ceramic phase within the pores at high temperature thus, densifying the preform. However, the competing effects of chemical kinetics and gas transport are known to result in non-uniform densification as the outer surfaces of the preform experience faster depositions compared to the core. Achieving spatial temperature control plays a key role improving the densification quality. Microwave heating can potentially create temperature inversion such that the core of the preform is hotter than the outer surface and subsequently, lead to improved densification while keeping the manufacturing times and costs low. In the present work, a computational modeling strategy has been developed that accounts for the key physical phenomena responsible for densification of porous preforms using microwave heating. Specifically, a chemical kinetics model has been formulated for Silicon Carbide (SiC) deposition from MTS/H2 precursor. The model is implemented in a pore-resolved reactive transport solver, Quilt. The CVI simulations are performed for several conditions. Initially, parameterized temperature control is imposed to identify optimum conditions for good densification quality at a fraction of processing time. Further, simulations of microwave heating of porous SiC preforms are performed using OpenFOAM. A strategy to achieve and enhance temperature inversion is identified. The resulting temperature profiles are used in the pore-resolved densification simulations to analyze the densification behavior. It is observed that the temperature inversion achieved by microwave heating leads to significant improvements in densification quality and at the same time, keeps the manufacturing time low.

36 MATERIALS SCIENCE↗

From Coherence to Function: Exploring the Connection in Chemical Systems

The role of quantum mechanical coherences or coherent superposition states in excited state processes has received considerable attention in the last two decades largely due to advancements in ultrafast laser spectroscopy. These coherence effects hold promise for enhancing the efficiency and robustness of functionally relevant processes, even when confronted with energy disorder and environmental fluctuations. Understanding coherence deeply drives us to unravel mechanisms and dynamics controlled by order and synchronization at a quantum mechanical level, envisioning optical control of coherence to enhance functions or create new ones in molecular and material systems. In this frontier, the interplay between electronic and vibrational dynamics, specifically the influence of vibrations in directing electronic dynamics, has emerged as the leading principle. Here, two energetically disparate quantum degrees of freedom work in-sync to dictate the trajectory of an excited state reaction. Moreover, with the vibrational degree being directly related to the structural composition of molecular or material systems, new molecular designs could be inspired by tailoring certain structural elements. In the realm of chemical kinetics, our understanding of the dynamics of chemical transformations is underpinned by fundamental theories, such as transition state theory, activated rate theory, and Marcus theory. These theories elucidate reaction rates by considering the energy barriers that must be overcome for reactants to transform into products. Those barriers are surmounted by the stochastic nature of energy gap fluctuations within reacting systems, emphasizing that the reaction coordinate, the pathway from reactants to products, is not rigidly defined by a specific vibrational motion but encompasses a diverse array of molecular motions. While less is known about the involvement of specific intramolecular vibrational modes, their significance in certain cases cannot be overlooked. In this Account, we summarize key experimental findings that offer deeper insights into the complex electronic–vibrational trajectories encompassing excited states afforded from state-of-the-art ultrafast laser spectroscopy in three exemplary processes: photoinduced electron transfer, singlet–triplet intersystem crossing, and intramolecular vibrational energy flow in molecular systems. We delve into the rapid decoherence, or loss of phase and amplitude correlations, of vibrational coherences along promoter vibrations during subpicosecond intersystem crossing dynamics in a series of binuclear platinum complexes. This rapid decoherence illustrates the vibration-driven reactive pathways from the Franck–Condon state to the curve crossing region. We also explore the generation of new vibrational coherences induced by impulsive reaction dynamics rather than by the laser pulse in these systems, which sheds light on specific energy dissipation pathways and thereby on the progression of the reaction trajectory in the vicinity of the curve crossing on the product side. Another property of vibrational coherences, amplitude, reveals how energy can flow from one vibration to another in the electronic excited state of a terpyridine–molybdenum complex hosting a nonreactive dinitrogen substrate. In conclusion, a slight change in vibrational energy triggers a quasi-resonant interaction, leading to constructive wavepacket interference and ultimately intramolecular vibrational redistribution from a Franck–Condon active terpyridine vibration to a dinitrogen stretching vibration, energizing the dinitrogen bond.

Electrical energy↗