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.
We discuss dynamical enhancement of entanglement in a driven Bose-Hubbard model and find an enhancement of two orders of magnitude from the ground state value which is robust against fluctuations in experimental parameters.
Despite being of fundamental importance and potential interest for topological quantum computing, spin-triplet superconductors remain rare in solid state materials after decades of research. In this work, we present a three-particle mechanism for spin-triplet superconductivity in multiband systems, where an effective attraction between doped electrons is produced from the Coulomb repulsion via a virtual interband transition involving a third electron [V. Crépel, L. Fu, Sci. Adv. 7, eabh2233 (2021)]. Our theory is analytically controlled by an interband hybridization parameter and explicitly demonstrated in doped band insulators with the example of an extended Hubbard model. Our theory of exciton-mediated pairing reveals how, as a matter of principle, a two-particle bound state can arise from the strong electron repulsion upon doping, opening a viable path to Bose–Einstein condensate (BEC)–Bardeen–Cooper–Schrieffer (BCS) physics in solid state systems. In light of this theory, we propose that recently discovered dilute superconductors such as ZrNCl, WTe 2 , and moiré materials can be spin-triplet and compare the expected consequences of our theory with experimental data.
Motivated by the discovery of superconductivity in the two-leg, quasi-one dimensional ladder compound, BaFe 2 S 3 we present a renormalization group study of electrons moving on a two leg, two orbital ladder, subjected to Hubbard repulsion U and Hund’s coupling J . In our calculations, we adopt tightbinding parameters obtained from ab-initio studies on this material. At incommensurate filling, the long wavelength analysis displays four phases as a function of 0 ≤ J/U < 1. We show that a fully gapped superconductor is stabilized at sufficiently large Hund’s coupling, the relative phases at the three Fermi points are “+,-,-”. By contrast, when the system is tuned to half filling, Umklapp scattering gives rise to Mott insulating phases. Furthermore, we discuss the general implications of our study for the broad class of iron-based superconductors.
The two-loop renormalization group (RG) calculation is considerably extended here for the two-dimensional (2D) fermionic effective field theory model, which includes only the so-called “hot spots” that are connected by the spin-density-wave (SDW) ordering wavevector on a Fermi surface generated by the 2D t−t{sup ′} Hubbard model at low hole doping. We compute the Callan–Symanzik RG equation up to two loops describing the flow of the single-particle Green’s function, the corresponding spectral function, the Fermi velocity, and some of the most important order-parameter susceptibilities in the model at lower energies. As a result, we establish that–in addition to clearly dominant SDW correlations–an approximate (pseudospin) symmetry relating a short-range incommensurated-wave charge order to the d-wave superconducting order indeed emerges at lower energy scales, which is in agreement with recent works available in the literature addressing the 2D spin-fermion model. We derive implications of this possible electronic phase in the ongoing attempt to describe the phenomenology of the pseudogap regime in underdoped cuprates.
Strongly correlated electron systems host a variety of poorly understood correlations in their high-temperature normal state. Unlike ordered phases defined by order parameters, regions of the normal state are often defined through unconventional properties such as strange metallic transport or spectroscopic pseudogaps. Characterizing the microscopic correlations in the normal state is necessary to elucidate mechanisms that lead to these properties and their connection to ground-state orders. Here, in this paper, we establish the presence of intertwined charge and spin stripes in the strange metal normal state of the Hubbard model using determinant quantum Monte Carlo calculations. The charge and spin density waves constituting the stripes are fluctuating and short ranged; yet they obey a mutual commensurability relation and remain microscopically interlocked, as evidenced through measurements of three-point spin-spin-hole correlation functions. Our findings demonstrate the ability of many-body numerical simulations to unravel the microscopic correlations that define quantum states of matter.
Abstract not provided.
We propose a way of obtaining effective low energy Hubbard-like model Hamiltonians from ab initio quantum Monte Carlo calculations for molecular and extended systems. The Hamiltonian parameters are fit to best match the ab initio two-body density matrices and energies of the ground and excited states, and thus we refer to the method as ab initio density matrix based downfolding. For benzene (a finite system), we find good agreement with experimentally available energy gaps without using any experimental inputs. For graphene, a two dimensional solid (extended system) with periodic boundary conditions, we find the effective on-site Hubbard U{sup ∗}/t to be 1.3 ± 0.2, comparable to a recent estimate based on the constrained random phase approximation. For molecules, such parameterizations enable calculation of excited states that are usually not accessible within ground state approaches. For solids, the effective Hamiltonian enables large-scale calculations using techniques designed for lattice models.
We have evaluated the successes and failures of the Hubbard-corrected density functional theory approach to study Mg doping of LiCoO{sub 2}. We computed the effect of the U parameter on the energetic, geometric, and electronic properties of two possible doping mechanisms: (1) substitution of Mg onto a Co (or Li) site with an associated impurity state and (2) formation of impurity-state-free complexes of substitutional Mg and point defects in LiCoO{sub 2}. We find that formation of impurity states results in changes on the valency of Co in LiCoO{sub 2}. Variation of the Co U shifts the energy of the impurity state, resulting in energetic, geometric, and electronic properties that depend significantly on the specific value of U. In contrast, the properties of the impurity-state-free complexes are insensitive to U. These results identify reasons for the strong dependence on the doping properties on the chosen value of U and for the overall difficulty of achieving agreement with the experimentally known energetic and electronic properties of doped transition metal oxides such as LiCoO{sub 2}.
We consider two-component one-dimensional quantum gases at special imbalanced commensurabilities which lead to the formation of multimer (multiparticle bound-states) as the dominant order parameter. Luttinger liquid theory supports a mode-locking mechanism in which mass (or velocity) asymmetry is identified as the key ingredient to stabilize such states. While the scenario is valid both in the continuum and on a lattice, the effects of umklapp terms relevant for densities commensurate with the lattice spacing are also mentioned. These ideas are illustrated and confronted with the physics of the asymmetric (mass-imbalanced) fermionic Hubbard model with attractive interactions and densities such that a trimer phase can be stabilized. Phase diagrams are computed using density-matrix renormalization group techniques, showing the important role of the total density in achieving the latter phase. The effective physics of the trimer gas is studied as well. Last, the effect of a parabolic confinement and the emergence of a crystal phase of trimers are briefly addressed. This model has connections with the physics of imbalanced two-component fermionic gases and Bose-Fermi mixtures as the latter gives a good phenomenological description of the numerics in the strong-coupling regime.
Simulating the Hubbard model is of great interest to a wide range of applications within condensed matter physics, however its solution on classical computers remains challenging in dimensions larger than one. The relative simplicity of this model, embodied by the sparseness of the Hamiltonian matrix, allows for its efficient implementation on quantum computers, and for its approximate solution using variational algorithms such as the variational quantum eigensolver. While these algorithms have been shown to reproduce the qualitative features of the Hubbard model, their quantitative accuracy in terms of producing true ground state energies and other properties, and the dependence of this accuracy on the system size and interaction strength, the choice of variational ansatz, and the degree of spatial inhomogeneity in the model, remains unknown. Here we present a rigorous classical benchmarking study, demonstrating the potential impact of these factors on the accuracy of the variational solution of the Hubbard model on quantum hardware, for systems with up to 32 qubits. We find that even when using the most accurate wavefunction ansätze for the Hubbard model, the error in its ground state energy and wavefunction plateaus for larger lattices, while stronger electronic correlations magnify this issue. Concurrently, spatially inhomogeneous parameters and the presence of off-site Coulomb interactions only have a small effect on the accuracy of the computed ground state energies. Our study highlights the capabilities and limitations of current approaches for solving the Hubbard model on quantum hardware, and we discuss potential future avenues of research.
Neural quantum states (NQSs) have emerged as a powerful ansatz for variational quantum Monte Carlo studies of strongly correlated systems. Here, we apply recurrent neural networks (RNNs) and autoregressive transformer neural networks to the Fermi-Hubbard and the (non-Hermitian) Hatano-Nelson-Hubbard models in one and two dimensions. In both cases, we observe that the convergence of the RNN ansatz is challenged when increasing the interaction strength. We present a physically motivated and easy-to-implement strategy for improving the optimization, namely, by ramping of the model parameters. Furthermore, we investigate the advantages and disadvantages of the autoregressive sampling property of both network architectures. For the Hatano-Nelson-Hubbard model, we identify convergence issues that stem from the autoregressive sampling scheme in combination with the non-Hermitian nature of the model. Our findings provide insights into the challenges of the NQS approach and make the first step towards exploring strongly correlated electrons using this ansatz. Published by the American Physical Society 2025
Using a separable many-body variational wavefunction, we formulate a self-consistent effective Hamiltonian theory for fermionic many-body system. The theory is applied to the two-dimensional (2D) Hubbard model as an example to demonstrate its capability and computational effectiveness. Most remarkably for the Hubbard model in 2D, a highly unconventional quadruple-fermion non-Cooper pair order parameter is discovered.
Here, we assess the performance of seven exchange–correlation functionals (some with and some without a Hubbard U correction) for their ability (i) to predict band gaps of silicon, diamond, and Li-ion battery cathode materials, (ii) to localize hole polarons and predict delithiation energies in Li-ion battery cathode materials, and (iii) to predict transition levels of charge carriers of doped silicon and diamond. Both local and hybrid exchange–correlation functionals were tested. The local functionals tend to underestimate band gaps and delocalize polarons. The hybrid functionals very often give a good description of both properties, but they may not be practical for calculations involving large unit cells, large ensembles, or dynamics, and therefore a local functional with a Hubbard U correction is often used (giving the method called DFT+U), where the value of a parameter U is adjusted according to the system and the property being investigated. Keeping in mind the importance of computational cost and the undesirability of having to adjust an empirical parameter for each system or property of interest, we recently developed a local functional, namely HLE17, to try to accurately predict band gaps and excitation energies, and we validated it using mostly main-group solids and molecules. Here we test the performance of HLE17 for its ability to predict band gaps and localize polarons in other solid-state materials, and we compare its performance to that of popular local functionals (PBE, PBEsol, and TPSS), a range-separated hybrid functional with screened exchange (HSE06), and DFT+U. We find that HLE17 predicts more accurate band gaps than other local functionals and can localize holes as polarons, which other local functionals usually fail to do, and for a number of cases it is comparable in performance quality to Hubbard-corrected functionals without the need for system-specific parametrization and to hybrid functionals without the high cost. Because HLE17 does not predict accurate lattice constants, we use the single-point method of quantum chemistry, where the geometry is optimized with one functional and the band gap is calculated with HLE17, or we perform calculations with the lattice constants obtained by TPSS and both the fractional intracell coordinates and the electronic structure obtained by HLE17 (a new method denoted HLE17\\TPSS). In particular, we performed calculations by HLE17//TPSS, HLE17//HSE06, HLE17//DFT+U, and HLE17\\TPSS, and these methods usually agree well with each other and give values similar to experiment.
We apply perturbation theory and cyclic spin permutation formalism to study the lowest energy states of the infinite-repulsion Hubbard model on finite fragments of n-leg ladder rings and show jump-wise behavior of the ground state spin S{sub 0} as a function of fragment parameters.
Among the various numerical techniques to study the physics of strongly correlated quantum many-body systems, the self-energy functional approach (SFA) has become increasingly important. In its previous form, however, SFA is not applicable to Bose-Einstein condensation or superfluidity. In this paper, we show how to overcome this shortcoming. To this end, we identify an appropriate quantity, which we term D, that represents the correlation correction of the condensate order parameter, as it does the self-energy for Green's function. An appropriate functional is derived, which is stationary at the exact physical realization of D and of the self-energy. Its derivation is based on a functional-integral representation of the grand potential followed by an appropriate sequence of Legendre transformations. The approach is not perturbative and, therefore, applicable to a wide range of models with local interactions. We show that the variational cluster approach based on the extended self-energy functional is equivalent to the ''pseudoparticle'' approach proposed in Phys. Rev. B 83, 134507 (2011). We present results for the superfluid density in the two-dimensional Bose-Hubbard model, which shows a remarkable agreement with those of quantum-Monte-Carlo calculations.
In this approach, we have employed the ab initio method of full-potential linearized plane waves plus local orbital based on the framework of density functional theory, within the spin generalized gradient approximation plus U (GGA+U), where U is the Hubbard correlation term for exchange-correlation potential. This method is used to predict the equilibrium structural parameters, the electronic structure, and the magnetic properties of Rh{sub 2}FeGa and Rh{sub 2}FeIn full-Heusler alloys. The structural properties result show that both Rh{sub 2}FeGa and Rh{sub 2}FeIn compounds are stable in ferromagnetic states. The spin-polarized electronic structure demonstrate that both full-Heusler compounds are half-metals with half-metallic band gap (E{sub HM}) of 0.326 eV for Rh{sub 2}FeGa alloy and 0.245 eV for Rh{sub 2}FeIn alloy, where the apparition of an indirect gap from Γ to X in minority spin is shown by these Heusler alloys. Furthermore, calculations of total magnetic moment give the value of 5 µ{sub B}, which agrees with Slater-Pauling rule; also, the atomic magnetic moment of both Rh and Fe transition metals have the same sign which depicts the ferromagnetic behavior.
In this work, we employ the fully anisotropic DFT+U+J approach with the PBEsol functional to investigate ground-state magnetic and electronic properties of bulk binary manganese oxides: MnO, Mn 3 O 4 , α-Mn 2 O 3 , and β-MnO 2 , in order of increasing Mn valence. The computed crystal structures, noncollinear magnetic ground states, and corresponding electronic structures are in good agreement with the experimental data and hybrid functional calculations available in the literature. We take into account the nonlinear core-valence interaction in our Mn pseudopotential designed by ourselves, as it has been proven to be important for transition-metal systems. Although the Hubbard U term is capable by itself of opening a band gap, the explicitly defined exchange parameter J plays an important role in improving the detailed electronic and noncollinear magnetic structure profiles. Appropriate band gaps are obtained with U values smaller than those used in previously reported calculations. Our results suggest that pseudopotential design together with DFT+U+J enables the acquisition of accurate properties of complex magnetic systems using a nonhybrid density functional.
In this work, we consider a multiband Hubbard model H m for Cu and O orbitals in Ba 2-x Sr x CuO 3+δ similar to the three-band model for two-dimensional cuprates. The hopping parameters are obtained from maximally localized Wannier functions derived from ab initio calculations. Using the cell perturbation method, we derive both a generalized t–J model H tJ and a one-band Hubbard model H H to describe the low-energy physics of the system. H tJ has the advantage of having a smaller relevant Hilbert space, facilitating numerical calculations, while additional terms should be included in H H to accurately describe the multiband physics of H m . Using H tJ and the density matrix renormalization group method, we calculate the wave-vector-resolved photoemission and discuss the relevant features in comparison with recent experiments. In agreement with previous calculations, we find that the addition of an attractive nearest-neighbor interaction of the order of the nearest-neighbor hopping shifts the weight from the 3k F to the holon-folding branch. Kinetic effects also contribute to this process.