SEARCH · Search NASA
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.
The existence of gaseous BeH 2 , predicted from simple valence theory, has still not been established experimentally. We have estimated the total energy of this molecule, assumed linear symmetric, as the sum of four negative terms; the calculated SCF energy, the SCF error to the Hartree-Fock limit, the correlation energy to the exact nonrelativistic energy, and the relativistic energy.
The packing geometry of macromolecules in complex mesophases is of key importance to self-organization in synthetic and biological soft materials. While approximate or heuristic models rely on often-untested assumptions about how flexible molecules “fit in” to distinct locations of complex assemblies, physical assemblies derive from ensembles of fluctuating conformations, obscuring the connection between mesophase geometry and the underlying arrangements. Here, we present an approach to extract and analyze features of molecular packing in diblock block copolymer (BCP) melts, a prototypical soft matter system, based on the statistical description of chain conformations in self-consistent field (SCF) theory. We show how average BCP chain trajectories in ordered morphologies can be analyzed from the SCF-derived orientational order parameter of chain segments. We use these extracted trajectories to analyze the features of local packing geometry, including chain bending and tilt, as well as the terminal boundaries that delineate distinct domains in ordered BCP morphologies. In conclusion, we illustrate this analysis by focusing on measurable features of packing frustration in 2D (columnar) and 3D (spherical and bicontinuous) morphologies, notably establishing an explicit link between chain conformations in complex morphologies and their medial geometry.
Here, density functional theory (DFT) is routinely used to make electronic structure predictions for high-throughput screening of materials and molecules for technologically relevant areas, like the identification of better catalysts, electronic materials, and drug discovery. However, the DFT formalism is limited by (a) its poor (quadratic-to-quartic) scaling, and (b) the need to perform repeated eigenvalue computations of the electronic Hamiltonian as part of its self-consistent field (SCF) iteration procedure to obtain the converged ground state electron density, ρ (r). Approaches that directly predict ρ (r) of a structure with high accuracy can accelerate conventional SCF calculations and can also be used in linearly scaling methods such as orbital-free DFT. To this end, we present a procedure to predict the ground state electron density of molecular and periodic three-dimensional systems directly from the atomic structure with a particular emphasis on physical interpretability. In our framework, ρ (r) is modeled using many-body correlation descriptors that accurately capture the effects of local atomic arrangements in the neighborhood of a grid point. Our use of a linear regression scheme to fit to charge density data enables transparent analysis of the relative contributions of various types of local atomic correlations. By systematically including increasingly complex correlations, our model is shown to accurately predict ρ (r) for a variety of chemically and electronically diverse systems — amorphous Ge, Al(001) slab, crystalline Ga 2 O 3 , molecular benzene, and polyethylene. We then demonstrate a symbolic regression-based protocol to construct easily computable, interpretable features from lower-order correlations that significantly improves our electron density predictions with effectively no increase in the computational cost.
In this work, we present a new stochastic extended Lagrangian molecular dynamics solution to charge equilibration that eliminates self-consistent field (SCF) calculations, thus eliminating the computational bottleneck in solving the charge distribution with standard SCF solvers. By formulating both charges and chemical potential as latent variables and introducing a holonomic constraint that satisfies charge conservation, the SC-XLMD method accurately reproduces thermodynamic, dynamic, and structural properties within the framework of ReaxFF for a bulk water system and highly reactive RDX molecules simulated at high temperature. The SC-XLMD method shows excellent computational performance and is available in the publicly available LAMMPS package.
We present the tensor hypercontraction difference self-consistent field (SCF) method, an approach that reduces the formal computational scaling of traditional naive self-consistent field methods from 𝑂(𝑁 4 ) to 𝑂(𝑁 3 ) with system size 𝑁. The scaling reduction is achieved by developing a new technique for constructing the tensor hypercontraction decomposition based on the fundamental approximation made in density fitting. Combining this scheme with the difference self-consistent field methodology, we achieve a method that enables 𝑂(𝑁 3 ) scaling SCF calculations with only 𝑂(𝑁 2 ) storage requirements. In conclusion, our proof-of-concept numerical tests demonstrate robust performance with errors in total energies below 8 × 10 –4 E h and with sub 1 kcal/mol errors for relative energies.
The photophysical and photochemical reactivity of the organometallic mono(cyclooctatetraenide)-Ce(III) complexes: [(C 8 H 8 )Ce(μ-X)(THF) 2 ] 2 X = O 3 SCF 3 - (1), Cl - (2), were studied. Here, in the course of these studies, a new polymorph of [(C 8 H 8 )Ce(μ-O 3 SCF 3 )(THF) 2 ] 2 (1) was described. Photoluminescence (PL) studies of 1 and 2 showed characteristic excitation and emission bands in the visible region, with bright green light emission under light irradiation from states corresponding to the 5d → 4f transition. Complexes 1 and 2 exhibit relatively long lifetime of 205.4 ± 0.2 and 145.8 ± 0.4 ns respectively. And the quantum yields (Φ) for 1 and 2 are 0.180 and 0.068 respectively. Electrochemical studies were performed on complexes 1 and 2 with a reversible Ce(III)/Ce(IV) redox couple recorded at $E_{1/2}$ = -1.53 V versus the Fc/Fc + for 2. Complex 1 shows an irreversible Ce(III/IV) oxidation wave. Complexes 1 and 2 revealed strongly-reducing, estimated, excited state reduction potentials ($E^*_{1/2}$) of -3.31 and -4.02 V versus the Fc/Fc + respectively, with small Stokes shifts of 0.12 eV. With their associated relatively long lifetimes, small Stokes shifts and large, negative $E^*_{1/2}$ values, both complexes were evaluated as potential photosensitisers for halogen atom transfer (XAT) using a test reaction of the dehalogenation of benzyl chloride.
The central importance of large-scale eigenvalue problems in scientific computation necessitates the development of massively parallel algorithms for their solution. Recent advances in dense numerical linear algebra have enabled the routine treatment of eigenvalue problems with dimensions on the order of hundreds of thousands on the world’s largest supercomputers. In cases where dense treatments are not feasible, Krylov subspace methods offer an attractive alternative due to the fact that they do not require storage of the problem matrices. However, demonstration of scalability of either of these classes of eigenvalue algorithms on computing architectures capable of expressing massive parallelism is non-trivial due to communication requirements and serial bottlenecks, respectively. In this work, we introduce the SISLICE method: a parallel shift-invert algorithm for the solution of the symmetric self-consistent field (SCF) eigenvalue problem. The SISLICE method drastically reduces the communication requirement of current parallel shift-invert eigenvalue algorithms through various shift selection and migration techniques based on density of states estimation and k-means clustering, respectively. This work demonstrates the robustness and parallel performance of the SISLICE method on a representative set of SCF eigenvalue problems and outlines research directions that will be explored in future work.
Accurate prediction of snow water equivalent (SWE) can be valuable for water resource managers. Recently, deep learning methods such as long short-term memory (LSTM) have exhibited high accuracy in simulating hydrologic variables and can integrate lagged observations to improve prediction, but their benefits were not clear for SWE simulations. Here we tested an LSTM network with data integration (DI) for SWE in the western United States to integrate 30-day-lagged or 7-day-lagged observations of either SWE or satellite-observed snow cover fraction (SCF) to improve future predictions. SCF proved beneficial only for shallow-snow sites during snowmelt, while lagged SWE integration significantly improved prediction accuracy for both shallow- and deep-snow sites. The median Nash–Sutcliffe model efficiency coefficient (NSE) in temporal testing improved from 0.92 to 0.97 with 30-day-lagged SWE integration, and root-mean-square error (RMSE) and the difference between estimated and observed peak SWE values d max were reduced by 41% and 57%, respectively. DI effectively mitigated accumulated model and forcing errors that would otherwise be persistent. Moreover, by applying DI to different observations (30-day-lagged, 7-day-lagged), we revealed the spatial distribution of errors with different persistent lengths. For example, integrating 30-day-lagged SWE was ineffective for ephemeral snow sites in the southwestern United States, but significantly reduced monthly-scale biases for regions with stable seasonal snowpack such as high-elevation sites in California. These biases are likely attributable to large interannual variability in snowfall or site-specific snow redistribution patterns that can accumulate to impactful levels over time for nonephemeral sites. These results set up benchmark levels and provide guidance for future model improvement strategies.
Understanding how host gene regulation responds to viral infection is essential for developing effective antiviral strategies. Emerging evidence suggests that host transcripts undergo dynamic chemical modifications to counteract viral invasion. Conversely, viruses that rely on nuclear transcription exploit host RNA methyltransferases to enhance mRNA export and translation. Orthopoxviruses, however, complete their entire replication cycle within compartmentalized cytoplasmic “factories” utilizing enzymes encoded by their large double-stranded viral DNA genomes. The dynamic interplay between host and poxviral epitranscriptome remains poorly characterized. Using a temporally resolved model of Vaccinia virus (VV) infection, we investigated host-virus interactions through transcriptome and N6-methyladenosine (m⁶A) epitranscriptome whole genome sequencing. We found that host m⁶A modifications respond rapidly to VV infection, preceding the delayed transcriptional changes that emerge at later stages. Early m⁶A signatures included key innate immunity factors as well as host genes involved in transcriptional regulation, post-transcriptional modification, and protein ubiquitination. Functional assays validated two host factors with early m⁶A modification changes that are essential for VV infection: a m⁶A reader, YTHDF1, and a component of the SCF E3 ubiquitin ligase complex, FBXO31. The m⁶A gain on YTHDF1 enhanced its protein expression and promoted efficient VV replication. In addition, we identified previously unrecognized roles of FBXO31 and the SCF E3 ligase complex in supporting VV infection. Temporal profiling of the m⁶A epitranscriptome reveals how VV exploits host post-transcriptional regulatory pathways, specifically m⁶A RNA modification and protein ubiquitination. These findings highlight critical host factors co-opted during poxvirus infection and identify potential targets for therapeutic intervention.
A numerical-variational method for performing self-consistent molecular calculations in the Hartree-Fock-Slater (HFS) model is presented. Molecular wavefunctions are expanded in terms of basis sets constructed from numerical HFS solutions of selected one-center atomlike problems. Binding energies and wavefunctions for the molecules are generated using a discrete variational method for a given molecular potential. In the self-consistent-charge (SCC) approximation to the complete self-consistent-field (SCF) method, results of a Mulliken population analysis of the molecular eigenfunctions are used in each iteration to produce 'atomic' occupation numbers. The simplest SCC potential is then obtained from overlapping spherical atomlike charge distributions. Molecular ionization energies are calculated using the transition-state procedure; results are given for CO, H2O, H2S, AlCl, InCl, and the Ni5O surface complex. Agreement between experimental and theoretical ionization energies for the free-molecule valence levels is generally within 1 eV. The simple SCC procedure gives a reasonably good approximation to the molecular potential, as shown by comparison with experiment, and with complete SCF calculations for CO, H2O, and H2S.
It is noted that for molecular calculations, the basis sets presented by Huzinaga (1977) need to be augmented by (1) 5p functions to describe the 5s - 5p near degeneracy; (2) a diffuse 4d function to provide for a balanced description of the 5s2 4dn, 5s1 4d(n+1), and 4d(n+2) states of the atom; and (3) a set of 4f functions to correlate the 4d functions. Here, the diffuse 4d function is similar in function to the diffuse 3d function for the first transition row elements recommended by Hay (1977). A table is included giving the optimized values for the diffuse 4d, the 5p, and 4f (STO exponent) functions. The diffuse 4d function and the 5p functions are optimized at the SCF level on the basis of the 5s1 4d(n+1) state (except for Pd, which is optimized for the 4d10 state) and the 5s1 5p1 4dn state, respectively. The table also gives the energies and the atomic symmetries for each of the SCF calculations.
The nature of the bonding of CO and NH3 ligands to Al is analyzed, and the intra-unit charge polarization and inter-unit donation for the interaction of ligands with metals are studied. The consequences of metal-to-ligand and ligand-to-metal charge transfer are separately considered by performing a constrained space orbital variation (CSOV) with the electrons of the metal member of the complex in the field of frozen ligand. The electrons of the metal atoms are then frozen in the relaxed distribution given by the CSOV SCF wave function and the ligand electrons are allowed to relax. Quantitative measures of the importance of inter-unit charge transfers and intra-unit polarization are obtained using results of SCF studies of Al4CO and Al4NH3 clusters chosen to simulate the adsorption of the ligands at an on-top side of the Al(111) surface. The electrostatic attraction of the effective dipole moments of the metal and ligand units makes an important contribution to the bond.
Research conducted during this period was focused on two main areas: (1) bonding in transition metal oxides; and (2) adsorption of CO on Al and Ni. In both of these theoretical studies a major interest was to obtain a better understanding of the nature of the bonding in transition metal containing systems. The studies used self consistent field (SCF), multi-configuration self cosistent field (MCSCF) and configuration interaction (CI) methods in the treatment of the transition metal oxides and only the SCF method in the adsorption studies. The reports of three principle investigators who contributed to this work during the tenure of the project are presented along with associated published papers.
Full configuration-interaction (FCI) calculations are performed at selected geometries for the 1-sigma(+) state of HF and the 2-B(1) and 2-A(1) states of NH2 using both DZ and DZP gaussian basis sets. Higher excitations become more important when the bonds are stretched and the self-consistent field (SCF) reference becomes a poorer zeroth-order description of the wave function. The complete active space SCF - multireference configuration-interaction (CASSCF-MRCI) procedure gives excellent agreement with the FCI potentials, especially when corrected with a multi-reference analog of the Davidson correction.
The present paper is concerned with an analysis of the bonding in Fe2(CO)9, taking into account the results of self-consistent-field (SCF) calculations. The basis sets used are the same as those employed in the study of Fe(CO)5 conducted by Bauschlicher and Bagus (1985). The Fe basis set starts with the 14s9p5d primitive set optimized by Wachters (1970), which is contracted to 8s4p5d. An SCF wave function, consisting of only the core electrons, was determined, and then electrons were slowly added to the virtual orbitals with the lowest eigenvalues. The conducted calculations show no direct Fe-Fe bonding, but rather bonding resulting from the bridging CO's. The bonding is found to involve both the 5 sigma and 2 pi(asterisk) orbitals of the bridging CO's.
Ab initio calculations at the self-consistent-field (SCF), singles-plus-doubles configuration-interaction (SDCI), and coupled-pair functional (CPF) level, are reported for the dipole moments and dipole derivatives of the X2Sigma(+) ground states of BeF, BeCl, MgF, MgCl, CaF, CaCl, and SrF. For comparison, analogous calculations are performed for the X1Sigma(+) state of KCl. The CPF results are found to be in remarkably better agreement with experiment than are the SCF and SDCI results. Apparently higher excitations are required to properly describe the radial extent along the bond axis of the remaining valence electron on the alkaline-earth metal.
The procedure employing bond functions recently suggested by Wright and Buenker has been applied to the N2 X 1 Sigma g + potential curve within the CAS SCF + MRSD CI treatment of electron correlation. The basis set used herein is identical to that employed by these authors in their SCF + CI calculations. The De and and the shape of the resulting potential curve, as judged by the computed vibrational levels, is not so accurate as would be expected from the results reported by Wright and Buenker (1984). The results indicate that using the CI superposition errors associated with bond functions to cancel basis set incompleteness depends on the treatment of the electron correlation.
Previously reported low-temperature pressure-broadening calculations (Green, 1985) for CO-He interacting via an SCF-CI potential are compared with new calculations in which the attractive part of the potential is either reduced by half or eliminated entirely. Results demonstrate that the attractive well is responsible for low-temperature enhancement of pressure-broadening cross sections and suggest that agreement with recent experimental values at 4 K (Messer and DeLucia, 1984) can be obtained by a modest reduction, probably within the expected uncertainty, in the attractive part of the SCF-CI potential.