SEARCH · Engineering Papers
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.
Controlling dimensionality and strain in actinide heterojunctions will provide unrivaled opportunities for exploring novel quantum phenomena. We discuss the promises, challenges, and synthesis routes for these actinide-bearing heterostructures with complex electron correlation for functional and energy materials.
Mott materials are archetypal quantum systems actively explored as next-generation electronic and photonic platforms, with potential applications spanning non-Von Neumann computing, robotics, energy storage, and microwave technologies. Among these, vanadium dioxide (VO 2 ) has emerged as one of the most intensively studied compounds, owing to its sharp, near-room-temperature insulator-to-metal phase transition. VO 2 also serves as a benchmark system for testing cutting-edge theories and experimental techniques. Here, we directly visualize the electrically driven transition dynamics in VO 2 using a microwave-driven, frequency-tunable pulsed transmission electron microscope that combines nanometer spatial and picosecond temporal resolution. Under high-frequency (MHz–GHz) excitation, we capture the ultrafast nucleation, propagation, and dissolution of metallic domains within an operating device over millions of reversible cycles. We observe the ultrafast formation of consistent metallic nuclei beneath the electrodes, followed by the propagation of a structural phase front at 4.54 nm/ns. Our experiments show that phonon-mediated structural recovery ultimately limits reversible switching of VO 2 at GHz frequencies, and that a tunable regime for reversible operation spans from kHz to GHz through device engineering. Beyond VO 2 , our approach provides a powerful framework for probing non-equilibrium structural transformations in correlated and functional materials under realistic electrical stimuli.
Magnetic field aligned filaments such as blobs and edge localized mode filaments carry significant amounts of heat and particles to the plasma facing components and they decrease their lifetime. The dynamics of these filaments determine at least a part of the heat and particle loads. These dynamics can be characterized by their translation and rotation. In this paper, we present an analysis method novel for fusion plasmas, which can estimate the angular velocity of the filaments on frame-by-frame time resolution. After pre-processing, the frames are two-dimensional (2D) Fourier-transformed, then the resulting 2D Fourier magnitude spectra are transformed to log-polar coordinates, and finally the 2D cross-correlation coefficient function (CCCF) is calculated between the consecutive frames. The displacement of the CCCF’s peak along the angular coordinate estimates the angle of rotation of the most intense structure in the frame. Further, the proposed angular velocity estimation method is tested and validated for its accuracy and robustness by applying it to rotating Gaussian-structures. The method is also applied to gas-puff imaging measurements of filaments in National Spherical Torus Experiment plasmas.
Here, we present a fast, asymptotically linear-scaling implementation of the perturbative quadruples energy correction in coupled-cluster theory using local natural orbitals. Our work follows the domain-based local pair natural orbital (DLPNO) approach previously applied to lower levels of excitations in coupled-cluster theory. Our DLPNO-CCSDT(Q) algorithm uses converged doubles and triples amplitudes from a preceding DLPNO-CCSDT computation to compute the quadruples amplitude and energy in the quadruples natural orbital (QNO) basis. We demonstrate the compactness of the QNO space, showing that more than 95% of the (Q) correction can be recovered using relatively loose natural orbital cutoffs, compared to the tighter cutoffs used in pair and triples natural orbitals at lower levels of coupled-cluster theory. We also highlight the accuracy of our algorithm in the computation of relative energies, which yields deviations of sub-kJ mol −1 in relative energy compared to the canonical CCSDT(Q). Timings are conducted on a series of growing linear alkanes (up to 10 carbons and 608 basis functions) and water clusters (up to 49 water molecules and 2842 basis functions) to establish the asymptotic linear-scaling of our DLPNO-(Q) algorithm.
The accurate prediction of band gaps and structural properties in periodic systems continues to be one of the central goals of electronic structure theory. However, band gaps obtained from popular exchange–correlation (XC) functionals (such as LDA and PBE) are severely underestimated partly due to the spurious self-interaction error (SIE) inherent to these functionals. In this work, we present a new formulation and implementation of Wannier function-derived Fermi–Löwdin (WFL) orbitals for correcting the SIE in periodic systems. Since our approach utilizes a variational minimization of the self-interaction energy with respect to the Wannier charge centers (WCC), it is computationally more efficient than the HSE hybrid functional and other self-interaction corrections that require a large number of transformation matrix elements. Here, calculations on several (17 in total) prototypical molecular solids, semiconductors, and wide-bandgap materials show that our WFL self-interaction correction approach gives better band gaps and bulk moduli compared to semilocal functionals, largely due to the partial removal of self-interaction errors.
Recently, transition-metal-based kagome metals have aroused much research interest as a novel platform to explore exotic topological quantum phenomena. Here we report on the synthesis, structure, and physical properties of a bilayer kagome lattice compound V 3 Sb 2 . The polycrystalline V 3 Sb 2 samples were synthesized by conventional solid-state-reaction method in a sealed quartz tube at temperatures below 850 °C. Measurements of magnetic susceptibility and resistivity revealed consistently a density-wave-like transition at T dw ≈ 160 K with a large thermal hysteresis, even though some sample-dependent behaviors were observed presumably due to the different preparation conditions. Upon cooling through T dw , no strong anomaly in lattice parameters and no indication of symmetry lowering were detected in powder x-ray diffraction measurements. This transition can be suppressed completely by applying hydrostatic pressures of about 1.8 GPa, around which no sign of superconductivity was observed down to 1.5 K. Specific-heat measurements revealed a relatively large Sommerfeld coefficient γ = 18.5 mJ∙mol –1 ∙K –2 , confirming the metallic ground state with moderate electronic correlations. Density functional theory calculations indicate that V 3 Sb 2 shows a non-trivial topological crystalline property. Furthermore, our study makes V 3 Sb 2 a new candidate of metallic kagome compound to study the interplay between density-wave-order, nontrivial band topology, and possible superconductivity.
Recent laser-shock experiments on an H–He mixture containing 11% helium (atomic fraction) have suggested the presence of an immiscibility region inside Jupiter. Reflectivity measurements were used as the primary diagnostic of H–He demixing, with discontinuities in the optical reflectivity proposed as a signature of phase separation under conditions relevant to Jupiter's interior. Here, we investigate shock-compressed H–He using ab initio molecular dynamics simulations with optical properties evaluated within the Kubo–Greenwood formalism. The equation of state and ionic configurations were obtained using the thermal Tr 2 SCANL meta-GGA exchange–correlation (XC) functional, while optical properties were computed using the recently developed RS-KDT0 range-separated thermal hybrid XC, which provides state-of-the-art accuracy for band-gap predictions in the warm dense matter regime. The calculated reflectivity shows overall good agreement with experimental measurements; however, no discontinuity is observed at elevated temperatures. Moreover, the reflectivity predictions for the mixed system are consistent with the experimental measurements in the temperature range where the mixture is inferred to be demixed. Furthermore, these results suggest that reflectivity alone may not provide a unique or sensitive diagnostic of H–He demixing at low helium concentrations under these conditions.
The most stringent laboratory-based experimental limits on the existence of submegaelectronvolt sterile neutrinos are currently set by decay spectroscopy of radioactive 7 Be embedded into superconducting sensors. The systematic uncertainties are dominated by the modeling of the electron shakeup and shakeoff spectra that are not based on state-of-the-art atomic theory and do not include electron correlations or relativistic effects. We have used the multiconfiguration Dirac-Fock formalism to obtain correlated wave functions ab initio and compute all single and double shake processes in the electron capture decay of atomic 7 Be . The simulations can explain some but not all of the observed spectral features, likely because the wave functions are modified by the Ta sensor material that the 7 Be is embedded into. The new models also show that the 𝐿/𝐾 electron capture ratio of 7 Be in Ta has previously been slightly underestimated revising the previous value of 0.070(7) to a new value of 0.0756(20).
A distinct feature of Hermitian quantum chaotic dynamics is the exponential increase of certain out-of-time-order correlation (OTOC) functions around the Ehrenfest time with a rate given by a Lyapunov exponent. Physically, the OTOCs describe the growth of quantum uncertainty that crucially depends on the nature of the quantum motion. Here, we employ the OTOC in order to provide a precise definition of dissipative quantum chaos. For this purpose, we compute analytically the Lyapunov exponent for the vectorized formulation of the large- q limit of a q -body Sachdev-Ye-Kitaev model coupled to a Markovian bath. These analytic results are confirmed by an explicit numerical calculation of the Lyapunov exponent for several values of q ≥ 4 based on the solutions of the Schwinger-Dyson and Bethe-Salpeter equations. We show that the Lyapunov exponent decreases monotonically as the coupling to the bath increases and eventually becomes negative at a critical value of the coupling signaling a transition to a dynamics which is no longer quantum chaotic. Therefore, a positive Lyapunov exponent is a defining feature of dissipative many-body quantum chaos. The observation of the breaking of the exponential growth for sufficiently strong coupling suggests that dissipative quantum chaos may require in certain cases a sufficiently weak coupling to the environment. Published by the American Physical Society 2024
Quantitative convergent beam electron diffraction (CBED) enables determination of aspherical valence electron distributions through refinement of low-order structure factors, which are highly sensitive to chemical bonding and charge density variations. However, conventional quantitative CBED (QCBED) requires solving a highly nonlinear inverse problem with many coupled parameters, and computationally intensive dynamical diffraction calculations, making it time-consuming and difficult to apply to complex systems. More broadly, reconstructing charge density and orbital electron distribution from diffraction data has long been a central challenge in both x-ray and electron crystallography. Here, in this study, we introduce an artificial-intelligence (AI)-based framework that replaces traditional refinement with a data-driven inverse solver. Using a large synthetic CBED dataset generated by Bloch-wave simulations, we train a conditional diffusion model to directly infer crystal structural parameters and multipole density formalism parameters, and hence valence electron distributions, from CBED patterns alone. By learning from forward simulations across realistic parameter space, the model effectively solves the inverse problem. Compared with direct regression approaches, the diffusion-based framework provides posterior parameter distributions for rigorous uncertainty quantification while preserving quantitative fidelity and reducing analysis time by orders of magnitude. By eliminating the need for external single-crystal x-ray diffraction data and complex nonlinear refinement, this approach enables practical, high-throughput, and in situ quantitative CBED, enabling real-time mapping of valence electron distributions and their correlation with functional responses in quantum and energy materials.
X-ray free-electron lasers (X-FELs) present new opportunities to study ultrafast lattice dynamics in complex materials. While the unprecedented source brilliance enables high fidelity measurement of structural dynamics, it also raises experimental challenges related to the understanding and control of beam-induced irreversible structural changes in samples that can ultimately impact the interpretation of experimental results. This is also important for designing reliable high performance X-ray optical components. In this work, X-FEL beam-induced lattice alterations are investigated by measuring the shot-to-shot evolution of near-Bragg coherent scattering from a single crystalline germanium sample. It is shown that X-ray photon correlation analysis of sequential speckle patterns measurements can be used to monitor the nature and extent of lattice rearrangements. Additionally, abrupt, irreversible changes are observed following intermittent high-fluence monochromatic X-ray pulses, thus revealing the existence of a threshold response to X-FEL pulse intensity.
ExchCXX is a modern C++ library for the evaluation of exchange-correlation (XC) functionals required for density functional theory (DFT) calculations. In addition to providing high-level wrappers for the ubiquitous Libxc library for the CPU evaluation of XC functionals, ExchCXX also provides locally developed implementations of a small subset of XC functionals which may be evaluated either on the host (CPU) or device (GPU, FPGA, etc).
Filamentous fungi degrade complex plant material to its monomeric building blocks, which have many biotechnological applications. Transcription factors play a key role in plant biomass degradation, but little is known about their interactions in the regulation of polysaccharide degradation. Here, we deepened the knowledge about the storage polysaccharide regulators AmyR and InuR in Aspergillus niger. AmyR controls starch degradation, while InuR is involved in sucrose and inulin utilization. In our study, the phenotypes of A. niger parental, ΔamyR, ΔinuR and ΔamyRΔinuR strains were assessed in both solid and liquid media containing sucrose or inulin as carbon source to evaluate the roles of AmyR and InuR and the effect of culture conditions on their functions. In correlation with previous studies, our data showed that AmyR has a minor contribution to sucrose and inulin utilization when InuR is active. In contrast, growth profiles and transcriptomic data showed that the deletion of amyR in the ΔinuR background strain resulted in more pronounced growth reduction on both substrates, mainly evidenced by data originating from solid cultures. Overall, our results show that submerged cultures do not always reflect the role of transcription factors in the natural growth condition, which is better represented on solid substrates. Importance: The type of growth has critical implications in enzyme production by filamentous fungi, a process that is controlled by transcription factors. Submerged cultures are the preferred setups in laboratory and industry and are often used for studying the physiology of fungi. In this study, we showed that the genetic response of A. niger to starch and inulin was highly affected by the culture condition, since the transcriptomic response obtained in a liquid environment did not fully match the behavior of the fungus in a solid environment. These results have direct implications in enzyme production and would help industry choose the best approaches to produce specific CAZymes for industrial purposes.
We develop methods to learn the correlation potential for a time-dependent Kohn-Sham (TDKS) system in one spatial dimension. We start from a low-dimensional two-electron system for which we can numerically solve the time-dependent Schr¨odinger equation; this yields electron densities suitable for training models of the correlation potential. We frame the learning problem as one of optimizing a least-squares objective subject to the constraint that the dynamics obey the TDKS equation. Applying adjoints, we develop efficient methods to compute gradients and thereby learn models of the correlation potential. Our results show that it is possible to learn values of the correlation potential such that the resulting electron densities match ground truth densities. We also show how to learn correlation potential functionals with memory, demonstrating one such model that yields reasonable results for trajectories outside the training set.
Correlations and measures of entanglement in ground-state wave functions of relativistic quantum field theories are spatially localized over length scales set by the mass of the field. We utilize this localization to design digital quantum circuits for preparing the ground states of lattice scalar quantum field theories. Controlled rotations that are exponentially localized in their position-space extent are found to provide exponentially convergent wave function fidelity. These angles scale with the correlation between sites and the classical two-point correlation function, as opposed to the more localized mutual information or the hyperlocalized (1×1)-site negativity. We anticipate that further investigations will uncover quantum circuit designs with controlled rotations dictated by the measures of entanglement. This work is expected to impact quantum simulations of systems of importance to nuclear physics, high-energy physics, and basic energy sciences research.
A search for collective effects inside jets produced in proton-proton collisions is performed via correlation measurements of charged particles using the CMS detector at the CERN LHC. The analysis uses data collected at a center-of-mass energy of s = 13 TeV , corresponding to an integrated luminosity of 138 fb − 1 . Jets are reconstructed with the anti- k T algorithm with a distance parameter of 0.8 and are required to have transverse momentum greater than 550 GeV and pseudorapidity | η jet | < 1.6 . Two-particle correlations among the charged particles within the jets are studied as functions of the particles’ azimuthal angle and pseudorapidity separations ( Δ ϕ * and Δ η * ) in a jet coordinate basis, where particles’ η * , ϕ * are defined relative to the direction of the jet. The correlation functions are studied in classes of in-jet charged-particle multiplicity up to N ch j ≈ 100 . Fourier harmonics are extracted from long-range azimuthal correlation functions to characterize azimuthal anisotropy for | Δ η * | > 2 . For low- N ch j jets, the long-range elliptic anisotropic harmonic, v 2 * , is observed to decrease with N ch j . This trend is well described by Monte Carlo event generators. However, a rising trend for v 2 * emerges at N ch j ≳ 80 , hinting at a possible onset of collective behavior, which is not reproduced by the models tested. This observation yields new insights into the dynamics of jet evolution in the vacuum. © 2024 CERN, for the CMS Collaboration 2024 CERN
We present a real-space method for computing the random phase approximation (RPA) correlation energy within Kohn–Sham density functional theory, leveraging the low-rank nature of the frequency-dependent density response operator. In particular, we employ a cubic-scaling formalism based on density functional perturbation theory that circumvents the calculation of the response function matrix, instead relying on the ability to compute its product with a vector through the solution of the associated Sternheimer linear systems. We develop a large-scale parallel implementation of this formalism using the subspace iteration method in conjunction with the spectral quadrature method while employing the Kronecker product-based method for the application of the Coulomb operator and the conjugate orthogonal conjugate gradient method for the solution of the linear systems. We demonstrate convergence with respect to key parameters and verify the method’s accuracy by comparing with plane-wave results. We show that the framework achieves good strong scaling to many thousands of processors, reducing the time to solution for a lithium hydride system with 128 electrons to around 150 s on 4608 processors.