Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Level density”

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,009 records · Page 56
All sourcesNASA NTRSDOE OSTI
NASA NTRS2002

Application of Accelerometer Data to Mars Odyssey Aerobraking and Atmospheric Modeling

Aerobraking was an enabling technology for the Mars Odyssey mission even though it involved risk due primarily to the variability of the Mars upper atmosphere. Consequently, numerous analyses based on various data types were performed during operations to reduce these risk and among these data were measurements from spacecraft accelerometers. This paper reports on the use of accelerometer data for determining atmospheric density during Odyssey aerobraking operations. Acceleration was measured along three orthogonal axes, although only data from the component along the axis nominally into the flow was used during operations. For a one second count time, the RMS noise level varied from 0.07 to 0.5 mm/s2 permitting density recovery to between 0.15 and 1.1 kg per cu km or about 2% of the mean density at periapsis during aerobraking. Accelerometer data were analyzed in near real time to provide estimates of density at periapsis, maximum density, density scale height, latitudinal gradient, longitudinal wave variations and location of the polar vortex. Summaries are given of the aerobraking phase of the mission, the accelerometer data analysis methods and operational procedures, some applications to determining thermospheric properties, and some remaining issues on interpretation of the data. Pre-flight estimates of natural variability based on Mars Global Surveyor accelerometer measurements proved reliable in the mid-latitudes, but overestimated the variability inside the polar vortex.

Tolson, R. H.↗
NASA NTRS2002

Subcutaneous administration of insulin-like growth factor (IGF)-II/IGF binding protein-2 complex stimulates bone formation and prevents loss of bone mineral density in a rat model of disuse osteoporosis

Elevated serum levels of insulin-like growth factor binding protein-2 (IGFBP-2) and a precursor form of IGF-II are associated with marked increases in bone formation and skeletal mass in patients with hepatitis C-associated osteosclerosis. In vitro studies indicate that IGF-II in complex with IGFBP-2 has high affinity for bone matrix and is able to stimulate osteoblast proliferation. The purpose of this study was to determine the ability of the IGF-II/IGFBP-2 complex to increase bone mass in vivo. Osteopenia of the femur was induced by unilateral sciatic neurectomy in rats. At the time of surgery, 14-day osmotic minipumps containing vehicle or 2 microg IGF-II+9 microg IGFBP-2/100g body weight/day were implanted subcutaneously in the neck. Bone mineral density (BMD) measurements were taken the day of surgery and 14 days later using a PIXImus small animal densitometer. Neurectomy of the right hindlimb resulted in a 9% decrease in right femur BMD (P<0.05 vs. baseline). This loss in BMD was completely prevented by treatment with IGF-II/IGFBP-2. On the control limb, there was no loss of BMD over the 14 days and IGF-II/IGFBP-2 treatment resulted in a 9% increase in left femur BMD (P<0.05). Bone histomorphometry indicated increases in endocortical and cancellous bone formation rates and in trabecular thickness. These results demonstrate that short-term administration of the IGF-II/IGFBP-2 complex can prevent loss of BMD associated with disuse osteoporosis and stimulate bone formation in adult rats. Furthermore, they provide proof of concept for a novel anabolic approach to increasing bone mass in humans with osteoporosis.

NASA Discipline Regulatory Physiology↗
DOE OSTI2020

Validation of gyrokinetic simulations in NSTX and projections for high-k turbulence measurements in NSTX-U

An extensive validation effort performed for a modest-beta NSTX NBI-heated H-mode discharge predicts that electron thermal transport can be entirely explained by electron-scale turbulence fluctuations driven by the electron temperature gradient mode (ETG), both in conditions of strong and weak ETG turbulence drive. Thermal power-balance estimates computed by TRANSP as well as the shape of the high-k density fluctuation wavenumber spectrum and the fluctuation level ratio between strongly driven and weakly driven ETG-turbulence conditions can be matched by nonlinear gyrokinetic simulations and a synthetic diagnostic for high-k scattering. Linear gyrokinetic simulations suggest that the ion-scale instability in the weak ETG condition is close to the critical threshold for the kinetic ballooning mode instability, and nonlinear ion-scale gyrokinetic simulations show that turbulence might be in a state reminiscent of a Dimits' shift regime, opening speculation on the role that ion-scale turbulence might play for the weak ETG condition. A simulation that matched all experimental constraints is chosen to project high-k turbulence spectra in NSTX-U, revealing that the new high-k system [R. Barchfeld et al., Rev. Sci. Instrum. 89, 10C114 (2018)] should be sensitive to density fluctuations from radially elongated streamer structures. In this study, two schemes are designed to characterize the radial and poloidal wavenumber dependence of the density fluctuation wavenumber power spectrum around the streamer peak, suggesting future high-k fluctuation measurements could be sensitive to an asymmetry in the k r spectrum introduced due to the presence of strong background flow shear.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1996

The Effects of Dissipation and Coarse Grid Resolution for Multigrid in Flow Problems

The objective of this paper is to investigate the effects of the numerical dissipation and the resolution of the solution on coarser grids for multigrid with the Euler equation approximations. The convergence is accomplished by multi-stage explicit time-stepping to steady state accelerated by FAS multigrid. A theoretical investigation is carried out for linear hyperbolic equations in one and two dimensions. The spectra reveals that for stability and hence robustness of spatial discretizations with a small amount of numerical dissipation the grid transfer operators have to be accurate enough and the smoother of low temporal accuracy. Numerical results give grid independent convergence in one dimension. For two-dimensional problems with a small amount of numerical dissipation, however, only a few grid levels contribute to an increased speed of convergence. This is explained by the small numerical dissipation leading to dispersion. Increasing the mesh density and hence making the problem over resolved increases the number of mesh levels contributing to an increased speed of convergence. If the steady state equations are elliptic, all grid levels contribute to the convergence regardless of the mesh density.

Eliasson, Peter↗
NASA NTRS2016

Ames S-32 O-16 O-18 Line List for High-Resolution Experimental IR Analysis

By comparing to the most recent experimental data and spectra of the SO2 628 ν1/ν3 bands (see Ulenikov et al., JQSRT 168 (2016) 29-39), this study illustrates the reliability and accuracy of the Ames-296K SO2 line list, which is accurate enough to facilitate such high-resolution spectroscopic analysis. The SO2 628 IR line list is computed on a recently improved potential energy surface (PES) refinement, denoted Ames-Pre2, and the published purely ab initio CCSD(T)/aug-cc-pVQZ dipole moment surface. Progress has been made in both energy level convergence and rovibrational quantum number assignments agreeing with laboratory analysis models. The accuracy of the computed 628 energy levels and line list is similar to what has been achieved and reported for SO2 626 and 646, i.e. 0.01-0.03 cm(exp −1) for bands up to 5500 cm(exp −1). During the comparison, we found some discrepancies in addition to overall good agreements. The three-IR-list based feature-by-feature analysis in a 0.25 cm(exp −1) spectral window clearly demonstrates the power of the current Ames line lists with new assignments, correction of some errors, and intensity contributions from varied sources including other isotopologues. We are inclined to attribute part of detected discrepancies to an incomplete experimental analysis and missing intensity in the model. With complete line position, intensity, and rovibrational quantum numbers determined at 296 K, spectroscopic analysis is significantly facilitated especially for a spectral range exhibiting such an unusually high density of lines. The computed 628 rovibrational levels and line list are accurate enough to provide alternatives for the missing bands or suspicious assignments, as well as helpful to identify these isotopologues in various celestial environments. The next step will be to revisit the SO2 828 and 646 spectral analyses.

Sulfur dioxide↗
NASA NTRS1989

On the populating mechanisms of the autoionising states of Al III ions produced in a Penning ionisation discharge

Spectra emitted by Al II and Al III in the 100-1900-A range from a Penning ionization discharge (PID) have been analyzed in order to study the populating mechanisms of the autoionizing levels of the Al III ion. Electron temperature and density estimates and Al III line intensities obtained using a collisional-radiative model are found to differ from those obtained experimentally, implying that the autoionizing levels of the Na I-like ion may be populated, in the PID plasma, by the inner-shell ionization of excited Mg I-like ions. It is suggested that such autoionizing lines could be used to search for populations inversions in the soft X-ray domain.

Finkenthal, M.↗
NASA NTRS2016

Gene Expression and Structural Skeletal Responses to Long-Duration Simulated Microgravity in Rats

Spaceflight has deleterious effects on skeletal structure and function, specifically causingprofound loss in bone mass, density, and strength, as well as changes in expression levels of genes related to oxidative stress [Hyeon et al., Smith et al.]. It is known that bone resorption remains elevated after spaceflight and that bone density and strength fail to recover completely even years following spaceflight [Smith et al., Carpenter et al.]. However, our current understanding of the signaling pathways and molecular mechanisms that control bone loss and that link oxidative stress, bone resorption, and mechanical unloading of skeletal tissue is incomplete. Here, we aim to examine skeletal responses to simulated long-duration spaceflight on bone loss using the ground-based hindlimb unloading (HU) model in adult (9 months old) male rats. We hypothesized that simulated microgravity leads to the temporal regulation of oxidative-defense genes and pro-osteoclastogenic factors, showing progression and eventual plateau during long-term unloading, and that transient changes at early timepoints in these pathways precede skeletal adaptations to long-duration unloading. We will identify oxidativestress and bone resorption-related changes using global gene expression analysis (Affymetrix arrays) for both acute (within 14 days) and long-term timepoints (90 days). We will also use quantitative PCR to examine changes in expression of genes related to oxidative metabolism (e.g. Nrf2, SOD-1), bone turnover (resorption and formation markers, e.g. TRAP, osteocalcin respectively, SOST), and osteoclastogenesis (e.g. RANKL, OPG) at both early and late timepoints. We will then use detailed microarchitectural and structural analysis through microcomputed tomography to relate gene expression changes with structural changes in bone, expecting that plateaus in gene expression correlate with long-term changes in bone microarchitecture.

Rael, Victoria E.↗
DOE OSTI2023

First large capsule implosions in a frustum-shaped hohlraum

Here, we report on the first indirect-drive implosions driven by a dual conical frustum-shaped hohlraum denoted “frustraum” and the experimental tuning campaigns leading up to two layered implosions. The campaign used 1.2 and 1.4 mm inner radius high density carbon (HDC) capsules and represented the largest HDC capsules to be imploded on the National Ignition Facility via indirect drive. Several techniques were successfully implemented to control the Legendre mode 2 capsule symmetry of the implosions, including changing the wall angle of the frustraum, which is not possible with cylindrical hohlraums. A mode 4 feature was observed and its implications for hotspot mix discussed. Two layered implosions were conducted with 1.2 mm inner radius capsules, the latter of which achieved the highest layered capsule absorbed energy on the National Ignition Facility using only 1.74 MJ of laser energy. The layered implosion results, along with generalized Lawson parameters, suggest that increasing the energy absorbed by the capsule at the expense of long coast times makes it more challenging to achieve ignition and that further reducing coast time (time between end of laser pulse and bang time) closer to the 1 ns level is warranted to improve the areal density and make it easier to achieve the hotspot temperature, alpha heating, and yield amplification required for ignition.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Design Rules of the Mixing Phase and Impacts on Device Performance in High-Efficiency Organic Photovoltaics

In nonfullerene acceptor- (NFA-) based solar cells, the exciton splitting takes place at both domain interface and donor/acceptor mixture, which brings in the state of mixing phase into focus. The energetics and morphology are key parameters dictating the charge generation, diffusion, and recombination. It is revealed that tailoringthe electronic properties of the mixing region by doping with larger-bandgap components could reduce the density of state but elevate the filling state level, leading to improved open-circuit voltage ( V OC ) and reduced recombination. The monomolecular and bimolecular recombinations are shown to be intercorrelated, which show a Gaussian-like relationship with V OC and linear relationship with short-circuit current density ( J SC ) and fill factor (FF). The kinetics of hole transfer and exciton diffusion scale with J SC similarly, indicating the carrier generation in mixing region and crystalline domain are equally important. From the morphology perspective, the crystalline order could contribute to V OC improvement, and the fibrillar structure strongly affects the FF. These observations highlight the importance of the mixing region and its connection with crystalline domains and point out the design rules to optimize the mixing phase structure, which is an effective approach to further improve device performance.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1989

Parallel Optical Random Access Memory (PORAM)

It is shown that the need to minimize component count, power and size, and to maximize packing density require a parallel optical random access memory to be designed in a two-level hierarchy: a modular level and an interconnect level. Three module designs are proposed, in the order of research and development requirements. The first uses state-of-the-art components, including individually addressed laser diode arrays, acousto-optic (AO) deflectors and magneto-optic (MO) storage medium, aimed at moderate size, moderate power, and high packing density. The next design level uses an electron-trapping (ET) medium to reduce optical power requirements. The third design uses a beam-steering grating surface emitter (GSE) array to reduce size further and minimize the number of components.

Alphonse, G. A.↗
DOE OSTI2025

Differentially Private Map Matching (DPMM) v1.0

Human mobility trajectories provide valuable information for developing mobility applications, as they contain diverse and rich information about the users. User mobility data is valuable for various applications such as intelligent transportation systems (ITS), commercial business models, and disease-spread models. However, such spatio-temporal traces may pose a threat to user privacy. GPS trajectories in their raw form are not suitable for transportation studies, as they require matching locations with nearest road links — a process called map-matching. This software implements a differential privacy (DP)-based map-matching algorithm, called DPMM, that generates link-level location trajectories in a privacy-preserving manner to protect users' origin destinations (OD) and travel paths. OD privacy is achieved by injecting Planar Laplace noise to the user OD GPS points. Travel-path privacy is provided with randomized travel path construction using exponential DP mechanism. The injected noise level is selected adaptively, by considering the link density of the location and the functional category of the localized links. For path privacy, our mechanism samples waypoints and selects candidate paths between waypoints. DPMM provides privacy effectively with respect to link density instead of other trajectory samples in the database compared to other privacy mechanisms. Compared to the different baseline models our DP-based privacy model offers closer query responses to the raw data in terms of individual and aggregate trajectory-level statistics with an average at absolute deviation from the baseline for individual statistics on ϵ = 1.0. Beyond individual trajectory statistics, the DPMM outperforms the other benchmark DP-based mechanisms on different aggregate statistics with up to 8x improvement in utility.

Peisert, Sean [Lawrence Berkeley National Laborato↗
DOE OSTI2025

Advancing Grid Resilience through Smart Charge Management: Findings from Maryland’s Pilot

This report presents research findings from a four-year Smart Charge Management (SCM) pilot program conducted by Maryland’s largest electric utilities—Baltimore Gas and Electric (BGE), Potomac Electric Power Company (Pepco), and Delmarva Power & Light (DPL)—to evaluate strategies for optimizing electric vehicle (EV) charging loads and enhancing grid stability. Supported by the U.S. Department of Energy (DOE), Argonne National Laboratory collaborated with all project partners and examined the effectiveness of Time-of-Use (TOU) and Load Balancing (LB) strategies in managing peak demand, deferring costly infrastructure upgrades, and reducing grid constraints at the feeder level. Using charging data from over 4,600 EV drivers, the study analyzed SCM’s impact on the distribution systems of BGE and Pepco, which consists of over 2000 feeders. Unlike prior research that focused on system-wide trends or synthetic feeders, this analysis offers granular, feeder-level insights based on real-world operational data. It highlights how transformer density, load profiles, and infrastructure constraints influence smart charging performance. Results show feeder-level conditions play a crucial role in SCM effectiveness, with most feeders benefiting more from LB, while TOU-based SCM may be sufficient for others. By 2035, LB reduced peak charging loads by 27% on average, compared to 23% under TOU-based SCM, though some feeders saw reductions exceeding 35%, while others experienced minimal impact. Feeders with higher transformer utilization and limited capacity benefited more from LB, which more effectively distributed charging demand during off-peak hours. Beyond reducing grid constraints, SCM offers long-term operational and financial benefits. By shifting EV charging demand strategically, utilities can optimize asset utilization, delay infrastructure investments, and enhance grid performance. In terms of infrastructure upgrade deferrals, at the feeder level, LB consistently reduced peak charging loads and resulting infrastructure upgrade costs, particularly in high EV enrollment areas, decreasing the number of overloaded transformers by up to 35%, while TOU-based SCM achieved 20-30% reductions depending on feeder characteristics. At the system level, LB has the potential to defer total upgrade costs by $\$$186 million for BGE, compared to $\$$159 million under TOU-based SCM. For Pepco, TOU-based SCM performed slightly better, deferring upgrade costs by $\$$30 million, compared to $\$$29 million under LB. Section 4.5 reviews some of the system differences between BGE and Pepco. However, as EV adoption scales, TOU-based SCM will introduce secondary peak charging loads, reinforcing the need for more advanced, adaptive SCM approaches to prevent new grid challenges. As EV adoption continues to grow, feeder-level managed charging strategies will be essential for mitigating grid stress, improving infrastructure efficiency, and maintaining energy affordability for consumers. This report provides critical insights for utilities, Public Utility Commissions (PUCs), and state agencies on the role of feeder-specific smart charging in infrastructure planning, policy development, and grid modernization. The findings underscore the importance of tailored, data-driven SCM solutions that align with local grid conditions, ensuring a resilient, cost-effective transition to increasing EV adoption while safeguarding distribution system performance.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2025

Field-level reconstruction from foreground-contaminated 21-cm maps

Current and upcoming 21-cm experiments will soon be able to map 21-cm spatial fluctuations in three dimensions for a wide range of redshifts. However, bright foreground contamination and the nature of radio interferometry create significant challenges, making it difficult to access rich cosmological information from the Fourier modes that lie within the “foreground wedge”. Here, in this work, we introduce two approaches aiming to reconstruct the full 21-cm density field, including the missing modes in the wedge: (a) a field-level inference under an effective field theory (EFT) framework; (b) a diffusion-based deep generative model trained on simulations. Under the EFT framework, we implement a fully differentiable forward model that maps the initial conditions of matter fluctuations to the observed, foreground-filtered 21-cm maps. This enables a gradient-based sampler to simultaneously sample the initial conditions and bias parameters, allowing a physically motivated mode reconstruction. Alternatively, we apply a variational diffusion model to perform 21-cm density reconstruction at the map level. Our model is trained on semi-numerical simulations over a wide range of astrophysical parameters. Our results from both approaches should provide improved cosmological constraints from the field level and also enable cross-correlation between experiments that have little or no overlapping modes.

cosmological perturbation theory↗
DOE OSTI2026

Physics basis for the reference flat-top plasma scenario in the ST–E1 fusion power plant

As part of the U.S. Department of Energy’s Milestone-Based Fusion Energy Development Program, Tokamak Energy has completed the pre-concept design of the ST–E1 fusion power plant. ST–E1 is envisaged to operate in two phases: a pilot plant phase, targeting sustained net power production of 300 - 500 MWe for a duration >1 hr, followed by a commercial power plant phase targeting steady-state operations and a normalised overnight capital cost of ⩽12 000 $\$$/kWe. The design process adopted was highly iterative, integrating all major plant systems and progressing in a phased fidelity approach. At the pre-conceptual stage, the emphasis has been on exploring the design space, identifying the main system-level trade-offs, and making the key decisions that define the overall plant concept, rather than optimising a single operating point. This paper, part of a focused collection detailing the ST–E1 pre-concept design, addresses the development of a series of reference flat-top plasma operating points for the pilot plant phase. A modelling workflow was established to develop and assess candidate plasma design points and explore key dependencies. The workflow includes integrated core plasma modelling, magnetohydrodynamic (MHD) stability assessment, equilibrium generation, scrape-off-layer and exhaust modelling, heating & current drive design and optimisation, and turbulent transport modelling. Using this framework, the impact of several key parameters on the flat-top operating space was investigated, including the density limit, core radiation fraction and divertor power loading, level of external heating and curent drive power and assumed pedestal characteristics. The MHD stability, controllability and micro-stability characteristics of these plasmas were also analysed. These investigations informed the definition of a set of fully non-inductive, flat-top reference operating points that satisfy the high-level ST–E1 mission, including a low and high density case, a case that is stable to resistive wall modes and a case with reduced divertor power loading.

ST–E1↗
NASA NTRS1986

Short-term variability of nitrogen dioxide in the winter stratosphere

Measurements of limb radiance from the Solar Mesosphere Explorer (SME) satellite are used to infer the NO2 density at the 10- and 16-mbar pressure levels from January 1 to March 31 of 1982. A photochemical-dynamical model is developed using the presently accepted chemistry of the NO(x) family. The dynamical model produces isentropic trajectories which simulate the history of air parcels. From the trajectories the photochemical model calculates NO2 densities, which are compared to those observed by SME. Although the model generally reproduces the spatial and temporal variations rather well, some disagreement was noted for conditions of exceptionally low temperatures. Further analysis indicates that the temperature sensitivity of the N2O5 photolysis cross sections may be overestimated at low temperatures.

Zawodny, J. M.↗
DOE OSTI2022

Coalescence of GaP on V-Groove Si

In recent years, better understanding and control over the formation of crystalline defects during the direct epitaxy of III-V semiconductors on Si substrates via metal organic vapor phase epitaxy (MOVPE) has enabled large gains in III-V-on-Si solar cell efficiency, pointing to pathway to lower-cost, high-performance III-V solar cells. However, such results have only been achieved on costly chemo-mechanically-polished (CMP) Si wafers. The use of V-groove nanopatterned Si substrates has demonstrated similarly high-crystalline-quality III-V-on-Si epitaxy, but also can be combined with lower-cost polishing techniques. Although they offer a potential cost advantage, growth on V-groove substrates adds challenges not present for epitaxy on planar wafers (the III-V material must be coalesced into a thin film after an initial nucleation stage). MOVPE growth conditions that promote highly facet-selective lateral growth needed for coalescence are generally actively avoided for conventional thin fim growth, so growth conditions need to be re-developed for V-groove-based epitaxy. Additionally, coalescence add complexity to the dislocation dynamics related to lattice relaxation, so strategies used to keep threading dislocation density (TDD) low on planar substrates need to be re-tuned for V-grooves. We have studied the morphological evolution and dislocation dynamics of GaP grown on V-groove Si by MOVPE. Growth conditions of V/III=5,000 and Tg=800 degrees C were uniquely found to produce extremely smooth coalesced thin films, with an RMS roughness of 0.2 nm measured by AFM. Additionally, at this growth condition, we identified two regimes of growth determined by the width of the SiNx cap (a remnant of the nanopatterning process) at the top of the V-grooves. For narrow caps, the GaP coalesces into a thin film, and for wide caps, the GaP evolves into {1 1 1}-faceted diamonds that do not coalesce. We suggest the influence of Si from the sidewalls of the V-grooves on the surface reconstruction of the GaP as the mechanism for this effect, with reflection difference spectroscopy (RDS) and Si doping experiments supporting this theory. In addition to morphology, the dislocation dynamics of the system were studied with electron channeling contrast imaging (ECCI) and transmission electron microscopy (TEM). The TDD of the coalesced GaP films was found to be 5 x 10^7 cm^-2 after coalescence via ECCI, a level still too high for high-quality solar cells. However, misfit dislocations crossing multiple grooves greater than 20 micrometers long were observed in ECCI, suggesting that the V-grooves do not block dislocation glide. TEM prior to and after coalescence was used to distinguish between dislocation creation driven by growth conditions and coalescence. Finally, strategies to reduce the dislocation density to levels acceptable for solar cells will be discussed.

GaP↗
DOE OSTI2025

Density Functional Tight Binding Insights into Plasmonic Silver–Platinum Nanoparticles and Alloys for Enhanced Photocatalysis

Developing accurate and efficient Slater-Koster (SK) tight-binding parameter sets is essential for quantum plasmonic studies of alloyed metal nanoparticles, as conventional time dependent density functional theory (TD-DFT) calculations are computationally prohibitive for larger clusters. In this work, we develop and validate density functional tight binding (DFTB) parameter sets for both ground state (GS-SK) and excited state (ES-SK) calculations to study the structural, electronic, and optical properties of silver (Ag), platinum (Pt), and Ag–Pt nanoalloys. Our investigation of the ground state properties demonstrates that the GS-SK parameters enable DFTB to closely reproduce the electronic structures of platinum clusters with diverse sizes and geometries – showing qualitative agreement with DFT for density of states (DOS) profiles and energy levels. The ES-SK parameters accurately describe excited-state properties compared to TD-DFT reference calculations, including the broad, featureless absorption profiles of Pt that are dominated by interband transitions. Using the ES-SK parameters within a real-time TD-DFTB framework, we compute size-dependent optical absorption spectra of Ag, Pt and Ag-Pt nanocubes containing up to 1099 atoms (size ∼4.18 nm). A detailed study of Ag–Pt and Pt-Ag core–shell nanoparticles shows quenching of the Ag plasmon resonance even at monolayer coverage for Ag-Pt, but not for Pt-Ag. We also show how to define submonolayer Ag-core Pt-shell cubic structures that have similar optical properties to those generated experimentally for much larger particles, which offers potential for describing plasmon-enhanced photocatalysis. Collectively, the GS-SK and ES-SK parameter sets provide an accurate, computationally efficient approach for modeling the complex optical and electronic behavior of noble–transition metal nanostructures and their alloys.

SPR↗
DOE OSTI2025

Efficient Modeling of Structural, Electronic, and Optical Properties of Silver and Gold Metal Nanoclusters and Alloys Using Optimized SCC-DFTB Parameters

Computation of optical properties using conventional time-dependent density functional theory (TD-DFT) is time-consuming and memory-intensive. In this study, we investigate the accuracy and efficiency of the density functional tight binding (DFTB) framework with newly optimized Slater–Koster (SK) parameters for modeling the structural, electronic properties, and absorption spectra of silver and gold nanoclusters and their alloys. Our investigation of the ground state (GS) properties demonstrates that the newly developed GS-SK parameters enable DFTB to closely approximate DFT-calculated bond lengths for octahedron, tetrahedron, icosahedra, and truncated octahedron with sizes Ag n /Au n (n = 19, 20, 38, 55), nanoclusters and Ag 20 /Au 20 nanoalloys, with a maximum deviation of approximately 0.15 Å. Formation energy results indicate that the GS-SK parameters can closely estimate changes in formation energies with alloy composition, and the comparison of electronic structures for Ag 20 , Au 20 , and AgAu alloy nanoclusters using the DFTB approximation reveals good agreement in the projected density of states (DOS) profiles and energy levels. A second set of SK parameters, ES-SK, has been developed to describe excited state (ES) properties, including the absorption spectra of silver octahedron Ag 19 , tetrahedral Ag n (n = 20, 56, 84), truncated octahedron Ag 38 , and icosahedra Ag 55 closed-shell clusters and their gold and alloy counterparts over a broad range of alloy compositions. This parametrization uses TD-DFTB calculations and fine-tunes the d and p eigenvalues by comparing them to reference absorption spectra from first-principles TD-DFT. This enables the generation of absorption spectra that closely match the reference spectra when plasmon excitation is dominant, as demonstrated by studying the plasmonic properties of icosahedral Ag n and Au n (n = 309 and 561) nanoparticles. This includes the rapid loss in plasmon quality when Au partially replaces Ag in alloy clusters. Furthermore, these results provide a foundation for addressing computational bottlenecks in plasmonics and with new prospects for applications in the quantum plasmonics for bimetallic alloys.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗