Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Single crystal model”

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 757 records · Page 42
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Single-ion properties of the transverse-field Ising model material CoNb 2 O 6

We report CoNb 2 O 6 is one of the few materials that is known to approximate the one-dimensional transverse-field Ising model (1D-TFIM) near its quantum critical point. It has been inferred that Co 2+ acts as a pseudospin 1/2 with anisotropic exchange interactions that are largely Ising-like, enabling the realization of the TFIM. However, the behavior of CoNb 2 O 6 is known to diverge from the ideal TFIM under transverse magnetic fields that are far from the quantum critical point, requiring the consideration of additional anisotropic, bond-dependent (Kitaev-like) terms in the microscopic pseudospin-1/2 Hamiltonian. These terms are expected to be controlled in part by single-ion physics, namely the wave function for the pseudospin-1/2 angular momentum doublet. Here, we present the results of both inelastic neutron scattering measurements and electron paramagnetic resonance spectroscopy on CoNb 2 O 6 , which elucidate the single-ion physics of Co 2+ in CoNb 2 O 6 for the first time. We find that the system is well-described by an intermediate spin-orbit coupled Hamiltonian, and the ground state is a well-isolated Kramers doublet with an anisotropic g tensor. We provide the approximate wave functions for this doublet, which we expect will be useful in theoretical investigations of the anisotropic exchange interactions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

CHARMM-GUI Nanomaterial Modeler for Modeling and Simulation of Nanomaterial Systems

Molecular modeling and simulation are invaluable tools for nanoscience that predict mechanical, physicochemical, and thermodynamic properties of nanomaterials and provide molecular-level insight into underlying mechanisms. However, building nanomaterial-containing systems remains challenging due to the lack of reliable and integrated cyberinfrastructures. Here we present Nanomaterial Modeler in CHARMM-GUI, a web-based cyberinfrastructure that provides an automated process to generate various nanomaterial models, associated topologies, and configuration files to perform state-of- the-art molecular dynamics simulations using most simulation packages. The nanomaterial models are based on the interface force field, one of the most reliable force fields (FFs). The transferability of nanomaterial models among the simulation programs was assessed by single-point energy calculations, which yielded 0.01% relative absolute energy differences for various surface models and equilibrium nanoparticle shapes. Three widely used Lennard-Jones (LJ) cutoff methods are employed to evaluate the compatibility of nanomaterial models with respect to conventional biomolecular FFs: simple truncation at r = 12 Å (12 cutoff), force-based switching over 10 to 12 Å (10-12 fsw), and LJ particle mesh Ewald with no cutoff (LJPME). The FF parameters with these LJ cutoff methods are extensively validated by reproducing structural, interfacial, and mechanical properties. As a result, we find that the computed density and surface energies are in good agreement with reported experimental results, although the simulation results increase in the following order: 10-12 fsw <12 cutoff < LJPME. Nanomaterials in which LJ interactions are a major component show relatively higher deviations (up to 4% in density and 8% in surface energy differences) compared with the experiment. Nanomaterial Modeler's capability is also demonstrated by generating complex systems of nanomaterial-biomolecule and nanomaterial-polymer interfaces with a combination of existing CHARMM-GUI modules. We hope that Nanomaterial Modeler can be used to carry out innovative nanomaterial modeling and simulations to acquire insight into the structure, dynamics, and underlying mechanisms of complex nanomaterial-containing systems.

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

Elucidating Processes Controlling Arctic Atmospheric Aerosol Sources, Aging, and Mixing States (Final Report)

Atmospheric aerosols play critical roles in the Earth’s energy budget, directly by scattering or absorbing solar and terrestrial radiation and indirectly by serving as seeds (nuclei) for cloud droplet and ice crystal formation and by depositing on snow and ice surface, thereby changing the surface albedo. These effects are dependent on aerosol particle size and chemical composition and impact the hydrological cycle as well. This project provided single-particle size and chemical composition measurements across the entire annual cycle in the high Arctic and in the Alaskan Arctic during fall – winter, addressing the most significant gaps in Arctic aerosol observational data. These needs were based on recent rapid sea ice loss across the entire Arctic, as well as the major annual delays in sea ice freeze-up during fall in the Chukchi Sea and increased wintertime sea ice fracturing in the Beaufort Sea, both off the North Slope of Alaska. Two DOE Atmospheric Radiation Measurement (ARM) field campaigns were conducted for atmospheric aerosol sampling. The Aerosols during the Polar Utqiagvik Night (APUN – ‘snow on ground’ in Iñupiaq) ARM field campaign at Utqagivik, Alaska was conducted from Oct. 28 – Dec. 22, 2018. Aerosol sizing instrumentation and a single-particle mass spectrometer were successfully deployed for size-resolved number concentration measurements and measurements of individual particle size and chemical composition, respectively. These results show the influence of locally-produced sea spray aerosol, with high cloud-forming potential, due to delayed sea ice freeze-up in the fall. During the 2019‐2020 international Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC) expedition, daily atmospheric aerosol particles were collected aboard the German icebreaker Polarstern in the Central Arctic from Nov. 2019 – Oct. 2020. Sea salt aerosol and marine organics were observed year-round during MOSAiC with varying morphologies and sources. These findings are important because most Arctic models do not include a sea spray aerosol source, despite this source increasing with declining sea ice extent. In addition to collecting new samples and data, this project also conducted further analysis of previously collected single-particle chemical composition measurements within the North Slope of Alaska oil fields and at Utqiaġvik, AK, during Aug. – Sep. 2015 and 2016 field campaigns. This work resulted in the discovery of chemical reactions of oil field combustion emissions occurring within fog droplets across the North Slope of Alaska and forming secondary aerosol, showing the impact of Arctic oil field emissions beyond black carbon aerosol and greenhouse gases. In addition, the distribution of chemical species across the aerosol population within the oil fields was quantified, using these data and a previously development framework. We also presented the first ambient evidence of the collision of two atmospheric particles resulting in formation of an organic-coated ammonium sulfate particle of marine origin, which has implications for cloud formation with declining sea ice extent. Overall, this project has elucidated connections between seawater biogeochemistry, resource extraction activities, atmospheric composition, clouds, and the energy budget of the Arctic region. The results of this project are expected to improve weather and sea ice forecasting for security and development in the Arctic and beyond.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Toward Understanding the Simulated Phase Partitioning of Arctic Single-Layer Mixed-Phase Clouds in E3SM

21 Mixed-phase clouds are widespread in the Arctic and have substantial impacts on the 22 surface energy budget, sea ice extent, and regional climate. Significant changes are found in the 23 modeled phase partitioning of Arctic mixed-phase clouds in the newly developed U.S. 24 Department of Energy (DOE) Energy Exascale Earth System Model (E3SM) Atmosphere Model 25 version 1 (EAMv1) compared to its predecessor (CAM5/EAMv0). The mixed-phase clouds 26 simulated by EAMv1 is overly dominated by super-cool liquid which is often underestimated by 27 CAM5. Compared to CAM5, EAMv1 has made several important changes in its cloud physical 28 parameterizations, including using CLUBB to unify shallow convection, PBL turbulence, and 29 cloud macrophysics, as well as updating the MG two-moment cloud microphysics scheme and 30 using the CNT heterogeneous ice nucleation scheme to replace the Meyers scheme for mixed- 31 phase clouds. To understand which of the newly-implemented physical parameterizations is 32 responsible for the change of modeled mixed-phase cloud phase partitioning in EAMv1, 33 sensitivity tests are conducted by utilizing the short-term hindcast approach and results are 34 compared to the observations from the ARM Mixed-Phase Arctic Cloud Experiment (M-PACE) 35 field campaign. 36 It is found that although each scheme contributes to the decrease of cloud ice water 37 content (IWC) and increase of cloud liquid water content (LWC), only the combined effect of 38 CLUBB and CNT can explain the significant underestimation of IWC. By examining the budgets 39 of cloud microphysical processes, CLUBB is found to play a dominant role due to the fact that 40 CLUBB does not diagnoses cloud ice phase water in its assumed joint PDF method and thus may 41 not provide sufficient cloud ice to initialize ice phase microphysical processes such as the growth 42 of ice crystals through vapor deposition. Although a more effective heterogeneous ice nucleation 43 may alleviate the underestimation of IWC, the ice phase process in CLUBB requires urgent 44 attention in the future model development.

Zhang, Meng↗
DOE OSTI2025

Composition dependence of atomic order in strain-relaxed, metastable GeSn alloys

Extended x-ray absorption fine structure (EXAFS) measurements of single-crystal Ge/GeSn radial heterostructure nanowires are used to examine the effects of composition on both short-range order (SRO) and longer-range disorder in GeSn alloys. GeSn has prompted significant interest because it can achieve a direct band gap for sufficient Sn concentrations beyond the equilibrium solid solubility limit in an all-group IV system. Short-range order in this material is particularly interesting as it has been predicted to affect the band gap independent of average composition or strain effects. By independently controlling the Sn composition and GeSn thickness during chemical vapor deposition of misfitting GeSn shells around ultrathin, elastically compliant, Ge core nanowires, the elastic misfit strain in the GeSn is minimized for Sn compositions over the studied range ≈Ge 0.96 Sn 0.04 to Ge 0.88 Sn 0.12 . The degree of SRO was found to decrease with increasing Sn composition. Additionally, damping of the EXAFS signal was observed as the Sn content increased, particularly for increasingly distant atomic shells about the absorbing atom, even for scattering paths not involving Sn atoms. This result is quantified as an increase in the mean-squared relative displacement parameters of the shells. These measurements reveal the accommodation of local strain due to the presence of the highly size-mismatched Sn atoms in the Ge diamond cubic lattice (≈14%), which may have effects on the band structure of the material in addition to the influence of short-range atomic order. Comparison among the nanowire samples allows for calculation of the topological rigidity parameter, a ∗∗ , for the first-neighbor bond lengths. Furthermore, these exhibit chemically distinct values for Ge-Ge, Ge-Sn, and Sn-Sn, and they are consistent with the value a ∗∗ = 0.75 ± 0.07 confirming the general applicability of the model to alloys with both large amounts of natural misfit strain and the potential for short-range order.

Crystal structure Semiconductors Transmission elec↗
DOE OSTI2021

Ground state in proximity to a possible Kitaev spin liquid: The undistorted honeycomb iridate Na x IrO 3 (0.60 ≤ x ≤ 0.80)

Here, we report the results of our study of a recently synthesized honeycomb iridate Na x IrO 3 (0.60 ≤ x ≤ 0.80). Single-crystal Na x IrO 3 adopts a honeycomb lattice noticeably without distortions and stacking disorder inherently existent in its sister compound Na 2 IrO 3 . The oxidation state of the Ir ion is a mixed valence state resulting from a majority Ir 5+ (5d 4 ) ion and a minority Ir 6+ (5d 3 ) ion. Na x IrO 3 is a Mott insulator likely with a predominant pseudospin =1 state. It exhibits an effective moment of 1.1 μB/Ir and a Curie-Weiss temperature of –19 K but with no discernible long-range order above 1 K. The physical behavior below 1 K features two prominent anomalies at T h = 0.9 K and T l = 0.12 K in both the heat capacity and AC magnetic susceptibility. Intermediate between T h and T l lies a pronounced temperature linearity of the heat capacity with a large slope of 77 mJ /mole K 2 , a feature expected for highly correlated metals but not at all for insulators. These results along with a comparison drawn with the honeycomb lattices Na 2 IrO 3 and (Na 0.2 Li 0.8 ) 2 IrO 3 point to an exotic ground state in proximity to a possible Kitaev spin liquid.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Metadynamics simulations reveal mechanisms of Na + and Ca 2+ transport in two open states of the channelrhodopsin chimera, C1C2

Cation conducting channelrhodopsins (ChRs) are a popular tool used in optogenetics to control the activity of excitable cells and tissues using light. ChRs with altered ion selectivity are in high demand for use in different cell types and for other specialized applications. However, a detailed mechanism of ion permeation in ChRs is not fully resolved. Here, we use complementary experimental and computational methods to uncover the mechanisms of cation transport and valence selectivity through the channelrhodopsin chimera, C1C2, in the high- and low-conducting open states. Electrophysiology measurements identified a single-residue substitution within the central gate, N297D, that increased Ca 2+ permeability vs. Na + by nearly two-fold at peak current, but less so at stationary current. We then developed molecular models of dimeric wild-type C1C2 and N297D mutant channels in both open states and calculated the PMF profiles for Na + and Ca 2+ permeation through each protein using well-tempered/multiple-walker metadynamics. Results of these studies agree well with experimental measurements and demonstrate that the pore entrance on the extracellular side differs from original predictions and is actually located in a gap between helices I and II. Cation transport occurs via a relay mechanism where cations are passed between flexible carboxylate sidechains lining the full length of the pore by sidechain swinging, like a monkey swinging on vines. In the mutant channel, residue D297 enhances Ca 2+ permeability by mediating the handoff between the central and cytosolic binding sites via direct coordination and sidechain swinging. We also found that altered cation binding affinities at both the extracellular entrance and central binding sites underly the distinct transport properties of the low-conducting open state. This work significantly advances our understanding of ion selectivity and permeation in cation channelrhodopsins and provides the insights needed for successful development of new ion-selective optogenetic tools.

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

Oxygen Adsorption, Subsurface Oxygen Layer Formation and Reaction with Hydrogen on Surfaces of a Pt–Rh Alloy Nanocrystal

The oxygen adsorption and its catalytic reaction with hydrogen on Pt–Rh single crystals were studied at the nanoscale by Field Emission Microscopy (FEM) and Field Ion Microscopy (FIM) techniques at 700 K. Both FEM and FIM use samples prepared as sharp tips, apexes of which mimic a single nanoparticle of catalyst considering their similar size and morphology. Oxygen adsorption on Pt-17.4 at.%Rh samples leads to the formation of subsurface oxygen, which is manifested in the field emission (FE) patterns: for O2 exposure of ~3 Langmuir (L), {113} planes appear bright in the emission pattern, while for higher oxygen doses, i.e. 84 L, the bright regions correspond to the high index planes between the {012} and {011} planes. Formation of subsurface oxygen is probably accompanied by a surface reconstruction of the nanocrystal. The subsurface oxygen can be effectively reacted off by subsequent exposure of the sample to hydrogen gas at 700 K. The hydrogenation reaction was observed as a sudden, eruptive change of the brightness seen on the FE pattern. This reaction resulted in the recovery of the initial field emission pattern characteristic of a clean tip, with {012} facets being the most visible. It was shown that the oxygen accumulation-reduction process is completely reversible. The obtained results indicate that the presence of subsurface species must be considered in the description of reactive processes on Pt–Rh catalysts.

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

Heterogeneous fatigue damage in a nickel-based single-crystal superalloy unraveled using correlative 3D X-ray technology

Nickel-based single-crystal (Ni-SX) superalloys under cyclic stress are susceptible to cracking at stress-concentration sites, eventually leading to low-cycle fatigue (LCF) failure. LCF cracks typically originate from intrinsic defects (e.g., voids and carbides) within solidified dendrites. However, systematic quantitative experimental analyses of defect-mediated local damage remain limited. To thoroughly understand the microscopic origins and evolution of LCF damage, correlated 3D mapping of dendrites across various regions is essential. Here, in this study, macroscale micro-computed tomography (μ-CT) was initially used to capture internal interdendritic secondary cracks within bulk DD413 superalloy after LCF testing at 760 °C. Subsequently, a multimodal methodology combining synchrotron 3D microdiffraction (3D-μXRD), high-resolution μ-CT, and electron microscopy was established. This approach allowed precise localization of internal damage zones near interdendritic secondary cracks and detailed mapping of the 3D correlated distributions of dendrites, defects, and residual stress/strain fields within these zones at submicron spatial resolution. Finally, the same approach was applied to specimens subjected to interrupted loading at approximately 40 % of the fatigue life to uncover the early damage states of dendrites. The dendrite cores (DCs) and interdendritic regions (IDs) exhibit microscale heterogeneous mechanical responses: nearly defect-free DCs accumulate local irreversible slip along specific slip systems to generate slip bands, while the IDs containing various defects accommodate local microplasticity through the activation of multiple slip systems around these defects. The local tensile stress near defects in the IDs exceeds that in the DC slip band regions by more than threefold, leading to the generation of local damage zones within the IDs. Chain-like defect distributions facilitate the interconnection of these local zones into a continuous damage region, further elevating the overall tensile stress in the IDs. Additionally, geometrically necessary dislocations alone are insufficient as indicators of LCF damage; both the internal stress state and its magnitude must be considered. These experimental results provide critical data and insights for the development of multi-physics fatigue models.

Localized deformation↗
DOE OSTI2020

Ice Nucleating Particles, Aerosols and Clouds over the Higher Latitude Southern Ocean

Improved understanding of aerosol-cloud-precipitation interactions is a key goal of the DOE-ASR program. Ice Nucleating Particles (INPs) are rare aerosol particles that trigger ice formation in clouds. By doing so they may initiate precipitation in clouds, and this impacts cloud lifetime. Nowhere may the role of INPs be manifested so dramatically as in the Southern Ocean (SO). This vast region was the focus of the DOE-ARM MARCUS (Measurement of Aerosols, Radiation and CloUds over the Southern Oceans) and MICRE (Macquarie Island Cloud and Radiation Experiment) campaigns. In the SO, a deficit of INPs has been a leading hypothesis for why liquid clouds persist, which underlays the bias of global climate models in predicting excess shortwave radiation reaching the ocean surface between 55°S and Antarctica. In this study, we used the 6-month acquisition of INP data on four MARCUS cruises from Hobart, Tasmania, to Antarctica, and the single-location full annual INP cycle observed in MICRE, to gain a holistic, quantitative picture of the INP number concentrations as a function of temperature and their variability over SO latitudes from 43°S to Antarctica. Augmenting initial processing of filter collections of aerosol particles during MARCUS and MICRE, additional freezing studies following thermal treatment (to removes INPs associated with microbes) and peroxide treatment (to remove all organics), as well as ionic chemistry and total organic carbon analysis, were applied to samples to improve time and space resolution, and differentiate marine from terrestrial contributions. We also performed Next Generation Sequencing to determine that bacterial composition can tag periods of enhanced and degraded marine organic INPs and the relative lack of occurrence of land-sourced bio-particles. Through these and other analyses, INP data were categorized into representative source types (e.g., marine versus terrestrially-influenced). The products of our analyses are being used to parameterize ice formation for use in numerical modeling studies to determine if, and when, ocean-derived INPs control the microphysical composition and radiative balance of SO clouds in a manner not presently captured by global models. The specific objectives and approaches proposed were largely accomplished, including: 1) Additional and value-added chemical and biological analyses of archived samples of aerosols were completed to categorize and quantify INP types and concentrations over the SO. 2) Analyses were completed to place INP data in meteorological and aerosol context. 3) Sea spray aerosol INPs were demonstrated to dominate the SO region in all but episodic events. 4) New parameterizations for INP sources relevant to the SO were constructed. 5) Collaborative modeling studies demonstrating the crucial role of INPs in determining cloud radiative and precipitation properties were conducted with MARCUS and MICRE partners. Collaborative publication submissions in this regard ensued within the two year study. This work will advance the science of interactions of aerosols, clouds and precipitation, to improve their representation in regional and global climate models. Results have informed representation of primary ice crystal nucleation, give inference to the role of secondary processes, and improve investigations of aerosol influences on clouds via ice nucleation over SO high latitudes and similar vast ocean regions. Application will improve representation of such interactions for clouds in both regional and global climate models. The data base and analyses methods applied will continue to serve research studies in the future.

58 GEOSCIENCES↗
DOE OSTI2024

Evidence of Dirac Quantum Spin Liquid in YbZn 2⁢ GaO 5

The emergence of a quantum spin liquid (QSL), a state of matter that can result when electron spins are highly correlated but do not become ordered, has been the subject of a considerable body of research in condensed matter physics [1,2]. Spin liquid states have been proposed as hosts for high-temperature superconductivity [3] and can host topological properties with potential applications in quantum information science [4]. The excitations of most quantum spin liquids are not conventional spin waves but rather quasiparticles known as spinons, whose existence is well established experimentally only in one-dimensional systems; the unambiguous experimental realization of QSL behavior in higher dimensions remains challenging. Here, in this work, we investigate the novel compound $YbZn_{2}GaO_{5}$, which hosts an ideal triangular lattice of effective spin-1/2 moments with no detectable inherent chemical disorder. Thermodynamic and inelastic neutron scattering measurements performed on high-quality single crystal samples of $YbZn_{2}GaO_{5}$ exclude the possibility of long-range magnetic ordering down to 0.06 K, demonstrate a quadratic power law for the specific heat and reveal a continuum of magnetic excitations in parts of the Brillouin zone. Both low-temperature thermodynamics and inelastic neutron scattering spectra suggest that $YbZn_{2}GaO_{5}$ is a U(1) Dirac QSL with spinon excitations concentrated at certain points in the Brillouin zone. We advanced these results by performing additional specific heat measurements under finite fields, further confirming the theoretical expectations for a Dirac QSL on the triangular lattice of $YbZn_{2}GaO_{5}$.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

New Mechanism for Yield Point Phenomena

Yield point phenomena (YPP) are widely attributed to discrete dislocation locking by solute atmospheres. An alternate YPP mechanism was recently suggested by simulations of Ta single crystals without any influence of solutes or discrete dislocations. The general meso-scale (GM) simulations consist of crystal plasticity (CP) plus accounting for internal stresses of geometrically necessary dislocation content. GM predicted the YPP while CP did not, suggesting a novel internal stress mechanism. The predicted YPP varied with crystal orientation and boundary conditions, contrary to expectations for a solute mechanism. The internal stress mechanism was probed by experimentally deforming oligocrystal Ta samples and comparing the results with independent GM simulations. Strain distributions of the experiments were observed with high-resolution digital image correlation. A YPP stress–strain response occurred in the 0–2% strain range in agreement with GM predictions. Shear bands appeared concurrent with the YPP stress–strain perturbation in agreement with GM predictions. At higher strains, the shear bands grew at progressively slower rates in agreement with GM predictions. It was concluded that the internal stress mechanism can account for the existence of YPP in a wide variety of materials including ones where interstitial-dislocation interactions and dislocation transient avalanches are improbable. Furthermore, the internal stress mechanism is a CP analog of various micro-scale mechanisms of discrete dislocations such as pile-up or bow-out. It may operate concurrently with strain aging, or either mechanism may operate alone. A suggestion was made for a future experiment to answer this question.

42 ENGINEERING↗
DOE OSTI2025

Interplay Between Metastability and Mechanically Induced Structural Instability in CsPbBr 3 Photovoltaic Perovskite

An atomic-level understanding of the underlying structural metastability is still absent in halide perovskite photovoltaic systems. Focusing on the model material CsPbBr 3 , the impact of mechanically induced atomic structure alteration is elucidated through structural modeling and X-ray diffraction measurements. Sudden Cs–Br bond breaking drives the system metastability, where the first-order transition arises from cation-halide bond strain relief by severing the corner connectivity of the PbBr 6 units. The pressure–volume and entropy terms govern the Gibbs free-energy landscape in halide perovskites. In conclusion, metastability is revealed as a critical factor limiting the performance of lead halide perovskites under extreme but natural environments.

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

OpenAWSEM with Open3SPN2: A fast, flexible, and accessible framework for large-scale coarse-grained biomolecular simulations

We present OpenAWSEM and Open3SPN2, new cross-compatible implementations of coarse-grained models for protein (AWSEM) and DNA (3SPN2) molecular dynamics simulations within the OpenMM framework. These new implementations retain the chemical accuracy and intrinsic efficiency of the original models while adding GPU acceleration and the ease of forcefield modification provided by OpenMM’s Custom Forces software framework. By utilizing GPUs, we achieve around a 30-fold speedup in protein and protein-DNA simulations over the existing LAMMPS-based implementations running on a single CPU core. We showcase the benefits of OpenMM’s Custom Forces framework by devising and implementing two new potentials that allow us to address important aspects of protein folding and structure prediction and by testing the ability of the combined OpenAWSEM and Open3SPN2 to model protein-DNA binding. The first potential is used to describe the changes in effective interactions that occur as a protein becomes partially buried in a membrane. We also introduced an interaction to describe proteins with multiple disulfide bonds. Using simple pairwise disulfide bonding terms results in unphysical clustering of cysteine residues, posing a problem when simulating the folding of proteins with many cysteines. We now can computationally reproduce Anfinsen’s early Nobel prize winning experiments by using OpenMM’s Custom Forces framework to introduce a multi-body disulfide bonding term that prevents unphysical clustering. Our protein-DNA simulations show that the binding landscape is funneled towards structures that are quite similar to those found using experiments. In summary, this paper provides a simulation tool for the molecular biophysics community that is both easy to use and sufficiently efficient to simulate large proteins and large protein-DNA systems that are central to many cellular processes. These codes should facilitate the interplay between molecular simulations and cellular studies, which have been hampered by the large mismatch between the time and length scales accessible to molecular simulations and those relevant to cell biology.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Diffusion of Hydrophilic to Hydrophobic Forms of Nile Red in Aqueous C 12 EO 10 Gels by Variable Area Fluorescence Correlation Spectroscopy

In this work, solute diffusion within lyotropic liquid crystal gels prepared from a series of water and decaethylene glycol monododecyl ether (C 12 EO 10 ) mixtures was explored by variable area fluorescence correlation spectroscopy. Aqueous C 12 EO 10 gels were prepared in concentrations ranging from 55:45 to 70:30 wt% of surfactant and water. Small angle X-ray scattering revealed that these gels comprise hexagonal mesophases of cylindrical micelles. Micelle spacing was found to decrease with increasing C 12 EO 10 concentration. Three different Nile red (NR) dyes were employed as model solutes and were separately doped into the gels at nanomolar levels. These include a hydrophilic form of NR incorporating an anionic sulfonate group (NRSO 3 - ), a hydrophobic form incorporating a fourteen-carbon alkane tail (NRC 14 ), and commercial NR as an intermediate case. FCS data acquired from the gels revealed that NRSO 3 - diffused primarily in 3D. Its diffusion coefficient exhibited a monotonic decrease with increasing gel concentration and micelle packing density, consistent with confinement of its motions as a result of its exclusion from the micelle cores. NRC 14 exhibited the smallest diffusion coefficient, most likely due to its larger size and enhanced interactions with the micelle cores. NR yielded an intermediate diffusion coefficient and the most anomalous behavior of the three dyes, attributable to its facile partitioning between core and corona regions, and greater participation by 1D diffusion. The results of these studies afford an improved understanding of molecular mass transport through soft-matter nanomaterials like those being developed for use in drug delivery and membrane based chemical separations.

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

Spin-strain interactions under hydrostatic pressure in α-RuCl 3

We investigate the effects of hydrostatic pressure on 𝛼−RuCl 3 , a prototypical material for the Kitaev spin model on a honeycomb lattice with a possible spin-liquid ground state. Using ultrasound measurements at pressures up to 1.16 GPa, we reveal significant modifications of the acoustic properties and the 𝐻−𝑇 phase diagram of this material. Hydrostatic pressure suppresses the three-dimensional magnetic order and induces a dimerization transition at higher pressures. At low pressures, the sound attenuation exhibits a linear temperature dependence, while above 0.28 GPa, it becomes nearly temperature independent, suggesting a shift in the phonon scattering regime dominated by Majorana fermions. These findings provide new insights into spin-strain interactions in Kitaev magnets and deliver a detailed characterization of the 𝐻−𝑇 phase diagram of 𝛼−RuCl 3 under hydrostatic pressure.

Focused ion beam↗
DOE OSTI2022

Investigating Local Carrier Dynamics in PERC Patterned CdTe Solar Cells

Remarkable progress has been achieved in CdTe photovoltaics (PVs) to further improve cell performance while reducing manufacturing costs. Researchers optimized the front contact that leads the short-circuit current density (J sc ) over 31mA/cm 2 (e.g., (Zn, Mg)O buffer layer and/or alloying CdTe with Se). Reese et al. reported an open-circuit voltage (V oc ) over 1 V with group-V doped single crystal, showing the feasibility of increasing V oc in CdTe PV that is fixed to less than 900 mV in most cases. Recent studies suggested that the V oc improvement in polycrystalline CdTe PVs requires well-passivated back contact. One possible strategy is to use CdMgTe alloy (E g > 1.8 eV), where the band-gap offset (ΔE CB ≈ 0.2 eV) reflects minority carrier electrons, thereby decreasing surface recombination. Another strategy is to utilize a stable Al 2 O 3 layer as an electron reflector; experimental works confirmed the improved performance with a conformal Al 2 O 3 layer. In this configuration, the precise control of Al 2 O 3 (≈1 nm) over the entire CdTe layer is essential because the J sc greatly depends on the tunneling current. While promising, both CdMgTe and Al 2 O 3 passivation on CdTe introduce unfavorable valence band offset that blocks the hole transport. Kephart et al. proposed a patterned passixc vation layer (≈ 20 nm Al 2 O 3 ) with point contacts (≈ 3 μm in diameter) that extract hole carriers. This configuration is similar to the passive emitter and rear contact (PERC) design that has been extensively studied in Si PVs. While remarkably high carrier lifetime and photoluminescence (> 10 X) were measured, consistent improvement in V oc and efficiency in this CdTe PERC were not yet observed. At this time, it is unclear how the surface potential is distributed in the presence of patterned Al 2 O 3 reflectors and how this additional field impacts the PV performance for Cu-doped and group-V doped CdTe PERC devices. To tackle this challenge, we develop robust nanofabrication and characterization platform to study CdTe PERC devices consisting of a patterned Al 2 O 3 layer on CdSe( 1-x )Te x (Cu-doped, GrV-doped; x = 0 to 1; CST). First Solar supplies high-quality CST materials. The dark and light I-Vs under 1-sun illumination are extracted for a series of each set of devices. We perform quantitative and qualitative PV analysis for a series of Cu-doped and Gr-V doped PERC devices, showing notably different V oc changes for GrV-doped compared to Cu-doped PERC devices. This project attempts to produce individual contacts on single grains and grain boundaries to measure lateral transport. As a model system, we use a CdSe( 1-x )Te x film on glass prepared by a colossal grain growth (CGG) technique at National Renewable Research Laboratories (NREL). Our preliminary results confirm that intergranular transport is hindered by the electrical barrier formed near grain boundaries. The PERC fabrication and the lateral transport measurement platform developed in this project can easily be applied to other types of advanced CdSe( 1-x )Te x architectures to better understand local photocarrier transport, in turn, providing the fundamental knowledge to improve the performance of CdTe PVs.

14 SOLAR ENERGY↗
DOE OSTI2021

Computationally Guided Discovery of Axis-Dependent Conduction Polarity in NaSnAs Crystals

Most electronic materials exhibit a single dominant charge carrier type, either holes or electrons, along all crystallographic directions. However, there are a small number of compounds, mostly metals, that exhibit simultaneous p-type and n-type conduction behavior along different crystallographic directions. We demonstrate that the experimental discovery of semiconductors with this axis-dependent conduction polarity can be facilitated by identifying a large anisotropy of either the electron or hole effective masses (m*) or both, providing the electron and hole masses dominate along different crystallographic directions. We calculated the layered semiconductor NaSnAs to have a lower electron m* in-plane than the cross-plane and a very large hole m* in-plane and small hole m* cross-plane. We established the growth of >3 mm-sized NaSnAs crystals via Sn flux and confirmed the band gap to be 0.65 eV, in agreement with theory. NaSnAs exhibits p-type thermopowers cross-plane and n-type thermopowers in-plane, confirming that the large anisotropy in the effective mass at the band edges is an excellent indicator for axis-dependent conduction polarity. Altogether, this work shows that the discovery of semiconductors with such a phenomenon can be accelerated by computationally evaluating the anisotropic curvatures of the band edges, paving the way for their future discovery and application.

30 DIRECT ENERGY CONVERSION↗