Publications

We solve the fully connected ferromagnetic XY model in a local Gaussian random magnetic field using the replica method. We work out the replica symmetric (RS) solution and prove its exactness by diagonalizing the Hessian around the RS saddle point and showing that all its eigenvalues are non-negative for all values of the temperature T and random field strength σ. As a consequence, there is no spin-glass phase with replica symmetry breaking in this model. Moreover, we show that the critical line σc(T) separating the paramagnetic and ferromagnetic phases is a strictly monotonic decreasing function of T and thus there is no reentrant paramagnetic phase at low temperature. Furthermore, all points on the critical line are second order phase transition critical points, and thus there are no first-order phase transitions for any value of σ and T. Finally, we show that the calculation we present here can be easily generalized beyond the XY model to the O(ν) model with an arbitrary number ν of spin components.

Correlated sampling has wide-ranging applications in Monte Carlo calculations. When branching random walks are involved, as commonly found in many algorithms in quantum physics and electronic structure, population control is typically not applied with correlated sampling due to technical difficulties. This hinders the stability and efficiency of correlated sampling. In this work, we study schemes for allowing birth/death in correlated sampling and propose an algorithm for population control. The algorithm can be realized in several variants depending on the application. One variant is a static method that creates a reference run and allows other correlated calculations to be added a posteriori. Another optimizes the population control for a set of correlated, concurrent runs dynamically. These approaches are tested in different applications in quantum systems, including both the Hubbard model and electronic structure calculations in real materials.

Click chemistry, which refers to chemical reactions that are fast, selective, and with high product yields, has become a powerful approach in organic synthesis and chemical biology. Due to the cytotoxicity of the transition metals employed in click chemistry reactions, a search for novel metal-free alternatives continues. Herein we demonstrate that an optical cavity can be utilized as a metal-free alternative in click chemistry cycloaddition reaction between cyanoacetylene and formylazide using the quantum electrodynamics coupled cluster (QED-CC) method. We show that by changing the molecular orientation with respect to the polarization of the cavity mode(s), the reaction can be selectively catalyzed to form a major 1,4-disubstituted or 1,5-disubstituted product. This work highlights that a cavity has the same effect on the investigated cycloaddition as the transition metal catalysts traditionally employed in click chemistry reactions. We expect our findings to further stimulate research in cavity-assisted click chemistry reactions.

We study the quantum Hall effect in a two-dimensional homogeneous electron gas coupled to a quantum cavity field. As initially pointed out by Kohn, Galilean invariance for a homogeneous quantum Hall system implies that the electronic center of mass (CM) decouples from the electron-electron interaction, and the energy of the CM mode, also known as Kohn mode, is equal to the single particle cyclotron transition. In this work, we point out that strong light-matter hybridization between the Kohn mode and the cavity photons gives rise to collective hybrid modes between the Landau levels and the photons. We provide the exact solution for the collective Landau polaritons and we demonstrate the weakening of topological protection at zero temperature due to the existence of the lower polariton mode which is softer than the Kohn mode. This provides an intrinsic mechanism for the recently observed topological breakdown of the quantum Hall effect in a cavity [Appugliese et al., Science 375, 1030-1034 (2022)]. Importantly, our theory predicts the cavity suppression of the thermal activation gap in the quantum Hall transport. Our work paves the way for future developments in the cavity control of quantum materials.

We investigate optical nonlinearities that are induced and enhanced due to the strong phonon resonance in hexagonal boron nitride. We predict and observe large sub-picosecond duration signals due to four-wave mixing (FWM) during resonant excitation. The resulting FWM signal allows for time-resolved observation of the crystal motion. In addition, we observe enhancements of third-harmonic generation with resonant pumping at the hBN transverse optical phonon. Phonon-induced nonlinear enhancements are also predicted to yield large increases in high-harmonic efficiencies.

Many stellarator coil design problems are plagued by multiple minima, where the locally optimal coil sets can sometimes vary substantially in performance. As a result, solving a coil design problem a single time with a local optimization algorithm is usually insufficient and better optima likely do exist. To address this problem, we propose a global optimization algorithm for the design of stellarator coils and outline how to apply box constraints to the physical positions of the coils. The algorithm has a global exploration phase that searches for interesting regions of design space and is followed by three local optimization algorithms that search in these interesting regions (a "global-to-local" approach). The first local algorithm (phase I), following the globalization phase, is based on near-axis expansions and finds stellarator coils that optimize for quasisymmetry in the neighborhood of a magnetic axis. The second local algorithm (phase II) takes these coil sets and optimizes them for nested flux surfaces and quasisymmetry on a toroidal volume. The final local algorithm (phase III) polishes these configurations for an accurate approximation of quasisymmetry. Using our global algorithm, we study the trade-off between coil length, aspect ratio, rotational transform, and quality of quasi-axisymmetry. The database of stellarators, which comprises almost 140,000 coil sets, is available online and is called QUASR, for "QUAsi-symmetric Stellarator Repository".

Inferring dynamical models from low-resolution temporal data continues to be a significant challenge in biophysics, especially within transcriptomics, where separating molecular programs from noise remains an important open problem. We explore a common scenario in which we have access to an adequate amount of cross-sectional samples at a few time-points, and assume that our samples are generated from a latent diffusion process. We propose an approach that relies on the probability flow associated with an underlying diffusion process to infer an autonomous, nonlinear force field interpolating between the distributions. Given a prior on the noise model, we employ score-matching to differentiate the force field from the intrinsic noise. Using relevant biophysical examples, we demonstrate that our approach can extract non-conservative forces from non-stationary data, that it learns equilibrium dynamics when applied to steady-state data, and that it can do so with both additive and multiplicative noise models.

The quantum Fourier transform (QFT) is a key component of many important quantum algorithms, most famously being the essential ingredient in Shor’s algorithm for factoring products of primes. Given its remarkable capability, one would think it can introduce large entanglement to qubit systems and would be difficult to simulate classically. While early results showed the QFT indeed has maximal operator entanglement, we show that this is entirely due to the bit reversal in the QFT. The core part of the QFT has Schmidt coefficients decaying exponentially quickly, and thus it can generate only a constant amount of entanglement regardless of the number of qubits. In addition, we show the entangling power of the QFT is the same as the time evolution of a Hamiltonian with exponentially decaying interactions, and thus a variant of the area law for dynamics can be used to understand the low entanglement intuitively. Using the low entanglement property of the QFT, we show that classical simulations of the QFT on a matrix product state with low bond dimension take time linear in the number of qubits, providing a potential speedup over the classical fast Fourier transform on many classes of functions. We demonstrate this speedup in test calculations on some simple functions. For data vectors of length 10^6 – 10^8, the speedup can be a few orders of magnitude.

To check the accuracy of Bayesian computations, it is common to use rank-based simulation-based calibration (SBC). However, SBC has drawbacks: The test statistic is somewhat ad-hoc, interactions are difficult to examine, multiple testing is a challenge, and the resulting p-value is not a divergence metric. We propose to replace the marginal rank test with a flexible classification approach that learns test statistics from data. This measure typically has a higher statistical power than the SBC rank test and returns an interpretable divergence measure of miscalibration, computed from classification accuracy. This approach can be used with different data generating processes to address likelihood-free inference or traditional inference methods like Markov chain Monte Carlo or variational inference. We illustrate an automated implementation using neural networks and statistically-inspired features, and validate the method with numerical and real data experiments.