Search NASA⌕ Search

SEARCH · Search NASA

Results for “SCF”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Generalized spin σ -SCF method

We introduce a generalization of the σ-SCF method to approximate noncollinear spin ground and excited single-reference electronic states by minimizing the Hamiltonian variance. The new method is based on the σ-SCF method, originally proposed by Ye et al. [J. Chem. Phys. 147, 214104 (2017)], and provides a prescription to determine ground and excited noncollinear spin states on an equal footing. Our implementation was carried out utilizing an initial simulated annealing stage followed by a mean-field iterative self-consistent approach to simplify the cumbersome search introduced by generalizing the spin degrees of freedom. The simulated annealing stage ensures a broad exploration of the Hilbert space spanned by the generalized spin single-reference states with random complex element-wise rotations of the generalized density matrix elements in the simulated annealing stage. The mean-field iterative self-consistent stage employs an effective Fockian derived from the variance, which is utilized to converge tightly to the solutions. This process helps us to easily find complex spin structures, avoiding manipulating the initial guess. As proof-of-concept tests, we present results for Hn (n = 3–7) planar rings and polyhedral clusters with geometrical spin frustration. We show that most of these systems have noncollinear spin excited states that can be interpreted in terms of geometric spin frustration. These states are not directly targeted by energy minimization methods, which are meant to converge to the ground state. This stresses the capability of the σ-SCF methodology to find approximate noncollinear spin structures as mean-field excited states.

Chemistry↗
DOE OSTI2026

SCF Framework, HF Stability, and RPA Correlation for Jordan–Wigner-Transformed Spin Hamiltonians on Arbitrary Coupling Topologies

Mapping spins to fermions via the Jordan–Wigner (JW) transformation can render mean-field (Hartree–Fock, HF) descriptions effective for strongly correlated spin systems. As established in recent work, the application of such approaches is not limited by the nonlocal structure of JW strings or by site ordering because string operators can be absorbed into Thouless rotations of a Slater determinant, and the variational optimization of a unitary Lie-algebraic similarity transformation removes any ordering dependence. Leveraging these ideas, we develop a self-consistent field (SCF) scheme that expresses the mean-field energy as a functional of the single-particle density matrix, providing an alternative to gradient-based optimization of Thouless parameters. We derive the analytical orbital Hessian to diagnose HF stability and compute the ground-state correlation energy through the random-phase approximation (RPA). Benchmark results for the XXZ and J 1 –J 2 model on one- and two-dimensional lattices demonstrate that RPA significantly improves mean-field accuracy.

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

Correction to “Frank–Kasper Phases of Diblock Copolymer Melts Studied with the DPD Model: SCF Results”

In a recent study, (1) we performed the self-consistent field (SCF) calculations of a dissipative particle dynamics (DPD) model of conformationally asymmetric diblock copolymer A-B melts to study the stability of various Frank–Kasper (FK) phases formed by such A-B melts. While all the equations in ref (1) are correct, we did not calculate the dimensionless (mean-field) Helmholtz free energies per chain βf c = βu c,χ + βu c,κ – s c /k B and the dimensionless stresses dβf c /dθ of ordered phases accurately. For example, the dimensionless (nonbonded) internal energy per chain due to the system compressibility βu c,κ = (N/2κ)∑ q βû 0 (|q|)(ϕ̂ A (q) + ϕ̂ B (q) – δ q,0 )(ϕ̂ A (−q) + ϕ̂ B (−q) – δ– q,0 ) should be calculated in the reciprocal space as given (note that βû 0 (|q|) is obtained analytically, and that since ϕ A (r) is real 𝜙̂ A (−𝐪)= $\overline{𝜙̂ A}$(𝐪), i.e., the complex conjugate of ϕ̂ A (q)); this corresponds to numerically evaluating βu c,κ = (N/2κ)∫dr(ϕ A (r) + ϕ B (r) – 1)ψ(r) with ψ(r) ≡ (∫dr′βu 0 (|r–r′|)(ϕ A (r′) + ϕ B (r′) – 1)/V in the real space via the (composite) trapezoidal rule. While the trapezoidal rule is less accurate than the Romberg integration (2) in general, the opposite occurs when the integration domain is periodic; this was pointed out by Matsen in ref (3) and explained in detail in refs (4and5). With our ignorance of this exceptional convergence property of the trapezoidal rule for periodic functions, we calculated all convolutions having the form of ψ(r) in the reciprocal space via the fast Fourier transforms (6) but those having the form of βu c,κ in the real space via Romberg integration; the latter leads to inaccurate results and unfortunately incorrect phase diagrams reported in ref (1), which are corrected here.

He, Juntong↗
DOE OSTI2021

CAS without SCF—Why to use CASCI and where to get the orbitals

The complete active space self-consistent field (CASSCF) method has seen broad adoption due to its ability to describe the electronic structure of both the ground and excited states of molecules over a broader swath of the potential energy surface than is possible with the simpler Hartree–Fock approximation. However, it also has a reputation for being unwieldy, computationally costly, and un-black-box. Here, we discuss a class of alternatives, complete active space configuration interaction (CASCI) methods, paying particular attention to their application to electronic excited states. The goal of this Perspective is fourfold. First, we argue that CASCI is not merely an approximation to CASSCF, in that it can be designed to have important qualitative advantages over CASSCF. Second, we present several insights drawn from our experience experimenting with different schemes for computing orbitals to be employed in CASCI. Third, we argue that CASCI is well suited for application to nanomaterials. Lastly, we reason that, with the rise in new low-scaling approaches for describing multireference systems, there is a greater need than ever to develop new methods for defining orbitals that provide an efficient and accurate description of both static correlation and electronic excitations in a limited active space.

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

C-terminal conformational changes in SCF-D3/MAX2 ubiquitin ligase are required for KAI2-mediated signaling

Karrikins (KARs) are bioactive molecules derived from burning vegetation. Plants perceive KARs through the hydrolase KARRIKIN INSENSITIVE 2 (KAI2) that interacts with the F-box protein ubiquitin ligase MORE AXILLARY GROWTH 2 (MAX2). MAX2 also plays a role in the perception and signal activation by Strigolactone (SL), a phytohormone controlling various developmental processes in plants. SL also acts as a rhizosphere signal to activate arbuscular mycorrhiza fungi that can be exploited by parasitic plants. kai2 knockouts exhibit distinct developmental defects and therefore KAI2 is hypothesized to perceive an unidentified endogenous ligand provisionally termed KAI2-Ligand (KL). Upon KAR/KL perception, the protein complex of KAI2-MAX2 targets SUPPRESSOR OF MAX2-1/2 (SMAX1)/SMXL2 for proteasomal degradation. Despite the identification of the key components KAI2, MAX2, and SMAX1 in KAR/KL signaling, their mode of interaction and regulation remains elusive. Recently, the regulatory function of the conformational switch of MAX2 C-terminal helix (CTH) in SL signaling has been demonstrated however its role in KAR/KL signaling remained unknown. Here we address the function of MAX2-CTH dynamics both in vitro and in planta and show that the central role of CTH is conserved between SL and KAR/KL signaling pathway.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Pseudodiagonalization Method for Accelerating Nonlinear Subspace Diagonalization in Density Functional Theory

In density functional theory, each self-consistent field (SCF) nonlinear step updates the discretized Kohn-Sham orbitals by solving a linear eigenvalue problem. The concept of pseudodiagonalization is to solve this linear eigenvalue problem approximately, and specifically utilizing a method involving a small number of Jacobi rotations that takes advantage of the good initial guess to the solution given by the approximation to the orbitals from the previous SCF iteration. The approximate solution to the linear eigenvalue problem can be very rapid, particularly for those steps near SCF convergence. Here, we adapt pseudodiagonalization to finite-temperature and metallic systems, where partially-occupied orbitals must be individually resolved with some accuracy. We apply pseudodiagonalization to the subspace eigenvalue problem that arises in Chebyshev-filtered subspace iteration. In tests on metallic and other systems for a range of temperatures, we show that pseudodiagonalization achieves similar rates of SCF convergence to exact diagonalization.

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

Revisiting Artifacts of Kohn–Sham Density Functionals for Biosimulation

We revisit the problem of unphysical charge density delocalization/fractionalization induced by the self-interaction error of common approximate Kohn–Sham (KS) density functional theory functionals on simulation of small to medium-sized proteins in a vacuum. Aside from producing unphysical electron densities and total energies, the vanishing of the HOMO–LUMO gap associated with the unphysical charge delocalization leads to an unphysical low-energy spectrum and catastrophic failure of most popular solvers for the KS self-consistent field (SCF) problem. We apply a robust quasi-Newton SCF solver to obtain solutions for some of these difficult cases. The anatomy of the charge delocalization is revealed by the natural deformation orbitals obtained from the density matrix difference between the Hartree–Fock and KS solutions; the charge delocalization not only can occur between charged fragments (such as in zwitterionic polypeptides) but also involves neutral fragments. The vanishing-gap phenomenon and troublesome SCF convergence are both attributed to the unphysical KS Fock operator eigenspectra of molecular fragments (e.g., amino acids or their side chains). Analysis of amino acid pairs suggests that the unphysical charge delocalization can be partially ameliorated by the use of some range-separated hybrid functionals but not by semilocal or standard hybrid functionals. Last, we demonstrate that solutions without the unphysical charge delocalization can be located even for semilocal KS functionals highly prone to such defects, but such solutions have non-Aufbau character and are unstable with respect to mixing of the non-overlapping “frontier” orbitals. Caution should be exercised when unexpectedly small (or vanishing) HOMO–LUMO gaps and atypical SCF convergence patterns (e.g., oscillatory) are observed in KS DFT simulations in any context (bio or otherwise).

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

Structural basis for dimerization quality control

Most quality control pathways target misfolded proteins to prevent toxic aggregation and neurodegeneration 1 . Dimerization quality control further improves proteostasis by eliminating complexes of aberrant composition 2 , but how it detects incorrect subunits remains unknown. Here we provide structural insight into target selection by SCF-FBXL17, a dimerization-quality-control E3 ligase that ubiquitylates and helps to degrade inactive heterodimers of BTB proteins while sparing functional homodimers. We find that SCF-FBXL17 disrupts aberrant BTB dimers that fail to stabilize an intermolecular β-sheet around a highly divergent β-strand of the BTB domain. Furthermore, complex dissociation allows SCF-FBXL17 to wrap around a single BTB domain, resulting in robust ubiquitylation. SCF-FBXL17 therefore probes both shape and complementarity of BTB domains, a mechanism that is well suited to establish quality control of complex composition for recurrent interaction modules.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Generating accurate density matrices on the tangent space of a Grassmann manifold

Interpolating a density matrix from a set of known density matrices is not a trivial task. This is because a linear combination of density matrices does not necessarily correspond to another density matrix. In this Communication, density matrices are examined as objects of a Grassmann manifold. Although this manifold is not a vector space, its tangent space is a vector space. As a result, one can map the density matrices on this manifold to their corresponding vectors in the tangent space and then perform interpolations on that tangent space. The resulting interpolated vector can be mapped back to the Grassmann manifold, which can then be utilized (1) as an optimal initial guess for a self-consistent field (SCF) calculation or (2) to derive energy directly without time-consuming SCF iterations. Such a promising approach is denoted as Grassmann interpolation (G-Int). The hydrogen molecule has been used to illustrate that the described interpolated method in this work preserves the essential attributes of a density matrix. For phosphorus mononitride and ferrocene, it was demonstrated numerically that reference points for the definition of the corresponding tangent spaces can be chosen arbitrarily. In addition, the interpolated density matrices provide a superior and essentially converged initial guess for an SCF calculation to make the SCF procedure itself unnecessary. Finally, this accurate, efficient, robust, and systematically improved G-Int strategy has been used for the first time to generate highly accurate potential energy surfaces with fine details for the difficult case, ferrocene.

Chemistry↗
DOE OSTI2023

The Grassmann interpolation method for spin-unrestricted open-shell systems

The recently reported Grassmann interpolation (G-Int) method [J. A. Tan and K. U. Lao, J. Chem. Phys. 158, 051101 (2023)] has been extended to spin-unrestricted open-shell systems. In contrast to closed-shell systems, where G-Int has to be performed only once since the α and β density matrices are the same, spin-unrestricted open-shell systems require G-Int to be performed twice—one for the α spin and another for the β spin density matrix. In this work, we tested the performance of G-Int to the carbon monoxide radical cation CO●+ and nickelocene complex, which have the doublet and triple ground states, respectively. We found that the Frobenius norm errors associated with the interpolations for the α and β spin density matrices are comparable for a given molecular geometry. These G-Int density matrices, when used as an initial guess for a self-consistent field (SCF) calculation, outperform the conventional SCF guess schemes, such as the superposition of atomic densities, purified superposition of atomic densities, core Hamiltonian, and generalized Wolfsberg–Helmholtz approximation. Depending on the desired accuracy, these G-Int density matrices can be used to directly evaluate the SCF energy without performing SCF iterations. In addition, the spin-unrestricted G-Int density matrices have been used for the first time to directly calculate the atomic charges using the Mulliken and ChElPG population analysis.

Chemistry↗
DOE OSTI2022

Evaluation of digestively resistant or soluble fibers, short- and medium-chain fatty acids, trace minerals, and antibiotics in nonchallenged nursery pigs on performance, digestibility, and intestinal integrity

Abstract Three experiments (EXP) were conducted to determine the effect of feed additives on performance, intestinal integrity, gastrointestinal volatile fatty acids (VFA), and energy and nutrient digestion in nonchallenged nursery pigs. In EXP 1, 480 pigs (6.36-kg body weight, BW) were placed into 96 pens with 5 pigs/pen, and allotted to 1 of 10 dietary treatments: 1) negative control containing no feed additive (NC), 2) NC + 44 mg chlortetracycline and 38.5 mg tiamulin/kg diet (CTsb), 3) NC + 5% resistant potato starch (RSpo), 4) NC + 5% soluble corn fiber (SCF), 5) NC + 5% sugar beet pulp (SBP), 6) NC + 0.30% fatty acid mix (FAM), 7) NC + 0.10% phytogenic blend of essential oils and flavoring compounds (PHY), 8) NC + 50 mg Cu and 1,600 mg zinc oxide/kg diet (CuZn), 9) NC + 5% resistant corn starch (RScn), and 10) NC + 0.05% β-glucan (BG) for 28 d. There was no impact of dietary treatment on BW gain or feed intake (P ≥ 0.22). Pigs fed diets containing SCF, CTsb, and RSpo resulted in microbial community differences compared to pigs fed the NC (P < 0.05). In EXP 2, 48 barrows (12.8 kg BW) were selected at the end of EXP 1 and fed the same dietary treatments they had previously received: 1) NC, 2) NC + 5% RScn, 3) NC + 5% SCF, and 4) NC + FAM for 8 d. There was no effect of feeding diets containing RScn, SCF, or FAM on in vivo intestinal permeability (P ≤ 0.21). Ileal or colon pH, concentrations of VFA did not differ due to dietary treatment (P ≥ 0.36), but pigs fed diets containing FAM resulted in a greater butyric acid concentration in the cecum compared to pigs fed the NC (P ≤ 0.05). In EXP 3, 156 pigs (6.11 kg BW) were placed into 52 pens with 3 pigs/pen and allotted to 1 of 4 dietary treatments arranged in a factorial manner: 1) NC, 2) NC + 5% RSpo, 3) NC + 0.30% FAM, and 4) NC + 5% RSpo + 0.30% FAM for 24 d. Feeding pigs diets containing RSpo did not affect BW gain (P = 0.91) while pigs fed diets containing FAM grew improved BW gain (P = 0.09). Colonic butyric acid concentrations were greater in pigs fed diets containing RSpo (P = 0.03), while pigs fed diets containing FAM exhibited reduced total VFA concentrations (P = 0.11). The results indicate that supplementing diets with digestively resistant but fermentable fibers, short- and medium-chain fatty acids, or antibiotics do not have a consistent effect, positive or negative, on markers of intestinal integrity or barrier function, intestinal VFA patterns, ATTD of energy and nutrients, or on pig performance.

Agriculture↗
DOE OSTI2024

Frank–Kasper Phases of Diblock Copolymer Melts: Self-Consistent Field Results of Two Commonly Used Models

We constructed phase diagrams of conformationally asymmetric diblock copolymer A-B melts using the polymer self-consistent field (SCF) calculations of both the dissipative particle dynamics chain (DPDC) model (i.e., compressible melts of discrete Gaussian chains with the DPD non-bonded potential) and the “standard” model (i.e., incompressible melts of continuous Gaussian chains with the Dirac δ-function non-bonded potential) in the χN-ε plane, where χN and ε characterize, respectively, the repulsion and conformational asymmetry between the A and B blocks, at the A-block volume fraction f = 0.2 and 0.3. Consistent with previous SCF calculations of the “standard” model, σ and A15 are the only stable Frank–Kasper (FK) phases among the five FK (i.e., σ, A15, C14, C15 and Z) phases considered. The stability of σ and A15 is due to their delicate balance between the energetic and entropic contributions to the Helmholtz free energy per chain of the system, which, within our parameter range, increases in the order of σ/A15, Z, and C14/C15. While in general the SCF phase diagrams of these two models are qualitatively consistent, A15 is not stable for the DPDC model at the copolymer chain length N = 10 and f = 0.3; any differences in the SCF phase diagrams are solely due to the differences between these two models.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗