Search NASA⌕ Search

SEARCH · Search NASA

Results for “Derived”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

The strip method for shape derivatives

Abstract A major challenge in shape optimization is the coupling of finite element method (FEM) codes in a way that facilitates efficient computation of shape derivatives. This is particularly difficult with multiphysics problems involving legacy codes, where the costs of implementing and maintaining shape derivative capabilities are prohibitive. The volume and boundary methods are two approaches to computing shape derivatives. Each has a major drawback: the boundary method is less accurate, while the volume method is more invasive to the FEM code. We introduce the strip method , which computes shape derivatives on a strip adjacent to the boundary. The strip method makes code coupling simple. Like the boundary method, it queries the state and adjoint solutions at quadrature nodes, but requires no knowledge of the FEM code implementations. At the same time, it exhibits the higher accuracy of the volume method. As an added benefit, its computational complexity is comparable to that of the boundary method, that is, it is faster than the volume method. We illustrate the benefits of the strip method with numerical examples.

Hardesty, Sean↗
DOE OSTI2022

Four-derivative corrections to minimal gauged supergravity in five dimensions

We study four-derivative corrections to pure N = 2, D = 5 gauged supergravity. In particular, we find that, up to field redefinitions, there is a single four-derivative superinvariant that one can add to the action, up to factors of the two-derivative action. Consequently, this selects out a unique set of coefficients for the four-derivative corrections. We confirm these coefficients (in the ungauged limit) on the BMPV solution.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Derivation of K-matrix reaction theory in a discrete basis formalism

The usual derivations of the S and K matrices for two-particle reactions proceed through the Lippmann–Schwinger equation with formal definitions of the incoming and outgoing scattering states. Here we present a simpler alternative derivation that is carried out completely in the Hamiltonian representation, using a discrete basis of configurations for the scattering channels as well as the quasi-bound configurations of the combined fragments. We use matrix algebra to derive an explicit expression for the K matrix in terms of the Hamiltonian of the internal states of the compound system and the coupling between the channels and the internal states. The formula for the K matrix includes explicitly a real dispersive shift matrix to the internal Hamiltonian that is easily computed in the formalism. That expression is applied to derive the usual form of the S matrix as a sum over poles in the complex energy plane. Some extensions and limitations of the discrete-basis Hamiltonian formalism are discussed in the concluding remarks and in the Appendix.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Subtle penetrant size effects on separation of carbon molecular sieve membranes derived from $\mathrm{6FDA:BPDA}$-$\mathrm{DAM}$ polyimide

Economically scalable carbon molecular sieve (CMS) hollow fiber membranes rely upon a tunable bimodal pore morphology to separate gas pairs using angstrom-level size discrimination. Freshly formed CMS membranes experience self-retarding physical aging, reflected by permeance losses and selectivity gains that stabilize over short times compared to the lifetime of the membrane. Self-retarding morphology rearrangements are of special interest here due to tightening of the largest size range of ultramicropores, thereby causing aging for polyimide-derived CMS. We report effects of such aging for two A/B penetrant pairs, C 3 H 6 /C 3 H 8 and CO 2 /CH 4 . These two pairs not only have different average sizes, $\overline{d}_{A/B}$, but also different size differences,Δd A/B , between the members in each pair. We focus primarily on the C 3 H 6 /C 3 H 8 pair, which is the most difficult of the two pairs to separate, and use some CO 2 /CH 4 results as a comparison case. We study CMS derived from 6FDA:BPDA-DAM (1:1) polyimide by pyrolysis at 550 °C, 600 °C and 650 °C. Here we suggest conditions and physical causes that expedite, retard or even suppress the aging process within CMS membranes. Finally, we suggest how analysis of CMS derived from other precursors and other penetrant pairs can be generalized by understanding the results for the 6FDA:BPDA-DAM (1:1) polyimide derived CMS.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Derivative-based SINDy (DSINDy): Addressing the challenge of discovering governing equations from noisy data

Recent advances in the field of data-driven dynamics allow for the discovery of ODE systems using state measurements. One approach, known as Sparse Identification of Nonlinear Dynamics (SINDy), assumes the dynamics are sparse within a predetermined basis in the states and finds the expansion coefficients through linear regression with sparsity constraints. This approach requires an accurate estimation of the state time derivatives, which is not necessarily possible in the high-noise regime without additional constraints. We present an approach called Derivative-based SINDy (DSINDy) that combines two novel methods to improve ODE recovery at high-noise levels. First, we denoise the state variables by applying a projection operator that leverages the assumed basis for the system dynamics. Second, we use a second order cone program (SOCP) to find the derivative and governing equations simultaneously. We derive theoretical results for the projection-based denoising step, which allow us to estimate the values of hyperparameters used in the SOCP formulation. This underlying theory helps limit the number of required user-specified parameters. Finally, we present results demonstrating that our approach leads to improved system recovery for the Van der Pol oscillator, the Duffing oscillator, the Rössler attractor, and the Lorenz 96 model.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Life-cycle energy use and greenhouse gas emissions of palm fatty acid distillate derived renewable diesel

This study aims to quantify life-cycle fossil energy use and greenhouse gas (GHG) emissions for palm fatty acid distillate (PFAD) derived renewable diesel (RD) taking into consideration different feedstock classifications that are applicable to PFAD (residue, byproduct, or coproduct), and incorporating updated data for key processes. Under the three classifications, the PFAD to RD pathway was modeled using the Greenhouse gases, Regulated Emissions, and Energy Use in Technologies (GREET®) model. PFAD-derived RD could reduce fossil energy consumption by 77%–88%, relative to petroleum diesel. GHG emissions are very sensitive to PFAD classification and coproduct handling methods. Considering the production of palm oil and PFAD and economic value, we maintain that PFAD should be treated as a byproduct in palm oil refineries. With this treatment, PFAD-derived RD could achieve 84% GHG emissions reductions, compared to the emissions of petroleum diesel. We also employed a substitution method to address the substitution of PFAD by other materials in the marketplace. Compared to coproduct allocation results, we found substituting PFAD by tallow, soy oil, barley, and canola oil results in lower GHG emissions. Due to high induced land-use change emissions associated with palm farming, if PFAD is treated as a coproduct with refined palm oil, PFAD-derived RD may not deliver GHG reductions. A sensitivity analysis identified key parameters such as palm fruit yield, oil extraction efficiency in oil mills, and energy use intensity for RD production affects LCA results significantly; future efforts to improve these parameters could result in further GHG reductions.

09 BIOMASS FUELS↗
DOE OSTI2021

Sorafenib derivatives-functionalized gold nanoparticles confer protection against tumor angiogenesis and proliferation via suppression of EGFR and VEGFR-2

Sorafenib is a multi-kinase inhibitor that has been highlighted as a tumor suppressor due to its anti-proliferative and anti-angiogenic properties, whereas the clinical application of Sorafenib is restricted by the side effects it may cause. The past decade has witnessed the development of a series of sorafenib derivatives to improve the clinical performance of sorafenib. Gold nanoparticles (AuNPs) have been widely utilized in drug delivery systems due to their unique properties, including biocompatible nature, simple preparation, and easy surface modification. Herein, this study is aimed to investigate the anti-tumor effects of new sorafenib derivatives-capped gold nanoparticles (AuNPs-New Sor) in tumor formation and metastasis as well as the underlying mechanisms. Initially, new sorafenib derivatives were constructed and combined with AuNPs to form AuNPs-New Sor, and the properties of synthesized AuNPs-New Sor were identified in a mouse model of tumorigenesis. The effect of AuNPs-New Sor on tumor vascular normalization was investigated by assessing vascular permeability and perfusion rate. Next, we evaluated the effect of AuNPs-New Sor on migration and viability of tumor cells and human umbilical vein endothelial cells (HUVECs) as well as on HUVEC angiogenesis in vitro. A melanoma mouse model was further established for in vivo substantiation of the anti-tumor effect of AuNPs-New Sor. According to the results, AuNPs could deliver new sorafenib derivatives into tumor tissues and downregulate the expression of epidermal growth factor receptor (EGFR) and vascular endothelial growth factor receptor-2 (VEGFR-2), thereby suppressing tumor migration, EMT, and angiogenesis in vitro. In addition, AuNPs-New Sor displayed competitive anti-tumor activities in vivo. Taken together, AuNPs-New Sor may attenuate tumor development and angiogenesis through downregulation of EGFR and VEGFR-2.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2021

MCL1 inhibition enhances the efficacy of docetaxel against airway-derived squamous cell carcinoma cells

Highlights: • Authentic squamous cell carcinoma (SqCC) cell lines of the lung are sparse. • HARA lung SqCC cells depend highly on MCL1 for survival in vitro and in vivo. • MCL1-silenced HARA cells decline in number upon docetaxel (DTX) exposure. • High MCL1 mRNA expression is linked with shorter survival in lung SqCC patients. • MCL1 depletion + DTX is highly effective against Detroit 562 pharyngeal SqCC cells. MCL1 is an anti-apoptotic BCL2 family member that is often overexpressed in various malignant tumors. However, few reports have described the role of MCL1 in squamous cell carcinoma (SqCC) derived from airways including the lung. In this study, we examined whether MCL1 could be a novel druggable target for airway-derived SqCC, for which effective molecular targeted drugs are unavailable. We searched the Kaplan–Meier Plotter database and found that high MCL1 mRNA expression was significantly associated with shorter survival in patients with lower airway (lung) or upper airway (head and neck) derived SqCC. We also explored the Expression Atlas database and learned that authentic lung SqCC cell lines expressing both TP63 and KRT5 mRNA were extremely sparse among the publicly available “lung SqCC cell lines”, with an exception being HARA cells. HARA cells were highly dependent on MCL1 for survival, and MCL1-depleted cells were not able to grow, and even declined in number, upon docetaxel (DTX) exposure in vitro and in vivo. Similar in vitro experimental findings, including those in a 3D culture model, were also obtained using Detroit 562 pharyngeal SqCC cells. These findings suggested that combined treatment with MCL1 silencing plus DTX appears highly effective against airway-derived SqCC.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2021

Metabolism of Syringyl Lignin-Derived Compounds in Pseudomonas putida Enables Convergent Production of 2-Pyrone-4,6-Dicarboxylic Acid

Valorization of lignin, an abundant component of plant cell walls, is critical to enabling the lignocellulosic bioeconomy. Biological funneling using microbial biocatalysts has emerged as an attractive approach to convert complex mixtures of lignin depolymerization products to value-added compounds. Ideally, biocatalysts would convert aromatic compounds derived from the three canonical types of lignin: syringyl (S), guaiacyl (G), and p-hydroxyphenyl (H). Pseudomonas putida KT2440 (hereafter KT2440) has been developed as a biocatalyst owing in part to its native catabolic capabilities but is not known to catabolize S-type lignin-derived compounds. Here, we demonstrate that syringate, a common S-type lignin-derived compound, is utilized by KT2440 only in the presence of another energy source or when vanAB was overexpressed, as syringate was found to be O-demethylated to gallate by VanAB, a two-component monooxygenase, and further catabolized via extradiol cleavage. Unexpectedly, the specificity (kcat/KM) of VanAB for syringate was within 25% that for vanillate and O-demethylation of both substrates was well-coupled to O2 consumption. However, the native KT2440 gallate-cleaving dioxygenase, GalA, was potently inactivated by 3-O-methylgallate. To engineer a biocatalyst to simultaneously convert S-, G-, and H-type monomers, we therefore employed VanAB from Pseudomonas sp. HR199, which has lower activity for 3MGA, and LigAB, an extradiol dioxygenase able to cleave protocatechuate and 3-O-methylgallate. This strain converted 93% of a mixture of lignin monomers to 2-pyrone-4,6-dicarboxylate, a promising bio-based chemical. Overall, this study elucidates a native pathway in KT2440 for catabolizing S-type lignin-derived compounds and demonstrates the potential of this robust chassis for lignin valorization.

2-pyrone-4,6-dicarboxylic acid↗
DOE OSTI2021

Controlling selectivity of modular microbial biosynthesis of butyryl-CoA-derived designer esters

Short-chain esters have broad utility as flavors, fragrances, solvents, and biofuels. Controlling selectivity of ester microbial biosynthesis has been an outstanding metabolic engineering problem. In this study, we enabled the de novo fermentative microbial biosynthesis of butyryl-CoA-derived designer esters (e.g., butyl acetate, ethyl butyrate, butyl butyrate) in Escherichia coli with controllable selectivity. Using the modular design principles, we generated the butyryl-CoA-derived ester pathways as exchangeable production modules compatible with an engineered chassis cell for anaerobic production of designer esters. We designed these modules derived from an acyl-CoA submodule (e.g., acetyl-CoA, butyryl-CoA), an alcohol submodule (e.g., ethanol, butanol), a cofactor regeneration submodule (e.g., NADH), and an alcohol acetyltransferase (AAT) submodule (e.g., ATF1, SAAT) for rapid module construction and optimization by manipulating replication (e.g., plasmid copy number), transcription (e.g., promoters), translation (e.g., codon optimization), pathway enzymes, and pathway induction conditions. To further enhance production of designer esters with high selectivity, we systematically screened various strategies of protein solubilization using protein fusion tags and chaperones to improve the soluble expression of multiple pathway enzymes. Finally, our engineered ester-producing strains could achieve 19-fold increase in butyl acetate production (0.64 g/L, 96% selectivity), 6-fold increase in ethyl butyrate production (0.41 g/L, 86% selectivity), and 13-fold increase in butyl butyrate production (0.45 g/L, 54% selectivity) as compared to the initial strains. Altogether, this study presented a generalizable framework to engineer modular microbial platforms for anaerobic production of butyryl-CoA-derived designer esters from renewable feedstocks.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

Synthesis and Properties of Bulk Mg 3 WN 4 in a Wurtzite-Derived Structure

Experimental synthesis of theoretically predicted materials with controlled elemental coordination environments can lead to the realization of useful properties, such as facile ion transport or ferroelectric switching. Among such materials are ternary nitrides in the Mg−W−N composition space, where several stable and metastable compounds have been predicted and synthesized recently in bulk and film forms. Here, we report for the first time on the bulk synthesis of Mg 3 WN 4 in a wurtzite-derived crystal structure via a solid-state metathesis reaction. In situ synchrotron powder X-ray diffraction shows how the ion exchange proceeds from Li 6 WN 4 + MgCl 2 precursors to Mg 3 WN 4 + LiCl products, with the reaction starting slowly near 380 °C and completing by 600 °C, including the presence of a competing disordered rocksalt-derived phase (Mg,W)N above 440 °C. The follow-up ex situ powder synthesis at 400 °C for 0.5 h with 10% excess MgCl 2 reveals the cation-ordered nature of the wurtzite-derived Mg 3 WN 4 structure with polar symmetry confirmed by second-harmonic generation measurements. Optical absorption spectra, chemical composition analysis, and electron microscopy imaging suggest that bulk Mg 3 WN 4 obtained via metathesis reaction is prone to defect formation. Overall, this study shows that selective ex situ synthesis of the phase-pure ternary nitrides, informed by in situ measurements, is possible by carefully controlling the thermal budget of the reaction and paves the way toward property characterization of the Mg 3 WN 4 wurtzite-derived material.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Recent Advances toward Efficient Calculation of Higher Nuclear Derivatives in Quantum Chemistry

In this article, we provide an overview of state-of-the-art techniques that are being developed for efficient calculation of second and higher nuclear derivatives of quantum mechanical (QM) energy. Calculations of nuclear Hessians and anharmonic terms incur high costs and memory and scale poorly with system size. Three emerging classes of methods—machine learning (ML), automatic differentiation (AD), and matrix completion (MC)—have demonstrated promise in overcoming these challenges. We illustrate studies that employ unsupervised ML methods to reduce the need for multiple Hessian calculations in dynamics simulations and those that utilize supervised ML to construct approximate potential energy surfaces and estimate Hessians and anharmonic terms at reduced cost. By extension, if electronic structure operations could be written in a manner similar to functions underlying ML methods, rapid differentiation or AD routines can be employed to inexpensively calculate higher arbitrary-order derivatives. While ML approaches are typically black-box, we describe methods such as compressed sensing (CS) and MC, which explicitly leverage problem-specific mathematical properties of higher derivatives such as sparsity and low-rank, to complete higher derivative information using only a small, incomplete sample. The three classes of methods facilitate reliable predictions of observables ranging from infrared spectra to thermal conductivity and constitute a promising way forward in accurately capturing otherwise intractable higher-order responses of QM energy to nuclear perturbations.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2023

R –8 Dispersion Interaction: Derivation and Application to the Effective Fragment Potential Method

The anisotropic and isotropic R –8 dispersion contributions (disp8) are derived and implemented within the framework of the effective fragment potential (EFP) method formulated with imaginary frequency-dependent Cartesian polarizability tensors distributed at the centroids of the localized molecular orbitals (LMOs). Two forms of damping functions, intermolecular overlap-based and Tang–Toennies, are extended for disp8. To obtain LMO polarizability tensors centered at LMO centroids, an origin-shifting transformation is derived and implemented for the dipole–octopole polarizability tensor and the quadrupole–quadrupole polarizability tensor. The analytic gradient is derived and implemented for the isotropic disp8 contribution. Relative to the previously implemented empirical EFP disp8 energy, the isotropic disp8 component of the interaction energy improves the overall agreement of the EFP dispersion energies with the symmetry-adapted perturbation theory (SAPT) benchmarks, reducing the mean absolute errors (MAEs) and mean absolute percentage errors for most of the databases examined in this work. Further, while the anisotropic disp8 can further enhance the accuracy of the EFP dispersion energy and yield smaller MAEs, significantly overbound dispersion energies are predicted by the anisotropic disp8 when the maximum element in the intermolecular overlap matrix is greater than 0.1, possibly due to the breakdown of the approximations made in the EFP dispersion derivation at a short range. For potential energy scan databases, the newly developed EFP dispersion model with isotropic disp8 yields the overall correct curvature and good agreement with SAPT benchmarks around equilibrium and longer but overestimates the dispersion interactions at a short range. While the overlap-based dispersion-damping functions produce better MAEs than Tang–Toennies damping functions, further improvement is needed to better screen the large attractive dispersion energies at a short range (overlap >0.1).

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

CO 2 Chemisorption Behavior in Conjugated Carbanion-Derived Ionic Liquids via Carboxylic Acid Formation

Superbase-derived task-specific ionic liquids (STSILs) represent one of the most attractive and extensively studied systems in carbon capture via chemisorption, in which the obtained CO 2 uptake capacity has a strong relationship with the basicity of the anions. High energy input in desorption and side reactions caused by the strong basicity of the anions are still unsolved issues. The development of other customized STSILs leveraging an alternative driving force to achieve efficient CO 2 chemisorption/desorption is highly desirable yet challenging. Here, in this work, carbanion-derived STSILs were developed for efficient CO 2 chemisorption via a carboxylic acid formation pathway. The STSIL with the deprotonated malononitrile molecule ([MN]) as the anion exhibited much higher CO 2 uptake capacity than the one derived from 2-methylmalononitrile ([MMN]). Notably, this trend was opposite to their basicity ([MN] < [MMN]). Detailed characterization of the products, supported by density functional theory simulations of spectra and calculations of the reaction energetics, demonstrated that carboxylic acid was formed upon reacting with CO 2 via proton transfer in [MN]-derived STSILs but not in the case of [MMN] due to lack of an α-H. The preference of the carboxylic acid product over carboxylate formation was driven by the extended conjugation among the central sp 2 carbon, the as-formed carboxylic acid, and the two nitrile groups. The achievements made in this work provide an alternative design principle of STSILs by leveraging the extended conjugation in the CO 2 -integrated product.

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

Biomass-Derived Polymeric Binders in Silicon Anodes for Battery Energy Storage Applications

The demand for portable electronic devices has increased rapidly during the past decade, and has driven a concordant growth in battery production. Since their development as a commercial energy storage solution in the 1990s, lithium-ion batteries (LIBs) have attracted significant attention in both science and industry due to their long cycle life, high energy density, low self-discharge rate, and high working voltage. Production of LIBs requires large amounts of a polymeric binder – commonly polyvinylidene difluoride (PVDF) – for processing and performance purposes. However, since this material is petrochemically-derived, it is far from “green” or sustainable. On the other hand, polymers and their building blocks are found widely throughout nature and can be renewably sourced from biomass at low cost; therefore, replacing PVDF with biomass-derived binders is a promising approach to reduce the environmental footprint of LIBs. Additionally, polymer binders play a critical role in next-generation battery performance. For instance, silicon (Si) is a promising high-capacity anode material for LIBs because of its high theoretical capacity (4200 mA h g –1 ), low working potential, and high abundance in Earth's crust. However, its huge volume change during charge/discharge tends to result in a shortened cycle life, since conventional binders interact only weakly with silicon's native surface and cannot maintain long-term integrity of the electrode. Naturally derived polymers have found better success in this role due to their high structural advantages. In this review, we summarize recent developments in silicon anode binders derived from various biomass sources, with a focus on polymer properties and their effect on battery performance. Further, we propose various perspectives based on our own assessment of these works, and provide a brief commentary on the future outlook of the field.

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

Controlling the toxicity of biomass-derived difunctional molecules as potential pharmaceutical ingredients for specific activity toward microorganisms and mammalian cells

A biomass-derived difuran compound, denoted as HAH (HMF–acetone–HMF), synthesized by aldol-condensation of 5-hydroxyfurfural (HMF) and acetone, can be partially hydrogenated to provide an electron-rich difuran compound (PHAH) for Diels–Alder reactions with maleimide derivatives. The nitrogen (N) site in the maleimide can be substituted by imidation with amine-containing compounds to control the hydrophobicity of the maleimide moiety in adducts of furans and maleimide by Diels–Alder reaction, denoted as norcantharimides (Diels–Alder adducts). The structural effects on the toxicity of various biomass-derived small molecules synthesized in this manner to regulate biological processes, defined as low molecular weight (≤1000 g mol -1 ) organic compounds, were investigated against diverse microbial and mammalian cell types. The biological toxicity increased when hydrophobic N-substitutions and C[double bond, length as m-dash]C bonds were introduced into the molecular structure. Among the synthesized norcantharamide derivatives, some compounds demonstrated pH-dependent toxicities against specific cell types. Reaction kinetics analyses of the norcantharimides in biological conditions suggest that this pH-dependent toxicity of norcantharimides could arise from retro Diels–Alder reactions in the presence of a Brønsted acid that catalyzes the release of an N-substituted maleimide, which has higher toxicity against fungal cells than the toxicity of the Diels–Alder adduct. These synthetic approaches can be used to design biologically-active small molecules that exhibit selective toxicity against various cell types (e.g., fungal, cancer cells) and provide a sustainable platform for production of prodrugs that could actively or passively affect the viability of infectious cells.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2020

On the proper derivation of the Floquet-based quantum classical Liouville equation and surface hopping describing a molecule or material subject to an external field

We investigate different approaches to derive the proper Floquet-based quantum–classical Liouville equation (F-QCLE) for laser-driven electron-nuclear dynamics. The first approach projects the operator form of the standard QCLE onto the diabatic Floquet basis and then transforms to the adiabatic representation. The second approach directly projects the QCLE onto the Floquet adiabatic basis. Both approaches yield a form that is similar to the usual QCLE with two modifications: (1) The electronic degrees of freedom are expanded to infinite dimension and (2) the nuclear motion follows Floquet quasi-energy surfaces. However, the second approach includes an additional cross derivative force due to the dual dependence on time and nuclear motion of the Floquet adiabatic states. Our analysis and numerical tests indicate that this cross derivative force is a fictitious artifact, suggesting that one cannot safely exchange the order of Floquet state projection with adiabatic transformation. Our results are in accord with similar findings by Izmaylov et al., [J. Chem. Phys. 140, 084104 (2014)] who found that transforming to the adiabatic representation must always be the last operation applied, although now we have extended this result to a time-dependent Hamiltonian. Finally, this paper and the proper derivation of the F-QCLE should lay the basis for further improvements of Floquet surface hopping.

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

A diagrammatic approach for automatically deriving analytical gradients of tensor hyper-contracted electronic structure methods

In this work, we introduce a diagrammatic approach to facilitate the automatic derivation of analytical nuclear gradients for tensor hyper-contraction (THC) based electronic structure methods. The automatically derived gradients are guaranteed to have the same scaling in terms of both operation count and memory footprint as the underlying energy calculations, and the computation of a gradient is roughly three times as costly as the underlying energy. The new diagrammatic approach enables the first cubic scaling implementation of nuclear derivatives for THC tensors fitted in molecular orbital basis (MO-THC). Furthermore, application of this new approach to THC-MP2 analytical gradients leads to an implementation, which is at least four times faster than the previously reported, manually derived implementation. Finally, we apply the new approach to the 14 tensor contraction patterns appearing in the supporting subspace formulation of multireference perturbation theory, laying the foundation for developments of analytical nuclear gradients and nonadiabatic coupling vectors for multi-state CASPT2.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗