SEARCH · Search NASA
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.
Single-electron acceleration in vacuum by exponentially chirped plane-wave and tightly focused Gaussian laser pulses is investigated. Numerical calculations and particle-in-cell simulations, employing the chirped fields, demonstrate that GeV net energy gain, over sub-millimeter excursion distances, is possible from PW-class femtosecond laser pulses.
We studied the LaFeSi phase using density functional theory based full-potential linearized augmented plane wave (FP-LAPW) method. Specifically we examined the effect of p-element substitution on the stability and magnetic properties focusing on pathways to induce Fe magnetic moments. We demonstrate that either partial or complete substitutions at the Si 2 c-site by several p-block elements lead to non-zero Fe moments. Our theoretical study shows that partially substituted LaFeSi 1-y Ga y and LaFeSi 1-y Al y are more thermodynamically stable than the fully substituted LaFeGa and LaFeAl. We also found that the P substitutions, either partial or complete, have the most negative formation energies, however the compounds containing phosphorus are non-magnetic. Our work highlights a pathway toward manipulations of the otherwise quenched Fe-moments in LaFeSi-related compounds and help in advancing the potential magnetic functionalities of LaFeX compounds.
Significance Density functionals with Hartree–Fock exchange have been widely used for a wide range of chemical applications, but the nonlocal character of exchange makes long-range exchange computationally expensive for solid-state calculations with periodic boundary conditions, and full exchange is nonphysical for condensed-phase systems. Here, we present a screened-exchange (SX) density functional, M06-SX, that is especially designed to have good accuracy for both solid-state physics and chemical applications with less computational cost than full Hartree–Fock exchange. The M06-SX functional gives accuracy comparable to functionals with full Hartree–Fock exchange for predicting chemical properties of molecules, while also being practical with good accuracy for plane wave calculations on band gaps and lattice constants of solids.
We embed the Penrose limit into the Weyl classical double copy. Thereby, we provide a lift of the double copy properties of plane wave spacetimes into black hole geometries and we open a novel avenue towards taking the classical double copy beyond statements about algebraically special backgrounds. In particular, the Penrose limit, viewed as the leading order Fermi coordinate expansion around a null geodesic, complements approaches leveraging asymptotic flatness such as the asymptotic Weyl double copy. Along the way, we show how our embedding of the Penrose limit within the Weyl double copy naturally fixes the functional ambiguity in the double copy for Petrov type N spacetimes. Here, we also highlight the utility of a spinorial approach to the Penrose limit. In particular, we use this spinorial approach to derive a simple analytical expression for arbitrary Penrose limits of four-dimensional, vacuum type D spacetimes.
The second-order optical susceptibility of semiconductors $\chi^{(2)}_{ijk}(-2\omega;\omega,\omega)$ finds application in metrology, spectroscopy, telecommunications, material characterization, and quantum information. Pioneering calculations of $\chi^{(2)}_{ijk}(-2\omega;\omega,\omega)$ utilized non-orthogonal Gaussian orbitals centered at atoms. That formulation transitioned into plane-wave-based algorithms as time went by. As of late, nevertheless, multiple tools for calculating optical susceptibilities have recast the problem using Wannier (i.e. localized) orbitals, making a comeback onto frameworks based on localized basis sets. Here, in this work, we present an approach for calculating $\chi^{(2)}_{ijk}(-2\omega;\omega,\omega)$ reliant on numerical pseudo-atomic orbitals (PAOs) within perturbation theory in the velocity gauge. Its salient feature is a calculation of ‘Slater–Koster-like’ two-center integrals of the momentum operator in between PAOs identified by symmetry. The approach was successfully tested on paradigmatic cubic silicon carbide (3C-SiC) and gallium arsenide, for which linear responses are contributed as well.
Highly coherent and powerful light sources capable of generating soft x-ray frequency combs are essential for high precision measurements and rigorous tests of fundamental physics. In this work, we derive the analytical conditions required for the emission of coherent radiation from an electron beam colliding with a laser pulse, modeled as a plane wave. These conditions are applied in a series of numerical simulations, where we show that a soft x-ray frequency comb can be produced if the electrons are regularly-spaced and sufficiently monoenergetic. High quality beams of this kind may be produced in the near future from laser-plasma interactions or linear accelerators. Furthermore, we highlight the advantageous role of employing few-cycle laser pulses in relaxing the stringent monoenergeticity requirements for coherent emission. The conditions derived here can also be used to optimize coherent emission in other frequency ranges, such as the terahertz domain.
Abstract Isolated attosecond pulses (IAPs) in the soft x-ray (SXR) region are highly desirable for attosecond time-resolved experiments. Here we identify a transient phase matching gating method for the generation of such IAPs with mid-infrared (MIR) lasers. This gating method works when a loosely focused, long-duration MIR Gaussian driving beam is slightly reshaped during its propagation in the gas medium under the ‘critical’ ionization condition. Quantitatively, the calculated coherence length of high harmonic is used to analyze the mechanism of the gating method, by using one-dimensional plane-wave beams and by comparing 2000 nm and 800 nm lasers. The robustness of the generation method is checked by varying laser’s carrier-envelope-phase and gas pressure. This gating method provides with an alternative route to efficiently produce tabletop ultrashort attosecond SXR light sources with the emerging MIR lasers.
In this article we investigate novel signatures of radiation reaction via the angular deflection of an electron beam colliding at 90 degrees with an intense laser pulse. Due to the radiation reaction effect, the electrons can be deflected towards the beam axis for plane wave backgrounds, which is not possible in the absence of radiation reaction effects. The magnitude and size of the deflection angle can be controlled by tailoring the laser pulse shapes. The effect is first derived analytically using the Landau–Lifshitz equation, which allows to determine the important scaling behavior with laser intensity and particle energy. We then move on to full scale 3D Monte Carlo simulations to verify the effect is observable with present day laser technology. We investigate the opportunities for an indirect observation of laser depletion in such side scattering scenarios.
Graph neural networks (GNNs) have been employed in materials research to predict physical and functional properties, and have achieved superior performance in several application domains over prior machine learning approaches. Recent studies incorporate features of increasing complexity such as Gaussian radial functions, plane wave functions, and angular terms to augment the neural network models, with the expectation that these features are critical for achieving a high performance. Here, we propose a GNN that adopts edge convolution where hidden edge features evolve during training and extensive attention mechanisms, and operates on simple graphs with atoms as nodes and distances between them as edges. As a result, the same model can be used for very different tasks as no other domain-specific features are used. With a model that uses no feature engineering, we achieve performance comparable with state-of-the-art models with elaborate features for formation energy and band gap prediction with standard benchmarks; we achieve even better performance when the dataset size increases. Although some domain-specific datasets still require hand-crafted features to achieve state-of-the-art results, our selected architecture choices greatly reduce the need for elaborate feature engineering and still maintain predictive power in comparison.
The quality of recent Short-Range Correlations/Color-Transparency (SRC/CT) Collaboration 𝐽/𝜓 photoproduction data off a 4 He target from Hall D at Jefferson Laboratory, combined with the feasibility of measuring the reaction close to the free-nucleon energy threshold, opens the door to using incoherent 𝐽/𝜓 photoproduction to explore a variety of interesting physics aspects. An example is the estimation of the 𝐽/𝜓 𝑝 scattering length |𝛼 𝐽/𝜓 𝑝 | on the bound proton, obtained using the vector meson dominance model. The new experiment E12–25–002, conducted by the SRC/CT Collaboration, was recently approved in Hall D with the GlueX spectrometer, which will significantly enhance the number of reconstructed 𝐽/𝜓 mesons. Using a plane-wave theoretical model to generate quasidata, we estimated that the experiment could achieve a result of |𝛼 𝐽/𝜓 𝑝 | = 3.08 ± 0.45 mfm. This value can be compared to an extraction based on the free proton data from the GlueX Collaboration, |𝛼 𝐽/𝜓 𝑝 | = 3.08 ± 0.55 mfm. A comparison between the two would provide an opportunity to evaluate the effects of medium modification in light nuclei.
Preconditioning is the most widely used and effective way for treating ill-conditioned linear systems in the context of classical iterative linear system solvers. We introduce a quantum primitive called fast inversion, which can be used as a preconditioner for solving quantum linear systems. The key idea of fast inversion is to directly block encode a matrix inverse through a quantum circuit implementing the inversion of eigenvalues via classical arithmetics. We demonstrate the application of preconditioned linear system solvers for computing single-particle Green's functions of quantum many-body systems, which are widely used in quantum physics, chemistry, and materials science. We analyze the complexities in three scenarios: the Hubbard model, the quantum many-body Hamiltonian in the plane-wave-dual basis, and the Schwinger model. We also provide a method for performing Green's function calculation in second quantization within a fixed-particle manifold and note that this approach may be valuable for simulation more broadly. Aside from solving linear systems, fast inversion also allows us to develop fast algorithms for computing matrix functions, such as the efficient preparation of Gibbs states. Furthermore, we introduce two efficient approaches for such a task, based on the contour-integral formulation and the inverse transform, respectively.
Here we present a technique to assess the focal volume of petawatt-class lasers at full power. Our approach exploits quantitative measurement of the angular distribution of electrons born in the focus via ionization of rarefied gas, which are accelerated forward and ejected ponderomotively by the field. We show that a bivariate (θ,Φ) angular distribution, which was obtained with image plates, not only enables the peak intensity to be extracted, but also reflects nonideality of the focal-spot intensity distribution. In our prototype demonstration at intensities of a few ×10 19 to a few ×10 20 W/cm 2 , an f/10 optic produced a focal spot in the paraxial regime. This allows a plane-wave parametrization of the peak intensity given by tanθ c =2/a 0 (a 0 being the normalized vector potential and θ c the minimum ejection angle) to be compared with our measurements. Qualitative agreement was found using an a 0 inferred from the pulse energy, pulse duration, and focal spot distribution with a modified parametrization, tanθ c =2η/a 0 ($η=2.02^{+0.26}_{–0.22}$). This highlights the need for (i) better understanding of intensity degradation due to focal-spot distortions and (ii) more robust modeling of the ejection dynamics. Using single-shot detection of electrons, we showed that while there is significant shot-to-shot variation in the number of electrons ejected at a given angular position, the average distribution scales with the pulse energy in a way that is consistent with that seen with the image plates. Finally, we note that the asymptotic behavior as θ→0° limits the usability of angular measurement. For 800 nm, this limit is at an intensity ~10 21 W/cm 2 .
We present ab initio calculations of electron capture and release coefficients as a function of temperature for Tl dopants in the widely used scintillator NaI, which is a soft ionic crystal. The modelled capture and release events occur by transitions mediated by multiple phonon absorption and emission between states with significantly different local geometries around the trapping site. We demonstrate that such transitions are not well-described by the normal harmonic approximation to the nuclear dynamics. We avoid the harmonic approximation by numerically solving the vibrational Schrödinger equation for the motion along/of a single effective phonon coordinate. The localized trapped state is not correctly described by semilocal density functionals, so we employ a hybrid functional tuned to reproduce the band gap and lattice constant of the host material. Our calculations, which combine plane-wave and cluster electronic structure methods in a novel way, are in reasonable agreement with available experimental data for detrapping and predict an unusual temperature dependence for trapping rates.
In this work, first-principles methods were used to understand magnetic flux trapping at vacancies, dislocations, and grain boundaries in high-purity superconducting niobium. Full-potential linear augmented plane-wave methods were applied in progressively greater complexity, starting at simple vacancies and extending to screw dislocations and tilt grain boundaries to analyze the effects of magnetic field on the superconducting state surrounding these defects. Density-functional theory calculations identified changes in electronic structure at the dislocation core and different types of symmetric tilt grain boundaries relative to bulk niobium. Electron redistribution enhanced nonparamagnetic effects and thus perturb superconductivity, resulting in local conditions suitable for premature flux penetration and subsequently flux pinning. Since the coherence length of superconducting niobium at 0 K is significantly larger than the lattice parameter, the effects of line and planar defects in niobium are predicted to be stronger for defect clusters than single defects in isolation, which is consistent with recent experimental observations. Controlling accumulation or depletion of charge at the defects, e.g., by segregation of an impurity atom, can mitigate these tendencies thus increasing the quality of superconducting niobium.
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.
Here, we present a computationally efficient approach to perform systematically convergent real-space all-electron Kohn-Sham density functional theory calculations for solids using an enriched finite element (FE) basis. The enriched FE basis is constructed by augmenting the classical FE basis with atom-centered numerical basis functions, comprising of atomic solutions to the Kohn-Sham problem. Notably, to improve the conditioning, we orthogonalize the enrichment functions with respect to the classical FE basis, without sacrificing the locality of the resultant basis. In addition to improved conditioning, this orthogonalization procedure also renders the overlap matrix block diagonal, greatly simplifying its inversion. Subsequently, we use a Chebyshev polynomial based filtering technique to efficiently compute the occupied eigenspace in each self-consistent field iteration. We demonstrate the accuracy and efficiency of the proposed approach on periodic unit cells and supercells. The benchmark studies show a staggering 130× speedup of the orthogonalized enriched FE basis over the classical FE basis. We also present a comparison of the orthogonalized enriched FE basis with the linearized augmented plane-wave + local orbitals basis, both in terms of accuracy and efficiency. Notably, we demonstrate that the orthogonalized enriched FE basis can handle large system sizes of ~10 000 electrons. Finally, we observe good parallel scalability of our implementation with 92% efficiency at 22× speedup for a system with 620 electrons
Here we present a two-step method specifically tailored for band structure calculation of the small-angle moiré-pattern materials which contain tens of thousands of atoms in a unit cell. In the first step, the self-consistent field calculation for the ground state is performed with the O(N) Krylov subspace method implemented in openmx. Second, the crystal momentum-dependent Bloch Hamiltonian and overlap matrix are constructed from the results obtained in the first step and only a small number of eigenvalues near the Fermi energy are solved with shift-invert and Lanczos techniques. By systematically tuning two key parameters, the cutoff radius for electron hopping interaction and the dimension of the Krylov subspace, we obtained the band structures for both rigid and corrugated twisted bilayer graphene structures down to the first magic angle (θ = 1.08°) with high enough accuracy at affordable costs. The band structures are in good agreement with those from tight-binding models, continuum models, plane-wave pseudopotential based ab initio calculations, and experimental observations. This method is also shown to be efficient in twisted double-bilayer graphene and bilayer WSe 2 . We think this two-step method can play a crucial role in other twisted two-dimensional materials, especially those with much more complex band structure and where the effective model is hard to construct.
The implementation of optimal strictly localized atomic orbitals basis for plutonium (Pu) using norm-conserving pseudopotential density-functional theory (DFT) is presented. The basis set was applied to the α, β, γ, δ, δ', and ε phases of Pu, δ-Pu surface, and δ-PuGa alloys. The computed properties of the Pu phases and δ-Pu surface were in good agreement with both available experimental data and prior DFT calculations based on plane-wave methodologies. Results for the δ-PuGa alloys were also in good agreement with experimental data. Finally, the reliability of the basis set was further demonstrated by using ab initio molecular dynamics to model the diffusion coefficient and activation barrier for atomic diffusion in a δ-PuGa alloy.