Search NASA⌕ Search

SEARCH · Search NASA

Results for “dependability”

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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Transformation of TiN to TiNO Films via In-Situ Temperature-Dependent Oxygen Diffusion Process and Their Electrochemical Behavior

Titanium oxynitride (TiNO) thin films represent a multifaceted material system applicable in diverse fields, including energy storage, solar cells, sensors, protective coatings, and electrocatalysis. This study reports the synthesis of TiNO thin films grown at different substrate temperatures using pulsed laser deposition. A comprehensive structural investigation was conducted by X-ray diffraction (XRD), X-ray photoelectron spectroscopy (XPS), Non-Rutherford backscattering spectrometry (N-RBS), and X-ray absorption spectroscopy (XAS), which facilitated a detailed analysis that determined the phase, composition, and crystallinity of the films. Structural control was achieved via temperature-dependent oxygen in-diffusion, nitrogen out-diffusion, and the nucleation growth process related to adatom mobility. The XPS analysis indicates that the TiNO films consist of heterogeneous mixtures of TiN, TiNO, and TiO2 phases with temperature-dependent relative abundances. The correlation between the structure and electrochemical behavior of the thin films was examined. The TiNO films with relatively higher N/O ratio, meaning less oxidized, were more electrochemically active than the films with lower N/O ratio, i.e., more oxidized films. Films with higher oxidation levels demonstrated enhanced crystallinity and greater stability under electrochemical polarization. These findings demonstrate the importance of substrate temperature control in tailoring the properties of TiNO film, which is a fundamental part of designing and optimizing an efficient electrode material.

Cherono, Sheilah↗
DOE OSTI2025

The Principle of Maximum Conformality Correctly Resolves the Renormalization-Scheme-Dependence Problem

In this paper, we clarify a serious misinterpretation and consequent misuse of the Principle of Maximum Conformality (PMC), which also can serve as a mini-review of PMC. In a recently published article, P. M. Stevenson has claimed that “the PMC is ineffective and does nothing to resolve the renormalization-scheme-dependence problem”, concluding incorrectly that the success of PMC predictions is due to the PMC being a “laborious, ad hoc, and back-door” version of the Principle of Minimal Sensitivity (PMS). We show that such conclusions are incorrect, deriving from a misinterpretation of the PMC and an overestimation of the applicability of the PMS. The purpose of the PMC is to achieve precise fixed-order pQCD predictions, free from conventional renormalization schemes and scale ambiguities. We demonstrate that the PMC predictions satisfy all the self-consistency conditions of the renormalization group and standard renormalization-group invariance; the PMC predictions are thus independent of any initial choice of renormalization scheme and scale. The scheme independence of the PMC is also ensured by commensurate scale relations, which relate different observables to each other. Moreover, in the Abelian limit, the PMC dovetails into the well-known Gell-Mann–Low framework, a method universally revered for its precision in QED calculations. Due to the elimination of factorially divergent renormalon terms, the PMC series not only attains a convergence behavior far superior to that of its conventional counterparts but also deftly curtails any residual scale dependence caused by the unknown higher-order terms. This refined convergence, coupled with its robust suppression of residual uncertainties, furnishes a sound and reliable foundation for estimating the contributions from unknown higher-order terms. Anchored in the bedrock of standard renormalization-group invariance, the PMC simultaneously eradicates the factorial divergences and eliminates superfluous systematic errors, which inversely provides a good foundation for achieving high-precision pQCD predictions. Consequently, owing to its rigorous theoretical underpinnings, the PMC is eminently applicable to virtually all high-energy hadronic processes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Polarization-dependent photoluminescence of Ce-implanted MgO and MgAl 2 O 4

Since the qubit's performance of solid-state spin centers depends highly on the host material, spin centers using new host materials may offer new qubit applications. We investigate the optical properties of Ce-implanted MgO and MgAl 2 O 4 as potential materials holding the optically accessible qubit. We find that the photoluminescence of Ce-implanted MgAl 2 O 4 is more than 10 times brighter than that of Ce-implanted MgO and observe polarization-dependent emission of Ce center in MgAl 2 O 4 with 2% at 4 K under 500 mT, suggesting that the properties required for initializing and reading the state of the spin qubit have been achieved.

42 ENGINEERING↗
DOE OSTI2025

DESI Emission-line Galaxies: Clustering Dependence on Stellar Mass and [O II] Luminosity

We measure the projected two-point correlation functions of emission-line galaxies (ELGs) from the Dark Energy Spectroscopic Instrument One-Percent Survey and model their dependence on stellar mass and [O II] luminosity. We select ∼180,000 ELGs with redshifts of 0.8 < z < 1.6, and define 27 samples according to cuts in redshift and both galaxy properties. Following a framework that describes the conditional [O II] luminosity–stellar mass distribution as a function of halo mass, we simultaneously model the clustering measurements of all samples at fixed redshift. Based on the modeling result, most ELGs in our samples are classified as central galaxies, residing in halos of a narrow mass range with a typical median of ∼10 12.2−12.4 h −1 M ⊙ . We observe a weak dependence of clustering amplitude on stellar mass, which is reflected in the model constraints and is likely a consequence of the 0.5 dex measurement uncertainty in the stellar mass estimates. The model shows a trend between galaxy bias and [O II] luminosity at high redshift (1.2 < z < 1.6) that is otherwise absent at lower redshifts.

Cosmology↗
DOE OSTI2026

Dependence of CCN closure relationship with organic fraction from two airborne field campaigns over mid-latitude land and ocean

This study investigates the relationship between measured and calculated cloud condensation nuclei (CCN) number concentration and its dependence with organic fraction utilizing aircraft observations from The Aerosol and Cloud Experiments in the Eastern North Atlantic (ACE-ENA, 2017–2018) and The Holistic Interactions of Shallow Clouds, Aerosols, and Land Ecosystems (HI-SCALE, 2016) campaigns, which represent midlatitude marine and continental environments, respectively. For the ACE-ENA marine region, aerosol and CCN concentrations were significantly higher in summer than in winter, whereas at continental site for HI-SCALE, aerosol and CCN concentrations showed no pronounced differences between spring and autumn. Using aerosol chemical composition and number size distribution data, CCN concentrations at various supersaturations are calculated based on Köhler theory and then compared with observations from CCN counter. The results show that CCN closure performs well at both sites with a slight overestimation, with mean closure ratio (CR) of 1.13 and 1.17, respectively. Further investigation reveals that CR at lower supersaturation perform better than that at higher supersaturation. The dependence of CR on organic mass fraction (MForg) varies by environment: for marine aerosols, CR decreases with increasing organic fraction at lower supersaturations, whereas continental aerosols exhibit a consistent overestimation, with CR decreasing as organic fraction increases at higher supersaturations. This study provides key insights into CCN characteristics over midlatitude marine and continental environments, emphasizing the necessity of incorporating size-resolved chemical composition and mixing states into future model parameterizations, and contributing to a better understanding of aerosol–cloud interactions.

ACE-ENA field campaign↗
DOE OSTI2026

DESI Data Release 2 ELGs: Property-dependent subsamples, imaging systematics, and clustering

Using emission-line galaxies (ELGs) from the Dark Energy Spectroscopic Instrument (DESI) Data Release 2, we evaluate a property-dependent correction to imaging systematics. We derive systematic weights following the same linear regression method used for other DESI tracers, but do so separately on ELG subsamples to provide a physically-informed alternative to the fiducial, neural-network-based approach. In doing so, we show that the deeper imaging in the Dark Energy Survey (DES) footprint leads to a higher overall number density but a lack of targets with extreme $g-r$ and $r-z$ colors. ELGs in the DES region also show a distinct redshift distribution when subsampled by position in the $g-r$ vs. $r-z$ plane. To address these effects, we implement a separate treatment of the DES footprint within the DESI catalog production pipeline, which is generally well-motivated and, in some cases, imperative for accurate clustering measurements. With DES treated separately, we find that property-dependent systematic weights further mitigate spurious clustering signal in $\sim$10% of subsamples, while the fiducial scheme remains optimal for the full sample.

Hagen, T. [Utah U.]↗
DOE OSTI2025

Resolution dependence of the turbulent atmospheric boundary layer in global storm-resolving climate simulations

The current generation of state-of-the-art global climate models are being run at increasingly higher horizontal resolutions, with the goal of resolving organised deep convection explicitly. How a kilometre-scale resolution impacts the representation of the atmospheric boundary layer is, however, not well known. Using statistical analysis on global fields as well as high-frequency data at selected locations, produced with the Integrated Forecasting System (IFS) model for the Next Generation Earth-system Models (nextGEMS) project, we investigate the horizontal resolution dependence of some boundary-layer processes. We find that a change in resolution from 9 to 2.8 km causes no substantial changes to boundary-layer properties and processes at most of the locations studied, although some global changes are detected that indicate circulation changes. Small changes to the boundary-layer depth and structure are found in the Tropics. The short simulation length and lack of data for optimal boundary-layer analysis limits the conclusions, especially in relation to the connection between the boundary layer and the atmosphere general circulation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Temperature-Dependent Mechanical Properties of Ni-Based Concentrated Alloys: Insights from First-Principles Calculations

The present work focuses on predicting temperature-dependent mechanical properties of Ni-based concentrated alloys Ni 18 Cr 10 Co 10 Fe 6 M 4 (abbreviated by X 44 M 4 , with M = Al, V, Mn, Fe, Nb, Mo, and W) using density functional theory (DFT). These predictions are based on shear (plastic) and elastic deformations, utilizing the special quasirandom structure (SQS), the phonon-based quasiharmonic approach (QHA), and the quasistatic approach. The resulting properties include coefficient of thermal expansion via QHA, ideal shear strength (τ IS ), and stable and unstable stacking fault energies (γ SF and γ US ) through pure alias shear deformation, and elastic constants (c ij ), bulk modulus (B 0 ), and shear modules (G 0 ) via elastic deformation. Notably, predicting accurate γ SF is challenging due to uncertainties that can exceed the γ SF values. τ IS and γ US exhibit a strong linear relationship, enabling the accurate prediction of γ US based on the precisely determined τ IS . All mechanical properties of X 44 M 4 decrease with increasing temperature, except for some γ SF cases such as X 44 M 4 with M = V, Mn, Fe, Mo, and W. Among the X 44 M 4 alloys, X 44 Nb 4 exhibits the lowest τ IS , γ US , and G 0 values, and the highest B 0 /G 0 ratio, while X 44 Mn 4 has the lowest B 0 and B 0 /G 0 ratio. We found that volume is a crucial descriptor for understanding and modeling mechanical properties (except B0 and maybe also γ SF ) affected by alloying elements and temperature. Ni-based dilute alloys (e.g., Ni 11 M 1 and Ni 31 M 1 ) and concentrated alloys (e.g., X 44 M 4 ) show similar trends in mechanical properties influenced by alloying elements and temperature, simplifying the analysis and design of Ni-based alloys.

Elastic properties↗
DOE OSTI2025

Coverage-dependent structures and thermodynamic stability of intercalated Gd layers beneath buffer-layer graphene on SiC(0001)

Electronic properties of two-dimensional (2D) materials are strongly influenced by their atomic arrangements, making the theoretically-aided characterization of experimentally-synthesized 2D structures crucial. Using first-principles density functional theory, we analyze nearly 200 configurations of intercalated Gd layers beneath buffer-layer graphene on SiC(0001) over a Gd coverage range of 0.01 < θ < 1.2. By fully relaxing selectively-constructed configurations at each coverage within a large, low-strain supercell, we determine the coverage dependence of the chemical potential for intercalated Gd structures. Thermodynamically-preferred configurations below θ ≈ 0.8 form single-atom-thick monolayers, while 3D-like or multilayer structures emerge beyond θ ≈ 0.9. Most structures are amorphous-like, including the configuration at the chemical potential minimum around θ ≈ 0.4. In contrast, a strongly stretched Gd(0001)-like monolayer at θ = 1/3 and a nearly perfect Gd(0001) monolayer at θ = 1 are significantly less favorable with 0.16 eV and 0.82 eV higher chemical potentials above the minimum, respectively. Furthermore, the graphene layer decoupled by intercalated Gd near the chemical potential minimum is significantly flatter compared to its morphology above intercalated 3D structures at higher coverages and nearly isolated Gd atoms in the lowest coverage region. In conclusion, these findings align with our experimental results and underscore the need for further research on this unique intercalated system, which holds significant potential for diverse applications.

36 MATERIALS SCIENCE↗
DOE OSTI2024

QRCODE: Massively parallelized real-time time-dependent density functional theory for periodic systems

We present a new software module, QRCODE (Quantum Research for Calculating Optically Driven Excitations), for massively parallelized real-time time-dependent density functional theory (RT-TDDFT) calculations of periodic systems in the open-source Qbox software package. Our approach utilizes a custom implementation of a fast Fourier transformation scheme that significantly reduces inter-node message passing interface (MPI) communication of the major computational kernel and shows impressive scaling up to 16,344 CPU cores. In addition to improving computational performance, QRCODE contains a suite of various time propagators for accurate RT-TDDFT calculations. As benchmark applications of QRCODE, we calculate the current density and optical absorption spectra of hexagonal boron nitride (h-BN) and photo-driven reaction dynamics of the ozone-oxygen reaction. We also calculate the second and higher harmonic generation of monolayer and multi-layer boron nitride structures as examples of large material systems. Our optimized implementation of RT-TDDFT in QRCODE enables large-scale calculations of real-time electron dynamics of chemical and material systems with enhanced computational performance and impressive scaling across several thousand CPU cores.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Containment strategy for subsurface H 2 storage based on time-dependent soft solids

Here, we propose an innovative containment strategy based on time-dependent soft solids, namely 2 wt% Laponite suspensions, to reinforce natural subsurface seals and engineer flow barriers, with an eye toward making H 2 subsurface storage scalable and geographically agnostic. This suspension can be injected at its initial low viscosity and elasticity into a porous medium, allowing for easy pumping and targeted delivery. Once inside the target zone, it matures into a soft solid with much higher viscosity and elasticity, acting as a potential flow barrier. We discuss the rheological properties of the suspensions and demonstrate that hydrogen does not adversely affect their microstructure, but rather increases the suspensions’ viscosity and elasticity. Moreover, the suspensions enhance the rock samples’ compressive strength, while hydrogen exposure increases their stiffness and ductility. The ability of rock samples saturated with the suspensions to contain higher injected gas pressures is enhanced by aging at higher temperatures.

08 HYDROGEN↗
DOE OSTI2025

A finite-strain rate- and pressure-dependent constitutive framework for analyzing shock compression behavior of cemented tungsten carbides to 100 GPa

In the present study a thermodynamically-consistent finite-strain rate-and-pressure-dependent constitutive framework is implemented to analyze the shock-compression behavior of cemented tungsten carbides to 100 GPa. Central to this framework is the use of logarithmic strain with a set of invariant basis that allow the Cauchy stress tensor to be expressed as a sum of three response terms that are mutually orthogonal, thus permitting a complete separation of the deviatoric and volumetric (pressure) response. An overstress viscoplasticity model that includes strain and strain rate hardening along with thermal softening is used to represent the deviatoric response, while a complete Mie-Grüneisen equation of state (EoS) is used to obtain the pressure response. Using this formulation, the shock-induced compression behavior of cemented tungsten carbide - obtained from planar plate impact experiments using a 30 mm powder gun to peak stresses of up to ~100 GPa - is analyzed to better understand the structure of the measured shock wave profiles and the associated in-material shock quantities. Of particular interest is the evolution of material inelasticity and strength, and temperature in the tungsten carbide samples during the shock compression process.

Cemented tungsten carbide↗
DOE OSTI2025

Geometry and size dependent microstructure and crack formation in Rene 41 superalloy fabricated by laser powder bed fusion

This study provides a systematic investigation into the size and geometry-dependent microstructural evolution and cracking susceptibility of LPBF-fabricated Rene 41. By coupling experimental microstructural analysis with thermal modeling, this research uniquely identifies the relationship between geometry-specific thermal histories, carbide coarsening, and liquation cracking. Four different geometries with varying thickness were fabricated with the identical process parameters. It was found that the grain size and morphology are not affected by the part size. However, the thinner parts exhibited coarser sub-grain structures compared to the thicker ones. The crack formation was observed for the parts with cross-sections smaller than 1 mm, whereas thicker parts had high density without any defects. Further, the cracks were observed in the interdendritic regions, suggesting that liquation cracking was the active micro-crack formation mechanism. The detailed microstructural analysis combined with a thermal finite element analysis proved that the heat extraction efficiency was lower for thinner parts causing a lower cooling rate and coarser carbides, making them more susceptible to constitutional liquation. Microhardness measurements were conducted for each geometry and correlated with the observed microstructural variations. The findings highlight the critical need for geometry-specific optimization of LPBF process parameters to mitigate cracking and achieve microstructural uniformity, offering valuable insights into the fabrication of complex, high-performance aerospace components.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Facet-Dependent Adsorption of Pb(II) on Hematite (001), (116), and (104) Surfaces

Hematite’s common (001) and (012) facets are frequently used in model studies of lead (Pb) adsorption behavior, but there is a lack of research on the high-energy facets, e.g. (104), present in nature. Also, few studies have attempted to connect the molecular details of facet-specific Pb adsorption to macroscopic uptake behavior. Here, to address these knowledge gaps, we investigated Pb(II) adsorption behaviors on facet-engineered hematite nanoparticles dominated by (001), (104), and (116). Adsorption experiments revealed significant variations in Pb(II) uptake among the three samples, with (001) demonstrating the highest capacity and (116) showing the best adsorption efficiency when normalized to specific surface area. Adsorption kinetics followed the pseudo-second-order model, indicating the adsorption process is governed mostly by chemisorption. Adsorption isotherms were well fitted by the Langmuir model, indicating uptake proceeds until roughly monolayer adsorption. Detailed characterization revealed Pb(II) was adsorbed as single atoms with complex inner-sphere binding modes that varied across different facets, indicating adsorption is both structurally and energetically facet-dependent. Co-adsorption experiments further demonstrated Cu 2+ , Zn 2+ , and humic acid significantly promoted Pb(II) adsorption. This study advances the understanding of hematite surface reactivity in controlling macroscopic wet adsorption behaviors, providing valuable insights into the environmental fate of Pb(II).

Contaminant removal↗
DOE OSTI2025

Excitation Laser Energy Dependence of the Gap-Mode TERS Spectra of WS 2 and MoS 2 on Silver

In this work we present a systematic study of the dependence of the gap mode tip-enhanced Raman scattering (TERS) response of the mono- and bi-layer WS 2 and MoS 2 on silver as a function of the excitation laser energy in a broad spectral range from 473nm to 830 nm. For this purpose, we collected consecutive TERS maps of the same area in the sample containing mono-and bi-layer regions with the same TERS probe with 6 different excitation lasers. To decrease the number of collected TERS maps, we used for the first time to the best of our knowledge, concurrent excitation and collection with two lasers simultaneously. We found that the E 2g /A 1g peak intensity ratio for the bilayer WS 2 @Ag and the ratio of the A’/A 1g peak intensity of the out-of-plane mode for the mono- and the bilayer, change in a significantly non-monotonous way as the excitation laser energy is swept from 1.58 eV to 2.62 eV. The former ratio increases at energies corresponding to A and B excitons (~2.0 eV and 2.4 eV correspondingly) in bilayer WS 2 . The absolute intensity of the A’ peak in the monolayer, and correspondingly the A’/A 1g ratio, is surprisingly high at lower excitation energies, but dips dramatically at the energy corresponding to the A exciton, being restored partially in between A and B excitons, but still showing the descending trend as the excitation laser energy increases. A somewhat similar picture was observed in mono- and bi-layers of MoS 2 @Ag, though the existing set of excitation lasers did not match the excitonic profile of this material as nicely as for the case of WS 2 . We attribute the observed behavior to the presence of intermediate (Fano resonance) or strong (Rabi splitting) coupling between the excitons in transition metal dichalcogenides (TMDs) and the plasmons in the tip-substrate nanocavity. This is akin to the so-called Fano (Rabi) transparency experimentally observed in far field scattering from transition metal dichalcogenides between two plasmonic metals. The possibility of the formation of intermediate/strong coupling between the excitonic resonances in TMDs and the nanocavity re-evaluates the role of various resonances in the gap-mode TERS and should become an important factor to be considered by TERS practitioners during planning the experiments. Lastly, based on observed phenomena and its explanation, we propose the “ideal” substrate for efficient TERS and tip enhanced photoluminescence (TEPL) measurements.

Raman spectroscopy↗
DOE OSTI2024

Time-dependent density-functional theory study on nonlocal electron stopping for inertial confinement fusion

Understanding laser–target coupling is of the utmost importance for achieving high performance in laser-direct-drive (LDD) inertial confinement fusion (ICF) experiments. Thus, accurate modeling of electron transport and deposition through ICF-relevant materials and conditions is necessary to quantify the total thermal conduction and ablation. The stopping range is a key transport quantity used in thermal conduction models; in this work, we review the overall role that the electron mean free path (MFP) plays in thermal conduction and hydrodynamic simulations. The currently used modified Lee–More model employs various physics approximations. We discuss a recent model that uses time-dependent density functional theory (TD-DFT) to eliminate these approximations in both the calculation of the electron stopping power and corresponding MFP in conduction zone polystyrene (CH) plasma. In general, the TD-DFT calculations showed a larger MFP (lower stopping power) than the standard modified Lee–More model. Using the TD-DFT results, an analytical model for the electron deposition range, λTD−DFT(ρ,T,K), was devised for CH plasmas between ρ=[0.05−1.05] g/cm3, kBT=[100−1000] eV. We implemented this model into LILAC, for simulations of a National Ignition Facility-scale LDD implosion and compared key physics quantities to ones obtained by simulations using the standard model. The implications of the obtained results and the path moving forward to calculate this same quantity in conduction-zone deuterium–tritium plasmas are further discussed, to hopefully close the understanding gap for laser target coupling in LDD-ICF simulations.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Jet Definition and Transverse-Momentum–Dependent Factorization in Semi-inclusive Deep-Inelastic Scattering

Using the color dipole picture of deep inelastic scattering (DIS) and the color glass condensate effective theory, we study semi-inclusive jet production in DIS at small x in the limit where the photon virtuality Q 2 is much larger than the transverse momentum squared P$^{2}_{⊥}$ of the produced jet. In this limit, the cross section is dominated by aligned jet configurations, that is, quark–antiquark pairs in which one of the fermions—the would-be struck quark in the Breit frame—carries most of the longitudinal momentum of the virtual photon. We show that physically meaningful jet definitions in DIS are such that the effective axis of the jet sourced by the struck quark is controlled by its virtuality rather than by its transverse momentum. For such jet definitions, we show that the next-to-leading order cross section admits factorization in terms of the (sea) quark transverse momentum dependent distribution, which in turn satisfies a universal Dokshitzer-Gribov-Lipatov-Altarelli-Parisi and Sudakov evolution.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Decision-Dependent Uncertainty-Aware Distribution System Planning Under Wildfire Risk

The interaction between power systems and wildfires can be dangerous and costly. Distribution grids can be liable for the outbreak of wildfires during extreme weather. In wildfire-prone areas, investment planning should consider the impact of operational actions on wildfire-related uncertainties affecting line failure likelihood. Here, in this case, endogenous-based uncertainty modeling should comprise the backbone of the investment planning model viz-a-viz the inability of standard exogenous-based uncertainty modeling. Therefore, we propose a decision-dependent uncertainty (DDU) aware methodology to optimize investment portfolios for distribution systems, considering that high power-flow levels in high-threat areas can ignite wildfires and increase line failure probability. The methodology identifies the best combination of upgrades (new lines, hardening existing lines, and placing switching devices). Methodologically, we propose a two-stage distributionally robust planning optimization problem with DDU that considers the distribution system's multiperiod operation. The first stage determines optimal switching actions and line investments, and the second stage evaluates the worst-case expected operational cost under a DDU framework designed to account for the endogenous impact of power-flow levels and hardening investment decisions in the line failure probabilities. An iterative method is tailored to handle the problem and numerical experiments demonstrate a more prepared grid to deal with wildfire risk.

Power systems investment planning↗