Search NASA⌕ Search

SEARCH · Search NASA

Results for “Zero point energy”

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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

First Solar Thermal Energy Planner (STEP 1) and Nationwide Industrial Heat and Power Analysis

The First Solar Thermal Energy Planner (STEP 1) and Nationwide Industrial Heat and Power Analysis (aka the STEP 1 Project or the Project) aimed to (1) developed a brand-new web tool that could provide decision support through free, rapid techno-economic analysis of behind-the-meter solar+storage systems for industrial process heat and (2) conduct high-level analyses of the cost-competitiveness of the same systems across the US in key sectors. The STEP 1 web tool collects key location, land availability, thermal load profile, proccess heat temperature and media, and other key parameters through an easy-to-use user interface (UI). The UI was designed to meet the user at their level of understanding by minimizing the number of required inputs as much as possible while including options for more nuanced inputs if the user desires. STEP 1 advises users on which solar thermal tehcnologies that could fit their needs based on the inputs provided (namely process media and temperature). The tool can model a wide range of solar thermal technologies including flat plate collectors, evacuated tubes, parabolic troughs, linear Fresnel, and molten salt towers all with corresponding thermal energy storage (TES) - solar PV with resistive heating and TES is also included. Once the parameters are collected, a nominal thermal energy production profile for the facility's location is generated using NREL's System Advisory Model (SAM) and then passed to a modified version of NREL's REopt platform to optimize the size (capacity) and dispatch of the solar+storage system to minimize lifecycle costs subject to energy balance, fuel and electricity rates, emission reductions goals, and other constraints. This entire process takes less than 20 minutes, is completely free, requires zero coding skills, and provides the user with a high-level assessment on the techno-economic feasability of deploying solar+storage systems for their energy needs. In addition to the development of the STEP 1 web tool, the project completed two complementary analyses focused on leveraging the backend code of STEP 1, public industrial facility locations and fuel consumption data, and key sector information to assess the economic opportunity of reducing fuel costs at various levels of capacity factor around the US. Both analyses found that there are key markets, locations, industrial sectors (namely food and beverage), and degrees of offset where solar thermal technologies could be cost-effectively deployed, highlighting a key market entry point for these technologies, and assessing deployment potential. In summary, the STEP 1 project improved the opportunities for solar+storage systems to expand into the industrial process heat market through breaking down barriers to assessing these technologies.

14 SOLAR ENERGY↗
DOE OSTI2021

Properties of giant dipole resonances within an extended pairing model with a focus on spectral statistics

In this paper, we report on a study of the spectral features associated with dipole resonances in medium mass nuclei (E ≤ 12 MeV), as revealed in the framework of the spd-interacting boson model. The effect of pairing correlations on the theory follows from solutions obtained through an application of the Bethe Ansatz Equation. In general, calculated spectra around the critical point of the vibrational to γ-soft transitions appears to approach that of a Gaussian Orthogonal Ensemble, while near the rotational and vibrational limits of the theory the spectra show more regular behavior. Specifically, the results reveal that the statistical features of the spectra are sensitive to the vector boson pairing strength, c p , in the transition region; that is, when c p is zero, or when the system approaches one of its dynamical symmetries limits the spectrum display regular features, while for stronger c p values, or when near to the critical phase transition region, the spectral feature show more chaotic behavior. Overall, our results indicate that the statistical features are governed by the interplay between dipole resonant energies, pairing correlations, and interactions between and among the single and vector bosons modes of the theory. As part of this work we also found out that chaoticity occurs when results were fit to a Berry-Robnik distribution. Throughout our analyses, we used experimentally known information about both positive and negative parity states. Our findings suggest that dipole resonances appear to be best-described by Poisson statistics for A ≈ 32-138 nuclei.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Value of Geothermal Energy Storage for Supply-Side and Demand-Side Applications

This report presents the results of a study examining the value potential for geothermal energy storage (GES), a long-duration energy storage resource that stores thermal and/or geomechanical energy in the subsurface. GES could benefit the overall U.S. power system by temporally shifting electricity generation (supply-side) or meeting building heating and cooling load (demand-side). This report analyzes supply-side and demand-side opportunities independently because of differences in applications and models. Currently there is significant uncertainty about the development costs for GES, with only a limited number of demonstration plants for electric energy storage and building heating and cooling storage developments. In this report, we estimate the value of supply-side and demand-side GES to the bulk power system in the contiguous United States. Because of the significant uncertainty about GES development costs, this analysis does not consider GES deployment costs but instead focuses on the value of GES to the U.S. electricity system. The estimated values of GES provide reference points for economically competitive commercial cost targets. Supply-side GES is modeled as part of an enhanced geothermal system (EGS) generation plant in NREL's Regional Energy Deployment System (ReEDS) capacity expansion model (Ho et al. 2021). In contrast to conventional geothermal plants, which generate constant power, EGS plants have unique features that may allow for in-reservoir energy storage for flexible generation. Demand-side GES for heating and cooling, including seasonal hot and cold storage and short-duration heat pump storage, is incorporated into a price-taker model using Cambium electricity marginal cost projections. To establish an upper bound for the value of GES, analysis focused on favorable scenarios for storage with high generation from zero marginal cost, variable renewable energy resources. High penetrations of variable renewable energy generation can increase hourly electricity price variability, which increases the value of temporal energy arbitrage for storage technologies like GES.

15 GEOTHERMAL ENERGY↗
DOE OSTI2022

Embedding a critical point in a hadron to quark-gluon crossover equation of state

Lattice QCD simulations have shown unequivocally that the transition from hadrons to quarks and gluons is a crossover when the baryon chemical potential is zero or small. Many model calculations predict the existence of a critical point at a value of the chemical potential where current lattice simulations are unreliable. We show how to embed a critical point in a smooth background equation of state so as to yield the critical exponents and critical amplitude ratios expected of a transition in the same universality class as the liquid-gas phase transition and the three-dimensional Ising model. There are only two independent critical exponents; the relations α + 2β + γ = 2 and β(δ-1) = γ arise automatically, as does a relation between the two critical amplitudes. The resulting equation of state has parameters that may be inferred by hydrodynamic modeling of heavy-ion collisions in the Beam Energy Scan II at the BNL Relativistic Heavy Ion Collider or in experiments at other accelerators.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Trace anomaly form factors from lattice QCD

The hadron mass can be obtained through the calculation of the trace of the energy-momentum tensor in the hadron which includes the trace anomaly and sigma terms. The anomaly due to conformal symmetry breaking is believed to be an important ingredient for hadron mass generation and confinement. In this work, we will present the calculation of the glue part of the trace anomaly form factors of the pion up to Q 2 ∼ 4.3 GeV 2 and the nucleon up to Q 2 ∼ 1 GeV 2 . The calculations are performed on a domain wall fermion ensemble with overlap valence quarks at seven valence pion masses varying from ∼ 250 to ∼ 540 MeV , including the unitary point ∼ 340 MeV . We calculate the radius of the glue trace anomaly for the pion and the nucleon from the z expansion. By performing a two-dimensional Fourier transform on the glue trace anomaly form factors in the infinite momentum frame with no energy transfer, we also obtain their spatial distributions for several valence quark masses. The results are qualitatively extrapolated to the physical valence pion mass with systematic errors from the unphysical sea quark mass, discretization effects in the renormalization sum rule, and finite-volume effects to be addressed in the future. We find the pion’s form factor changes sign, as does its spatial distribution, for light quark masses. This explains how the trace anomaly contribution to the pion mass approaches zero toward the chiral limit. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Superconductivity from Luttinger surfaces: Emergent Sachdev-Ye-Kitaev physics with infinite-body interactions

The pairing of two electrons on a Fermi surface due to an infinitesimal attraction between them always results in a superconducting instability at zero temperature (T = 0). The equivalent question of pairing instability on a Luttinger surface (LS)—190a contour of zeros of the propagator—instead leads to a quantum critical point (QCP) that separates a non-Fermi liquid and superconductor. A surprising and little understood aspect of pair fluctuations at this QCP is that their thermodynamics maps to that of the Sachdev-Ye-Kitaev (SYK) model in the strong coupling limit. Here, we offer a simple justification for this mapping by demonstrating that (i) LS models share the reparametrization symmetry of the q → ∞ SYK model with q-body interactions close to the LS, and (ii) the enforcement of gauge invariance results in a $\frac{1}{√τ}$ (τ ~ T -1 ) behavior of the fluctuation propagator near the QCP, as is the feature of SYK conformal Green functions, but leaves the overall form of the free energy map robust.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

A hybrid analytical and numerical model for cross-over and performance decay in a unit cell vanadium redox flow battery

Developing an accurate and efficient model for cross-over is critical for improving the long-term performance of redox flow batteries (RFBs). Here, this work presents a hybrid analytical and numerical model that combines a two-dimensional analytical solution to the active species, a one-dimensional analytical model for cross-over mechanisms, and a zero-dimensional numerical model for outlet concentrations of reactants. By comparing with experiment over 41 cycles (ca. 144 h), the model reported a mean voltage difference of 0.0089 V, mean time difference of 14 s per cycle (of 3.5 h), maximum relative difference for capacity and energy of 1.34% and 1.63%, respectively. The predicted mean concentrations for V 2+ , V 3+ , and VO$_{2}^{+}$ in membrane are 65%~77% of measured values from the literature. Upon validation, the model reproduced behaviors in electrolyte imbalance similar to those observed in experiments and numerical models, and revealed the control of cross-over, self-discharge, stoichiometry of side reactions, and Coulombic efficiency on electrolyte imbalance. The model also demonstrates excellent computational efficiency for simulating 41 cycles (around 37,000 points) within 3~4 s. The demonstrated efficiency and accuracy for predicting cross-over and its impacts on voltage, capacity&energy decay, membrane concentrations, and electrolyte imbalance makes it a reliable tool for optimizing RFBs’ long-term performance.

25 ENERGY STORAGE↗
DOE OSTI2026

Unveiling in-gap states and Majorana zero modes in superconductor-topological insulator bilayer model

Interfaces between topological insulators and superconductors (SCs) are promising platforms for realizing Majorana zero modes (MZMs) via the superconducting proximity effect. We introduce a bilayer model consisting of the surface states of a three-dimensional topological insulator (3DTI) coupled to an 𝑠-wave superconductor and systematically study the role of interlayer tunneling strength (𝑡⊥) motivated by the recent growth of the Fe⁡(Te,Se)/Bi 2 ⁢Te 3 heterostructure. We find that increasing 𝑡⊥ shifts the proximity-induced (PrI) gap minimum away from the Γ point, giving rise to momentum-selective interference patterns that manifest as spatial oscillations in the in-gap states. We introduce an antidot with a magnetic vortex in the model and investigate the resulting in-gap states, including MZMs and Caroli–de Gennes–Matricon (CdGM) modes. With increasing hybridization strength, the energy separation between MZMs and CdGM states increases, enhancing the isolation of MZMs. Importantly, in the strong hybridization limit, the leading CdGM separation remains large in spite of the decrease in the PrI gap. Spin- and spatially resolved wavefunction analysis reveals angular momentum asymmetries absent in conventional 𝑠-wave systems. A direct comparison with a standalone 𝑠-wave superconductor confirms the emergence of distinct 𝑝-wave-like features in the bilayer geometry. Our results provide experimentally relevant predictions for tuning the stability of MZMs and their differentiation from the CdGM modes in SC-3DTI heterostructures and offer a theoretical framework for probing unconventional superconductivity in engineered topological systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Polyakov model in ’t Hooft flux background: a quantum mechanical reduction with memory

We construct a compactification of Polyakov model on T 2 × $\mathbb{R}$ down to quantum mechanics which remembers non-perturbative aspects of field theory even at an arbitrarily small area. Standard compactification on small T 2 × $\mathbb{R}$ possesses a unique perturbative vacuum (zero magnetic flux state), separated parametrically from higher flux states, and the instanton effects do not survive in the Born-Oppenheimer approximation. By turning on a background magnetic GNO flux in co-weight lattice corresponding to a non-zero ’t Hooft flux, we show that N-degenerate vacua appear at small torus, and there are N - 1 types of flux changing instantons between them. We construct QM instantons starting with QFT instantons using the method of replicas. For example, SU(2) gauge theory with flux reduces to the double-well potential where each well is a fractional flux state. Despite the absence of a mixed anomaly, the vacuum structure of QFT and the one of QM are continuously connected. We also compare the quantum mechanical reduction of the Polyakov model with the deformed Yang-Mills, by coupling both theories to TQFTs. In particular, we compare the mass spectrum for dual photons and energy spectrum in the QM limit. We give a detailed description of critical points at infinity in the semi-classical expansion, and their role in resurgence structure.

nonperturbative effects↗
DOE OSTI2023

Long-pulse high-performance H-mode plasmas achieved on EAST

A record duration of a 310 s H-mode plasma (H 98y2 ~ 1.3, n e /n GW ~ 0.7, f BS > 50%) has been recently achieved on experimental advanced superconducting tokamak (EAST) with metal walls, exploiting the device's improved long-pulse capabilities. The experiment demonstrates good control of tungsten concentration, core/edge MHD stability, and particle and heat exhaust with an ITER-like tungsten divertor and zero injected torque, establishing a milestone on the path to steady-state long-pulse high-performance scenarios in support of ITER and CFETR. Important synergistic effects are leveraged toward this result, which relies purely on radio frequency (RF) powers for heating and current drive (H&CD). On-axis electron cyclotron heating enhances the H&CD efficiency from lower hybrid wave injection, increasing confinement quality and enabling fully non-inductive operation at high density (n e /n GW ~ 70%) and high poloidal beta (β P ~ 2.5). A small-amplitude grassy edge localized mode regime facilitates the RF power coupling to the H-mode edge and reduces divertor sputtering/erosion. The high energy confinement quality (H 98y2 ~ 1.3) is achieved with the experimental and simulated results pointing to the strong effect of Shafranov shift on turbulence. Transport analysis suggests that trapped electron modes dominate in the core region during the record discharge. The detailed physics processes (RF synergy, core-edge integration, confinement properties, etc.) of the steady-state operation will be illustrated in the content. In the future, EAST will aim at accessing more relevant dimensionless parameters to develop long-pulse high-performance plasma toward ITER and CFETR steady-state advanced operation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Parity-odd four-point correlation function from the DESI data release 1 luminous red galaxy sample

The parity-odd four-point function provides a unique probe of fundamental symmetries and potential new physics in the large-scale structure of the Universe. We present measurements of the parity-odd four-point function using the Dark Energy Spectroscopic Instrument (DESI) DR1 luminous red galaxy (LRG) sample and assess its detection significance. Our analysis considers both auto- and cross-correlations, using two complementary approaches to the covariance: (i) the full analytic covariance matrix applied to the uncompressed data vector, and (ii) a compressed data vector combined with a hybrid covariance matrix constructed from simulations and analytic estimates. When using the full analytic covariance matrix without corrections, we observe apparent auto-correlation signals with significance up to 4⁢𝜎. However, this excess is also consistent with a mismatch between the statistical fluctuations estimated from the simulations and those present in the real data. Our findings therefore suggest that the parity-odd signal in the current DESI DR1 LRG sample is consistent with zero. We note, however, that the low completeness of this sample may have a non-negligible impact on the detection sensitivity. Future data releases with improved completeness will be crucial for further investigation.

Hou, Jiamin [Ludwig-Maximilians-Universität; Unive↗
DOE OSTI2020

Ab Initio Phase Diagram of Tungsten

The phase diagram of tungsten (W) to a pressure (P) of 2500 GPa is investigated using a comprehensive ab initio approach that includes (i) the calculation of the zero temperature (T) free energies (enthalpies) of different solid structures, (ii) the quantum molecular dynamics simulation of the melting curves of different solid structures, (iii) the derivation of the analytic form for the solid-solid phase transition boundary, and (iv) the simulations of the solidification of liquid W into the final solid states on both sides of the solid-solid phase transition boundary, in order to confirm the corresponding analytic form. There are two solid structures confirmed to be present on the phase diagram of W, the ambient body-centered cubic (bcc) and the high-pressure double hexagonal close-packed (dhcp). At T = 0, the bcc-dhcp transition occurs at 1060 GPa, and the transition boundary has a positive slope dT/dP : the bcc-dhcp-liquid triple point is at (P, T) = (1675 GPa, 23680 K).

36 MATERIALS SCIENCE↗
DOE OSTI2020

Iron-based binary ferromagnets for transverse thermoelectric conversion

Thermoelectric generation using the anomalous Nernst effect (ANE) has great potential for application in energy harvesting technology because the transverse geometry of the Nernst effect should enable efficient, large-area and flexible coverage of a heat source. For such applications to be viable, substantial improvements will be necessary not only for their performance but also for the associated material costs, safety and stability. In terms of the electronic structure, the anomalous Nernst effect (ANE) originates from the Berry curvature of the conduction electrons near the Fermi energy. To design a large Berry curvature, several approaches have been considered using nodal points and lines in momentum space. Here we perform a high-throughput computational search and find that 25 percent doping of aluminium and gallium in alpha iron, a naturally abundant and low-cost element, dramatically enhances the ANE by a factor of more than ten, reaching about 4 and 6 microvolts per kelvin at room temperature, respectively, close to the highest value reported so far. The comparison between experiment and theory indicates that the Fermi energy tuning to the nodal web—a flat band structure made of interconnected nodal lines—is the key for the strong enhancement in the transverse thermoelectric coefficient, reaching a value of about 5 amperes per kelvin per metre with a logarithmic temperature dependence. We have also succeeded in fabricating thin films that exhibit a large ANE at zero field, which could be suitable for designing low-cost, flexible microelectronic thermoelectric generators.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Computational discovery of two-dimensional rare-earth iodides: promising ferrovalley materials for valleytronics

Two-dimensional Ferrovalley materials with intrinsic valley polarization are rare but highly promising for valley-based nonvolatile random access memory and valley filter devices. These ferromagnetic materials exhibit valleys at or near the Fermi level with intrinsic magnetism. The strong coupling between magnetism and spin–orbit coupling induces intrinsic valley polarization. Using Kinetically Limited Minimization, an unconstrained crystal structure prediction algorithm, and prototype sampling based on first-principles calculations, we have discovered new Ferrovalley materials, rare-earth iodides RI 2 , where R is a rare-earth element belonging to Sc, Y, or La-Lu, and I is Iodine. The rare-earth iodides are layered and demonstrate either 2H, 1T, or 1T d phase as the ground state in bulk, analogous to transition metal dichalcogenides (TMDCs). The calculated exfoliation energy of monolayers (MLs) is comparable to that of graphene and TMDCs, suggesting possible experimental synthesis. The MLs in the 2H phase exhibit ferromagnetism due to unpaired electrons in d and f orbitals. Throughout the rare-earth series, d bands have valley polarization at K and $\overline{K}$ points in the Brillouin zone in the vicinity of the Fermi level. Large intrinsic valley polarization in the range of 15–143 meV without external stimuli is observed in these Ferrovalley materials, which can be enhanced further by applying an in-plane bi-axial strain. These valleys can selectively be probed and manipulated for information storage and processing, potentially offering superior performance beyond conventional electronics and spintronics. Here we further show that the 2H ferromagnetic phase of RI 2 MLs possesses non-zero Berry curvature and exhibits anomalous valley Hall effect with considerable anomalous Hall conductivity. Our work will incite exploratory synthesis of the predicted Ferrovalley materials and their application in valleytronics and beyond.

2D materials↗
DOE OSTI2025

DAmodel: hierarchical Bayesian modelling of DA white dwarfs for spectrophotometric calibration

We use hierarchical Bayesian modelling to calibrate a network of 32 all-sky faint DA white dwarf (DA WD) spectrophotometric standards (⁠16.5 < V , 19.5⁠) alongside three CALSPEC standards, from 912 Å to 32 μm. The framework is the first of its kind to jointly infer photometric zero points and WD parameters (surface gravity log g⁠, effective temperature T eff ⁠, extinction A V ⁠, dust relation parameter R V ) by simultaneously modelling both photometric and spectroscopic data. We model panchromatic Hubble Space Telescope Wide Field Camera 3 (HST/WFC3) UVIS and IR photometry, HST/STIS UV spectroscopy, and ground-based optical spectroscopy to sub-per cent precision. Photometric residuals for the sample are the lowest yet yielding < 0.004 mag RMS on average from the UV to the NIR, achieved by jointly inferring time-dependent changes in system sensitivity and WFC3/IR count-rate nonlinearity. Our GPU-accelerated implementation enables efficient sampling via Hamiltonian Monte Carlo, critical for exploring the high-dimensional posterior space. The hierarchical nature of the model enables population analysis of intrinsic WD and dust parameters. Inferred spectral energy distributions from this model will be essential for calibrating the James Webb Space Telescope as well as next-generation surveys, including Vera Rubin Observatory’s Legacy Survey of Space and Time and the Nancy Grace Roman Space Telescope.

methods: statistical↗
DOE OSTI2023

On falling film evaporator – A review of mechanisms and critical assessment of correlation on a horizontal tube bundle with updated development

Optimization of the energy efficiency and reduction of the greenhouse gas emission for the heat pump system is imperative to meet the net-zero goals at 2050. Chillers with falling film evaporator design not only possess better system performance but also contain less refrigerant inventory. Hence, accurate prediction of the evaporator performance is pivotal especially when charged with low-GWP refrigerant. Here, the study reviews the correlations for falling film evaporator with a horizontal tube bundle configuration. The major efforts of this study include four tasks: (a) literature review of the experimental studies and available empirical correlations; (b) comprehensive discussion of the falling film evaporation heat transfer mechanism; (c) development of a new rationally based correlation based on available literature; and (d) comparison of different correlations based on the existing data. The collected data includes 4114 data points from 8 sources, 6 refrigerants (R-600a, R-290, R-245fa, R-134a, R-1234ze(E), R-123), 5 types of the tubes (Plain, Turbo-GII-HP, GEWA-B, Low-fin, High-Flux), liquid Weber number from 2.2 × 10 –6 to 0.7, imposed vapor Weber number from 0 to 37.2, heat flux from 2.5 to 151.5 kW/m 2 , and film Reynolds number from 1 to 3159.8. The new correlation gives a MAD of 28.3%, and an R 2 of 0.86. Yet, the developed correlation considers various heat transfer mechanisms, including the transition point from falling film evaporation to the nucleate boiling, local evaporation, dry-out, mist flow, imposed flow, and enhanced tube effects.

42 ENGINEERING↗
DOE OSTI2025

Measurement of the parity-violating asymmetry in the 𝑁 → Δ transition at low 𝑄 2

Here, we report the measurement of the parity-violating asymmetry in the 𝑁 → Δ transition via the 𝑒 − + 𝑝 → 𝑒 − + Δ + reaction at two different kinematic points with low four-momentum transfer 𝑄 2 . Measurements were made with incident electron beam energies of 0.877 and 1.16 GeV, corresponding to 𝑄 2 values of 0.0111 and 0.0208 (GeV/c) 2 , respectively. These measurements put constraints on a low-energy constant in the weak Lagrangian, 𝑑 Δ , corresponding to a parity-violating electric-dipole transition matrix element. This matrix element has been shown to be large in the strangeness-changing channel, via weak hyperon decays such as Σ + → 𝑝⁢𝛾. The measurements reported here constrain 𝑑 Δ in the strangeness-conserving channel. The final asymmetries were −0.65±1.00⁢(stat.)±1.02⁢(syst.) ppm (parts per million) for 0.877 GeV and −3.59±0.82⁢(stat.)±1.33⁢(syst.) ppm for 1.16 GeV. With these results we deduce a small value for 𝑑 Δ , consistent with zero, in the strangeness-conserving channel, in contrast to the large value for 𝑑 Δ previously reported in the strangeness-changing channel.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Majorana Edge Modes in Isolated Wires

Topological superconductors are believed to host exotic quasiparticle excitations known as Majorana zero modes (MZMs), with much of the evidence based on BCS mean-field theory. The direct application of mean-field arguments is tenuous in finite, isolated systems relevant in some experiments. Here, we develop a new correlation-based method for identifying MZMs in interacting, number-conserving systems. Using the density matrix renormalization group, we study fermion number-conserving models with long-range interactions, which under periodic boundary conditions exhibit robust topological and nontopological superconductivity, tuned by the strength of interaction [Ortiz et al., Phys. Rev. Lett. 113, 267002 (2014)]. We find evidence that, on the topological side, Majorana edge modes appear in open chains, manifesting as the vanishing of the energy splitting between odd- and even-parity ground states with increasing system size. Additionally, off-diagonal two-point correlation functions show nonlocal, parity-dependent edge effects. These correlations reveal the spatial structure of Majorana modes in the many-body wave function. We show that the correlation diagnostic applies broadly, including to short-range interacting models, where topological superconductivity is more fragile due to the absence of a bulk excitation gap.

Thomas-Markarian, Jaden↗