Publications

In a previous work [Phys. Rev. Lett. 123, 047202 (2019)] a translationally invariant framework called quantum-electrodynamical Bloch (QED-Bloch) theory was introduced for the description of periodic materials in homogeneous magnetic fields and strongly coupled to the quantized photon field in the optical limit. For such systems, we show that QED-Bloch theory predicts the existence of fractal polaritonic spectra as a function of the cavity coupling strength. In addition, for the energy spectrum as a function of the relative magnetic flux we find that a terahertz cavity can modify the standard Hofstadter butterfly. In the limit of no quantized photon field, QED-Bloch theory captures the well-known fractal spectrum of the Hofstadter butterfly and can be used for the description of 2D materials in strong magnetic fields, which are of great experimental interest. As a further application, we consider Landau levels under cavity confinement and show that the cavity alters the quantized Hall conductance and that the Hall plateaus are modified as σ

A deeper theoretical understanding of driven-dissipative interacting systems and their nonequilibrium phase transitions is essential both to advance our fundamental physics understanding and to harness technological opportunities arising from optically controlled quantum many-body states. This paper provides a numerical study of dynamical phases and the transitions between them in the nonequilibrium steady state of the prototypical two-dimensional Heisenberg antiferromagnet with drive and dissipation. We demonstrate a nonthermal transition that is characterized by a qualitative change in the magnon distribution, from subthermal at low drive to a generalized Bose-Einstein form including a nonvanishing condensate fraction at high drive. A finite-size analysis reveals static and dynamical critical scaling at the transition, with a discontinuous slope of the magnon number versus driving field strength and critical slowing down at the transition point. Implications for experiments on quantum materials and polariton condensates are discussed.

Interaction of light with an excited free-electron gas is a fundamental process spanning a large variety of fields in physics. The advent of femtosecond laser pulses and extreme-ultraviolet sources allowed one to put theoretical models to the test. Recent experimental and theoretical investigations of nonequilibrium aluminum, which is considered to be a good real-world representation of an ideal free-electron metal, showed that, despite significant progress, the transient hot-electron/cold-ion state is not well understood. In particular, the role of plasmon broadening, screening, and electron degeneracy remains unclear. Here, we experimentally investigate the free-free opacity in aluminum on the few-femtosecond timescale at laser intensities close to the damage threshold. Few-femtosecond time resolution allows us to track the purely electronic contribution to nonequilibrium absorption and unambiguously separate it from the slower lattice contribution. We support the experiments with ab initio calculations and a nearly free electron model in the Sommerfeld expansion. We find that the simplest independent-particle model with a fixed band structure is sufficient to explain the experimental findings without the need to include changes in screening or electron scattering, contrasting previous observations in 3d transition metals. We further find that electronic heating of a free-electron gas shifts the spectral weight of the absorption to higher photon energies, and we are able to distinguish the influence of the population change and the chemical potential shift based on the comparison of ab initio calculations to a simplified free-electron model. Our findings provide a benchmark for further investigations and modeling of dense nonequilibrium plasma under even more extreme conditions.

Physical theories that depend on many parameters or are tested against data from many different experiments pose unique challenges to statistical inference. Many models in particle physics, astrophysics and cosmology fall into one or both of these categories. These issues are often sidestepped with statistically unsound ad hoc methods, involving intersection of parameter intervals estimated by multiple experiments, and random or grid sampling of model parameters. Whilst these methods are easy to apply, they exhibit pathologies even in low-dimensional parameter spaces, and quickly become problematic to use and interpret in higher dimensions. In this article we give clear guidance for going beyond these procedures, suggesting where possible simple methods for performing statistically sound inference, and recommendations of readily-available software tools and standards that can assist in doing so. Our aim is to provide any physicists lacking comprehensive statistical training with recommendations for reaching correct scientific conclusions, with only a modest increase in analysis burden. Our examples can be reproduced with the code publicly available at Zenodo.

Epitopes that bind simultaneously to all human alleles of Major Histocompatibility Complex class II (MHC II) are considered one of the key factors for the development of improved vaccines and cancer immunotherapies. To engineer MHC II multiple-allele binders, we developed a protocol called PanMHC-PARCE, based on the unsupervised optimization of the epitope sequence by single-point mutations, parallel explicit-solvent molecular dynamics simulations and scoring of the MHC II-epitope complexes. The key idea is accepting mutations that not only improve the affinity but also reduce the affinity gap between the alleles. We applied this methodology to enhance a Plasmodium vivax epitope for multiple-allele binding. In vitro rate-binding assays showed that four engineered peptides were able to bind with improved affinity toward multiple human MHC II alleles. Moreover, we demonstrated that mice immunized with the peptides exhibited interferon-gamma cellular immune response. Overall, the method enables the engineering of peptides with improved binding properties that can be used for the generation of new immunotherapies.

The distance to the Galactic center R0 is a fundamental parameter for understanding the Milky Way, because all observations of our Galaxy are made from our heliocentric reference point. The uncertainty in R0 limits our knowledge of many aspects of the Milky Way, including its total mass and the relative mass of its major components, and any orbital parameters of stars employed in chemo-dynamical analyses. While measurements of R0 have been improving over a century, measurements in the past few years from a variety of methods still find a wide range of R0 being somewhere within 8.0 to 8.5kpc. The most precise measurements to date have to assume that Sgr A∗ is at rest at the Galactic center, which may not be the case. In this paper, we use maps of the kinematics of stars in the Galactic bar derived from APOGEE DR17 and Gaia EDR3 data augmented with spectro-photometric distances from the \texttt{astroNN} neural-network method. These maps clearly display the minimum in the rotational velocity vT and the quadrupolar signature in radial velocity vR expected for stars orbiting in a bar. From the minimum in vT, we measure R0=8.23±0.12kpc. We validate our measurement using realistic N-body simulations of the Milky Way. We further measure the pattern speed of the bar to be Ωbar=40.08±1.78kms−1kpc−1. Because the bar forms out of the disk, its center is manifestly the barycenter of the bar+disc system and our measurement is therefore the most robust and accurate measurement of R0 to date.

Parametric computer-aided design (CAD) tools are the predominant way that engineers specify physical structures, from bicycle pedals to airplanes to printed circuit boards. The key characteristic of parametric CAD is that design intent is encoded not only via geometric primitives, but also by parameterized constraints between the elements. This relational specification can be viewed as the construction of a constraint program, allowing edits to coherently propagate to other parts of the design. Machine learning offers the intriguing possibility of accelerating the design process via generative modeling of these structures, enabling new tools such as autocompletion, constraint inference, and conditional synthesis. In this work, we present such an approach to generative modeling of parametric CAD sketches, which constitute the basic computational building blocks of modern mechanical design. Our model, trained on real-world designs from the SketchGraphs dataset, autoregressively synthesizes sketches as sequences of primitives, with initial coordinates, and constraints that reference back to the sampled primitives. As samples from the model match the constraint graph representation used in standard CAD software, they may be directly imported, solved, and edited according to downstream design tasks. In addition, we condition the model on various contexts, including partial sketches (primers) and images of hand-drawn sketches. Evaluation of the proposed approach demonstrates its ability to synthesize realistic CAD sketches and its potential to aid the mechanical design workflow.

Observations of the Milky Way's low-α disk show that at fixed metallicity, [Fe/H], several element abundance, [X/Fe], correlate with age, with unique slopes and small scatters around the age-[X/Fe] relations. In this study, we turn to simulations to explore the age-[X/Fe] relations for the elements C, N, O, Mg, Si, S, and Ca that are traced in a FIRE-2 cosmological zoom-in simulation of a Milky Way-like galaxy, m12i, and understand what physical conditions give rise to the observed age-[X/Fe] trends. We first explore the distributions of mono-age populations in their birth and current locations, [Fe/H], and [X/Fe], and find evidence for inside-out radial growth for stars with ages < 7 Gyr. We then examine the age-[X/Fe] relations across m12i's disk and find that the direction of the trends agree with observations, apart from C, O, and Ca, with remarkably small intrinsic scatters, σint (0.01-0.04 dex). This σint measured in the simulations is also metallicity-dependent, with σint ≈ 0.025 dex at [Fe/H]=-0.25 dex versus σint ≈ 0.015 dex at [Fe/H]=0 dex, and a similar metallicity dependence is seen in the GALAH survey for the elements in common. Additionally, we find that σint is higher in the inner galaxy, where stars are older and formed in less chemically-homogeneous environments. The age-[X/Fe] relations and the small scatter around them indicate that simulations capture similar chemical enrichment variance as observed in the Milky Way, arising from stars sharing similar element abundances at a given birth place and time.

Uncertain feedback processes in galaxies affect the distribution of matter, currently limiting the power of weak lensing surveys. If we can identify cosmological statistics that are robust against these uncertainties, or constrain these effects by other means, then we can enhance the power of current and upcoming observations from weak lensing surveys such as DES, Euclid, the Rubin Observatory, and the Roman Space Telescope. In this work, we investigate the potential of the electron density auto-power spectrum as a robust probe of cosmology and baryonic feedback. We use a suite of (magneto-)hydrodynamic simulations from the CAMELS project and perform an idealized analysis to forecast statistical uncertainties on a limited set of cosmological and physically-motivated astrophysical parameters. We find that the electron number density auto-correlation, measurable through either kinematic Sunyaev-Zel'dovich observations or through Fast Radio Burst dispersion measures, provides tight constraints on Ωm and the mean baryon fraction in intermediate-mass halos, f¯bar. By obtaining an empirical measure for the associated systematic uncertainties, we find these constraints to be largely robust to differences in baryonic feedback models implemented in hydrodynamic simulations. We further discuss the main caveats associated with our analysis, and point out possible directions for future work.

Modern surveys of gravitational microlensing events have progressed to detecting thousands per year. Surveys are capable of probing Galactic structure, stellar evolution, lens populations, black hole physics, and the nature of dark matter. One of the key avenues for doing this is studying the microlensing Einstein radius crossing time distribution (tE). However, systematics in individual light curves as well as over-simplistic modeling can lead to biased results. To address this, we developed a model to simultaneously handle the microlensing parallax due to Earth's motion, systematic instrumental effects, and unlensed stellar variability with a Gaussian Process model. We used light curves for nearly 10,000 OGLE-III and IV Milky Way bulge microlensing events and fit each with our model. We also developed a forward model approach to infer the timescale distribution by forward modeling from the data rather than using point estimates from individual events. We find that modeling the variability in the baseline removes a source of significant bias in individual events, and previous analyses over-estimated the number of long timescale (tE>100 days) events due to their over simplistic models ignoring parallax effects and stellar variability. We use our fits to identify hundreds of events that are likely black holes.