• We present new observations of the binary neutron star merger GW170817 at $\Delta t\approx 220-260$ days post-merger, at radio (Karl G.~Jansky Very Large Array; VLA), X-ray ({\it Chandra X-ray Observatory}) and optical ({\it Hubble Space Telescope}; HST) wavelengths. These observations provide the first evidence for a turnover in the X-ray light curve, mirroring a decline in the radio emission at $\gtrsim5\sigma$ significance. The radio to X-ray spectral energy distribution exhibits no evolution into the declining phase. Our full multi-wavelength dataset is consistent with the predicted behavior of our previously-published models of a structured jet expanding into a low density circumbinary medium. If future observations continue to track our predictions, we expect the radio and X-ray emission to remain detectable until $\sim 1000$ days post-merger.
  • We report deep Chandra, HST and VLA observations of the binary neutron star event GW170817 at $t<160$ d after merger. These observations show that GW170817 has been steadily brightening with time and might have now reached its peak, and constrain the emission process as non-thermal synchrotron emission where the cooling frequency $\nu_c$ is above the X-ray band and the synchrotron frequency $\nu_m$ is below the radio band. The very simple power-law spectrum extending for eight orders of magnitude in frequency enables the most precise measurement of the index $p$ of the distribution of non-thermal relativistic electrons $N(\gamma)\propto \gamma^{-p}$ accelerated by a shock launched by a NS-NS merger to date. We find $p=2.17\pm0.01$, which indicates that radiation from ejecta with $\Gamma\sim3-10$ dominates the observed emission. While constraining the nature of the emission process, these observations do \emph{not} constrain the nature of the relativistic ejecta. We employ simulations of explosive outflows launched in NS ejecta clouds to show that the spectral and temporal evolution of the non-thermal emission from GW170817 is consistent with both emission from radially stratified quasi-spherical ejecta traveling at mildly relativistic speeds, \emph{and} emission from off-axis collimated ejecta characterized by a narrow cone of ultra-relativistic material with slower wings extending to larger angles. In the latter scenario, GW170817 harbored a normal SGRB directed away from our line of sight. Observations at $t\le 200$ days are unlikely to settle the debate as in both scenarios the observed emission is effectively dominated by radiation from mildly relativistic material.
  • We propose a data-driven filtered reduced order model (DDF-ROM) framework for the numerical simulation of fluid flows. The novel DDF-ROM framework consists of two steps: (i) In the first step, we use explicit ROM spatial filtering of the nonlinear PDE to construct a filtered ROM. This filtered ROM is low-dimensional, but is not closed (because of the nonlinearity in the given PDE). (ii) In the second step, we use data-driven modeling to close the filtered ROM, i.e., to model the interaction between the resolved and unresolved modes. To this end, we use a quadratic ansatz to model this interaction and close the filtered ROM. To find the new coefficients in the closed filtered ROM, we solve an optimization problem that minimizes the difference between the full order model data and our ansatz. We emphasize that the new DDF-ROM is built on general ideas of spatial filtering and optimization and is independent of (restrictive) phenomenological arguments. We investigate the DDF-ROM in the numerical simulation of a 2D channel flow past a circular cylinder at Reynolds number $Re=100$. The DDF-ROM is significantly more accurate than the standard projection ROM. Furthermore, the computational costs of the DDF-ROM and the standard projection ROM are similar, both costs being orders of magnitude lower than the computational cost of the full order model. We also compare the new DDF-ROM with modern ROM closure models in the numerical simulation of the 1D Burgers equation. The DDF-ROM is more accurate and significantly more efficient than these ROM closure models.
  • Spatial filtering has been central in the development of large eddy simulation reduced order models (LES-ROMs) and regularized reduced order models (Reg-ROMs), In this paper, we perform a numerical investigation of spatial filtering. To this end, we consider one of the simplest Reg-ROMs, the Leray ROM (L-ROM), which uses ROM spatial filtering to smooth the flow variables and decreases the amount of energy aliased to the lower index ROM basis functions. We also propose a new form of ROM differential filter and use it as a spatial filter for the L-ROM. We investigate the performance of this new form of ROM differential filter in the numerical simulation of a flow past a circular cylinder at a Reynolds number $Re=760$.
  • We propose a calibrated filtered reduced order model (CF-ROM) framework for the numerical simulation of general nonlinear PDEs that are amenable to reduced order modeling. The novel CF-ROM framework consists of two steps: (i) In the first step, we use explicit ROM spatial filtering of the nonlinear PDE to construct a filtered ROM. This filtered ROM is low-dimensional, but is not closed (because of the nonlinearity in the given PDE). (ii) In the second step, we use a calibration procedure to close the filtered ROM, i.e., to model the interaction between the resolved and unresolved modes. To this end, we use a linear or quadratic ansatz to model this interaction and close the filtered ROM. To find the new coefficients in the closed filtered ROM, we solve an optimization problem that minimizes the difference between the full order model data and our ansatz. Although we use a fluid dynamics setting to illustrate how to construct and use the CF-ROM framework, we emphasize that it is built on general ideas of spatial filtering and optimization and is independent of (restrictive) phenomenological arguments. Thus, the CF-ROM framework can be applied to a wide variety of PDEs.
  • We demonstrated a laser depth imaging system based on the time-correlated single-photon counting technique, which was incorporated with a low-jitter superconducting nanowire single-photon detector (SNSPD), operated at the wavelength of 1550 nm. A sub-picosecond time-bin width was chosen for photon counting, resulting in a discrete noise of less than one/two counts for each time bin under indoor/outdoor daylight conditions, with a collection time of 50 ms. Because of the low-jitter SNSPD, the target signal histogram was significantly distinguishable, even for a fairly low retro-reflected photon flux. The depth information was determined directly by the highest bin counts, instead of using any data fitting combined with complex algorithms. Millimeter resolution depth imaging of a low-signature object was obtained, and more accurate data than that produced by the traditional Gaussian fitting method was generated. Combined with the intensity of the return photons, three-dimensional reconstruction overlaid with reflectivity data was realized.
  • We study the expected number of real zeros for random linear combinations of orthogonal polynomials. It is well known that Kac polynomials, spanned by monomials with i.i.d. Gaussian coefficients, have only $(2/\pi + o(1))\log{n}$ expected real zeros in terms of the degree $n$. On the other hand, if the basis is given by Legendre (or more generally by Jacobi) polynomials, then random linear combinations have $n/\sqrt{3} + o(n)$ expected real zeros. We prove that the latter asymptotic relation holds universally for a large class of random orthogonal polynomials on the real line, and also give more general local results on the expected number of real zeros.
  • The Stokes-Einstein relation for the self-diffusion coefficient of a spherical particle suspended in an incompressible fluid is an asymptotic result in the limit of large Schmidt number, that is, when momentum diffuses much faster than the particle. When the Schmidt number is moderate, which happens in most particle methods for hydrodynamics, deviations from the Stokes-Einstein prediction are expected. We study these corrections computationally using a recently-developed minimally-resolved method for coupling particles to an incompressible fluctuating fluid in both two and three dimensions. We find that for moderate Schmidt numbers the diffusion coefficient is reduced relative to the Stokes-Einstein prediction by an amount inversely proportional to the Schmidt number in both two and three dimensions. We find, however, that the Einstein formula is obeyed at all Schmidt numbers, consistent with linear response theory. The numerical data is in good agreement with an approximate self-consistent theory, which can be used to estimate finite-Schmidt number corrections in a variety of methods. Our results indicate that the corrections to the Stokes-Einstein formula come primarily from the fact that the particle itself diffuses together with the momentum. Our study separates effects coming from corrections to no-slip hydrodynamics from those of finite separation of time scales, allowing for a better understanding of widely observed deviations from the Stokes-Einstein prediction in particle methods such as molecular dynamics.
  • Single-layer FeSe films with extremely expanded in-plane lattice constant of 3.99A are fabricated by epitaxially growing FeSe/Nb:SrTiO3/KTaO3 heterostructures, and studied by in situ angle-resolved photoemission spectroscopy. Two elliptical electron pockets at the Brillion zone corner are resolved with negligible hybridization between them, indicating the symmetry of the low energy electronic structure remains intact as a free-standing single-layer FeSe, although it is on a substrate. The superconducting gap closes at a record high temperature of 70K for the iron based superconductors. Intriguingly, the superconducting gap distribution is anisotropic but nodeless around the electron pockets, with minima at the crossings of the two pockets. Our results put strong constraints on the current theories, and support the coexistence of both even and odd parity spin-singlet pairing channels as classified by the lattice symmetry.
  • Transition metal oxide hetero-structure has great potential for multifunctional devices. However, the degraded physical properties at interface, known as dead-layer behavior, present a main obstacle for device applications. Here we present the systematic study of the dead-layer behavior in La0.67Sr0.33MnO3 thin film grown on SrTiO3 substrate with ozone assisted molecular beam epitaxy. We found that the evolution of electric and magnetic properties as a function of thickness shows a remarkable resemblance to the phase diagram as a function of doping for bulk materials, providing compelling evidences of the hole depletion in near interface layers that causes dead-layer. Detailed electronic and surface structure studies indicate that the hole depletion is due to the intrinsic oxygen vacancy formation. Furthermore, we show that oxygen vacancies are partly caused by interfacial electric dipolar field, and thus by doping-engineering at the single-atomic-layer level, we demonstrate the dead-layer reduction in films with higher interfacial hole concentration.
  • The recent development of techniques for laser-driven shock compression of hydrogen has opened the door to the experimental determination of its behavior under conditions characteristic of stellar and planetary interiors. The new data probe the equation of state (EOS) of dense hydrogen in the complex regime of pressure ionization. The structure and evolution of dense astrophysical bodies depend on whether the pressure ionization of hydrogen occurs continuously or through a ``plasma phase transition'' (PPT) between a molecular state and a plasma state. For the first time, the new experiments constrain predictions for the PPT. We show here that the EOS model developed by Saumon and Chabrier can successfully account for the data, and we propose an experiment that should provide a definitive test of the predicted PPT of hydrogen. The usefulness of the chemical picture for computing astrophysical EOS and in modeling pressure ionization is discussed.
  • Raman scattering from the acoustic phonons confied in si nanocrystals is investigated by a microscopic lattice-dynamical calculation. Polarized and depolarized spectra are found to be rather different, indicated from different acoustic phonons. The polarized scattering is more pronounced than the depolarized one. The calculated results are compared with the recent experimental data. The matrix effects are discussed.