Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “wave function”

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 559 records · Page 31
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Cohesion and excitations of diamond-structure silicon by quantum Monte Carlo: Benchmarks and control of systematic biases

We have carried out quantum Monte Carlo (QMC) calculations of silicon crystal focusing on the accuracy and systematic biases that affect the electronic structure characteristics. The results show that 64 and 216 atom supercells provide an excellent consistency for extrapolated energies per atom in the thermodynamic limit for ground, excited, and ionized states. We have calculated the ground state cohesion energy with both systematic and statistical errors below ≈ 0.05 eV. The ground state exhibits a fixed-node error of only 1.3(2)% of the correlation energy, suggesting an unusually high accuracy of the corresponding single-reference trial wave function. We obtain a very good agreement between optical and quasiparticle gaps that affirms the marginal impact of excitonic effects. Our most accurate results for band gaps differ from the experiments by about 0.2 eV. This difference is assigned to a combination of residual finite-size and fixed-node errors. In conclusion, we have estimated the crystal Fermi level referenced to vacuum that enabled us to calculate the edges of valence and conduction bands in agreement with experiments.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Fractional chiral hinge insulator

We propose and study a wave function describing an interacting three-dimensional fractional chiral hinge insulator (FCHI) constructed by Gutzwiller projection of two noninteracting second-order topological insulators with chiral hinge modes at half filling. We use large-scale variational Monte Carlo computations to characterize the model states via the entanglement entropy and charge-spin fluctuations. We show that the FCHI possesses fractional chiral hinge modes characterized by a central charge c = 1 and Luttinger parameter K = 1/2, like the edge modes of a Laughlin 1/2 state. The bulk and surface topology is characterized by the topological entanglement entropy (TEE) correction to the area law. While our computations indicate a vanishing bulk TEE, we show that the gapped surfaces host an unconventional two-dimensional topological phase. In a clear departure from the physics of a Laughlin 1/2 state, we find a TEE per surface compatible with (ln √2 )/2, half that of a Laughlin 1/2 state. This value cannot be obtained from topological quantum field theory for purely two-dimensional systems. For the sake of completeness, we also investigate the topological degeneracy.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Entanglement transitions from restricted Boltzmann machines

The search for novel entangled phases of matter has lead to the recent discovery of a new class of “en- tanglement transitions,” exemplified by random tensor networks and monitored quantum circuits. Most known examples can be understood as some classical ordering transitions in an underlying statistical mechanics model, where entanglement maps onto the free-energy cost of inserting a domain wall. In this paper we study the possibility of entanglement transitions driven by physics beyond such statistical mechanics mappings. Motivated by recent applications of neural-network-inspired variational Ansätze, we investigate under what conditions on the variational parameters these Ansätze can capture an entanglement transition. We study the entanglement scaling of short-range restricted Boltzmann machine (RBM) quantum states with random phases. For uncorrelated random phases, we analytically demonstrate the absence of an entanglement transition and reveal subtle finite-size effects in finite-size numerical simulations. Introducing phases with correlations decaying as 1/r α in real space, we observe three regions with a different scaling of entanglement entropy depending on the exponent α. We study the nature of the transition between these regions, finding numerical evidence for critical behavior. Furthermore, our work establishes the presence of long-range correlated phases in RBM-based wave functions as a required ingredient for entanglement transitions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Interaction effects in graphene in a weak magnetic field

A weak perpendicular magnetic field, B, breaks the chiral symmetry of each valley in the electron spectrum of graphene, preserving the overall chiral symmetry in the Brillouin zone. Here we explore the consequences of this symmetry breaking for the interaction effects in graphene. In particular, we demonstrate that the electron-electron interaction lifetime acquires an anomalous B dependence. Also, the ballistic zero-bias anomaly, δv(ω), where ω is the energy measured from the Fermi level, emerges at a weak B and has the form δv(B) ~ B 2 /ω 2 . Temperature dependence of the magnetic-field corrections to the thermodynamic characteristics of graphene is also anomalous. We discuss experimental manifestations of the effects predicted. The microscopic origin of the B-field sensitivity is an extra phase acquired by the electron wave function resulting from the chirality-induced pseudospin precession.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Trions in twisted bilayer graphene

The strong-coupling phase diagram of magic-angle twisted bilayer graphene (TBG) predicts a series of exact one-particle charge ±1 gapped excitations on top of the integer-filled ferromagnetic ground states. Finite-size exact diagonalization studies showed that these are the lowest charge ±1 excitations in the system (for 10 nm screening length), with the exception of charge +1 at filling -1 in the chiral limit. In the current paper we show that this “trion bound state,” a 3-particle, charge 1 excitation of the insulating ferromagnetic ground state of the projected Hamiltonian of TBG, is the lowest charge +1 overall excitation at ν =-1, and also for some large (≈20 nm) screening lengths at ν =-2 in the chiral limit and with very small binding energy. At other fillings, we show that trion bound states do exist, but only for momentum ranges that do not cover the entire moiré Brillouin zone. The trion bound states (at different momenta) exist for finite parameter range w 0 /w 1 but they all disappear in the continuum far below the realistic values of w 0 /w 1 = 0.8. Thus we find the conditions for the existence of the trion bound state, a good variational wave function for it, and investigate its behavior for different screening lengths, at all integer fillings, on both the electron and hole sides.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Effect of emitters on quantum state transfer in coupled cavity arrays

Over the last decade, conditions for perfect state transfer in quantum spin chains have been discovered and their experimental realizations addressed, as have their extensions to more complex geometries of coupled cavity-emitter arrays. In this paper we further consider such studies and situations in which quantum state transfer can occur with high fidelity, even when the cavity-cavity coupling rates and cavity-emitter interaction rates are comparable. This is accomplished through the development and use of a Monte Carlo approach to the inverse eigenvalue problem, which allows the determination of coupling rates which optimize quantum state transfer fidelity and subsequent time evolution of the polariton wave function through exact diagonalization of the resulting Jaynes-Cummings-Hubbard Hamiltonian. The effect of inhomogeneous emitter locations is also evaluated. Furthermore, our key results include the demonstration that our methodology can be used successfully to establish Hamiltonian parameters for high-fidelity state transfer in more general lattice geometries and excitation number sectors, and also a determination of the effects of fluctuations in those parameters about their optimal values.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2022

Twistronics of Janus transition metal dichalcogenide bilayers

Twisted multilayers of two-dimensional (2D) materials are an increasingly important platform for investigating quantum phases of matter, and in particular, strongly correlated electrons. The moiré pattern introduced by the relative twist between layers creates effective potentials of long wavelength, leading to electron localization. However, in contrast to the abundance of 2D materials, few twisted heterostructures have been studied until now. Here, in this study, we develop a first-principles continuum theory to study the electronic bands introduced by moire patterns of twisted Janus transition metal dichalcogenide (TMD) homo- and heterobilayers. The model includes lattice relaxation, stacking-dependent effective mass, and Rashba spin-orbit coupling. We then perform a high-throughput generation and characterization of DFT-extracted continuum models for more than a hundred possible combinations of materials and stackings. Our model predicts that the moiré physics and emergent symmetries depend on chemical composition, vertical layer orientation, and twist angle, so that the miniband wave functions can form triangular, honeycomb, and kagome networks. Rashba spin-orbit effects, peculiar of these systems, can dominate the moiré bandwidth at small angles. Our work enables detailed investigation of Janus twisted heterostructures, allowing the discovery and control of novel electronic phenomena.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Influence of local symmetry on lattice dynamics coupled to topological surface states

Herein we investigate coupled electron-lattice dynamics in the topological insulator Bi 2 Te 3 with time-resolved photoemission and time-resolved x-ray diffraction. It is well established that coherent phonons can be launched by optical excitation, but selection rules generally restrict these modes to zone-center wave vectors and Raman-active branches. We find that the topological surface state couples to additional modes, including a continuum of surface-projected bulk modes from both Raman and infrared branches, with possible contributions from surface-localized modes when they exist. Our calculations show that this surface vibrational spectrum occurs naturally as a consequence of the translational and inversion symmetries broken at the surface, without requiring the splitting-off of surface-localized phonon modes. The generality of this result suggests that coherent phonon spectra are useful by providing unique fingerprints for identifying surface states in more controversial materials. These effects may also expand the phase space for tailoring surface state wave functions via ultrafast optical excitation.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Realization of flat bands by lattice intercalation in kagome metals

Recently there has been intense interest in kagome metals, which are expected to host flat bands (FBs). However, the observed “FBs” are not flat over the whole two-dimensional Brillouin zone and overlap strongly with other bands. In fact, the FB does not truly exist in a default d-orbital kagome lattice, and the conditions for its existence in kagome metals are unknown. Here, based on tight-binding model analyses of the interplay between orbital and lattice symmetry, we establish such conditions. Additionally, we show that for a single d-orbital kagome lattice assuming large crystal field splitting (CFS), only the d z 2 orbital gives rise to a FB, while d xy , d x 2 -y 2 , d x⁢z , and d yz orbitals can only produce a FB with a rotated d-orbital basis so that they conform with the underlying kagome lattice symmetry. Most importantly, we demonstrate that both conditions of d-orbital rotation and large CFS can be ideally satisfied by intercalating the kagome lattice with a hexagonal sublattice without disrupting the destructive interference of FB wave function. Furthermore, we propose layered metalorganic frameworks as promising candidate kagome metals to realize FBs.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Importance of nonuniform Brillouin zone sampling for ab initio Bethe-Salpeter equation calculations of exciton binding energies in crystalline solids

Excitons are prevalent in semiconductors and insulators, and their binding energies are critical for optoelectronic applications. The state-of-the-art method for first-principles calculations of excitons in extended systems is the ab initio GW-Bethe-Salpeter equation (BSE) approach, which can require a fine sampling of reciprocal space to accurately resolve solid-state exciton properties. Here, in this study, we show, for a range of semiconductors and insulators, that the commonly employed approach of uniformly sampling the Brillouin zone can lead to underconverged exciton binding energies, as impractical grid sizes are required to achieve adequate convergence. We further show that nonuniform sampling of the Brillouin zone, focused on the region of reciprocal space where the exciton wave function resides, enables efficient rapid numerical convergence of exciton binding energies at a given level of theory. We propose a well-defined convergence procedure, which can be carried out at relatively low computational cost and which in some cases leads to a correction of previous best theoretical estimates by almost a factor of 2, qualitatively changing the predicted exciton physics. These results call for the adoption of nonuniform sampling methods for ab initio GW-BSE calculations and for revisiting previously computed values for exciton binding energies of many systems.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Random insights into the complexity of two-dimensional tensor network calculations

Projected entangled pair states (PEPS) offer memory-efficient representations of some quantum many-body states that obey an entanglement area law and are the basis for classical simulations of ground states in two-dimensional (2d) condensed matter systems. However, rigorous results show that exactly computing observables from a 2d PEPS state is generically a computationally hard problem. Yet approximation schemes for computing properties of 2d PEPS are regularly used, and empirically seen to succeed, for a large subclass of (“not too entangled”) condensed matter ground states. Adopting the philosophy of random matrix theory, in this work, we analyze the complexity of approximately contracting a 2d random PEPS by exploiting an analytic mapping to an effective replicated statistical mechanics model that permits a controlled analysis at a large bond dimension. Through this statistical-mechanics lens, we argue that (i) although approximately sampling wave-function amplitudes of random PEPS faces a computational-complexity phase transition above a critical bond dimension, and (ii) one can generically efficiently estimate the norm and correlation functions for any finite bond dimension. Furthermore, these results are supported numerically for various bond-dimension regimes. It is an important open question whether the above results for random PEPS apply more generally also to PEPS representing physically relevant ground states.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Protected fermionic zero modes in periodic gauge fields

It is well known that macroscopically normalizable zero-energy wave functions of spin-$\frac{1}{2}$ particles in a two-dimensional inhomogeneous magnetic field are spin-polarized and exactly calculable with degeneracy equaling the number of flux quanta linking the whole system. Here, extending this argument to massless Dirac fermions subjected to magnetic fields that have zero net flux but are doubly periodic in real space, we show that there exist only two Bloch-normalizable zero-energy eigenstates, one for each spin flavor. This result is immediately relevant to graphene multilayer systems subjected to doubly periodic strain fields, which at low energies enter the Hamiltonian as periodic pseudogauge vector potentials. Furthermore, we explore various related settings including nonlinearly dispersing band structure models and systems with singly periodic magnetic fields.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Second harmonic generation in topological insulators under quantizing magnetic fields

We theoretically investigate the second harmonic generation (SHG) of topological insulator surface states in a perpendicular magnetic field. Our theory is based on the microscopic expression of the second-order magneto-optical conductivity developed from the density matrix formalism, taking into account hexagonal warping effects on the surface states' band structure. Using numerically exact Landau-level energies and wave functions including hexagonal warping, we calculate the spectrum of SHG conductivities under normal incidence for different values of magnetic field and chemical potential. The imaginary parts of the SHG conductivities show prominent resonant peaks corresponding to one-photon and two-photon inter-Landau-level transitions. Treating the hexagonal warping term perturbatively, these transitions are clarified analytically within a perturbation theory from which approximate selection rules for the allowable optical transitions for SHG are determined. Furthermore, our results show extremely high SHG susceptibility that is easily tunable by magnetic field and doping level for topological surface states in the far-infrared regime, exceeding that of many conventional nonlinear materials. This work highlights the key role of hexagonal warping effects in generating second-order optical responses and provides insights into the nonlinear magneto-optical properties of the topological insulators.

Landau levels↗
DOE OSTI2020

Covariant spectator theory of n p scattering: Deuteron form factors

The deuteron form factors are calculated using two model wave functions obtained from the 2007 CST high precision fits to $np$ scattering data. Included in the calculation are a new class of isoscalar $np$ interaction currents which are automatically generated by the nuclear force model used in these fits. If the nuclear model WJC2 is used, a precision fit ($\chi^2$/datum $\eqsim1$) to the Sick Global Analysis (GA) of all $ed$ elastic scattering data can be obtained by adjusting the unknown off-shell nucleon form factors $F_3(Q^2)$ (discussed before) and $F_4(Q^2)$ (introduced in this paper), and predicting the high $Q^2$ behavior of the neutron charge form factor $G_{En}(Q^2)$ well beyond the region where it has been measured directly. Relativistic corrections, isoscalar interaction currents, and off-shell effects are defined, discussed, and their size displayed. A rationale for extending $ed$ elastic scattering measurements to higher $Q^2$ is presented.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Two-neutron knockout as a probe of the composition of states in Mg 22 , Al 23 , and Si 24

Simpson and Tostevin proposed that the width and shape of exclusive parallel momentum distributions of the A - 2 residue in direct two-nucleon knockout reactions carry a measurable sensitivity to the nucleon single-particle configurations and their couplings within the wave functions of exotic nuclei. We further report here on the first benchmarks and use of this new spectroscopic tool. Exclusive parallel momentum distributions for states in the neutron-deficient nuclei Mg 22 , Al 23 , and Si 24 populated in such direct two-neutron removal reactions were extracted and compared to predictions combining eikonal reaction theory and shell-model calculations. For the well-known Mg 22 and Al 23 nuclei, measurements and calculations were found to agree, supporting the dependence of the parallel momentum distribution width on the angular momentum composition of the shell-model two-neutron amplitudes. In Si 24 , a level at 3439(9) keV, of relevance for the important Al 23 ( p , γ ) Si 24 astrophysical reaction rate, was confirmed to be the 2 2 + state, whereas the 4 1 + state, expected to be strongly populated in two-neutron knockout, was not observed. This puzzle is resolved by theoretical considerations of the Thomas-Ehrman shift, which also indicates that a previously reported 3471-keV state in Si 24 is, in fact, the ( 0 2 + ) level with one of the largest experimental mirror-energy shifts ever observed.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Possibility of 14 C cluster as a building block of medium-mass nuclei

The possibility of the 14 C cluster being a basic building block of medium-mass nuclei is discussed. Although α cluster structures have been widely discussed in the light N ≈ Z mass region, the neutron-to-proton ratio deviates from unity in the nuclei near β-stability line and in neutron-rich nuclei. Thus, more neutron-rich objects with N > Z could become the building blocks of cluster structures in such nuclei. The 14 C nucleus is strongly bound and can be regarded as such a candidate. In addition, the path to the lowest shell-model configuration at short relative distances is closed for the 14 C+ 14 C structure, contrary to the case of the 12 C+ 12 C structure; this allows an appreciable separation distance between the 14 C clusters. The recent development of the antisymmetrized quasicluster model (AQCM) allows us to utilize a j j-coupling shell-model wave function for each cluster in a simplified way. The AQCM results for the 14 C+ 14 C structure in 28 Mg are compared with the ones of cranked relativistic mean field (CRMF) calculations. Although theoretical frameworks of these two models are quite different, they give similar results for the nucleonic densities and rotational properties of the structure under investigation. The existence of a linear chain three- 14 C cluster structure in 42 Ar has also been predicted in AQCM. Furthermore, these results confirm the role of the 14 C cluster as a possible building block of cluster structures in medium-mass nuclei.

20 ≤ A ≤ 38↗
DOE OSTI2020

Basis light-front quantization for a chiral nucleon-pion Lagrangian

We present the first application of the basis light-front quantization method to a simple chiral model of the nucleon-pion system as a relativistic bound state for the physical proton. In this work, the light-front mass-squared matrix of the nucleon-pion system is obtained within a truncated basis. The mass and the corresponding light-front wave function (LFWF) of the proton are computed by numerical diagonalization of the resulting mass-squared matrix. With the boost-invariant LFWF, we calculate the probability density distribution of the pion's longitudinal momentum fraction and the Dirac form factor of the proton.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Shape evolution of neutron-rich 106,108,110 Mo isotopes in the triaxial degree of freedom

Background: Neutron-rich nuclei with mass number between 100 and 110 attract much attention, since several kinds of shapes, such as spherical, prolate, oblate, and triaxial shapes, are predicted. In particular, for neutron-rich Mo isotopes, different models predict different magnitudes and rigidity of triaxial deformation. Previous interpretations of experimental results based solely on low-lying $2^+_2$ states are insufficient to distinguish between the rigid triaxial shape, $\gamma$ vibration, or $\gamma$-soft rotor. Purpose: The shape evolution of 106 Mo, 108 Mo, and 110 Mo is investigated through their $2^+_1$-state lifetimes, decay-branching ratios of the $2^+_2$ state, and energies of the low-lying collective excited states with $K^π = 0^+, 2^+$, and $4^+$. Method: $\beta$-delayed $\gamma$-ray spectroscopy was employed for neutron-rich Nb and Zr isotopes produced at the RIKEN RI Beam Factory to populate excited states in 106 Mo, 108 Mo, and 110 Mo. The EUroball-RIKEN Cluster Array was used for high-resolution $\gamma$-ray detection and lifetimes of the $2^+_1$ states were determined using the UK fast-timing array of LaBr 3 (Ce) detectors. Results: New $\gamma$-ray transitions and levels are reported, including newly assigned $0^+_2$ states in 108,110 Mo. Quadrupole deformations were obtained for 106,108,110 Mo from their $2^+_1$ energies and lifetimes. The $\beta$-delayed neutron-emission probabilities of 108 Nb and 110 Nb were determined by examining the $\gamma$ rays of their respective daughter decays. Conclusions: In this work, the even-odd energy staggering in the $2^+_2$ band was compared with typical patterns of the $\gamma$-vibrational band, rigid triaxial rotor, and $\gamma$-soft rotor. The very small even-odd staggering of 106 Mo, 108 Mo, and 110 Mo favors a $\gamma$-vibrational band assignment. The kinematic moment of inertia for the $2^+_2$ band showed a trend similar to the ground-state band, which is as expected for the $\gamma$-vibrational band. Beyond-mean-field calculations employing the constrained Hartree-Fock-Bogoliubov and local quasiparticle-random-phase approximation method using the SLy5 + T interaction reproduced the ground and $2^+_2$ bands in 106 Mo and 108 Mo. The collective wave functions are consistent with the interpretation of the $2^+_2$ band as the $\gamma$-vibrational band of the prolate shape. However, the staggering pattern observed in 110 Mo differs from the one suggested in the calculations which predict a $\gamma$-soft rotor. There was no experimental indication of the oblate shape or the $\gamma$-soft rotor predicted in these Mo isotopes.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗