Search NASA⌕ Search

SEARCH · Search NASA

Results for “Lattice Integrable Models”

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 235 records · Page 13

Calculations of detuning with amplitude for the McMillan electron lens in the Fermilab Integrable Optics Test Accelerator (IOTA)

The McMillan electron lens will be one of the experimental implementations of nonlinear integrable lattices in the Fermilab Integrable Optics Test Accelerator (IOTA). We describe the physics of the nonlinear McMillan lens and calculate the tune dependence with amplitude from the experimental parameters. Some of the implications for the design of experiments in IOTA are discussed.

43 PARTICLE ACCELERATORS↗

Thermodynamics of Dual Doping in Quantum Dots

Dual doping is a powerful way to tailor the properties of semiconductor quantum dots (QDs) arising out of host–dopant and dopant–dopant interactions. Nevertheless, it has seldom been explored due to a variety of thermodynamic challenges, such as the differential bonding strength and diffusion constant within the host matrix that integrates with the host in dissimilar ways. This work discusses the challenges involved in administering them within the constraints of one host under similar conditions of temperature, time, and chemical parameters such as solubility and reactivity using CoPt-doped CdS QDs as a model system. In addition, the various forces in play, such as Kirkendall diffusion, solid- and liquid-state diffusion, hard acid soft base interaction with the host, and the effect of lattice strain due to lattice mismatch, are studied to understand the feasibility of the core to doped transformation. Furthermore, these findings suggest a potential approach for manipulating the properties of semiconductors by dual doping engineering.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

Rolling moments in a trailing vortex flow field

Pressure distributions are presented which were measured on a wing in close proximity to a tip vortex of known structure generated by a larger, upstream semispan wing. Overall loads calculated by integration of these pressures are checked by independent measurements made with an identical model mounted on a force balance. Several conventional methods of wing analysis are used to predict the loads on the following wing. Strip theory is shown to give uniformly poor results for loading distribution, although predictions of overall lift and rolling moment are sometimes acceptable. Good results are obtained for overall coefficients and loading distribution by using linearized pressures in vortex-lattice theory in conjunction with a rectilinear vortex. The equivalent relation from reverse-flow theory that can be used to give economic predictions for overall loads is presented.

Mcmillan, O. J.↗

Rolling moments in a trailing vortex flow field

Pressure distributions are presented which were measured on a wing in close proximity to a tip vortex of known structure generated by a larger, upstream semispan wing. Overall loads calculated by integration of these pressures are checked by independent measurements made with an identical model mounted on a force balance. Several conventional methods of wing analysis are used to predict the loads on the following wing. Strip theory is shown to give uniformly poor results for loading distribution, although predictions of overall lift and rolling moment are sometimes acceptable. Good results are obtained for overall coefficients and loading distribution by using linearized pressures in vortex-lattice theory in conjunction with a rectilinear vortex. The equivalent relation from reverse-flow theory that can be used to give economic predictions for overall loads is presented.

Mcmillan, O. J.↗

Assessment of current state of the art in modeling techniques and analysis methods for large space structures

Advances in continuum modeling, progress in reduction methods, and analysis and modeling needs for large space structures are covered with specific attention given to repetitive lattice trusses. As far as continuum modeling is concerned, an effective and verified analysis capability exists for linear thermoelastic stress, birfurcation buckling, and free vibration problems of repetitive lattices. However, application of continuum modeling to nonlinear analysis needs more development. Reduction methods are very effective for bifurcation buckling and static (steady-state) nonlinear analysis. However, more work is needed to realize their full potential for nonlinear dynamic and time-dependent problems. As far as analysis and modeling needs are concerned, three areas are identified: loads determination, modeling and nonclassical behavior characteristics, and computational algorithms. The impact of new advances in computer hardware, software, integrated analysis, CAD/CAM stems, and materials technology is also discussed.

Noor, A. K.↗

Particle-hole symmetric slave-boson method for the mixed valence problem

We introduce an analytic slave-boson method for treating the finite-𝑈 Anderson impurity model. Our approach introduces two bosons to track both 𝑄 ⇌ 𝑄 ± 1 valence fluctuations and reduces to a single symmetric 𝑠 boson in the effective action, encoding all the high-energy atomic physics information in the boson's kinematics, while the low-energy part of the action remains unchanged across finite-𝑈, infinite-𝑈, and Kondo limits. We recover the infinite-𝑈- and Kondo-limit actions from our approach and show that the Kondo resonance already develops in the normal state when the slave boson has yet to condense. We show that the slave rotor and 𝑠 boson have the same algebraic structure, and we establish a unified functional integral framework connecting the 𝑠-boson and slave-rotor representations for the single-impurity Anderson model.

Anderson impurity model↗

Probing configurational disorder in ZnGeN 2 using cluster-based Monte Carlo

ZnGeN 2 is sought as a semiconductor with comparable lattice constant to GaN and tunable band gap for integration in optoelectronic devices. Configurational disorder on the cation sublattice of ZnGeN 2 can strongly modify the electronic structure compared to the ordered material, and both ordered and disordered forms of ZnGeN 2 are candidates for light-emitting diodes and other emitters. The nonisovalent character of the disordered species (Zn 2+ and Ge 4+ ) subjects the cation ordering to strong short-range order effects. To model these effects, we use Monte Carlo (MC) simulations utilizing a cluster expansion to approximate formation enthalpy. Representative disordered configurations in 1024-atom supercells are relaxed in density functional theory calculations. From the MC structures, we extract a short-range order parameter (the N-cation coordination motif), and two long-range order parameters (Bragg-Williams and stretching parameters), and examine their correlations. We perform a thermodynamic integration to determine the mixing entropy and free energy. ZnGeN 2 exhibits a first-order phase transition with pronounced discontinuities in enthalpy and entropy, as well as in the structural order parameters. We discuss the relationship between the effective temperature used in the MC simulation and the growth temperatures in experiment in relation to the crossover from the nonequilibrium to the equilibrium growth regime. This work expands on current models of site disorder in ZnGeN 2 and provides atomic structure models with a systematic variation of the degree of cation disorder.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Establishing a neutronics design and equilibrium cycle analysis for the I 2 S-LWR reactor with UO 2 and U 3 Si 2 Fuel

The Integral Inherently Safe (I 2 S-LWR) reactor was developed not only to enhance passive reactor safety performance, but also to improve reactor economics. This paper describes the steady-state and depletion analysis of the I 2 S-LWR core using both UO 2 and U 3 Si 2 fuel designs. The development and validation of the core models are described, to include the generation of HELIOS cross sections and comparison with the WEC lattice code PARAGON and the Monte Carlo code SERPENT. The few group homogenized cross sections were then prepared for the U.S. NRC PARCS core simulator using the GenPMAXS code. The multicycle capability in PARCS was used to determine the equilibrium core condition for three batch cores loaded with the UO 2 and then with the U 3 Si 2 fuel designs. Comparing the results from the WEC code ANC at beginning, middle, and end of cycle showed good agreement. Core conditions at beginning and end of cycle were then used to determine the converged steady-state conditions of a coupled RELAP5/PARCS model, which provided the initial condition for various transient accident simulations which are reported in a companion paper.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Nanoparticle self-assemblies with modern complexity

Thanks to decades of tireless efforts, nanoparticle assemblies have reached at an extremely high level of controllability, sophistication, and complexity, with new insights provided by integration with graph theory, cutting-edge characterization, and machine learning (ML)-based computation and modeling, as well as with ever-diversifying applications in energy, catalysis, biomedicine, optics, electronics, magnetics, organic biosynthesis, and quantum technology. Nanoparticle assemblies can be crystalline, known as superlattices or supracrystals. Their assembly entails a transition from disorder—dispersed nanoparticles—to order, which can be achieved through classical nucleation pathways or nonclassical pathways via prenucleation precursors or particle aggregation. Further, the periodic lattices allow facile manipulations of electrons, phonons, photons, and even spins, leading to advanced device components and metamaterials. Meanwhile, aperiodic assemblies out of nanoparticles, such as gels, networks, and amorphous solids, also start to attract attentions. Despite the loss of periodicity, symmetry-lowering or symmetry-breaking three-dimensional (3D) structures emerge with unique properties, such as chiroptical activity, topological mechanical strength, and quantum entanglement. Real-space imaging such as electron microscopy and X-ray based tomography methods are utilized to characterize these complex structures, while mathematical tools such as graph theories are in need to describe such complex structures. This issue aims to provide a timely review of the efforts in this greatly broadened materials design space including experiment, simulation, theory, and applications. Nine top experts (and their teams) from four countries deliver six review papers, summarizing fundamental mechanistic understandings of nanoparticle assemblies, highlighted with the developments of state-of-the-art in situ characterization tools and ML-assisted reverse engineering, and newly emergent applications of nanoarchitectures.

36 MATERIALS SCIENCE↗

Enhancing Kokkos with OpenACC

C++ template metaprogramming has emerged as a prominent approach for achieving performance portability in heterogeneous computing. Kokkos represents a notable paradigm in this domain, offering programmers a suite of high-level abstractions for generic programming while deferring much of the device-specific code generation and optimization to the compiler through template specializations. Kokkos furnishes a range of device-specific code specializations across multiple back ends, including CUDA and HIP. Diverging from conventional back ends, the OpenACC implementation presents a high-level, multicompiler, multidevice, and directive-based programming model. This paper presents recent advancements in the OpenACC back end for Kokkos (i.e., KokkACC) and focuses on its integration into the Kokkos ecosystem, exploration of automatic device selection capabilities to enhance productivity, and performance evaluation on modern hardware such as NVIDIA H100 GPUs. The study includes implementation details and a thorough performance assessment across various computational benchmarks, including minibenchmarks (AXPY and DOT product), miniapps (LULESH, MiniFE, and SNAP-LAMMPS), and a scientific kernel based on the lattice Boltzmann method.

Valero Lara, Pedro↗

Ginzburg-Landau theory of the bcc-liquid interface kinetic coefficient

We extend the Ginzburg-Landau (GL) theory of atomically rough bcc-liquid interfaces outside of equilibrium. We use this extension to derive an analytical expression for the kinetic coefficient, which is the proportionality constant μ(nˆ) between the interface velocity along a direction nˆ normal to the interface and the interface undercooling. The kinetic coefficient is expressed as a spatial integral along the normal direction of a sum of gradient square terms corresponding to different nonlinear density wave profiles. Anisotropy arises naturally from the dependence of those profiles on the angles between the principal reciprocal lattice vectors K → i and nˆ. Values of the kinetic coefficient for the (100), (110), and (111) interfaces are compared quantitatively to the prediction of linear Mikheev-Chernov (MC) theory and previous molecular dynamics (MD) simulation studies of crystallization kinetics for a classical model of Fe. Additional MD simulations are carried out here to compute the relaxation time of density waves in the liquid in order to make this comparison free of fit parameters. The GL theory predicts an expression for μ similar to the MC theory but yields a better agreement with MD simulations for both its magnitude and anisotropy due to a fully nonlinear description of density wave profiles across the solid-liquid interface. In particular, the overall magnitude of μ predicted by GL theory is an order of magnitude larger than predicted by the MC theory. GL theory is also used to derive an inverse relation between μ and the solid-liquid interfacial free energy. Here, the general methodology used here to derive an expression for μ(nˆ) also applies to amplitude equations derived from the phase-field-crystal model, which only differ from GL theory by the choice of cubic and higher order nonlinearities in the free-energy density.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Experimental Assessment of Elastic Modulus vs. Relative Density for Stretch- and Bend-Dominated Lattices

This project was designed to study the possibility of using structural properties of lattices to replicate the material properties of certain hard to manufacture designs and use topology optimization to determine the lattice type and density required to mimic these properties. The integration of additively manufactured lattice structures with topology optimization highlights the need for well characterized mechanical properties and uncertainty analyses to insure these optimized structures respond as predicted. Stereolithographically printed octet and rhombic dodecahedron lattices were manufactured at 10%, 25% and 65% density by volume. As a separate task, yet integrated into this work, finite element analysis (FEA) was used to predict the printed lattice’s mechanical properties, which were then compared to our experimental results. After comparing FEA and experimentallymeasured elastic moduli, it was determined that the FEA provides highly reliable predictions for the modulus of these printed lattice structures. These lattices also exhibited greater tensile stiffness than that of the solid material demonstrating the flexibility that lattices provide to designing parts with designer structural properties. The accurate printing and reliable modeling of these lattices will enable topology optimization of complex parts from well-characterized rhombic dodecahedron and octet lattice structures of varying densities.

36 MATERIALS SCIENCE↗

APOLLO: a facility-scale differentiable virtual accelerator for Fermilab

As the design complexity of modern accelerators grows, there is more interest in using advanced simulations that have fast execution time or yield additional insights like gradients. The FAST/IOTA facility has been working on implementing and experimentally validating an end-to-end digital twin that is both fast and gradient-aware, allowing for rapid prototyping of new software and experiments with minimal beam time costs. Our framework integrates physics and ML codes for linac and ring simulation through a set of generic interfaces between surrogate and physics-based sections. To reproduce device inputs and outputs, system state is exposed as a deterministic discrete event simulator. Because Fermilab is undergoing control system transition, both EPICS and ACNET frontends are supported. Recently, we have begun transitioning to a new community lattice standard, PALS, as well as developing standardized infrastructure for data ingest and normalization to prepare for model calibration during FAST proton injector commissioning. We discuss implementation details as well as challenges, and future plans to extend modelling to main complex proton accelerators like PIPII and Booster.

Kuklev, Nikita [Fermilab]↗

Decoding THz‐Driven Dynamic Fingerprints of Ferroelectric Nanotwin Networks

Ultrafast polarization dynamics in ferroelectrics are of considerable interest for high-speed tunable dielectrics and electro-optics. Extended domain wall networks formed in ferroelectric twin nanodomains can support collective dynamics in the terahertz regime but require techniques that track polarization and strain evolution driven by ultrafast stimulus. Here, we use multi-modal probing of THz-pulse-driven excitations in PbTiO 3 /SrTiO 3 superlattices by combining X-ray free electron laser measurements that directly tracks lattice changes, with optical second harmonic generation that tracks the electronic potential coupled with the lattice potential. Dynamical phase-field modeling enables fingerprinting of these collective modes as superpositions of domain “breathing” through wall oscillations and polarization “rotations” with still walls. Ultrafast domain wall motion at 0.1–0.5 THz is observed at practical fields of 100 kV/cm with wall velocities of >4000 m/s, approaching typical speed of sound in PbTiO 3 . A unique “charging” mode is discovered that can electrically charge and discharge domain walls on ∼4 ps time scale thus dynamically tuning wall conductivity. Integrated experimental and theoretical fingerprinting of the dynamical landscape presented here enables ultrafast control of ferroics for high-speed microelectronics and optical applications.

THz dynamics↗

Advancing simulations of coupled electron and phonon nonequilibrium dynamics using adaptive and multirate time integration

Electronic structure calculations in the time domain provide a deeper understanding of nonequilibrium dynamics in materials. The real-time Boltzmann equation (rt-BTE), used in conjunction with accurate interactions computed from first principles, has enabled reliable predictions of coupled electron and lattice dynamics. However, the timescales and system sizes accessible with this approach are still limited, with two main challenges being the different timescales of electron and phonon interactions and the cost of computing collision integrals. As a result, only a few examples of these calculations exist, mainly for two-dimensional (2D) materials. Here we leverage adaptive and multirate time integration methods to achieve a major step forward in solving the coupled rt-BTEs for electrons and phonons. Relative to conventional (non-adaptive) time-stepping, our approach achieves a 10x speedup for a target accuracy, or greater accuracy by 3–6 orders of magnitude for the same computational cost, enabling efficient calculations in both 2D and bulk materials. This efficiency is showcased by computing the coupled electron and lattice dynamics in graphene up to ~100 ps, as well as modeling ultrafast lattice dynamics and thermal diffuse scattering maps in bulk materials (silicon and gallium arsenide). In addition to improved efficiency, our adaptive method can resolve the characteristic rates of different physical processes, thus naturally bridging different timescales. This enables simulations of longer timescales and provides a framework for modeling multiscale dynamics of coupled degrees of freedom in matter. Our work opens new opportunities for quantitative studies of nonequilibrium physics in materials, including driven lattice dynamics with phonons coupled to electrons, spin, and other degrees of freedom.

Yao, Jia [California Institute of Technology (CalT↗

Competition of lattice and spin excitations in the temperature dependence of spin-wave properties

The interplay of magnons and phonons can induce strong temperature variations in the magnetic exchange interactions, leading to changes in the magnetothermal response. This is a central mechanism in many magnetic phenomena, and in the new field of Spin Caloritronics, which focuses on the combination of heat and spin currents. Boson model systems have previously been developed to describe the magnon-phonon coupling but, until recently, studies rely on empirical parameters. In this paper, we propose a first-principles approach to describe the dependence of the magnetic exchange integrals on phonon renormalization, leading to changes in the magnon dispersion as a function of temperature. The temperature enters into the spin dynamics (by introducing fluctuations) as well as in the magnetic exchange itself. Depending on the strength of the coupling, these two temperatures may or may not be equilibrated, yielding different regimes. We test our approach in typical and well-known ferromagnetic materials: Ni, Fe, and Permalloy. Finally, we compare our results to recent experiments on the spin-wave stiffness, and discuss departures from Bloch’s law and parabolic dispersion.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Mesoscopic Elastic Distortions in GaAs Quantum Dot Heterostructures

Quantum devices formed in high-electron-mobility semiconductor heterostructures provide a route through which quantum mechanical effects can be exploited on length scales accessible to lithography and integrated electronics. The electrostatic definition of quantum dots in semiconductor heterostructure devices intrinsically involves the lithographic fabrication of intricate patterns of metallic electrodes. The formation of metal/semiconductor interfaces, growth processes associated with polycrystalline metallic layers, and differential thermal expansion produce elastic distortion in the active areas of quantum devices. Understanding and controlling these distortions present a significant challenge in quantum device development. In this paper, we report synchrotron X-ray nanodiffraction measurements combined with dynamical X-ray diffraction modeling that reveal lattice tilts with a depth-averaged value up to 0.04° and strain on the order of 10 –4 in the two-dimensional electron gas (2DEG) in a GaAs/AlGaAs heterostructure. Elastic distortions in GaAs/AlGaAs heterostructures modify the potential energy landscape in the 2DEG due to the generation of a deformation potential and an electric field through the piezoelectric effect. Finally, the stress induced by metal electrodes directly impacts the ability to control the positions of the potential minima where quantum dots form and the coupling between neighboring quantum dots.

36 MATERIALS SCIENCE↗

Ergodic and nonergodic many-body dynamics in strongly nonlinear lattices

The study of nonlinear oscillator chains in classical many-body dynamics has a storied history going back to the seminal work of Fermi et al. [Los Alamos Scientific Laboratory Report No. LA-1940, 1955 (unpublished)]. Here, we introduce a family of such systems which consist of chains of N harmonically coupled particles with the nonlinearity introduced by confining the motion of each individual particle to a box or stadium with hard walls. The stadia are arranged on a one-dimensional lattice but they individually do not have to be one dimensional, thus permitting the introduction of chaos already at the lattice scale. For the most part we study the case where the motion is entirely one dimensional. We find that the system exhibits a mixed phase space for any finite value of N . Computations of Lyapunov spectra at randomly picked phase space locations and a direct comparison between Hamiltonian evolution and phase space averages indicate that the regular regions of phase space are not significant at large system sizes. While the continuum limit of our model is itself a singular limit of the integrable sinh Gordon theory, we do not see any evidence for the kind of nonergodicity famously seen in the work of Fermi et al. Finally, we examine the chain with particles confined to two-dimensional stadia where the individual stadium is already chaotic and find a much more chaotic phase space at small system sizes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗