Search NASA⌕ Search

SEARCH · Search NASA

Results for “Structural simulation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,135 records · Page 63
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Solar cell contacts: quantifying the impact of interfacial layers on selectivity, recombination, charge transfer, and V oc

Interfacial layers (IFLs) are ubiquitous in solar cells, but their precise impact on carrier transfer rates and the relation of these rates to performance metrics and the concepts of selectivity and recombination is lacking. We report the use of a well-defined interdigitated back-contact (IBC) silicon solar cell to determine the precise role of such IFLs. We characterize the action of the common IFL spiro-OMeTAD by making it a third contact to the IBC cell. This architecture creates three solar cells within a single structure that, with numerical simulation, provide the exchange current densities (i.e., charge transfer rates) for electrons ( J 0n ) and holes ( J 0p ) and the quasi Fermi-level splitting in the absorber, which measures the balance of generation and recombination. Further, we describe the relation of V oc to contact recombination, the asymmetry in electron/hole collection rates at a single contact (contact selectivity), and the asymmetry in collection rates of the same carrier at separate contacts (carrier selectivity). Relative to bare gold, neat spiro-OMeTAD reduces J 0n and J 0p (their geometric mean ( J 0n J 0p ) 0.5 decreases by 10 4 ), decreasing contact recombination. Addition of the common dopant Li-TFSI and air increase J 0p / J 0n by 10 6 with little effect on ( J 0n J 0p ) 0.5 , increasing contact hole selectivity. The significant increase in V oc observed by introducing spiro-OMeTAD/Li-TFSI IFLs into the cells studied, however, is due to an increase in the carrier selectivity rather than the contact selectivity or recombination of the spiro-OMeTAD-modified contact. Operando measurements further show voltage-dependent changes in the J 0 s, demonstrating that spiro-OMeTAD contributes to hysteresis. Broader context : Solar cells operate by photogenerating excess charge carriers in an absorber material and, in competition with recombination, asymmetrically extracting them at so-called carrier-selective contacts, one that ideally collects only electrons and the other, only holes. Particularly in emerging photovoltaics such as perovskites, thin layers of organic semiconductors or related materials are introduced between the absorber and contact to improve power conversion efficiency. In terms of interfacial charge transfer, a prevailing view is that such interfacial layers improve performance by helping block the collection of the undesired carrier, considered a form of recombination. We use a novel platform to study the simultaneous impact of spiro-OMeTAD, a common interfacial layer, on the collection of electrons, the collection of holes, and on the recombination of electrons and holes. We quantitatively demonstrate that spiro-OMeTAD layers indeed passivate the contact toward recombination, but that larger improvements in the open-circuit voltage, a key cell metric, can come not only from this but also from increasing the asymmetry of the collection of electrons in the system. Further, operando measurements show transient changes in the properties of spiro-OMeTAD which suggest that it contributes to hysteresis phenomena commonly observed in perovskite and other solar cells.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Advancing density functional tight-binding method for large organic molecules through equivariant neural networks

Semi-empirical quantum-mechanical (QM) methods have become valuable tools for studying complex (bio)molecular systems due to their balance between computational efficiency and accuracy. A key aspect of these methods is their parameterization, which not only governs the reliability of the results but also provides an opportunity to enhance their overall performance. In our previous work [J. Phys. Chem. Lett., 2021, 11, 16], we advanced the third-order semi-empirical density functional tight-binding (DFTB3) method for computing multiple properties of small molecules by developing the machine learning (ML) potential NN rep to bridge the gap between DFTB3 electronic components and those of the hybrid DFT-PBE0 functional. To overcome the limitations of NN rep , we introduce the EquiDTB framework, which leverages physics-inspired equivariant neural networks (NN) to parameterize scalable and transferable many-body Δ TB potentials, replacing the standard pairwise DFTB repulsive potential. This advancement extends the applicability of our ML-corrected DFTB approach to larger molecules and non-covalent systems (including only C, N, O, and H atoms), going beyond the chemical space represented in the training QM datasets. The enhanced performance of EquiDTB over the standard TB methods is demonstrated by the accurate computation of the atomic forces of S66x8 molecular dimers, as well as their interaction energies. Moreover, EquiDTB can be effectively employed to explore the potential energy surfaces of large and flexible drug-like molecules—for example, to determine the minimum energy path between isomers, analyze structural transitions during dynamical simulations, compute vibrational modes, and investigate energetic rankings. The performance for single molecules slightly decreases when the DFTB electronic energy is reduced to first-order but remains superior to standard TB methods. Our work thus demonstrates that an optimal integration of an equivariant NN with QM datasets can advance the DFTB method while maintaining high efficiency, paving the way for reliable (bio)molecular simulations.

Medrano Sandonas, Leonardo [Technische Universität↗
DOE OSTI2017

Irradiation- vs. vitrification-induced disordering: The case of α-quartz and glassy silica

Irradiation and vitrification can both result in the disordering of minerals. However, it remains unclear whether these effects are comparable or if the glassy state represents an upper limit for irradiation-induced disordering. By reactive molecular dynamics simulations, we compare the structure of irradiated quartz to that of glassy silica. We show that although they share some degree of similarity, the structure of irradiated quartz and glassy silica differs from each other, both at the short- (<3 Å) and the medium-range (>3 Å and <10 Å). In particular, the atomic network of irradiated quartz is found to comprise coordination defects, edge-sharing units, and large rings, which are absent from glassy silica. These results highlight the different nature of irradiation- and vitrification-induced disordering.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2018

Alloying effects on superionic conductivity in lithium indium halides for all-solid-state batteries

Alloying of anions is a promising engineering strategy for tuning ionic conductivity in halide-based inorganic solid electrolytes. We explain the alloying effects in Li 3 InBr 6-x Cl x , in terms of strain, chemistry, and microstructure, using first-principles molecular dynamics simulations and novel electronic structure analysis. We find that strain and bond chemistry can be tuned through alloying and affect the activation energy and maximum diffusivity coefficient. The similar conductivities of the x=3 and x=6 compositions can be understood by assuming the alloy separates into Brrich and Cl-rich regions. Phase-separation increases diffusivity at the interface and in the expanded Cl-region, suggesting microstructure effects are critical.

25 ENERGY STORAGE↗
DOE OSTI2021

Reduced bandwidth Compton photons from a laser-plasma accelerator using tailored plasma channels

It has been demonstrated experimentally that laser plasma accelerators can produce multi-100 MeV electron bunches with a few percent energy spread, and from these electrons, multi-MeV quasi-monoenergetic photons have been demonstrated based on Compton up-scattering from a counter-propagating laser. This offers the potential of a high-quality, narrow-bandwidth, compact, photon source with broad application. The bandwidth of the resulting photons depends directly on the distribution of the electron bunch and is limited, in particular, by the bunch divergence (i.e., the spread in transverse velocity angle). At the same time, the ability to decelerate electrons after scattering is important to source deployment. We describe a series of plasma structures that expand and then collimate the electron bunch, reducing its divergence and thus reducing the bandwidth of the scattered photons while enabling both high performance scattering and deceleration. Furthermore, these plasma structures are demonstrated in simulations of the accelerator system, showing the potential to reach few-percent photon spread which is important for applications using nuclear resonance fluorescence.

42 ENGINEERING↗
DOE OSTI2022

Molecular-based analysis of nanoparticle solvation: Classical density functional approach

Proper statistical mechanics understanding of nanoparticle solvation processes requires an accurate description of the molecular structure of the solvent. Achieving this goal with standard molecular dynamics (MD) simulation methods is challenging due to large length scales. An alternative approach to this problem can be formulated using classical density functional theory (cDFT), where a full configurational description of the positions of all the atoms is replaced by collective atomic site densities in the molecule. Using an example of the negatively charged silica-like system in an aqueous polar environment represented by a two-site water model, we demonstrate here that cDFT can reproduce MD data at a fraction of the computational cost. An important implication of this result is the ability to understand how the solvent molecular features may affect the system’s properties at the macroscopic scale. A concrete example highlighted in this work is the analysis of nanoparticle interactions with sizes of up to 100 nm in diameter.

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

RANGE: A robust adaptive nature-inspired global explorer of potential energy surfaces

With the growing demand for realistic representations of chemical structures and the advent of exascale computing, the intelligent sampling of potential energy surfaces and efficient identification of global minima have become more essential but also more feasible. Building on prior studies demonstrating the efficiency of the Artificial Bee Colony (ABC) swarm intelligence algorithm, we report a hybrid metaheuristic framework that integrates the adaptive exploration capabilities of ABC coupled with the exploitation strengths of genetic algorithms (GA) in a scalable, Python-based implementation. The resulting tool, RANGE (Robust Adaptive Nature-inspired Global Explorer), provides seamless interfaces to multiple potential energy evaluators, either directly or via widely used Python libraries, and is designed for high-performance computing environments. We describe the implementation details of RANGE and evaluate its performance, relative to ABC- or GA-alone based algorithms, on a variety of chemical systems, including molecular clusters and heterogeneous surfaces. In conclusion, our results demonstrate RANGE’s efficiency, robustness, and broad applicability in addressing challenging global optimization problems in computational chemistry and materials science.

Algorithms and data structure↗
DOE OSTI2017

Study of blood plasma optical properties in mice grafted with Ehrlich carcinoma in the frequency range 0.1–1.0 THz

In the course of in vitro studies of blood of laboratory animals with progressing Ehrlich carcinoma, we have revealed the change of the blood plasma optical properties in the THz range, which can be used for developing the express diagnostics of the presence of oncological diseases. An applied software package is elaborated that allows the phantoms of biological samples having a complex structure to be numerically simulated and the parameters of the electromagnetic wave reflected from these samples in the THz frequency range to be calculated. (paper)

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2022

Bivalent recognition of fatty acyl-CoA by a human integral membrane palmitoyltransferase

S-acylation, also known as palmitoylation, is the most abundant form of protein lipidation in humans. This reversible posttranslational modification, which targets thousands of proteins, is catalyzed by 23 members of the DHHC family of integral membrane enzymes. DHHC enzymes use fatty acyl-CoA as the ubiquitous fatty acyl donor and become autoacylated at a catalytic cysteine; this intermediate subsequently transfers the fatty acyl group to a cysteine in the target protein. Protein S-acylation intersects with almost all areas of human physiology, and several DHHC enzymes are considered as possible therapeutic targets against diseases such as cancer. These efforts would greatly benefit from a detailed understanding of the molecular basis for this crucial enzymatic reaction. Here, in this study, we combine X-ray crystallography with all-atom molecular dynamics simulations to elucidate the structure of the precatalytic complex of human DHHC20 in complex with palmitoyl CoA. The resulting structure reveals that the fatty acyl chain inserts into a hydrophobic pocket within the transmembrane spanning region of the protein, whereas the CoA headgroup is recognized by the cytosolic domain through polar and ionic interactions. Biochemical experiments corroborate the predictions from our structural model. We show, using both computational and experimental analyses, that palmitoyl CoA acts as a bivalent ligand where the interaction of the DHHC enzyme with both the fatty acyl chain and the CoA headgroup is important for catalytic chemistry to proceed. This bivalency explains how, in the presence of high concentrations of free CoA under physiological conditions, DHHC enzymes can efficiently use palmitoyl CoA as a substrate for autoacylation.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Stripping away ion hydration shells in electrical double-layer formation: Water networks matter

Significance For centuries the double layer at the solid/electrolyte interface has been a central concept in electrochemistry. Today, it is still crucial for virtually all renewable energy storage and conversion technologies. Here, the double-layer formation is probed by THz spectroscopy with ultrabright synchrotron light as a source. Our results capture the molecular details of double-layer formation at positively/negatively charged Au electrodes for an NaCl electrolyte. We reveal a contrasting response applying positive versus negative bias, which is dictated by the interfacial water network and rationalized by accompanying molecular dynamics simulations and electronic-structure calculations. While Na + is directly attracted toward the negatively charged electrode, stripping of the Cl − hydration shell is observed only at larger potential values.

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

Stable calcium-free myocilin olfactomedin domain variants reveal challenges in differentiating between benign and glaucoma-causing mutations

Nonsynonymous gene mutations can be beneficial, neutral, or detrimental to the stability, structure, and biological function of the encoded protein, but the effects of these mutations are often not readily predictable. For example, the β-propeller olfactomedin domain of myocilin (mOLF) exhibits a complex interrelationship among structure(s), stability, and aggregation. Numerous mutations within mOLF are linked to glaucoma; the resulting variants are less stable, aggregation-prone, and sequestered intracellularly, causing cytotoxicity. In this study, we report the first stable mOLF variants carrying substitutions in the calcium-binding site that exhibit solution characteristics indistinguishable from those of glaucoma variants. Crystal structures of these stable variants at 1.8–2.0-Å resolution revealed features that we could not predict by molecular dynamics simulations, including loss of loop structure, helix unwinding, and a blade shift. Double mutants that combined a stabilizing substitution and a selected glaucoma-causing single-point mutant rescued in vitro folding and stability defects. In the context of full-length myocilin, secretion of stable single variants was indistinguishable from that of the WT protein, and the double mutants were secreted to varying extents. In summary, our finding that mOLF can tolerate particular substitutions that render the protein stable despite a conformational switch emphasizes the complexities in differentiating between benign and glaucoma-causing variants and provides new insight into the possible biological function of myocilin.

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

Novel challenges in tracking self-interacting dark matter subhalos

Cosmological N-body simulations are among the primary tools for studying structure formation in the Universe. Analyses of these simulations critically depend on accurately identifying and tracking dark matter subhalos over time. In recent years, several new algorithms have been developed to improve the accuracy and consistency of subhalo tracking in cold dark matter (CDM) simulations. These algorithms should be revisited in the context of new physics beyond gravity, which can modify the evolution and final properties of subhalo populations. In this work, we apply the particle-tracking-based subhalo finder Symfind to velocity-dependent self-interacting dark matter (SIDM) simulations with large cross section amplitudes to assess the performance of particle-tracking methods beyond the CDM paradigm. We find that the core-particle-tracking technique, which is key to the success of these algorithms in CDM, does not always yield accurate results in SIDM. In particular, the interplay between dark matter self-interactions and tidal stripping can cause the diffusion of core particles to larger radii, leading particle-tracking-based algorithms to prematurely lose track of SIDM subhalos. For massive core-expansion subhalos and core-collapse subhalos that experience close or repeated pericentric passages, a significant fraction of core particles can be lost, and particle-tracking-based finders such as Symfind offer no clear advantage over traditional methods that rely on identifying phase-space overdensities. On the other hand, for subhalos with large pericentric distances or fewer, more distant passages, Symfind tends to outperform. These differences depend sensitively on the cross section amplitude and turnover velocity of the underlying SIDM model. We therefore recommend a hybrid approach that leverages the strengths of both techniques to produce complete and robust catalogs of core-expansion and core-collapse SIDM subhalos.

cosmological simulations↗
DOE OSTI2022

Development and construction status of the Mu2e electromagnetic calorimeter mechanical structures

Abstract The Mu2e experiment at Fermilab will search for the CLFV neutrino-less coherent conversion of a muon into an electron in the field of an aluminum nucleus. The observation of this process would be the evidence of physics beyond the Standard Model. Mu2e comprises a straw-tracker, an electromagnetic calorimeter and an external veto for cosmic rays. The calorimeter provides electron identification, a fast trigger and aids track reconstruction. It is a state-of-the-art crystal calorimeter and employs 1340 pure CsI crystals readout by UV-extended SiPM and fast electronics. The design consists of two identical annular disks positioned at the relative distance of 70 cm downstream the target.The hostile Mu2e conditions (total ionizing dose of 12 krad and a neutron fluence of 5 × 10 10 n/cm 2 @ 1 MeVeq (Si)/y, 1 T magnetic field and vacuum level of 10 −4 Torr) posed tight constraints on the mechanical structures and materials choice. The support structure of the two crystal matrices employs two aluminum hollow rings and parts made of open-cell vacuum-compatible carbon fiber. SiPMs and front-end electronics for each crystal are assembled in one mechanical unit inserted in a machined copper holder. The units are supported by a plate made of vacuum-compatible material. The plate integrates the cooling system made of a network of copper lines flowing a low temperature fluid and placed in thermal contact with the copper holders. The DAQ is hosted in aluminum crates positioned on the lateral surface of the disks. The crates also integrate the DAQ electronics cooling system. We review the constraints on the calorimeter structures design, the development of all the structural components, including the simulations that have determined the materials and technological choices and the specifications of the cooling station, components production and quality assurance tests, the procedures for detector assembly, transportation and installation in the experimental area.

Instruments & Instrumentation↗
DOE OSTI2020

Bipolar planetary nebulae from outflow collimation by common envelope evolution

ABSTRACT The morphology of bipolar planetary nebulae (PNe) can be attributed to interactions between a fast wind from the central engine and the dense toroidal-shaped ejecta left over from common envelope (CE) evolution. Here we use the 3D hydrodynamic adaptive mesh refinement (AMR) code AstroBEAR to study the possibility that bipolar PN outflows can emerge collimated even from an uncollimated spherical wind in the aftermath of a CE event. The output of a single CE simulation via the smoothed particle hydrodynamics (SPH) code phantom serves as the initial conditions. Four cases of winds, all with high enough momenta to account for observed high momenta pre-PN outflows, are injected spherically from the region of the CE binary remnant into the ejecta. We compare cases with two different momenta and cases with no radiative cooling versus application of optically thin emission via a cooling curve to the outflow. Our simulations show that in all cases highly collimated bipolar outflows result from deflection of the spherical wind via the interaction with the CE ejecta. Significant asymmetries between the top and bottom lobes are seen in all cases. The asymmetry is strongest for the lower momentum case with radiative cooling. While real post-CE winds may be aspherical, our models show that collimation via ‘inertial confinement’ will be strong enough to create jet-like outflows even beginning with maximally uncollimated drivers. Our simulations reveal detailed shock structures in the shock-focused inertial confinement (SFIC) model and develop a lens-shaped inner shock that is a new feature of SFIC-driven bipolar lobes.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

Sparse Bayesian mass mapping with uncertainties: peak statistics and feature locations

ABSTRACT Weak lensing convergence maps – upon which higher order statistics can be calculated – can be recovered from observations of the shear field by solving the lensing inverse problem. For typical surveys this inverse problem is ill-posed (often seriously) leading to substantial uncertainty on the recovered convergence maps. In this paper we propose novel methods for quantifying the Bayesian uncertainty in the location of recovered features and the uncertainty in the cumulative peak statistic – the peak count as a function of signal-to-noise ratio (SNR). We adopt the sparse hierarchical Bayesian mass-mapping framework developed in previous work, which provides robust reconstructions and principled statistical interpretation of reconstructed convergence maps without the need to assume or impose Gaussianity. We demonstrate our uncertainty quantification techniques on both Bolshoi N-body (cluster scale) and Buzzard V-1.6 (large-scale structure) N-body simulations. For the first time, this methodology allows one to recover approximate Bayesian upper and lower limits on the cumulative peak statistic at well-defined confidence levels.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

Highly dispersive magnons with spin-gap-like features in the frustrated ferromagnetic S = 1 2 chain compound Ca 2 Y 2 Cu 5 O 10 detected by inelastic neutron scattering

In this paper, we report inelastic neutron scattering experiments in Ca 2 Y 2 Cu 5 O 10 and map out the full one-magnon dispersion which extends up to a record value of 53 meV for frustrated ferromagnetic (FM) edge-sharing CuO 2 chain (FFESC) cuprates. A homogeneous spin-1/2 chain model with a FM nearest-neighbor (NN), an antiferromagnetic (AFM) next-nearest-neighbor (NNN) inchain, and two diagonal AFM interchain couplings (ICs) analyzed within linear spin-wave theory (LSWT) reproduces well the observed strong dispersion along the chains and a weak one perpendicularly. The ratio α=| J a 2 / J a 1 | of the FM NN and the AFM NNN couplings is found as ~0.23, close to the critical point α c =1/4 which separates ferromagnetically and antiferromagnetically correlated spiral magnetic ground states in single chains, whereas α c >0.25 for coupled chains is considerably upshifted even for relatively weak IC. Although the measured dispersion can be described by homogeneous LSWT, the scattering intensity appears to be considerably reduced at ~11.5 and ~28 meV. The gaplike feature at 11.5 meV is attributed to magnon-phonon coupling whereas based on density matrix renormalization group simulations of the dynamical structure factor the gap at 28 meV is considered to stem partly from quantum effects due to the AFM IC. Another contribution to that gap is ascribed to the intrinsic superstructure from the distorting incommensurate pattern of CaY cationic chains adjacent to the CuO 2 ones. It gives rise to nonequivalent CuO 4 units and Cu-O-Cu bond angles Φ and a resulting distribution of all exchange integrals. The J 's fitted by homogeneous LSWT are regarded as average values. The record value of the FM NN integral J 1 =24 meV among FFESC cuprates can be explained by a nonuniversal Φ(≠90°) and Cu-O bond length dependent anisotropic mean direct FM Cu-O exchange $\overline{K}_{pd}$~120 meV, similar to a value of 105 meV for Li 2 CuO 2 , in accord with larger values for La 2 CuO 4 and CuGeO 3 (~110 meV) reported by Braden et al . [Phys. Rev. B 54, 1105 (1996)] phenomenologically. Lastly, enhanced $K_{pd}$ values are also needed to compensate a significant AFM J dd ≥ 6 meV from the dd channel, generic for FFESC cuprates but ignored so far.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Superlubric-pinned transition of a two-dimensional solid dusty plasma under a periodic triangular substrate

Here, the superlubric-pinned transition in the depinning dynamics of a two-dimensional (2D) solid dusty plasma modulated by 2D triangular periodic substrates is investigated using Langevin dynamical simulations. When the lattice structure of the 2D solid dusty plasma perfectly matches the triangular substrate, two distinctive pinned and moving ordered states are observed as the external uniform driving force gradually increases from zero. When there is a mismatch between the lattice structure and the triangular substrate, however, on shallow substrates, it is discovered that all of the particles can slide freely on the substrate even when the applied driving force is tiny. This is a typical example of superlubricity, which is caused by the competition between the substrate-particle and particle-particle interactions. If the substrate depth increases further, as the driving force increases from zero, there are three dynamical states consisting of the pinned state, the disordered plastic flow state, and the moving ordered state. In an underdense system, where there are fewer particles than potential well minima, it is found that the occurrence of the three different dynamical states is controlled by the depth of the substrate, which is quantitatively characterized using the average mobility.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗