Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Density matrix methods”

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 1,171 records · Page 65
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Comparing the Effects of Nb, Pb, Y, and La Replacement on the Structural, Electrical, and Magnetic Characteristics of Bi-Based Superconductors

In this study, the effects of Pb, Nb, La, and Y replacements were investigated on Bi-based superconducting materials. In preparing the samples, we used a method called solid-state reaction method. The patterns of the X-ray diffraction of all samples indicated presence of Bi-2212 and Bi-2223 phases. The results obtained from XRD revealed that with increase of the melting point of substation elements, the Bi-2223 phase decreased while the Bi-2212 phase and impurity phases of samples grew. From the electrical resistivity measurements using the four-probe method, it was found that sample A with Pb and sample B with La replacements had the maximum and minimum critical temperatures of 111.4 and 81.6 K, respectively. Based on hysteresis loop (M–H) measurement using Bean’s model, estimation of critical current density (J{sub c}) showed that sample A with Pb and sample B with La substitution had the maximum and minimum values respectively. These results may be due to the melting point of these elements with values of 888, 1512, 2315, and 2425 {sup ∘}C for PbO, Nb{sub 2}O{sub 5}, La{sub 2}O{sub 3}, and Y{sub 2}O{sub 3}, respectively. These elements were replaced by Bi{sub 2}O{sub 3} with a melting point of 817 {sup ∘} C. Further, the samples were prepared at the temperature of 845 {sup ∘}C. It seems at this temperature, these elements not only dissolve within the main matrix and participate in the formation of the Bi-2212 phase during the sintering process but they also participate in the development of the variety of the impurity phases as confirmed by XRD results.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2013

Recent development of Monte Carlo shell model and its application to no-core calculations

One of the major challenges in nuclear theory is to reproduce and to predict nuclear structure from ab initio calculations with realistic nuclear forces. As the current limitation of direct diagonalization of Hamiltonian matrices by Lanczos iteration method is around the order of matrix dimensionality 10 10 in shell-model calculations, it is difficult to access heavier nuclei beyond the p shell with sufficiently large basis spaces. It is possible to overcome this difficulty by utilizing efficient approximate methods to reproduce full ab initio solutions with good precision and quantified uncertainties. Following the major success of the Monte Carlo shell model (MCSM) with an assumed inert core in the sd- and pf-shell regions and also by recent developments in the MCSM algorithm, the no-core MCSM is expected to be one of the most powerful tools to meet these conditions. We have performed benchmark calculations in the p-shell region. Results of energies are compared with those in the full configuration interaction and no-core full configuration methods. These are found to be consistent with each other within quoted uncertainties when they could be quantified. We also compare and discuss the radial density of the helium-4 ground state extracted from the MCSM and FCI many-body wave functions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS2012

NASA Tech Briefs, July 2012

Topics covered include: Instrument Suite for Vertical Characterization of the Ionosphere-Thermosphere System; Terahertz Radiation Heterodyne Detector Using Two-Dimensional Electron Gas in a GaN Heterostructure; Pattern Recognition Algorithm for High-Sensitivity Odorant Detection in Unknown Environments; Determining Performance Acceptability of Electrochemical Oxygen Sensors; Versatile Controller for Infrared Lamp and Heater Arrays; High-Speed Scanning Interferometer Using CMOS Image Sensor and FPGA Based on Multifrequency Phase-Tracking Detection; Ultra-Low-Power MEMS Selective Gas Sensors; Compact Receiver Front Ends for Submillimeter-Wave Applications; Dynamically Reconfigurable Systolic Array Accelerator; Blocking Losses With a Photon Counter; Motion-Capture-Enabled Software for Gestural Control of 3D Mod; Orbit Software Suite; CoNNeCT Baseband Processor Module Boot Code SoftWare (BCSW); Trajectory Software With Upper Atmosphere Model; ALSSAT Version 6.0; Employing a Grinding Technology to Assess the Microbial Density for Encapsulated Organisms; Demonstration of Minimally Machined Honeycomb Silicon Carbide Mirrors; Polyimide Aerogel Thin Films; Nanoengineered Thermal Materials Based on Carbon Nanotube Array Composites; Composite Laminate With Coefficient of Thermal Expansion Matching D263 Glass; Robust Tensioned Kevlar Suspension Design; Focal Plane Alignment Utilizing Optical CMM; Purifying, Separating, and Concentrating Cells From a Sample Low in Biomass; Virtual Ultrasound Guidance for Inexperienced Operators; Beat-to-Beat Blood Pressure Monitor; Non-Contact Conductivity Measurement for Automated Sample Processing Systems; An MSK Radar Waveform; Telescope Alignment From Sparsely Sampled Wavefront Measurements Over Pupil Subapertures; Method to Remove Particulate Matter from Dusty Gases at Low Pressures; Terahertz Quantum Cascade Laser With Efficient Coupling and Beam Profile; Measurement Via Optical Near-Nulling and Subaperture Stitching; 885-nm Pumped Ceramic Nd:YAG Master Oscillator Power Amplifier Laser System; Airborne Hyperspectral Imaging System; Heat Shield Employing Cured Thermal Protection Material Blocks Bonded in a Large-Cell Honeycomb Matrix; and Asymmetric Supercapacitor for Long-Duration Power Storage.

Source record↗
DOE OSTI2014

Daubechies wavelets for linear scaling density functional theory

We demonstrate that Daubechies wavelets can be used to construct a minimal set of optimized localized adaptively contracted basis functions in which the Kohn-Sham orbitals can be represented with an arbitrarily high, controllable precision. Ground state energies and the forces acting on the ions can be calculated in this basis with the same accuracy as if they were calculated directly in a Daubechies wavelets basis, provided that the amplitude of these adaptively contracted basis functions is sufficiently small on the surface of the localization region, which is guaranteed by the optimization procedure described in this work. This approach reduces the computational costs of density functional theory calculations, and can be combined with sparse matrix algebra to obtain linear scaling with respect to the number of electrons in the system. Calculations on systems of 10 000 atoms or more thus become feasible in a systematic basis set with moderate computational resources. Further computational savings can be achieved by exploiting the similarity of the adaptively contracted basis functions for closely related environments, e.g., in geometry optimizations or combined calculations of neutral and charged systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2016

Mesh Convergence Requirements for Composite Damage Models

The ability of the finite element method to accurately represent the response of objects with intricate geometry and loading renders the finite element method as an extremely versatile analysis technique for structural analysis. Finite element analysis is routinely used in industry to calculate deflections, stress concentrations, natural frequencies, buckling loads, and much more. The method works by discretizing complex problems into smaller, simpler approximations that are valid over small uniform domains. For common analyses, the maximum size of the elements that can be used is often be determined by experience. However, to verify the quality of a solution, analyses with several levels of mesh refinement should be performed to ensure that the solution has converged. In recent years, the finite element method has been used to calculate the resistance of structures, and in particular that of composite structures. A number of techniques such as cohesive zone modeling, the virtual crack closure technique, and continuum damage modeling have emerged that can be used to predict cracking, delaminations, fiber failure, and other composite damage modes that lead to structural collapse. However, damage models present mesh refinement requirements that are not well understood. In this presentation, we examine different mesh refinement issues related to the representation of damage in composite materials. Damage process zone sizes and their corresponding mesh requirements will be discussed. The difficulties of modeling discontinuities and the associated need for regularization techniques will be illustrated, and some unexpected element size constraints will be presented. Finally, some of the difficulties in constructing models of composite structures capable of predicting transverse matrix cracking will be discussed. It will be shown that to predict the initiation and propagation of transverse matrix cracks, their density, and their saturation may require models that are significantly more refined than those that have been contemplated in the past.

Davila, Carlos G.↗
DOE OSTI2013

Study on generating of thermal neutron scattering cross sections for LiH

LiH is designated as a promising moderator and shielding material because of its low density, high melting point and large fraction of H atoms. However, lack of the thermal neutron cross sections of LiH makes numerical calculation deviate from experimental data to some extent. As a result, it is necessary to study LiH thermal kernel effect. The phonon property of LiH has been investigated by first-principles calculations using the plane-wave pseudo potential method with CASTEP code. The scattering law and the thermal neutron scattering cross sections for Li and H have been generated using this distribution. The results have been compared with zirconium hydride data. The GASKET and NJOY/LEAPR codes have been used in the calculation of scattering law, whose results have been compared with the reference; the discrepancy mainly comes from phonon spectrums and its expansion. LEAPR had the capability to compute scattering through larger energy and momentum transfers than GASKET did. By studying LiH phonon spectrum and constructing the model of LiH thermal kernel and scattering matrix, the ACE format LiH thermal neutron cross sections for MCNP software could be made and used for reactor Neutronics calculation. (authors)

36 MATERIALS SCIENCE↗
DOE OSTI2024

Microstructural and Hall–Petch Analysis of Additively Manufactured Ferritic Alloy Using 2507 Duplex Stainless Steel Powder

The powder bed fusion–laser beam (PBF-LB) process, a method of additive manufacturing (AM), was used to print duplex stainless steel (DSS) using commercial-grade 2507 powders. While conventionally processed DSS has a two-phase microstructure consisting of 50% austenite and 50% ferrite, the PBF-LB-printed 2507 alloy was nearly 100% ferrite. Optimal processing conditions that minimized porosity were determined to be 290 W laser power and 1000 mm/s scan speed, and grain size, texture, and phases were characterized as a function of laser power and scan speed. Grain size increased with increasing laser power but decreased with increasing scan speed. A <100> texture diminished with increasing scan speed from 1000 mm/s to 1400 mm/s. No austenite phase was detected. Transmission electron microscopy (TEM) and scanning transmission electron microscopy (STEM) characterization revealed nanoscale chromium nitride precipitates in the ferritic matrix (incoherent hexagonal close-packed (HCP) precipitates at grain boundaries and coherent body-centered cubic (BCC) precipitates within the grains) and a high density of tangled dislocations. Tensile tests of as-printed alloys showed a yield strength of 570 MPa, an ultimate tensile strength of 756 MPa, and an elongation to failure of 10%. The tensile properties were analyzed based on the observed microstructure considering grain size, nanoscale precipitates, and the high density of dislocations.

Crystallography↗
DOE OSTI2022

Pairwise connected tensor network representation of path integrals

It has been recently shown how the tensorial nature of real-time path integrals (PIs) involving the Feynman-Vernon influence functional can be utilized with matrix product states, taking advantage of the finite length of the bath-induced memory. Tensor networks (TNs) promise to provide a unified language to express the structure of a PI. A generalized TN specifically incorporating the pairwise interaction structure of the influence functional and its invariance with respect to the average forward-backward position or the sojourn value in the form of the blip representation is derived and implemented. This pairwise connected TNPI (PC-TNPI) is illustrated through applications to typical spin-boson problems and explorations of the differences caused by the exact form of the spectral density. The storage and performance scalings are reported, showing the compactness of the representation and the efficiency of the contraction process. Finally, taking advantage of the compressed representation, the viability of using PC-TNPI for simulating multistate problems is demonstrated. The PC-TNPI structure can be shown to yield other TN algorithms currently in use. Consequently, it should be possible to use it as a starting point for deriving other optimized procedures.

36 MATERIALS SCIENCE↗
DOE OSTI2016

TU-FG-BRB-04: A New Optimization Method for Pre-Treatment Patient-Specific Stopping-Power by Combining Proton Radiography and X-Ray CT

Purpose: The relative stopping power (RSP) uncertainty is the largest contributor to the range uncertainty in proton therapy. The purpose of this work is to develop a robust and systematic method that yields accurate patient specific RSPs by combining pre-treatment X-ray CT and daily proton radiography. Methods: The method is formulated as a penalized least squares optimization (PLSO) problem min(|Ax-B|). The matrix A represents the cumulative path-length crossed in each material computed by calculating proton trajectories through the X-ray CT. The material RSPs are denoted by x and B is the pRad, expressed as water equivalent thickness. The equation is solved using a convex-conic optimizer. Geant4 simulations were made to assess the feasibility of the method. RSP extracted from the Geant4 materials were used as a reference and the clinical HU-RSP curve as a comparison. The PLSO was first tested on a Gammex RMI-467 phantom. Then, anthropomorphic phantoms of the head, pelvis and lung were studied and resulting RSPs were evaluated. A pencil beam was generated in each phantom to evaluate the proton range accuracy achievable by using the optimized RSPs. Finally, experimental data of a pediatric head phantom (CIRS) were acquired using a recently completed experimental pCT scanner. Results: Numerical simulations showed precise RSP (<0.75%) for Gammex materials except low-density lung (LN-300) (1.2%). Accurate RSP have been obtained for the head (µ=−0.10%, 1.5σ=1.12%), lung (µ=−0.33, 1.5σ=1.02%) and pelvis anthropomorphic phantoms (µ=0.12, 1.5σ=0,99%). The range precision has been improved with an average R80 difference to the reference (µ±1.5σ) of −0.20±0.35%, −0.47±0.92% and −0.06±0.17% in the head, lung and pelvis phantoms respectively, compared to the 3.5% clinical margin. Experimental HU-RSP curve have been produced on the CIRS pediatric head. Conclusion: The proposed PLSO with prior knowledge X-ray CT shows promising potential (R80 σ<1.0% over all sites) to decrease the range uncertainty.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2011

I-NERI-2007-004-K, Development and Characterization of New High-Level Waste Forms for Achieving Waste Minimization from Pyroprocessing

Work describe in this report represents the final year activities for the 3-year International Nuclear Energy Research Initiative (I-NERI) project: Development and Characterization of New High-Level Waste Forms for Achieving Waste Minimization from Pyroprocessing. Used electrorefiner salt that contained actinide chlorides and was highly loaded with surrogate fission products was processed into three candidate waste forms. The first waste form, a high-loaded ceramic waste form is a variant to the CWF produced during the treatment of Experimental Breeder Reactor-II used fuel at the Idaho National Laboratory (INL). The two other waste forms were developed by researchers at the Korean Atomic Energy Research Institute (KAERI). These materials are based on a silica-alumina-phosphate matrix and a zinc/titanium oxide matrix. The proposed waste forms, and the processes to fabricate them, were designed to immobilize spent electrorefiner chloride salts containing alkali, alkaline earth, lanthanide, and halide fission products that accumulate in the salt during the processing of used nuclear fuel. This aspect of the I-NERI project was to demonstrate 'hot cell' fabrication and characterization of the proposed waste forms. The outline of the report includes the processing of the spent electrorefiner salt and the fabrication of each of the three waste forms. Also described is the characterization of the waste forms, and chemical durability testing of the material. While waste form fabrication and sample preparation for characterization must be accomplished in a radiological hot cell facility due to hazardous radioactivity levels, smaller quantities of each waste form were removed from the hot cell to perform various analyses. Characterization included density measurement, elemental analysis, x-ray diffraction, scanning electron microscopy and the Product Consistency Test, which is a leaching method to measure chemical durability. Favorable results from this demonstration project will provide additional options for fission product immobilization and waste management associated the electrochemical/pyrometallurgical processing of used nuclear fuel.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2016

Chebyshev polynomial filtered subspace iteration in the discontinuous Galerkin method for large-scale electronic structure calculations

The Discontinuous Galerkin (DG) electronic structure method employs an adaptive local basis (ALB) set to solve the Kohn-Sham equations of density functional theory in a discontinuous Galerkin framework. The adaptive local basis is generated on-the-fly to capture the local material physics and can systematically attain chemical accuracy with only a few tens of degrees of freedom per atom. A central issue for large-scale calculations, however, is the computation of the electron density (and subsequently, ground state properties) from the discretized Hamiltonian in an efficient and scalable manner. We show in this work how Chebyshev polynomial filtered subspace iteration (CheFSI) can be used to address this issue and push the envelope in large-scale materials simulations in a discontinuous Galerkin framework. We describe how the subspace filtering steps can be performed in an efficient and scalable manner using a two-dimensional parallelization scheme, thanks to the orthogonality of the DG basis set and block-sparse structure of the DG Hamiltonian matrix. The on-the-fly nature of the ALB functions requires additional care in carrying out the subspace iterations. We demonstrate the parallel scalability of the DG-CheFSI approach in calculations of large-scale twodimensional graphene sheets and bulk three-dimensional lithium-ion electrolyte systems. In conclusion, employing 55 296 computational cores, the time per self-consistent field iteration for a sample of the bulk 3D electrolyte containing 8586 atoms is 90 s, and the time for a graphene sheet containing 11 520 atoms is 75 s.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Experimental evaluation of localized creep deformation in grade 91 steel weldments

Spatially resolved measurement of localized creep deformation in heterogeneous creep resistant steel weldments is crucial but challenging for lifetime assessments of critical steam components in power plants. Here, experimental approaches were established to quantitatively evaluate commonly observed localized creep deformation in multi-pass Grade 91 steel weldments. An in-situ digital image correlation (DIC) system was utilized with a creep testing frame to monitor and measure both full-field strain and localized strain accumulation across the weldments during long-term creep testing at elevated temperatures. The in situ DIC method measured not only the creep deformation behavior of the weld metal, heat affected zone (HAZ), and base metal, but also creep strain evolution for each sub-region within the HAZ itself, including coarse-grained HAZ, fine-grained HAZ (FGHAZ), and intercritical HAZ (ICHAZ). The DIC results revealed that local creep strain in the ICHAZ reached up to 90% strain before final rupture, whereas nominal creep strain measured by the standard extensometer of the tested cross-weld specimen was below 10%, indicative of Type IV cracking of the Grade 91 weld. Microstructural analyses revealed that the faster creep degradation/deformation in the HAZ was caused mainly by accelerated matrix grain recrystallization/growth and a reduced pinning effect from the segregated and coarsened precipitates in the FGHAZ and ICHAZ. The ultimate creep rupture occurred in the ICHAZ owing to its lowest creep resistance induced by the largest recrystallized grain size, the lowest fraction of coincidence site lattice, and the lowest local strain energies/dislocation densities.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Phase-field based Multiscale Modeling of Heterogeneous Solid Electrolytes: Applications to Nanoporous Li 3 PS 4

Modeling the effective ion conductivities of heterogeneous solid electrolytes typically involves the use of a computer-generated microstructure consisting of randomly or uniformly oriented fillers in a matrix. But, the structural features of the filler/matrix interface, which critically determine the interface ion conductivity and the microstructure morphology, have not been considered during the microstructure generation. In using nanoporous β-Li 3 PS 4 electrolyte as an example, we develop a phase-field model that enables generating nanoporous microstructures of different porosities and connectivity patterns based on the depth and the energy of the surface (pore/electrolyte interface), both of which are predicted through density functional theory (DFT) calculations. Room-temperature effective ion conductivities of the generated microstructures are then calculated numerically, using DFT-estimated surface Li-ion conductivity (3.14×10 -3 S/cm) and experimentally measured bulk Li-ion conductivity (8.93×10 -7 S/cm) of β-Li 3 PS 4 as the inputs. We also use the generated microstructures to inform effective medium theories to rapidly predict the effective ion conductivity via analytical calculations. Furthemore, when porosity approaches the percolation threshold, both the numerical and analytical methods predict a significantly enhanced Li-ion conductivity (1.74×10 -4 S/cm) that is in good agreement with experimental data (1.64×10 -4 S/cm). The present phase-field based multiscale model is generally applicable to predict both the microstructure patterns and the effective properties of heterogeneous solid electrolytes.

36 MATERIALS SCIENCE↗
NASA NTRS2014

Development and Characterization of Embedded Sensory Particles Using Multi-Scale 3D Digital Image Correlation

A method for detecting fatigue cracks has been explored at NASA Langley Research Center. Microscopic NiTi shape memory alloy (sensory) particles were embedded in a 7050 aluminum alloy matrix to detect the presence of fatigue cracks. Cracks exhibit an elevated stress field near their tip inducing a martensitic phase transformation in nearby sensory particles. Detectable levels of acoustic energy are emitted upon particle phase transformation such that the existence and location of fatigue cracks can be detected. To test this concept, a fatigue crack was grown in a mode-I single-edge notch fatigue crack growth specimen containing sensory particles. As the crack approached the sensory particles, measurements of particle strain, matrix-particle debonding, and phase transformation behavior of the sensory particles were performed. Full-field deformation measurements were performed using a novel multi-scale optical 3D digital image correlation (DIC) system. This information will be used in a finite element-based study to determine optimal sensory material behavior and density.

Cornell, Stephen R.↗
DOE OSTI2021

Thermal History of Matrix Forsterite Grains from Murchison Based on High-resolution Tomography

Protoplanetary disks are dust- and gas-rich structures surrounding protostars. Depending on the distance from the protostar, this dust is thermally processed to different degrees and accreted to form bodies of varying chemical compositions. The primordial accretion processes occurring in the early protoplanetary disk such as chondrule formation and metal segregation are not well understood. One way to constrain them is to study the morphology and composition of forsteritic grains from the matrix of carbonaceous chondrites. Here, we present high-resolution ptychographic X-ray nanotomography and multimodal chemical microtomography (X-ray diffraction and X-ray fluorescence) to reveal the early history of forsteritic grains extracted from the matrix of the Murchison CM2.5 chondrite. Here,. the 3D electron density maps revealed, at unprecedented resolution (64 nm), spherical inclusions containing Fe–Ni, very little silica-rich glass and void caps (i.e., volumes where the electron density is consistent with conditions close to vacuum) trapped in forsterite. The presence of the voids along with the overall composition, petrological textures, and shrinkage calculations is consistent with the grains experiencing one or more heating events with peak temperatures close to the melting point of forsterite (~2100 K), and subsequently cooled and contracted, in agreement with chondrule-forming conditions.

58 GEOSCIENCES↗
DOE OSTI2019

Scalable electric-field-assisted fabrication of vertically aligned carbon nanotube membranes with flow enhancement

In this work, we report a solution-based approach for fabricating vertically aligned carbon-nanotube (VACNT) membranes, which improves on the unaligned nature of CNT mixed-matrix membranes. This approach addresses important challenges regarding scalability and incorporation of different types and sizes of nanotube as pores in a membrane. In a new, solvent-deposition approach, CNTs are dispersed and then aligned and concentrated via electro-deposition from an organic solvent, using a composite AC + DC electric field to form aligned CNT arrays. This enables VACNT number densities that are two orders-of-magnitude higher than possible from direct suspension of CNTs in the liquid oligomer. The solvent is then replaced by a UV-curable oligomer, which is selectively cured to form macroscopic membranes of highly controlled thickness. After plasma etching to open pores, the VACNT membranes are shown to be permeable, and to have pore sizes consistent with flow through internal channels of the CNTs. In a first for field-aligned, solution-fabricated VACNT membranes, the pores show substantial, 100–300 × , gas-flow enhancement compared to theory, similar to membranes fabricated via less scalable methods (i.e., by CVD growth of already aligned CNT forests). The solution-based, electric-field-assisted approach aligns nanotubes of different sizes and types grown by any means, and is scalable to the efficient fabrication of large-area VACNT membranes for a variety of important applications.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Theory of x-ray scattering from laser-driven electronic systems

We describe, within the framework of quantum electrodynamics, an interaction between a nonresonant hard x-ray pulse and an electronic system in the presence of a temporally periodic laser field driving electron dynamics in this system. Here we apply Floquet theory to describe the laser-driven electronic system, and then obtain the scattering probability of an arbitrary nonresonant x-ray pulse from such a system employing the density-matrix formalism. We show that the scattering probability can be connected to the time-dependent electron density of the driven electronic system only under certain conditions, in particular, if the bandwidth of the probe x-ray pulse is sufficiently narrow to spectroscopically resolve transitions to different final states. A special focus is laid on application of the theory to laser-driven crystals in a strongly nonperturbative regime. We show how the time-dependent electron density of a crystal can be reconstructed from energy-resolved scattering patterns. This is illustrated by a calculation of a diffraction signal from a driven MgO crystal.

36 MATERIALS SCIENCE↗
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↗