Publications

Predictions of the mean and covariance matrix of summary statistics are critical for confronting cosmological theories with observations, not least for likelihood approximations and parameter inference. Accurate estimates require running costly N-body and hydrodynamics simulations. Approximate solvers, or surrogates, greatly reduce the computational cost but introduce biases, especially in the non-linear regime of structure growth. We propose ‘CARPool Bayes’ to solve the inference problem for both the means and covariances using a combination of simulations and surrogates. Our approach allows incorporating prior information for the mean and covariance. We derive closed-form solutions for maximum a posteriori covariance estimates that are efficient Bayesian shrinkage estimators, guarantee positive semidefiniteness, and can optionally leverage analytical covariance approximations. We discuss choices of the prior and propose a procedure for obtaining optimal prior hyperparameter values with a small set of test simulations. We test our method by estimating the covariances of clustering statistics of GADGET-IIIN-body simulations at redshift z = 0.5 using surrogates from a 100–1000× faster particle-mesh code. Taking the sample covariance from 15 000 simulations as the truth, and using an empirical Bayes prior with diagonal blocks, our estimator produces nearly identical Fisher matrix contours for ΛCDM parameters using only 15 simulations of the non-linear dark matter power spectrum. In this case, the number of simulations is so small that the sample covariance is degenerate. We show cases where even with a naïve prior our method improves the estimate. Our framework is applicable to a wide range of cosmological problems where fast surrogates are available.

We use angle-resolved photoemission to map the Fermi surface and quasiparticle dispersion of bulk-like thin films of SrMoO

We theoretically study the THz-induced high-order harmonic generation (HHG) and nonlinear electric transport in graphene based on the quantum master equation with the relaxation time approximation. To obtain microscopic insight into the phenomena, we compare the results of the fully dynamical calculations with those under a quasi-static approximation, where the electronic system is approximated as a nonequilibrium steady state. As a result, we find that the THz-induced electron dynamics in graphene can be accurately modeled with the nonequilibrium steady-state at each instance. The population distribution analysis further clarifies that the THz-induced HHG in graphene originates from the reduction of effective conductivity due to a large displacement of electrons in the Brillouin zone. By comparing the present nonequilibrium picture with a thermodynamic picture, we explore the role of the nonequilibrium nature of electron dynamics on the extremely nonlinear optical and transport phenomena in graphene.

Over a decade of work has culminated in the consensus that one-dimensional systems, subject to sufficiently large disorder, fail to thermalize and possess an extensive set of local integrals of motion. In this Letter, I will provide numerical evidence for the contrary.

The early dark energy (EDE) scenario aims to increase the value of the Hubble constant (H0) inferred from cosmic microwave background (CMB) data over that found in ΛCDM, via the introduction of a new form of energy density in the early universe. The EDE component briefly accelerates cosmic expansion just prior to recombination, which reduces the physical size of the sound horizon imprinted in the CMB. Previous work has found that non-zero EDE is not preferred by Planck CMB power spectrum data alone, which yield a 95% confidence level (CL) upper limit fEDE99.7% CL: fEDE=0.091+0.020−0.036, with H0=70.9+1.0−2.0 km/s/Mpc (both 68% CL). From a model-selection standpoint, we find that EDE is favored over ΛCDM by these data at roughly 3σ significance. In contrast, a joint analysis of the full Planck and ACT data yields no evidence for EDE, as previously found for Planck alone. We show that the preference for EDE in ACT alone is driven by its TE and EE power spectrum data. The tight constraint on EDE from Planck alone is driven by its high-ℓ TT power spectrum data. Understanding whether these differing constraints are physical in nature, due to systematics, or simply a rare statistical fluctuation is of high priority. The best-fit EDE models to ACT and Planck exhibit coherent differences across a wide range of multipoles in TE and EE, indicating that a powerful test of this scenario is anticipated with near-future data from ACT and other ground-based experiments.

Precise Gaia measurements of positions, parallaxes, and proper motions provide an opportunity to calculate 3D positions and 2D velocities (i.e., 5D phase-space) of Milky Way stars. Where available, spectroscopic radial velocity (RV) measurements provide full 6D phase-space information, however there are now and will remain many stars without RV measurements. Without an RV it is not possible to directly calculate 3D stellar velocities; however, one can infer 3D stellar velocities by marginalizing over the missing RV dimension. In this paper, we infer the 3D velocities of stars in the Kepler field in Cartesian Galactocentric coordinates (vx, vy, vz). We directly calculate velocities for around a quarter of all Kepler targets, using RV measurements available from the Gaia, LAMOST, and APOGEE spectroscopic surveys. Using the velocity distributions of these stars as our prior, we infer velocities for the remaining three quarters of the sample by marginalizing over the RV dimension. The median uncertainties on our inferred vx, vy, and vz velocities are around 4, 18, and 4 km s−1, respectively. We provide 3D velocities for a total of 148,590 stars in the Kepler field. These 3D velocities could enable kinematic age-dating, Milky Way stellar population studies, and other scientific studies using the benchmark sample of well-studied Kepler stars. Although the methodology used here is broadly applicable to targets across the sky, our prior is specifically constructed from and for the Kepler field. Care should be taken to use a suitable prior when extending this method to other parts of the Galaxy.

Polaritons in hyperbolic van der Waals materials---where principal axes have permittivities of opposite signs---are light-matter modes with unique properties and promising applications. Isofrequency contours of hyperbolic polaritons may undergo topological transitions from open hyperbolas to closed ellipse-like curves, prompting an abrupt change in physical properties. Electronically-tunable topological transitions are especially desirable for future integrated technologies but have yet to be demonstrated. In this work, we present a doping-induced topological transition effected by plasmon-phonon hybridization in graphene/α-MoO

Recent work has revealed that the light curves of hydrogen-poor (Type I) superluminous supernovae (SLSNe), thought to be powered by magnetar central engines, do not always follow the smooth decline predicted by a simple magnetar spin-down model. Here we present the first systematic study of the prevalence and properties of "bumps" in the post-peak light curves of 34 SLSNe. We find that the majority (44-76%) of events cannot be explained by a smooth magnetar model alone. We do not find any difference in supernova properties between events with and without bumps. By fitting a simple Gaussian model to the light-curve residuals, we characterize each bump with an amplitude, temperature, phase, and duration. We find that most bumps correspond with an increase in the photospheric temperature of the ejecta, although we do not see drastic changes in spectroscopic features during the bump. We also find a moderate correlation (ρ≈0.5; p≈0.01) between the phase of the bumps and the rise time, implying that such bumps tend to happen at a certain "evolutionary phase," (3.7±1.4)trise. Most bumps are consistent with having diffused from a central source of variable luminosity, although sources further out in the ejecta are not excluded. With this evidence, we explore whether the cause of these bumps is intrinsic to the supernova (e.g., a variable central engine) or extrinsic (e.g., circumstellar interaction). Both cases are plausible, requiring low-level variability in the magnetar input luminosity, small decreases in the ejecta opacity, or a thin circumstellar shell or disk.

The search for planets orbiting other stars has recently expanded to include stars from galaxies outside the Milky Way. With the TESS and Gaia surveys, photometric and kinematic information can be combined to identify transiting planet candidates of extragalactic origin. Here, 1,080 low-luminosity red giant branch stars observed by Gaia and TESS with kinematics suggesting a high likelihood of extragalactic origin were searched for planet transits. Transit injection-recovery tests were performed to measure the sensitivity of the TESS data and completeness of the transit search. Injected signals of planets larger than Jupiter with orbital periods of 10 days or less were recovered in ≈44 percent of cases. Although no planet transits were detected in this sample, we find an upper limit on planet occurrence of 0.52 percent for hot Jupiters, consistent with previous studies of planet occurrence around similar host stars. As stars in the halo tend to be lower metallicity, and short period giant planet occurrence tends to be strongly correlated with stellar metallicity, we predict that relative to the Galactic disk population, a smaller fraction of halo stars will host planets detectable by transit surveys. Thus, applying the known planet occurrence trends to potential planet detection around halo stars, we predict ≳7,000 stars must be searched with similar cadence and precision as the stars studied here before a detection of a planet of extragalactic origin is likely. This may be possible with future data releases from the TESS and Gaia missions.

Many active biological particles, such as swimming microorganisms or motor-proteins, do work on their environment by going though a periodic sequence of shapes. Interactions between particles can lead to the phase-synchronization of their duty cycles. Here we consider collective dynamics in a suspension of such active particles coupled through hydrodynamics. We demonstrate that the emergent non-equilibrium states feature stationary patterned flows and robust unidirectional pumping states under confinement. Moreover the phase-synchronized state of the suspension exhibits spatially robust chimera patterns in which synchronized and phase-isotropic regions coexist within the same system. These findings demonstrate a new route to pattern formation and could guide the design of new active materials.