Search NASA⌕ Search

SEARCH · Search NASA

Results for “finite-size scaling”

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 55 records · Page 3

Neutral band gap of carbon by quantum Monte Carlo methods

We present a method of calculating the energy gap of a charge-neutral excitation using only ground-state calculations. We report Quantum Monte Carlo calculations of Γ→ Γ and Γ → X particle-hole excitation energies in diamond carbon. We analyze the finite-size effect and find the same 1/L decay rate as that in a charged excitation, where L is the linear extension of the supercell. This slow decay is attributed to the delocalized nature of the excitation in supercells too small to accommodate excitonic binding effects. At larger system sizes, the apparent 1/L decay crosses over to a 1/L 3 behavior. Estimation of the scale of exciton binding can be used to correct finite-size effects of neutral gaps.

36 MATERIALS SCIENCE↗

Electrostatic Kelvin-Helmholtz instability in a radially injected plasma cloud

An elecrostatic finite-sized particle simulation model is used to study the early time-scale phenomena asociated with an impulsively injected plasma, expanding radially, normal to a strong ambient magnetic field. Results of the simulation show the early formation of a radial polarization electric field due to ion-electron charge separation which causes electrons to drift in the azimuthal direction. Velocity shear within this motion gives rise to a Kelvin-Helmholtz (dioxotron) instability, creating a radial fluted pattern as the plasma expands. In the nonlinear stage of instability, azimuthal variations in the electric field are observed which cause electrons to drift across the magnetic field, thereby reducing the space charge created by energetic ion expansion. The results of a linear stability analysis reveal that for a decreasing amount of charge separation in the plasma, the number of unstable azimuthal modes at maximum growth rate increases. Also, as the thickness of the electron ring increases, a fewer number of unstable modes develops for a given amount of charge separation. It is concluded that electrons become unstable at the very early time scale, within an ion gyroperiod, and that unstable electrons lead to modification of the ion dynamics.

Sydora, R. D.↗

Trajectory sampling and finite-size effects in first-principles stopping power calculations

Abstract Real-time time-dependent density functional theory (TDDFT) is presently the most accurate available method for computing electronic stopping powers from first principles. However, obtaining application-relevant results often involves either costly averages over multiple calculations or ad hoc selection of a representative ion trajectory. We consider a broadly applicable, quantitative metric for evaluating and optimizing trajectories in this context. This methodology enables rigorous analysis of the failure modes of various common trajectory choices in crystalline materials. Although randomly selecting trajectories is common practice in stopping power calculations in solids, we show that nearly 30% of random trajectories in an FCC aluminum crystal will not representatively sample the material over the time and length scales feasibly simulated with TDDFT, and unrepresentative choices incur errors of up to 60%. We also show that finite-size effects depend on ion trajectory via “ouroboros” effects beyond the prevailing plasmon-based interpretation, and we propose a cost-reducing scheme to obtain converged results even when expensive core-electron contributions preclude large supercells. This work helps to mitigate poorly controlled approximations in first-principles stopping power calculations, allowing 1–2 order of magnitude cost reductions for obtaining representatively averaged and converged results.

36 MATERIALS SCIENCE↗

Multiscale simulation of spatially correlated microstructure via a latent space representation

When deformation gradients act on the scale of the microstructure of a part due to geometry and loading, spatial correlations and finite-size effects in simulation cells cannot be neglected. We propose a multiscale method that accounts for these effects using a variational autoencoder to encode the structure–property map of the stochastic volume elements making up the statistical description of the part. In this paradigm the autoencoder can be used to directly encode the microstructure or, alternatively, its latent space can be sampled to provide likely realizations. Furthermore, we demonstrate the method on three examples using the common additively manufactured material AlSi10Mg in: (a) a comparison with direct numerical simulation of the part microstructure, (b) a push forward of microstructural uncertainty to performance quantities of interest, and (c) a simulation of functional gradation of a part with stochastic microstructure.

Elastoplasticity↗

Itinerant ferromagnetism in transition metal dichalcogenide moiré superlattices

Moiré materials are artificial crystals formed at van der Waals heterojunctions that have emerged as a highly tunable platform that is able to realize much of the rich quantum physics of electrons in atomic scale solids, and in several cases even new quantum phases of matter. Here we use finite-size exact diagonalization methods to explore the physics of single-band itinerant electron ferromagnetism in semiconductor moiré materials. As a result, we predict where ferromagnetism is likely to occur in triangular-lattice moiré systems, and where it is likely to yield the highest Curie temperatures.

36 MATERIALS SCIENCE↗

Finite-size effects on long-range correlations: implications for analyzing DNA sequences

We analyze the fluctuations in the correlation exponents obtained for noncoding DNA sequences. We find prominent sample-to-sample variations as well as variations within a single sample in the scaling exponent. To determine if these fluctuations may result from finite system size, we generate correlated random sequences of comparable length and study the fluctuations in this control system. We find that the DNA exponent fluctuations are consistent with those obtained from the control sequences having long-range power-law correlations. Finally, we compare our exponents for the DNA sequences with the exponents obtained from power-spectrum analysis and correlation-function techniques, and demonstrate that the original "DNA-walk" method is intrinsically more accurate due to reduced noise.

NASA Discipline Cardiopulmonary↗

Cavity electrodynamics of van der Waals heterostructures

Van der Waals heterostructures host many-body quantum phenomena that are tunable in situ using electrostatic gates. Their constituent two-dimensional materials and gates can naturally form plasmonic self-cavities, confining light in standing waves of current density due to finite-size effects. The plasmonic resonances of typical graphite gates fall in the gigahertz to terahertz range, corresponding to the same microelectronvolt to millielectronvolt energy scale as the phenomena in van der Waals heterostructures that they electrically control. This raises the possibility that the built-in cavity modes of graphite gates are relevant for shaping the low-energy physics of these heterostructures. However, probing these cavity-coupled electrodynamics is challenging as devices are notably smaller than the diffraction limit at the relevant wavelengths. Here we report on the intrinsic cavity conductivity of gate-tunable graphene heterostructures. As the carrier density is tuned, we observe coupling and spectral weight transfer between graphene and graphite plasmonic cavity modes in the ultrastrong coupling regime. We present an analytical model to describe the results and provide general principles for cavity design. Our findings show that intrinsic cavity effects are important for understanding the low-energy electrodynamics of van der Waals heterostructures and open a pathway for useful functionality through cavity control.

Electronic properties and materials↗

Generalized design principles for hydrodynamic electron transport in anisotropic metals

Interactions of charge carriers with lattice vibrations, or phonons, play a critical role in unconventional electronic transport of metals and semimetals. Recent observations of phonon-mediated collective electron flow in bulk semimetals, termed electron hydrodynamics, present new opportunities in the search for strong electron-electron interactions in high carrier density materials. Here we present the general transport signatures of such a second-order scattering mechanism, along with analytical limits at the Eliashberg level of theory. We study electronic transport, using ab initio calculations, in finite-size channels of semimetallic ZrSiS and TaAs 2 with and without topological band crossings, respectively. The order of magnitude separation between momentum-relaxing and momentum-conserving scattering length scales across a wide temperature range make both of them promising candidates for further experimental observation of electron hydrodynamics. More generally, our calculations suggest that the hydrodynamic transport regime does not, to first order, rely on the topological nature of the bands. Finally, we discuss general design principles guiding future search for hydrodynamic candidates, based on the analytical formulation and our ab initio predictions. Further, we find that systems with strong electron-phonon interactions, reduced electronic phase space, and suppressed phonon-phonon scattering at temperatures of interest are likely to feature hydrodynamic electron transport. We predict that layered and/or anisotropic semimetals composed of half-filled d shells and light group V/VI elements with lower crystal symmetry are promising candidates to observe hydrodynamic phenomena in the future.

36 MATERIALS SCIENCE↗

Numerical simulation of nonaxisymmetric adiabatic accretion flow

The hydrodynamics of gas flow past a finite-sized gravitating central object is studied in the adiabatic approximation in two spatial dimensions. The flow morphology is shown to depend on the ratio (epsilon) of the axisymmetric accretion radius to the length scale of the density variation. For small to intermediate epsilon values the flows are highly time dependent, and for large epsilon values a quasi-steady state is obtained.

Fryxell, B. A.↗

Scale setting of SU⁡(𝑁) Yang–Mills theory, topology and large-𝑁 volume independence

We set the scale of SU⁡(𝑁) Yang-Mills theories for 𝑁 =3, 5, 8 and in the large-𝑁 limit via gradient flow, as a first step towards the computation of the large-𝑁 Λ-parameter using step scaling. We adopt twisted boundary conditions to achieve large-𝑁 volume reduction and the Parallel Tempering on Boundary Conditions algorithm to tame topological freezing. This setup allows accurate determinations of the gradient-flow scales down to lattice spacings as fine as ∼0.025 fm for all the explored values of 𝑁, a regime that has never been reached with ergodic algorithms. Moreover, we are able to precisely estimate the finite-size systematics related to topological freezing, and to show the suppression of finite-volume effects expected by virtue of large-𝑁 twisted volume reduction.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Molecular concentration field design using closed-form steady-state solutions

Control over spatial concentration fields represents a fundamental challenge in designing synthetic biological systems and programmable soft materials. While nature creates morphogen gradients that orchestrate complex developmental processes, synthetic approaches have largely relied on empirical optimization and computationally intensive simulations. Here, we present an analytical framework for steady-state concentration fields generated by finite-sized localized sources in diffusion–degradation systems and derive closed-form solutions for one-, two-, and three-dimensional geometries. By expressing these solutions in dimensionless form, we show that gradient steepness and spatial structure are organized by the Thiele modulus, which captures the competition between diffusion and degradation length scales. The analysis reveals distinct design regimes: in degradation-dominated systems, gradient shape is governed by exponential decay and becomes dimension-independent, whereas in diffusion-dominated systems, gradient magnitude and extent follow dimension-dependent power-law scaling. Building on these results, we introduce a quantitative design strategy that uses threshold-based criteria to program concentration ranges by tuning physically accessible parameters, most directly the production rate, while holding transport and degradation properties fixed. Comparisons with numerical solutions and reported experimental systems demonstrate consistency with the predicted scaling behavior. Together, this work provides a generalizable and physically transparent framework for designing steady-state concentration fields in synthetic biological and soft matter systems, enabling predictive control of gradient-mediated organization without reliance on extensive numerical optimization.

Kim, Dong Woo [Johns Hopkins University, Baltimore↗

On nonsteady accretion in stellar wind-fed X-ray sources

The paper reports on the results of hydrodynamical simulations of nonaxisymmetric gas flow past a finite-sized gravitating object. The asymmetry in the flow is assumed to be due to a transverse density gradient. The solutions are found to be highly time-dependent and exhibit phases in which a disk forms. The direction of circulation of the disk alternates between the direct and counter sense on time scales of approximately 1 hr. While the disk is present, the specific angular momentum is high, and the mass accretion rate, low. The angular momentum capture oscillates about a mean which secularly increases. Such temporal behavior may be relevant to those X-ray pulsing systems (e.g., Vela X-1) which exhibit fluctuations in the X-ray luminosity and pulse period on time scales much shorter than the orbital period of the binary system.

Taam, Ronald E.↗

A ductility metric for refractory-based multi-principal-element alloys

We propose a quantum-mechanical dimensionless metric, the local-lattice distortion (LLD), as a reliable predictor of ductility in refractory multi-principal-element alloys (RMPEAs). The LLD metric is based on electronegativity differences in localized chemical environments and combines atomic-scale displacements due to local lattice distortions with a weighted average of valence-electron count. To evaluate the effectiveness of this metric, we examined body-centered cubic (bcc) refractory alloys that exhibit ductile-to-brittle behavior. Our findings demonstrate that local-charge behavior can be tuned via composition to enhance ductility in RMPEAs. With finite-sized cell effects eliminated, the LLD metric accurately predicted the ductility of arbitrary alloys, which compares well with existing tensile-elongation experiments. To validate further, we qualitatively evaluated the ductility of two refractory RMPEAs, i.e., NbTaMoW and Mo 72 W 13 Ta 10 Ti 2.5 Zr 2.5 , through the observation of crack formation under indentation, again showing excellent agreement with LLD predictions. Additionally, a comparative study of three refractory alloys provides further insights into the electronic-structure origin of ductility in refractory RMPEAs. This proposed metric enables rapid and accurate assessment of ductility behavior in the vast RMPEA composition space.

36 MATERIALS SCIENCE↗

A Dipolar Chiral Spin Liquid on the Breathed Kagome Lattice

Continuous control over lattice geometry, when combined with long-range interactions, offers a powerful yet underexplored tool to generate highly frustrated quantum spin systems. By considering long-range dipolar antiferromagnetic interactions on a breathed Kagome lattice, we demonstrate how these tools can be leveraged to stabilize a chiral spin liquid. We support this prediction with large-scale density-matrix renormalization group calculations and explore the surrounding phase diagram, identifying a route to adiabatic preparation via a locally varying magnetic field. At the same time, we identify the relevant low-energy degrees of freedom in each unit cell, providing a complementary language to study the chiral spin liquid. Finally, we carefully analyze its stability and signatures in finite-sized clusters, proposing direct, experimentally viable measurements of the chiral edge mode in both Rydberg atom and ultracold polar molecule arrays.

Machado↗

Resolving Nonequilibrium Shape Variations among Millions of Gold Nanoparticles

Nanoparticles, exhibiting functionally relevant structural heterogeneity, are at the forefront of cutting-edge research. Now, high-throughput single-particle imaging (SPI) with X-ray free-electron lasers (XFELs) creates opportunities for recovering the shape distributions of millions of particles that exhibit functionally relevant structural heterogeneity. To realize this potential, three challenges have to be overcome: (1) simultaneous parametrization of structural variability in real and reciprocal spaces; (2) efficiently inferring the latent parameters of each SPI measurement; (3) scaling up comparisons between 10 5 structural models and 10 6 XFEL-SPI measurements. Here, we describe how we overcame these three challenges to resolve the nonequilibrium shape distributions within millions of gold nanoparticles imaged at the European XFEL. These shape distributions allowed us to quantify the degree of asymmetry in these particles, discover a relatively stable “shape envelope” among nanoparticles, discern finite-size effects related to shape-controlling surfactants, and extrapolate nanoparticles’ shapes to their idealized thermodynamic limit. Ultimately, these demonstrations show that XFEL SPI can help transform nanoparticle shape characterization from anecdotally interesting to statistically meaningful.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Stability, metallicity, and magnetism in niobium silicide nanofilms

Modern superconducting qubits based on two-dimensional (2D) transmons typically involve the growth of Nb thin films on high-resistivity Si substrates. Since imperfections at the Nb-Si heterointerface have been implicated as a source of two-level systems that limit quantum coherence times, detailed characterization and understanding of niobium silicide interfacial layers are critical to improving superconducting qubit technology. While bulk binary intermetallic niobium silicide phases are well understood, the thermodynamic phase stability and properties of ultrathin niobium silicides, such as those found at the Nb-Si heterointerface in 2D transmons, have not yet been explored. Here, we report finite-sized effects for ultrathin niobium silicide films using density functional theory calculations and predict nanoscale stabilization of Nb 6 Si 5 over the bulk α-Nb 5 Si 3 phase. This result is consistent with our experimental observations of a niobium silicide interfacial layer between a sputtered Nb thin film and the underlying Si substrate. Furthermore, our calculations show that Nb 6 Si 5 nanofilms are nonmagnetic, making them superior to nanofilms of α-Nb 5 Si 3 that exhibit antiferromagnetic correlations detrimental to long coherence times in superconducting qubits. Furthermore, by providing atomic-scale insight into niobium silicide nanofilms, this paper can help guide ongoing efforts to optimize Nb-Si heterointerfaces for long coherence times in superconducting qubits.

36 MATERIALS SCIENCE↗

Bridging the time scale in exascale computing of chemical systems (Final Technical Report)

This report summarizes the work carried out with support of the United States Department of Energy under Award DE-SC0019441. The theme of this project was to develop and apply methods that allowed for the acceleration of atomistic calculations, particularly in challenging areas such as multiphase systems, electrified interfaces, uncertainty estimation, and applications requiring chemical accuracy, which tend to be applications where simulation time is severely bottlenecked by the computational time requirements. Much of the focus was on the application of emerging machine-learning methodologies, although a wide range of methodologies were employed. This report has two major sections. The first focuses on the methodological advances themselves. Within this part, we report a number of major advances, a few examples of which are described here. We report the first machine-learning scheme for the acceleration of electronically grand-canonical calculations (that is, those applicable to electrochemistry). We report new methods of performing transfer learning, in which physics-based priors can be used to provide predictions, often with uncertainty estimates, of images well outside of training sets; we also offer ways to fine-tune these transfer-learning models. We provide a new systematic means to generate and apply minimal training data sets to very large (10,000’s of atoms) systems, with only small training sets appropriate for electronic structure. We developed new methodologies to integrate surface vibrations into surface adsorption calculations. We made advances to the applicability of diffusion Monte Carlo methods to allow (learned) force prediction, finite-size error correction, and force-free means of searching for transition states. We integrated machine-learned atomistic predictions into mechanism generation codes. Additionally, we released new software including AmpTorch, a modernized version of our original atomistic machine-learning code Amp. The second part of this report focuses on the scientific applications that accompanied, and were often enabled by, the methodological advances described earlier. A few examples follow, but full details are in the individual chapters of the report. For example, we developed a general theory of phonon-induced friction on molecular adsorbates. We showed fundamentally how solvent influences the adsorption and desorption process and how it differs from the processes typically involved at the solid–gas interface, making aqueous-phase and electrocatalysis different from traditional thermocatalysis. We examined how metal–insulator and magnetic transitions can be probed, and accelerated exciton dynamics via Frenkel Hamiltonian parameters. We showed that the nearsighted force-training approach, developed within this project, can predict both the stability and reactivity of large nanoparticles, and can also lead to insights on catalyst coverage on binding energies and entropies. These applied studies, which generally integrated with our method development, allowed us to push forward the theoretical understanding of several reaction classes.

08 HYDROGEN↗

Finite-size subthermal regime in disordered SU ( N ) -symmetric Heisenberg chains

SU⁡(N) symmetry is incompatible with the many-body localized (MBL) phase, even when strong disorder is present. However, recent studies have shown that finite-size SU⁡(2) systems exhibit nonergodic, subthermal behavior, characterized by the breakdown of the eigenstate thermalization hypothesis, and by the excited eigenstates entanglement entropy that is intermediate between area and volume law. In this paper, we extend previous studies of the SU⁡(2)-symmetric disordered Heisenberg model to larger systems, using the time-dependent density matrix renormalization group (tDMRG) method. We simulate quench dynamics from weakly entangled initial states up to long times, finding robust subthermal behavior at stronger disorder. Although we find an increased tendency towards thermalization at larger system sizes, the subthermal regime persists at intermediate time scales, nevertheless, and therefore should be accessible experimentally. At weaker disorder, we observe signatures of thermalization; however, entanglement entropy exhibits slow sublinear growth, in contrast to conventional thermalizing systems. Furthermore, we study dynamics of the SU⁡(3)-symmetric disordered Heisenberg model. Similarly, strong disorder drives the system into subthermal regime, albeit thermalizing phase is broader compared to the SU⁡(2) case. Finally, our findings demonstrate the robustness of the subthermal regime in spin chains with non-Abelian continuous symmetry, and are consistent with eventual thermalization at large system sizes and long time scales, suggested by previous studies.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗