Search NASA⌕ Search

SEARCH · Search NASA

Results for “real-space 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.

126 records · Page 7

Exotic phases in finite-density ℤ 3 theories

Lattice ℤ 3 theories with complex actions share many key features with finite- density QCD including a sign problem and $\mathcal{CK}$ symmetry. Complex ℤ 3 spin and gauge models exhibit a generalized Kramers-Wannier duality mapping them onto chiral ℤ 3 spin and gauge models, which are simulatable with standard lattice methods in large regions of parameter space. The Migdal-Kadanoff real-space renormalization group (RG) preserves this duality, and we use it to compute the approximate phase diagram of both spin and gauge ℤ 3 models in dimensions one through four. Chiral ℤ 3 spin models are known to exhibit a Devil’s Flower phase structure, with inhomogeneous phases that can be thought of as ℤ 3 analogues of chiral spirals. Out of the large class of models we study, we find that only chiral spin models and their duals have a Devil’s Flower structure with an infinite set of inhomogeneous phases, a result we attribute to Elitzur’s theorem. We also find that different forms of the Migdal-Kadanoff RG produce different numbers of phases, a violation of the expectation for universal behavior from a real-space RG. We discuss extensions of our work to ℤ N models, SU(N) models and nonzero temperature.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Framework of compressive sensing and data compression for 4D-STEM

Four-dimensional Scanning Transmission Electron Microscopy (4D-STEM) is a powerful technique for high-resolution and high-precision materials characterization at multiple length scales, including the characterization of beam-sensitive materials. However, the field of view of 4D-STEM is relatively small, which in absence of live processing is limited by the data size required for storage. Furthermore, the rectilinear scan approach currently employed in 4D-STEM places a resolution- and signal-dependent dose limit for the study of beam sensitive materials. Improving 4D-STEM data and dose efficiency, by keeping the data size manageable while limiting the amount of electron dose, is thus critical for broader applications. Here we introduce a general method for reconstructing 4D-STEM data with subsampling in both real and reciprocal spaces at high fidelity. The approach is first tested on the subsampled datasets created from a full 4D-STEM dataset, and then demonstrated experimentally using random scan in real-space. The same reconstruction algorithm can also be used for compression of 4D-STEM datasets, leading to a large reduction (100 times or more) in data size, while retaining the fine features of 4D-STEM imaging, for crystalline samples.

4D-STEM↗

Earth Mover’s Distance as a Metric to Evaluate the Extent of Charge Transfer in Excitations Using Discretized Real-Space Densities

This paper presents a novel theoretical measure, μ EMD , based on the earth mover's distance (EMD), for quantifying the density shift caused by electronic excitations in molecules. As input, the EMD metric uses only the discretized ground- and excited-state electron densities in real space, rendering it compatible with almost all electronic structure methods used to calculate excited states. The EMD metric is compared against other popular theoretical metrics for describing the extent of electron-hole separation in a wide range of excited states (valence, Rydberg, charge transfer, etc.). Further, the results showcase the EMD metric's effectiveness across all excitation types and suggest that it is useful as an additional tool to characterize electronic excitations. The study also reveals that μ EMD can function as a promising diagnostic tool for predicting the failure of pure exchange-correlation functionals. Specifically, we show statistical relationships among the functional-driven errors, the exact exchange content within the functional, and the magnitude of μ EMD values.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Simulation of 24,000 Electron Dynamics: Real-Time Time-Dependent Density Functional Theory (TDDFT) with the Real-Space Multigrids (RMG)

Here, we present the theory, implementation, and benchmarking of a real-time time-dependent density functional theory (RT-TDDFT) module within the RMG code, designed to simulate the electronic response of molecular systems to external perturbations. Our method offers insights into nonequilibrium dynamics and excited states across a diverse range of systems, from small organic molecules to large metallic nanoparticles. Benchmarking results demonstrate excellent agreement with established TDDFT implementations and showcase the superior stability of our time integration algorithm, enabling long-term simulations with minimal energy drift. The scalability and efficiency of RMG on massively parallel architectures allow for simulations of complex systems, such as plasmonic nanoparticles with thousands of atoms. Future extensions, including nuclear and spin dynamics, will broaden the applicability of this RT-TDDFT implementation, providing a powerful toolset for studies of photoactive materials, nanoscale devices, and other systems where real-time electronic dynamics is essential.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Mechanisms of Skyrmion and Skyrmion Crystal Formation from the Conical Phase

Real-space topological magnetic structures such as skyrmions and merons are promising candidates for information storage and transport. However, the microscopic mechanisms that control their formation and evolution are still unclear. In this work, using in situ Lorentz transmission electron microscopy, we demonstrate that skyrmion crystals (SkXs) can nucleate, grow, and evolve from the conical phase in the same ways that real nanocrystals form from vapors or solutions. More intriguingly, individual skyrmions can also “reproduce” by division in a mitosis-like process that allows them to annihilate SkX lattice imperfections, which is not available to crystals made of mass-conserving particles. Combined string method and micromagnetic calculations show that competition between repulsive and attractive interactions between skyrmions governs particle-like SkX growth, but nonconservative SkX growth appears to be defect mediated. Our results provide insights toward manipulating magnetic topological states by applying established crystal growth theory, adapted to account for the new process of skyrmion mitosis.

36 MATERIALS SCIENCE↗

Bridging the gap between molecules and materials in quantum chemistry with localized active spaces

The number of materials that “bridge the gap” between single molecules and extended solids, such as metal-organic frameworks and organic semiconductors, has been increasing. Consequently, there is a growing need for modeling approaches that effectively integrate the real-space molecular perspective employed by computational chemists and the reciprocal-space dispersive perspective employed by computational physicists. Here, we propose the localized active space (LAS) approach as a promising method to successfully bridge this gap. The LAS approach extends the active space concept from multiconfigurational methods such as complete active space self-consistent field theory to multiple molecular fragments via a product-form wave function ansatz. Here, we apply this method to solid state phenomena by treating each unit cell as a fragment with different sets of local quantum numbers (e.g., charge and excitation number). State interaction between these LAS states (LASSI) thus provides a comprehensive basis for the study of charge and energy transfer, meeting and surpassing the capabilities of single-reference fragmentation approaches such as constrained density functional theory (cDFT). Most centrally, we show how combining this LASSI approach with multiconfigurational pair-density functional theory (MC-PDFT) provides an elegant and efficient method to compute band structures that capture multiconfigurational character. We apply the LASSI band structure approach to the computation of band gaps in stretched hydrogen chain, polyacetylene, and bulk nickel oxide (NiO), finding good or excellent quantitative agreement with reference values in all cases. Additionally, we use the LAS basis in one-dimensional model systems to demonstrate its ability to treat difficult solid-state phenomena such as exciton transfer and excitation at p-n junctions.

method development↗

Data-driven discovery of dynamics from time-resolved coherent scattering

Coherent X-ray scattering (CXS) techniques are capable of interrogating dynamics of nano- to mesoscale materials systems at time scales spanning several orders of magnitude. However, obtaining accurate theoretical descriptions of complex dynamics is often limited by one or more factors—the ability to visualize dynamics in real space, computational cost of high-fidelity simulations, and effectiveness of approximate or phenomenological models. In this work, we develop a data-driven framework to uncover mechanistic models of dynamics directly from time-resolved CXS measurements without solving the phase reconstruction problem for the entire time series of diffraction patterns. Our approach uses neural differential equations to parameterize unknown real-space dynamics and implements a computational scattering forward model to relate real-space predictions to reciprocal-space observations. This method is shown to recover the dynamics of several computational model systems under various simulated conditions of measurement resolution and noise. Moreover, the trained model enables estimation of long-term dynamics well beyond the maximum observation time, which can be used to inform and refine experimental parameters in practice. Finally, we demonstrate an experimental proof-of-concept by applying our framework to recover the probe trajectory from a ptychographic scan. Our proposed framework bridges the wide existing gap between approximate models and complex data.

36 MATERIALS SCIENCE↗

QMCPACK: Advances in the development, efficiency, and application of auxiliary field and real-space variational and diffusion quantum Monte Carlo

We review recent advances in the capabilities of the open source ab initio Quantum Monte Carlo (QMC) package QMCPACK and the workflow tool Nexus used for greater efficiency and reproducibility. The auxiliary field QMC (AFQMC) implementation has been greatly expanded to include k-point symmetries, tensor-hypercontraction, and accelerated graphical processing unit (GPU) support. These scaling and memory reductions greatly increase the number of orbitals that can practically be included in AFQMC calculations, increasing the accuracy. Advances in real space methods include techniques for accurate computation of bandgaps and for systematically improving the nodal surface of ground state wavefunctions. Results of these calculations can be used to validate application of more approximate electronic structure methods, including GW and density functional based techniques. To provide an improved foundation for these calculations, we utilize a new set of correlation-consistent effective core potentials (pseudopotentials) that are more accurate than previous sets; these can also be applied in quantum-chemical and other many-body applications, not only QMC. These advances increase the efficiency, accuracy, and range of properties that can be studied in both molecules and materials with QMC and QMCPACK.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A super-resolution technique to analyze single-crystal inelastic neutron scattering measurements using direct-geometry chopper spectrometers

Direct-geometry time-of-flight chopper neutron spectroscopy is instrumental in studying dynamics in liquid, powder, and single crystal systems. We report here that real-space techniques in optical imagery can be adapted to obtain reciprocal-space super resolution dispersion for phonon or magnetic excitations from single-crystal neutron spectroscopy measurements. The procedure to reconstruct super-resolution energy dispersion of excitations relies on an accurate determination of the momentum and energy-dependent point spread function and a dispersion correction technique inspired by an image disparity calculation technique commonly used in stereo imaging. Here by, applying these methods to spinwave dispersion data from a virtual neutron experiment demonstrates ~5-fold improvement over nominal energy resolution.

47 OTHER INSTRUMENTATION↗

Dark energy survey year 3 results: cosmological constraints from the analysis of cosmic shear in harmonic space

ABSTRACT We present cosmological constraints from the analysis of angular power spectra of cosmic shear maps based on data from the first three years of observations by the Dark Energy Survey (DES Y3). Our measurements are based on the pseudo-Cℓ method and complement the analysis of the two-point correlation functions in real space, as the two estimators are known to compress and select Gaussian information in different ways, due to scale cuts. They may also be differently affected by systematic effects and theoretical uncertainties, making this analysis an important cross-check. Using the same fiducial Lambda cold dark matter model as in the DES Y3 real-space analysis, we find ${S_8 \equiv \sigma _8 \sqrt{\Omega _{\rm m}/0.3} = 0.793^{+0.038}_{-0.025}}$, which further improves to S8 = 0.784 ± 0.026 when including shear ratios. This result is within expected statistical fluctuations from the real-space constraint, and in agreement with DES Y3 analyses of non-Gaussian statistics, but favours a slightly higher value of S8, which reduces the tension with the Planck 2018 constraints from 2.3σ in the real space analysis to 1.5σ here. We explore less conservative intrinsic alignments models than the one adopted in our fiducial analysis, finding no clear preference for a more complex model. We also include small scales, using an increased Fourier mode cut-off up to $k_{\rm max}={5}\, {h}\, {\rm Mpc}^{-1}$, which allows to constrain baryonic feedback while leaving cosmological constraints essentially unchanged. Finally, we present an approximate reconstruction of the linear matter power spectrum at present time, found to be about 20 per cent lower than predicted by Planck 2018, as reflected by the lower S8 value.

79 ASTRONOMY AND ASTROPHYSICS↗

Boosting H I -Galaxy Cross-Clustering Signal through Higher-Order Cross-Correlations

After reionization, neutral hydrogen (${\rm H\, \small {I}}$) traces the large-scale structure (LSS) of the Universe, enabling ${\rm H\, \small {I}}$ intensity mapping (IM) to capture the LSS in 3D and constrain key cosmological parameters. We present a new framework utilizing higher-order cross-correlations to study ${\rm H\, \small {I}}$ clustering around galaxies, tested using real-space data from the IllustrisTNG300 simulation. This approach computes the joint distributions of k-nearest neighbor (kNN) optical galaxies and the ${\rm H\, \small {I}}$ brightness temperature field smoothed at relevant scales (the kNN-field framework), providing sensitivity to all higher-order cross-correlations, unlike two-point statistics. To simulate ${\rm H\, \small {I}}$ data from actual surveys, we add random thermal noise and apply a simple foreground cleaning model, filtering out Fourier modes of the brightness temperature field with k ∥ < k min,∥ . Under current levels of thermal noise and foreground cleaning, typical of a Canadian Hydrogen Intensity Mapping Experiment (CHIME)-like survey, the ${\rm H\, \small {I}}$-galaxy cross-correlation signal in our simulations, using the kNN-field framework, is detectable at >30σ across r = [3, 12] h –1 Mpc. In contrast, the detectability of the standard two-point correlation function (2PCF) over the same scales depends strongly on the foreground filter: a sharp k ∥ filter can spuriously boost detection to 8σ due to position-space ringing, whereas a less sharp filter yields no detection. Nonetheless, we conclude that kNN-field cross-correlations are robustly detectable across a broad range of foreground filtering and thermal noise conditions, suggesting their potential for enhanced constraining power over 2PCFs.

79 ASTRONOMY AND ASTROPHYSICS↗

Many-body electronic structure of LaScO 3 by real-space quantum Monte Carlo

We present real space quantum Monte Carlo (QMC) calculations of the scandate LaSc0 3 that proved to be challenging for traditional electronic structure approaches due to strong correlation effects resulting in inaccurate band gaps from DFT and GW methods when compared with existing experimental data. Besides calculating an accurate QMC band gap corrected for supercell size biases and in agreement with numerous experiments, we also predict the cohesive energy of the crystal using the standard fixed-node QMC without any empirical or non-variational parameters. We show that promotion (optical) gap and fundamental gap agree with each other illustrating a clear absence of significant excitonic effects in the ideal crystal. We obtained these results in perfect consistency in two independent tracks that employ different basis sets (plane wave vs. localized gaussians), different codes for generating orbitals (QUANTUM ESPRESSO vs. CRYSTAL), different QMC codes (QMCPACK vs. QWALK) and different high-accuracy pseudopotentials (ccECPs vs. Troullier-Martins) presenting the maturity and consistency of QMC methodology and tools for studies of strongly correlated problems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Real-space representation of the quasiparticle self-consistent GW self-energy and its application to defect calculations

The quasiparticle self-consistent (QS) GW (G for Green's function, W for screened Coulomb interaction) approach incorporates the corrections of the quasiparticle energies from their Kohn-Sham density functional theory (DFT) eigenvalues by means of an energy-independent and Hermitian self-energy matrix usually given in the basis set of the DFT eigenstates. By expanding these into an atom-centered basis set (specifically here the linearized muffin-tin orbitals) a real space representation of the self-energy corrections becomes possible. In this work, We show that this representation is relatively short-ranged. This offers opportunities to construct the self-energy of a complex system from parts of the system by a cut-and-paste method. Specifically for a point defect, represented in a large supercell, the self-energy can be constructed from those of the host and a smaller defect-containing cell. The self-energy of the periodic host can be constructed simply from a GW calculation for the primitive cell. We show for the case of the As Ga in GaAs that the defect part can already be well represented by a minimal eight-atom cell and allows us to construct the self-energy for a 64-atom cell in good agreement with direct QSGW calculations for the large cell. Using this approach to an even larger 216-atom cell shows the defect band approaches an isolated defect level. The calculations also allow us to identify a second defect band which appears as a resonance near the conduction band minimum. The results on the extracted defect levels agree well with Green's function calculations for an isolated defect and with experimental data.

36 MATERIALS SCIENCE↗

Dielectric screening and vacancy formation for large neutral and charged Si n H m ( n > 1500 ) nanocrystals using real-space pseudopotentials

A commonly used procedure for computing the properties of defects in crystalline materials is to consider a large supercell that includes the defect of interest. This is a straightforward technique as standard energy band codes can be used for such computations. For neutral defects, the only impediment of such an approach is to avoid defect-defect interactions between adjoining cells. However, this procedure can be complex if the defect of interest is charged as the system at large contains Coulombic divergences. Moreover, some have recently argued that the conventional definition of formation energies for charged defects cannot be reconciled with statistical mechanics. Here, we focus on an alternative approach. We consider large nanocrystals wherein a charged defect can be placed. Since the system is confined, a charged defect within the nanocrystal does not result in a Coulombic divergence. The chief impediment is computational, i.e., while no defect-defect or Coulombic divergences are present, the nanocrystal must be sufficiently large to allow the system to properly replicate a bulklike configuration. With the development of new algorithms and hardware advances, computations for systems of sufficient size to address this issue are feasible. In particular, we solve the Kohn-Sham equation in real space using pseudopotential-density-functional theory for large silicon nanocrystals, which contain thousands of atoms. Further, we focus on (i) the screening of a point charge and (ii) the formation of a charged vacancy in hydrogen-terminated silicon nanocrystals. This approach allows us to examine the role of quantum confinement in addition to exploring the bulk limit. Comparisons to other methods confirm the viability of this approach.

36 MATERIALS SCIENCE↗

Illustrated formalisms for total scattering data: a guide for new practitioners

The total scattering method is the simultaneous study of both the real- and reciprocal-space representations of diffraction data. While conventional Bragg-scattering analysis (employing methods such as Rietveld refinement) provides insight into the average structure of the material, pair distribution function (PDF) analysis allows for a more focused study of the local atomic arrangement of a material. Generically speaking, a PDF is generated by Fourier transforming the total measured reciprocal-space diffraction data (Bragg and diffuse) into a real-space representation. However, the details of the transformation employed and, by consequence, the resultant appearance and weighting of the real-space representation of the system can vary between different research communities. As the worldwide total scattering community continues to grow, these subtle differences in nomenclature and data representation have led to conflicting and confusing descriptions of how the PDF is defined and calculated. This paper provides a consistent derivation of many of these different forms of the PDF and the transformations required to bridge between them. Some general considerations and advice for total scattering practitioners in selecting and defining the appropriate choice of PDF in their own research are presented. This contribution aims to benefit people starting in the field or trying to compare their results with those of other researchers.

36 MATERIALS SCIENCE↗

CERES : a cryo-EM re-refinement system for continuous improvement of deposited models

The field of electron cryomicroscopy (cryo-EM) has advanced quickly in recent years as the result of numerous technological and methodological developments. This has led to an increase in the number of atomic structures determined using this method. Recently, several tools for the analysis of cryo-EM data and models have been developed within the Phenix software package, such as phenix.real_space_refine for the refinement of atomic models against real-space maps. Also, new validation metrics have been developed for low-resolution cryo-EM models. To understand the quality of deposited cryo-EM structures and how they might be improved, models deposited in the Protein Data Bank that have map resolutions of better than 5 Å were automatically re-refined using current versions of Phenix tools. The results are available on a publicly accessible web page (https://cci.lbl.gov/ceres). The implementation of a Cryo-EM Re-refinement System ( CERES ) for the improvement of models deposited in the wwPDB, and the results of the re-refinements, are described. Based on these results, contents are proposed for a `cryo-EM Table 1', which summarizes experimental details and validation metrics in a similar way to `Table 1' in crystallography. The consistent use of robust metrics for the evaluation of cryo-EM models and data should accompany every structure deposition and be reported in scientific publications.

59 BASIC BIOLOGICAL SCIENCES↗

Imaging Nanoscale Energy Transport and Conversion with Ultrafast Electron Microscopy (Final Technical Report)

Light-matter interactions are ubiquitous in nature and reside at the heart of innumerable technologies. The cascade of processes that occur when a material absorbs a photon of light are exceedingly complex and are interwoven in both space and time, rendering precise determination of the atomic-scale and ultrafast mechanisms immensely challenging. The advent of methods for generating short pulses of light several decades ago led to major advances in understanding the initial moments of light absorption and the resultant effects, though directly interrogating the response of the atoms within the material continued to prove challenging. More recently, methods for generating ultrashort pulses of X-rays and fast electrons have opened the way to probing photoinduced structural dynamics of a wide range of matter in multiple phases. An especially promising laboratory-scale method is ultrafast electron microscopy (UEM), wherein the modalities of conventional transmission electron microscopes (imaging, diffraction, spectroscopy) are extended into the femtosecond temporal regime. Here, the ultrafast imaging and diffraction modalities of UEM were used to study the transient structural aspects of photoexcitation of semiconducting materials (e.g., spatially-resolved electron-phonon coupling, excitation and emission of acoustic phonons, and discrete nanoscale scattering processes). The project had three main objectives: (1) determination of the excitation mechanisms of dense, hypersonic charge-carrier waves and the spatially mediated means by which they couple to the lattice via coherent phonon emission, especially with UEM imaging, wherein effects of nanoscale structural and morphological features on the coupling and relaxation dynamics are expected to dictate nucleation sites and preferred wave vectors, (2) elucidation of photoinduced acoustic-phonon seeding, emergence, propagation, and decay over nanoscale crystal regions and especially with respect to local strain fields and atomic-scale disruptions in lattice order, and (3) interwoven with the first two were efforts aimed at realizing combined angstrom-femtosecond spatiotemporal imaging with UEM. The outcomes and impacts of this project were the generation of new knowledge with respect to fundamental light-matter interactions and, especially, the spatially-mediated excitation and evolution of the structural response of materials following coherent photoexcitation. Importantly, the spatial and temporal resolutions of the UEM imaging modalities used are well-suited for such studies and enable spatially-resolved mechanisms to be determined as a function of atomic order, structural features, and morphology. In addition, ultrafast crystallographic measurements were used to correlate real- and reciprocal-space dynamics in order to determine atomic-scale preferential wave vectors and ultrafast scattering mechanisms, especially as dictated by specimen boundary conditions. The obtained results, as detailed in peer-reviewed publications and presentations, illustrate the importance of ultrafast, angstrom-scale real-space imaging for developing a comprehensive understanding of energy transport and conversion in materials.

47 OTHER INSTRUMENTATION↗

Power spectrum, correlation function, and tests for luminosity bias in the CfA redshift survey

We describe and apply a method for directly computing the power spectrum for the galaxy distribution in the extension of the Center for Astrophysics Redshift Survey. Tests show that our technique accurately reproduces the true power spectrum for k greater than 0.03 h Mpc(exp -1). The dense sampling and large spatial coverage of this survey allow accurate measurement of the redshift-space power spectrum on scales from 5 to approximately 200 h(exp -1) Mpc. The power spectrum has slope n approximately equal -2.1 on small scales (lambda less than or equal 25 h(exp -1) Mpc) and n approximately -1.1 on scales 30 less than lambda less than 120 h(exp -1) Mpc. On larger scales the power spectrum flattens somewhat, but we do not detect a turnover. Comparison with N-body simulations of cosmological models shows that an unbiased, open universe CDM model (OMEGA h = 0.2) and a nonzero cosmological constant (CDM) model (OMEGA h = 0.24, lambda(sub zero) = 0.6, b = 1.3) match the CfA power spectrum over the wavelength range we explore. The standard biased CDM model (OMEGA h = 0.5, b = 1.5) fails (99% significance level) because it has insufficient power on scales lambda greater than 30 h(exp -1) Mpc. Biased CDM with a normalization that matches the Cosmic Microwave Background (CMB) anisotropy (OMEGA h = 0.5, b = 1.4, sigma(sub 8) (mass) = 1) has too much power on small scales to match the observed galaxy power spectrum. This model with b = 1 matches both Cosmic Background Explorer Satellite (COBE) and the small-scale power spect rum but has insufficient power on scales lambda approximately 100 h(exp -1) Mpc. We derive a formula for the effect of small-scale peculiar velocities on the power spectrum and combine this formula with the linear-regime amplification described by Kaiser to compute an estimate of the real-space power spectrum. Two tests reveal luminosity bias in the galaxy distribution: First, the amplitude of the pwer spectrum is approximately 40% larger for the brightest 50% of galaxies in volume-limited samples that have M(sub lim) greater than M*. This bias in the power spectrum is independent of scale, consistent with the peaks-bias paradigm for galaxy formation. Second, the distribution of local density around galaxies shows that regions of moderate and high density contain both very bright (M less than M* = -19.2 + 5 log h) and fainter galaxies, but that voids preferentially harbor fainter galaxies (approximately 2 sigma significance level).

Park, Changbom↗