SEARCH · Engineering Papers

Results for “PHONONS”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 829 records · Page 46

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Unraveling excitation energy transfer assisted by collective behaviors of vibrations

Abstract We investigate how collective behaviors of vibrations such as cooperativity and interference can enhance energy transfer in a nontrivial way, focusing on an example of a donor–bridge–acceptor trimeric chromophore system coupled to two vibrational degrees of freedom. Employing parameters selected to provide an overall uphill energy transfer from donor to acceptor, we use numerical calculations of dynamics in a coupled exciton–vibration basis, together with perturbation-based analytics and calculation of vibronic spectra, to identify clear spectral features of single- and multi-phonon vibrationally-assisted energy transfer (VAET) dynamics, where the latter include up to six-phonon contributions. We identify signatures of vibrational cooperation and interference that provide enhancement of energy transfer relative to that obtained from VAET with a single vibrational mode. We observe a phononic analogue of two-photon absorption, as well as a novel heteroexcitation mechanism in which a single phonon gives rise to simultaneous excitation of both the trimeric system and the vibrational degrees of freedom. The impacts of vibrations and of the one- and two-phonon VAET processes on the energy transfer are seen to be quite different in the weak and strong site–vibration coupling regimes. In the weak coupling regime, two-phonon processes dominate, whereas in the strong coupling regime up to six-phonon VAET processes can be induced. The VAET features are seen to be enhanced with increasing temperature and site–vibration coupling strength, and are reduced in the presence of dissipation. We analyze the dependence of these phenomena on the explicit form of the chromophore–vibration couplings, with comparison of VAET spectra for local and non-local couplings.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2021

Anharmonic theory of superconductivity in the high-pressure materials

Electron-phonon superconductors at high pressures have displayed the highest values of critical superconducting temperature Tc on record, now rapidly approaching room temperature. Despite the importance of high-P superconductivity in the quest for room-temperature superconductors, a mechanistic understanding of the effect of pressure and its complex interplay with phonon anharmonicity and superconductivity is missing, as numerical simulations can bring only system-specific details, clouding out key players controlling the physics. In this work, we develop a minimal model of electron-phonon superconductivity under an applied pressure which takes into account the anharmonic decoherence of the optical phonons. We find that T c behaves nonmonotonically as a function of the ratio Γ/ω 0 , where Γ is the optical phonon damping and ω 0 is the optical phonon energy at zero pressure and momentum. Optimal pairing occurs for a critical ratio Γ/ω 0 when the phonons are on the verge of decoherence (“diffusonlike” limit). Our framework gives insights into recent experimental observations of T c as a function of pressure in the complex BCS material TlInTe 2 .

36 MATERIALS SCIENCE↗

DOE OSTI2024

Thermal transport in Yb-based 1-2-20 materials

Abstract Lattice thermal-transport properties of heavy-fermion YbT 2 Zn 20 -based (T = Co, Rh, and Ir) 1-2-20 compounds are calculated with the single-mode relaxation-time approximation and the full solution of the linearized phonon Boltzmann transport equation from first-principles anharmonic phonon calculations. We predict low lattice thermal conductivity κ L with a maximum value of ∼ 24.16 · K at 85.14 K, 137.29 · K at 9.67 K, and 23.55 · K at 67.50 K and ∼ 12.02, 40.04, and 10.30 · K at room temperature for YbCo 2 Zn 20 , YbRh 2 Zn 20 , and YbIr 2 Zn 20 , respectively. Based on the analysis of the cumulative κ L as a function of the phonon mean free path and the frequency-dependent joint density of states at various temperatures, we attribute the low κ L to a rattling mode with an average characteristic rattling frequency that enhanced phonon scattering processes, which reduced the phonon mean free paths, suppressed the phonon lifetime, and enhanced the probability of three-phonon scattering events. The predicted low κ L , especially in YbCo 2 Zn 20 and YbIr 2 Zn 20 , makes them promising candidate materials for thermoelectric applications and thermal management.

Physics↗

DOE OSTI2025

Crystalline Order Yet Glass-Like Heat Transport Driven by Hidden Local Distortions as the Structural Origin of Ultralow Thermal Conductivity in AgErTe 2

Given their rich chemical diversity and the interplay among the p-, d-, and f-orbitals of chalcogens, transition metals, and lanthanides, respectively, rare-earth transition-metal chalcogenides exhibit a wide variety of structural, magnetic, and transport phenomena. As a result, they form a particularly appealing platform for investigating structure–property relationships, emergent electronic and magnetic behaviors, and thermal transport. Here we investigate AgErTe 2 as a model system to understand phonon-glass behavior in ordered crystalline solids, which establishes the design principles for thermal barrier coatings and next-generation thermoelectrics. The local bonding asymmetry and lattice softness suppress the inherently low lattice thermal conductivity, resembling the characteristics of amorphous materials. This suppression is significantly influenced by local off-centering of Ag atoms, which breaks lattice periodicity while maintaining global crystallinity. The presence of antibonding states just below the Fermi level, arising from Ag 4d and Te 5p orbital interactions, leads to lattice softening and destabilizes ideal tetrahedral coordination, resulting in a pseudo Jahn–Teller distortion. Furthermore, the coexistence of weaker, more polarizable Ag–Te bonds and stronger Er–Te bonds creates a complex vibrational landscape enriched with low-frequency modes and enhanced phonon scattering. A pronounced disparity in interatomic force constants gives rise to highly localized, low-energy optical phonons linked to Ag rattling. These flat vibrational modes exhibit strong coupling with transverse acoustic phonons, resulting in ultrashort phonon lifetimes and mean free paths approaching interatomic distances. These features collectively enhance phonon scattering across a broad range of length and energy scales. This work offers a framework for engineering suppressed thermal conductivity in crystalline systems without the introduction of alloying elements.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Optimal carrier concentration for FeSb 2 colossal thermopower

Crystals of FeSb 2 correlated narrow-gap semiconductor host colossal thermopower values. By tuning the impurity level here, we demonstrate that electron-phonon scattering that transfers phonon momentum to electrons is efficient only for certain optimal carrier concentration in the low-mobility band. Phonon drag acting on such states in crystals with high phonon mean free path enhances thermopower to colossal values, whereas for different carrier concentration, dominant thermal transport mechanism is electronic diffusion. This highlights the dual nature of correlated in-gap states that take part in the phonon drag but also reduce phonon mean free path.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2022

Direct visualization of polaron formation in the thermoelectric SnSe

SnSe is a layered material that currently holds the record for bulk thermoelectric efficiency. The primary determinant of this high efficiency is thought to be the anomalously low thermal conductivity resulting from strong anharmonic coupling within the phonon system. Here we show that the nature of the carrier system in SnSe is also determined by strong coupling to phonons by directly visualizing polaron formation in the material. We employ ultrafast electron diffraction and diffuse scattering to track the response of phonons in both momentum and time to the photodoping of free carriers across the bandgap, observing the bimodal and anisotropic lattice distortions that drive carrier localization. Relatively large (18.7 Å), quasi-one-dimensional (1D) polarons are formed on the 300-fs timescale with smaller (4.2 Å) 3D polarons taking an order of magnitude longer (4 ps) to form. This difference appears to be a consequence of the profoundly anisotropic electron–phonon coupling in SnSe, with strong Fröhlich coupling only to zone-center polar optical phonons. These results demonstrate a high density of polarons in SnSe at optimal doping levels. Strong electron-phonon coupling is critical to the thermoelectric performance of this benchmark material and, potentially, high performance thermoelectrics more generally.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Extracting correlation effects from momentum-resolved electron energy loss spectroscopy: Synergistic origin of the dispersion kink in Bi 2.1 Sr 1.9 CaCu 2 O 8 + x

We employ momentum-resolved electron energy loss spectroscopy (M-EELS) on Bi 2.1 Sr 1.9 CaCu 2 O 8 + x to resolve the issue of the kink feature in the electron dispersion widely observed in the cuprates. To this end, we utilize the GW approximation to relate the density response function measured in M-EELS to the self-energy, isolating contributions from phonons, electrons, and the momentum dependence of the effective interaction to the decay rates. The phononic contributions, present in the M-EELS spectra due to electron-phonon coupling, lead to kink features in the corresponding single-particle spectra at energies between 40 and 80 meV, independent of the doping level. We find that a repulsive interaction constant in momentum space is able to yield the kink attributed to phonons in ARPES. Hence, our analysis of the M-EELS spectra points to local repulsive interactions as a factor that enhances the spectroscopic signatures of electron-phonon coupling in cuprates. We conclude that the strength of the kink feature in cuprates is determined by the combined action of electron-phonon coupling and electron-electron interactions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2021

Thermodynamic properties of ε -Fe with thermal electronic excitation effects on vibrational spectra

The thermodynamic properties of hexagonal-close-packed iron (ε–Fe) are essential for investigating the internal structure and dynamic properties of planetary cores. Despite their importance to planetary sciences, experimental investigations of ε–Fe at relevant conditions are still challenging. Therefore, ab initio calculations are crucial to elucidating the thermodynamic properties of this system. Here, we use a free energy calculation scheme based on the phonon gas model compatible with temperature-dependent phonon frequencies. We investigate the effects of electronic thermal excitations, which introduces a temperature dependence on phonon frequencies, and the implication for the thermodynamic properties of ε–Fe at extreme pressure (P) and temperature (T) conditions. We disregard phonon-phonon interactions, i.e., anharmonicity and their effect on phonon frequencies. Nevertheless, the current scheme is also applicable to T -dependent anharmonic frequencies. We conclude that the impact of thermal electronic excitations on vibrational properties is not significant up to ~4000 K at 200 GPa but should not be ignored at higher temperatures or pressures. However, the static free energy F st must always include the effect of thermal excitation fully in a continuum of T. Furthermore, our results for isentropic equations of state show good agreement with data from recent ramp compression experiments up to 1400 GPa conducted at the National Ignition Facility.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2025

Temperature-dependent lattice dynamics in polycrystalline ZrN

Zirconium nitride holds promise for improving nuclear reactor components, where thermal conductivity is important, highlighting the need for a detailed understanding of its phonons. However, the phonon dispersion derived from previous inelastic neutron scattering data is in disagreement with modern theoretical calculations, pointing to the necessity for further experimental validation. Here, in this study, we use inelastic neutron scattering experiments on polycrystalline ZrN at 6, 100, 300, 400, and 600 K, and density functional theory (DFT) calculations to study the momentum-resolved scattering and the phonon density of states. Both experimental and theoretical results reveal a phonon cutoff energy near 65 meV, in stark disagreement with previous single-crystal neutron scattering measurements indicating a cutoff near 75 meV. The 65 meV cutoff energy is, however, in reasonable agreement with previous measurements on ZrN 0.9 . Consistency between our theoretical and experimental results confirms that ZrN is a wide phonon band gap system and that DFT provides a reliable description of its phonons.

DFT↗

DOE OSTI2025

Atomic dynamics in 𝑀⁢Cr⁢𝑋 2 (𝑀=Ag, Cu; 𝑋 = S, Se) across magnetic and superionic transitions

Here, a systematic study of atomic dynamics and thermal properties of the family of layered chalcogenide compounds 𝑀⁢Cr⁢𝑋 2 (𝑀= Ag, Cu; 𝑋 = S, Se) was performed, including neutron and x-ray scattering, thermal characterization, and first-principles simulations. In all compounds, we observe a breakdown of specific phonon modes across the superionic phase transition, for phonons whose eigenvectors exhibit large contributions of mobile ions. In particular, the nondispersive portions of transverse acoustic (TA) branches at short-wavelengths and the low-energy optical phonons with large contributions from Ag + or Cu + become severely damped in the superionic phase. However, well-defined quasiparticles persist in the superionic state for long-wavelength TA phonons. In the case of AgCrS 2 , the coupling of lattice dynamics with its antiferromagnetic transition was also investigated. The magnetic ordering couples with the monoclinic–rhombohedral structural transition, and the Cr 3+ spin arrangement strongly affects the phonon dispersions. We qualitatively reproduce the magnetic and nuclear components of the INS measurement for antiferromagnetic AgCrS 2 by combining models of spin-waves and spin-polarized first-principles phonon simulations. Quasielastic magnetic fluctuations persist in the paramagnetic phase up to high temperature, but are clearly distinguished from the nuclear component through their momentum dependence. Finally, we report measurements of the thermal properties of the selenide compounds and find good agreement with our DFT simulations.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Evaluation of Thermal Neutron Scattering Cross Section of Uranium Silicide with Ab Initio Lattice Dynamics

Uranium silicide (U 3 Si 2 ) is a candidate material for the high-density nuclear fuel in commercial light water reactors [1], [2]. Its higher uranium density, 11.3 g-U/cm3, compared to that of uranium dioxide (UO 2 ), 9.7 g-U/cm3, can improve the performance of a nuclear reactor while using low enriched uranium (LEU) and diversify the choice of cladding materials [1]–[3]. It also has a higher thermal conductivity than UO 2 , which can reduce the thermal stress on the material caused by a temperature gradient across the fuel pellet and provide a larger margin for some postulated accidents [1], [2], [4]. Furthermore, compared to U3Si, another high-density fuel candidate, it has better resistance to in-pile swelling due to less irradiation-induced rapid amorphization [1], [3]. Corresponding to its importance in nuclear engineering, many previous studies have reported the properties of U3Si2. Experiments showed that U 3 Si 2 is a paramagnetic (PM) metal, where a slight linear increase in magnetic susceptibility was measured with increasing temperature [5], [6]. In addition, thermodynamic quantities such as thermal expansion coefficient, heat capacity, and thermal conductivity were experimentally determined over a wide temperature range [1], [7], [8]. In several computational studies, ab initio atomistic simulations based on density functional theory (DFT) were performed to calculate various properties including elastic constants, electronic density of states (DOS), and phonon dispersion curves [9]–[12]. Nevertheless, thermal neutron scattering cross sections, which are critical to the prediction of the parameters in reactor physics that are ultimately related to reactor criticality, have not yet been evaluated for U3Si2. The scattering cross section can be calculated from the phonon DOS, or the energy spectrum of lattice vibrations, of the crystalline system [13], [14]. However, there is also no experimental data available for the phonon DOS of U 3 Si 2 . While some computational studies reported the phonon DOS and/or dispersion curves from ab initio simulations [9]–[12], the accuracy cannot be guaranteed because it is unclear whether the spin-polarization behavior of PM U 3 Si 2 was properly described. In the present study, the thermal neutron scattering cross section for U 3 Si 2 is evaluated for the first time by calculating the phonon DOS for U3Si2 from ab initio lattice dynamics (AILD) simulations based on DFT. First, U 3 Si 2 is modeled based on the experimental structure, and AILD simulations are performed on the modeled U3Si2 to optimize the structure. Next, AILD simulations are performed for supercells with atomic displacement to calculate Hellmann-Feynman forces. Based on the calculated forces, partial phonon DOSs for U and Si are obtained, and the thermal neutron scattering law (TSL) for U 3 Si 2 is finally evaluated. To verify the accuracy of the calculations in the present study, the calculation results are compared with experimental data on the structure and heat capacity of U3Si2 [1], [7], [8], [15].

Geometry Optimization↗

DOE OSTI2024

Manipulating Spin–Lattice Coupling in Layered Magnetic Topological Insulator Heterostructure via Interface Engineering

Induced magnetic order in a topological insulator (TI) can be realized either by depositing magnetic adatoms on the surface of a TI or engineering the interface with epitaxial thin film or stacked assembly of 2D van der Waals (vdW) materials. Herein, the observation of spin-phonon coupling in the otherwise non-magnetic TI Bi 2 Te 3 is reported, due to the proximity of FePS 3 (an antiferromagnet (AFM), T N ≈ 120 K), in a vdW heterostructure framework. Temperature-dependent Raman spectroscopic studies reveal deviation from the usual phonon anharmonicity originated from spin-lattice coupling at the Bi 2 Te 3 /FePS 3 interface at/below 60 K in the peak position (self-energy) and linewidth (lifetime) of the characteristic phonon modes of Bi 2 Te 3 (106 and 138 cm –1 ) in the stacked heterostructure. The Ginzburg-Landau (GL) formalism, where the respective phonon frequencies of Bi 2 Te 3 couple to phonons of similar frequencies of FePS 3 in the AFM phase, is adopted to understand the origin of the hybrid magneto-elastic modes. At the same time, the reduction of characteristic TN of FePS 3 from 120 K in isolated flakes to 65 K in the heterostructure, possibly due to the interfacial strain, which leads to smaller Fe-S-Fe bond angles as corroborated by computational studies using density functional theory (DFT). Besides, inserting hexagonal boron nitride within Bi 2 Te 3 /FePS 3 stacking regains the anharmonicity in Bi 2 Te 3 . As a result, controlling interfacial spin-phonon coupling in stacked heterostructure can have potential application in surface code spin logic devices.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Dislocation-Limited Thermal Conductivity in LiF: Revisiting Perturbative Models

We revisit models to describe phonon-dislocation interactions and resulting dislocation-limited thermal conductivity (k), particularly as applied to describe the temperature-dependent k behavior of LiF with significant dislocation densities. Coupling semiempirical models of phonon scattering from dislocation strain fields provided by Klemens and Carruthers with density functional theory description of the phonons and input parameters, we find that both models significantly overpredict the measured k of LiF, as found previously. However, more direct application of the quantum perturbation theory description, from which the models were derived, gives significantly stronger phonon-dislocation scattering and strongly underpredicts the measured k data. We revisit the derivation of phonon-dislocation-strain field scattering and provide numerical details regarding its implementation. This work demonstrates that quantum perturbative calculations, now accessible to modern computational architectures, may provide a reasonable description of dislocation-limited k in materials as governed by phonon-strain field interactions provided that questions regarding the relevant strain field range can be resolved.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Microscopic mechanism of unusual lattice thermal transport in TlInTe 2

We investigate the microscopic mechanism of ultralow lattice thermal conductivity (κ l ) of TlInTe 2 and its weak temperature dependence using a unified theory of lattice heat transport, that considers contributions arising from the particle-like propagation as well as wave-like tunneling of phonons. While we use the Peierls–Boltzmann transport equation (PBTE) to calculate the particle-like contributions (κ l (PBTE)), we explicitly calculate the off-diagonal (OD) components of the heat-flux operator within a first-principles density functional theory framework to determine the contributions (κ l (OD)) arising from the wave-like tunneling of phonons. At each temperature, T, we anharmonically renormalize the phonon frequencies using the self-consistent phonon theory including quartic anharmonicity, and utilize them to calculate κ l (PBTE) and κ l (OD). With the combined inclusion of κ l (PBTE), κ l (OD), and additional grain-boundary scatterings, our calculations successfully reproduce the experimental results. Our analysis shows that large quartic anharmonicity of TlInTe 2 (a) strongly hardens the low-energy phonon branches, (b) diminishes the three-phonon scattering processes at finite T, and (c) recovers the weaker than T –1 decay of the measured κ l .

36 MATERIALS SCIENCE↗

DOE OSTI2022

Anomalous thermal transport under high pressure in boron arsenide

High pressure represents extreme environments and provides opportunities for materials discovery. Thermal transport under high hydrostatic pressure has been investigated for more than 100 years and all measurements of crystals so far have indicated a monotonically increasing lattice thermal conductivity. Here we report in situ thermal transport measurements in the newly discovered semiconductor crystal boron arsenide, and observe an anomalous pressure dependence of the thermal conductivity. We use ultrafast optics, Raman spectroscopy and inelastic X-ray scattering measurements to examine the phonon bandstructure evolution of the optical and acoustic branches, as well as thermal conductivity under varied temperatures and pressures up to 32 gigapascals. Using atomistic theory, we attribute the anomalous high-pressure behaviour to competitive heat conduction channels from interactive high-order anharmonicity physics inherent to the unique phonon bandstructure. Our study verifies ab initio theory calculations and we show that the phonon dynamics—resulting from competing three-phonon and four-phonon scattering processes—are beyond those expected from classical models and seen in common materials. In conclusion, this work uses high-pressure spectroscopy combined with atomistic theory as a powerful approach to probe complex phonon physics and provide fundamental insights for understanding microscopic energy transport in materials of extreme properties.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2023

Advanced Modeling and Simulation Methods for Evaluation of Thermal Neutron Scattering Materials

With the rise of interest in thermal neutron scattering data for advanced reactor, criticality safety, and shielding applications, new experimental data are required for evaluation of new materials or for re-evaluation (or validations) of previously evaluated materials. New experimental data are evaluated in a three-step process: (1) computing the phonon characteristics, (2) computing the dynamic structure factor (DSF) from the data, and (3) using the experimental setup to simulate the experimental data. All three steps have challenges, ranging from the need for a sufficiently general material simulation code—a processing code that can compute the corresponding DSF—to having a detailed layout of the instrument/beamline/facility where the data were measured. Whereas phonon characteristics of materials can be computed using various methods (molecular dynamics, density functional theory, etc.), a high-fidelity computation of the DSF and the simulation of the experiment based on the DSF is vital to the accuracy of the evaluation. The latter two steps can be achieved by using the two corresponding code systems developed by instrument scientists at the Spallation Neutron Source (SNS) at Oak Ridge National Laboratory: (1) OCLIMAX, a program that calculates the dynamic structure factor from DFT and MD simulation results, and (2) MCViNE, a Monte Carlo neutron ray-tracing program designed to simulate neutron scattering experiments. Recently, polyethylene and yttrium hydride were measured at the Wide Angular-Range Chopper (ARCS) and SEQUOIA instrument stations of the SNS. These experiments are simulated using the density functional theory code, the Cambridge Serial Total Energy Package (CASTEP), to compute its phonon characteristics (eigenvalues/vectors and PDOS), which is then processed using OCLIMAX to yield the DSF, and finally the data at each instrument station are simulated by the MCViNE for comparison to the measured data for evaluation. For comparison to conventional evaluation methods, the scattering data processed from OCLIMAX are compared against those processed from the LEAPR module of NJOY, and the results from MCViNE simulations are compared against previously used simplified beamline models implemented in the Monte Carlo N-Particle (MCNP) code.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2024

Chemical bonding dictates drastic critical temperature difference in two seemingly identical superconductors

Though YB 6 and LaB 6 share the same crystal structure, atomic valence electron configuration, and phonon modes, they exhibit drastically different phonon-mediated superconductivity. YB 6 superconducts below 8.4 K, giving it the second-highest critical temperature of known borides, second only to MgB 2 . LaB 6 does not superconduct until near-absolute zero temperatures (below 0.45 K), however. Though previous studies have quantified the canonical superconductivity descriptors of YB 6 ’s greater Fermi-level (E f ) density of states and higher electron–phonon coupling (EPC), the root of this difference has not been assessed with full detail of the electronic structure. Through chemical bonding, we determine low-lying, unoccupied 4f atomic orbitals in lanthanum to be the key difference between these superconductors. These orbitals, which are not accessible in YB 6 , hybridize with π B–B bonds and bring this π-system lower in energy than the σ B–B bonds otherwise at E f . This inversion of bands is crucial: the optical phonon modes we show responsible for superconductivity cause the σ-orbitals of YB 6 to change drastically in overlap, but couple weakly to the π-orbitals of LaB 6 . These phonons in YB 6 even access a crossing of electronic states, indicating strong EPC. No such crossing in LaB 6 is observed. Finally, a supercell (the M k-point) is shown to undergo Peierls-like effects in YB 6 , introducing additional EPC from both softened acoustic phonons and the same electron-coupled optical modes as in the unit cell. Overall, we find that LaB 6 and YB 6 have fundamentally different mechanisms of superconductivity, despite their otherwise near-identity.

Science & Technology - Other Topics↗

DOE OSTI2020

Charge-density-wave melting in the one-dimensional Holstein model

We study the Holstein model of spinless fermions, which at half filling exhibits a quantum phase transition from a metallic Tomonaga-Luttinger liquid phase to an insulating charge-density-wave (CDW) phase at a critical electron-phonon coupling strength. In our work, we focus on the real-time evolution starting from two different types of initial states that are CDW ordered: (i) ideal CDW states with and without additional phonons in the system and (ii) correlated ground states in the CDW phase. We identify the mechanism for CDW melting in the ensuing real-time dynamics and show that it strongly depends on the type of initial state. We focus on the far-from-equilibrium regime and emphasize the role of electron-phonon coupling rather than dominant electronic correlations, thus complementing a previous study of photoinduced CDW melting [H. Hashimoto and S. Ishihara, Phys. Rev. B 96, 035154 (2017)]. The numerical simulations are performed by means of matrix-product-state based methods with a local basis optimization (LBO). Within these techniques, one rotates the local (bosonic) Hilbert spaces adaptively into an optimized basis that can then be truncated while still maintaining a high precision. In this work, we extend the time-evolving block decimation (TEBD) algorithm with LBO, previously applied to single-polaron dynamics, to a half-filled system. We demonstrate that in some parameter regimes, a conventional TEBD method without LBO would fail. Furthermore, we introduce and use a ground-state density-matrix renormalization group method for electron-phonon systems using local basis optimization. In our examples, we account for up to M ph =40 bare phonons per site by working with O(10) optimal phonon modes.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗