Search NASA⌕ Search

SEARCH · Search NASA

Results for “gradient structure”

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

Gradient-based optimization of complex nanoparticle heterostructures enabled by deep learning on heterogeneous graphs

Applications of deep learning (DL) to design nanomaterials are hampered by a lack of suitable data representations and training data. Here, in this study, we report efforts to overcome these limitations and leverage DL to optimize the nonlinear optical properties of core–shell upconverting nanoparticles (UCNPs). UCNPs, which have applications in fields such as biosensing, super-resolution microscopy and three-dimensional printing, can emit visible and ultraviolet light from near-infrared excitations. We report a large-scale dataset of UCNP emission spectra based on accurate but expensive kinetic Monte Carlo simulations (N > 6,000) and use these data to train a heterogeneous graph neural network using a physically motivated representation of UCNP nanostructure. Applying gradient-based optimization on the trained graph neural network, we identify structures with 6.5× higher predicted emission under 800-nm illumination than any UCNP in our training set. Our work reveals design principles for UCNP heterostructures and presents a roadmap for DL-based inverse design of nanomaterials.

Sivonxay, Eric [Lawrence Berkeley National Laborat↗
DOE OSTI2021

The effects of laser polarization and wavelength on injection dynamics of a laser wakefield accelerator

Here, we investigate the effects of laser polarization and wavelength on electron injection dynamics in a laser wakefield accelerator. During the ionization process, electrons gain residual momentum and kinetic energy via above threshold ionization, which has a strong dependence on laser polarization. A circularly polarized laser pulse results in a much higher residual momentum and kinetic energy gain for the ionized electrons compared with the linearly polarized case. This residual momentum results in particle injection because of the sensitivity of particle trapping to the initial conditions and enhanced the total injected beam charge in both experiments and particle-in-cell simulations. Due to the strong correlation of above threshold ionization with laser wavelength, in this work we extended the investigation to long wavelength (up to 20 μm) drive pulses using particle-in-cell simulations. Owing to the gain in kinetic energy, it may be expected that the charge trapped would consistently increase for circular polarization with increasing laser wavelength, but this was not observed. Instead, there are oscillations with wavelength in the relative trapped charge between linear and circular polarization cases, which arise because of ionization and heating effects on the plasma. Our studies highlight the complex interplay between several different physical effects, including injection regimes—above threshold ionization assisted injection, wave-breaking injection by carrier-envelope-phase effects and ionization injection—ionization gradient induced laser pulse evolution, and thermal modifications to the wake structure that need considering when extrapolating laser wakefield acceleration to different wavelength regimes.

Ma, Y. (ORCID:0000000239194856)↗
DOE OSTI2024

Investigating instabilities in magnetized low-pressure capacitively coupled RF plasma using particle-in-cell (PIC) simulations

The effect of a uniform magnetic field on particle transport in low-pressure radio frequency (RF) capacitively coupled plasma (CCP) has been studied using a particle-in-cell model. Three distinct regimes of plasma behavior can be identified as a function of the magnetic field. In the first regime at low magnetic fields, asymmetric plasma profiles are observed within the CCP chamber due to the effect of $\overrightarrow{E}$ × $\overrightarrow{B}$ drift. As the magnetic field increases, instabilities develop and form self-organized spoke-shaped structures that are distinctly seen within the bulk plasma closer to the sheath. In this second regime, the spoke-shaped coherent structures rotate inside the plasma chamber in the -$\overrightarrow{E}$ × $\overrightarrow{B}$ direction, where $\overrightarrow{E}$ × $\overrightarrow{B}$ are the DC electric and magnetic field vectors, respectively, and the DC electric field exists in the sheath and pre-sheath regions. The spoke rotation frequency is in the megahertz range. As the magnetic field strength increases further, the rotating coherent spokes continue to exist near the sheath. The coherent structures are, however, accompanied by new small-scale incoherent structures originating and moving within the bulk plasma region away from the sheath. This is the third regime of plasma behavior. The threshold values of the magnetic field between these regimes were found not to vary with changing plasma reactor geometry (e.g., area ratio between ground and powered electrodes) or the use of an external capacitor between the RF-powered electrode and the RF source. The threshold values of the magnetic field between these regimes shift toward higher values with increasing gas pressure. Analysis of the results indicates that the rotating structures are due to the lower hybrid instability driven by density gradients and electron-neutral collisions. This paper provides guidance on the upper limit of the magnetic field for instability-free operation in low-pressure CCP-based semiconductor deposition and etch systems that use the external magnetic field for plasma uniformity control.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Equatorial Electronic Structure in the Uranyl Ion: Cs 2 UO 2 Cl 4 and Cs 2 UO 2 Br 4

Electric field gradient (EFG) tensors in the equatorial plane of the linear UO 2 2+ ion have been measured by nuclear magnetic resonance (NMR) and nuclear quadrupole resonance (NQR) experiments and computed by relativistic Kohn–Sham methods with and without environment embedding for Cs 2 UO 2 Cl 4 and Cs 2 UO 2 Br 4 . This approach expands the possibilities for probing the electronic structure in uranyl complexes beyond the strongly covalent U–O bonds. The combined analyses find that one of the two largest principal EFG tensor components at the halogen sites points along the U–X bond (X = Cl, Br), and the second is parallel to the UO 2 2+ ion; in Cs 2 UO 2 Cl 4 , the components are nearly equal in magnitude, whereas in Cs 2 UO 2 Br 4 , due to short-range bromide–cesium interactions, the equatorial component is dominant for one pair of Br sites and the axial component is larger for the second pair. Additionally, the directions and relative magnitudes of the field gradient principal axes are found to be sensitive to the σ and π electron donation by the ligands and the model of the environment. Chlorine-35 NQR spectra of 235 U-depleted and 235 U-enriched Cs 2 UO 2 Cl 4 exhibited no uranium-isotope-dependent shift, but the resonance of the depleted sample displayed a 58% broader line width.

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

High-throughput characterization of Ag–V–O nanostructured thin-film materials libraries for photoelectrochemical solar water splitting

Ag–V–O thin-film materials libraries, with both composition (Ag 22-77 V 23-78 O x ) and thickness (123–714 nm) gradients were fabricated using combinatorial reactive magnetron co-sputtering aiming on establishing relations between composition, structure, and functional properties. As-deposited libraries were annealed in air at 300 °C for 10 h. High-throughput characterization methods of composition, structure and functional properties were used to identify photoelectrochemically active regions. The phases AgV 6 O 15 , Ag 2 V 4 O 11 , AgVO 3 , and Ag 4 V 2 O 7 were observed throughout the composition gradient. The photoelectrochemical properties of Ag–V–O films are dependent on composition and morphology. An enhanced photocurrent density (~300–554 μA/cm 2 ) was obtained at 30 to 45 at.% Ag along the thickness gradient. Thin films of these compositions show a nanowire morphology, which is an important factor for the enhancement of photoelectrochemical performance. The photoelectrochemically active regions were further investigated by high-throughput synchrotron-X-ray diffraction and transmission electron microscopy (Ag 32 V 68 O x ) which confirmed the presence of Ag 2 V 4 O 11 as the dominating phase along with the minor phases AgV 6 O 15 and AgVO 3 . This enhanced photoactive region shows bandgap values of ~2.30 eV for the direct and ~1.87 eV for the indirect bandgap energies. Finally, the porous nanostructured films improve charge transport and are hence of interest for photoelectrochemical water splitting.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Stabilizing high-Ni cathodes with gradient surface Ti-enrichment

High-Ni cathodes are being intensely pursued worldwide for electric vehicles and other energy-dense applications due to their high capacity and low cost. However, structural instabilities during electrochemical cycling and when subjected to thermal treatment have been the major issues hindering their practical deployment. We here report a rational design of coating-integrated-into-synthesis protocol for fabricating surface Ti-enriched LiNi 0.8 Mn 0.1 Co 0.1 O 2 (NMC811#Ti) material. The coating to intermediates is crucial to obtain high structural ordering, both in the bulk and surface of high-Ni cathodes, and the Ti substitute has a unique tri-valence (Ti 3+ ) in a gradient surface distribution. Further, the simulations of projected density of states in the atomistic understanding further certify significantly enhanced stability of lattice oxygen for the NMC811 through such a Ti 3+ -based structure reinforcement. Consequently, the NMC811#Ti cathode delivers a high capacity up to 200mAhg -1 at 0.1 C, along with superior stabilities during air-storage and thermal treatment (up to 297°C at the fully charged state under differential scanning calorimetric measurements). The corresponding NMC811#Ti||graphite full cell exhibits a desired 83.6% capacity retention after 1000 cycles at 0.5 C in a voltage range of 2.8–4.3V. This work demonstrates a delicate surface reinforcement to stabilize high-Ni cathodes for long-life and safe lithium-ion batteries.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Novel Concepts for High Gradient Acceleration (Final Technical Report)

We have conducted an intensive, pioneering program to demonstrate novel concepts for achieving high gradient acceleration at frequencies from the conventional microwave bands up to the millimeter wave /THz bands. High gradient accelerators hold the promise of smaller and less costly accelerators for applications ranging from the largest scale accelerators used for discovery science down to the smallest accelerators used for industrial, homeland security and medical applications. The research consisted of two major research thrusts: 1.) Structure-based wakefield accelerator (SWFA) research in collaboration with the Argonne Wakefield Accelerator research group and 2.) Millimeter Wave / THz high gradient acceleration in collaboration with SLAC. The specific goals of the research program were: Design novel metallic metamaterial structures that increase the beam-wave coupling for the accelerator mode and reduce the effect of high order modes; test novel metamaterial structures to achieve higher output power, > 1 GW at X-Band (11.7 GHz), in test at the Argonne Wakefield Accelerator (AWA); determine the experimental breakdown threshold for nanosecond-scale pulses at X-Band in testing at the AWA; test a 110 GHz accelerator structure with a field emission electron gun, built at SLAC, using pulses from a 1 MW, 110 GHz gyrotron; design, build and test a 110 GHz quasi-optical, resonant-ring pulse compressor to compress microsecond pulses from the 1 MW gyrotron into > 20 MW, 5 ns output pulses for accelerator structure testing. The proposed research program built on our successes in our research program including: Generation of 510 MW, 2.1 ns (FWHM) pulses at 11.7 GHz from a metallic metamaterial structure in test at the Argonne Wakefield Accelerator using a train of eight 65 MeV electron bunches spaced at 1.3 GHz with a total charge of 280 nC. The metamaterial structure consisted of 100 copper unit cells each consisting of a “wagon-wheel” plate and a spacer plate with a total structure length of 0.2 m. The 510 MW pulse generated an on-axis wakefield of 130 MV/m that could be used to accelerate a trailing witness bunch. Demonstration of coupling of an unprecedented rf power level of 575kW into a 110 GHz accelerator structure using a quasi-optical setup. The standing structure consisted of a central copper cavity located between two matching cavities fed by a TM01 mode. The 6 ns input pulses were sliced from 3 microsecond pulses from the gyrotron using a laser-driven silicon switch. We obtained an unprecedented high gradient up to 230MV/m corresponding to a peak surface electric field of more than 520 MV/m.

43 PARTICLE ACCELERATORS↗
DOE OSTI2026

Analytic Nuclear Gradients Including Oriented External Electric Fields in a Molecule-Fixed Frame

Electric-field-assisted chemistry has attracted much attention in recent years, particularly in the context of oriented external electric fields for controlling molecular structure and reactivity. Such fields have been explored in a wide range of applications, including switching materials, nanoparticles, controllable catalysts, medicines, and clinical therapies. However, the determination of fixed fields in the laboratory frame becomes ineffective for flexible molecules, as conformational changes can significantly alter the relative orientation between the applied field and molecular structure. In this work, we propose two molecular reference frames─the principal axis frame and the local reference frame─to define oriented electric fields within the molecular framework. These coordinate systems powerfully eliminate ambiguities in the relative orientation between the applied field and the molecule. Analytic nuclear gradients in the presence of external electric fields are derived and implemented, with an initial application to field-dependent geometry optimizations of cis - and trans -formanilide. Analysis of the resulting field-induced equilibrium structures reveals distinct structural responses, validating the accuracy and robustness of the proposed formalism. The analytic gradient framework enables systematic investigations of molecular properties and reactivity under arbitrarily oriented electric fields, opening new opportunities for computational modeling and rational design in electric-field-controlled chemistry.

electric fields↗
DOE OSTI2022

A deep learning-enhanced framework for multiphysics joint inversion

Joint inversion has drawn considerable attention due to the availability of multiple geophysical data sets, ever-increasing computational resources, the development of advanced algorithms, and its ability to reduce inversion uncertainty. A key issue of joint inversion is to develop effective strategies to link different geophysical data in a unified mathematical framework, in which the information obtained from different models can complement each other. We have developed a deep learning-enhanced joint inversion framework to simultaneously reconstruct different physical models by fusing different types of geophysical data. Traditionally, structure similarity constraints are pursued by joint inversion algorithms using manually crafted formulations (e.g., cross gradient). The constraint is constructed by a deep neural network (DNN) during the learning process. The framework is designed to combine the DNN and a traditional independent inversion workflow and improve the joint inversion result iteratively. The network can be easily extended to incorporate multiphysics without structural changes. Numerical experiments on the joint inversion of 2D DC resistivity data and seismic traveltime are used to validate our method. In addition, this learning-based framework demonstrates excellent generalization abilities when tested on data sets using different geologic structures. It also can handle different sensing configurations and nonconforming discretization.

Geochemistry & Geophysics↗
DOE OSTI2021

Nucleon-pair coupling scheme in Elliott's SU(3) model

Elliott's SU(3) model is at the basis of the shell-model description of rotational motion in atomic nuclei. Here we demonstrate that SU(3) symmetry can be realized in a truncated shell-model space if constructed in terms of a sufficient number of collective S, D, G,...pairs (i.e., with angular momentum zero, two, four,...) and if the structure of the pairs is optimally determined either by a conjugate-gradient minimization method or from a Hartree-Fock intrinsic state. We illustrate the procedure for six protons and six neutrons in the pf (sdg) shell and exactly reproduce the level energies and electric quadrupole properties of the ground-state rotational band with SDG (SDGI) pairs. The SD-pair approximation without significant renormalization, on the other hand, cannot describe the full SU(3) collectivity. A mapping from Elliott's fermionic SU(3) model to systems with s, d, g,... bosons provides insight into the existence of a decoupled collective subspace in terms of S, D, G,... pairs.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Full Concentration Gradient-Tailored Li-Rich Layered Oxides for High-Energy Lithium-Ion Batteries

Lithium-rich layered oxides (LLOs) are prospective cathode materials for next-generation lithium-ion batteries (LIBs), but severe voltage decay and energy attenuation with cycling still hinder their practical applications. In this study, a series of full concentration gradient-tailored agglomerated-sphere LLOs are designed with linearly decreasing Mn and linearly increasing Ni and Co from the particle center to the surface. The gradient-tailored LLOs exhibit noticeably reduced voltage decay, enhanced rate performance, improved cycle stability, and thermal stability. Without any material modifications or electrolyte optimizations, the gradient-tailored LLO with medium-slope shows the best electrochemical performance, with a very low average voltage decay of 0.8 mV per cycle as well as a capacity retention of 88.4% within 200 cycles at 200 mA g –1 . These excellent findings are due to spinel structure suppression, electrochemical stress optimization, and Jahn-Teller effect inhibition. Further investigation shows that the gradient-tailored LLO reduces the thermal release percentage by as much as about 41% when the battery is charged to 4.4 V. This study provides an effective method to suppress the voltage decay of LLOs for further practical utilization in LIBs and also puts forward a bulk-structure design strategy to prepare better electrode materials for different rechargeable batteries.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Optical response and band structure of LiCoO 2 including electron-hole interaction effects

The optical response functions and band structures of LiCoO 2 are studied at different levels of approximation, from density functional theory (DFT) in the generalized gradient approximation (GGA) to quasiparticle self-consistent QSGW (with G for Green's function and W for screened Coulomb interaction) without and with ladder diagrams ($QSG\hat{W}$) and the Bethe Salpeter Equation (BSE) approach. The QSGW method is found to strongly overestimate the band gap and electron-hole or excitonic effects are found to be important. They lower the quasiparticle gap by only about 11~\% but the lowest energy peaks in absorption are found to be excitonic in nature. The contributions from different band to band transitions and the relation of excitons to band-to-band transitions are analyzed. Finally, the excitons are found to be strongly localized. Finally, a comparison to experimental data is presented.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Microbial and Environmental Processes Shape the Link between Organic Matter Functional Traits and Composition

Dissolved organic matter (DOM) is a large and complex mixture of molecules that fuels biogeochemical reaction in virtually all ecosystems on Earth. However, the relative importance of deterministic and stochastic processes in structuring DOM composition remains poorly characterized. Here we develop a framework for partitioning molecular composition based on key molecular traits, including lability vs. recalcitrance and activity vs. inactivity. Within this framework, we examine the ecological processes governing the assembly of DOM fractions by deploying aquatic microcosms on mountainsides that span gradients of temperature and nutrient loading in subtropical and subarctic ecosystems. Across study regions, deterministic and stochastic processes primarily structure active and inactive fractions, respectively. However, recalcitrant molecules are more deterministically assembled than labile molecules in the inactive fraction. Deterministic processes leading to variable selection generally exhibit more variation across the energy supply gradient for inactive fractions, and their importance increases with energy supply for recalcitrant molecules in both active and inactive fractions. Together, our results indicate that active and inactive fractions of DOM assemblages are structured by contrasting ecological processes, and their recalcitrant components are sensitive to global change. In conclusion, our framework opens new avenues to understand the assembly and turnover of DOM in a changing world, which can be used to predict carbon cycling at local to global scales.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Structural and dynamic properties of solvated hydroxide and hydronium ions in water from ab initio modeling

Predicting the asymmetric structure and dynamics of solvated hydroxide and hydronium in water from ab initio molecular dynamics (AIMD) has been a challenging task. The difficulty mainly comes from a lack of accurate and efficient exchange–correlation functional in elucidating the amphiphilic nature and the ubiquitous proton transfer behaviors of the two ions. By adopting the strongly constrained and appropriately normed (SCAN) meta-generalized gradient approximation functional in AIMD simulations, we systematically examine the amphiphilic properties, the solvation structures, the electronic structures, and the dynamic properties of the two water ions. In particular, we compare these results to those predicted by the PBE0-TS functional, which is an accurate yet computationally more expensive exchange–correlation functional. We demonstrate that the general-purpose SCAN functional provides a reliable choice for describing the two water ions. Specifically, in the SCAN picture of water ions, the appearance of the fourth and fifth hydrogen bonds near hydroxide stabilizes the pot-like shape solvation structure and suppresses the structural diffusion, while the hydronium stably donates three hydrogen bonds to its neighbors. We apply a detailed analysis of the proton transfer mechanism of the two ions and find the two ions exhibit substantially different proton transfer patterns. In conclusion, our AIMD simulations indicate that hydroxide diffuses more slowly than hydronium in water, which is consistent with the experimental results.

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

Evaluation of thermal energy storage for central receiver concentrating solar power plants under charging cycles through structural assessment

Electricity production by concentrated solar power (CSP) systems has stood out among energy transition alternatives due to their large generation capacity (between 30 MW and 400 MW) and high capacity factor (80%) when incorporating thermal storage systems (TES). Although TES tanks in operation worldwide are made of austenitic steels with remarkable mechanical and anticorrosive properties, operational failures have been reported due to low-cycle fatigue, creep, and the abrupt release of residual stresses generated during tank manufacturing. A critical operation for salt tanks is the initial daily charging operation, given the low fluid level and the thermal gradients between the inventory and the salt entering the tank. This study presents a structural simulation of a 39.6 m-diameter molten salt tank for central receiver CSP plants, accounting for temperature and pressure variations during charging and different sparger ring inlet configurations. The initial pre-stressed condition of the floor, resulting from the manufacture of welded plates, is used as a boundary condition. The fatigue and creep life analyses were performed in accordance with the ASME BPVC and API standards. It was found that the tank's floor is in a less favourable mechanical condition than the wall, due to its initial manufacturing condition and the maximum temperature differences during charging. A critical zone was identified in an area affected by residual stress on the floor and under the sparger ring. The baseline sparger ring configuration, with 52 2-inch-diameter orifices and a flow direction of 90 degrees , results in an average floor stress of approximately 25 MPa and creep damage after 1.2 years of tank operation.

14 SOLAR ENERGY↗
DOE OSTI2024

Coupled Cluster Theory for Nonadiabatic Dynamics: Nuclear Gradients and Nonadiabatic Couplings in Similarity Constrained Coupled Cluster Theory

Coupled cluster theory is one of the most accurate electronic structure methods for predicting ground and excited state chemistry. However, the presence of numerical artifacts at electronic degeneracies, such as complex energies, has made it difficult to apply the method in nonadiabatic dynamics simulations. While it has already been shown that such numerical artifacts can be fully removed by using similarity constrained coupled cluster (SCC) theory [J. Phys. Chem. Lett. 2017, 8(19), 4801–4807], simulating dynamics requires efficient implementations of gradients and nonadiabatic couplings. Here, we present an implementation of nuclear gradients and nonadiabatic derivative couplings at the similarity constrained coupled cluster singles and doubles (SCCSD) level of theory, thereby making possible nonadiabatic dynamics simulations using a coupled cluster theory that provides a correct description of conical intersections between excited states. We present a few numerical examples that show good agreement with literature values and discuss some limitations of the method.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2025

Radiation-induced bowing of SiC/SiC composites under neutron flux gradients—integral experimental data for model validation

Here, the radiation-induced swelling of SiC and its composites, including strong dependencies on temperature and dose, can drive significant lateral bowing in the presence of temperature and/or dose gradients. In recent years, simulations have been performed to assess the extent of bowing in SiC composite light-water reactor (LWR) fuel cladding and boiling water reactor (BWR) channel boxes. However, to date, no integral experimental data exist to validate these models. This work provides the first experimental bowing evaluation of three ∼380 mm long SiC composite specimens irradiated under varying neutron dose gradients (∼50°C–60°C, 0.03–0.06 dpa): two tubes (∼9.8 mm diameter) and a miniature BWR channel box (∼30 mm square). The measured radiation-induced length swelling (∼0.3%–0.7% linear) was consistently 10%–21% higher than values obtained from 3D finite element structural analyses with inputs from 3D radiation transport calculations. This discrepancy could be at least partially explained by differences in dose rate (∼10 -8 dpa/s) compared to the literature data (∼10-6 dpa/s) used to establish the dose-to-swelling correlations in the model. Nevertheless, the modeled bowing magnitudes (<2 mm) obtained from finite element analyses and simple analytical equations were within the bounds of the experimental measurements for all specimens. With improved confidence in the ability to predict the structural response and measure the macroscopic deformations, future experiments will target transient bowing under neutron flux gradients at representative LWR temperatures and assess whether grid spacers can mitigate the tens of millimeters of bowing that would otherwise be expected in ∼4 m long LWR components.

bowing↗
DOE OSTI2021

Analytical derivatives of the individual state energies in ensemble density functional theory. II. Implementation on graphical processing units (GPUs)

Conical intersections control excited state reactivity, and thus, elucidating and predicting their geometric and energetic characteristics are crucial for understanding photochemistry. Locating these intersections requires accurate and efficient electronic structure methods. Unfortunately, the most accurate methods (e.g., multireference perturbation theories such as XMS-CASPT2) are computationally challenging for large molecules. The state-interaction state-averaged restricted ensemble referenced Kohn–Sham (SI-SA-REKS) method is a computationally efficient alternative. The application of SI-SA-REKS to photochemistry was previously hampered by a lack of analytical nuclear gradients and nonadiabatic coupling matrix elements. We have recently derived analytical energy derivatives for the SI-SA-REKS method and implemented the method effectively on graphical processing units. We demonstrate that our implementation gives the correct conical intersection topography and energetics for several examples. Furthermore, our implementation of SI-SA-REKS is computationally efficient, with observed sub-quadratic scaling as a function of molecular size. This demonstrates the promise of SI-SA-REKS for excited state dynamics of large molecular systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗