Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “local density approximations”

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 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Effect of Composition and Local Environment on CO 2 Adsorption on Nickel and Magnesium Oxide Solid Solutions

CO 2 adsorption energies on the (100) surfaces of a nickel oxide doped with Mg, magnesium oxide doped with Ni, and their 50:50 solid solution were calculated using density functional theory. The composition and atomic arrangement of the adsorption site was varied to understand how the local environment affects CO 2 adsorption and the basicity of the surfaces at the atomic level. The dispersive and electronic contributions to the adsorption energies were quantified, and the results indicate that the variation of the adsorption energy with adsorption site configuration and metal composition is dominated by electronic interactions. Interestingly, for magnesium oxide doped with nickel, a single substitution can create stronger CO 2 binding sites, which implies stronger basic sites, even though nickel oxide is less basic than magnesium oxide. The effect of double substitution at the binding site can be reasonably approximated by summing the effects of single substitutions. Furthermore, this work provides guidance for the preparation of metal oxides with tailored Lewis basicity.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

A charged particle transport approximation for thick and thin plasmas

In this work, we will construct a simple method to be utilized in the interpretation of the properties of a thin plasma from a study of the transport of charged particles through it. We do so by first demonstrating how charged particle fluxes encode local plasma information in their spectra for sufficiently thick plasmas and then propose how to extend that analysis to the thin plasma limit where the extensive geometry cannot be neglected through arguments from scale separation. We provide a numerical treatment of the transport problem and demonstrate that it is in good agreement with the approximation we present. Finally, we utilize both the numerical solution and our novel approximation to study the impact of the extensive scale of the plasma on the yield and flux normalized high-energy neutron spectra resulting from the upscattering of charged fuel ions. Using this analysis, we show that (after controlling for fusion yield) for the same uniform densities and temperatures, larger plasmas have higher magnitude but softer reaction-in-flight (RIF) neutron spectra relative to smaller plasmas. This is because larger plasmas retain more knocked-on suprathermal ions within their bulk and can downscatter them to lower average energies, while smaller plasmas allow a larger fraction of high energy particles to “range out” of the system prior to substantial downscattering or inducing an RIF reaction.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Mott localization in the van der Waals crystal CrI 3 : A GGA+DMFT study

Using the generalized gradient approximation plus dynamical mean-field theory (GGA+DMFT) we confirm the importance of multi-orbital dynamical correlations in determining the paramagnetic insulating state of CrI 3 . Furthermore, while the ferromagnetic phase reveals weak electronic correlation effects due to strong spin-orbital polarization, the Mott insulating state of paramagnetic CrI 3 crystal is shown to be driven by the interplay between orbital-dependent one-electron lineshape and multi-orbital electronic interactions. To probe the paramagnetic Mott insulating state we performed x-ray absorption spectroscopy (XAS) measurements for the two structural phases of CrI 3 . Our study is relevant to understanding the orbital-selective electronic structure reconstruction of Mott insulators and should be applicable to other van der Waals bonded materials from bulk to the ultrathin limit.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Density functional theory of water with the machine-learned DM21 functional

The delicate interplay between functional-driven and density-driven errors in density functional theory (DFT) has hindered traditional density functional approximations (DFAs) from providing an accurate description of water for over 30 years. Recently, the deep-learned DeepMind 21 (DM21) functional has been shown to overcome the limitations of traditional DFAs as it is free of delocalization error. To determine if DM21 can enable a molecular-level description of the physical properties of aqueous systems within Kohn–Sham DFT, we assess the accuracy of the DM21 functional for neutral, protonated, and deprotonated water clusters. Furthermore, we find that the ability of DM21 to accurately predict the energetics of aqueous clusters varies significantly with cluster size. Additionally, we introduce the many-body MB-DM21 potential derived from DM21 data within the many-body expansion of the energy and use it in simulations of liquid water as a function of temperature at ambient pressure. We find that size-dependent functional-driven errors identified in the analysis of the energetics of small clusters calculated with the DM21 functional result in the MB-DM21 potential systematically overestimating the hydrogen-bond strength and, consequently, predicting a more ice-like local structure of water at room temperature.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Semiclassical theory and the Koopman-van Hove equation *

Abstract The phase space Koopman-van Hove (KvH) equation can be derived from the asymptotic semiclassical analysis of partial differential equations. Semiclassical theory yields the Hamilton–Jacobi equation for the complex phase factor and the transport equation for the amplitude. These two equations can be combined to form a nonlinear semiclassical version of the KvH equation in configuration space. There is a natural injection of configuration space solutions into phase space and a natural projection of phase space solutions onto configuration space. Hence, every solution of the configuration space KvH equation satisfies both the semiclassical phase space KvH equation and the Hamilton–Jacobi constraint. For configuration space solutions, this constraint resolves the paradox that there are two different conserved densities in phase space. For integrable systems, the KvH spectrum is the Cartesian product of a classical and a semiclassical spectrum. If the classical spectrum is eliminated, then, with the correct choice of Jeffreys–Wentzel–Kramers–Brillouin (JWKB) matching conditions, the semiclassical spectrum satisfies the Einstein–Brillouin–Keller quantization conditions which include the correction due to the Maslov index. However, semiclassical analysis uses different choices for boundary conditions, continuity requirements, and the domain of definition. For example, use of the complex JWKB method allows for the treatment of tunneling through the complexification of phase space. Finally, although KvH wavefunctions include the possibility of interference effects, interference is not observable when all observables are approximated as local operators on phase space. Observing interference effects requires consideration of nonlocal operations, e.g. through higher orders in the asymptotic theory.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Investigation of local distortion effects on X-ray absorption of ferroelectric perovskites from first principles simulations

Understanding the role of ferroelectric polarization in modulating the electronic and structural properties of crystals is critical for advancing these materials for overcoming various technological and scientific challenges. However, due to difficulties in performing experimental methods with the required resolution, or in interpreting the results of methods therein, the nanoscale morphology and response of these surfaces to external electric fields has not been properly elaborated. Here, in this work, we investigate the effect of ferroelectric polarization and local distortions in a BaTiO 3 perovskite, using two widely used computational approaches which treat the many-body nature of X-ray excitations using different philosophies, namely the many-body, delta-self-consistent-field determinant (mb-ΔSCF) and the Bethe–Salpeter equation (BSE) approaches. We show that in agreement with our experiments, both approaches consistently predict higher excitations of the main peak in the O–K edge for the surface with upward polarization. However, the mb-ΔSCF approach mostly fails to capture the L 2,3 separations at the Ti–L edge, due to the absence of spin–orbit coupling in Kohn–Sham density functional theory (KS-DFT) at the generalized gradient approximation level. On the other hand, and most promising, we show that application of the GW/BSE approach successfully reproduces the experimental XAS, both the relative peak intensities as well as the L 2,3 separations at the Ti–L edges upon ferroelectric switching. Thus simulated XAS is shown to be a powerful method for capturing the nanoscale structure of complex materials, and we underscore the need for many-body perturbation approaches, with explicit consideration of core-hole and multiplet effects, for capturing the essential physics in these systems.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2025

Libration of hydroxyl groups in layered aluminum (oxy)hydroxides and other material analogs: insights from inelastic neutron scattering and theory

We analyzed the hydroxyl librational signatures of five structurally related aluminum (oxy)hydroxides, using inelastic neutron scattering (INS) and plane-wave lattice dynamics simulations. A clear trend across these aluminum-containing phases illustrates the relationship between hydrogen bonding, local atomic structure, and the spectral location and profile of the librational bands. The INS spectra have been compared to previous optical spectroscopy and computational studies, highlighting the complementary nature of the INS technique. Taking into account other structurally or chemically related material analogs, we have identified a correlation between a blueshift (to higher energy) of the upper librational band edge and the geometry of the hydrogen bond interactions, mirroring (with opposite correlation) the well-known redshift in the intramolecular O–H stretching energy with increasing hydrogen bond strength. For hydroxyl groups that do not participate in hydrogen bonding effectively, the bending librations occur at lower energies and hybridize with metal–oxygen lattice modes. Standard density functional theory approximations, including dispersion corrections, struggle to correctly predict vibrational frequencies of motions dominated by H but perform well for metal–oxygen modes, allowing us to make detailed mode assignments in several cases, including a demonstration of how layer-to-layer disorder in boehmite hydrogen bond orientations is reflected in the sharp but minor low energy peaks (at ∼70–80 meV) of the INS spectrum.

Wang, Hsiu-Wen [Oak Ridge National Laboratory (ORN↗
DOE OSTI2020

An explicit scheme to enforce charge conservation in transient Particle-in-Cell simulations with Maxwell-Boltzmann electrons

The Maxwell-Boltzmann electron model is a very popular approximation of the electronic behavior in electrostatic plasmas, adopted in all cases where the electron dynamics is governed only by a balance between electric and pressure forces. In such a model, the electron dynamics is reduced to solving for a nonlinear Poisson problem together with an additional expression enforcing charge conservation. In this work we derive an expression of charge conservation which can be conveniently applied to explicit schemes to update the reference Boltzmann density, particularly useful when simulating plasma sheaths. The scheme can be equally applied to both steady-state and transient problems. Two examples are shown, a steady-state plasma sheath, and a radio-frequency magnetic presheath. Our proposed scheme allows to enforce global charge conservation locally in time, and can thus be applied to the simulation of transient phenomena.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Spatial correlations of entangled polymer dynamics

The spatial correlations of entangled polymer dynamics are examined by molecular dynamics simulations and neutron spin-echo spectroscopy. Due to the soft nature of topological constraints, the initial spatial decays of intermediate scattering functions of entangled chains are, to the first approximation, surprisingly similar to those of an unentangled system in the functional forms. However, entanglements reveal themselves as a long tail in the reciprocal-space correlations, implying a weak but persistent dynamic localization in real space. Furthermore, comparison with a number of existing theoretical models of entangled polymers suggests that they cannot fully describe the spatial correlations revealed by simulations and experiments. In particular, the strict one-dimensional diffusion idea of the original tube model is shown to be flawed. The dynamic spatial correlation analysis demonstrated in this work provides a useful tool for interrogating the dynamics of entangled polymers. Lastly, the failure of the investigated models to even qualitatively predict the spatial correlations of collective single-chain density fluctuations points to a possible critical role of incompressibility in polymer melt dynamics.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2022

Dynamic aspects of graphene deformation and fracture from approximate density functional theory

Graphene is one of the most intriguing two-dimensional carbon materials. Its mechanical strength and failure are key concerns for materials engineering and applications. Despite the success of fracture mechanics, the mechanism of how pristine materials fail remains an elusive problem. While many theoretical studies based on molecular dynamics using empirical forcefields have tried to address this question, atomic-scale mechanics are not clearly understood. Especially, a widely employed bond-breaking approach based on the critical bond length has not been rigorously tested. Here, utilizing molecular dynamics simulations with density functional based tight binding, we investigate how the failure of the pristine material initiates. The Wiberg bond order (W BO ) to estimate the change of chemical bonds shows a transition from sp 2 (W BO ~ 1.33) to sp 3 (WBO < 1.0) during the deformation. However, it reveals that a single threshold value for either the WBO or bond length is insufficient to decide failure of pristine material without free edges or defects. Instead, collective behaviors of the local atomic group govern the fracture initiation of pristine graphene. Overall, our study provides dynamic mechanical responses based on quantum mechanics, which have not been captured by empirical forcefields, opening opportunities to design properties by precisely coupling the mechanics and quantum mechanics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Theory of Plasmonic Hot-Carrier Generation and Relaxation

Hot-carrier (HC) generation from (localized) surface plasmon decay has recently attracted much attention due to its promising applications in physical, chemical, materials, and energy science. However, the detailed mechanisms of plasmonic HC generation, relaxation, and trapping are less studied. In this work, we developed and applied a quantum-mechanical model and coupled master equation method to study the generation of HCs from plasmon decay and their following relaxation processes with different mechanisms treated on equal footing. First, a quantum-mechanical model for HC generation is developed. Its connection to existing semiclassical models and time-dependent density functional theory (TDDFT) is discussed. Second, the relaxation and lifetimes of HCs are investigated in the presence of electron–electron and electron–phonon interactions. A GW-like approximation is introduced to account for the electron–electron scattering. The numerical simulations on the Jellium nanoparticles with a size up to 1.6 nm demonstrate the electron–electron scattering and electron–phonon scattering dominate different time scale in the relaxation dynamics. We also generalize the model to study the extraction of HCs to attached molecules. The quantum yield of extracting HCs for other applications is found to be size-dependent. In general, the smaller size of NP improves the quantum yield, which is in agreement with recent experimental measurements. Even though we demonstrate this newly developed theoretical formalism with Jellium model, the theory applies to any other atomistic models.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Cosmological quasiparticles and the cosmological collider

The interplay between cosmology and strongly coupled dynamics can yield transient spectral features that vanish at late times, but which may leave behind phenomenological signatures in the spectrum of primordial fluctuations. Of particular interest are strongly coupled extensions of the standard model featuring approximate conformal invariance. In flat space, the spectral density for a scalar operator in a conformal field theory is characterized by a continuum with scaling law governed by the dimension of the operator, and is otherwise featureless. Anti–de Sitter/conformal field theory arguments suggest that for large 𝑁, in an inflationary background with Hubble rate 𝐻, this continuum is gapped. We demonstrate that there can be additional peak structures that become sharp and particlelike at phenomenologically interesting regions in parameter space, and we estimate their contribution to cosmological observables. We find phenomena that are potentially observable in future experiments that are unique to these models, including displaced oscillatory features in the squeezed limit of the bispectrum. These particles can be either fundamental, and localized to a UV brane, or composite at the Hubble scale, 𝐻, and bound to a horizon in the bulk of the 5D geometry. We comment on how stabilization of conformal symmetry breaking vacua can be correlated with these spectral features and their phenomenology.

conformal field theory↗
DOE OSTI2023

A New Plasma Radar Concept for Simultaneous Magnetic and Density Measurements

An innovative, compact 288GHz interferometer has been fabricated, tested, installed and successfully demonstrated on the LAPD-U magnetized plasma at UCLA. The system takes advantage of frequency modulated (FM) radar techniques to deliver a compact heterodyne system. In addition, the reflected power from the source is taken advantage of to eliminate the need for additional quasi-optical components. Electron density in LAPD-U plasma has recently been increased substantially thereby requiring a higher frequency/shorter wavelength interferometer to avoid deleterious refractive effects. This system satisfies those needs. The system uses a 96GHz varactor tuned Gunn oscillator which passes to a passive tripler. This tripler has ~3% conversion efficiency. The 288GHz radiation is then coupled to free space using a so-called dual-mode or Pickett horn. The output 288GHz beam is then coupled to an aspheric lens manufactured from low-loss, high-density polyethylene. This lens is employed to collimate the emerging beam. Small axial adjustment of the lens position can also be used to create a slowly focusing beam so as to optimize the measured signal. In addition, up-down or side-to-side adjustment of the lens can be utilized to steer the beam vertically or horizontally – again to optimize alignment. The propagating beam passes through a beam splitter and then through a water-free, bubble-free fused quartz window into the LAPD-U vacuum vessel. The beam-splitter is a thin sheet of G10 which reflects a small fraction of the incident power (~5 %) towards a zero-bias detector optimized for the frequency range from 220 to 300GHz. Note that waveguides at this frequency have dimensions of ~0.9mm x 0.45mm and so have very large conductive losses. This drives the use of quasi-optical propagation. The detector requires no DC bias and is very responsive (> 1V/mW into 1MΩ). Radiation is coupled to the detector via a similar lens-horn arrangement used for the launch. This reflected beam acts as the local oscillator or reference millimeter-wave beam for the detector. The remainder of the launched source beam then enters the LAPD-U vacuum vessel and passes through the plasma at the mid-plane until reaching the opposing port which is closed off with an aluminum flange. This flange is used as a mirror to retroreflect the incident 288GHz beam back along its path. The retroreflected beam exits the input port but does NOT couple directly into the zero-bias detector. Instead, the majority of the return power continues towards the 288GHz source. As mentioned above the transmitted beam enters the source a second time. This would appear undesirable. However, at these frequencies multipliers are highly non-linear elements which results in a significant portion of the return beam (~20%) re-emerging from the multiplier and horn and then coupling via the G10 beam-splitter to the zero-bias detector. This approach eliminated the need for a second quasi-optical beam-splitter. The system is extremely compact measuring approximately 28 inches x 20 inches. The above did not explain how heterodyne operation was achieved. As mentioned above the Gunn oscillator is able to be varactor tuned. This allows a low voltage to be applied to control the operating frequency of the Gunn oscillator. During heterodyne operation a sawtooth shaped voltage is applied to the varactor at 750kHz using an 80MHz Arbitrary Waveform Generator (AWG). This voltage changes the Gunn frequency linearly during the up-sweep which is then reset abruptly at the sawtooth crash to be immediately followed by another linear sweep. Passage through the 288GHz multiplier triples the frequency change experienced by the electromagnetic wave. These frequency changes are small – tens of megahertz. This FM radar approach results in the launched electromagnetic wave frequencies at the detector for the reference and plasma wave to be different. The approximately 10 ns delay propagation delay for the plasma beam results in the local oscillator and plasma beams NOT having an identical frequency – there is in fact a fixed difference frequency. The frequency tuning level of the Gunn oscillator is then adjusted so that there is ONE cycle of this difference frequency during each linear ramp. During the sawtooth crash or downward re-sweep this one cycle replays in reverse but on a very fast timescale. The process then repeats. Low-pass filtering eliminates the fast re-sweep to leave a pure sine wave heterodyne signal. When the plasma is present it introduces a phase delay in the sine wave (caused by the extremely small Doppler shift resulting from the optical path length change). Of course, to measure this phase change we need a reference. This is simply obtained from an arbitrary waveform generator which provides a synchronized output pulse train which again is low pass filtered to obtain a 750kHz sinusoidal voltage reference for the interferometer. The interferometer was installed on LAPD-U where it has worked reliably and has established that electron densities exceeding 1x10 13 cm -3 are routinely achieved. In addition, the system sensitivity was able to easily observe density fluctuation at frequencies up to 50kHz. FM Radar techniques have enabled a full demonstration of a compact, sensitive, high frequency (288GHz/1mm) heterodyne interferometer.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

How Many Bursts Does It Take to Form a Core at the Center of a Galaxy?

We present a novel method for systematically assessing the impact of the central potential fluctuations associated with bursty outflows on the structures of dark matter halos for classical and ultrafaint dwarf (UFD) galaxies. Specifically, we use dark-matter-only simulations augmented with a manually added massive particle that modifies the central potential and approximately accounts for a centrally concentrated baryonic component. This approach enables precise control over the magnitude, frequency, and timing of rapid outflow events. We demonstrate that this method can reproduce the established result of core formation for systems that undergo multiple episodes of bursty outflows. In contrast, we also find that equivalent models involving only single (or a small number of) burst episodes do not form cores with the same efficacy. This is important because many UFDs in the local Universe are observed to have tightly constrained star formation histories that are best described by a single early burst of star formation. Using a suite of cosmological zoom-in simulations, we identify the regimes in which single bursts can and cannot form a cored density profile. Our results suggest that it may be difficult to form cores in UFD-mass systems with a single early burst, regardless of its magnitude.

cold dark matter↗
DOE OSTI2025

Simulations of Quantum Approximate Optimization Algorithm on HPC-QC Integrated Systems

The Quantum Approximate Optimization Algorithm (QAOA) has emerged as a promising tool for accelerating optimization processes in the Noisy Intermediate-Scale Quantum (NISQ) era. Compared to classical methods, QAOA efficiently solves optimization problems, often formulated as Quadratic Unconstrained Binary Optimization (QUBO) problems. Classical quantum simulators are crucial for evaluating quantum algorithms due to limited quantum resources. However, QAOA's performance can vary with different simulation methods. This study analyzes QAOA's performance using various quantum simulators (e.g., density _matrix, statevector, and matrix_product_state) and demonstrates the benefits of HPC-QC integrated systems in solving QUBO problems on an active learning workflow. By simulating QAOA on dense, large-matrix QUBO problems, we evaluate accuracy and problem-solving time. We also assess QAOA's performance on local computers and HPC-QC inte-grated systems, using Oak Ridge Leadership Computing Facility (OLCF)'s Frontier supercomputer with local Qiskit Aer and remote IBM Quantum simulators.

Kim, Seongmin [ORNL] (ORCID:0000000159063004)↗
DOE OSTI2022

Applying Thomson scattering to diagnosing turbulent density and velocity fluctuations in a gas-puff z-pinch

The electron plasma wave feature (EPW) in the time-resolved Thomson scattering spectrum is used to obtain the local electron density in imploding high energy density gas-puff z-pinch plasmas. The optical setup was optimized to allow the relatively weak EPW feature obtained from 1 MA imploding neon gas-puff z-pinches to be seen above the continuum emission as well as the brighter ion acoustic wave (IAW) feature. Using a frequency-doubled Nd:YLF laser (E = 10 J, λ = 526.5 nm, Δt= 2.3 ns, spot size ~ 250 μm) and two visible light streak cameras, we determined the average electron density in the imploding plasma sheath 40 ns prior to stagnation to be n e = 2.5 × 10 18 /cm 3 . At pinch time, it reached n e = 1.7 × 10 19 /cm 3 . Here, the electron temperature during implosion measured via the IAW (approximately 50 eV) was four times lower than the implosion electron temperature measured via the EPW (approximately 200 eV), assuming that neither feature is affected by turbulent fluctuations in the plasma. In order for the electron temperatures inferred from the EPW and IAW spectral features to be self-consistent, we find that it is necessary to include velocity fluctuations in the analysis of the IAW feature peaks and corresponding density fluctuation in the peak widths of the EPW feature.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

FINETUNA: fine-tuning accelerated molecular simulations

Abstract Progress towards the energy breakthroughs needed to combat climate change can be significantly accelerated through the efficient simulation of atomistic systems. However, simulation techniques based on first principles, such as density functional theory (DFT), are limited in their practical use due to their high computational expense. Machine learning approaches have the potential to approximate DFT in a computationally efficient manner, which could dramatically increase the impact of computational simulations on real-world problems. However, they are limited by their accuracy and the cost of generating labeled data. Here, we present an online active learning framework for accelerating the simulation of atomic systems efficiently and accurately by incorporating prior physical information learned by large-scale pre-trained graph neural network models from the Open Catalyst Project. Accelerating these simulations enables useful data to be generated more cheaply, allowing better models to be trained and more atomistic systems to be screened. We also present a method of comparing local optimization techniques on the basis of both their speed and accuracy. Experiments on 30 benchmark adsorbate-catalyst systems show that our method of transfer learning to incorporate prior information from pre-trained models accelerates simulations by reducing the number of DFT calculations by 91%, while meeting an accuracy threshold of 0.02 eV 93% of the time. Finally, we demonstrate a technique for leveraging the interactive functionality built in to Vienna ab initio Simulation Package (VASP) to efficiently compute single point calculations within our online active learning framework without the significant startup costs. This allows VASP to work in tandem with our framework while requiring 75% fewer self-consistent cycles than conventional single point calculations. The online active learning implementation, and examples using the VASP interactive code, are available in the open source FINETUNA package on Github.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Dynamical cluster approximation study of electron localization in the extended Hubbard model

We perform a detailed study of the phase transitions and mechanisms of electron localization in the extended Hubbard model using the dynamical cluster approximation on a 2×2 cluster. We explore the interplay of charge order and Mott physics. We find that a nearest-neighbor Coulomb interaction V causes “screening” effects close to the Mott phase transition, pushing the phase boundary to larger values of U. We also demonstrate the different effects of V on correlations in metallic and insulating regimes, and document the different correlation aspects of charge order and Mott states.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗