Temperature-dependent phonon linewidths and shifts in bismuth from inelastic neutron scattering
Study of temperature-dependent phonon linewidths and shifts in bismuth from inelastic neutron scattering.
36 MATERIALS SCIENCE↗
SEARCH · Search NASA
Results for “first-principles”
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.
Predicting Open Quantum Dynamics with Data-Informed Quantum-Classical Dynamics
We introduce a data-informed quantum-classical dynamics (DIQCD) approach for predicting the evolution of an open quantum system. The equation of motion in DIQCD is a Lindblad equation with a flexible, time-dependent Hamiltonian that can be optimized to fit sparse and noisy data from local observations of an extensive open quantum system. We demonstrate the accuracy and efficiency of DIQCD for both experimental and simulated quantum devices. We show that DIQCD can predict entanglement dynamics of ultracold molecules (calcium fluoride) in optical tweezer arrays. DIQCD also successfully predicts carrier mobility in organic semiconductors (rubrene) with accuracy comparable to nearly exact numerical methods.
Lindblad equation↗
Kaon gluon parton distribution and momentum fraction from 2+1+1 lattice QCD with high statistics
We present a high-statistics lattice-QCD determination of the kaon gluon parton distribution function and gluon momentum fraction. We use clover valence fermion action to take 1,296,640 kaon-correlator measurements on a highly improved staggered quark ensemble with 𝑎 ≈ 0.12 fm and 310-MeV pion mass generated by the MILC Collaboration. A detailed investigation into the impact of gauge-link smearing on the gluonic matrix elements indicates that five steps of hypercubic smearing offer an effective balance between signal quality and preservation of long-distance physics. We report a nonperturbatively renormalized kaon gluon momentum fraction of Math output error at 𝜇 = 2 GeV in the Math output error scheme. Using reduced pseudo-Ioffe-time-distribution matrix elements and pseudo-parton-distribution-function (PDF) matching, we extract the kaon gluon PDF and compare with the prediction from the Dyson-Schwinger equation and with the pion PDF obtained from the same ensemble.
First-principles calculations↗
Magnetic semiconductors from functionalized Cr 2 C MXenes
Here, we report an ab initio investigation of functionalized and 3𝑑-electron-doped Cr 2 C MXenes. Upon functionalization, the Cr 2 C becomes chemically, dynamically, and mechanically stable, and it exhibits magnetic semiconducting behavior. Cr 2 CF 2 stands out as a wide band gap semiconductor, possessing superexchange interaction mediated by F atoms within the layer; however, the applied strain transforms it from an indirect to a direct band gap semiconductor. Strong spin-phonon coupling found in Cr 2 CH 2 is supported by the distorted Cr spin density due to the hydrogen environment. Two magnon branches, associated with two sublattice spins, are found in the ferromagnetic Cr 2 CO 2 and antiferromagnetic Cr 2 CF 2 . Depending on the types of 3𝑑-electron dopants and functionalization, Cr 2 C MXenes (except for Cr 2 CO 2 ) change from an indirect band gap magnetic semiconductor to different states of electronic and magnetic matter, including an exotic direct band gap magnetic bipolar semiconductor. In addition, we reveal a band inversion between the two highest valence bands in the Fe-doped Cr 2 CCl 2 .
2D materials↗
Development of an ab initio learned model of electron deposition range in deuterium-tritium plasmas through time-dependent density functional theory calculations and machine learning
Accurate hydrodynamic modeling for laser-direct-drive (LDD) inertial-confinement-fusion (ICF) relies on precise calculations of the electron thermal conduction in all target materials. The nonlocal stopping range of electrons in ICF plasmas directly influences thermal conduction; yet, no first principles model exists for the electron mean free path in the conduction-zone regime. This work utilized time-dependent stochastic density-functional theory (TD-sDFT) to calculate the electron stopping power in deuterium-tritium (DT) plasmas at (ρ, T) conditions relevant to the conduction zone and the compressed shell in ICF. Using a combination of our TD-sDFT data and already established analytical models, we developed and trained an artificial neural network to create a global model for the nonlocal electron deposition range, λ E . We compared our machine-learning (ML) based model for λ E to the currently-used modified-Lee-More model in LDD radiation-hydrodynamic codes, such as lilac, and saw an overall decrease in the deposition range. To understand the effects of λ E on LDD ICF implosion dynamics, we implemented the ML-based model into lilac; specifically, we looked at designs consistent with a current experiment on the OMEGA laser and for a newly designed LDD-ICF target for the future OMEGA-Next facility. In both cases, we saw an overall drop in predicted ablation pressure, peak areal density, and neutron yield due to the reduced thermal conduction (smaller λ E ) in DT plasmas. Comparisons with the experiment on OMEGA are also made.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Electronic coherence evolution at the nearly commensurate-incommensurate CDW boundary of 1T−TaS 2
Transition-metal dichalcogenides host a variety of charge-density-wave phases that couple lattice, charge, and correlation effects. In 1𝑇−TaS 2 , the commensurate and nearly commensurate states are well characterized, yet the transition near 350 K into the incommensurate phase has lacked direct momentum-resolved insight. In this work, we use temperature-dependent angle-resolved photoemission spectroscopy to track the electronic structure across this transition. We observe a suppression of quasiparticle spectral weight at the Brillouin-zone center, coincident with the transport anomaly, but without clear evidence of a full band-gap opening. The transition appears to involve momentum-dependent redistribution of spectral weight, consistent with a loss of coherence that reshapes the Fermi surface while leaving conduction dispersions largely intact. These results suggest that the nearly commensurate–incommensurate transition may not align with a conventional metal-insulator transition picture, but rather as an electronic reconstruction driven by loss of coherence. Our work provides new microscopic insight into the resistivity anomaly near room temperature and may guide design principles for collective electronic switching in transition-metal dichalcogenides.
36 MATERIALS SCIENCE↗
Comprehensive analytical model of the dynamic 𝑍 pinch
In this work we present an analytical 1D axisymmetric model describing the evolution of the dynamic 𝑍 pinch. This model is capable of predicting the trajectories of the imploding sheath's magnetic piston and preceding shock front, along with the velocity, pressure, density, and magnetic field profiles, for any time-dependent current, spatially varying initial density profile, and weak initial axial field. The implosion is divided into stages, with each stage described by a set of coupled ordinary differential equations derived from the ideal MHD equations. Comparisons with experimental data from the COBRA pulsed-power facility are quite promising and imply this model could prove useful in designing and analyzing future pulsed-power experiments.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
𝐴𝑏 initio density-matrix approach to exciton coherence: Phonon scattering, Coulomb interactions, and radiative recombination
Relaxation processes following light excitation in semiconductors are key in materials-based quantum technology applications. These processes are broadly studied in atomically thin transition-metal dichalcogenides, quasi-two-dimensional excitonic semiconductors in which atomistic design allows for tunable excited-state properties, such as relaxation lifetimes and photoinduced coherence. In this work, we present a density-matrix-based approach to compute exciton relaxation within a many-body ab initio perspective. We expand our previously developed Lindblad density-matrix formalism to capture multichannel electron-hole pair relaxation processes, including phonon and Coulomb scattering as well as radiative recombination, and we study their effect on the time-resolved excited-state propagation. Using monolayer MoSe 2 as a prototypical example, we examine many-body effects on the time-dependent dynamics of photoactive excitations, exploring how the electron-hole pair interactions are reflected in variations of the excitation energy, spectral signature, and state coherence. In conclusion, our method supplies a detailed understanding of exciton relaxation mechanisms in realistic materials, offering a previously unexplored pathway to study excited-state dynamics in semiconductors from first principles.
Band structure methods↗
Measuring the Debye energy in superconductors via two-electron photoemission spectroscopy
In this work, we demonstrate theoretically that double angle resolved photoemission spectroscopy (2eARPES) can directly probe the existence of Cooper pairs away from the Fermi surface, and can thus provide insight into the characteristic energy scale around the Fermi surface, the Debye energy, in which electrons are bound into Cooper pairs. To this end, we compute the photoelectron counting rate 𝑃 (2) in two different types of unconventional superconductors, a 𝑑 𝑥 2 −𝑦 2 -wave superconductor, and a topological superconductor with a broken time-reversal symmetry. We show that 𝑃 (2) provides insight into the relative strength of intra- and interband pairing in multiband systems, as well as into the spin polarization of the bands.
75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
Anisotropy, frustration, and saddle point in the twisted kagome compound ErPdPb
The kagome lattice, with its inherent geometric frustration, provides a rich platform for exploring intriguing magnetic phenomena and topological electronic structures. In reduced-symmetry structures, such as twisted kagome systems involving rare earth elements, additional anisotropy can arise, enabling intriguing properties including spin-ice states, magnetocaloric effects, noncollinear magnetic ordering, and the anomalous Hall effect. Here, we report the synthesis of single crystals of ErPdPb, which features a twisted kagome lattice net of Er atoms within the hexagonal ZrNiAl-type structure, and we investigate its magnetic, electronic, and thermal properties. The material exhibits a highly anisotropic, correlated magnetic state below 2.7 K, as evidenced by magnetic, transport, and heat capacity measurements. Density functional theory (DFT) calculations indicate strong easy-axis anisotropy, consistent with experiment and crystal-field expectations, as well as quasi-one-dimensional bands and a spin-split saddle point at the zone center. Here, the coexistence of competing magnetic interactions along different crystallographic directions suggests an inherent degree of frustration. ErPdPb thus provides a promising platform for exploring the interplay of frustration, anisotropy, and electronic structure in a twisted kagome lattice.
Antiferromagnets↗
Single-crystalline orthorhombic GdAlGe as a rare-earth magnetic Dirac nodal-line metal
Crystal engineering is a method for discovering new quantum materials and phases, which may be achieved using external pressure or strain. Chemical pressure is unique in that it generates internal pressure perpetually to the lattice. As an example, GdAlSi from the rare-earth (𝑅) 𝑅Al𝑋 (𝑋=Si or Ge) family of Weyl semimetals is considered. Replacing Si with the larger isovalent element Ge creates sufficiently large chemical pressure to induce a structural transition from the tetragonal structure of GdAlSi, compatible with a Weyl semimetallic state, to an orthorhombic phase in GdAlGe, resulting in an inversion-symmetry-protected nodal-line metal. We find that GdAlGe hosts an antiferromagnetic ground state with two successive orderings, at 𝑇 N 1=35K and 𝑇 N 2=30K. In-plane isothermal magnetization shows a magnetic field induced metamagnetic transition at 6.2 T for 2 K. Furthermore, electron-hole compensation gives rise to a large magnetoresistance of ∼100% at 2 K and 14 T. Angle-resolved photoemission spectroscopy measurements and density functional theory calculations reveal a Dirac-like linear band dispersion over an exceptionally large energy range of ∼1.5eV with a high Fermi velocity of ∼10 6 m/s, a rare feature not observed in any magnetic topological materials.
75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
Ferroelectric phase transition in group-IV monochalcogenides from an equivariant machine learned force field
Group-IV monochalcogenides are a class of layered ferroelectric semiconductors that have demonstrated spontaneous intrinsic polarization above room temperature. Here, in this study, we use the multi-atomic cluster expansion (MACE) machine learning architecture to train and test a force field capable of modeling the structural properties and second-order ferroelectric-to-paraelectric phase transition in a Group-IV monochalcogenide, GeSe. The model captures the double-well potential energy surface associated with the onset of macroscopic polarization in bulk GeSe within 12.5 meV/atom, as well as near-equilibrium properties like the phonon dispersion. The development of this quantitatively accurate force field enables long-time molecular dynamics simulations, which predict the critical temperature of the ferroelectric-to-paraelectric phase transition in bulk GeSe to be T c = 600 K. This study demonstrates the capabilities of equivariant force-fields to accurately describe phenomena associated with structural symmetry breaking.
ferroelectricity↗
Engineering 2D Square Lattice Hubbard Models in 90° Twisted GeX/SnX (X =S, Se) Moiré Superlattices
Because of the large-period superlattices emerging in moiré two-dimensional (2D) materials, electronic states in such systems exhibit low energy flat bands that can be used to simulate strongly correlated physics in a highly tunable setup. While many investigations have thus far focused on moiré flat bands and emergent correlated electron physics in triangular, honeycomb, and quasi-one-dimensional lattices, tunable moiré realizations of square lattices subject to strong correlations remain elusive. Here, in this work, we propose a feasible scheme to construct moiré square lattice systems by twisting two or more layers of 2D materials in a rectangular lattice by 90°. We demonstrate the concept with twisted GeX/SnX (X =S, Se) moiré superlattices and calculate their electronic structures from first principles. We show that the lowest conduction flat band in these systems can be described by a square lattice Hubbard model with parameters which can be controlled by varying the choice of host materials, number of layers, and external electric fields. In particular, twisted double bilayer GeSe realizes a square lattice Hubbard model with strong frustration due to the next-nearest-neighbor hopping that could host unconventional superconductivity, in close analogy to the Hubbard model for copper-oxygen planes of cuprate high-temperature superconductors. The presented scheme uses 90° twisted 2D materials with rectangular unit cells as a promising platform for realizing the physical phenomena of square lattice Hubbard models, establishing a new route for studying its rich phase diagram of magnetism, charge order, and unconventional superconductivity in a highly tunable setting.
2-dimensional systems↗
Experimental confirmation of the magnetic ordering transition induced by an electronic structure change in the metallic triangular antiferromagnet Co 1/3 TaS 2
We report angle-resolved photoemission spectroscopy (ARPES) studies combined with DFT+DMFT calculations to confirm that the magnetic ordering vector transition from Q=(1/2,0,0) to Q=(1/3,0,0) in the metallic triangular antiferromagnets Co 1/3±𝜀 TaS 2 (𝜀 ≈ 0.007) is induced by the electronic structure change in the system. The ARPES-measured Fermi surface (FS) maps of Co 0.325 TaS 2 show two hexagonal and one circular holelike FSs around Γ, which matches well with the triple- Q state by taking into account the contribution of nesting vectors occurring between Co 3𝑑 and Ta 5𝑑 orbitals. In the case of Co 0.340 TaS 2 , a new electron pocket around K appears, and the FS geometry changes as a result of the correlation effect of Co 4 S 18 tripods forming in the system. The magnetic susceptibility calculations based on the charge-self-consistent DFT+DMFT band structures and the random phase approximation indicate that the most stable magnetic ordering vector (1/2,0,0) split into (1/6,0,0) and (1/2,0,0), which is consistent with the magnetic phase transition around 𝑥=1/3 in Co 𝑥 TaS 2 .
Fermi surface↗
Impact of absorption due to zero-field splitting on loss in dielectrics: A case study in sapphire
The coherence times of superconducting qubits are limited by loss mechanisms, whose microscopic origins have remained elusive. We propose a mechanism caused by transitions between zero-field-split states of paramagnetic impurities or defects. We derive the absorption cross section for a magnetic dipole transition and apply it to calculate the bulk loss tangent. For Cr, Fe, and V impurities in sapphire, we find loss tangents at 4.5 GHz in the range of 10 −9 –10 −8 , comparable to the loss measured in experiments. This value suggests that magnetic loss may be a limiting factor in the coherence times of superconducting qubits.
36 MATERIALS SCIENCE↗
Fokker-Planck Equation Governing the Distribution of Walkers in Auxiliary-Field Quantum Monte Carlo
Auxiliary-field quantum Monte Carlo (AFQMC) is typically formulated as an open-ended random walk in an overcomplete space of Slater determinants, implemented through a Langevin equation. However, the explicit form of the underlying Fokker-Planck equation governing the walker population distribution has remained unknown. Here, in this Letter, we derive the Fokker-Planck equation for AFQMC and propose a novel numerical scheme to solve it. The solution of the Fokker-Planck equation reveals the wave function actually sampled by the AFQMC algorithm. Interestingly, we find that even when the exact ground state is used as a guiding wave function in constrained path AFQMC, contrary to the common assumption, the wave function sampled by AFQMC is not exact. Beyond clarifying several fundamental aspects of AFQMC, the availability of a Fokker-Planck equation formulation opens new avenues for systematically improving its accuracy, which we outline in this Letter.
Monte Carlo methods↗
Stability, electronic quantum states, and magnetic interactions of Er 3+ ions in Ga 2 O 3
Here, we report an ab initio study of phase stability, defect formation, electronic structure, and multiple magnetic, Dzyaloshinskii-Moriya, optical, hyperfine, and crystal field interactions in erbium (Er)-doped wide band gap 𝛼- and 𝛽-gallium oxides (Ga 2 O 3 ), critically important to make a foundation for both optoelectronic and quantum information applications. The chemical, structural, mechanical, and dynamical stabilities of the pristine phases are confirmed from respective negative formation energies, negative cohesive energies, favorable elastic constants, and positive phonon frequencies. The phonon dispersions indicate that the Ga-O bonds are uniform in the 𝛼-phase, while they vary in the 𝛽-phase due to the anisotropic polyhedral movement. The defect formation energy analysis confirms that both Er-doped 𝛼- and 𝛽−Ga 2 O 3 prefer Er 3+ (neutral) state. The underestimated band gaps of the pristine phases from standard density functional theory (DFT) calculations as compared to experimental values are corrected by employing the hybrid functional calculations, resulting in the indirect band gaps of 5.21 eV in 𝛼−Ga 2 O 3 and 4.94 eV in 𝛽−Ga 2 O 3 . The site preference energy analysis indicates partial occupation of Er in the octahedral site of Ga. The anisotropic nature of hyperfine tensor coefficients of Er is similar in both phases, which may be due to the occupation of Er in the same octahedral Ga site. On the other hand, the calculated magnetic exchange interaction between two Er dopants is negative for 𝛼 and positive for 𝛽, indicating an antiferromagnetic (AFM) ground state in the former and a ferromagnetic (FM) ground state in the latter. Large values of Dzyaloshinskii-Moriya interactions (DMIs) are obtained along the 𝑥 direction in the 𝛼 and along the 𝑦 direction in the 𝛽. The large DMI may support exotic magnetic textures, a promising direction for spintronic applications. The analysis of dielectric constants and refractive indices of both pristine and Er-doped phases shows a good agreement with available experimental values. The calculated optical anisotropy is slightly higher in 𝛽 than those in 𝛼, which is due to the involvement of lower symmetry in 𝛽. The crystal field coefficients (CFCs) calculated from DFT are used to analyze 4𝑓 multiplets and 4𝑓 −4𝑓 transitions. Thus calculated lowest energy level of the first excited state to the lowest energy level of the ground state is about 1.53 µm, which is in a good agreement with available experiments, and it falls within the quantum telecommunication wavelength range.
3-dimensional systems↗
Electronic and optical properties of transition metal impurities in GaN and AlN
The optical properties of Cr(III) impurities are well studied and exploited in 𝛼−Al 2 O 3 , which enables ruby lasers. Photoluminescence spectra resembling the famous 1.78 eV emission of Cr in 𝛼−Al 2 O 3 have been observed in other wide-band-gap host materials, in particular AlN. However, it is unclear whether Cr would produce the same emission in the nitrides as in 𝛼−Al 2 O 3 , especially given the different local environment (tetrahedral vs octahedral coordination). We investigate the electronic and optical properties of transition metal substitutional impurities (Cr, Co, Ni, and Cu) in wurtzite GaN and AlN from first principles, utilizing hybrid density functional theory. We calculate the optical transition energies using a constrained occupation approach that accounts for the multideterminant nature of the excited states. Our results indicate that substitutional Cr(III) is not responsible for the emission near 1.78 eV observed in the nitrides. Among the examined impurities, only substitutional Co(II) produces emission near 1.78 eV. Based on the favorable properties of Co(II) (low formation energy, small Huang-Rhys factor, and an 𝑆 = 3/2 ground state) we propose it as a promising single-photon emitter or spin qubit. We additionally assign the PL peak observed near 1.05 eV in GaN and AlN to the 2 𝐸→ 4 𝐴 2 transition of the Ni(III) substitution.
36 MATERIALS SCIENCE↗