Publications

One-dimensional crystals of passively-driven particles in microfluidic channels exhibit collective vibrational modes reminiscent of acoustic ‘phonons’. These phonons are induced by the long-range hydrodynamic interactions among the particles and are neutrally stable at the linear level. Here, we analyze the effect of particle activity – self-propulsion – on the emergence and stability of these phonons. We show that the direction of wave propagation in active crystals is sensitive to the intensity of the background flow. We also show that activity couples, at the linear level, transverse waves to the particles' rotational motion, inducing a new mode of instability that persists in the limit of large background flow, or, equivalently, vanishingly small activity. We then report a new phenomenon of phonons switching back and forth between two adjacent crystals in both passively-driven and active systems, similar in nature to the wave switching observed in quantum mechanics, optical communication, and density stratified fluids. These findings could have implications for the design of commercial microfluidic systems and the self-assembly of passive and active micro-particles into one-dimensional structures.

Observational tests of stellar and Galactic chemical evolution call for the joint knowledge of a star's physical parameters, detailed element abundances, and precise age. For cool main-sequence (MS) stars the abundances of many elements can be measured from spectroscopy, but ages are very hard to determine. The situation is different if the MS star has a white dwarf (WD) companion and a known distance, as the age of such a binary system can then be determined precisely from the photometric properties of the cooling WD. As a pilot study for obtaining precise age determinations of field MS stars, we identify nearly one hundred candidates for such wide binary systems: a faint WD whose GPS1 proper motion matches that of a brighter MS star in Gaia/TGAS with a good parallax (σϖ/ϖ≤0.05). We model the WD's multi-band photometry with the BASE-9 code using this precise distance (assumed to be common for the pair) and infer ages for each binary system. The resulting age estimates are precise to ≤10% (≤20%) for 42 (67) MS-WD systems. Our analysis more than doubles the number of MS-WD systems with precise distances known to date, and it boosts the number of such systems with precise age determination by an order of magnitude. With the advent of the Gaia DR2 data, this approach will be applicable to a far larger sample, providing ages for many MS stars (that can yield detailed abundances for over 20 elements), especially in the age range 2 to 8\,\Gyr, where there are only few known star clusters.

The generation of high-order harmonics from atomic and molecular gases enables the production of high-energy photons and ultrashort isolated pulses. Obtaining efficiently similar photon energy from solid-state systems could lead, for instance, to more compact extreme ultraviolet and soft x-ray sources. We demonstrate from ab initio si- mulations that it is possible to generate high-order harmonics from free-standing monolayer materials, with an energy cutoff similar to that of atomic and molecular gases. In the limit in which electrons are driven by the pump laser perpendicularly to the monolayer, they behave qualitatively the same as the electrons responsible for high- harmonic generation (HHG) in atoms, where their trajectories are described by the widely used semiclassical model, and exhibit real-space trajectories similar to those of the atomic case. Despite the similarities, the first and last steps of the well-established three-step model for atomic HHG are remarkably different in the two-dimensional materials from gases. Moreover, we show that the electron-electron interaction plays an important role in harmonic gener- ation from monolayer materials because of strong local-field effects, which modify how the material is ionized. The recombination of the accelerated electron wave packet is also found to be modified because of the infinite extension of the material in the monolayer plane, thus leading to a more favorable wavelength scaling of the harmonic yield than in atomic HHG. Our results establish a novel and efficient way of generating high-order harmonics based on a solid-state device, with an energy cutoff and a more favorable wavelength scaling of the harmonic yield similar to those of atomic and molecular gases. Two-dimensional materials offer a unique platform where both bulk and atomic HHG can be investigated, depending on the angle of incidence. Devices based on two-dimensional materials can extend the limit of existing sources.

Estimates of the Hubble constant, H0, from the distance ladder and the cosmic microwave background (CMB) differ at the ∼3-σ level, indicating a potential issue with the standard ΛCDM cosmology. Interpreting this tension correctly requires a model comparison calculation depending on not only the traditional `n-σ' mismatch but also the tails of the likelihoods. Determining the form of the tails of the local H0 likelihood is impossible with the standard Gaussian least-squares approximation, as it requires using non-Gaussian distributions to faithfully represent anchor likelihoods and model outliers in the Cepheid and supernova (SN) populations, and simultaneous fitting of the full distance-ladder dataset to correctly propagate uncertainties. We have developed a Bayesian hierarchical model that describes the full distance ladder, from nearby geometric anchors through Cepheids to Hubble-Flow SNe. This model does not rely on any distributions being Gaussian, allowing outliers to be modeled and obviating the need for arbitrary data cuts. Sampling from the ∼3000-parameter joint posterior using Hamiltonian Monte Carlo, we find H0 = (72.72 ± 1.67) kms−1Mpc−1 when applied to the outlier-cleaned Riess et al. (2016) data, and (73.15±1.78) kms−1Mpc−1 with SN outliers reintroduced. Our high-fidelity sampling of the low-H0 tail of the distance-ladder likelihood allows us to apply Bayesian model comparison to assess the evidence for deviation from ΛCDM. We set up this comparison to yield a lower limit on the odds of the underlying model being ΛCDM given the distance-ladder and Planck XIII (2016) CMB data. The odds against ΛCDM are at worst 10:1 or 7:1, depending on whether the SNe outliers are cut or modeled, or 60:1 if an approximation to the Planck Int. XLVI (2016) likelihood is used.

An important but missing component in the application of the kernel independent fast multipole method (KIFMM) is the capability for flexibly and efficiently imposing singly, doubly, and triply periodic boundary conditions. In most popular packages such periodicities are imposed with the hierarchical repetition of periodic boxes, which may give an incorrect answer due to the conditional convergence of some kernel sums. Here we present an efficient method to properly impose periodic boundary conditions using a near-far splitting scheme. The near-field contribution is directly calculated with the KIFMM method, while the far-field contribution is calculated with a multipole-to-local (M2L) operator which is independent of the source and target point distribution. The M2L operator is constructed with the far-field portion of the kernel function to generate the far-field contribution with the downward equivalent source points in KIFMM. This method guarantees the sum of the near-field & far-field converge pointwise to results satisfying periodicity and compatibility conditions. The computational cost of the far-field calculation observes the same O(N) complexity as FMM and is designed to be small by reusing the data computed by KIFMM for the near-field. The far-field calculations require no additional control parameters, and observes the same theoretical error bound as KIFMM. We present accuracy and timing test results for the Laplace kernel in singly periodic domains and the Stokes velocity kernel in doubly and triply periodic domains.

We study the large source asymptotics of the generating functional in quantum field theory using the holographic renormalization group, and draw comparisons with the asymptotics of the Hopf characteristic function in fractal geometry. Based on the asymptotic behavior, we find a correspondence relating the Weyl anomaly and the fractal dimension of the Euclidean path integral measure. We are led to propose an equivalence between the logarithmic ultraviolet divergence of the Shannon entropy of this measure and the integrated Weyl anomaly, reminiscent of a known relation between logarithmic divergences of entanglement entropy and a central charge. It follows that the information dimension associated with the Euclidean path integral measure satisfies a c-theorem.

Two-dimensional materials equipped with strong spin–orbit coupling can display novel electronic, spintronic, and topological properties originating from the breaking of time or inversion symmetry. A lot of interest has focused on the valley degrees of freedom that can be used to encode binary information. By performing ab initio time-dependent density functional simulation on MoS2, here we show that the spin is not only locked to the valley momenta but strongly coupled to the optical E′′ phonon that lifts the lattice mirror symmetry. Once the phonon is pumped so as to break time-reversal symmetry, the resulting Floquet spectra of the phonon-dressed spins carry a net out-of-plane magnetization (≈0.024μB for single-phonon quantum) even though the original system is non-magnetic. This dichroic magnetic response of the valley states is general for all 2H semiconducting transition-metal dichalcogenides and can be probed and controlled by infrared coherent laser excitation.

The compositions of stars are a critical diagnostic tool for many topics in astronomy such as the evolution of our Galaxy, the formation of planets, and the uniqueness of the Sun. Previous spectroscopic measurements indicate a large intrinsic variation in the elemental abundance patterns of stars with similar overall metal content. However, systematic errors arising from inaccuracies in stellar models are known to be a limiting factor in such studies, and thus it is uncertain to what extent the observed diversity of stellar abundance patterns is real. Here we report the abundances of 30 elements with precisions of 2% for 79 Sun-like stars within 100 parsecs. Systematic errors are minimized in this study by focusing on solar twin stars and performing a line-by-line differential analysis using high-resolution, high-signal-to-noise spectra. We resolve [X/Fe] abundance trends in galactic chemical evolution at precisions of 10^-3 dex Gyr^-1 and reveal that stars with similar ages and metallicities have nearly identical abundance patterns. Contrary to previous results, we find that the ratios of carbon-to-oxygen and magnesium-to-silicon are homogeneous to within 10% throughout the solar neighborhood, implying that exoplanets may exhibit much less compositional diversity than previously thought. Finally, we demonstrate that the Sun has a subtle deficiency in refractory material relative to ~95% of solar twins, suggesting a possible signpost for planetary systems like our own.

Particle picking is a crucial first step in the computational pipeline of single-particle cryo-electron microscopy (cryo-EM). Selecting particles from the micrographs is difficult especially for small particles with low contrast. As high-resolution reconstruction typically requires hundreds of thousands of particles, manually picking that many particles is often too time-consuming. While semi-automated particle picking is currently a popular approach, it may suffer from introducing manual bias into the selection process. In addition, semi-automated particle picking is still somewhat time-consuming. This paper presents the APPLE (Automatic Particle Picking with Low user Effort) picker, a simple and novel approach for fast, accurate, and fully automatic particle picking. While our approach was inspired by template matching, it is completely template-free. This approach is evaluated on publicly available datasets containing micrographs of β-galactosidase and keyhole limpet hemocyanin projections.