🔬 materials science

Wavefunction-Free Approach for Predicting Nonlinear Responses in Weyl Semimetals

This paper introduces a wavefunction-free approach that achieves a 1000-fold computational speedup in predicting nonlinear responses like the circular photogalvanic effect in Weyl semimetals, enabling the efficient discovery of materials such as Ta $_3$ S $_2$ with significantly enhanced photocurrents.

Original authors: Mohammad Yahyavi, Ilya Belopolski, Yuanjun Jin, Yilin Zhao, Jinyang Ni, Naizhou Wang, Yi-Chun Hung, Zi-Jia Cheng, Tyler A. Cochran, Tay-Rong Chang, Wei-bo Gao, Su-Yang Xu, Jia-Xin Yin, Qiong Ma, Md Sh

Published 2026-02-02

📖 4 min read☕ Coffee break read

CC BY 4.0

Original authors: Mohammad Yahyavi, Ilya Belopolski, Yuanjun Jin, Yilin Zhao, Jinyang Ni, Naizhou Wang, Yi-Chun Hung, Zi-Jia Cheng, Tyler A. Cochran, Tay-Rong Chang, Wei-bo Gao, Su-Yang Xu, Jia-Xin Yin, Qiong Ma, Md Shafayat Hossain, Arun Bansil, Naoto Nagaosa, Guoqing Chang

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Imagine you are trying to predict how a new type of super-material will react when hit by a laser. In the world of quantum physics, these materials are like complex mazes made of invisible energy. Traditionally, to figure out how they react, scientists had to map out every single path a tiny particle (an electron) could take through this maze. This process is like trying to count every single grain of sand on a beach to predict how the tide will move. It is incredibly accurate, but it takes so much time and computer power that it's almost impossible to test many different materials quickly.

This paper introduces a "shortcut" that changes the game. Instead of counting every grain of sand, the researchers found a way to predict the outcome just by looking at the shape of the beach itself.

Here is the breakdown of their discovery:

1. The Problem: The "Wavefunction" Bottleneck

To understand how these materials work, scientists usually have to calculate something called a "wavefunction." Think of the wavefunction as a detailed, 3D blueprint of every possible move an electron can make. Calculating this for a whole material is like trying to solve a massive, multi-layered jigsaw puzzle where the pieces keep changing shape. It takes supercomputers days or weeks to solve just one puzzle. This makes it very hard to screen thousands of materials to find the best ones for future technology.

2. The Solution: The "Recipe" Shortcut

The researchers realized that for a specific type of particle called a "Weyl fermion" (which acts like a massless, super-fast messenger in these materials), you don't need the full blueprint. You only need a few key numbers from the material's "recipe" (called Hamiltonian parameters).

They developed a new mathematical formula that skips the complex blueprint entirely. Instead of asking, "What is the electron doing right now?" they ask, "What are the basic rules of the game?"

The Analogy: Imagine you want to know how fast a car will go. The old way was to build a full simulation of the engine, the tires, the road friction, and the wind for every single second of the drive. The new way is just to look at the car's horsepower and the slope of the hill. If you know those two numbers, you can predict the speed instantly without simulating the whole drive.

3. The Result: A 1,000-Fold Speed Boost

By using this shortcut, the team made calculations 1,000 times faster.

Real-world example: In their study, calculating the response for a material called TaAs took their old method 106 seconds (and huge amounts of memory). Their new method did the same job in a fraction of that time.
This speed allows scientists to test hundreds of materials in the time it used to take to test just one.

4. The Discovery: Finding the "Super-Material"

Because they could now test materials so quickly, they screened a long list of known "Weyl semimetals" (materials with these special particles).

They found that a material called Ta3S2 is a hidden champion. It produces a "photocurrent" (electricity generated by light) that is 10 times stronger than the previous record-holder, TaAs.
Even better, they discovered that if you squeeze (strain) this material in a specific direction, it could become 10 times stronger again. This is like finding a car that not only goes fast but can be tuned to go even faster just by adjusting the suspension.

5. Why It Matters (According to the Paper)

The paper claims this approach does two main things:

Efficiency: It allows scientists to rapidly screen and optimize materials for "nonlinear" effects (where light and electricity interact in complex ways). This is crucial for developing next-generation devices like ultra-fast wireless communication and advanced sensors.
Understanding: It provides a clearer, simpler way to understand why these materials behave the way they do. Instead of getting lost in complex math, scientists can now see that these effects are caused by simple "tilts" and "warping" in the material's structure, breaking a fundamental symmetry in physics.

Summary

The authors have built a "fast-forward" button for quantum material discovery. By replacing a slow, complex calculation with a simple formula based on basic material properties, they found a new material (Ta3S2) that is vastly superior for generating electricity from light. They also showed that this method works for other types of physics problems, not just light, making it a powerful new tool for designing the quantum technologies of the future.

Technical Summary: Wavefunction-Free Approach for Predicting Nonlinear Responses in Weyl Semimetals

Problem Statement
While Density Functional Theory (DFT) has revolutionized the prediction of material ground states by bypassing intractable many-body wavefunction calculations, predicting nonlinear electromagnetic responses remains a significant computational bottleneck. Current methods for evaluating nonlinear responses—critical for next-generation quantum devices—rely heavily on explicit integration of intricate Bloch wavefunctions across the Brillouin zone (BZ). This "brute-force" approach is computationally expensive and resource-intensive, hindering the high-throughput screening of materials with tailored nonlinear properties. Specifically, calculating the Circular Photogalvanic Effect (CPGE) and Berry-curvature dipole (BCD) in topological materials like Weyl semimetals requires evaluating complex sums over bands and momenta involving wavefunction matrix elements, often requiring terabytes of memory and millions of seconds of CPU time.

Methodology
The authors propose a "wavefunction-free" approach that circumvents explicit wavefunction integrations by leveraging the one-to-one correspondence between the free parameters of Weyl fermion Hamiltonians and their associated nonlinear responses. Instead of calculating the full wavefunction geometry, the method derives closed-form empirical formulas directly from Hamiltonian parameters (velocities, tilting terms, warping terms, and chemical potential).

The methodology proceeds as follows:

Hamiltonian Analysis: The authors analyze a high-symmetry Weyl fermion Hamiltonian ( $H_0$ ) and introduce perturbations, specifically a tilting term ( $H_{p1}$ ) and a warping term ( $H_{p2}$ ), to emulate the low-symmetry conditions found in real materials (e.g., TaAs) that break momentum-space balance and yield net nonlinear currents.
Derivation of Empirical Formulas: By performing systematic numerical simulations on model Hamiltonians, the authors identify the functional dependencies of the nonlinear response tensors on the Hamiltonian parameters. They derive analytical expressions for the transverse CPGE ( $\beta_{ab}$ ) and the BCD ( $D_{ab}$ ) that depend solely on parameters such as anisotropic velocities ( $\nu_i$ ), tilting velocities ( $\nu^t_i$ ), the warping coefficient ( $\delta_i$ ), and the energy window defined by the chemical potential ( $\mu$ ).
Validation: The derived formulas are validated against conventional wavefunction-based calculations for both Type-I ( $W_T < 1$ ) and Type-II ( $W_T > 1$ ) Weyl semimetals. The approach is then applied to real materials by extracting Hamiltonian parameters from DFT band structures and fitting them to the empirical formulas, bypassing the need for Wannier function generation and subsequent wavefunction integration.

Key Contributions and Results

Computational Efficiency: The primary contribution is a reduction in computational cost by approximately three orders of magnitude (1000-fold). For the material TaAs, the conventional method required 106 seconds and over 3 TB of CPU memory for the final calculation step, whereas the wavefunction-free approach reduced the process to two steps (DFT and parameter fitting), with the most time-consuming step being the initial DFT calculation (103 seconds).
Discovery of High-Performance Materials: Applying this efficient screening method to a library of known Weyl semimetals (including TaAs, TaP, NbP, Ta3S2, etc.), the authors identified that Ta3S2 exhibits transverse CPGE photocurrents an order of magnitude greater than those in TaAs.
Strain Engineering: The study demonstrates that the nonlinear response can be further enhanced by manipulating Hamiltonian parameters via external strain. For Ta3S2, applying a -5% strain along the [110] plane was predicted to increase the CPGE amplitude by a factor of 10. The formulas also reveal that tuning the velocity ratio ( $\nu_x/\nu_y$ ) allows for the anisotropic tuning of CPGE.
Generalization to Transport: The approach was extended beyond optical effects to the nonlinear Hall effect induced by the Berry-curvature dipole (BCD). A wavefunction-free formula for BCD was derived, showing that contributions from Weyl fermions dominate over trivial pockets, even when pocket sizes are comparable.
Physical Insight: The work provides a transparent physical interpretation of nonlinear responses, linking them directly to the breaking of "generalized" Lorentz symmetry in Weyl fermions via tilting and warping terms, rather than the complex geometry of Bloch wavefunctions.

Significance
The paper claims that this framework paves the way for the substantially more efficient screening and optimization of nonlinear electromagnetic properties in topological quantum materials. By eliminating the dependence on explicit wavefunctions, the method enables the rapid identification of promising material candidates for next-generation quantum devices operating in the terahertz regime. Furthermore, the work establishes a general framework that bridges high-energy physics concepts (Weyl fermions) with condensed matter applications, offering clear physical insights into the origins of nonlinear responses. The authors suggest that with advances in synthesis, these strategies can guide the rational design of novel quantum materials capable of achieving significantly larger photocurrents.