Real-time cumulant approach for charge-transfer satellites in x-ray photoemission spectra
Not Available
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.
Not Available
We studied the lattice vibrations of two inter-penetrating atomic sublattices via the Frenkel-Kontorova (FK) model of a linear chain of harmonically interacting atoms subjected to an on-site potential, using the technique of thermodynamic Green's functions based on quantum field-theoretical methods. General expressions were deduced for the phonon frequency-wave-vector dispersion relations, number density, and energy of the FK model system. In addition, as the application of the theory, we investigated in detail cases of linear chains with various periods of the on-site potential of the FK model. Some unusual but interesting features for different amplitudes of the on-site potential of the FK model are discussed. In the commensurate structure, the phonon spectrum always starts at a finite frequency, and the gaps of the spectrum are true ones with a zero density of modes. In the incommensurate structure, the phonon spectrum starts from zero frequency, but at a non-zero wave vector; there are some modes inside these gap regions, but their density is very low. In our approximation, the energy of a higher-order commensurate state of the one-dimensional system at a finite temperature may become indefinitely close to the energy of an incommensurate state. This finding implies that the higher-order incommensurate-commensurate transitions are continuous ones and that the phase transition may exhibit a “devil's staircase” behavior at a finite temperature.
Using state-of-the-art many-body Green's function calculations based on the GW plus cumulant approach, we analyze the properties of plasmon satellites in the electron spectral function resulting from electron-plasmon interactions in one-, two-, and three-dimensional systems. Specifically, we show how their dispersion relation, line shape, and linewidth are related to the properties of the constituent electrons and plasmons. In addition, to gain insight into the many-body processes giving rise to the formation of plasmon satellites, we connect the GW plus cumulant approach to a many-body wave-function picture of electron-plasmon interactions and introduce the coupling-strength-weighted electron-plasmon joint density states as a powerful concept for understanding plasmon satellites.
By first-principles calculations, we find that Ba 2 X (X = Si, Ge) hosts a topological semimetal phase with one nodal ring in the k x = 0 plane, which is protected by the glide mirror symmetry when spin-orbit coupling (SOC) is ignored. The corresponding drumheadlike surface flat band appears on the (100) surface in surface Green's function calculation. Furthermore, a topological-semimetal-to-insulator transition is found. The nodal line semimetal would evolve into a topological insulator as SOC is turned on. The topologically protected metallic surface states emerge around the Γ point, which could be tuned into the topologically trivial insulator state by more than 3% hydrostatic strain. These results reveal a category of materials showing quantum phase transition between topological semimetal and insulator, and tunability through elastic strain engineering.
High resolution inelastic x-ray scattering experiments are widely used to study the electronic and chemical properties of materials under a range of conditions, from ambient temperature to the warm dense matter regime. Here, we use the real-space multiple scattering (RSMS) Green's function formalism coupled with density functional theory molecular dynamics (DFT-MD) to study thermal effects on the Compton profile (CP) of disordered systems. The RSMS method is advantageous for calculations of highly disordered, aperiodic systems because it places no restriction on symmetry. As a test, we apply our approach to thermally disordered Be, Li, and Si in both liquid and solid phases. We find good agreement with experimental and other theoretical results, showing that the real-space multiple scattering approach coupled with DFT-MD is an efficient and reliable method for calculating the CP of disordered systems at finite temperatures.
Single-atom gating, achieved by manipulation of adatoms on a surface, has been shown in experiments to allow precise control over superposition of electronic states in quantum corrals. Using a Green's function approach, we demonstrate theoretically that such atom gating can also be used to control the coupling between magnetic degrees of freedom in these systems. Atomic gating enables control not only on the direct interaction between magnetic adatoms, but also over superpositions of many-body states which can then control long distance interactions. We illustrate this effect by considering the competition between direct exchange between magnetic impurities and the Kondo screening mediated by the host electrons, and how this is affected by gating. These results suggest that both magnetic and nonmagnetic single-atom gating may be used to investigate magnetic impurity systems with tailored interactions, and may allow the control of entanglement of different spin states.
The effect of disorder on lattice vibrational modes has been a topic of interest for several decades. In this article, we employ a Green's function based approach, namely, the dynamical cluster approximation (DCA), to investigate phonons in mass-disordered systems. Detailed benchmarks with previous exact calculations are used to validate the method in a wide parameter space. An extension of the method, namely, the typical medium DCA (TMDCA), is used to study Anderson localization of phonons in three dimensions. We show that, for binary isotopic disorder, lighter impurities induce localized modes beyond the bandwidth of the host system, while heavier impurities lead to a partial localization of the low-frequency acoustic modes. For a uniform (box) distribution of masses, the physical spectrum is shown to develop long tails comprising mostly localized modes. The mobility edge separating extended and localized modes, obtained through the TMDCA, agrees well with results from the transfer matrix method. A reentrance behavior of the mobility edge with increasing disorder is found that is similar to, but somewhat more pronounced than, the behavior in disordered electronic systems. Our work establishes a computational approach, which recovers the thermodynamic limit, is versatile and computationally inexpensive, to investigate lattice vibrations in disordered lattice systems.
In this work, we present a theoretical framework and implementation details for self-energy embedding theory (SEET) with the GW approximation for the treatment of weakly correlated degrees of freedom and configuration interactions solver for handing the strongly correlated degrees. On a series of molecular examples, for which the exact results are known within a given basis, we demonstrate that SEET(CI/GW) is a systematically improvable and well controlled method capable of giving accurate results and well behaved causal self-energies and Green's functions. We compare the theoretical framework of SEET(CI/GW) to that of GW+DMFT and comment on differences between these to approaches that aim to treat both strongly and weakly correlated simultaneously.
In this article, we present a two-dimensional lattice model for quantum spin-1/2 for which the low-energy limit is governed by four flavors of strongly interacting Majorana fermions. We study this low-energy effective theory using two alternative approaches. The first consists of a mean-field approximation. The second consists of a random phase approximation (RPA) for the single-particle Green's functions of the Majorana fermions built from their exact forms in a certain one-dimensional limit. The resulting phase diagram consists of two competing chiral phases, one with Abelian and the other with non-Abelian topological order, separated by a continuous phase transition. Remarkably, the Majorana fermions propagate in the two-dimensional bulk, as in the Kitaev model for a spin liquid on the honeycomb lattice. We identify the vison fields, which are mobile (they are static in the Kitaev model) domain walls propagating along only one of the two space directions.
Following the recent success of realizing exciton-polariton condensates in cavities, we examine the hybridization of cavity photons with the closest analog of excitons within a superconductor, states called Bardasis-Schrieffer modes. Although these modes do not typically couple linearly to light, one can engineer a coupling with an externally imposed supercurrent, leading to the formation of hybridized Bardasis-Schrieffer-polariton states, which we obtain both as poles of the bosonic Green's function and through the derivation of an effective Hamiltonian picture for the model. Here, these new excitations have nontrivial overlap with both the original photon states and d-wave superconducting fluctuations. We conjecture that a phase-coherent density of these objects could produce a finite d-wave component of the superconducting order parameter—an s ± id superconducting state.
The vibrational thermal conductances (G) across GaN-AlN interfaces are computed using a nonequilibrium Green's function formalism in the harmonic limit with bulk and interfacial interatomic force constants (IFCs) fully from density functional theory. Several numerical methods and supercell configurations are employed to examine the sensitivity of G to variances of IFCs. In particular, the effects of supercell size, the enforcement of symmetry constraints, and truncation of IFCs near the interface, and atomic relaxation on phonon transmission and conductance are explored. Our fully first-principles calculations are compared with common approximations and measured G values inferred from thermal conductivity measurements for GaN-AlN superlattices. Here, our calculated value, G~300MWm -2 K -1 , is nearly half that from measurements. This discrepancy is critically analyzed in terms of the physical assumptions of the calculations and the derivation of the experimental values. This work provides guidelines to determine “physically correct” sets of interfacial IFCs from first principles for thermal conductance calculations using minimal computational resources. It also contributes toward developing predictive calculations and a more complete picture of thermal conduction across interfaces, a step toward first-principles multiscale thermal transport.
Veselago lens focusing in graphene p-n junction is promising for realizations of new generation electron optics devices. However, the effect of the strain-induced Aharonov-Bohm interference in a p-n junction has not been discussed before. We provide an experimentally feasible setup based on the Veselago lens in which the presence of strain can result in both the valley-dependent Lorentz force and Aharonov-Bohm interference. Inparticular, by employing the Green’s function and tight-binding methods, we study the strain induced by dislocations and line defects in a p-n junction and show how the resulting Aharonov-Bohm phase and interference can be detected. Furthermore, for a different strain configuration, e.g., corresponding to corrugated graphene, we find strong signatures of valley splitting induced by the fictitious magnetic field. Our proposal can be useful for mapping elastic deformations and defects, and for studying valley-dependent effects in graphene.
We study if commonly used nucleon-nucleon effective interactions, obtained from fitting the properties of cold nuclear matter and of finite nuclei, can properly describe the hot dense nuclear matter produced in intermediate-energy heavy-ion collisions. We use two representative effective interactions, i.e., an improved isospin- and momentum-dependent interaction with its isovector part calibrated by the results from the ab initio nonperturbative self-consistent Green's function (SCGF) approach with chiral forces, and a Skyme-type interaction fitted to the equation of state of cold nuclear matter from chiral effective many-body perturbation theory and the binding energy of finite nuclei. In the mean-field approximation, we evaluate the equation of state and the single-nucleon potential for nuclear matter at finite temperatures and compare them to those from the SCGF approach. We find that the improved isospin- and momentum-dependent interaction reproduces reasonably well the SCGF results due to its weaker momentum dependence of the mean-field potential than in the Skyrme-type interaction. Furthermore, our study thus indicates that effective interactions with the correct momentum dependence of the mean-filed potential can properly describe the properties of hot dense nuclear matter and are thus suitable for use in transport models to study heavy-ion collisions at intermediate energies.
Monoenergetic muon neutrinos at 236 MeV are readily produced in intense medium-energy proton facilities (≳2–3GeV) when a positive kaon decays at rest (KDAR; K + →μ + ν μ ). These neutrinos provide a unique opportunity to both study the neutrino interaction and probe the nucleus with a monoenergetic weak-interaction-only tool. We present cross-section calculations for quasielastic scattering of these 236-MeV neutrinos off 12 C and 40 Ar, paying special attention to low-energy aspects of the scattering process. Our model takes the description of the nucleus in a mean-field approach as the starting point, where we solve Hartree-Fock equations using a Skyrme type nucleon-nucleon interaction. Thereby, we introduce long-range nuclear correlations by means of a continuum random-phase approximation (CRPA) framework where we solve the CRPA equations using a Green's function method. The model successfully describes (e,e') data on 12 C and 40 Ca in the kinematic region that overlaps with the KDAR ν μ phase space. In addition to these results, we present future prospects for precision KDAR cross-section measurements and applications of our calculations in current and future experiments that will utilize these neutrinos.
Searches for neutrinoless-double beta decay rates are crucial in addressing questions within fundamental symmetries and neutrino physics. The rates of these decays depend not only on unknown parameters associated with neutrinos, but also on nuclear properties. In order to reliably extract information about the neutrino, one needs an accurate treatment of the complex many-body dynamics of the nucleus. Neutrinoless-double beta decays take place at momentum transfers on the order of 100 MeV/$c$ and require both nuclear electroweak vector and axial current matrix elements. Muon capture, a process in the same momentum transfer regime, has readily available experimental data to validate these currents. In this work, we present results of {\it ab initio} calculations of partial muon capture rates for $^3$He and $^6$Li nuclei using variational and Green's Function Monte Carlo computational methods. We estimate the impact of the three-nucleon interactions, the cutoffs used to regularize two-nucleon ($2N$) interactions, and the energy range of $2N$ scattering data used to fit these interactions.
We perform a first-principles calculation of optical potentials for nucleon elastic scattering off medium-mass isotopes. Fully based on a saturating chiral Hamiltonian, the optical potentials are derived by folding nuclear density distributions computed with self-consistent Green's function theory with a nucleon-nucleon t matrix computed with a consistent chiral interaction. The dependence on the folding interaction as well as the convergence of the target densities are investigated. Numerical results are presented and discussed for differential cross sections and analyzing powers, with focus on elastic proton scattering off calcium and nickel isotopes. Our optical potentials generally show a remarkable agreement with the available experimental data for laboratory energies in the range 65–200 MeV. We study the evolution of the scattering observables with increasing proton-neutron asymmetry by computing theoretical predictions of the cross section and analyzing power over the calcium and nickel isotopic chains. Published by the American Physical Society 2024
Background: Ab initio many-body methods have been developed over the past 10 yr to address closed-shell nuclei up to mass A≈130 on the basis of realistic two- and three-nucleon interactions. A current frontier relates to the extension of those many-body methods to the description of open-shell nuclei. Several routes to address open-shell nuclei are currently under investigation, including ideas that exploit spontaneous symmetry breaking. Purpose: Singly open-shell nuclei can be efficiently described via the sole breaking of U(1) gauge symmetry associated with particle-number conservation as a way to account for their superfluid character. While this route was recently followed within the framework of self-consistent Green's function theory, the goal of the present work is to formulate a similar extension within the framework of coupled cluster theory. Methods: We formulate and apply Bogoliubov coupled cluster (BCC) theory, which consists of representing the exact ground-state wave function of the system as the exponential of a quasiparticle excitation cluster operator acting on a Bogoliubov reference state. Equations for the ground-state energy and the cluster amplitudes are derived at the singles and doubles level (BCCSD) both algebraically and diagrammatically. The formalism includes three-nucleon forces at the normal-ordered two-body level. The first BCC code is implemented in m scheme, which will permit the treatment of doubly open-shell nuclei via the further breaking of SU(2) symmetry associated with angular momentum conservation. Results: Proof-of-principle calculations in an N max =6 spherical harmonic oscillator basis for 16,18 O and 18 Ne in the BCCD approximation are in good agreement with standard coupled cluster results with the same chiral two-nucleon interaction, while 20 O and 20 Mg display underbinding relative to experiment. The breaking of U(1) symmetry, monitored by computing the variance associated with the particle-number operator, is relatively constant for all five nuclei, in both the Hartree-Fock-Bogoliubov and BCCD approximations. Conclusions: The newly developed many-body formalism increases the potential span of ab initio calculations based on single-reference coupled cluster techniques tremendously, i.e., potentially to reach several hundred additional midmass nuclei. The new formalism offers a wealth of potential applications and further extensions dedicated to the description of ground and excited states of open-shell nuclei. Short-term goals include the implementation of three-nucleon forces at the normal-ordered two-body level. Midterm extensions include the approximate treatment of triples corrections and the development of the equation-of-motion methodology to treat both excited states and odd nuclei. Long-term extensions include exact restoration of U(1) and SU(2) symmetries.