Search NASA⌕ Search

SEARCH · Search NASA

Results for “Orientation gradients”

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 145 records · Page 8

Reducing heat load density with asymmetric and inclined double-crystal monochromators: principles and requirements revisited

Asymmetric double-crystal monochromators (aDCMs) and inclined DCMs (iDCMs) can significantly expand the X-ray beam footprint and consequently reduce the heat load density and gradient. Based on rigorous dynamical theory calculations, the major principles and properties of aDCMs and iDCMs are presented to guide their design and development, particularly for fourth-generation synchrotrons. In addition to the large beam footprint, aDCMs have very large bandwidths (up to ∼10 eV) and angular acceptance, but the narrow angular acceptance of the second crystal requires precise control of the relative orientations and strains. Based on Fourier coupled-wave diffraction theory calculations, it is rigorously proved that the iDCM has almost the same properties as the conventional symmetric DCM, including the efficiency, angular acceptance, bandwidth, tuning energy range and sensitivity to misalignment. The exception is that, for the extremely inclined geometry that can achieve very large footprint expansion, the iDCM has (beneficially) a larger bandwidth and wider angular acceptance. Inclined diffraction has the `rho-kick effect' that can be cancelled by the second reflection of the iDCM (even with misalignment), except that inhomogeneous strains may cause non-uniform rho-kick angles. At present, fabrication/mounting-induced strains pose low risk since they can be controlled to <0.5 µrad over large areas. The only uncertain challenge is the thermally induced strains, yet it is estimated that these strains are naturally lowered by the large footprint and may be further mitigated by optimized cryogenic cooling to the 1–2 µrad level. Overall, aDCMs and iDCMs have more stringent requirements than normal DCMs, but they are feasible schemes in practice.

asymmetric monochromator↗

Simulation and Analysis on Reactor Pressure Vessel (RPV) subjected to Pressurized Thermal Shock (PTS) under SBLOCA scenario by using Cardinal to support the fracture mechanics analyses

The structural components that comprise nuclear reactors and their supporting structures are subjected to harsh operating environments that can challenge their integrity, especially after exposure for extended durations or under accident condition. As one of the most significant components of a Reactor, the Reactor Pressure Vessel (RPV) is exposed to an aggressive environment during the operation time (e.g. more than 40 years). Ageing degradation mechanisms (e.g. thermos-fatigue) could grow initial defects up to a critical size, increasing the susceptibility to failure in the RPV. The conventional methods are mostly based on simple crack and structure geometries. Very limited studies consider the real conditions of the RPV subjected to a thermal shock due to a Loss of Coolant Accident (LOCA). During a LOCA event, the most severe conditions take place when the emergency core cooling (ECC) water is injected inside the cold legs filled initially with hotter water and/or steam. The rapid cooling of the down-comer and the internal RPV surface followed probably by re-pressurization of the RPV causes large temperature gradients and variation of pressure which induces thermal-mechanical stresses. In order to develop the model for integrity assessment of a reactor pressure vessel (RPV) subjected to pressurized thermal shock (PTS), a multi-physics simulation, which includes the thermo-hydraulic, thermo-mechanical and fracture mechanics analyses is necessary. The multi-physics simulations are performed using Cardinal, a wrapping of the GPU-oriented spectral element Computational Fluid Dynamics (CFD) code NekRS and other multi-physics sub-modules within the MOOSE framework. Cardinal now fully supports MOOSE stochastic perturbations of NekRS models with varying boundary conditions, initial conditions, material properties, and any other quantity which is defined by a kernel (such as coefficients in a momentum source model). The implementation is designed in a flexible manner so that scalar values are sent from MOOSE into a user scratch space in NekRS, which can then be applied for any purpose within the NekRS case files (both on the host and device). When modeling PTS, several factors can impact the results significantly. In this report, the impacts of the geometry of the model, Reynolds number and buoyancy effect are investigated. Two geometry, i.e., a simplified model and a realistic RPV model, with both laminar and turbulent flow condition are adopted for the PTS simulation with and without buoyancy effect. The purpose of the investigation is to understand the impact of these factors on the prediction of temperature history of RPV. The accurate prediction on the temperature evolution, which will be exported to Grizzly code for further analyses on the progression of aging mechanisms and their effects on the integrity of RPV structures, is very crucial. Based on the understanding of these factors, a more sophisticated model is built to analysis the PTS under SBLOCA scenario. A literature survey is conducted to pick the SBLOCA scenario for the multi-physics simulation. The analysis helps to explain the form and the transformation of the cold plum when the ECC is activated under SBLOCA. This model can be can be applied to study the PTS effect for different RPV configurations. The results can help to assess structural component degradation for advanced reactors.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Effect of microstructure on fatigue crack propagation in additive manufactured nickel-based superalloy Haynes 282: an experiment and crystal plasticity study

Haynes 282 is a γ' precipitation-strengthened nickel-based superalloy known for its exceptional high-temperature creep resistance and excellent fabricability. Recent advancements in powder-bed fusion-based additive manufacturing (PBF-AM) have enabled the fabrication of Haynes 282 with accurate, site-specific control of the grain orientation and morphology at the microscale. This ability opens up new avenues for microstructure design, to improve the material’s fatigue crack resistance and service life. This paper investigates the fatigue crack growth behavior of hybrid microstructure Haynes 282 fabricated via PBF-AM. Previous experiments revealed a higher crack propagation rate in the coarse columnar-grained microstructural regions when compared against fine-grained areas. Here, a strain gradient crystal plasticity model was adapted to study the fracture-related mechanical fields at the crack tip in the two microstructures. The simulation results showed a consistent influence of grain structure and texture on crack propagation, as was seen in the experiment. The model analysis revealed higher crack propagation driving force along crack direction in the coarse-grained sharply textured microstructure and higher driving force for crack kinking in fine-grained more diffusely textured microstructure, which is ascribed to the combined effect of yield stress, hardening rate, texture and grain morphology. The presented modeling approach will facilitate the development of the AM-based accurate microstructure design by deepening the fundamental study in AM-specific microstructure-properties relations.

36 MATERIALS SCIENCE↗

Simulations of a hypersonic turbulent boundary layer over wavy surfaces

Here, we conduct large-eddy simulations of a Mach 5.84 cold wall turbulent boundary layer over one-dimensional wavy walls with varying amplitudes and wavelengths. Across all wall topologies, a series of alternating shock and expansion waves is shown to influence the entire boundary layer, and generate repeating wave patterns in the turbulent stresses, dispersive stresses, and turbulent kinetic energy budget. The series of alternating shocks and expansions imposes repeating adverse and favourable pressure gradients across the wavy wall, and at sufficient wall amplitude, triggers flow separation in the trough of the wave. Flow separation is demonstrated to influence the behaviour of wall pressure fluctuations over the wavy wall. In attached flows, the prominent frequencies are consistent with integral-scale boundary layer turbulence, whereas in separated flows, a two-decade frequency range is present, akin to two-dimensional shock–boundary layer interactions. Counter-rotating streamwise-oriented structures are observed on the windward side of the wave, which diminish over the wave crest. A conditional analysis demonstrates that these structures are present in the upstream boundary layer, and are amplified with increasing wall amplitude. An examination of the Görtler number and probability density function (PDF) of the fluctuating lateral wall shear stress demonstrates the strong correlation between a large Görtler number and growth of the PDF tail density, suggesting that the amplification of the counter-rotating streamwise-oriented structures are linked to centrifugal instabilities in regions of streamline concavity.

boundary layers↗

Wind Shear and Wind Veer Effects on Wind Turbines

This chapter highlights key contributions to the scientific literature on the sources of wind shear and wind veer in the atmospheric boundary layer, observations of shear and veer, and the effects of shear and veer on wind turbine power production, wind turbine wake evolution, and wind turbine loads. As wind turbines have grown larger, they encounter deeper and more complicated regions of the atmosphere. Over this height, profiles of wind speed shear and wind direction veer play a quantifiable role. Changes in the wind speed and wind direction across the vertical extent of a wind turbine rotor disk modify the inflow vector on the blades of the turbine, thereby affecting the magnitude and orientation of the lift and drag forces of the blade's airfoil. These changes can affect the power production and loads on large modern turbines, as well as the evolution of the wake that could affect a downwind turbine.

atmospheric stability↗

Application of electron backscatter diffraction techniques to quantify effects of aging on sub-grain and spatial heterogeneity in NMC cathodes

Identification and evaluation of structural heterogeneity in spent cathode materials is crucial to developing appropriate remediation strategies for novel recycling processes. Native heterogeneities may be exacerbated during the cell's operational lifetime, as sub-particle-scale variations induce anisotropic expansion and contraction upon cycling. Structural transformations resulting from repeated cycling and calendar aging predominantly occur at the secondary particle surface and at the grain boundaries (GBs) between primary particles. However, the diffusion and stress build up around and across GBs are poorly understood. In this study, electron backscatter diffraction (EBSD) is employed to track sub-grain lattice structure across a statistically relevant number of Li(Ni0.33Mn0.33Co0.33)O2 (NMC-111) particles. Specifically, differences in lattice misorientation – measured as the deviation from the grain's average orientation – are tracked as a function of position within the electrode (near current collector, middle, and near separator) and as a function of electrochemical cycling. Further, a novel method of structural analysis is developed, offering insight into sub-grain diffusion behavior by comparing lattice misorientation near the grain boundary versus in the grain bulk. The present results suggest that electrode-scale spatial heterogeneity in lattice structure is induced by initial manufacturing conditions, and that radial gradients in lattice misorientation evolve at the primary particle scale with repeated electrochemical cycling.

25 ENERGY STORAGE↗

Plasma confinement state classification via FPP relevant microwave diagnostics

We present a parsimonious and robust machine learning approach for identifying plasma confinement states in fusion power plants (FPPs) where reliable identification of the low-confinement and high-confinement regimes is critical for safe and efficient operation. Unlike research-oriented devices, FPPs must operate with a severely constrained set of diagnostics. To address this challenge, we demonstrate that a minimalist model, using only electron cyclotron emission (ECE) signals, can achieve accurate and reliable state classification. ECE provides electron temperature profiles without the engineering or survivability issues of in-vessel probes, making it a primary candidate for FPP-relevant diagnostics. Our framework employs ECE as input, extracts features using radial basis functions, and applies a gradient boosting classifier, achieving a test accuracy of 96% (correct predictions). Robustness analysis and feature importance analyzes confirm the approach’s reliability. These results demonstrate that state-of-the-art performance is attainable from a restricted diagnostic set, paving the way for minimalist yet resilient plasma control architectures for FPPs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Shear-Induced Conformations of Salt-Free Polyelectrolytes in Semidilute Solutions

Here, we present coarse-grained molecular dynamics simulations of salt-free polyelectrolyte chains in semidilute solutions under simple shear flow, with full hydrodynamic interactions and explicit dipolar solvent. At equilibrium, chain orientation statistics follow a pseudo-Voigt distribution, and the structural correlation length and chain end-to-end vector autocorrelation function exhibit scaling behavior consistent with theoretical predictions for polyelectrolytes. Under shear, chains transition from coiled to stretched states and the end-to-end vector autocorrelation function reveals oscillatory dynamics at high Weissenberg numbers. Analysis of the gyration tensor and shear strain distributions identified three distinct chain populations with directional alignment along and against the shear gradient. Compared with their neutral polymer counterparts, polyelectrolytes exhibit stronger shear thinning and enhanced chain alignment under the same Weissenberg number, which is attributed to electrostatic interactions and shear-induced counterion release. These findings provide molecular insight into the distinct flow response of charged polymers and have implications for tailoring the rheological properties of polyelectrolyte-based materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Optimizing Magnetic Force Microscopy Resolution and Sensitivity to Visualize Nanoscale Magnetic Domains

Magnetic force microscopy (MFM) enables mapping local magnetic fields across a sample surface with nanoscale resolution. To perform MFM, an atomic force microscopy (AFM) probe whose tip has been magnetized vertically (i.e., perpendicular to the probe cantilever) is oscillated at a fixed height above the sample surface. The resultant shifts in the oscillation phase or frequency, which are proportional to the magnitude and sign of the vertical magnetic force gradient at each pixel location, are then tracked and mapped. Although the spatial resolution and sensitivity of the technique increases with decreasing lift height above the surface, this seemingly straightforward path to improved MFM images is complicated by considerations such as minimizing topographical artifacts due to shorter range van der Waals forces, increasing the oscillation amplitude to further improve sensitivity, and the presence of surface contaminants (in particular water due to humidity under ambient conditions). In addition, due to the orientation of the probe's magnetic dipole moment, MFM is intrinsically more sensitive to samples with an out-of-plane magnetization vector. Here, high-resolution topographical and magnetic phase images of single and bicomponent nanomagnet artificial spin-ice (ASI) arrays obtained in an inert (argon) atmosphere glovebox with <0.1 ppm O 2 and H 2 O are reported. Further, optimization of lift height and drive amplitude for high resolution and sensitivity while simultaneously avoiding the introduction of topographical artifacts is discussed, and detection of the stray magnetic fields emanating from either end of the nanoscale bar magnets (~250 nm long and <100 nm wide) aligned in the plane of the ASI sample surface is shown. Likewise, using the example of a Ni-Mn-Ga magnetic shape memory alloy (MSMA), MFM is demonstrated in an inert atmosphere with magnetic phase sensitivity capable of resolving a series of adjacent magnetic domains each ~200 nm wide.

47 OTHER INSTRUMENTATION↗

Enhanced mechanical property evaluation using innovative data analytics capability

This report focuses on efforts to improve Multiphysics Object-Oriented Simulation Environ-ment (MOOSE) for mechanical property evaluation using data analytics. These efforts involveimprovements to the stochastic tools module (STM) for stochastic simulations of MOOSE mul-tiphysics model and the development of inverse optimization capabilities. The report gives anoverview of the STM and describes recent updates to its core capabilities and theory on itsreduced-order model (ROM) schemes. Examples are also provided showing the impact of theseupdates and exhibits the usefulness of ROMs. An overview of the gradient based inverse opti-mziation algorithms are given along with examples of their application to source identification.Inverse optimization will provide a new methodology in STM for fitting model parameters toexperimental data.

97 MATHEMATICS AND COMPUTING↗

Tuning Organic Semiconductor Packing and Morphology through Non-equilibrium Solution Processing (Final Report)

Organic semiconductors (OSCs) are a promising candidate to produce low-cost, large area, and flexible electronics. There has previously been successful development of organic field effect transistors (OFETs), photovoltaics (OPV), and bioelectronics using OSCs. Solution processing of these materials allows for the fabrication of large-area devices in the kinetic crystallization regime. The charge transport capabilities have been shown to depend on the morphology and molecular packing of the OSC within thin films. Our hypothesis is that solution processing conditions may significantly impact OSC morphology and as a result the charge transport. Therefore, this work has focused on obtaining a fundamental understanding of the morphology and molecular packing of OSCs for optimal device performance. Through our work, we have gained a better understanding of how these structural conditions were influenced by the solution processing conditions. Our studies have led to a more systematic understanding of the various parameters that impact OSC morphology. Solution processing leads to non-equilibrium films, thus allowing for the formation of diverse morphologies that are inaccessible by other fabrication methods. Previously, our group has focused on tuning the morphology of small-molecular OSCs [53, 55-57]. However, little work had been devoted to polymer OSCs. There was a lack of detailed studies characterizing the solution-state of polymer OSCs in terms of their conformation, degree of entanglement, polymer aggregation, and chain relaxation dynamics. The solution state properties are also likely to be influenced by the rigidity and molecular weight of the polymer OSCs investigated. Thus, we have focused our attention on gaining a better insight of polymer OSC films and solution processing methods. Through this proposal, our approach is to investigate the correlation between solution-state properties and the morphology of the resulting polymer OSC films. We worked three specific aims: investigate the effects of 1) polymer OSC solution-state properties on final film morphology; 2) molecular additives on the solution-state and final film properties; 3) controlled pre-aggregation in the solution-state in the final film morphology. More rigid and planar polymer backbones should promote interchain charge transport and more efficient interchain hopping between polymer chains. Through tailoring the polymer backbone, polymer sidechains, and molecular weight, we expected the altered solution-state properties to affect the final film morphology. In addition, reducing the entanglements and promoting chain alignment will likely prevent charge carrier trapping through conformation disorders. We thus studied different mechanisms, such as molecular additives, to reduce entanglements in solution. Devices fabricated from these solutions were expected to have improved charge transport abilities. In addition to tailoring the solution-state characteristics of the polymer OSCs, we investigated the effect of solution processing methods on the molecular packing and morphology of polymer OSC films. Two main solution processing techniques, e.g., spin-coating and solution shearing, were employed to fabricate OFETs. Spin coating was employed to prepare OSC films, which produce isotropic films. This technique creates several parameters to tune such as spin coating speed, acceleration, and time. Additionally, our research group developed the solution shearing method, which consists of the solution initially sandwiched between two plates. By sliding the top plate, the solution front is exposed and drying begins. This technique allows for the creation of aligned large crystalline domains. Solution shearing consists of different processing parameters which can affect the final morphology, such as shearing speed, substrate temperature and temperature gradient, distance between both plates, and the tilt angle of the top plate. Due to the complexity of the polymer systems investigated, numerous techniques were employed to characterize the polymer OSC solution state and films. All the materials were characterized using the DOE supported synchrotron X-ray scattering facilities at the Stanford Synchrotron Radiation Lightsource (SSRL). Using grazing incidence X-ray scattering (GIXS) and Near Edge X-ray Absorption Fine Structure (NEXAFS) techniques, the crystalline structure and molecular orientations of the thin films were measured. Optical absorption (UV-Vis) spectroscopy was employed to determine the aggregation state in solution and films of the polymer systems. Polarized UV-Vis also allowed for the determination of the relative degree of polymer chain alignment for solution sheared films. Additionally, various other techniques were used to investigate other properties within the film, such as atomic force microscopy (AFM) and solution rheology. Finally, the device performance is quantified through the fabrication and characterization of OFETs, which will highlight the effects of morphology on charge transport.

36 MATERIALS SCIENCE↗

Massively parallel axisymmetric fluid model for streamer discharges

A highly parallelizable fluid plasma simulation tool based upon the first-order drift-diffusion equations is discussed. Atmospheric pressure plasmas have densities and gradients that require small element sizes in order to accurately simulate the plasm resulting in computational meshes on the order of millions to tens of millions of elements for realistic size plasma reactors. To enable simulations of this nature, parallel computing is required and must be optimized for the particular problem. Here, a finite-volume, electrostatic drift-diffusion implementation for low-temperature plasma is discussed. The implementation is built upon the Message Passing Interface (MPI) library in C++ using Object Oriented Programming. The underlying numerical method is outlined in detail and benchmarked against simple streamer formation from other streamer codes. Electron densities, electric field, and propagation speeds are compared with the reference case and show good agreement. Convergence studies are also performed showing a minimal space step of approximately 4 μm required to reduce relative error to below 1% during early streamer simulation times and even finer space steps are required for longer times. Additionally, strong and weak scaling of the implementation are studied and demonstrate the excellent performance behavior of the implementation up to 100 million elements on 1024 processors. Lastly, different advection schemes are compared for the simple streamer problem to analyze the influence of numerical diffusion on the resulting quantities of interest.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A multiphase flow model of water droplets dielectrophoretic-induced air dehumidification phenomena

Air humidity in indoor spaces plays a critical role in human comfort and health. Dehumidification systems are used for building humidity controls, but they can take significant energy consumption, especially in geographic locations with high outdoor humidity and warm climates. Consequently, there is a growing demand for innovative dehumidification processes that consume minimal energy. Dielectrophoretic air dehumidification represents one such promising approach. However, it has not garnered significant attention due to the absence of engineering models and simulation tools capable of evaluating its performance and limitations at large-scale airflows. A new numerical multiphase CFD model, which is also experimentally validated, is developed in a customized Reacting Foam solver based on OpenFOAM® version 9. The newly developed model seeks to decrease substantial energy consumption and lower costs by leveraging the dielectrophoretic phenomenon to regulate moisture levels in the air. The solver integrates a hybrid Eulerian-Lagrangian framework to track the droplet's trajectory and growth rate while solving the continuum equations for the moist air. An electrospray produces electrically charged droplets, which grow during their in-flight trajectories as water vapor condenses onto their surfaces. The role of electrostatic forces in promoting vapor condensation within a high-gradient electrical field is investigated, and the dielectrophoretic vapor nucleation process on charged water droplets is discussed. The CFD model was validated against results from the literature and from proof-of-concept experiments conducted by the authors, which showed a 2 % air dehumidification with a single electrospray and airflow rate of 5 cubic feet per minute. The simulation results indicated that augmenting the number of electrically charged spray droplets increased the dehumidification of the air to 25 %. The initial mean droplet diameter, the orientation of the injector and relative humidity significantly influence the assessment of dehumidification. As a result, scaling up this approach to larger airflow volumes is identified as a potential future research direction.

42 ENGINEERING↗

Direct current magnetic Hall probe technique for measurement of field penetration in thin film superconductors for superconducting radio frequency resonators

Superconducting Radio Frequency (SRF) cavities used in particle accelerators are typically formed from or coated with superconducting materials. Currently, high purity niobium is the material of choice for SRF cavities that have been optimized to operate near their theoretical field limits. This brings about the need for significant R & D efforts to develop next generation superconducting materials that could outperform Nb and keep up with the demands of new accelerator facilities. To achieve high quality factors and accelerating gradients, the cavity material should be able to remain in the superconducting Meissner state under a high RF magnetic field without penetration of quantized magnetic vortices through the cavity wall. Therefore, the magnetic field at which vortices penetrate a superconductor is one of the key parameters of merit of SRF cavities. Techniques to measure the onset of magnetic field penetration on thin film samples need to be developed to mitigate the issues with the conventional magnetometry measurements that are strongly influenced by the film orientation and shape and edge effects. In this work, we report the development of an experimental setup to measure the field of full flux penetration through films and multi-layered superconductors. Our system combines a small superconducting solenoid that can generate a magnetic field of up to 500 mT at the sample surface and three Hall probes to detect the full flux penetration through the superconductor. We report this setup can be used to study alternative materials that could potentially outperform niobium, as well as superconductor–insulator–superconductor (SIS) multilayer coatings on niobium.

47 OTHER INSTRUMENTATION↗

Atomic-Scale Surface Studies of Bulk Metallic Glasses. Final Report

Bulk metallic glasses (BMGs) are of both scientific and technological interest because the absence of periodic atomic arrangements provides them with unique physical, chemical, and mechanical properties. Their high strength, superior elasticity, and an ability to be easily formed into virtually unlimited shapes with feature sizes from centimeters to Angstroms by blow molding and thermoplastic forming makes them an attractive choice for more and more practical applications and products. Due to their complex internal structure, however, experiments that yield insight into their exact atomic arrangements have been scarce. As a result, glass physics is one of the last remaining unexplored fields of materials science despite the scientific and technological importance of glasses in our daily lives, and the question how to characterize and control matter away from equilibrium, as glasses are, was listed as one of five Grand Challenges in a recent DoE report. The main reason for the slow progress in glass physics is the lack of experimental tools that enable access to atomic-scale structural and behavioral information for disordered materials. Such atomic-scale knowledge is mandatory to establish structure-property relationships that ultimately could allow to custom-design alloys featuring specific desired characteristics. With no such relationships available, theory development in glass remains basic and the few that exist are often untested. The aim of this research was to enable progress in our understanding of BMGs by developing a new approach that will allow a meaningful application of local surface science methods to specially prepared BMG samples to obtain a wealth of quantitative information on their atomic arrangements. Key was the availability of specially prepared samples whose surfaces feature large atomically flat terraces despite being entirely amorphous, which we have produced from a Pt 57.5 Cu 14.7 Ni 5.3 P 22.5 alloy (‘Pt-BMG’) both under ambient conditions as well as in ultrahigh vacuum using a unique setup that has been specially developed within this grant. Our approach starts with the in-situ preparation of oxide crystals that are terminated by large terraces, from which exact mirror images out of BMG will be produced using thermoplastic forming (TPF); for the research within this grant, we have successfully used (001)-oriented SrTiO 3 single crystal surfaces as well as (100)-, (110)-, and (111)-oriented single crystals made from LaAlO 3 . Since the resulting BMG replicas display all features of the original crystal with sub-Angstrom fidelity, thereby mimicking the original crystal’s termination by atomically flat terraces without being crystalline themselves, they are ideally suited for further investigation. The following atomic-scale local studies were then carried within this proposal: (i) high-resolution surface imaging and local spectroscopy using scanning probe microscopy, which showed disordered atom-like features and revealed changes in the gradient of the local surface potential on a 1-2 nm length scale; (ii) characterization of atomic-scale plastic flow with affected volumes as low as 1000 atoms, which showed local hardness near or above the theoretically predicted maximum and, once plastic deformation was initiated, homogeneous flow of the atoms involved; (iii) characterization of surface relaxation processes and the onset of crystallization induced by annealing, which showed that upon heating over the material’s glass transition temperature, the surface rearranges and relaxes towards a more stable, denser packed glass, which increases on-terrace surface roughness, while surface tension smoothens step edges; and (iv) studies that investigate the dependence of the material’s mechanical properties and structure on processing parameters, revealing that relaxed glasses get denser, harder, and more elastic. In combination, this information allows to combine structural models with mechanical properties and preparation history, thereby facilitating the development of preparation-structure-property relationships for metallic glasses. With the availability of such information, bulk metallic glasses can be further optimized to be used in more and more applications in industry.

36 MATERIALS SCIENCE↗

Measurements of Magnetic Field Penetration of Materials for Superconducting Radiofrequency Cavities

Superconducting Radio Frequency (SRF) cavities used in particle accelerators are typically formed from or coated with superconducting materials. Currently high purity niobium is the material of choice for SRF cavities which have been optimized to operate near their theoretical field limits. This brings about the need for significant R&D efforts to develop next generation superconducting materials which could outperform Nb and keep up with the demands of new accelerator facilities. To achieve high quality factors and accelerating gradients, the cavity material should be able to remain in the superconducting Meissner state under high RF magnetic field without penetration of quantized magnetic vortices through the cavity wall. Therefore, the magnetic field at which vortices penetrate in a supercon- ductor is one of the key parameters of merit of SRF cavities. Techniques to measure the onset of magnetic field penetration on thin film samples need to be developed to mitigate the issues with the conventional magnetometry measurements which are strongly influenced by the film orientation and shape and edge effects. The applied magnetic field is also needed to be parallel to the one side of the superconductor to resemble the magnetic field profile at the surface of the SRF cavities operating at fundamental accelerating mode. In this work we report the development of an experimental setup called Magnetic Field Penetration (MFP) magnetometer to measure the field of full flux penetration through bulk, thin films and multilayered superconductors. Our system combines a small superconducting solenoid which can generate the magnetic field up to 500 mT at the sample surface and three Hall probes to detect the full flux penetration through the planner superconductor with 2 inch diameter. This setup was used to study alternative materials which could potentially outperform niobium, as well as SIS multilayer coatings on niobium.

Senevirathne, Iresha↗

Mesoscale structural gradients in human tooth enamel

The outstanding mechanical and chemical properties of dental enamel emerge from its complex hierarchical architecture. An accurate, detailed multiscale model of the structure and composition of enamel is important for understanding lesion formation in tooth decay (dental caries), enamel development (amelogenesis) and associated pathologies (e.g., amelogenesis imperfecta or molar hypomineralization), and minimally invasive dentistry. Although features at length scales smaller than 100 nm (individual crystallites) and greater than 50 µm (multiple rods) are well understood, competing field of view and sampling considerations have hindered exploration of mesoscale features, i.e., at the level of single enamel rods and the interrod enamel (1 to 10 µm). Here, we combine synchrotron X-ray diffraction at submicrometer resolution, analysis of crystallite orientation distribution, and unsupervised machine learning to show that crystallographic parameters differ between rod head and rod tail/interrod enamel. This variation strongly suggests that crystallites in different microarchitectural domains also differ in their composition. Thus, we use a dilute linear model to predict the concentrations of minority ions in hydroxylapatite (Mg 2+ and CO 3 2- /Na + ) that plausibly explain the observed lattice parameter variations. While differences within samples are highly significant and of similar magnitude, absolute values and the sign of the effect for some crystallographic parameters show interindividual variation that warrants further investigation. In conclusion, by revealing additional complexity at the rod/interrod level of human enamel and leaving open the possibility of modulation across larger length scales, these results inform future investigations into mechanisms governing amelogenesis and introduce another feature to consider when modeling the mechanical and chemical performance of enamel.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Kohn-Sham Solver (KSSOLV) v2.0

KSSOLV is a MATLAB toolbox for solving Kohn-Sham density functional theory based electronic structure eigenvalue problems. It uses an object oriented features of MATLAB to represent atom, molecules, wavefunctions and Hamiltonians and their operations. It is designed to make it easier for users to prototype and test new algorithms for solving the Kohn-Sham problem. KSSOLV2.0 contains significant improvement over the original KSSOLV described in a paper published in ACM Transaction on Mathematical Software (attached). In addition to performing ground state calculation for small molecules, it can also perform geometry optimization for both molecules and solids. It uses standard pseudopotentials and implements local density approximation, generalized gradient approximation and hybrid functionals. Future releases will also include time-dependent DFT and post DFT calculations such as the GW quasi-particle energy calculation and Bethe-Salpeter equation solver for optical absorption.

Yang, Chao↗