Search NASA⌕ Search

SEARCH · Search NASA

Results for “Real Analytic Functions”

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 631 records · Page 35
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Relationships between Atomic Level Surface Structure and Stability/Activity of Platinum Surface Atoms in Aqueous Environments

The development of alternative energy systems for clean production, storage and conversion of energy is strongly dependent on our ability to understand, at atomic-molecular-levels, functional links between activity and stability of electrochemical interfaces. Whereas structure-activity relationships are rapidly evolving, the corresponding structure-stability relationships are still missing. Primarily, this is because there is no adequate experimental approach capable of monitoring in situ stability of well-defined single crystals. Here, by blending the power of Inductively Coupled Plasma-Mass Spectrometer (ICP-MS) connected to a stationary probe to measure in situ and real time dissolution rates of surface atoms (at above 0.4 pg cm-2s-1 levels) and a rotating disk electrode method for monitoring simultaneously the kinetic rates of electrochemical reactions in a single unite, it was possible to establish almost “atom-by-atom” the structure-stability-activity relationships for platinum single crystals in both acidic and alkaline environments. Furthermore, we found that the degree of stability is strongly dependent on the coordination of surface atoms (less coordinated yields less stable), the nature of covalent (adsorption of hydroxyl, oxygen atoms and halides species), and non-covalent interactions (interactions between hydrated Li cations and surface oxide), the thermodynamic driving force for Pt complexation (Pt ion speciation in solution) and the nature of the electrochemical reaction (the oxygen reduction/evolution and CO oxidation reactions). Consequently, these findings are opening new opportunities for elucidating key fundamental descriptors that govern both activity and stability trends, that ultimately, will assist to develop real energy conversion and storage systems.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Simulating Electron Dynamics of Complex Molecules with Time-Dependent Complete Active Space Configuration Interaction

In this paper, time-dependent electronic structure methods are growing in popularity as tools for modeling ultrafast and/or nonlinear processes, for computing spectra, and as the electronic structure component of mean-field molecular dynamics simulations. Time-dependent configuration interaction (TD-CI) offers several advantages over the widely used real-time time-dependent density functional theory: namely, that it correctly models Rabi oscillations; it offers a spin-pure description of open-shell systems; and a hierarchy of TD-CI methods can be defined that systematically approach the exact solution of the time-dependent Schrodinger equation (TDSE). In this work, we present a novel TD-CI approach that extends TD-CI to large complete active-space configuration expansions. Such extension is enabled by use of a direct configuration interaction approach that eliminates the need to explicitly build, store, or diagonalize the Hamiltonian matrix. Graphics processing unit (GPU) acceleration enables fast solution of the TDSE even for large active spaces—up to 12 electrons in 12 orbitals (853776 determinants) in this work. A symplectic split operator propagator yields long-time norm conservation. We demonstrate the applicability of our approach by computing the response of a large molecule with a strongly correlated ground state, decacene (C 42 H 24 ), to various pulses (δ-function, transform limited, chirped). Our simulations predict that chirped pulses can be used to induce dipole-forbidden transitions. Simulations of decacene using the 6-31G(d) basis set and a 12 electrons/12 orbitals active space took 20.1 h to propagate for 100 fs with a 1 attosecond time step on a single NVIDIA K40 GPU. Convergence with respect to time step is found to depend on the property being computed and the chosen active space.

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

Semi-analytic models for the CANDELS survey: comparison of predictions for intrinsic galaxy properties

We compare the predictions of three independently developed semi-analytic galaxy formation models (SAMs) that are being used to aid in the interpretation of results from the CANDELS survey. These models are each applied to the same set of halo merger trees extracted from the 'Bolshoi' high-resolution cosmological N-body simulation and are carefully tuned to match the local galaxy stellar mass function using the powerful method of Bayesian Inference coupled with Markov Chain Monte Carlo or by hand. The comparisons reveal that in spite of the significantly different parameterizations for star formation and feedback processes, the three models yield qualitatively similar predictions for the assembly histories of galaxy stellar mass and star formation over cosmic time. Comparing SAM predictions with existing estimates of the stellar mass function from z = 0-8, we show that the SAMs generally require strong outflows to suppress star formation in low-mass halos to match the present-day stellar mass function, as is the present common wisdom. However, all of the models considered produce predictions for the star formation rates (SFRs) and metallicities of low-mass galaxies that are inconsistent with existing data. The predictions for metallicity-stellar mass relations and their evolution clearly diverge between the models. We suggest that large differences in the metallicity relations and small differences in the stellar mass assembly histories of model galaxies stem from different assumptions for the outflow mass-loading factor produced by feedback. Importantly, while more accurate observational measurements for stellar mass, SFR and metallicity of galaxies at 1 < z < 5 will discriminate between models, the discrepancies between the constrained models and existing data of these observables have already revealed challenging problems in understanding star formation and its feedback in galaxy formation. The three sets of models are being used to construct catalogs of mock galaxies on light cones that have the same geometry as the CANDELS survey, which should be particularly useful for quantifying the biases and uncertainties on measurements and inferences from the real observations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Generating mock galaxy catalogues for flux-limited samples like the DESI Bright Galaxy Survey

ABSTRACT Accurate mock galaxy catalogues are crucial to validate analysis pipelines used to constrain dark energy models. We present a fast HOD-fitting method which we apply to the AbacusSummit simulations to create a set of mock catalogues for the DESI Bright Galaxy Survey, which contain r-band magnitudes and $(g-r)$ colours. The halo tabulation method fits HODs for different absolute magnitude threshold samples simultaneously, preventing unphysical HOD crossing between samples. We validate the HOD fitting procedure by fitting to real-space clustering measurements and galaxy number densities from the MXXL BGS mock, which was tuned to the SDSS and GAMA surveys. The best-fitting clustering measurements and number densities are mostly within the assumed errors, but the clustering for the faint samples is low on large scales. The best-fitting HOD parameters are robust when fitting to simulations with different realizations of the initial conditions. When varying the cosmology, trends are seen as a function of each cosmological parameter. We use the best-fitting HOD parameters to create cubic box and cut sky mocks from the AbacusSummit simulations, in a range of cosmologies. As an illustration, we compare the ${}^{0.1}M_r\lt -20$ sample of galaxies in the mock with BGS measurements from the DESI one-percent survey. We find good agreement in the number densities, and the projected correlation function is reasonable, with differences that can be improved in the future by fitting directly to BGS clustering measurements. The cubic box and cut-sky mocks in different cosmologies are made publicly available.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1999

Sonic Boom Prediction and Minimization of the Douglas Reference OPT5 Configuration

Conventional CFD methods and grids do not yield adequate resolution of the complex shock flow pattern generated by a real aircraft geometry. As a result, a unique grid topology and supersonic flow solver was developed at Northrop Grumman based on the characteristic behavior of supersonic wave patterns emanating from the aircraft. Using this approach, it was possible to compute flow fields with adequate resolution several body lengths below the aircraft. In this region, three-dimensional effects are diminished and conventional two-dimensional modified linear theory (MLT) can be applied to estimate ground pressure signatures or sonic booms. To accommodate real aircraft geometries and alleviate the burdensome grid generation task, an implicit marching multi-block, multi-grid finite-volume Euler code was developed as the basis for the sonic boom prediction methodology. The Thomas two-dimensional extrapolation method is built into the Euler code so that ground signatures can be obtained quickly and efficiently with minimum computational effort suitable to the aircraft design environment. The loudness levels of these signatures can then be determined using a NASA generated noise code. Since the Euler code is a three-dimensional flow field solver, the complete circumferential region below the aircraft is computed. The extrapolation of all this field data from a cylinder of constant radius leads to the definition of the entire boom corridor occurring directly below and off to the side of the aircraft's flight path yielding an estimate for the entire noise "annoyance" corridor in miles as well as its magnitude. An automated multidisciplinary sonic boom design optimization software system was developed during the latter part of HSR Phase 1. Using this system, it was found that sonic boom signatures could be reduced through optimization of a variety of geometric aircraft parameters. This system uses a gradient based nonlinear optimizer as the driver in conjunction with a computationally efficient Euler CFD solver (NIIM3DSB) for computing the three-dimensional near-field characteristics of the aircraft. The intent of the design system is to identify and optimize geometric design variables that have a beneficial impact on the ground sonic boom. The system uses a simple wave drag data format to specify the aircraft geometry. The geometry is internally enhanced and analytic methods are used to generate marching grids suitable for the multi-block Euler solver. The Thomas extrapolation method is integrated into this system, and hence, the aircraft's centerline ground sonic boom signature is also automatically computed for a specified cruise altitude and yields the parameters necessary to evaluate the design function. The entire design system has been automated since the gradient based optimization software requires many flow analyses in order to obtain the required sensitivity derivatives for each design variable in order to converge on an optimal solution. Hence, once the problem is defined which includes defining the objective function and geometric and aerodynamic constraints, the system will automatically regenerate the perturbed geometry, the necessary grids, the Euler solution, and finally the ground sonic boom signature at the request of the optimizer.

Siclari, Michael J.↗
DOE OSTI2025

Simulation of 24,000 Electron Dynamics: Real-Time Time-Dependent Density Functional Theory (TDDFT) with the Real-Space Multigrids (RMG)

Here, we present the theory, implementation, and benchmarking of a real-time time-dependent density functional theory (RT-TDDFT) module within the RMG code, designed to simulate the electronic response of molecular systems to external perturbations. Our method offers insights into nonequilibrium dynamics and excited states across a diverse range of systems, from small organic molecules to large metallic nanoparticles. Benchmarking results demonstrate excellent agreement with established TDDFT implementations and showcase the superior stability of our time integration algorithm, enabling long-term simulations with minimal energy drift. The scalability and efficiency of RMG on massively parallel architectures allow for simulations of complex systems, such as plasmonic nanoparticles with thousands of atoms. Future extensions, including nuclear and spin dynamics, will broaden the applicability of this RT-TDDFT implementation, providing a powerful toolset for studies of photoactive materials, nanoscale devices, and other systems where real-time electronic dynamics is essential.

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

Characterizing Fluorocarbon Assisted Atomic Layer Etching of Si Using Cyclic Ar/C 4 F 8 and Ar/CHF 3 Plasma

With the increasing interest in establishing directional etching methods capable of atomic scale resolution for fabricating highly scaled electronic devices, the need for development and characterization of atomic layer etching (ALE) processes, or generally etch processes with atomic layer precision, is growing. In this work, a flux-controlled cyclic plasma process is used for etching of SiO 2 and Si at the Angstrom-level. This is based on steady-state Ar plasma, with periodic, precise injection of a fluorocarbon (FC) precursor (C 4 F 8 and CHF 3 ), and synchronized, plasma-based Ar+ ion bombardment [D. Metzler et al., J Vac Sci Technol A 32, 020603 (2014), and D. Metzler et al., J Vac Sci Technol A 34, 01B101 (2016)]. For low energy Ar+ ion bombardment conditions, physical sputter rates are minimized, whereas material can be etched when FC reactants are present at the surface. This cyclic approach offers a large parameter space for process optimization. Etch depth per cycle, removal rates, and self-limitation of removal, along with material dependence of these aspects, were examined as a function of FC surface coverage, ion energy, and etch step length using in situ real time ellipsometry. The deposited FC thickness per cycle is found to have a strong impact on etch depth per cycle of SiO 2 and Si, but is limited with regard to control over material etching selectivity. Ion energy over the 20 to 30 eV range strongly impacts material selectivity. The choice of precursor can have a significant impact on the surface chemistry and chemically enhanced etching. CHF 3 has a lower FC deposition yield for both SiO 2 and Si, and also exhibits a strong substrate dependence of FC deposition yield, in contrast to C4F 8 . The thickness of deposited FC layers using CHF 3 is found to be greater for Si than for SiO 2 . X-ray photoelectron spectroscopy was used to study surface chemistry. When thicker FC films of 11 Å are employed, strong changes of FC film chemistry during a cycle are seen whereas the chemical state of the substrate varies much less. On the other hand, for FC film deposition of 5 Å for each cycle, strong substrate surface chemical changes are seen during an etching cycle. The nature of this cyclic etching with periodic deposition of thin FC films differs significantly from conventional etching with steady-state FC layers since surface conditions change strongly throughout each cycle.

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

Stochastic density functional theory: Real- and energy-space fragmentation for noise reduction

Stochastic density functional theory (sDFT) is becoming a valuable tool for studying ground-state properties of extended materials. The computational complexity of describing the Kohn-Sham orbitals is replaced by introducing a set of random (stochastic) orbitals leading to linear and often sub-linear scaling of certain ground-state observables at the account of introducing a statistical error. Schemes to reduce the noise are essential, for example, for determining the structure using the forces obtained from sDFT. Recently, we have introduced two embedding schemes to mitigate the statistical fluctuations in the electron density and resultant forces on the nuclei. Both techniques were based on fragmenting the system either in real space or slicing the occupied space into energy windows, allowing for a significant reduction in the statistical fluctuations. For chemical accuracy, further reduction of the noise is required, which could be achieved by increasing the number of stochastic orbitals. However, the convergence is relatively slow as the statistical error scales as 1 / N &#x3C7; according to the central limit theorem, where N χ is the number of random orbitals. In this paper, we combined the embedding schemes mentioned above and introduced a new approach that builds on overlapped fragments and energy windows. The new approach significantly lowers the noise for ground-state properties, such as the electron density, total energy, and forces on the nuclei, as demonstrated for a G-center in bulk silicon.

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

Predicting Pulsed-Laser Deposition SrTiO 3 Homoepitaxy Growth Dynamics Using High-Speed Reflection High-Energy Electron Diffraction

Pulsed-laser deposition (PLD) is a powerful technique for growing complex oxides with controlled stoichiometry. To understand growth dynamics therein, it is common to leverage in situ spectroscopies, such as reflection high-energy electron diffraction (RHEED), to monitor surface crystallinity. Most commercial systems rely on video-rate cameras operating at 60-120 Hz that lack sufficient temporal resolution to capture growth dynamics at practical deposition frequencies. Here, a high-speed platform to record in situ dynamics via RHEED at >500 Hz is implemented. An open-source analysis package is designed to fit diffraction spots to 2D Gaussians, allowing single-pulse surface reconstruction kinetics extraction. Using homoepitaxially deposited (001)-oriented SrTiO 3 as a model system, we demonstrate how high-speed RHEED can provide real-time insight into growth processes obscured by slower acquisition systems. By fitting the single-pulse intensity to a set of exponential functions, we observe changes in the characteristic decay time and mechanism correlated to the substrate step width and surface termination. We observe distinct surface effects, with diffraction intensity decaying on lower-energy TiO 2 -terminated surfaces and stabilizing on SrO- or mixed-terminated surfaces. Similarly, using an exponential model, the extracted characteristic time of adatom deposition decreases with increased density of bonding sites associated with mixed termination and narrower step widths. Ultimately, this work shows how increasing RHEED temporal resolution can uncover new insights into growth processes, with practical implications for the design and control of PLD processes. This experimental platform provides new capabilities to enable data-driven machine learning analysis and autonomous control systems to enhance the complexity and fecundity of PLD.

(SrO)↗
DOE OSTI2025

Autonomous Flow Electrochemistry for Accelerated Catalyst Discovery

Our objective is to develop an Autonomous Chemical Experimentation (ACE) platform that accelerates discovery of new catalytic transformations and other energy-relevant chemical reactions and processes. We intentionally designed ACE to be highly modular, both with respect to its rapid deployment to different chemistries and experimental workflows as well as incorporation of a wide range of different AI algorithms. In addition to the development of the core software architecture, initial efforts were made to incorporate Large Language Models to provide human-interpretable reasoning of the optimizer’s actions, and to develop a user-friendly graphical interface for experimental researchers. ACE was demonstrated using a flow electrocatalysis platform containing an inline FTIR spectrometer for real-time analysis and quantification of the reaction outcome. Human-in-the-loop experiments were performed in which a human researcher conducted an experiment using electrode potentials suggested by ACE, then fed the spectral data back to ACE for decision making. After confirming the successful function of the optimizer, efforts were next directed to automation of the hardware and performed full autonomy tests using three reactions: catalytic oxidation of formate, catalytic oxidation of cyclohexanol, and oxidation of hydroquinone. These studies confirm that ACE can close the loop between reaction execution, analysis, and optimization. They also reveal that more improved product detection methods will be essential for ACE to make well-informed decisions for reactions with low conversions.

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

Development and Optimization of Normal Pulse Voltammetry for Real Time Analysis of Electrorefiner Salt with High Concentrations of MgCl{sub 2} and UCl{sub 3}

Electro-refiners can accumulate large holdup of fissile actinides from used nuclear fuel including plutonium and will need to be well safeguarded in a commercial facility. Idaho National Laboratory has operated the Fuel Conditioning Facility, an engineering scale demonstration facility, with two electro-refiners for two decades using mass tracking models supplemented with destructive analysis for the requirements of domestic safeguards. This approach is known to be insufficient of meet international safeguards timeliness requirements. Since there are non-nuclear weapons states that are interested in commercial electro-refiners, there is a need to develop a real time concentration measurement technique. In addition to being useful for safeguards, such a real time measurement would be functional for maintaining optimal process control. Normal pulse voltammetry (NPV) is a powerful electrochemical analysis method that has been identified to have great potential for measuring concentrations of various ions in a given salt mixture. While cyclic voltammetry (CV), chronoamperometry (CA), and chronopotentiometry (CP) result in progressive growth of the working electrode surface area, NPV potentially can circumvent this problem via periodically stripping deposited metals from the electrode (WE) surface. With each successive pulse, the potential is dropped to a more negative level. The corresponding current is a summation of partial currents from each of the reducing species. The ultimate objective of this work is to develop correlations between limiting currents measured via NPV with known concentrations. This is potentially a more accurate method than CV and other methods, because it should keep the WE surface area from growing significantly. Previous work on using NPV for actinide concentration determination in eutectic LiCl-KCl has been reported by M. Iizuka et al. However, only low uranium concentrations were studied in that paper. The maximum concentration of uranium tested was 2.48 wt%. Limiting currents from NPV were found to be proportional to their concentrations up to about 1.7 wt%. They were unable to maintain double plateaus (U and Pu) at uranium concentrations above 1.69 w%. This issue was resolved during our high concentration study. This result represents a significant breakthrough for optimizing voltammetry for analysis of actinides in molten salt.

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

Building the Scale Model TELBOT Master Arm - 19585

For centuries, the art and science of modeling have gone hand in hand with the advancements of technology. Among the most impressive examples of the past modeling were ship models constructed for the British Royal Navy. Referred to as Admiralty models, they were principally constructed during the 18. and 19. centuries to depict proposed warship design. To this day, models are important in the field of engineering, where analytical modeling of a new design needs to be verified. In our project the old practice of building models will be fortified by the most advanced software and 3-D printing technology and used to manufacture a fully working Telbot Master Arm scaled down to 0.4-0.5 of its original size. The Telbot is a robotic arm used in different fields like the nuclear and chemical industry, often in combination with other products of our company. The main features of the Telbot's sleek construction are: 1) the high load capacity design; 2) unlimited rotation of all axes; 3) no external or internal wiring of the arm; 4) unlimited, fast and precise arm movements; 5) smooth and sealed arm bodies and joints; 6) high repeat accuracy; 7) resistance to radiation, and it can handle precise loads from 5 kg to 150 kg. Presently, to demonstrate such impressive design features and movement possibilities to its customers, HWM has to build a 3-D virtual reality simulation of the environment where Telbot is proposed to be used and then the arm simulates its movement in this environment. In addition, though HWM does have a showroom for Telbot, due to its weight and size it is not possible to take the full size Telbot on the road to show off its capabilities to customers. The challenge is to build a scaled down full working version of the Telbot robotic arm that can be moved in every direction and angle. The model will be based on the redeveloped T300 Telbot, with 3 joints and 7 free to move pivots. It will be at a 0.4 or 0.5 scale. The outer shell will be 3D printed. The sensors will forward the position and angle to HWM virtual reality software. With these features, any changes made on the model would get displayed in the software one-to-one. The result will be a Telbot Master arm model which can be easily transported; Waellischmiller would be able to present the Telbot not just virtually, but as a real product with all its functions and versatility and show it all over the world. The extensive construction of the virtual 3-D stimulation can be eliminated and the course of movement to reach specific positions could get simulated faster and easier. It can also be used internally and externally for training. This paper will address the challenges, materials used and manufacturing process of building a fully-working scaled-down version of the Telbot Master Arm. (authors)

3D PRINTING↗
DOE OSTI2018

Importance of sulfate radical anion formation and chemistry in heterogeneous OH oxidation of sodium methyl sulfate, the smallest organosulfate

Organosulfates are important organosulfur compounds present in atmospheric particles. While the abundance, composition, and formation mechanisms of organosulfates have been extensively investigated, it remains unclear how they transform and evolve throughout their atmospheric lifetime. To acquire a fundamental understanding of how organosulfates chemically transform in the atmosphere, this study investigates the heterogeneous OH radical-initiated oxidation of sodium methyl sulfate (CH 3 SO 4 Na) droplets, the smallest organosulfate detected in atmospheric particles, using an aerosol flow tube reactor at a high relative humidity (RH) of 85%. Aerosol mass spectra measured by a soft atmospheric pressure ionization source (direct analysis in real time, DART) coupled with a high-resolution mass spectrometer showed that neither functionalization nor fragmentation products are detected. Instead, the ion signal intensity of the bisulfate ion (HSO 4 - ) has been found to increase significantly after OH oxidation. We postulate that sodium methyl sulfate tends to fragment into a formaldehyde (CH 2 O) and a sulfate radical anion (SO 4 .- ) upon OH oxidation. The formaldehyde is likely partitioned back to the gas phase due to its high volatility. The sulfate radical anion, similar to OH radical, can abstract a hydrogen atom from neighboring sodium methyl sulfate to form the bisulfate ion, contributing to the secondary chemistry. Kinetic measurements show that the heterogeneous OH reaction rate constant, k , is (3.79 ± 0.19) × 10 -13 cm 3 molecule -1 s -1 with an effective OH uptake coefficient, γ eff , of 0.17 ± 0.03. While about 40% of sodium methyl sulfate is being oxidized at the maximum OH exposure (1.27 × 10 12 molecule cm -3 s), only a 3% decrease in particle diameter is observed. This can be attributed to a small fraction of particle mass lost via the formation and volatilization of formaldehyde. Overall, we firstly demonstrate that the heterogeneous OH oxidation of an organosulfate can lead to the formation of sulfate radical anion and produce inorganic sulfate. Fragmentation processes and sulfate radical anion chemistry play a key role in determining the compositional evolution of sodium methyl sulfate during heterogeneous OH oxidation.

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

A divide-conquer-recombine algorithmic paradigm for large spatiotemporal quantum molecular dynamics simulations

We introduce an extension of the divide-and-conquer (DC) algorithmic paradigm called divide-conquer-recombine (DCR) to perform large quantum molecular dynamics (QMD) simulations on massively parallel supercomputers, in which interatomic forces are computed quantum mechanically in the framework of density functional theory (DFT). In DCR, the DC phase constructs globally informed, overlapping local-domain solutions, which in the recombine phase are synthesized into a global solution encompassing large spatiotemporal scales. For the DC phase, we design a lean divide-and-conquer (LDC) DFT algorithm, which significantly reduces the prefactor of the O(N) computational cost for N electrons by applying a density-adaptive boundary condition at the peripheries of the DC domains. Our globally scalable and locally efficient solver is based on a hybrid real-reciprocal space approach that combines: (1) a highly scalable real-space multigrid to represent the global charge density; and (2) a numerically efficient plane-wave basis for local electronic wave functions and charge density within each domain. Hybrid space-band decomposition is used to implement the LDC-DFT algorithm on parallel computers. A benchmark test on an IBM Blue Gene/Q computer exhibits an isogranular parallel efficiency of 0.984 on 786 432 cores for a 50.3 × 10{sup 6}-atom SiC system. As a test of production runs, LDC-DFT-based QMD simulation involving 16 661 atoms is performed on the Blue Gene/Q to study on-demand production of hydrogen gas from water using LiAl alloy particles. As an example of the recombine phase, LDC-DFT electronic structures are used as a basis set to describe global photoexcitation dynamics with nonadiabatic QMD (NAQMD) and kinetic Monte Carlo (KMC) methods. The NAQMD simulations are based on the linear response time-dependent density functional theory to describe electronic excited states and a surface-hopping approach to describe transitions between the excited states. A series of techniques are employed for efficiently calculating the long-range exact exchange correction and excited-state forces. The NAQMD trajectories are analyzed to extract the rates of various excitonic processes, which are then used in KMC simulation to study the dynamics of the global exciton flow network. This has allowed the study of large-scale photoexcitation dynamics in 6400-atom amorphous molecular solid, reaching the experimental time scales.

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

Ab-Initio Investigation of Finite Size Effects in Rutile Titania Nanoparticles with Semilocal and Nonlocal Density Functionals

In this work, we employ hybrid and generalized gradient approximation (GGA) level density functional theory (DFT) calculations to investigate the convergence of surface properties and electronic gap of rutile titania nanoparticles with particle size. The surface energies and electronic gaps are calculated for cuboidal particles with minimum dimension ranging from 3.7 Angstrom (24 atoms) to 10.3 Angstrom (384 atoms) using a highly-parallel real-space DFT code to enable hybrid level DFT calculations of larger nanoparticles than are typically practical. We deconvolute the geometric and electronic finite size effects in surface energy, and evaluate the influence of defects on electronic gap and density of states (DOS). The electronic finite size effects in surface energy vanish when the minimum length scale of the nanoparticles becomes greater than 10 Angstrom. We show that this length scale is consistent with a computationally efficient numerical analysis of the characteristic length scale of electronic interactions. The surface energy of nanoparticles having minimum dimension beyond this characteristic length can be approximated using slab calculations that account for the geometric defects. In contrast, the finite size effects on the electronic gap and DOS is highly dependent on the shape and size of these particles. Furthermore, the DOS for cuboidal particles and more realistic particles constructed using the Wulff algorithm reveal that defect states within the electronic gap play a key role in determining the eigen value distribution of nanoparticles and the electronic gap does not converge to the bulk limit for the particle sizes investigated.

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

Two-electron Rabi oscillations in real-time time-dependent density-functional theory

We investigate the Rabi oscillations of electrons excited by an applied electric field in several simple molecular systems using time-dependent configuration interaction (TDCI) and real-time time-dependent density-functional theory (RT-TDDFT) dynamics. While the TDCI simulations exhibit the expected single-electron Rabi oscillations at a single resonant electric field frequency, Rabi oscillations in the RT-TDDFT simulations are a two-electron process. The existence of two-electron Rabi oscillations is determined both by full population inversion between field-free molecular orbitals and the behavior of the instantaneous dipole moment during the simulations. Furthermore, the Rabi oscillations in RT-TDDFT are subject to an intensity threshold of the electric field, below which Rabi oscillations do not occur and above which the two-electron Rabi oscillations occur at a broad range of frequencies. It is also shown that at field intensities near the threshold intensity, the field frequency predicted to induce Rabi oscillations by linear response TDDFT only produces detuned Rabi oscillations. Instead, the field frequency that yields the full two-electron population inversion and Rabi oscillation behavior is shown to be the average of single-electron transition frequencies from the ground S{sub 0} state and the doubly-excited S{sub 2} state. The behavior of the two-electron Rabi oscillations is rationalized via two possible models. The first model is a multi-photon process that results from the electric field interacting with the three level system such that three level Rabi oscillations may occur. The second model suggests that the mean-field nature of RT-TDDFT induces paired electron propagation.

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

Direct observation of ring-opening dynamics in strong-field ionized selenophene using femtosecond inner-shell absorption spectroscopy

Femtosecond extreme ultraviolet transient absorption spectroscopy is used to explore strong-field ionization induced dynamics in selenophene (C 4 H 4 Se). The dynamics are monitored in real-time from the viewpoint of the Se atom by recording the temporal evolution of element-specific spectral features near the Se 3d inner-shell absorption edge (~58 eV). The interpretation of the experimental results is supported by first-principles time-dependent density functional theory calculations. The experiments simultaneously capture the instantaneous population of stable molecular ions, the emergence and decay of excited cation states, and the appearance of atomic fragments. The experiments reveal, in particular, insight into the strong-field induced ring-opening dynamics in the selenophene cation, which are traced by the emergence of non-cyclic molecules as well as the liberation of Se + ions within an overall time scale of approximately 170 fs. In this study, we propose that both products may be associated with dynamics on the same electronic surfaces but with different degrees of vibrational excitation. The time-dependent inner-shell absorption features provide direct evidence for a complex relaxation mechanism that may be approximated by a two-step model, whereby the initially prepared, excited cyclic cation decays within τ 1 = 80 ± 30 fs into a transient molecular species, which then gives rise to the emergence of bare Se + and ring-open cations within an additional τ 2 = 80 ± 30 fs. The combined experimental and theoretical results suggest a close relationship between σ* excited cation states and the observed ring-opening reactions. In conclusion, the findings demonstrate that the combination of femtosecond time-resolved core-level spectroscopy with ab initio estimates of spectroscopic signatures provide new insights into complex, ultrafast photochemical reactions such as ring-opening dynamics in organic molecules in real-time and with simultaneous sensitivity for electronic and structural rearrangements.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗