Publications

We present the Julia package PauliStrings ( this https URL ) for quantum many-body simulations, which performs fast operations on the Pauli group by encoding Pauli strings in binary. All of the Pauli string algebra is encoded into low-level logic operations on integers, and is made efficient by various truncation methods which allow for systematic extrapolation of the results. We illustrate the effectiveness of our package by (i) performing Heisenberg time evolution through direct numerical integration and (ii) by constructing a Liouvillian Krylov space. We benchmark the results against tensor network methods, and we find our package performs favorably. In addition, we show that this representation allows for easy encoding of any geometry. We present results for chaotic and integrable spin systems in 1D as well as some examples in 2D. Currently, the main limitations are the inefficiency of representing non-trivial pure states (or other low-rank operators), as well as the need to introduce dissipation to probe long-time dynamics.

'Strange' metals that do not follow the predictions of Fermi liquid theory are prevalent in materials that feature superconductivity arising from electron interactions. In recent years, it has been hypothesized that spatial randomness in electron interactions must play a crucial role in strange metals for their hallmark linear-in-temperature (T) resistivity to survive down to low temperatures where phonon and Umklapp processes are ineffective, as is observed in experiments. However, a clear picture of how this happens has not yet been provided in a realistic model free from artificial constructions such as large-N limits and replica tricks. We study a realistic model of two-dimensional metals with spatially random antiferromagnetic interactions in a non-perturbative regime, using numerically exact high-performance large-scale hybrid Monte Carlo and exact averages over the quenched spatial randomness. Our simulations reproduce strange metals' key experimental signature of linear-in-T resistivity with a 'planckian' transport scattering rate Γ

We simulate high-pressure hydrogen in its liquid phase close to molecular dissociation using a machine-learned interatomic potential. The model is trained with density functional theory (DFT) forces and energies, with the Perdew-Burke-Ernzerhof (PBE) exchange-correlation functional. We show that an accurate NequIP model, an E(3)-equivariant neural network potential, accurately reproduces the phase transition present in PBE. Moreover, the computational efficiency of this model allows for substantially longer molecular dynamics trajectories, enabling us to perform a finite-size scaling (FSS) analysis to distinguish between a crossover and a true first-order phase transition. We locate the critical point of this transition, the liquid-liquid phase transition (LLPT), at 1200-1300 K and 155-160 GPa, a temperature lower than most previous estimates and close to the melting transition.

Thanks to its low or negative surface electron affinity and chemical inertness, diamond is attracting broad attention as a source material of solvated electrons produced by optical excitation of the solid–liquid interface. Unfortunately, its wide bandgap typically imposes the use of wavelengths in the ultraviolet range, hence complicating practical applications. Here, we probe the photocurrent response of water surrounded by single-crystal diamond surfaces engineered to host shallow nitrogen-vacancy (NV) centers. We observe clear signatures of diamond-induced photocurrent generation throughout the visible range and for wavelengths reaching up to 594 nm. Experiments as a function of laser power suggest that NV centers and other coexisting defects─likely in the form of surface traps─contribute to carrier injection, though we find that NVs dominate the system response in the limit of high illumination intensities. Given our growing understanding of near-surface NV centers and adjacent point defects, these results open new perspectives in the application of diamond–liquid interfaces to photocarrier-initiated chemical and spin processes in fluids.

Identifying the rank of species in a complex ecosystem is a difficult task, since the rank of each species invariably depends on the interactions stipulated with other species through the adjacency matrix of the network. A common ranking method in economic and ecological networks is to sort the nodes such that the layout of the reordered adjacency matrix looks maximally nested with all nonzero entries packed in the upper left corner, called Nestedness Maximization Problem (NMP). Here we solve this problem by defining a suitable cost-energy function for the NMP which reveals the equivalence between the NMP and the Quadratic Assignment Problem, one of the most important combinatorial optimization problems, and use statistical physics techniques to derive a set of self-consistent equations whose fixed point represents the optimal nodes’ rankings in an arbitrary bipartite mutualistic network. Concurrently, we present an efficient algorithm to solve the NMP that outperforms state-of-the-art network-based metrics and genetic algorithms. Eventually, our theoretical framework may be easily generalized to study the relationship between ranking and network structure beyond pairwise interactions, e.g. in higher-order networks.

Optically bright emitters in hexagonal boron nitride (hBN) often acting as a source of a single-photon are mostly attributed to point-defect centers, featuring localized intra-bandgap electronic states. Although vacancies, anti-sites, and impurities have been proposed as candidates, the exact physical and chemical nature of most hBN single-photon emitters (SPEs) within the visible region are still up for debate. Combining site-specific high-angle annular dark-field imaging (HAADF) with electron energy loss spectroscopy (EELS), we resolve and identify a few carbon substitutions among neighboring hBN hexagons, all within the same sample region, from which typical defect emission is observed. Our experimental results are further supported by first-principles calculations, through which the stability and possible optical transitions of the proposed carbon-defect complex are assessed. The presented correlation between optical emission and defects provides valuable information toward the controlled creation of emitters in hBN, highlighting carbon complexes as another probable cause of its visible SPEs.

Interactions between light and matter allow the realization of out-of-equilibrium states in quantum solids. In particular, nonlinear phononics is one of the efficient approaches to realizing the stationary electronic state in non-equilibrium. Herein, by using extended ab initio molecular dynamics, we identify that long-lived light-driven quasi-stationary geometry could stabilize the topological nature in the material family of HgTe compounds. We show that coherent excitation of the infrared-active phonon mode results in a distortion of the atomic geometry with a lifetime of several picoseconds. We show that four Weyl points are located exactly at the Fermi level in this non-equilibrium geometry, making it an ideal long-lived metastable Weyl semimetal. We propose that such a metastable topological phase can be identified by photoelectron spectroscopy of the Fermi arc surface states or ultrafast pump-probe transport measurements of the nonlinear Hall effect.

We theoretically investigate high-order harmonic generation (HHG) in graphene under mid-infrared (MIR) and terahertz (THz) fields based on a quantum master equation. Numerical simulations show that MIR-induced HHG in graphene can be enhanced by a factor of 10 for fifth harmonic and a factor of 25 for seventh harmonic under a THz field with a peak strength of 0.5 MV/cm by optimizing the relative angle between the MIR and THz fields. To identify the origin of this enhancement, we compare the fully dynamical calculations with a simple thermodynamic model and a nonequilibrium population model. The analysis shows that the enhancement of the high-order harmonics mainly results from a coherent coupling between MIR- and THz-induced transitions that goes beyond a simple THz-induced population contribution.

We develop new sub-optimality bounds for gradient descent (GD) that depend on the conditioning of the objective along the path of optimization, rather than on global, worst-case constants. Key to our proofs is directional smoothness, a measure of gradient variation that we use to develop upper-bounds on the objective. Minimizing these upper-bounds requires solving implicit equations to obtain a sequence of strongly adapted step-sizes; we show that these equations are straightforward to solve for convex quadratics and lead to new guarantees for two classical step-sizes. For general functions, we prove that the Polyak step-size and normalized GD obtain fast, path-dependent rates despite using no knowledge of the directional smoothness. Experiments on logistic regression show our convergence guarantees are tighter than the classical theory based on -smoothness.

One of the core problems in mean-field control and mean-field games is to solve the corresponding McKean–Vlasov forward-backward stochastic differential equations (MV-FBSDEs). Most existing methods are tailored to special cases in which the mean-field interaction only depends on expectation or other moments and thus are inadequate to solve problems when the mean-field interaction has full distribution dependence. In this paper, we propose a novel deep learning method for computing MV-FBSDEs with a general form of mean-field interactions. Specifically, built on fictitious play, we recast the problem into repeatedly solving standard FBSDEs with explicit coefficient functions. These coefficient functions are used to approximate the MV-FBSDEs’ model coefficients with full distribution dependence, and are updated by solving another supervising learning problem using training data simulated from the last iteration’s FBSDE solutions. We use deep neural networks to solve standard BSDEs and approximate coefficient functions in order to solve high-dimensional MV-FBSDEs. Under proper assumptions on the learned functions, we prove that the convergence of the proposed method is free of the curse of dimensionality (CoD) by using a class of integral probability metrics previously developed in [J. Han, R. Hu, and J. Long, Stochastic Process. Appl., 164 (2023), pp. 242–287]. The proved theorem shows the advantage of the method in high dimensions. We present the numerical performance in high-dimensional MV-FBSDE problems, including a mean-field game example of the well-known Cucker–Smale model, the cost of which depends on the full distribution of the forward process.