Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “density functional theory development”

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,297 records · Page 72
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Formation of 6H-Ba 3 Ce 0.75 Mn 2.25 O 9 during Thermochemical Reduction of 12R-Ba 4 CeMn 3 O 12 : Identification of a Polytype in the Ba(Ce,Mn)O 3 Family

The resurgence of interest in a hydrogen economy and the development of hydrogen-related technologies has initiated numerous research and development efforts aimed at making the generation, storage, and transportation of hydrogen more efficient and affordable. Solar thermochemical hydrogen production (STCH) is a process that potentially exhibits numerous benefits such as high reaction efficiencies, tunable thermodynamics, and continued performance over extended cycling. Although CeO 2 has been the de facto standard STCH material for many years, more recently 12R-Ba 4 CeMn 3 O 12 (BCM) has demonstrated enhanced hydrogen production at intermediate H 2 /H 2 O conditions compared to CeO 2 , making it a contender for large-scale hydrogen production. However, the thermo-reduction stability of 12R-BCM dictates the oxygen partial pressure (pO 2 ) and temperature conditions optimal for cycling. In this study, we identify the formation of a 6H-BCM polytype at high temperature and reducing conditions, experimentally and computationally, as a mechanism and pathway for 12R-BCM decomposition. 12R-BCM was synthesized with high purity and then controllably reduced using thermogravimetric analysis (TGA). Synchrotron X-ray diffraction (XRD) data is used to identify the formation of a 6H-Ba 3 Ce 0.75 Mn 2.25 O 9 (6H-BCM) polytype that is formed at 1350 degrees C under strongly reducing pO 2 . Density functional theory (DFT) total energy and defect calculations show a window of thermodynamic stability for the 6H-polytype consistent with the XRD results. These data provide the first evidence of the 6H-BCM polytype and could provide a mechanistic explanation for the superior water-splitting behaviors of 12R-BCM.

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

Ab initio study of the structure and properties of amorphous silicon hydride from accelerated molecular dynamics simulations

This paper presents a large-scale ab initio simulation study of amorphous silicon hydride (a-Si 1-x H x ) with an emphasis on the structure and properties of the material across a range of hydrogen concentration by combining accelerated molecular dynamics (MD) simulations with first-principles density-functional calculations. The accelerated MD scheme relied on classical metadynamics, which enabled the development of 2500+ high-quality structural models of a-Si 1-x H x , with system sizes ranging from 150 to 6000 atoms and hydrogen concentrations vary from 6 to 20 at. %. The resulting amorphous networks were found to be completely free from any coordination defects and that they all exhibited a pristine band-gap in their electronic spectrum. The microstructural properties of hydrogen distributions were examined with an emphasis on the presence of isolated and clustered environments of hydrogen atoms. The results were compared with experimental data obtained from X-ray diffraction, infrared spectroscopy and nuclear magnetic resonance studies.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Accurate prediction of global-density-dependent range-separation parameters based on machine learning

In this work, we develop an accurate and efficient XGBoost machine learning model for predicting the global-density-dependent range-separation parameter, ωGDD, for long-range corrected functional (LRC)-ωPBE. This ωGDDML model has been built using a wide range of systems (11 466 complexes, ten different elements, and up to 139 heavy atoms) with fingerprints for the local atomic environment and histograms of distances for the long-range atomic correlation for mapping the quantum mechanical range-separation values. The promising performance on the testing set with 7046 complexes shows a mean absolute error of 0.001 117 a0−1 and only five systems (0.07%) with an absolute error larger than 0.01 a0−1, which indicates the good transferability of our ωGDDML model. In addition, the only required input to obtain ωGDDML is the Cartesian coordinates without electronic structure calculations, thereby enabling rapid predictions. LRC-ωPBE(ωGDDML) is used to predict polarizabilities for a series of oligomers, where polarizabilities are sensitive to the asymptotic density decay and are crucial in a variety of applications, including the calculations of dispersion corrections and refractive index, and surpasses the performance of all other popular density functionals except for the non-tuned LRC-ωPBE. Finally, LRC-ωPBE (ωGDDML) combined with (extended) symmetry-adapted perturbation theory is used in calculating noncovalent interactions to further show that the traditional ab initio system-specific tuning procedure can be bypassed. The present study not only provides an accurate and efficient way to determine the range-separation parameter for LRC-ωPBE but also shows the synergistic benefits of fusing the power of physically inspired density functional LRC-ωPBE and the data-driven ωGDDML model.

Chemistry↗
DOE OSTI2025

Dynamical structure factors of warm dense matter from time-dependent orbital-free and mixed-stochastic-deterministic density functional theory

Abstract We present the first calculations of the inelastic part of the dynamical structure factor (DSF) for warm dense matter (WDM) using time-dependent orbital-free density functional theory (TD-OF-DFT) and mixed-stochastic-deterministic (mixed) Kohn Sham TD-DFT (KS TD-DFT). WDM is an intermediate phase of matter found in planetary cores and laser-driven experiments, where the accurate calculation of the DSF is critical for interpreting x-ray Thomson scattering measurements. Traditional TD-DFT methods, while highly accurate, are computationally expensive, motivating the exploration of TD-OF-DFT and mixed TD-KS-DFT as more efficient alternatives. We applied these methods to experimentally measured WDM systems, including solid-density aluminum and beryllium, compressed beryllium, and carbon–hydrogen mixtures. Our results show that TD-OF-DFT requires a dynamical kinetic energy potential in order to qualitatively capture the plasmon response. Additionally, it struggles with capturing bound electron contributions. In contrast, mixed TD-KS-DFT offers greater accuracy in distinguishing bound and free electron effects, aligning well with experimental data, though at a higher computational cost. This study highlights the trade-offs between computational efficiency and accuracy, demonstrating that TD-OF-DFT remains a valuable tool for rapid scans of parameter space, while mixed TD-KS-DFT should be preferred for high-fidelity simulations. Our findings provide insight into the future development of DFT methods for WDM and suggest potential improvements for TD-OF-DFT.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Atomic-Level Structure of Mesoporous Hexagonal Boron Nitride Determined by High-Resolution Solid-State Multinuclear Magnetic Resonance Spectroscopy and Density Functional Theory Calculations

Mesoporous hexagonal boron nitride (p-BN) has received significant attention over the last decade as a promising candidate for water cleaning/pollutant removal and hydrogen storage applications. In this work, high-resolution solid-state NMR spectroscopy and plane-wave density-functional theory (DFT) calculations are used to obtain an atomic-level description of p-BN. 1 H– 15 N or 1 H– 14 N heteronuclear (HETCOR) correlation experiments recorded with either conventional NMR at room temperature or dynamic nuclear polarization surface-enhanced spectroscopy (DNP-SENS) at ca. 100 K reveal NB 2 H, NBH 2 , NBH 3 + species residing on the edges of BN sheets. Ultra-high field 35.2 T 11 B NMR spectroscopy was used to resolve 11 B NMR signals from BN 3 , BN 2 O x (OH) 1–x (x = 0–1), BNO x (OH) 2–x (x = 0–2), BO x (OH) 3–x (x = 0–3), and BO x (OH) 4–x – (x = 0–4). Importantly, 2D 11 B dipolar double-quantum–single-quantum homonuclear correlation spectra reveal that many pore/defect sites are composed of boron oxide/hydroxide clusters connected to the BN framework through BN 2 O units. 1D and 2D 11 B{ 15 N} HETCOR NMR experiments, in addition to plane-wave DFT calculations of nine different structural models, further confirm the assignment of all NMR signals. The detailed structure determination of the pore and edge/defect sites within p-BN should further enable the rational design and development of next-generation p-BN-based materials. In addition, the techniques outlined here should be applicable to determine structure within other porous and/or boron-based materials.

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

Two Spacers, One Perovskite: Integrating Ruddlesden–Popper and Dion–Jacobson Halide Perovskites

Hybrid organic–inorganic perovskites are a rapidly developing class of materials due to their desirable properties for optoelectronic applications such as their ease of synthesis, solution-processable film formation, tunable band gap, strong photoluminescence, good charge carrier mobilities, and high defect tolerance. Here, we present a novel 2D perovskite motif that seamlessly integrates the structural elements from both Ruddlesden–Popper and Dion–Jacobson halide perovskites. We demonstrate the incorporation of two different organic spacer cations in an ordered manner in 5 novel 2D perovskite iodide materials. Both a cyclic diammonium cation 3-(aminomethyl)piperidinium (3AMP) and a linear alkyl monoammonium cation (Cn = C n H 2n+4 N, n = 4–8) are present in distinct alternating layers, making the new series (Cn) 2 (3AMP)[PbI 4 ] 2 . The crystallization of these materials was optimized through careful temperature control to obtain precise crystal structures via single-crystal X-ray diffraction (XRD) which confirmed the presence of the two cations in distinct layers. The influence of the two spacers on optical properties including the band gaps and photoluminescence spectra are found to more closely resemble (3AMP)PbI 4 than (Cn) 2 PbI 4 , which can be attributed to the amount of distortion imposed by the 3AMP spacer on the lead iodide layers. The findings are supported by density functional theory calculations. The strong photoelectric response of solution-processed thin films shows the potential of these materials in photodetectors or photovoltaics. This unprecedented amalgamation of RP–DJ in (Cn) 2 (3AMP)[PbI 4 ] 2 in a structurally ordered fashion suggests a potential vastly underexplored phase space of 2D perovskites in which there are chemically different spacers in distinct layers of the structure, providing an additional parameter to tune perovskite properties.

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

Zentropy Theory for Transformative Functionalities of Magnetic and Superconducting Materials

The proposed research developed the zentropy theory through applications to complex magnetic materials and superconductors under the hypothesis that the emergent properties of complex magnetic materials and superconductors can be predicted by statistical mechanics of ergodic microstates with their partition functions computed from DFT-predicted free energies. The key objective is to develop approaches to systematically determine the types and number of microstates and the supercell size in DFT-based calculations through convergency of macroscopic functionalities, with the incorporation of our mixed-space approach accounting for the interactions between periodic supercells. In addition to use scientific intuitions to guide the design of important microstates, the key innovation of the proposed research is to integrate the domain knowledge and the material-property-descriptor database (MPDD) with 4 million microstates, which is supported by our deep neural network machine learning models (SIPFENN: structure-informed prediction of formation energy using neural networks) and integrated with our high throughput DFT Tool Kit (DFTTK). For complex magnetic materials, one of the objectives is to develop approaches to calculate short-range ordering from the statistical distribution of each microstate. For superconductors, the divergency of quasiparticle effective mass at a quantum critical point will be investigated, and the superconducting and non-superconducting microstates will be delineated through analysis of electronic band structure, density of states, charge density, and Fermi surface.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Modulating the Electron Affinity of Small Bipyridyl Molecules on Single Gold Nanoparticles for Plasmon-Driven Electron Transfer

Developing controlled platforms for plasmon-driven chemistry is of great importance in catalytic reactions at the nanoscale. We report anion radical formation for five bipyridyl complexes of varying degrees of electron affinity utilizing optically fo-cused intraband (594 nm) and interband (532 nm) pump excitation of single gold nanoparticles. The surface-enhanced Raman scattering (SERS) of anion radicals for the five non-resonant adsorbed molecules 2,2’-bipyridine (22BPY), 4,4’-bipyridine (44BPY), trans-1,2-bis(4-pyridyl)ethylene (BPE), 1,2-bis(4-pyridyl)acetylene (BPA), and 1,2-bis(4-pyridyl)ethane (BPEt) were detected using localized surface-plasmon resonance (LSPR) excitation with 785 nm. The electron affinity of the five bipyridyl complexes were determined using electrochemistry. Molecules with low electron affinity experienced high-er instances of radical anion formation under a plasmon-coupled intraband electron transfer excitation (594 nm) whereas molecules with high electron affinity showed a preference for anion radical formation under direct interband electron trans-fer excitation (532 nm). The lowest unoccupied molecular orbital (LUMO) energy levels for low electron affinity surface-bound molecules (22BPY, BPEt) are on average ~0.43 eV higher for high electron affinity surface-bound molecules (BPA, BPE, 44BPY) as calculated using time-dependent density functional theory, elucidating the importance of plasmon coupling to energy levels that facilitate charge transfer pathways. We also show the ability to ‘activate’ high vs low electron affinity single nanoparticles with the choice of pump excitation wavelength. The findings show the complex interplay between molecular electron affinity, orbital overlap with the density of states of the plasmonic metal, and excitation energetics of the pump laser wavelength. Potential applications of this work include enhanced control over molecular scale catalysis, biosensor design, and solar energy capture.

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

Neural network potential from bispectrum components: A case study on crystalline silicon

In this article, we present a systematic study on developing machine learning force fields (MLFFs) for crystalline silicon. While the main-stream approach of fitting a MLFF is to use a small and localized training set from molecular dynamics simulations, it is unlikely to cover the global features of the potential energy surface. Additionally, to remedy this issue, we used randomly generated symmetrical crystal structures to train a more general Si-MLFF. Furthermore, we performed substantial benchmarks among different choices of material descriptors and regression techniques on two different sets of silicon data. Our results show that neural network potential fitting with bispectrum coefficients as descriptors is a feasible method for obtaining accurate and transferable MLFFs.

36 MATERIALS SCIENCE↗
NASA NTRS1995

The Isolated Bubble Regime in Pool Nucleate Boiling

We consider an isolated bubble boiling regime in which vapour bubbles are intermittently produced at a prearranged set of nucleation site on an upward facing overheated wall plane. In this boiling regime, the bubbles depart from the wall and move as separate entities. Except in the matter of rise velocity, the bubbles do not interfere and are independent of one another. However, the rise velocity is dependent on bubble volume concentration in the bulk. Heat transfer properties specific to this regime cannot be described without bubble detachment size, and we apply our previously developed dynamic theory of vapour bubble growth and detachment to determine this size. Bubble growth is presumed to be thermally controlled. Two limiting cases of bubble evolution are considered: the one in which buoyancy prevails in promoting bubble detachment and the one in which surface tension prevails. We prove termination of the isolated regime of pool nucleate boiling to result from one of the four possible causes, depending on relevant parameters values. The first cause consists in the fact that the upward flow of rising bubbles hampers the downward liquid flow, and under certain conditions, prevents the liquid from coming to the wall in an amount that would be sufficient to compensate for vapour removal from the wall. The second cause is due to the lateral coalescence of growing bubbles that are attached to their corresponding nucleation sites, with ensuing generation of larger bubbles and extended vapour patches near the wall. The other two causes involve longitudinal coalescence either 1) immediately in the wall vicinity, accompanied by the establishment of the multiple bubble boiling regime, or 2) in the bulk, with the formation of vapour columns. The longitudinal coalescence in the bulk is shown to be the most important cause. The critical wall temperature and the heat flux density associated with isolated bubble regime termination are found to be functions of the physical and operating parameters and are discussed in detail.

Buyevich, Y. A.↗
DOE OSTI2025

Predicting nonequilibrium Green’s function dynamics and photoemission spectra via nonlinear integral operator learning

Understanding the dynamics of nonequilibrium quantum many-body systems is an important research topic in a wide range of fields across condensed matter physics, quantum optics, and high-energy physics. However, numerical studies of large-scale nonequilibrium phenomena in realistic materials face serious challenges due to intrinsic high-dimensionality of quantum many-body problems and the absence of time-invariance. The nonequilibrium properties of many-body systems can be described by the dynamics of the correlator, or the Green's function of the system, whose time evolution is given by a high-dimensional system of integro-differential equations, known as the Kadanoff–Baym equations (KBEs). The time-convolution term in KBEs, which needs to be recalculated at each time step, makes it difficult to perform long-time numerical simulation. In this paper, we develop an operator-learning framework based on recurrent neural networks (RNNs) to address this challenge. We utilize RNNs to learn the nonlinear mapping between Green's functions and convolution integrals in KBEs. By using the learned operators as a surrogate model in the KBE solver, we obtain a general machine-learning scheme for predicting the dynamics of nonequilibrium Green's functions. Besides significant savings per each time step, the new methodology reduces the temporal computational complexity from $O(N_t^3)$ to $O(N_t)$ where N t is the number of steps taken in a simulation, thereby making it possible to study large many-body problems which are currently infeasible with conventional KBE solvers. Through various numerical examples, we demonstrate the effectiveness of the operator-learning based approach in providing accurate predictions of physical observables such as the reduced density matrix and time-resolved photoemission spectra. Moreover, our framework exhibits clear numerical convergence and can be easily parallelized, thereby facilitating many possible further developments and applications.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Kinetic and X-ray Absorption Spectroscopic Analysis of Catalytic Redox Cycles over Highly Uniform Polymetal Oxo Clusters

Metal–organic framework materials (MOFs) offer an opportunity for investigating catalytic properties of polymetal oxo clusters that are highly well defined and uniform in nature, in contrast to other classes of catalysts that may exhibit a propensity toward active site heterogeneity. We report herein a kinetic and X-ray absorption spectroscopy (XAS) analysis of the two-electron oxidation of CO over divalent metal sites in MIL-100(M = Fe, Cr) (MIL = Materials of Institut Lavoisier) materials carrying μ 3 -oxo bridged trimers, and connect observations about the kinetic relevance of redox steps to density functional theory (DFT) predictions published previously. The high degree of uniformity evident from in situ titration measurements leads to a congruence in mechanistic inferences made from steady-state catalytic, transient stoichiometric, isotopic exchange, and isotopic tracer data that all point to a sequential mechanism comprised of separate oxidation and reduction half-cycle steps conjoined by an active oxygen intermediate. In situ XAS data reinforce mechanistic conclusions derived from kinetic analysis, and suggest that the active oxygen intermediate may be more appropriately characterized as an iron-oxyl (Fe 3+ –O – ) rather than an iron-oxo (Fe 4+ =O 2– ) species. The Cr analogue of MIL-100 exhibits contrasting rate features that can be rationalized using an identical sequence of steps as MIL-100(Fe), but with a highly dissimilar set of kinetic parameters that can also be validated using transient stoichiometric experiments. The larger coverages of active oxygen intermediates on MIL-100(Cr) are consistent with predictions from prior DFT studies that suggest more stable and less reactive active oxygen species for metals with lower d-electron counts, and point to metal identity as a lever for precisely controlling the kinetic relevance of oxidation and reduction half-cycles in catalytic redox sequences over polymetal oxo clusters. Importantly, the results presented point to the utility of developing broadly applicable structure–catalytic property relationships over MOF nodes specifically, and highly uniform catalysts more generally.

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

Prediction of low-Z collinear and noncollinear antiferromagnetic compounds having momentum-dependent spin splitting even without spin-orbit coupling

The over 60 years old Rashba-Dresselhaus effect predicts spin-orbit coupling (SOC) induced momentum-dependent spin splitting and spin polarization in materials with noncentrosymmetric structures. Strong SOC induced effects usually require high-atomic number (Z) elements such as rare-earth elements. It has recently been pointed out that antiferromagnets could hold SOC-independent spin splitting and spin polarization. In the present work we develop the spatial and magnetic symmetry conditions enabling such antiferromagnet (AFM)-induced spin splitting, dividing the 1651 magnetic space groups into seven different spin splitting prototypes (SST-1 to SST-7). This analysis places the physics of AFM spin splitting (SST-4) within the broader context of symmetry conditions that enable the more familiar forms of spin splitting, such as ferromagnetic Zeeman effect (SST-5), nonmagnetic no spin splitting (SST-6), and the nonmagnetic Rashba and Dresselhaus effects (SST-7). The AFM-induced spin splitting and spin polarization do not necessarily require breaking of inversion symmetry or the presence of SOC, hence can exist even in centrosymmetric, low-Z light element compounds, considerably broadening the material base for spin manipulations. We use the “inverse design” approach of first formulating the target property (here, spin splitting in low-Z compounds not restricted to low symmetry structures), then derive the enabling physical design principles—the magnetic symmetry conditions—to search realizable compounds that satisfy these a priori design principles. This process uncovers 422 magnetic space groups (160 centrosymmetric and 262 noncentrosymmetric) that could hold AFM-induced, SOC-independent spin splitting and spin polarization. We then search for stable compounds following such enabling symmetries. We investigate the electronic and spin structures of some selected prototype compounds by density functional theory (DFT) and find spin textures that are different than the traditional Rashba-Dresselhaus patterns and exist even in the absence of SOC effect. We provide the DFT results for all antiferromagnetic spin splitting prototypes (SST-1, SST-2, SST-3, SST-4), and concentrate on two limits of SST-4 that are particularly unusual: When spin splitting is momentum dependent (just like the Rashba effect) but is enabled in antiferromagnets even in the absence of SOC in the Hamiltonian. This includes examples of (a) centrosymmetric SST-4A compounds (e.g., orthorhombic LaMnO 3 illustrating collinear AFM, as well as cubic NiS 2 illustrating noncollinear AFM) and (b) noncentrosymmetric SST-4B compounds (e.g., rhombohedral MnTiO 3 illustrating collinear AFM and hexagonal ScMnO 3 illustrating noncollinear AFM). The symmetry design principles outlined here, along with their transformation into an inverse design material search approach and DFT verification, could open the way to their experimental examination.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Correlated Displacement of Dynamic Elastic Dipoles Produces Nonclassical Electrostriction in Zr-Doped Ceria

By combining experimental data with density functional theory-based ab initio molecular dynamics modeling, this work provides evidence that nonclassical electrostriction in isovalent Zr-doped ceria is due to the correlated anharmonic motion of dynamic elastic dipoles associated with multiple [ZrO 8 ]-local bonding units with a high Zr concentration (Zr 0.1 Ce 0.9 O 2 ). Introduction of 0.5 mol % trivalent or divalent codopants (Sc, Yb, La, or Ca) reduces the longitudinal electrostriction strain coefficient by more than a factor of 10, produces a 3-fold decrease in the relative dielectric permittivity, and increases the elastic modulus. Since these changes depend neither on the radius nor on the valency of the codopant, we conclude that the responsible species are charge-compensating oxygen vacancies (VO). For trivalent dopants (Do 0.005 Zr 0.1 Ce 0.895 O 1.9975 ), oxygen vacancies are present at a concentration ratio 1:40 with respect to Zr, giving, for random distribution, a characteristic interaction distance of ≤2.3 unit cells (1.2 nm). Oxygen vacancies participate in [ZrO 7 -V O ] local bonding units, disrupting the correlated dynamic displacements of the connected [ZrO 8 ]-local bonding units. Finally, such correlated motion of dynamic elastic dipoles may also explain the exponential increase in the longitudinal electrostriction strain coefficient with an increase in Zr concentration to <0.2 mole fraction and must be taken into account for further development of nonclassical electrostrictors based on Zr-doped ceria.

36 MATERIALS SCIENCE↗
DOE OSTI2020

High-Resolution Three-Dimensional Sculpting of Two-Dimensional Graphene Oxide by E-Beam Direct Write

On demand switchable additive/subtractive patterning of 2D nanomaterials is an essential capability for developing new concepts of functional nanomaterials and their device realizations. Traditionally, this is performed via a multi-step process using photoresist coating and patterning by conventional photo or electron beam lithography, which is followed by bulk dry/wet etch or deposition. This limits the range of functionalities and structural topologies that can be achieved, as well as increases the complexity, cost, and possibility of contamination which are significant barriers to device fabrication from highly sensitive 2D materials. Focused electron beam induced processing (FEBIP) enables a material chemistry/site-specific, high-resolution multi-mode atomic scale processing and provides unprecedented opportunities for direct-write, single-step surface patterning of 2D nanomaterials with an in-situ imaging capability. It allows for realizing a rapid multi-scale/multi-mode approach, ranging from an atomic scale manipulation (e.g., via targeted defect introduction as an active site) to a large-area surface modification on nano and micro scales, including patterned doping and material removal/deposition with 2D (in-plane)/3D (out-of-plane) control. In this work, we report on a new capability of FEBIP for nanoscale patterning of graphene oxide via removal of oxygenated carbon moieties with no use of reactive gas required for etching complemented by carbon atom deposition using focused electron beam. The mechanism of experimentally observed phenomena is explored using the density functional theory (DFT) calculations, revealing that interactions of e-beam liberated reactive oxygen radicals with carbon atoms on graphene basal plane lead to creation of atomic vacancies in the material. The reaction by-products are volatile carbon di-oxides, which are dissociated and volatilized from the graphene oxide surface functional groups by interactions with an energetic focused electron beam. Along with selective “subtractive” patterning of graphene oxide, the same electron beam with increased irradiation doses can deposit out-of-plane 3D carbon nanostructures on top of or around the 2D etched pattern, thus forming a hybrid 2D/3D nanocomposite with feature control down to a few nanometers. Furthermore, this in-operando dual nanofabrication capability of FEBIP is unmatched by any other nanopatterning techniques and opens a new design window for forming 2D/3D complex nanostructures and functional nanodevices.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Alloying multiple halide perovskites on the same sublattice in search of stability and target band gaps

Single-component halide perovskites (HPs) rarely satisfy all the necessary criteria for optoelectronic applications, such as achieving an optimal band gap while maintaining high chemical and structural stability. Alloying halide perovskites has emerged as a promising strategy, not only to enhance stability but also to fine-tune their electronic and optical properties. In this work, we explore multiple degrees of freedom in alloy design, considering different substitution sublattices sites (A, B, or X in ABX3 perovskites), various chemical species (isovalent and hetero-valent elements), and multi-component compositions on a given sublattice. Using first-principles calculations based on density functional theory (DFT), we investigate how compositional variations influence the electronic (band gap) and structural properties (mixing enthalpy) of HP alloys. Our approach employs the polymorphous cell model, allowing full local relaxation which breaks local symmetry while preserving global cubic symmetry—an essential framework for accurately modeling HPs. Our results reveal that X-site mixing (halogen substitution) primarily affects the valence band maximum, allowing target band gap engineering. Additionally, variations in halogen radii introduce internal strain through octahedral distortions, influencing the mixing enthalpy. A-site substitution, while not directly contributing to the band edge states, modifies structural stability via volume effects, indirectly impacting the band gap. B-site alloying plays a dominant role in band gap modulation, leading to either positive or negative band gap bowing. Specifically, isovalent B-site mixing (Sn–Pb) induces strong positive bowing, where the alloy band gap is smaller than the average gap of parent compounds, whereas hetero-valent mixing (Cd–Pb) results in pronounced negative bowing. As an aside, we investigate the competition between the excess energy of disordered alloys vs. that of long-range ordered double perovskites of the same compositions, seeking examples of ordered phases emerging from disordered alloys. Furthermore, our findings provide fundamental insights into the electronic and structural behavior of HP alloys, offering valuable design principles for the development of stable and efficient materials for next-generation photovoltaic and optoelectronic devices.

14 SOLAR ENERGY↗
DOE OSTI2025

Development of Parameters for the Particle Size Distribution of TATB

Laser light scattering (LLS), manual counting of scanning electron microscopy (SEM) images, and reduction of SEM images using ImageJ (open source program) to determine Feret (caliper) diameters were applied to determine particle size distribution (PSD) of four preparations of 1,3,5‐triamino‐2,4,6‐trinitrobenzene (TATB) and yttria‐stabilized zirconia (YSZ), an SEM certified standard. Mie theory was used to reduce the LLS data. The spherical nature of the YSZ made it a good candidate for LLS. Variations in n , the refractive index, and iκ, the imaginary component, produced very little change in the PSD. However, changing the carrier liquid from H 2 O to a 40% aqueous sucrose solution, thereby changing the carrier refractive index, n 0 , substantially affected the PSD. The Mie complex refractive indices for the YSZ were n = 2.200, iκ = 0.100, with a 40% aqueous sucrose solution, n 0 = 1.400. The triclinic crystal structure of TATB made refractive index determinations more difficult, so a study was conducted varying Mie parameters and comparing them to the same data reduced using the Fraunhofer theory. Changing the n and κ parameters produced PSD with a small concentration of particles less than 1 µm in size or none in this range. SEM images, Feret data, manual counting, and Fraunhofer data reduction indicate particles less than 1 µm are probably < 5% in concentration. The final selection of Mie parameters for TATB was n = 2.283, iκ = 0.1, and suspension medium, n 0 = 1.330. Finally, computations, using density functional theory produced similar parameters.

Feret diameter↗