Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “crystal structure prediction”

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 343 records · Page 19

Compact Radiative Control Structures for Millimeter Astronomy

We have designed, fabricated, and tested compact radiative control structures, including antireflection coatings and resonant absorbers, for millimeter through submillimeter wave astronomy. The antireflection coatings consist of micromachined single crystal silicon dielectric sub-wavelength honeycombs. The effective dielectric constant of the structures is set by the honeycomb cell geometry. The resonant absorbers consist of pieces of solid single crystal silicon substrate and thin phosphorus implanted regions whose sheet resistance is tailored to maximize absorption by the structure. We present an implantation model that can be used to predict the ion energy and dose required for obtaining a target implant layer sheet resistance. A neutral density filter, a hybrid of a silicon dielectric honeycomb with an implanted region, has also been fabricated with this basic approach. These radiative control structures are scalable and compatible for use large focal plane detector arrays.

Brown, Ari D.↗

Topology optimization of 3D photonic crystals with complete bandgaps

The design of photonic crystals with complete bandgaps has recently received considerable research focus for numerous reasons. This work leverages well-known nonlinear programming techniques to alleviate the non-smoothness caused by degenerate eigenvalues such that topology optimization problems can be solved with the open-source IPOPT software. A fully-vectorial plane wave expansion technique is used with an iterative eigensolver to efficiently predict dispersion properties of candidate structures. Nonlinear programming is employed to solve the inverse problem of designing three-dimensional periodic structures that exhibit complete two-dimensional (2D) and three-dimensional (3D) photonic bandgaps. Mesh refinement is performed to alleviate the large computational burden of designing and analyzing photonic crystals, and a periodic density filter is implemented to impose a minimum feature size for manufacturability considerations.

97 MATHEMATICS AND COMPUTING↗

Orbital processing of high-quality CdTe compound semiconductors

CdZnTe crystals were grown in one-g and in micro-g for comparative analysis. The two micro-g crystals were grown in the Crystal Growth Furnace during the First United States Microgravity Laboratory mission (USML-1). The samples were analyzed for chemical homogeneity, structural perfection, and optoelectronic performance (infrared transmission). Fourier Transform Infrared (FTIR) transmission of both ground and flight materials showed that the infrared transmission was close to theoretical, 63% versus 66%, suggesting that the material was close to the stochiometric composition during both the ground and flight experiments. Infrared microscopy confirmed that the principal precipitates were Te and their size (1-10 microns) and density suggested that the primary flight and ground base samples experienced similar cooling rates. Macrosegregation was predicted, using scaling analysis, to be low even in one-g crystals and this was confirmed experimentally, with nearly diffusion controlled growth achieved even in the partial mixing regime on the ground. Radial segregation was monitored in the flight samples and was found to vary with fraction solidified, but was disturbed due to the asymmetric grvitational and thermal fields experienced by the flight samples. The flight samples, however, were found to be much higher in structural perfection than the ground samples produced in the same furnace under identical growth conditions except for the gravitational level. Rocking curve widths were found to be substantially reduced, from 20/35 (one-g) to 9/20 (micro-g) for the best regions of the crystals. The full width at half maximum (FWHM) of 9 arc seconds is as good as the best reported terrestrially for this material. The ground samples were found to have a fully developed mosaic structure consisting of subgrains, whereas the flight sample dislocations were discrete and no mosaic substructure was evident. The defect density was reduced from 50-100,000 (one-g) to 500-25000 EPD (micro-g). These results were confirmed using rocking curve analysis, synchrotron topography, and etch pit analysis. The low dislocation density is thought to have resulted from the near-absence of hydrostatic pressure which allowed the melt to solidify with minimum or no wall contact, resulting in very low stress being exerted on the crystal during growth or during post-solidification cooling.

Larson, David J., Jr.↗

Polyethylene Blends for Improved Oxygen Barrier: Processing-Dependent Microstructure and Gas Permeability

Here this work demonstrates a greater than expected enhancement of oxygen barrier properties in linear low-density polyethylene (LLDPE)-based materials by blending LLDPE with high-density polyethylene (HDPE). The films made by melt pressing the LLDPE/HDPE blends had a greater reduction in oxygen permeability coefficients (PO 2 ) than predicted using common permeability reduction models, i.e., the harmonic average model and zero-permeability nanofiller model. The reduction of PO 2 was attributed to the presence of spherulite crystal structures, as revealed by atomic force microscopy combined with infrared spectroscopy (AFM-IR). The LLDPE matrix exhibited significant spherulite formation even at a relatively low addition of HDPEs, which likely formed tortuous pathways for diffusing oxygen molecules. Transport results from melt-pressed films contrast with the results from films with similar compositions prepared by film blowing, which did not show barrier enhancement beyond expectation. AFM-IR revealed that the blown films lacked spherulite crystals likely due to stretching in the machine direction followed by rapid cooling. These findings demonstrate the role of processing in controlling microstructures and thus the oxygen barrier performance. This work offers the possibility of achieving easily recyclable LLDPE-based packaging materials by simple blending of polyethylenes with different crystalline content.

36 MATERIALS SCIENCE↗

Prospects for n -type doping of (Al x Ga 1– x ) 2 O 3 alloys

Here we systematically explore the properties of group-IV (C, Si, Ge, and Sn) and transition metal (Hf, Zr, and Ta) dopants substituting on the cation site in (Al x Ga 1– x ) 2 O 3 (AlGO) alloys using first-principles calculations with a hybrid functional. In Ga 2 O 3 , each of these dopants acts as a shallow donor. In Al 2 O 3 , they are deep defects characterized by the formation of either DX centers or positive- U (+/0) levels. Combining our calculations of dopant charge-state transition levels with information of the AlGO alloy band structure, we estimate the critical Al composition at which each dopant transitions from being a shallow to a deep donor. We identify Si to be the most efficient dopant to achieve n -type conductivity in high Al-content AlGO alloys, acting as a shallow donor over the entire predicted stability range for AlGO solid solution alloys.

36 MATERIALS SCIENCE↗

Identification of Small-Molecule Inhibitors of Fibroblast Growth Factor 23 Signaling via In Silico Hot Spot Prediction and Molecular Docking to α-Klotho

Fibroblast growth factor 23 (FGF23) is a therapeutic target for treating hereditary and acquired hypophosphatemic disorders, such as X-linked hypophosphatemic (XLH) rickets and tumor-induced osteomalacia (TIO), respectively. FGF23-induced hypophosphatemia is mediated by signaling through a ternary complex formed by FGF23, the FGF receptor (FGFR), and α-Klotho. Currently, disorders of excess FGF23 are treated with an FGF23-blocking antibody, burosumab. Small-molecule drugs that disrupt protein/protein interactions necessary for the ternary complex formation offer an alternative to disrupting FGF23 signaling. In this study, the FGF23:α-Klotho interface was targeted to identify small-molecule protein/protein interaction inhibitors since it was computationally predicted to have a large fraction of hot spots and two druggable residues on α-Klotho. We further identified Tyr433 on the KL1 domain of α-Klotho as a promising hot spot and α-Klotho as an appropriate drug-binding target at this interface. Subsequently, we performed in silico docking of ~5.5 million compounds from the ZINC database to the interface region of α-Klotho from the ternary crystal structure. Following docking, 24 and 20 compounds were in the final list based on the lowest binding free energies to α-Klotho and the largest number of contacts with Tyr433, respectively. Herein, five compounds were assessed experimentally by their FGF23-mediated extracellular signal-regulated kinase (ERK) activities in vitro, and two of these reduced activities significantly. Both these compounds were predicted to have favorable binding affinities to α-Klotho but not have a large number of contacts with the hot spot Tyr433. ZINC12409120 was found experimentally to disrupt FGF23:α-Klotho interaction to reduce FGF23-mediated ERK activities by 70% and have a half maximal inhibitory concentration (IC50) of 5.0 ± 0.23 μM. Molecular dynamics (MD) simulations of the ZINC12409120:α-Klotho complex starting from in silico docking poses reveal that the ligand exhibits contacts with residues on the KL1 domain, the KL1–KL2 linker, and the KL2 domain of α-Klotho simultaneously, thereby possibly disrupting the regular function of α-Klotho and impeding FGF23:α-Klotho interaction. ZINC12409120 is a candidate for lead optimization.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Candidate ferroelectrics via ab initio high-throughput screening of polar materials

Ferroelectrics are a class of polar and switchable functional materials with diverse applications, from microelectronics to energy conversion. Computational searches for new ferroelectric materials have been constrained by accurate prediction of the polarization and switchability with electric field, properties that, in principle, require a comparison with a nonpolar phase whose atomic-scale unit cell is continuously deformable from the polar ground state. For most polar materials, such a higher-symmetry nonpolar phase does not exist or is unknown. Here, we introduce a general high-throughput workflow that screens polar materials as potential ferroelectrics. We demonstrate our workflow on 1978 polar structures in the Materials Project database, for which we automatically generate a nonpolar reference structure using pseudosymmetries, and then compute the polarization difference and energy barrier between polar and nonpolar phases, comparing the predicted values to known ferroelectrics. Focusing on a subset of 182 potential ferroelectrics, we implement a systematic ranking strategy that prioritizes candidates with large polarization and small polar-nonpolar energy differences. To assess stability and synthesizability, we combine information including the computed formation energy above the convex hull, the Inorganic Crystal Structure Database id number, a previously reported machine learning-based synthesizability score, and ab initio phonon band structures. To distinguish between previously reported ferroelectrics, materials known for alternative applications, and lesser-known materials, we combine this ranking with a survey of the existing literature on these candidates through Google Scholar and Scopus databases, revealing ~130 promising materials uninvestigated as ferroelectric. Our workflow and large-scale high-throughput screening lays the groundwork for the discovery of novel ferroelectrics, revealing numerous candidates materials for future experimental and theoretical endeavors.

36 MATERIALS SCIENCE↗

An asynchronous parallel high-throughput model calibration framework for crystal plasticity finite element constitutive models

Crystal plasticity finite element model (CPFEM) is a powerful numerical simulation in the integrated computational materials engineering toolboxes that relates microstructures to homogenized materials properties and establishes the structure–property linkages in computational materials science. However, to establish the predictive capability, one needs to calibrate the underlying constitutive model, verify the solution and validate the model prediction against experimental data. Bayesian optimization (BO) has stood out as a gradient-free efficient global optimization algorithm that is capable of calibrating constitutive models for CPFEM. Here in this paper, we apply a recently developed asynchronous parallel constrained BO algorithm to calibrate phenomenological constitutive models for stainless steel 304 L, Tantalum, and Cantor high-entropy alloy.

304L stainless steel↗

Reversible structural transformations in supercooled liquid water from 135 to 245 K

Water has many anomalous properties compared to “simple” liquids, and these anomalies are typically enhanced in supercooled water.1-3 While numerous models have been proposed, including the liquid-liquid critical point,4,5 the singularity-free scenario,6 and the stability limit conjecture,1 a molecular-level understanding remains elusive. The main difficulty in determining which, if any, of these models is correct is the limited amount of data in the relevant temperature and pressure ranges. For water at ambient pressures, which is the focus of this work, data is largely missing from 160 – 232 K (“No Man’s Land”) due to rapid crystallization.2,3 Whether rapid crystallization is just an experimental obstacle, or a fundamental problem signaling the inability of water to thermally equilibrate prior to crystallization is also a major unanswered question.5,7 Here, we investigate the structural transformations of transiently-heated, supercooled water with nanosecond time resolution using infrared vibrational spectroscopy. The experiments demonstrate three key results. First, water’s structure relaxes from its initial configuration to a “steady-state” configuration prior to the onset of crystallization over a wide temperature range. Second, water’s steady-state structure can be reproduced by a linear combination of two, temperature-independent structures that correspond to a “high-temperature liquid” and a “low-temperature liquid.” Third, the observed structural changes are reversible over the full temperature range. Taken together, these results show that supercooled water can equilibrate prior to crystallization for temperatures from the homogeneous nucleation temperature, TH ~232 K,3 down to the glass transition temperature (Tg ~ 136 K). Second, the results provide support for the hypothesis that supercooled water can be described as a mixture of two, structurally-distinct, interconvertible liquids from 135 K to 245 K.5,8-18

Kringle, Loni M.↗

Ethane and methane at high pressures: Structure and stability

We have performed a combined experimental and theoretical study of ethane and methane at high pressures of up to 120 GPa at 300 K using x-ray diffraction and Raman spectroscopies and the USPEX ab initio evolutionary structural search algorithm, respectively. For ethane, we have determined the crystallization point, for room temperature, at 2.7 GPa and also the low pressure crystal structure (phase A). This crystal structure is orientationally disordered (plastic phase) and deviates from the known crystal structures for ethane at low temperatures. Moreover, a pressure induced phase transition has been identified, for the first time, at 13.6 GPa to a monoclinic phase B, the structure of which is solved based on good agreement with the experimental results and theoretical predictions. For methane, our x-ray diffraction measurements are in agreement with the previously reported high-pressure structures and equation of state (EOS). Finally, we have determined the EOSs of ethane and methane, which provides a solid basis for the discussion of their relative stability at high pressures.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

High-throughput multimodal exploration of a nanocrystalline Cu-Ag library

Sputter-deposited, nanocrystalline Cu-Ag thin films produced across a broad compositional and deposition-parameter space were evaluated to unravel the process-structure-property relationships important for creating hard, conductive electrical contacts and coatings. Combinatorial deposition involving pulsed direct current magnetron sputtering of elemental targets enabled swift examination of nearly the full range of alloy compositions and a relevant portion of deposition atomistics. Several high-throughput characterization modalities were employed to evaluate the chemistry, structure, and properties of the films. The resultant hardness, modulus, film density, crystal texture, and resistivity were analyzed in terms of key deposition characteristics (incident atom kinetic energy and incidence angle) predicted by binary-collision, kinematic Monte Carlo simulations. The study revealed improved hardness, parabolic resistivity dependence on composition, and compositional and process dependencies of film tarnishing. The results are discussed in the context of variations in microstructure and film density. Transmission electron microscopy and X-ray diffraction demonstrate several forms of compositional variation including solute segregation to grain boundaries as well as periodic, intragranular compositional modulations. Annealing of a Cu-rich alloy film exhibiting grain boundary segregation showed that this as-deposited, compositional variation is not stable above 100 °C. Finally, the Cu-Ag system is shown to have potential for hard, conductive, tarnish-resistant and room temperature-stable nanocrystalline thin films across the composition space.

36 MATERIALS SCIENCE↗

The gating mechanism of the large mechanosensitive channel MscL

The mechanosensitive channel of large conductance, MscL, is a ubiquitous membrane-embedded valve involved in turgor regulation in bacteria. The crystal structure of MscL from Mycobacterium tuberculosis provides a starting point for analysing molecular mechanisms of tension-dependent channel gating. Here we develop structural models in which a cytoplasmic gate is formed by a bundle of five amino-terminal helices (S1), previously unresolved in the crystal structure. When membrane tension is applied, the transmembrane barrel expands and pulls the gate apart through the S1-M1 linker. We tested these models by substituting cysteines for residues predicted to be near each other only in either the closed or open conformation. Our results demonstrate that S1 segments form the bundle when the channel is closed, and crosslinking between S1 segments prevents opening. S1 segments interact with M2 when the channel is open, and crosslinking of S1 to M2 impedes channel closing. Gating is affected by the length of the S1-M1 linker in a manner consistent with the model, revealing critical spatial relationships between the domains that transmit force from the lipid bilayer to the channel gate.

Non-NASA Center↗

Transferring predictions of formation energy across lattices of increasing size*

In this study, we show the transferability of graph convolutional neural network (GCNN) predictions of the formation energy of the nickel-platinum solid solution alloy across atomic structures of increasing sizes. The original dataset was generated with the large-scale atomic/molecular massively parallel simulator using the second nearest-neighbor modified embedded-atom method empirical interatomic potential. Geometry optimization was performed on the initially randomly generated face centered cubic crystal structures and the formation energy has been calculated at each step of the geometry optimization, with configurations spanning the whole compositional range. Using data from various steps of the geometry optimization, we first trained our open-source, scalable implementation of GCNN called HydraGNN on a lattice of 256 atoms, which accounts well for the short-range interactions. Using this data, we predicted the formation energy for lattices of 864 atoms and 2048 atoms, which resulted in lower-than-expected accuracy due to the long-range interactions present in these larger lattices. We accounted for the long-range interactions by including a small amount of training data representative for those two larger sizes, whereupon the predictions of HydraGNN scaled linearly with the size of the lattice. Therefore, our strategy ensured scalability while reducing significantly the computational cost of training on larger lattice sizes.

36 MATERIALS SCIENCE↗

Fabrication and Probabilistic Fracture Strength Prediction of High-Aspect-Ratio Single Crystal Silicon Carbide Microspecimens With Stress Concentration

Single crystal silicon carbide micro-sized tensile specimens were fabricated with deep reactive ion etching (DRIE) in order to investigate the effect of stress concentration on the room-temperature fracture strength. The fracture strength was defined as the level of stress at the highest stressed location in the structure at the instant of specimen rupture. Specimens with an elliptical hole, a circular hole, and without a hole (and hence with no stress concentration) were made. The average fracture strength of specimens with a higher stress concentration was larger than the average fracture strength of specimens with a lower stress concentration. Average strength of elliptical-hole, circular-hole, and without-hole specimens was 1.53, 1.26, and 0.66 GPa, respectively. Significant scatter in strength was observed with the Weibull modulus ranging between 2 and 6. No fractographic examination was performed but it was assumed that the strength controlling flaws originated from etching grooves along the specimen side-walls. The increase of observed fracture strength with increasing stress concentration was compared to predictions made with the Weibull stress-integral formulation by using the NASA CARES/Life code. In the analysis isotropic material and fracture behavior was assumed - hence it was not a completely rigorous analysis. However, even with these assumptions good correlation was achieved for the circular-hole specimen data when using the specimen data without stress concentration as a baseline. Strength was over predicted for the elliptical-hole specimen data. Significant specimen-to-specimen dimensional variation existed in the elliptical-hole specimens due to variations in the nickel mask used in the etching. To simulate the additional effect of the dimensional variability on the probabilistic strength response for the single crystal specimens the ANSYS Probabilistic Design System (PDS) was used with CARES/Life.

Nemeth, Noel N.↗

Establishing a process-structure-property-performance framework for SLS additive manufacturing through integrated multiscale modeling

This study presents a comprehensive suite of high-fidelity computational models that integrate multiscale and multiphysics simulations to capture the full Selective Laser Sintering (SLS) additive manufacturing process—from initial melting and solidification to mechanical response under external loads. Process simulations are linked with mechanical analysis through Representative Volume Elements (RVEs), establishing a process-structure–property-performance framework. The interaction between laser light and polyamide 12 (PA12) powder is modeled, accounting for laser characteristics and the optical, thermal, and geometrical properties of the powder. The heat source is incorporated into a heat transfer model, coupled with crystallization kinetics and densification models to predict material density and crystallinity. The porosity distribution from the densification model and crystallinity interpolated from experimental data are used to construct the RVEs. A multi-mechanism constitutive model is then calibrated using mechanical tests to predict the stress–strain response. Simulation results show good agreement with experimental data in terms of porosity, crystallinity, and mechanical performance when sufficient laser power (62 W or higher) is used. This research supports the inverse design of 3D-printed structures by introducing a high-fidelity framework that combines multiscale and multiphysics modeling with experimental calibration for predictive and performance-driven additive manufacturing.

SLS↗

Observation of an Unexpected n -Type Semiconducting Behavior in the New Ternary Zintl Phase Eu 3 InAs 3

The ternary arsenides Eu 3 InAs 3 and Sr 3 InAs 3 have been obtained by reactions of the elements in In flux at 1373 K. Structure elucidation by single-crystal X-ray diffraction reveals that Eu 3 InAs 3 and Sr 3 InAs 3 adopt the same orthorhombic structure (space group Pnma, Z = 4, Ca 3 AlAs 3 structure type) with unit cell parameters a = 12.9179(9) Å, b = 4.3990(3) Å, c = 13.9337(10) Å and a = 13.0218(11) Å, b = 4.4364(4) Å, c = 14.1339(12) Å, respectively. The structure consists of linear chains of corner-sharing InAs4 tetrahedra, [InAs 2 As 2/2 ] 6– , and Eu 2+ /Sr 2+ cations. Therefore, both Eu 3 InAs 3 and Sr 3 InAs 3 are valence-precise Zintl phases. As expected from the closed-shell electronic configurations, semiconducting behavior is confirmed by resistivity measurements on single crystals for both and by electronic band structure calculations for Sr 3 InAs 3 . The temperature dependence of resistivity and the computational work are in agreement that Eu 3 InAs 3 and Sr 3 InAs 3 are intrinsic semiconductors with narrow band gaps. Thermopower measurement on single-crystalline samples of Eu 3 InAs 3 shows that in the whole measured temperature range, from 300 to 700 K, the values for the Seebeck coefficient are negative. The observation of a negative Seebeck coefficient with very large absolute value (>400–500 μV K –1 at 700 K) is unexpected among the Zintl phases and suggestive that electrons are the majority charge carriers. Such a rare, n-type charge transport in an undoped compound such as Eu 3 InAs 3 , a material that has not been purposely optimized, could indicate native “defect” chemistry, and not extrinsic doping, as a reason for the unusual behavior. Here, a possible explanation involves a mixed-valent Eu 2+ /Eu 3+ state, which might be inferred from the measured effective paramagnetic moment of 7.2 μB per Eu atom, which is lower than the theoretically predicted value for free-ion moment of 7.9 μ B /Eu.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗