We present a proof of concept for a spectrally selective thermal midIR source based on nanopatterned graphene (NPG) with a typical mobility of CVDgrown graphene (up to 3000
There is a growing recognition that subsurface fluid injection can produce not only earthquakes, but also aseismic slip on faults. A major challenge in understanding interactions between injectionrelated aseismic and seismic slip on faults is identifying aseismic slip on the field scale, given that most monitored fields are only equipped with seismic arrays. We present a modeling workflow for evaluating the possibility of aseismic slip, given observational constraints on the spatialtemporal distribution of microseismicity, injection rate, and wellhead pressure. Our numerical model simultaneously simulates discrete offfault microseismic events and aseismic slip on a main fault during fluid injection. We apply the workflow to the 2012 Enhanced Geothermal System injection episode at Cooper Basin, Australia, which aimed to stimulate a watersaturated granitic reservoir containing a highly permeable (
 Award ID(s):
 1947448
 NSFPAR ID:
 10380300
 Publisher / Repository:
 Nature Publishing Group
 Date Published:
 Journal Name:
 Scientific Reports
 Volume:
 12
 Issue:
 1
 ISSN:
 20452322
 Format(s):
 Medium: X
 Sponsoring Org:
 National Science Foundation
More Like this

Abstract ), ensuring scalability to large areas. For that, we solve the electrostatic problem of a conducting hyperboloid with an elliptical wormhole in the presence of an$$\hbox {cm}^2\,\hbox {V}^{1}\,\hbox {s}^{1}$$ ${\text{cm}}^{2}\phantom{\rule{0ex}{0ex}}{\text{V}}^{1}\phantom{\rule{0ex}{0ex}}{\text{s}}^{1}$inplane electric field. The localized surface plasmons (LSPs) on the NPG sheet, partially hybridized with graphene phonons and surface phonons of the neighboring materials, allow for the control and tuning of the thermal emission spectrum in the wavelength regime from to 12$$\lambda =3$$ $\lambda =3$ m by adjusting the size of and distance between the circular holes in a hexagonal or square lattice structure. Most importantly, the LSPs along with an optical cavity increase the emittance of graphene from about 2.3% for pristine graphene to 80% for NPG, thereby outperforming stateoftheart pristine graphene light sources operating in the nearinfrared by at least a factor of 100. According to our COMSOL calculations, a maximum emission power per area of$$\upmu$$ $\mu $ W/$$11\times 10^3$$ $11\times {10}^{3}$ at$$\hbox {m}^2$$ ${\text{m}}^{2}$ K for a bias voltage of$$T=2000$$ $T=2000$ V is achieved by controlling the temperature of the hot electrons through the Joule heating. By generalizing Planck’s theory to any grey body and deriving the completely general nonlocal fluctuationdissipation theorem with nonlocal response of surface plasmons in the random phase approximation, we show that the coherence length of the graphene plasmons and the thermally emitted photons can be as large as 13$$V=23$$ $V=23$ m and 150$$\upmu$$ $\mu $ m, respectively, providing the opportunity to create phased arrays made of nanoantennas represented by the holes in NPG. The spatial phase variation of the coherence allows for beamsteering of the thermal emission in the range between$$\upmu$$ $\mu $ and$$12^\circ$$ ${12}^{\circ}$ by tuning the Fermi energy between$$80^\circ$$ ${80}^{\circ}$ eV and$$E_F=1.0$$ ${E}_{F}=1.0$ eV through the gate voltage. Our analysis of the nonlocal hydrodynamic response leads to the conjecture that the diffusion length and viscosity in graphene are frequencydependent. Using finitedifference time domain calculations, coupled mode theory, and RPA, we develop the model of a midIR light source based on NPG, which will pave the way to graphenebased optical midIR communication, midIR color displays, midIR spectroscopy, and virus detection.$$E_F=0.25$$ ${E}_{F}=0.25$ 
Abstract We present the first unquenched latticeQCD calculation of the form factors for the decay
at nonzero recoil. Our analysis includes 15 MILC ensembles with$$B\rightarrow D^*\ell \nu $$ $B\to {D}^{\ast}\ell \nu $ flavors of asqtad sea quarks, with a strange quark mass close to its physical mass. The lattice spacings range from$$N_f=2+1$$ ${N}_{f}=2+1$ fm down to 0.045 fm, while the ratio between the light and the strangequark masses ranges from 0.05 to 0.4. The valence$$a\approx 0.15$$ $a\approx 0.15$b andc quarks are treated using the Wilsonclover action with the Fermilab interpretation, whereas the light sector employs asqtad staggered fermions. We extrapolate our results to the physical point in the continuum limit using rooted staggered heavylight meson chiral perturbation theory. Then we apply a modelindependent parametrization to extend the form factors to the full kinematic range. With this parametrization we perform a joint latticeQCD/experiment fit using several experimental datasets to determine the CKM matrix element . We obtain$$V_{cb}$$ ${V}_{\mathrm{cb}}$ . The first error is theoretical, the second comes from experiment and the last one includes electromagnetic and electroweak uncertainties, with an overall$$\left V_{cb}\right = (38.40 \pm 0.68_{\text {th}} \pm 0.34_{\text {exp}} \pm 0.18_{\text {EM}})\times 10^{3}$$ $\left({V}_{\mathrm{cb}}\right)=(38.40\pm 0.{68}_{\text{th}}\pm 0.{34}_{\text{exp}}\pm 0.{18}_{\text{EM}})\times {10}^{3}$ , which illustrates the tensions between the experimental data sets, and between theory and experiment. This result is in agreement with previous exclusive determinations, but the tension with the inclusive determination remains. Finally, we integrate the differential decay rate obtained solely from lattice data to predict$$\chi ^2\text {/dof} = 126/84$$ ${\chi}^{2}\text{/dof}=126/84$ , which confirms the current tension between theory and experiment.$$R(D^*) = 0.265 \pm 0.013$$ $R\left({D}^{\ast}\right)=0.265\pm 0.013$ 
Abstract It has been recently established in David and Mayboroda (Approximation of green functions and domains with uniformly rectifiable boundaries of all dimensions.
arXiv:2010.09793 ) that on uniformly rectifiable sets the Green function is almost affine in the weak sense, and moreover, in some scenarios such Green function estimates are equivalent to the uniform rectifiability of a set. The present paper tackles a strong analogue of these results, starting with the “flagship degenerate operators on sets with lower dimensional boundaries. We consider the elliptic operators associated to a domain$$L_{\beta ,\gamma } = {\text {div}}D^{d+1+\gamma n} \nabla $$ ${L}_{\beta ,\gamma}=\text{div}{D}^{d+1+\gamma n}\nabla $ with a uniformly rectifiable boundary$$\Omega \subset {\mathbb {R}}^n$$ $\Omega \subset {R}^{n}$ of dimension$$\Gamma $$ $\Gamma $ , the now usual distance to the boundary$$d < n1$$ $d<n1$ given by$$D = D_\beta $$ $D={D}_{\beta}$ for$$D_\beta (X)^{\beta } = \int _{\Gamma } Xy^{d\beta } d\sigma (y)$$ ${D}_{\beta}{\left(X\right)}^{\beta}={\int}_{\Gamma}{Xy}^{d\beta}d\sigma \left(y\right)$ , where$$X \in \Omega $$ $X\in \Omega $ and$$\beta >0$$ $\beta >0$ . In this paper we show that the Green function$$\gamma \in (1,1)$$ $\gamma \in (1,1)$G for , with pole at infinity, is well approximated by multiples of$$L_{\beta ,\gamma }$$ ${L}_{\beta ,\gamma}$ , in the sense that the function$$D^{1\gamma }$$ ${D}^{1\gamma}$ satisfies a Carleson measure estimate on$$\big  D\nabla \big (\ln \big ( \frac{G}{D^{1\gamma }} \big )\big )\big ^2$$ $D\nabla (ln(\frac{G}{{D}^{1\gamma}})){}^{2}$ . We underline that the strong and the weak results are different in nature and, of course, at the level of the proofs: the latter extensively used compactness arguments, while the present paper relies on some intricate integration by parts and the properties of the “magical distance function from David et al. (Duke Math J, to appear).$$\Omega $$ $\Omega $ 
Abstract The search for neutrino events in correlation with gravitational wave (GW) events for three observing runs (O1, O2 and O3) from 09/2015 to 03/2020 has been performed using the Borexino dataset of the same period. We have searched for signals of neutrinoelectron scattering and inverse betadecay (IBD) within a time window of
s centered at the detection moment of a particular GW event. The search was done with three visible energy thresholds of 0.25, 0.8 and 3.0 MeV. Two types of incoming neutrino spectra were considered: the monoenergetic line and the supernovalike spectrum. GW candidates originated by merging binaries of black holes (BHBH), neutron stars (NSNS) and neutron star and black hole (NSBH) were analyzed separately. Additionally, the subset of most intensive BHBH mergers at closer distances and with larger radiative mass than the rest was considered. In total, followups of 74 out of 93 gravitational waves reported in the GWTC3 catalog were analyzed and no statistically significant excess over the background was observed. As a result, the strongest upper limits on GWassociated neutrino and antineutrino fluences for all flavors ($$\pm \, 1000$$ $\pm \phantom{\rule{0ex}{0ex}}1000$ ) at the level$$\nu _e, \nu _\mu , \nu _\tau $$ ${\nu}_{e},{\nu}_{\mu},{\nu}_{\tau}$ have been obtained in the 0.5–5 MeV neutrino energy range.$$10^9{}10^{15}~\textrm{cm}^{2}\,\textrm{GW}^{1}$$ ${10}^{9}{10}^{15}\phantom{\rule{0ex}{0ex}}{\text{cm}}^{2}\phantom{\rule{0ex}{0ex}}{\text{GW}}^{1}$ 
Abstract Finite volume, weighted essentially nonoscillatory (WENO) schemes require the computation of a smoothness indicator. This can be expensive, especially in multiple space dimensions. We consider the use of the simple smoothness indicator
, where$$\sigma ^{\textrm{S}}= \frac{1}{N_{\textrm{S}}1}\sum _{j} ({\bar{u}}_{j}  {\bar{u}}_{m})^2$$ ${\sigma}^{\text{S}}=\frac{1}{{N}_{\text{S}}1}{\sum}_{j}{({\overline{u}}_{j}{\overline{u}}_{m})}^{2}$ is the number of mesh elements in the stencil,$$N_{\textrm{S}}$$ ${N}_{\text{S}}$ is the local function average over mesh element$${\bar{u}}_j$$ ${\overline{u}}_{j}$j , and indexm gives the target element. Reconstructions utilizing standard WENO weighting fail with this smoothness indicator. We develop a modification of WENOZ weighting that gives a reliable and accurate reconstruction of adaptive order, which we denote as SWENOZAO. We prove that it attains the order of accuracy of the large stencil polynomial approximation when the solution is smooth, and drops to the order of the small stencil polynomial approximations when there is a jump discontinuity in the solution. Numerical examples in one and two space dimensions on general meshes verify the approximation properties of the reconstruction. They also show it to be about 10 times faster in two space dimensions than reconstructions using the classic smoothness indicator. The new reconstruction is applied to define finite volume schemes to approximate the solution of hyperbolic conservation laws. Numerical tests show results of the same quality as standard WENO schemes using the classic smoothness indicator, but with an overall speedup in the computation time of about 3.5–5 times in 2D tests. Moreover, the computational efficiency (CPU time versus error) is noticeably improved.