Search NASA⌕ Search

SEARCH · Search NASA

Results for “Work Function”

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,657 records · Page 92
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Ensuring thermodynamic consistency with invertible coarse-graining

Coarse-grained models are a core computational tool in theoretical chemistry and biophysics. A judicious choice of a coarse-grained model can yield physical insights by isolating the essential degrees of freedom that dictate the thermodynamic properties of a complex, condensed-phase system. The reduced complexity of the model typically leads to lower computational costs and more efficient sampling compared with atomistic models. Designing “good” coarse-grained models is an art. Generally, the mapping from fine-grained configurations to coarse-grained configurations itself is not optimized in any way; instead, the energy function associated with the mapped configurations is. In this work, we explore the consequences of optimizing the coarse-grained representation alongside its potential energy function. We use a graph machine learning framework to embed atomic configurations into a low-dimensional space to produce efficient representations of the original molecular system. Because the representation we obtain is no longer directly interpretable as a real-space representation of the atomic coordinates, we also introduce an inversion process and an associated thermodynamic consistency relation that allows us to rigorously sample fine-grained configurations conditioned on the coarse-grained sampling. We illustrate that this technique is robust, recovering the first two moments of the distribution of several observables in proteins such as chignolin and alanine dipeptide.

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

NEON AOP Survey of Upper East River CO Watersheds: Waveform LiDAR Binary Data

The waveform Light Detection and Ranging (LiDAR) data in this package were generated through a National Ecological Observatory Network Airborne Observation Platform (NEON AOP) acquisition over watersheds of interest surrounding Crested Butte, Colorado. The remote sensing imagery acquired by the NEON AOP supports an interdisciplinary project on hydrology, biogeochemistry, and ecosystem functioning in a snow-dominated headwater environment. These waveform LiDAR data enable spatially continuous estimation of vegetation structure parameters to facilitate analyses of the major environmental drivers of structural and compositional variability. The package contains 97 compressed file archives in 7-zip (.7z) format, each corresponding to one acquisition flightpath. Within each .7z archive is a set of constituent files describing properties of the LiDAR waveforms, such as return intensity, geolocation, outgoing pulse and other behavior of the sensor and signals. Once downloaded, the files must first be unzipped using the widely distributed command-line software utility 7z, using the command '7z x \[filename\].7z \[target_directory\]'. All files within the .7z archives can be opened in IDL, MatLab, or the open-source R statistical computing environment. Further details about the data package are in the attached user guide (neon_aop_crbu_waveformlidar_userguide.pdf). This work was supported by the Watershed Function Science Focus Area at Lawrence Berkeley National Laboratory funded by the US Department of Energy, Office of Science, Biological and Environmental Research under Contract No. DE-AC02-05CH11231.

2018 NEON and 2025 CHESS Campaigns↗
DOE OSTI2019

Learning nonlinear level sets for dimensionality reduction in function approximation

We developed a Nonlinear Level-set Learning (NLL) method for dimensionality reduction in high-dimensional function approximation with small data. This work is motivated by a variety of design tasks in real-world engineering applications, where practitioners would replace their computationally intensive physical models (e.g., high-resolution fluid simulators) with fast-to-evaluate predictive machine learning models, so as to accelerate the engineering design processes. There are two major challenges in constructing such predictive models: (a) high-dimensional inputs (e.g., many independent design parameters) and (b) small training data, generated by running extremely time-consuming simulations. Thus, reducing the input dimension is critical to alleviate the over-fitting issue caused by data insufficiency. Existing methods, including sliced inverse regression and active subspace approaches, reduce the input dimension by learning a linear coordinate transformation; our main contribution is to extend the transformation approach to a nonlinear regime. Specifically, we exploit reversible networks (RevNets) to learn nonlinear level sets of a high-dimensional function and parameterize its level sets in low-dimensional spaces. A new loss function was designed to utilize samples of the target functions' gradient to encourage the transformed function to be sensitive to only a few transformed coordinates. The NLL approach is demonstrated by applying it to three 2D functions and two 20D functions for showing the improved approximation accuracy with the use of nonlinear transformation, as well as to an 8D composite material design problem for optimizing the buckling-resistance performance of composite shells of rocket inter-stages.

Zhang, Guannan↗
DOE OSTI2017

eQE: An open‐source density functional embedding theory code for the condensed phase

Abstract In this work, we present the main features and algorithmic details of a novel implementation of the frozen density embedding (FDE) formulation of subsystem density functional theory (DFT) that is specifically designed to enable ab initio molecular dynamics (AIMD) simulations of large‐scale condensed‐phase systems containing 1000s of atoms. This code (available at http://eqe.rutgers.edu ) has been given the moniker of embedded Quantum ESPRESSO (eQE) as it is a generalization of the open‐source Quantum ESPRESSO (QE) suite of programs. The strengths of eQE reside in a hierarchical parallelization scheme that allows for an efficient and fully self‐consistent treatment of the electronic structure (via the addition of an additional DIIS extrapolation layer) while simultaneously exploiting the inherent symmetries and periodicities in the system (via sampling of subsystem‐specific first Brillouin zones and utilization of subsystem‐specific basis sets). While bulk liquids and molecular crystals are two classes of systems that exemplify the utility of the FDE approach (as these systems can be partitioned into weakly interacting subunits), we show that eQE has significantly extended this regime of applicability by outperforming standard semilocal Kohn–Sham DFT (KS‐DFT) for large‐scale heterogeneous catalysts with quite different layer‐specific electronic structure and intrinsic periodicities. eQE features very favorable strong parallel scaling for a model system of bulk liquid water composed of 256 water molecules, which allows for a significant decrease in the overall time to solution when compared to KS‐DFT. We show that eQE achieves speedups greater than one order of magnitude ( ) when performing AIMD simulations of such large‐scale condensed‐phase systems as: (1) molecular liquids via bulk liquid water represented by 1024 independent water molecules (3072 atoms with a 25.3× speedup over KS‐DFT), (2) polypeptide/biomolecule solvation via ( gly ) 6 solvated in (H 2 O) 395 (1230 atoms with a 38.6× speedup over KS‐DFT), and (3) molecular crystals via a 3 × 3 × 3 periodic supercell of pentacene (1940 atoms with a 12.0× speedup over KS‐DFT). These results represent a significant improvement over the current state‐of‐the‐art and now enable subsystem DFT‐based AIMD simulations of realistically sized condensed‐phase systems of interest throughout chemistry, physics, and materials science.

Genova, Alessandro↗
DOE OSTI2018

Isoprenylcysteine carboxylmethyltransferase regulates ovarian cancer cell response to chemotherapy and Ras activation

Highlights: • Icmt upregulation is a common feature in ovarian cancer. • Icmt expression level varies in ovarian cancer tissues and cell lines. • Icmt expression is associated with ovarian cancer cell chemoresistance. • Icmt inhibition sensitizes ovarian cancer to chemotherapy in vitro and in vivo. • Icmt inhibition suppresses Ras and its downstream signaling pathways. Inhibition of isoprenylcysteine carboxylmethyltransferase (Icmt), which catalyzes the final step of oncoproteins' prenylation, targets growth and survival of various cancers. In this work, we systematically studied the expression, functions and molecular signaling of Icmt in ovarian cancer. We show that the upregulation of Icmt expression is a common feature in patients with epithelial ovarian cancer regardless of age and disease stage. In line with the observations in ovarian cancer patients, a panel of epithelial ovarian cancer cell lines also demonstrates the significant increase on Icmt transcript and protein levels than normal ovarian epithelial cells. In addition, ovarian cancer cell lines with higher Icmt levels are more resistant to chemotherapeutic agents. We further show that Icmt inhibition by siRNA or inhibitor cysmethynil suppresses growth and induces apoptosis in ovarian cancer cells. Importantly, Icmt inhibition significantly augments chemotherapeutic agent's efficacy in vitro and in vivo, demonstrating the translational potential of Icmt inhibition in ovarian cancer. Mechanistically, we show that Ras activation is a critical effector of Icmt in ovarian cancer cells. Using cell culturing system, mouse model and patient samples, our work is the first to demonstrate the essential roles of Icmt in ovarian cancer via Ras signaling, particularly on its response to chemotherapy. Our findings suggest that Icmt inhibition is a promising therapeutic strategy to overcome chemoresistance in ovarian cancer, in particular, those patients with high Icmt expression.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2022

Emergent Magnetic States and Tunable Exchange Bias at 3d Nitride Heterointerfaces

Abstract Interfacial magnetism stimulates the discovery of giant magnetoresistance (MR) and spin–orbital coupling across the heterointerfaces, facilitating the intimate correlation between spin transport and complex magnetic structures. Over decades, functional heterointerfaces composed of nitrides have seldom been explored due to the difficulty in synthesizing high‐quality nitride films with correct compositions. Here, the fabrication of single‐crystalline ferromagnetic Fe 3 N thin films with precisely controlled thicknesses is reported. As film thickness decreases, the magnetization dramatically deteriorates, and the electronic state changes from metallic to insulating. Strikingly, the high‐temperature ferromagnetism is maintained in a Fe 3 N layer with a thickness down to 2 u.c. (≈8 Å). The MR exhibits a strong in‐plane anisotropy; meanwhile, the anomalous Hall resistivity reverses its sign when the Fe 3 N layer thickness exceeds 5 u.c. Furthermore, a sizable exchange bias is observed at the interfaces between a ferromagnetic Fe 3 N and an antiferromagnetic CrN. The exchange bias field and saturation moment strongly depend on the controllable bending curvature using the cylinder diameter engineering technique, implying the tunable magnetic states under lattice deformation. This work provides a guideline for exploring functional nitride films and applying their interfacial phenomena for innovative perspectives toward practical applications.

Jin, Qiao↗
DOE OSTI2016

Remote Retrieval, Disruption and Processing of Legacy Nuclear Waste - 16179

The UK National Nuclear Laboratory (NNL) is currently undertaking, in conjunction with Sellafield Ltd, a series of high profile trials of robotic equipment on a full scale test facility at their Workington Laboratory in England. These trials are designed to underpin and substantiate the processes proposed for the remote retrieval, disruption and processing of legacy Miscellaneous Beta Gamma Waste (MBGW) from a number of different donor plants on the Sellafield site, including some of the older fuel storage ponds and storage silos. These processes will be conducted inside the new Box Encapsulation Plant (BEP) which is due to be commissioned at Sellafield by 2020. These processes will subsequently enable Sellafield Ltd to safely recover legacy Intermediate Level Waste (ILW) from the older storage facilities on site and process it into a safer form which is more suitable for long term storage and then ultimately final disposal in a Geological Disposal Facility. As part of this legacy waste processing within the BEP plant, MBGW waste items will have to be remotely removed from a variety of different types of import skips. This waste will then be transferred, using robots, into the export container, which is a 3 m{sup 3} box liner, where the processed waste will be flood grouted to produce an encapsulated waste package. As part of this processing, any MBGW items consisting of sealed containers will have to be remotely disrupted before being grout encapsulated. This disruption is in order to (a) release any trapped liquors or gases from the containers and (b) ensure full grout infiltration of the container contents, thus minimising the amount of free liquors and voids within the final waste product. The NNL are conducting this technical work to confirm the viability of the proposed concept, processes and equipment. The developed solution consists of a linear layout of processing equipment with a pair of heavy duty, standard industrial, 'commercially off the shelf' (COTS) robots positioned around a Waste Handling Table. These robots have a 500 kg payload capacity and are operated remotely to handle all of the various types of MBGW waste, including the heavy and bulky items. The robots can be operated in manual tele-operation mode as well as pre-programmed automatic mode and have a range of primary and secondary tooling in order to perform the waste processing and disruption tasks as well as undertake general house-keeping and cleaning duties within the cell. The main tools include two types of hydraulic grab for the main waste recovery and handling tasks and a hydraulic shear for the waste disruption activities, but a number of other tools have also been designed and built for ancillary activities, such as a set of general purpose gripper jaws, drilling tool, bolt removal tool, pump tool, spray tool and a number of different housekeeping tools etc. The robot trials were started in November 2014, and have attracted a significant degree of high level stakeholder interest, including the UK Nuclear Decommissioning Authority (NDA), Office for Nuclear Regulation (ONR), Environment Agency (EA) and Radioactive Waste Management Limited (RWML) witnessing a number of the critical demonstration trials. Robots of the type selected are readily available, robust, reliable and easy to use. The trials have been very successful in proving the retrieval, processing and disruption operations, reducing the burden in the legacy waste facilities. Utilising robots to perform these tasks will lower the probability of producing non-compliant waste packages. In addition the robots can increase the consolidation of incoming waste skips into outgoing 3 m{sup 3} box liners, significantly increasing packing efficiencies and reducing the overall number of export containers produced thus lowering the decommissioning and waste storage/disposal costs. The trials have also enabled an assessment of the optimum way to operate the robots, as well as their reliability and maintainability issues. Engineering work associated with recovering failed robots and removing them from a typical Sellafield plant environment and then subsequently replacing them with functioning robots has also been undertaken and further work is now planned. The success of the trials, completed safely to a very challenging timescale, has enabled the use of robots in BEP to be endorsed as a viable option by the key stakeholders. Following a robust decision making process, this was a clear cut decision. The technical development and testing of the robot technology is currently continuing in order to investigate areas such as robot reliability, failure recoveries, human factors, ease of waste identification, optimisation of throughput rates, nuclearisation modifications etc. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2017

Quasiparticle and hybrid density functional methods in defect studies: An application to the nitrogen vacancy in GaN

Defects in semiconductors can play a vital role in the performance of electronic devices, with native defects often dominating the electronic properties of the semiconductor. Understanding the relationship between structural defects and electronic function will be central to the design of new high-performance materials. In particular, it is necessary to quantitatively understand the energy and lifetime of electronic states associated with the defect. In this work, we apply first-principles density functional theory (DFT) and many-body perturbation theory within the GW approximation to understand the nature and energy of the defect states associated with a charged nitrogen vacancy on the electronic properties of gallium nitride (GaN), as a model of a well-studied and important wide gap semiconductor grown with defects. We systematically investigate the sources of error associated with the GW approximation and the role of the underlying atomic structure on the predicted defect state energies. Additionally, analysis of the computed electronic density of states (DOS) reveals that there is one occupied defect state 0.2 eV below the valence band maximum and three unoccupied defect states at energy of 0.2–0.4 eV above the conduction band minimum, suggesting that this defect in the +1 charge state will not behave as a carrier trap. Furthermore, we compare the character and energy of the defect state obtained from GW and DFT using the HSE approximate density functional and find excellent agreement. This systematic study provides a more complete understanding of how to obtain quantitative defect energy states in bulk semiconductors.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2015

Ab initio Bogoliubov coupled cluster theory for open-shell nuclei

Background: Ab initio many-body methods have been developed over the past 10 yr to address closed-shell nuclei up to mass A≈130 on the basis of realistic two- and three-nucleon interactions. A current frontier relates to the extension of those many-body methods to the description of open-shell nuclei. Several routes to address open-shell nuclei are currently under investigation, including ideas that exploit spontaneous symmetry breaking. Purpose: Singly open-shell nuclei can be efficiently described via the sole breaking of U(1) gauge symmetry associated with particle-number conservation as a way to account for their superfluid character. While this route was recently followed within the framework of self-consistent Green's function theory, the goal of the present work is to formulate a similar extension within the framework of coupled cluster theory. Methods: We formulate and apply Bogoliubov coupled cluster (BCC) theory, which consists of representing the exact ground-state wave function of the system as the exponential of a quasiparticle excitation cluster operator acting on a Bogoliubov reference state. Equations for the ground-state energy and the cluster amplitudes are derived at the singles and doubles level (BCCSD) both algebraically and diagrammatically. The formalism includes three-nucleon forces at the normal-ordered two-body level. The first BCC code is implemented in m scheme, which will permit the treatment of doubly open-shell nuclei via the further breaking of SU(2) symmetry associated with angular momentum conservation. Results: Proof-of-principle calculations in an N max =6 spherical harmonic oscillator basis for 16,18 O and 18 Ne in the BCCD approximation are in good agreement with standard coupled cluster results with the same chiral two-nucleon interaction, while 20 O and 20 Mg display underbinding relative to experiment. The breaking of U(1) symmetry, monitored by computing the variance associated with the particle-number operator, is relatively constant for all five nuclei, in both the Hartree-Fock-Bogoliubov and BCCD approximations. Conclusions: The newly developed many-body formalism increases the potential span of ab initio calculations based on single-reference coupled cluster techniques tremendously, i.e., potentially to reach several hundred additional midmass nuclei. The new formalism offers a wealth of potential applications and further extensions dedicated to the description of ground and excited states of open-shell nuclei. Short-term goals include the implementation of three-nucleon forces at the normal-ordered two-body level. Midterm extensions include the approximate treatment of triples corrections and the development of the equation-of-motion methodology to treat both excited states and odd nuclei. Long-term extensions include exact restoration of U(1) and SU(2) symmetries.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Historical Contingency in Microbial Resilience to Hydrologic Perturbations

Development of reliable biogeochemical models requires a mechanistic consideration of microbial interactions with hydrology. Microbial response to and its recovery after hydrologic perturbations (i.e., resilience) is a critical component to understand in this regard, but generally difficult to predict because the impacts of future events can be dependent on the history of perturbations (i.e., historical contingency). Fundamental issues underlying this phenomenon include how microbial resilience to hydrologic perturbations is influenced by historical contingency and how their relationships vary depending on the characteristics of microbial functions. To answer these questions, we considered a simple microbial community composed of two species that redundantly consume a common substrate but specialize in producing distinct products and developed a continuous flow reactor model where the two species grow with trade-offs along the flow rate. Simulations of this model revealed that (1) the history of hydrologic perturbations can lead to the shifts in microbial populations, which consequently affect the community's functional dynamics, and (2) while historical contingency in resilience was consistently predicted for all microbial functions, it was more pronounced for specialized functions, compared to the redundant function. As a signature of historical contingency, our model also predicted the emergence of hysteresis in the transitions across conditions, a critical aspect that can affect transient formation of intermediate compounds in biogeochemistry. This work presents microbial growth traits and their functional redundancy or specialization as fundamental factors that control historical contingencies in resilience.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

Soil physical and chemical measurements for topsoils collected during NEON campaign in East River, CO (06/14/2018-06/28/2018)

The package is part of the DOE Watershed Function Science Focus Area (SFA) project and includes soil physical and chemical measurements from topsoils collected at the East River, Colorado, in conjunction with the National Ecological Observatory Network (NEON) Airborne Observation Platform (AOP) survey conducted in June 2018. The soil measurements include soil bulk density, soil volumetric water content, soil microbial biomass C (Carbon), N (Nitrogen) and C:N (C to N ratio), soil DNA yield, soil total extractable organic C, soil total extractable N, soil extractable nitrate, soil extractable ammonium, soil dissolved inorganic N, soil dissolved organic N, soil pH, soil TOC400 (total organic carbon at 400°C), soil ROC (residual oxidizable carbon), soil TIC (total inorganic carbon), soil TOC (total organic carbon), soil TC (total carbon), soil N, soil OM (organic matter) loss on ignition. Additional associated site metadata can be found in the ESS-DIVE package 10.15485/1618130. The dataset includes (1) 2018_NEON_soil_physical_chemical_measurements.csv: soil physical and chemical measurements indexed by soil sample IGSNs; (2) samples.csv: sample metadata file used to register International Generic Sample Numbers (IGSNs); (3) flmd.csv: file level metadata file; and (4) dd.csv: data dictionary file. This work was supported by the Watershed Function Science Focus Area at Lawrence Berkeley National Laboratory funded by the US Department of Energy, Office of Science, Biological and Environmental Research under Contract No. DE-AC02-05CH11231.

2018 NEON and 2025 CHESS (Catchment Hydrology and ↗
DOE OSTI2018

Microstructure and fracture properties of open-cell porous Ti-6Al-4V with high porosity fabricated by electron beam melting

Highlights: • The Ti-6Al-4V foams with 80% high porosity are prepared by electron beam melting. • The ligament width appears to play the major role for mechanical properties. • The work of fracture is a linear function of ligament width and roughness/width ratio. - Abstract: The open-cell porous Ti-6Al-4V structure with high porosity, intended to be applied as replacement for human cancellous bone, were fabricated by electron beam melting (EBM). Computer aided design (CAD) was applied to design porous structures using the same unit cell with different unit cell sizes from 2.5 to 4 mm, different ligament widths from 600 to 900 μm, and different pore sizes from 1200 to 1800 μm, in order to achieve high porosity of 80% in avoiding the stress shielding effect. In comparison with the CAD designs and the EBM samples, there were minor discrepancies in terms of pore size and ligament width, mainly a result of melting pool. The measured data on the Young's modulus and yield strength of the EBM porous samples can be predicted by the Gibson and Ashby model. All samples with high porosity were found to match well with cancellous bone, with Young's modulus of 2 GPa and yield stress of 31 MPa, effective to diminish the risk of stress shielding effect. For porous EBM sample with high 80% porosity, the work of fracture can increase from 15.9 to 47.6 kJ/m{sup 2} with increasing ligament width from 600 to 900 μm.

36 MATERIALS SCIENCE↗
DOE OSTI2026

FAD-Toolset (Floating Array Design Toolset) [SWR-26-056]

The Floating Array Design (FAD) Toolset is a collection of tools for modeling and designing arrays of floating offshore structures. It was originally designed for floating wind systems but has applicability for many offshore applications. A core part of the FAD Toolset is the floating array model, which serves as a high-level library for efficiently modeling a floating array, such as a floating wind array. It combines site condition information and a description of the floating array design, and contains functions for evaluating the array's behavior considering the site conditions. For example, it combines information about site soil conditions, mooring line loads, and an array's anchor characteristics to estimate the holding capacity of each anchor. The library works in conjunction with the tools RAFT, MoorPy, and FLORIS to model floating platforms, wind turbines, mooring systems, power cables, and array wakes respectively. Layered on top of the floating array model is a set of design tools that can be used for algorithmically adjusting or optimizing parts of the a floating array. Specific tools existing for mooring lines, shared mooring systems, dynamic power cables, static power cable routing, and overall array layout. These capabilities work with the design representation and evaluation functions in the floating array model, and they can be applied by users in various combinations to suit different purposes. In addition to standalone uses of the FAD Toolset, a coupling has been made with Ard, (https://github.com/NLRWindSystems/Ard) a sophisticated and flexible wind farm optimization tool. This coupling allows Ard to use certain mooring system capabilities from FAD to perform layout optimization of floating wind farms with Ard's more advanced layout optimization capabilities. The FAD Toolset works with the IEA Wind Task 49 Ontology (https://github.com/IEAWindTask49/Ontology), which provides a standardized format for describing floating wind farm sites and designs. See example use cases in our examples folder (https://github.com/NLRWindSystems/FAD-Toolset/blob/main/examples/README.md) For working with the library, it is important to understand the floating array model structure, which is described more here: https://github.com/NLRWindSystems/FAD-Toolset/blob/main/fad/README.md.

Sirkis, Leah [National Laboratory of the Rockies (↗
DOE OSTI2020

Loschmidt-amplitude wave function spectroscopy and the physics of dynamically driven phase transitions

In this work, we introduce the Loschmidt amplitude as a powerful tool to perform spectroscopy of generic many-body wave functions and use it to interrogate the wave function obtained after ramping the transverse field quantum Ising model through its quantum critical point. Previous results are confirmed and a more complete understanding of the population of defects and of the effects of magnon-magnon interaction or finite-size corrections is obtained. The influence of quantum coherence is clarified.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2019

Photocatalytic Facet Selectivity in BiVO 4 Nanoparticles: Polaron Electronic Structure and Thermodynamic Stability Considerations for Photocatalysis

Selective charge separation among different crystal facets of a semiconductor is an intriguing phenomenon for which there is no firm and full theoretical foundation currently. In this work, we report on a density functional theory + U characterization of band alignment and electron and hole polaron stabilities among the (010), (110), and (011) facets of bismuth vanadate BiVO 4 (BVO). Computation-derived band alignment indicates that the conduction band minima are at nearly the same level among the three facets but that the valence band maxima exhibit a shift. We also modeled electron and hole polarons as localized electrons and holes on vanadium and oxygen, respectively, and determined their relative stabilities from a “bulk” region to a surface region. Calculated stabilities reveal similar stability profiles across the various facets, with electron polarons most stable when localized on subsurface V atoms and hole polarons most stable on surface O atoms. Calculations indicate a small stability preference for electron polarons toward the (011) facet and for hole polarons toward the (110) facet, whereas, experimentally, interfacial reduction is observed to take place selectively on the (010) facet and oxidation on the (110) and (011) facets. Facet selectivity could be occurring on the basis of thermodynamics (electron or holes showing a stronger affinity for some facets over others) or kinetics (electron or hole transport and/or redox processes being more efficient toward/on some facets over others) or a combination of both. This work establishes that thermodynamic stability alone is not responsible for the observed facet selectivity in BVO. Furthermore, we surmise that polaron transport kinetics and interfacial redox kinetics are likely to have a role in facet selectivity in BVO. These issues will be the subject of future publications.

14 SOLAR ENERGY↗
DOE OSTI2015

Non-linear non-local molecular electrodynamics with nano-optical fields

The interaction of optical fields sculpted on the nano-scale with matter may not be described by the dipole approximation since the fields may vary appreciably across the molecular length scale. Rather than incrementally adding higher multipoles, it is advantageous and more physically transparent to describe the optical process using non-local response functions that intrinsically include all multi-poles. In this work, we present a semi-classical approach for calculating non-local response functions based on the minimal coupling Hamiltonian. The first, second, and third order response functions are expressed in terms of correlation functions of the charge and the current densities. Finally, this approach is based on the gauge invariant current rather than the polarization, and on the vector potential rather than the electric and magnetic fields.

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

On extending Kohn-Sham density functionals to systems with fractional number of electrons

Here, we analyze four ways of formulating the Kohn-Sham (KS) density functionals with a fractional number of electrons, through extending the constrained search space from the Kohn-Sham and the generalized Kohn-Sham (GKS) non-interacting $\mathcal{v}$-representable density domain for integer systems to four different sets of densities for fractional systems. In particular, these density sets are (I) ensemble interacting N -representable densities, (II) ensemble non-interacting N -representable densities, (III) non-interacting densities by the Janak construction, and (IV) non-interacting densities whose composing orbitals satisfy the Aufbau occupation principle. By proving the equivalence of the underlying first order reduced density matrices associated with these densities, we show that sets (I), (II), and (III) are equivalent, and all reduce to the Janak construction. Moreover, for functionals with the ensemble v-representable assumption at the minimizer, (III) reduces to (IV) and thus justifies the previous use of the Aufbau protocol within the (G)KS framework in the study of the ground state of fractional electron systems, as defined in the grand canonical ensemble at zero temperature. By further analyzing the Aufbau solution for different density functional approximations (DFAs) in the (G)KS scheme, we rigorously prove that there can be one and only one fractional occupation for the Hartree Fock functional, while there can be multiple fractional occupations for general DFAs in the presence of degeneracy. This has been confirmed by numerical calculations using the local density approximation as a representative of general DFAs. This work thus clarifies important issues on density functional theory calculations for fractional electron systems.

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

Action-angle formulation of generalized, orbit-based, fast-ion diagnostic weight functions

Due to the usually complicated and anisotropic nature of the fast-ion distribution function, diagnostic velocity-space weight functions, which indicate the sensitivity of a diagnostic to different fast-ion velocities, are used to facilitate the analysis of experimental data. Additionally, when velocity-space weight functions are discretized, a linear equation relating the fast-ion density and the expected diagnostic signal is formed. In a technique known as velocity-space tomography, many measurements can be combined to create an ill-conditioned system of linear equations that can be solved using various computational methods. However, when velocity-space weight functions (which by definition ignore spatial dependencies) are used, velocity-space tomography is restricted, both by the accuracy of its forward model and also by the availability of spatially overlapping diagnostic measurements. In this work, we extend velocity-space weight functions to a full 6D generalized coordinate system and then show how to reduce them to a 3D orbit-space without loss of generality using an action-angle formulation. Moreover, we show how diagnostic orbit-weight functions can be used to infer the full fast-ion distribution function, i.e., orbit tomography. In depth derivations of orbit weight functions for the neutron, neutral particle analyzer, and fast-ion D-α diagnostics are also shown

74 ATOMIC AND MOLECULAR PHYSICS↗