Search

We present an efficient linear-scaling algorithm for evaluating the analytical force and stress contributions derived from the exact-exchange energy, a key component in hybrid functional calculations. The algorithm, working equally well for molecular and periodic systems, is formulated within the framework of numerical atomic orbital (NAO) basis sets and takes advantage of the localized resolution-of-identity (LRI) technique for treating the two-electron Coulomb repulsion integrals. The linear-scaling behavior is realized by fully exploiting the sparsity of the expansion coefficients resulting from the strict locality of the NAOs and the LRI ansatz. Our implementation is massively parallel, and enables efficient structural relaxation based on hybrid density functionals for bulk materials containing thousands of atoms. In this work, we will present a detailed description of our algorithm and benchmark the performance of our implementation using illustrating examples. By optimizing the structures of the pristine and doped halide perovskite material CsSnI$_3$ with different functionals, we find that in the presence of lattice strain, hybrid functionals provide a more accurate description of the stereochemical expression of the lone pair.

With the fast developments of high-performance computing, first-principles methods based on quantum mechanics play a significant role in materials research, serving as fundamental tools for predicting and analyzing various properties of materials. However, the inherent complexity and substantial computational demands of first-principles algorithms, such as density functional theory, limit their use in larger systems. The rapid development of heterogeneous computing, particularly General-Purpose Graphics Processing Units (GPGPUs), has heralded new prospects for enhancing the performance and cost-effectiveness of first-principles algorithms. We utilize GPGPUs to accelerate the electronic structure algorithms in Atomic-orbital Based Ab-initio Computation at USTC (ABACUS), a first-principles computational package based on the linear combination of atomic orbitals (LCAO) basis set. We design algorithms on GPGPU to efficiently construct and diagonalize the Hamiltonian of a given system, including the related force and stress calculations. The effectiveness of this computational acceleration has been demonstrated through calculations on twisted bilayer graphene with the system size up to 10,444 atoms.

Simulating strongly-correlated quantum many-body systems at finite temperatures is a significant challenge in computational physics. In this work, we present a scalable finite-temperature tensor network algorithm for two-dimensional quantum many-body systems. We employ the (fermionic) projected entangled pair state (PEPS) to represent the vectorization of the quantum thermal state and utilize a stochastic reconfiguration method to cool down the quantum states from infinite temperature. We validate our method by benchmarking it against the 2D antiferromagnetic Heisenberg model, the $J_1$-$J_2$ model, and the Fermi-Hubbard model, comparing physical properties such as internal energy, specific heat, and magnetic susceptibility with results obtained from stochastic series expansion (SSE), exact diagonalization, and determinant quantum Monte Carlo (DQMC).

LibRPA is a software package designed for efficient calculations of random phase approximation (RPA) electron correlation energies from first principles using numerical atomic orbital (NAOs). Leveraging a localized resolution of identity (LRI) technique, LibRPA achieves $O(N^2)$ or better scaling behavior, making it suitable for large-scale calculation of periodic systems. Implemented in C++ and Python with MPI/OpenMP parallelism, LibRPA integrates seamlessly with NAO-based density functional theory (DFT) packages through flexible file-based and API-based interfaces. In this work, we present the theoretical framework, algorithm, software architecture, and installation and usage guide of LibRPA. Performance benchmarks, including the parallel efficiency with respect to the computational resources and the adsorption energy calculations for H$_2$O molecules on graphene, demonstrate its nearly ideal scalability and numerical reliability. LibRPA offers a useful tool for RPA-based calculations for large-scale extended systems., Comment: 29 pages, 4 figures

The recent discovery of high-temperature superconductivity in high-pressurized La$_3$Ni$_2$O$_{7-\delta}$ has garnered significant attention. Using density functional theory, we investigate the magnetic properties of La$_3$Ni$_2$O$_{7-\delta}$ at ambient pressure. Our calculations suggest that with $\delta=0$, the double spin stripe phase is favored as the magnetic ground state. Oxygen vacancies may effectively turn nearest Ni spins into \textit{charge} sites. Consequently, with moderate $\delta$ values, our theoretical magnetic ground state exhibits characteristics of both double spin stripe and spin-charge stripe configurations, providing a natural explanation to reconcile the seemingly contradictory experimental findings that suggest both the configurations as candidates for the spin-density-wave phase. With higher $\delta$ values, we anticipate the ground state to become a spin-glass-like noncollinear magnetic phase with only short-range order. The oxygen vacancies are expected to significantly impact the magnetic excitations and the transition temperatures $T_{SDW}$. Notably, the magnetic ordering also induces concomitant charge ordering and orbital ordering, driven by spin-lattice coupling under the low symmetry magnetic order. We further offer a plausible explanation for the experimental observations that the measured $T_{SDW}$ appears insensitive to the variation of samples and the lack of direct evidence for long-range magnetic ordering.

Projected Entangled Pair States (PEPS) are recognized as a potent tool for exploring two-dimensional quantum many-body systems. However, a significant challenge emerges when applying conventional PEPS methodologies to systems with periodic boundary conditions (PBC), attributed to the prohibitive computational scaling with the bond dimension. This has notably restricted the study of systems with complex boundary conditions. To address this challenge, we have developed a strategy that involves the superposition of PEPS with open boundary conditions (OBC) to treat systems with PBC. This approach significantly reduces the computational complexity of such systems while maintaining their translational invariance and the PBC. We benchmark this method against the Heisenberg model and the $J_1$-$J_2$ model, demonstrating its capability to yield highly accurate results at low computational costs, even for large system sizes. The techniques are adaptable to other boundary conditions, including cylindrical and twisted boundary conditions, and therefore significantly expands the application scope of the PEPS approach, shining new light on numerous applications.

We propose both a theoretical and a methodological framework to address a critical challenge in applying deep learning to physical systems: the reconciliation of non-linear expressiveness with SO(3)-equivariance in predictions of SO(3)-equivariant quantities. Inspired by covariant theory in physics, we present a solution by exploring the mathematical relationships between SO(3)-invariant and SO(3)-equivariant quantities and their representations. We first construct theoretical SO(3)-invariant quantities derived from the SO(3)-equivariant regression targets, and use these invariant quantities as supervisory labels to guide the learning of high-quality SO(3)-invariant features. Given that SO(3)-invariance is preserved under non-linear operations, the encoding process for invariant features can extensively utilize non-linear mappings, thereby fully capturing the non-linear patterns inherent in physical systems. Building on this, we propose a gradient-based mechanism to induce SO(3)-equivariant encodings of various degrees from the learned SO(3)-invariant features. This mechanism can incorporate non-linear expressive capabilities into SO(3)-equivariant representations, while theoretically preserving their equivariant properties as we prove, establishing a strong foundation for regressing complex SO(3)-equivariant targets. We apply our theory and method to the electronic-structure Hamiltonian prediction tasks, experimental results on eight benchmark databases covering multiple types of systems and challenging scenarios show substantial improvements on the state-of-the-art prediction accuracy of deep learning paradigm. Our method boosts Hamiltonian prediction accuracy by up to 40% and enhances downstream physical quantities, such as occupied orbital energy, by a maximum of 76%.

Deep learning for predicting the electronic-structure Hamiltonian of quantum systems necessitates satisfying the covariance laws, among which achieving SO(3)-equivariance without sacrificing the non-linear expressive capability of networks remains unsolved. To navigate the harmonization between equivariance and expressiveness, we propose a deep learning method synergizing two distinct categories of neural mechanisms as a two-stage encoding and regression framework. The first stage corresponds to group theory-based neural mechanisms with inherent SO(3)-equivariant properties prior to the parameter learning process, while the second stage is characterized by a non-linear 3D graph Transformer network we propose, featuring high capability on non-linear expressiveness. The novel combination lies in the point that, the first stage predicts baseline Hamiltonians with abundant SO(3)-equivariant features extracted, assisting the second stage in empirical learning of equivariance; and in turn, the second stage refines the first stage's output as a fine-grained prediction of Hamiltonians using powerful non-linear neural mappings, compensating for the intrinsic weakness on non-linear expressiveness capability of mechanisms in the first stage. Our method enables precise, generalizable predictions while capturing SO(3)-equivariance under rotational transformations, and achieves state-of-the-art performance in Hamiltonian prediction on six benchmark databases.

In materials without inversion symmetry, Berry curvature dipole (BCD) arises from the uneven distribution of Berry curvature in momentum space. This leads to nonlinear anomalous Hall effects even in systems with preserved time-reversal symmetry. A key goal is to engineer systems with prominent BCD near the Fermi level. Notably, TaAs, a type-I Weyl semimetal, exhibits substantial Berry curvature but a small BCD around the Fermi level. In this study, we employed first-principles methods to comprehensively investigate the BCD in TaAs. Our findings reveal significant cancellation effects not only within individual Weyl points but crucially, among distinct Weyl point pairs in bulk TaAs. We propose a strategic approach to enhance the BCD in TaAs by employing a layer-stacking technique. This greatly amplifies the BCD compared to the bulk material. By tuning the number of slab layers, we can selectively target specific Weyl point pairs near the Fermi level, while quantum confinement effects suppress contributions from other pairs, mitigating cancellation effects. Especially, the BCD of an 8-layer TaAs slab surpasses the bulk value near the Fermi level by orders of magnitude.

High-energy, few-cycle laser pulses are essential for numerous applications in the fields of ultrafast optics and strong-field physics, due to their ultrafast temporal resolution and high peak intensity. In this work, different from the traditional hollow-core fibers and multiple thin solid plates, we represent the first demonstration of the octave-spanning supercontinuum broadening by utilizing multiple ultrathin liquid films (MTLFs) as the nonlinear media. The continuum covers a range from 380 to 1050 nm, corresponding to a Fourier transform limit pulse width of 2.5 fs, when 35 fs Ti:sapphire laser pulse is applied on the MTLFs. The output pulses are compressed to 3.9 fs by employing chirped mirrors. Furthermore, a continuous high-order harmonic spectrum up to the 33rd order is realized by subjecting the compressed laser pulses to interact with Kr gas. The utilization of flowing water films eliminates permanent optical damage and enables wider and stronger spectrum broadening. Therefore, this MTLFs scheme provides new solutions for the generation of highly efficient femtosecond supercontinuum and nonlinear pulse compression, with potential applications in the fields of strong-field physics and attosecond science., Comment: 9 pages, 5 figures

The random phase approximation (RPA) as formulated as an orbital-dependent, fifth-rung functional within the density functional theory (DFT) framework offers a promising approach for calculating the ground-state energies and the derived properties of real materials. Its widespread use to large-size, complex materials is however impeded by the significantly increased computational cost, compared to lower-rung functionals. The standard implementation exhibits an $\mathcal{O}(N^4)$-scaling behavior with respect to system size $N$. In this work, we develop a low-scaling RPA algorithm for periodic systems, based on the numerical atomic orbital (NAO) basis-set framework and a localized variant of the resolution of identity (RI) approximation. The rate-determining step for RPA calculations -- the evaluation of non-interacting response function matrix, is reduced from $\mathcal{O}(N^4)$ to $\mathcal{O}(N^2)$ by just exploiting the sparsity of the RI expansion coefficients, resultant from localized RI (LRI) scheme and the strict locality of NAOs. The computational cost of this step can be further reduced to linear scaling if the decay behavior of the Green's function in real space can be further taken into account. Benchmark calculations against existing $\textbf k$-space based implementation confirms the validity and high numerical precision of the present algorithm and implementation. The new RPA algorithm allows us to readily handle three-dimensional, closely-packed solid state materials with over 1000 atoms. The algorithm and numerical techniques developed in this work also have implications for developing low-scaling algorithms for other correlated methods to be applicable to large-scale extended materials.

We present PYATB, a Python package designed for computing band structures and related properties of materials using the ab initio tight-binding Hamiltonian. The Hamiltonian is directly obtained after conducting self-consistent calculations with first-principles packages using numerical atomic orbital (NAO) bases, such as ABACUS. The package comprises three modules: Bands, Geometric, and Optical. In the Bands module, one can calculate essential properties of band structures, including the partial density of states (PDOS), fat bands, Fermi surfaces, and Weyl/Dirac points. The band unfolding method is utilized to obtain the energy band spectra of a supercell by projecting the electronic structure of the supercell onto the Brillouin zone of the primitive cell. With the Geometric module, one can compute the Berry phase and Berry curvature-related quantities, such as electric polarization, Wilson loops, Chern numbers, and anomalous Hall conductivities. The Optical module offers a range of optical property calculations, including optical conductivity and nonlinear optical responses, such as shift current and Berry curvature dipole.

Hybrid density functional calculation is indispensable to accurate description of electronic structure, whereas the formidable computational cost restricts its broad application. Here we develop a deep equivariant neural network method (named DeepH-hybrid) to learn the hybrid-functional Hamiltonian from self-consistent field calculations of small structures, and apply the trained neural networks for efficient electronic-structure calculation by passing the self-consistent iterations. The method is systematically checked to show high efficiency and accuracy, making the study of large-scale materials with hybrid-functional accuracy feasible. As an important application, the DeepH-hybrid method is applied to study large-supercell Moir\'{e} twisted materials, offering the first case study on how the inclusion of exact exchange affects flat bands in the magic-angle twisted bilayer graphene.

The bulk photovoltaic effect (BPVE) refers to the phenomenon of generating photocurrent or photovoltage in homogeneous noncentrosymmetric materials under illumination, and the intrinsic contribution to the BPVE is known as the shift current effect. We calculate the shift current conductivities of the ferroelectric SnTe monolayer using first-principles methods. We find that the monolayer SnTe has giant shift-current conductivity near the valley points. More remarkably, the linear optical absorption coefficient at this energy is very small, and therefore leads to an enormous Glass coefficient that is four orders of magnitude larger than that of BaTiO$_3$. The unusual shift-current effects are further investigated using a three-band model. We find that the giant shift current conductivities and Glass coefficient are induced by the nontrivial energy band geometries near the valley points, where the shift-vector diverges. This is a prominent example that the band geometry can play essential roles in the fundamental properties of solids.

The antiferromagnetic (AFM) topological insulator MnBi$_{2}$Te$_{4}$ was believed to have a topological surface state (TSS) with large band gap due to the ferromagnetic (FM) order on the surface, and be able to host the long-sought axion states. However, recent angle-resolved photoemission spectroscopy (ARPES) experiments indicate that the TSS is gapless, contradicting the theoretical predictions. Meanwhile, several experiments have suggested that there is robust out-of-plane FM order on the surface of MnBi$_{2}$Te$_{4}$. To understand these seemingly contradictory results, we carry out comprehensive first-principles calculations to investigate the interplay between the surface magnetism and the TSS. Our calculations provide direct evidence that in a wide range of parameters, the (nearly) gapless TSS can coexist with the surface FM order, therefore solving the paradox of the surface magnetism and the gapless TSS. We further show that proximity effects can be a promising route to open the gap in the TSS of MnBi$_{2}$Te$_{4}$. Our research deepens the understanding of the relationship between surface magnetism and TSS.

The quantum spin liquid (QSL) state has been searched intensively in Kitaev-like materials, such as the Iridium oxides $A_2$IrO$_3$ and $\alpha$-RuCl$_3$. The half-filled Kitaev-Hubbard model with bond dependent hopping terms is used to describe the Kitaev-like materials, which is calculated using the state-of-the-art fermionic projected entangled pair states (fPEPS) method. We find a QSL phase near the Mott insulator transition, which has a strong first-order transition to the semi-metal phase with the decrease of Hubbard $U$. We suggest that a promising routine to find the QSL is to find the Iridium oxides that are near the Mott insulator transitions.

Hybrid density functional (HDF) approximations usually deliver higher accuracy than local and semilocal approximations to the exchange-correlation functional, but this comes with drastically increased computational cost. Practical implementations of HDFs inevitably involve numerical approximations -- even more so than their local and semilocal counterparts due to the additional numerical complexity arising from treating the exact-exchange component. This raises the question regarding the reproducibility of the HDF results yielded by different implementations. In this work, we benchmark the numerical precision of four independent implementations of the popular Heyd-Scuseria-Ernzerhof (HSE) range-separated HDF on describing key materials' properties, including both properties derived from equations of states (EOS) and band gaps of 20 crystalline solids. We find that the energy band gaps obtained by the four codes agree with each other rather satisfactorily. However, for lattice constants and bulk moduli, the deviations between the results computed by different codes are of the same order of magnitude as the deviations between the computational and experimental results. On the one hand, this means that the HSE functional is rather accurate for describing the cohesive properties of simple insulating solids. On the other hand, this also suggests that the numerical precision achieved with current major HSE implementation is not sufficiently high to unambiguously assess the physical accuracy of HDFs. It is found that the pseudopotential treatment of the core electrons is a major factor that contributes to this uncertainty., Comment: 25 pages, 4 figures, 1 table

We investigate the lead-free halide double perovskites (HDPs) Cs$ _2BB'X_6$ ($B$=Ag, Na; $B'$=In, Bi; $X$=Cl, Br) via first-principles calculations. We find that both the geometric and electric structures of the HDPs obtained by the Heyd-Scuseria-Ernzerhof (HSE) hybrid functional are much better than those of the Perdew-Burke-Ernzerhof (PBE) functional. Importantly, we find that the electronic structures of DHPs are very sensitive to their geometries, especially the $B$-$X$ bond lengths. As a consequence, the electronic structures calculated by the HSE functional using the PBE optimized geometries may still significantly underestimate the band gaps, whereas the calculations on the HSE optimized geometries provide much more satisfactory results. The sensitivity of the band gaps of the DHPs to their geometries opens a promising path for the band structure engineering via doping and alloying. This work therefore provides an useful guideline for further improvement of HDPs materials.

The shift current effect, in materials lacking inversion symmetry, may potentially allow the performance of photovoltaics to surpass the Shockley-Queisser limit for traditional p-n junction-based photovoltaics. Although the shift-current effect has been studied from first-principles via second-order perturbation theory, an understanding of the dynamics of hot carriers is still lacking. We investigate the dynamics of the shift current in monolayer WS$_{2}$ via real-time propagation time-dependent density functional theory (rt-TDDFT). We find that the shift current can be generated within 10 - 20 fs after turning on the lights and dissipates within approximately a few tens of femtoseconds after turning off the lights. This property can be used for ultrafast photon detection. This work provides an important step toward understanding the dynamics of shift-current effects, which is crucial for device applications.

Characterizing an isolated attosecond pulse (IAP) is essential for its potential applications. A complete characterization of an IAP ultimately requires the determination of its electric field in both time and space domains. However, previous methods, like the widely-used RABBITT and attosecond streaking, only measure the temporal profile of the attosecond pulse. Here we demonstrate an all-optical method for the measurement of the space-time properties of an IAP. By introducing a non-collinear perturbing pulse to the driving field, the process of IAP generation is modified both spatially and temporally, manifesting as a spatial and a frequency modulation in the harmonic spectrum. By using a FROG-like retrieval method, the spatio-spectral phases of the harmonic spectrum are faithfully extracted from the induced spatio-spectral modulations, which allows a thoroughgoing characterization of the IAP in both time and space. With this method, the spatio-temporal structures of the IAP generated in a two-color driving field in both the near- and far-field are fully reconstructed, from which a weak spatio-temporal coupling in the IAP generation is revealed. Our approach overcomes the limitation in the temporal measurement in conventional in situ scheme, providing a reliable and holistic metrology for IAP characterization., Comment: 18 pages, 5 figures

We perform high-throughput first-principles computations to search the high Curie temperature ($T_{\rm C}$) two-dimensional ferromagnetic (2DFM) materials. We identify 79 2DFM materials and calculate their $T_{\rm C}$, in which Co$_2$F$_2$ has the highest $T_{\rm C}$=541K, well above the room temperature. The 79 2DFM materials are classified into different structural prototypes according to their structural similarity. We perform sure independence screening and sparsifying operator (SISSO) analysis to explore the relation between $T_{\rm C}$ and the material structures. The results suggest that the 2DFM materials with shorter distance between the magnetic atoms, larger local magnetic moments and more neighboring magnetic atoms are more likely to have higher $T_{\rm C}$.

Whether noisy quantum devices without error correction can provide quantum advantage over classical computers is a critical issue of current quantum computation. In this work, the random quantum circuits, which are used as the paradigm model to demonstrate quantum advantage, are simulated with depolarizing noise on experiment relevant two-dimensional architecture. With comprehensive numerical simulation and theoretical analysis, we find that the maximum achievable operator entanglement entropy, which indicates maximal simulation cost, has area law scaling with the system size for constant noise rate. On the other hand, we also find that the maximum achievable operator entanglement entropy has power law scaling with the noise rate for fixed system size, and the volume law scaling can be obtained only if the noise rate decreases when system size increase.

The combination of density functional theory with dynamical mean-field theory (DFT+DMFT) has become a powerful first-principles approach to tackle strongly correlated materials in condensed matter physics. The wide use of this approach relies on robust and easy-to-use implementations, and its implementation in various numerical frameworks will increase its applicability on the one hand and help crosscheck the validity of the obtained results on the other. In the work, we develop a formalism within the linear combination of numerical atomic orbital (NAO) basis set framework, which allows for merging NAO-based DFT codes with DMFT quantum impurity solvers. The formalism is implemented by interfacing two NAO-based DFT codes with three DMFT impurity solvers, and its validity is testified by benchmark calculations for a wide range of strongly correlated materials, including 3\textit{d} transition metal compounds, lanthanides, and actinides. Our work not only enables DFT+DMFT calculations using popular and rapidly developing NAO-based DFT code packages, but also facilitates the combination of more advanced beyond-DFT methodologies available in this codes with the DMFT machinery., Comment: 18 pages, 11 figures

For decades, people are developing efficient numerical methods for solving the challenging quantum many-body problem, whose Hilbert space grows exponentially with the size of the problem. However, this journey is far from over, as previous methods all have serious limitations. The recently developed deep learning methods provide a very promising new route to solve the long-standing quantum many-body problems. We report that a deep learning based simulation protocol can achieve the solution with state-of-the-art precision in the Hilbert space as large as $2^{1296}$ for spin system and $3^{144}$ for fermion system , using a HPC-AI hybrid framework on the new Sunway supercomputer. With highly scalability up to 40 million heterogeneous cores, our applications have measured 94% weak scaling efficiency and 72% strong scaling efficiency. The accomplishment of this work opens the door to simulate spin models and Fermion models on unprecedented lattice size with extreme high precision., Comment: Massive ground state optimizations of CNN-based wave-functions for $J1$-$J2$ model and $t$-$J$ model carried out on a heterogeneous cores supercomputer

One of the most remarkable characteristics of topological materials is that they have special surface states, which are determined by the topological properties of their bulk materials. The angle resolved photoemission spectroscopy (ARPES) is a powerful tool to explore the surface states, which allows to further investigate the topological phase transitions. However, it is very difficult to compare the first-principle calculated band structures to the ARPES results,when the systems are doped or alloyed, because the band structures are heavily folded.We develop an efficient band unfolding method based on numerical atomic orbitals (NAOs). We apply this method to study the surface states of the non-magnetically and magnetically doped topological insulators Bi$_2$Se$_3$ and the topological crystalline insulators Pb$_{1-x}$Sn$_{x}$Te.

We present a formulation and implementation of the DFT+\textit{U} method within the framework of linear combination of numerical atomic orbitals (NAO). Our implementation not only enables single-point total energy and electronic-structure calculations but also provides access to atomic forces and stresses, hence allowing for full structure relaxations of periodic systems. Furthermore, our implementation allows one to deal with non-collinear spin texture, with the spin-orbit coupling (SOC) effect treated self-consistently. The key aspect behind our implementation is a suitable definition of the correlated subspace when multiple atomic orbitals with the same angular momentum are used, and this is addressed via the "Mulliken charge projector" constructed in terms of the first (most localized) atomic orbital within the $d/f$ angular momentum channel. The important Hubbard $U$ and Hund $J$ parameters can be estimated from a screened Coulomb potential of the Yukawa type, with the screening parameter either chosen semi-empirically or determined from the Thomas-Fermi screening model. Benchmark calculations are performed for four late transition metal monoxide bulk systems, i.e., MnO, FeO, CoO, and NiO, and for the 5$d$-electron compounds IrO$_2$. For the former type of systems, we check the performance of our DFT+$U$ implementation for calculating band gaps, magnetic moments, electronic band structures, as well as forces and stresses; for the latter, the efficacy of our DFT+$U$+SOC implementation is assessed. Systematic comparisons with available experimental results, and especially with the results from other implementation schemes are carried out, which demonstrate the validity of our NAO-based DFT+$U$ formalism and implementation., Comment: 19 pages, 7 figures, 7 tables

Efficient numerical methods are promising tools for delivering unique insights into the fascinating properties of physics, such as the highly frustrated quantum many-body systems. However, the computational complexity of obtaining the wave functions for accurately describing the quantum states increases exponentially with respect to particle number. Here we present a novel convolutional neural network (CNN) for simulating the two-dimensional highly frustrated spin-$1/2$ $J_1-J_2$ Heisenberg model, meanwhile the simulation is performed at an extreme scale system with low cost and high scalability. By ingenious employment of transfer learning and CNN's translational invariance, we successfully investigate the quantum system with the lattice size up to $24\times24$, within 30 million cores of the new generation of sunway supercomputer. The final achievement demonstrates the effectiveness of CNN-based representation of quantum-state and brings the state-of-the-art record up to a brand-new level from both aspects of remarkable accuracy and unprecedented scales., Comment: This paper reviews the technical details of calculating spin-1/2 J1-J2 model using convolutional neural network ansatz on Sunway Supercomputer

The anisotropic-exchange spin-1/2 model on a triangular lattice has been used to describe the rare-earth chalcogenides, which may have exotic ground states. We investigate the quantum phase diagram of the model by using the projected entangled pair state (PEPS) method, and compare it to the classical phase diagram. Besides two stripe-ordered phase, and the 120$^\circ$ state, there is also a multi-\textbf{Q} phase. We identify the multi-\textbf{Q} phase as a $Z_{2}$ vortex state. No quantum spin liquid state is found in the phase diagram, contrary to the previous DMRG calculations.

Algorithms to simulate the ring-exchange models using the projected entangled pair states (PEPS) are developed. We generalize the imaginary time evolution (ITE) method to optimize PEPS wave functions for the models with ring-exchange interactions. We compare the effects of different approximations to the environment. To understand the numerical instability during the optimization, we introduce the ``singularity'' of a PEPS and develop a regulation procedure that can effectively reduce the singularity of a PEPS. We benchmark our method with the toric code model, and obtain extremely accurate ground state energies and topological entanglement entropy. We also benchmark our method with the two-dimensional cyclic ring exchange model, and find that the ground state has a strong vector chiral order. The algorithms can be a powerful tool to investigate the models with ring interactions. The methods developed in this work, e.g., the regularization process to reduce the singularity can also be applied to other models.

We present a derivation of the full formula to calculate the Berry curvature on non-orthogonal numerical atomic orbital (NAO) bases.Because usually, the number of NAOs is larger than that of the Wannier bases, we use a orbital contraction method to reduce the basis sizes, which can greatly improve the calculation efficiency without significantly reducing the calculation accuracy. We benchmark the formula by calculating the Berry curvature of ferroelectric BaTiO$_3$ and bcc Fe, as well as the anomalous Hall conductivity (AHC) for Fe. The results are in excellent agreement with the finite difference and previous results in the literature. We find that there are corrections terms to the Kubo formula of the Berry curvature. For the full NAO base, the differences between the two methods are negligibly small, but for the reduced bases sets, the correction terms become larger, which may not be neglected in some cases. The formula developed in this work can readily be applied to the non-orthogonal generalized Wannier functions.

We investigate the magnetic interactions in triangular rare-earth delafossites materials NaYbO$_2$ and NaYbS$_2$ via first-principles calculations. The calculated Curie-Weiss temperatures are in good agreement with experiments. We perform classical Monte Carlo simulations of the two compounds using the extracted exchange parameters. We find that if only the nearest neighbor interactions are considered, the magnetic ground states of NaYbO$_2$ and NaYbS$_2$ are a stripe and a planar 120\degree~ N\'{e}el state, respectively. The simulated transition temperatures are much higher than the lowest experimental temperatures, where no magnetic ordering was observed. However, we show by adding suitable second neighbor interactions, the {\it classical} magnetic ground state of NaYbO$_2$ becomes to the $Z_2$ vortex phase, and the simulated specific heat $C_v$ are very similar to the experimental observations, with no obvious phase transition down to the extremely low temperature.

We experimentally explore the fingerprint of the microscopic electron dynamics in second-order harmonic generation (SHG). It is shown that the interbond electron hopping induces a novel source of nonlinear polarization and plays an important role even when the driving laser intensity is 2 orders of magnitude lower than the characteristic atomic field. Our model predicts anomalous anisotropic structures of the SHG yield contributed by the interbond electron hopping, which is identified in our experiments with ZnO crystals. Moreover, a generalized second-order susceptibility with an explicit form is proposed, which provides a unified description in both the weak and strong field regimes. Our work reveals the nonlinear responses of materials at the electron scale and extends the nonlinear optics to a previously unexplored regime, where the nonlinearity related to the interbond electron hopping becomes dominant. It paves the way for realizing controllable nonlinearity on an ultrafast time scale.

Neural networks have been used as variational wave functions for quantum many-particle problems. It has been shown that the correct sign structure is crucial to obtain the high accurate ground state energies. In this work, we propose a hybrid wave function combining the convolutional neural network (CNN) and projected entangled pair states (PEPS), in which the sign structures are determined by the PEPS, and the amplitudes of the wave functions are provided by CNN. We benchmark the ansatz on the highly frustrated spin-1/2 $J_1$-$J_2$ model. We show that the achieved ground energies are competitive to state-of-the-art results.

We present an efficient, linear-scaling implementation for building the (screened) Hartree-Fock exchange (HFX) matrix for periodic systems within the framework of numerical atomic orbital (NAO) basis functions. Our implementation is based on the localized resolution of the identity approximation by which two-electron Coulomb repulsion integrals can be obtained by only computing two-center quantities -- a feature that is highly beneficial to NAOs. By exploiting the locality of basis functions and efficient prescreening of the intermediate three- and two-index tensors, one can achieve a linear scaling of the computational cost for building the HFX matrix with respect to the system size. Our implementation is massively parallel, thanks to a MPI/OpenMP hybrid parallelization strategy for distributing the computational load and memory storage. All these factors add together to enable highly efficient hybrid functional calculations for large-scale periodic systems. In this work we describe the key algorithms and implementation details for the HFX build as implemented in the ABACUS code package. The performance and scalability of our implementation with respect to the system size and the number of CPU cores are demonstrated for selected benchmark systems up to 4096 atoms.

Intrinsic magnetic van der Waals (vdW) materials have attracted much attention, especially Fe$_{3}$GeTe$_{2}$ (FGT), which exhibits highly tunable properties. However, despite vast efforts, there are still several challenging issues to be resolved. Here using a first-principles linear-response approach, we carry out comprehensive investigation of both bulk and monolayer FGT. We find that although the magnetic exchange interactions in FGT are frustrated, our Monte Carlo simulations agree with the total energy calculations, confirming that the ground state of bulk FGT is indeed ferromagnetic (FM). A tiny electron doping reduces the magnetic frustration, resulting in significant increasing of the Curie temperature. We also calculate the spin-wave dispersion, and find a small spin gap as well as a nearly flat band in the magnon spectra. These features can be used to compare with the future neutron scattering measurement and finally clarify the microscopic magnetic mechanism in this two-dimensional family materials.

In strong laser fields, the electron density in solids can show up tidal motions caused by the laser ``pulling'' on the electrons, which forms an electron density tide (EDT). However, the strong field processes in solids are always explained by the single-active-electron (SAE) model and the fluctuation of background electrons is neglected previously. Here, we demonstrate the strong field induced EDT effect and propose a model for revealing its role in high harmoinc generation (HHG) from solids. We show that the EDT effect induces an additional polarization current beyond the SAE approximation. It gives a new mechanism for HHG and leads to new anisotropic structures, which are experimentally observed with MgO. Our experiment indicates that the EDT effect becomes more obvious with increasing the laser intensity. Our work establishes the bridge between the HHG and the microscopic EDT in solids induced by the strong laser field, which paves the way to probe the ultrafast electron dynamics.

Understanding the retention and recycling of hydrogen isotopes in liquid metal plasma-facing materials such as liquid Li, Sn, and Li-Sn are of fundamental importance in designing magnetically confined fusion reactors. We perform first-principles molecules dynamics simulations of liquid Li-Sn slab with inserted D atoms to provide microscopic insights into the interactions of D with Li-Sn liquid metal. We observe evaporation of D$_2$ and LiD molecules out of the Li-Sn slabs. With detailed analysis, we unveil a cooperative process of forming D$_2$ molecules in liquid Li-Sn, where Li atoms act as catalytic centers to trap a D atom before another D comes nearby to form a molecule, and the surplus charges are transferred from D$_2$ to nearby Sn atoms. Furthermore, we predict a temperature window in which D$_2$ molecules can escape to vacuum, while LiD molecules cannot. The above findings deepen our understanding of interactions between hydrogen isotopes and Li-Sn liquid metal.