Search NASA⌕ Search

SEARCH · Search NASA

Results for “talbot”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

An Algorithm for Atom-Centered Lossy Compression of the Atomic Orbital Basis in Density Functional Theory Calculations

Large atomic-orbital (AO) basis sets of at least triple and preferably quadruple-ζ (QZ) size are required to adequately converge Kohn–Sham density functional theory (DFT) calculations toward the complete basis set limit. However, incrementing the cardinal number by one nearly doubles the AO basis dimension, and the computational cost scales as the cube of the AO dimension, so this is very computationally demanding. Here, in this work, we develop and test a threshold-based natural atomic orbital (NAO) scheme in which ϵ-NAOs are obtained as eigenfunctions of atomic blocks of the density matrix in a one-center orthogonalized representation. This enables compression of the AO basis that is optimal for a given threshold, 10 –ϵ , by discarding NAOs with occupation numbers below that threshold. Extensive pilot test calculations using the Hartree–Fock functional and taking the converged density matrix as input suggest that a threshold of 10 –5 can yield a compression factor (ratio of AO to compressed ϵ-NAO dimension) between 2.5 and 4.5 for the QZ pc-3 basis. The errors in relative energies are typically less than 0.1 kcal/mol when the compressed basis is used instead of the uncompressed basis. Between 10 and 100 times smaller errors (i.e., usually less than 0.01 kcal/mol) can be obtained with a threshold 10 –7 , while the compression factor is typically between 2 and 2.5.

basis sets↗
DOE OSTI2024

The Role of Spin–Orbit Coupling in the Linear Absorption Spectrum and Intersystem Crossing Rate Coefficients of Ruthenium Polypyridyl Dyes

The successful use of molecular dyes for solar energy conversion requires efficient charge injection, which in turn requires the formation of states with sufficiently long lifetimes (e.g., triplets). The molecular structure elements that confer this property can be found empirically, however computational predictions using ab initio electronic structure methods are invaluable to identify structure-property relations for dye sensitizers. The primary challenge for simulations to elucidate the electronic and nuclear origins of these properties is a spin-orbit interaction which drives transitions between electronic states. Here, in this work, we present a computational analysis of the spin-orbit corrected linear absorption cross sections and intersystem crossing rate coefficients for a derivative set of phosphonated tris(2,2'-bipyridine)ruthenium(2+) dye molecules. After sampling the ground state vibrational distributions, the predicted linear absorption cross sections indicate that the mixture between singlet and triplet states plays a crucial role in defining the line shape of the metal-to-ligand charge transfer bands in these derivatives. Additionally, an analysis of the intersystem crossing rate coefficients suggests that transitions from the singlet into the triplet manifolds are ultrafast with rate coefficients on the order of 10 13 s -1 for each dye molecule.

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

Unusual Base Pair between Two 2-Thiouridines and Its Implication for Nonenzymatic RNA Copying

2-Thiouridine (s 2 U) is a nucleobase modification that confers enhanced efficiency and fidelity both on modern tRNA codon translation and on nonenzymatic and ribozyme-catalyzed RNA copying. We have discovered an unusual base pair between two 2-thiouridines that stabilizes an RNA duplex to a degree that is comparable to that of a native A:U base pair. High-resolution crystal structures indicate similar base-pairing geometry and stacking interactions in duplexes containing s 2 U:s 2 U compared to those with U:U pairs. Notably, the C=O···H–N hydrogen bond in the U:U pair is replaced with a C=S···H–N hydrogen bond in the s 2 U:s 2 U base pair. The thermodynamic stability of the s 2 U:s 2 U base pair suggested that this self-pairing might lead to an increased error frequency during nonenzymatic RNA copying. However, competition experiments show that s 2 U:s 2 U base-pairing induces only a low level of misincorporation during nonenzymatic RNA template copying because the correct A:s 2 U base pair outcompetes the slightly weaker s 2 U:s 2 U base pair. In addition, even if an s 2 U is incorrectly incorporated, the addition of the next base is greatly hindered. This strong stalling effect would further increase the effective fidelity of nonenzymatic RNA copying with s 2 U. Our findings suggest that s 2 U may enhance the rate and extent of nonenzymatic copying with only a minimal cost in fidelity.

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

Enhancing Hydrogen Production from Bioenergy Crops via Photoreforming

Photoreforming perennial bioenergy crops (willow, Miscanthus , and poplar) has the potential to produce H 2 with reduced environmental impacts. To understand the compositional effects of the biomass on the average rate of H 2 production over the first 30 min of reaction ( r H 2 ), the r H 2 values of model biomass component (i.e., cellulose, hemicellulose, and lignin) mixtures were compared with those from the raw biomass. The higher cellulose or hemicellulose content in multicomponent mixtures resulted in higher r H 2 , whereas lignin reduced the hydrogen production rate. However, with raw biomass, the ratio of biomass components alone did not determine the r H 2 via photoreforming, with rates of hydrogen production for different varieties of willow ranging between 1.9 μmol h -1 and 12.3 μmol h -1 , 11.8 μmol h -1 for a poplar, and 6.8 μmol h -1 for a miscanthus biomass. In addition, comparable r H 2 values of raw poplar and its extracted cellulose via an IonoSolv treatment indicated the possibility of using raw biomass materials without delignification for generating H 2 via photoreforming. Importantly, r H 2 was positively correlated with the interaction between water and the biomass, as assessed by NMR relaxation via an examination of the T 1 /T 2 ratio. A stronger water-biomass interaction resulted in a higher r H 2 . Genetic modification of biomass has been suggested as a putative way to improve the r H 2 of biomass with an enhanced interaction with water. This research enhances the understanding of factors influencing H 2 production from lignocellulosic biomass by photoreforming and supports the breeding and management of perennial biomass crops to maximize H 2 yields while minimizing land area requirements.

biomass↗
DOE OSTI2021

Non-structural carbohydrates mediate seasonal water stress across Amazon forests

Non-structural carbohydrates (NSC) are major substrates for plant metabolism and have been implicated in mediating drought-induced tree mortality. Despite their significance, NSC dynamics in tropical forests remain little studied. We present leaf and branch NSC data for 82 Amazon canopy tree species in six sites spanning a broad precipitation gradient. During the wet season, total NSC (NSC T ) concentrations in both organs were remarkably similar across communities. However, NSC T and its soluble sugar (SS) and starch components varied much more across sites during the dry season. Notably, the proportion of leaf NSC T in the form of SS (SS:NSC T ) increased greatly in the dry season in almost all species in the driest sites, implying an important role of SS in mediating water stress in these sites. This adjustment of leaf NSC balance was not observed in tree species less-adapted to water deficit, even under exceptionally dry conditions. Thus, leaf carbon metabolism may help to explain floristic sorting across water availability gradients in Amazonia and enable better prediction of forest responses to future climate change.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Ocean current patterns drive the worldwide colonization of eelgrass (Zostera marina)

Currents are unique drivers of oceanic phylogeography and thus determine the distribution of marine coastal species, along with past glaciations and sea-level changes. Here we reconstruct the worldwide colonization history of eelgrass (Zostera marina L.), the most widely distributed marine flowering plant or seagrass from its origin in the Northwest Pacific, based on nuclear and chloroplast genomes. We identified two divergent Pacific clades with evidence for admixture along the East Pacific coast. Two west-to-east (trans-Pacific) colonization events support the key role of the North Pacific Current. Time-calibrated nuclear and chloroplast phylogenies yielded concordant estimates of the arrival of Z. marina in the Atlantic through the Canadian Arctic, suggesting that eelgrass-based ecosystems, hotspots of biodiversity and carbon sequestration, have only been present there for ~243 ky (thousand years). Mediterranean populations were founded ~44 kya, while extant distributions along western and eastern Atlantic shores were founded at the end of the Last Glacial Maximum (~19 kya), with at least one major refuge being the North Carolina region. The recent colonization and five- to sevenfold lower genomic diversity of the Atlantic compared to the Pacific populations raises concern and opportunity about how Atlantic eelgrass might respond to rapidly warming coastal oceans.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Basin-wide variation in tree hydraulic safety margins predicts the carbon balance of Amazon forests

Abstract Tropical forests face increasing climate risk 1,2 , yet our ability to predict their response to climate change is limited by poor understanding of their resistance to water stress. Although xylem embolism resistance thresholds (for example, $$\varPsi $$ Ψ 50 ) and hydraulic safety margins (for example, HSM 50 ) are important predictors of drought-induced mortality risk 3–5 , little is known about how these vary across Earth’s largest tropical forest. Here, we present a pan-Amazon, fully standardized hydraulic traits dataset and use it to assess regional variation in drought sensitivity and hydraulic trait ability to predict species distributions and long-term forest biomass accumulation. Parameters $$\varPsi $$ Ψ 50 and HSM 50 vary markedly across the Amazon and are related to average long-term rainfall characteristics. Both $$\varPsi $$ Ψ 50 and HSM 50 influence the biogeographical distribution of Amazon tree species. However, HSM 50 was the only significant predictor of observed decadal-scale changes in forest biomass. Old-growth forests with wide HSM 50 are gaining more biomass than are low HSM 50 forests. We propose that this may be associated with a growth–mortality trade-off whereby trees in forests consisting of fast-growing species take greater hydraulic risks and face greater mortality risk. Moreover, in regions of more pronounced climatic change, we find evidence that forests are losing biomass, suggesting that species in these regions may be operating beyond their hydraulic limits. Continued climate change is likely to further reduce HSM 50 in the Amazon 6,7 , with strong implications for the Amazon carbon sink.

Science & Technology - Other Topics↗
DOE OSTI2021

Structure of the master regulator Rns reveals an inhibitor of enterotoxigenic Escherichia coli virulence regulons

Enteric infections caused by the gram-negative bacteria enterotoxigenic Escherichia coli (ETEC), Vibrio cholerae, Shigella flexneri, and Salmonella enterica are among the most common and affect billions of people each year. These bacteria control expression of virulence factors using a network of transcriptional regulators, some of which are modulated by small molecules as has been shown for ToxT, an AraC family member from V. cholerae. In ETEC the expression of many types of adhesive pili is dependent upon the AraC family member Rns. We present here the 3 Å crystal structure of Rns and show it closely resembles ToxT. Rns crystallized as a dimer via an interface similar to that observed in other dimeric AraC’s. Furthermore, the structure of Rns revealed the presence of a ligand, decanoic acid, that inhibits its activity in a manner similar to the fatty acid mediated inhibition observed for ToxT and the S. enterica homologue HilD. Together, these results support our hypothesis that fatty acids regulate virulence controlling AraC family members in a common manner across a number of enteric pathogens. Furthermore, for the first time this work identifies a small molecule capable of inhibiting the ETEC Rns regulon, providing a basis for development of therapeutics against this deadly human pathogen.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2022

Dynamic signatures of electronically nonadiabatic coupling in sodium hydride: a rigorous test for the symmetric quasi-classical model applied to realistic, ab initio electronic states in the adiabatic representation

Sodium hydride (NaH) in the gas phase presents a seemingly simple electronic structure making it a potentially tractable system for the detailed investigation of nonadiabatic molecular dynamics from both computational and experimental standpoints. The single vibrational degree of freedom, as well as the strong nonadiabatic coupling that arises from the excited electronic states taking on considerable ionic character, provides a realistic chemical system to test the accuracy of quasi-classical methods to model population dynamics where the results are directly comparable against quantum mechanical benchmarks. Here, using a simulated pump–probe type experiment, this work presents computational predictions of population transfer through the avoided crossings of NaH via symmetric quasi-classical Meyer–Miller (SQC/MM), Ehrenfest, and exact quantum dynamics on realistic, ab initio potential energy surfaces. The main driving force for population transfer arises from the ground vibrational level of the D 1 Σ + adiabatic state that is embedded in the manifold of near-dissociation C 1 Σ + vibrational states. When coupled through a sharply localized first-order derivative coupling most of the population transfers between t = 15 and t = 30 fs depending on the initially excited vibronic wavepacket. While quantum mechanical effects are expected due to the reduced mass of NaH, predictions of the population dynamics from both the SQC/MM and Ehrenfest models perform remarkably well against the quantum dynamics benchmark. Additionally, an analysis of the vibronic structure in the nonadiabatically coupled regime is presented using a variational eigensolver methodology.

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

Compositional and structural control in LLZO solid electrolytes

Garnet-based solid-state electrolytes (SSEs) represent a promising class of materials for next-generation batteries with improved safety and performance. However, lack of control over the composition and crystal structure of the well-known Li 7 La 3 Zr 2 O 12 (LLZO) garnet material has led to poor reproducibility with a wide range of ionic conductivities reported in the literature. In this study, the role of precursor homogeneity in controlling the compositional and structural evolution of Al-doped LLZO is explored. A novel solution-based synthesis approach is employed to demonstrate enhanced atomic-scale mixing of the starting materials in comparison to conventional solid-state preparation methods. Through this technique, it is shown that the stability and formation temperature of the highly conductive cubic phase is directly impacted by the spatial distribution of the doping element and reactant species in the precursor mixture. Precursor homogeneity was also an important factor in mitigating the formation of unwanted secondary impurities. These findings can be used to guide the synthesis of SSEs with reproducible material characteristics and enhanced electrolytic performance.

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

Accelerating and stabilizing the convergence of vibrational self-consistent field calculations via the direct inversion of the iterative subspace (vDIIS) algorithm

The vibrational self-consistent field (VSCF) method yields anharmonic states and spectra for molecular vibrations, and it serves as the starting point for more sophisticated correlated-vibration methods. Convergence of the iterative, non-linear optimization in VSCF calculations can be erratic or altogether unsuccessful, particularly for chemical systems involving low-frequency motions. In this work, a vibrational formulation of the Direct Inversion of the Iterative Subspace method of Pulay is presented and investigated. This formulation accounts for distinct attributes of the vibrational and electronic cases, including the expansion of each single-mode vibrational wavefunction in its own basis set. The resulting Direct Inversion of the Iterative Subspace method is shown to substantially accelerate VSCF convergence in all convergent cases as well as rectify many cases where Roothaan-based methods fail. Performance across systems ranging from small, rigid molecules to weakly bound molecular clusters is investigated in this analysis.

Chemistry↗
DOE OSTI2023

Pitfalls in the n -mode representation of vibrational potentials

Simulations of anharmonic vibrational motion rely on computationally expedient representations of the governing potential energy surface. The n-mode representation (n-MR)—effectively a many-body expansion in the space of molecular vibrations—is a general and efficient approach that is often used for this purpose in vibrational self-consistent field (VSCF) calculations and correlated analogues thereof. In the present analysis, a lack of convergence in many VSCF calculations is shown to originate from negative and unbound potentials at truncated orders of the n-MR expansion. For cases of strong anharmonic coupling between modes, the n-MR can both dip below the true global minimum of the potential surface and lead to effective single-mode potentials in VSCF that do not correspond to bound vibrational problems, even for bound total potentials. The present analysis serves mainly as a pathology report of this issue. Furthermore, this insight into the origin of VSCF non-convergence provides a simple, albeit ad hoc, route to correct the problem by “painting in” the full representation of groups of modes that exhibit these negative potentials at little additional computational cost. Somewhat surprisingly, this approach also reasonably approximates the results of the next-higher n-MR order and identifies groups of modes with particularly strong coupling. The method is shown to identify and correct problematic triples of modes—and restore SCF convergence—in two-mode representations of challenging test systems, including the water dimer and trimer, as well as protonated tropine.

Chemistry↗
DOE OSTI2022

The symmetric quasi-classical model using on-the-fly time-dependent density functional theory within the Tamm–Dancoff approximation

The primary computational challenge when simulating nonadiabatic ab initio molecular dynamics is the unfavourable compute costs of electronic structure calculations with molecular size. Simple electronic structure theories, like time-dependent density functional theory within the Tamm–Dancoff approximation (TDDFT/TDA), alleviate this cost for moderately sized molecular systems simulated on realistic time scales. Although TDDFT/TDA does have some limitations in accuracy, an appealing feature is that, in addition to including electron correlation through the use of a density functional, the cost of calculating analytic nuclear gradients and nonadiabatic coupling vectors is often computationally feasible even for moderately sized basis sets. Here in this work, some of the benefits and limitations of TDDFT/TDA are discussed and analysed with regard to its applicability as a ‘back-end’ electronic structure method for the symmetric quasi-classical Meyer–Miller model (SQC/MM). In order to investigate the benefits and limitations of TDDFT/TDA, SQC/MM is employed to predict and analyse a prototypical example of excited-state hydrogen transfer in gas-phase malonaldehyde. Then, the ring-opening dynamics of selenophene are simulated, which highlight some of the deficiencies of TDDFT/TDA. Additionally, some new algorithms are proposed that speed up the calculation of analytic nuclear gradients and nonadiabatic coupling vectors for a set of excited electronic states.

molecular dynamics↗
DOE OSTI2025

Hydrology controls thermokarst and alters carbon cycling and methane emissions in peatlands near the southern limit of permafrost

Permafrost peatlands store vast amounts of frozen carbon across northern landscapes. When ground ice melts, surface subsidence produces thermokarst landforms that expand wetlands at the edges of permafrost plateaus. Thermokarst represents an accelerating climate feedback, but uncertainties remain about how ground ice, hydrology, and vegetation interact to shape landscape change and carbon fluxes. We extended the process-based model ecosys to simulate thermokarst dynamics in laterally coupled 2D transects at a well-characterized boreal peatland site in Canada’s Northwest Territories. After benchmarking against site observations, we varied ground ice content and hydrologic boundary conditions across ranges typical near the southern permafrost limit. Simulations revealed distinct degradation regimes governed by the elevation difference between the frost table and the external water table. Rates of lateral retreat, the thaw-driven encroachment of wetlands into adjacent plateaus, ranged from 0 to >2 m yr −1 under identical weather forcing, consistent with observations and highlighting the strong role of WT and ground ice. Simulated vegetation dynamics indicate that black spruce mortality cannot be explained by anoxia alone, pointing to additional stressors such as root damage, pathogens, or physical destabilization. Despite large hydrologic shifts, net ecosystem CO 2 exchange remained a slight sink after collapse, while methane (CH 4 ) emissions rose by one to two orders of magnitude. As a result, lateral retreat substantially increases the greenhouse warming potential of permafrost peatlands (1.7 million km 2 in area), with simulated emissions of 0.1–10 Mt CO 2 -eq decade −1 depending on hydrology and retreat rates. These results underscore the need to account for both ground ice and hydrologic dynamics when assessing thermokarst-driven climate feedbacks.

carbon cycling↗