Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “matrix completion”

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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Putting AlphaFold models to work with phenix.process_predicted_model and ISOLDE

AlphaFold has recently become an important tool in providing models for experimental structure determination by X-ray crystallography and cryo-EM. Large parts of the predicted models typically approach the accuracy of experimentally determined structures, although there are frequently local errors and errors in the relative orientations of domains. Importantly, residues in the model of a protein predicted by AlphaFold are tagged with a predicted local distance difference test score, informing users about which regions of the structure are predicted with less confidence. AlphaFold also produces a predicted aligned error matrix indicating its confidence in the relative positions of each pair of residues in the predicted model. The phenix.process_predicted_model tool downweights or removes low-confidence residues and can break a model into confidently predicted domains in preparation for molecular replacement or cryo-EM docking. These confidence metrics are further used in ISOLDE to weight torsion and atom–atom distance restraints, allowing the complete AlphaFold model to be interactively rearranged to match the docked fragments and reducing the need for the rebuilding of connecting regions.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Hydro-mechanical measurements and x-ray imaging of sheared crystalline rock fractures from EGS Collab experiments 1 & 2

Measuring hydro-mechanical properties of natural fractures is a prerequisite for optimizing hydraulic stimulation design and well placement. We completed experiments to characterize shear on natural fractures in schist, amphibolite, and rhyolite specimens drilled from EGS Collab Project’s field sites at the Sanford Underground Research Facility (SURF) in South Dakota. A triaxial direct shear method and coupled x-ray imaging were used to perform hydroshearing and mechanical shearing at the site’s in-situ stress conditions. This produced simultaneous measurements of fracture and matrix strength, permeability, stress-dependent aperture, dilation, and friction strength. Our results identified that only a subset of the natural fractures was weak enough for hydroshearing. Generally, hydroshearing increases fracture permeability by a factor of 10 or more and the enhancement is retainable over time. However, the shear slip does not always result in permeability enhancement. High content of phyllosilicates was found to associate with exceptionally weak fractures that also exhibited poor or even negative enhancement after stimulation. Combining our measurements with site data, we can predict that most observable fractures at the two EGS Collab sites do not meet the criteria for hydroshearing before tensile opening. In some cases, the visible fractures are low permeability and as strong as the adjacent rock. To induce hydroshearing before tensile opening, injection must target known weak and favorably oriented fractures with confirmed pre-existing permeability.

58 GEOSCIENCES↗
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 OSTI2023

Nucleation Selectivity and Lateral Coalescence of GaAs over Graphene on Ge(111)

Here we use epitaxial lateral overgrowth (ELO) to produce semimetallic graphene nanostructures embedded in a semiconducting GaAs matrix for potential applications in plasmonics, THz generation and detection, and tunnel junctions in multijunction solar cells. We show that (1) the combination of low sticking coefficient and fast surface diffusion on graphene enhances nucleation selectivity at exposed regions of the substrate and (2) high growth temperatures favor efficient lateral overgrowth, coalescence, and planarization of epitaxial GaAs films over the graphene nanostructures. Our work provides a more complete understanding of ELO using graphene masks, as opposed to more conventional dielectric masks, and enables new types of metal/semiconductor nanocomposites.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Total metals & anion concentration data; Slate River floodplain, Crested Butte, CO; May 2020-September 2020

This data package includes processed and undiluted measurements for metal and anion concentrations from pore water (groundwater) samples from the Slate River floodplain of Crested Butte, CO, a focus field site for the SLAC Floodplain Hydro-Biogeochemistry SFA. The data was generated as part of the work targeting the overarching research question for the SLAC SFA: How do ubiquitous subsurface interfaces mediate molecular-scale biogeochemical processes and groundwater quality in floodplains and watersheds? Samples were collected between May and September of 2020. These measurements were all recorded at the Arizona Laboratory for Emerging Contaminants (ALEC) at the University of Arizona located in Tucson, AZ. Groundwater samples were extracted from a network of installed rhizon (Rhizosphere Research Products, part no. 19.60.21F, 0.6 micrometer mesh size) and piezometer wells within the river floodplain. All water samples were shaded from sun exposure during extraction from the subsurface and preserved at 4C until measured at ALEC.Analysis by ICP-MS:Measurements for total metals were made on the Agilent 7700x ICP-MS (for total metals) – Agilent Technologies, Santa Clara, CA.The analytical QA/QC protocol was adapted from US EPA Method 200.8 for analysis by ICP-MS. Calibration standards were prepared from multi-element stock solution (Sigma-Aldrich Multielement standard solution for ICP, St. Louis, MO) using matrix matched to sample solutions (either 2% HCl or HNO3 from AriStar Plus,grade acids from VWR Scientific). Calibration curves include at least 7 points with correlation coefficients > 0.995. The QC protocol includes a continuing calibration blank (CCB), a continuing calibration verification (CCV) solution and at least one quality control sample (QCS) to be analyzed just after calibration and again after every 12 samples and at the completion of the run. The QCS solutions are from an independent source, such as NIST SRM 1643e - Trace Elements in Water, or QCS solutions from High Purity Standards (Charleston, SC). Acceptable QC responses must be between 90 and 110% of the certified value. An internal standard (Rh) is added via on-line addition into the sample line using a mixing tee.Analysis by Ion Chromatography (Anions):The protocol follows Method 4110 in Standard Methods for Examination of Water and Wastewater.The instrument used is the Thermo Scientific Dionex ICS-6000 using AS+AG22 column set for anion analysis with isocratic method using sodium carbonate eluent. Detection is by chemical suppression of eluent conductivity. Quality control solutions and mixed analyte standards purchased from Inorganic Ventures, Christiansburg, VA.All files are in csv format.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

LLNL measurements of thermally irradiated HEU sample and saltwater-matrixed HEU sample

Two scintillation vials, kindly provided by the PNNL team, were received in the LLNL radiochemistry building on June 10, 2024; one containing approximately 2 mL dissolved thermally-irradiated HEU, and the other containing dried salt from irradiated seawater. A quick, semiquantitative screening count of the salt showed that the activation product activities were very low, with the highest activity being ~1 Bq 24 Na. The dissolved HEU solution was transferred to a Prindle vial (LLNL standard counting geometry) and weighed. Approximately 2 mL 3 M HNO 3 was added so that the solution completely covered the bottom of the vial and Gamanal software would be able to accurately generate an efficiency curve. This sample, called “PNNLFP24” is the irradiated HEU solution received, gravimetrically diluted by a factor of 1.7356 ± 0.0003, and was sent for quick gamma counting prior to further modification. Meanwhile the salt was quantitatively transferred to a 250 mL polyethylene bottle and dissolved in 160 mL 3 M HNO 3 . After counting, 50 µL PNNLFP24 was aliquoted for Resonance Ionization Mass Spectrometry (RIMS) analysis, and then 10.0 mL of the salt solution (6.03% of the total salt, gravimetrically) was added to the remaining PNNLFP24 solution to make “PNNLFPSW24”. After weighing, two 100 µL aliquots were removed for RIMS and DES, and 3 mL was aliquoted for microfluidic chemistry. The remaining “PNNLFPSW24” solution was weighed and proceeded for singles and coincidence counting on the MCBOS system.

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

Evaluation of two-particle properties within finite-temperature self-consistent one-particle Green’s function methods: Theory and application to GW and GF2

One-particle Green’s function methods can model molecular and solid spectra at zero or non-zero temperatures. One-particle Green’s functions directly provide electronic energies and one-particle properties, such as dipole moment. However, the evaluation of two-particle properties, such as $\langle$S 2 $\rangle$ and $\langle$N 2 $\rangle$, can be challenging because they require a solution of the computationally expensive Bethe–Salpeter equation to find two-particle Green’s functions. We demonstrate that the solution of the Bethe–Salpeter equation can be completely avoided. Applying the thermodynamic Hellmann–Feynman theorem to self-consistent one-particle Green’s function methods, we derive expressions for two-particle density matrices in a general case and provide explicit expressions for GF2 and GW methods. Such density matrices can be decomposed into an antisymmetrized product of correlated one-electron density matrices and the two-particle electronic cumulant of the density matrix. Cumulant expressions reveal a deviation from ensemble representability for GW, explaining its known deficiencies. We analyze the temperature dependence of $\langle$S 2 $\rangle$ and $\langle$N 2 $\rangle$ for a set of small closed-shell systems. Interestingly, both GF2 and GW show a non-zero spin contamination and a non-zero fluctuation of the number of particles for closed-shell systems at the zero-temperature limit.

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

SQMBox: Interfacing a semiempirical integral library to modular ab initio electronic structure enables new semiempirical methods

Ab initio and semiempirical electronic structure methods are usually implemented in separate software packages or use entirely different code paths. As a result, it can be time-consuming to transfer an established ab initio electronic structure scheme to a semiempirical Hamiltonian. Here we present an approach to unify ab initio and semiempirical electronic structure code paths based on a separation of the wavefunction ansatz and the needed matrix representations of operators. With this separation, the Hamiltonian can refer to either an ab initio or semiempirical treatment of the resulting integrals. We built a semiempirical integral library and interfaced it to the GPU-accelerated electronic structure code TeraChem. Equivalency between ab initio and semiempirical tight-binding Hamiltonian terms is assigned according to their dependence on the one-electron density matrix. The new library provides semiempirical equivalents of the Hamiltonian matrix and gradient intermediates, corresponding to those provided by the ab initio integral library. This enables the straightforward combination of semiempirical Hamiltonians with the full pre-existing ground and excited state functionality of the ab initio electronic structure code. We demonstrate the capability of this approach by combining the extended tight-binding method GFN1-xTB with both spin-restricted ensemble-referenced Kohn–Sham and complete active space methods. We also present a highly efficient GPU implementation of the semiempirical Mulliken-approximated Fock exchange. The additional computational cost for this term becomes negligible even on consumer-grade GPUs, enabling Mulliken-approximated exchange in tight-binding methods for essentially no additional cost.

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

Physi-sorption of H 2 on pure and boron-doped graphene monolayers: A dispersion-corrected DFT study

Carbon based materials are of interest as potential candidates for H 2 storage. Earlier work has been inconclusive on the effect of boron doping on the energy of H 2 binding. However, earlier work has been inconclusive on the definitive effect of boron doping on the energy of H 2 binding, i.e. isosteric heats of adsorption (Qst). In this work, we completed a systematic DFT study to evaluate this effect, and found that doping graphene with boron provides only minor enhancement in H 2 binding. More importantly, the presence of the electron deficient boron into a graphene ring introduces a defect, such as terminal hydrogen or distortion from planarity, which creates hydrogen adsorption sites with slightly increased Qst. The increase is from ~ 5 kJ/mol H 2 for the pure carbon matrix to ~6 - 7 kJ/mol for the boron doped system. The more strongly bond H 2 is located near the defect and shows little direct interaction with the boron. Most significant enhancement is found in systems where H 2 is confined between layers at a distance of about 7Å. In this case, the H 2 binding nearly doubles, to ~10 kJ/mol for both pure, undistorted graphene, and the 2% boron doped system, containing tetrahedral C atom distortion. Interestingly, at higher doping levels of boron, we found that the Qst decreases compared to the all carbon materials. These finding suggest that interplanar nanoconfinement may be more effective way to enhance H 2 binding than boron doping.

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

Fuel performance analysis of fully-resolved TRISO compact

The TRi-structural ISOtropic (TRISO) fuel multilayered coating structure offers multiple barriers to fission product release, enhancing safety and performance. The heterogeneous nature of TRISO fuel compacts, comprising thousands of randomly distributed coated fuel particles embedded in a graphite matrix, creates intricate stress fields and thermal gradients that cannot be accurately modeled using simplified one-dimensional or homogenized approaches. Consequently, three-dimensional modeling enables the prediction of fuel compact dimensional changes, internal pressure buildup, and fission product transport pathways under diverse irradiation and thermal conditions. This capability facilitates detailed analysis of particle-to-particle interactions, matrix cracking mechanisms, and the statistical distribution of coating failures, which directly impact fuel performance and safety margins. This capability is particularly critical for advanced reactors, such as high-temperature gas-cooled reactors and other Generation IV reactor designs where TRISO fuel operates at elevated temperatures and burn-up levels. This work introduces a novel method to generate an optimized packing of TRISO compacts and a complete 3D mesh with random distribution of TRISO particles, which are discretized into each coating component layer.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Nonperturbative random matrix model of N = 2 JT supergravity

It is shown how to nonperturbatively define a random matrix model that captures key physics of N = 2 Jackiw-Teitelboim supergravity, going well beyond the perturbative topological expansion defined recently by Turiaci and Witten. A decomposition into an infinite family of certain multicritical models is derived, leading to the definition of a nonlinear ordinary differential equation from which the physics may be computed. Bogomol’nyi-Prasad-Sommerfield (BPS) states are naturally described by the model. The nonperturbative completions of the spectral densities for non-BPS multiplets are readily extracted. Published by the American Physical Society 2024

Johnson, Clifford V. (ORCID:0000000189645830)↗
DOE OSTI2022

Radiation induced amorphization of carbides in additively manufactured and conventional ferritic-martensitic steels: In-situ experiments on extraction replicas

In-situ irradiations using 1 MeV Kr 2+ ions in a transmission electron microscope were performed on extraction replica samples containing M 23 C 6 carbides and MX carbonitrides from two ferritic-martensitic (FM) steels: (i) 9%Cr-1%W-TaV based Eurofer97 and (ii) 9%Cr-1%Mo-VNb (all in wt.%) based additively manufactured (AM) Grade91 steel. The irradiations were performed between 100 and 773 K up to a maximum dose of ~2.4 displacements per atom (dpa). The M 23 C 6 carbides are highly susceptible to radiation induced amorphization (RIA), while the MX type nanoprecipitates are highly amorphization resistant across the entire irradiation temperature range. Between 100 and 423 K, RIA of M 23 C 6 carbides occurs very rapidly with critical amorphization doses ranging between ~0.35 and 0.9 dpa, increasing to higher doses at 573 K. Complete amorphization of the M 23 C 6 carbides up to doses of ~2.4 dpa is not possible at 773 K under the present irradiation conditions. The critical amorphization dose of M 23 C 6 carbides increases nearly exponentially with irradiation temperature. The critical temperature for the crystalline-to-amorphous phase transformation (T c→a ) of M 23 C 6 carbides irradiated as replica samples, i.e. without the surrounding metal matrix, was estimated to be ~812 K. Comparing the present results with neutron irradiation data on bulk samples reveals a decrease in the critical amorphization dose for M 23 C 6 and an increase of T c→a , highlighting the effect of dose rate on the amorphization behavior that is qualitatively consistent with literature on other non-metals. Changes in the minor chemistry of M 23 C 6 , such as presence or absence of W, V, Nb, Mo, seems to have little effect on the amorphization behavior. Here, in-situ irradiations on extraction replica samples provide a novel pathway to explore radiation tolerance of nanoprecipitates in nuclear structural materials

36 MATERIALS SCIENCE↗
DOE OSTI2020

Superallowed 0 + → 0 + nuclear β decays: 2020 critical survey, with implications for V ud and CKM unitarity

A new critical survey of all half-life, decay-energy and branching-ratio measurements related to 23 superallowed 0 + → 0 + β decays is presented. Included are 222 individual measurements of comparable precision obtained from 174 published references. Compared with our last survey in 2015, we have added results from 28 new publications and eliminated an approximately equal number whose results have been superseded by much more precise modern data. We obtain world-average ft values for each of the 21 transitions that have a complete set of data, then apply radiative and isospin-symmetry-breaking corrections to extract “corrected” Ft values. Fifteen of these Ft values now have a precision of 0.3% or better and all take the same value within statistics, as expected from conservation of the vector current. Their average, Ft¯, when combined with the muon lifetime, yields the up-down quark-mixing element of the Cabibbo-Kobayashi-Maskawa matrix, V ud = 0.97373 ± 0.00031. This is lower than our 2015 result by one standard deviation and its uncertainty is increased by 50%. This is a consequence, not of any shifts in the experimental data, but of new calculations for the radiative corrections. The lower V ud value now leads to greater tension in the top-row test of unitarity in the CKM matrix. Updates in experimental data have independently led to a factor-of-two tighter limit being set on the possible existence of a scalar interaction. In conclusion, the new limit on Fierz interference is b F ≤ 0.0033 at the 90% confidence level.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

A greedy algorithm for computing eigenvalues of a symmetric matrix with localized eigenvectors

Here, we present a greedy algorithm for computing selected eigenpairs of a large sparse matrix $H$ that can exploit localization features of the eigenvector. When the eigenvector to be computed is localized, meaning only a small number of its components have large magnitudes, the proposed algorithm identifies the location of these components in a greedy manner, and obtains approximations to the desired eigenpairs of $H$ by computing eigenpairs of a submatrix extracted from the corresponding rows and columns of $H$. Even when the eigenvector is not completely localized, the approximate eigenvectors obtained by the greedy algorithm can be used as good starting guesses to accelerate the convergence of an iterative eigensolver applied to $H$. We discuss a few possibilities for selecting important rows and columns of $H$ and techniques for constructing good initial guesses for an iterative eigensolver using the approximate eigenvectors returned from the greedy algorithm. We demonstrate the effectiveness of this approach with examples from nuclear quantum many-body calculations and many-body localization studies of quantum spin chains.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Solving quasi-free and quadratic Lindblad master equations for open fermionic and bosonic systems

The dynamics of Markovian open quantum systems are described by Lindblad master equations. For fermionic and bosonic systems that are quasi-free, i.e. with Hamiltonians that are quadratic in the ladder operators and Lindblad operators that are linear in the ladder operators, we derive the equation of motion for the covariance matrix. This determines the evolution of Gaussian initial states and the steady states, which are also Gaussian. Using ladder super-operators (a.k.a. third quantization), here we show how the Liouvillian can be transformed to a many-body Jordan normal form which also reveals the full many-body spectrum. Extending previous work by Prosen and Seligman, we treat fermionic and bosonic systems on equal footing with Majorana operators, shorten and complete some derivations, also address the odd-parity sector for fermions, give a criterion for the existence of bosonic steady states, cover non-diagonalizable Liouvillians also for bosons, and include quadratic systems. In extension of the quasi-free open systems, quadratic open systems comprise additional Hermitian Lindblad operators that are quadratic in the ladder operators. While Gaussian states may then evolve into non-Gaussian states, the Liouvillian can still be transformed to a useful block-triangular form, and the equations of motion for k-point Green’s functions form a closed hierarchy. Based on this formalism, results on criticality and dissipative phase transitions in such models are discussed in a companion paper.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Linear and Nonlinear Solvers for Simulating Multiphase Flow within Large-Scale Engineered Subsurface Systems

Simulation of multiphase flow in the subsurface is well-known to be computationally challenging. While there have been many studies that have explored approaches to overcoming these challenges, they often utilize relatively simple case studies. In this paper, we focus on the unique numerical challenges posed by modeling large-scale engineered subsurface systems, characterized by discrete features embedded in a heterogeneous natural subsurface setting. The man-made features such as shafts, tunnels, and barriers often cause multiple challenges in modeling the domain for multiphase porous media flow. This flow scenario can have a wide range of applications such as nuclear waste repositories, enhanced recovery of a petroleum reservoir, geothermal engineering, and carbon sequestration. An example of these severe numerical challenges is the case of performance assessment (PA) for Waste Isolation Pilot Plant (WIPP), the only operating deep geological repository in the US, which simulates extreme material properties of bedded salt rock formation and extreme contrast due to open excavation next to the formation. The models have extremes not only of permeability and porosity but also of the constitutive models needed for multiphase flow; additionally, they have process models like salt creep closure reducing porosity over time, fracturing in clay and anhydrite interbeds of the bedded salt, gas generation from the waste materials, and unintentional human borehole intrusions in some scenarios. Numerical simulations require the solution of coupled systems of nonlinear PDEs; in our work, we use the open-source simulator PFLOTRAN which is based on Finite Volume discretization. The solution of the nonlinear equations requires use of the Newton-Raphson iteration at each time step, which entails the solution of the linearized Jacobian system at each iteration. The effects of all the processes (i.e., large number of unknowns, highly nonlinear constitutive relations, large contrasts in material properties in short distances) lead to an ill-conditioned Jacobian matrix that severely challenges traditional linear solver, i.e., stabilized biconjugate gradient with block Jacobi incomplete LU preconditioner (BCGS-ILU) leading to non-convergence for traditional Newton-Raphson nonlinear solver causing unacceptably long computation time for each model. This paper presents linear solvers such as constrained pressure residual (CPR) two-stage preconditioner with alternate-block-factorization (ABF) and quasi- implicit pressure and explicit saturation (QIMPES) decouplers and flexible generalized residual solver (FGMRES). The new general-purpose nonlinear solver, Newton trust-region dogleg Cauchy (NTRDC), is also introduced to resolve extreme nonlinearities in the models. We demonstrate the effectiveness of each method relative to the default BCGS-Newton solver. The two best cases had nearly 50 times speed-up and achieved completion of a simulation in 14 hours that never completed due to non-convergence with the default solver. We also investigate the strong scalability of each method and discuss some of the deficiencies found for Block Jacobi preconditioner using parallel domain decomposition, and node packing effects of modern processor architecture.

Preconditioner, Nonlinear, Porous media, Multiphas↗
DOE OSTI2022

Randomized Cholesky Preconditioning for Graph Partitioning Applications

Graph partitioning has emerged as an area of interest due to its use in various applications in computational research. One way to partition a graph is to solve for the eigenvectors of the corresponding graph Laplacian matrix. This project focuses on the eigensolver LOBPCG and the evaluation of a new preconditioner: Randomized Cholesky Factorization (rchol). This proconditioner was tested for its speed and accuracy against other well-known preconditioners for the method. After experiments were run on several known test matrices, rchol appears to be a better preconditioner for structured matrices. This research was sponsored by National Nuclear Security Administration Minority Serving Institutions Internship Program (NNSA-MSIIP) and completed at host facility Sandia National Laboratories. As such, after discussion of the research project itself, this report contains a brief reflection on experience gained as a result of participating in the NNSA-MSIIP.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Hybrid quantum-classical simulations of magic angle spinning dynamic nuclear polarization in very large spin systems

Solid-state nuclear magnetic resonance can be enhanced using unpaired electron spins with a method known as dynamic nuclear polarization (DNP). Fundamentally, DNP involves ensembles of thousands of spins, a scale that is difficult to match computationally. This scale prevents us from gaining a complete understanding of the spin dynamics and applying simulations to design sample formulations. We recently developed an ab initio model capable of calculating DNP enhancements in systems of up to ~1000 nuclei; however, this scale is insufficient to accurately simulate the dependence of DNP enhancements on radical concentration or magic angle spinning (MAS) frequency. We build on this work by using ab initio simulations to train a hybrid model that makes use of a rate matrix to treat nuclear spin diffusion. We show that this model can reproduce the MAS rate and concentration dependence of DNP enhancements and build-up time constants. Finally, we then apply it to predict the DNP enhancements in core–shell metal-organic-framework nanoparticles and reveal new insights into the composition of the particles’ shells.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗