BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Wenbin Lu

Wenbin Lu contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
astro-ph.HEastro-ph.SRMachine LearningMethodologyastro-ph.GAastro-ph.IMApplicationsDistributed, Parallel, and Cluster Computingmath.STNetworking and Internet Architecturephysics.data-anphysics.flu-dynphysics.plasm-phStatistics Theory

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
30works
0followers
14topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

30 published item(s)

preprint2023arXiv

Mining the Factor Zoo: Estimation of Latent Factor Models with Sufficient Proxies

Latent factor model estimation typically relies on either using domain knowledge to manually pick several observed covariates as factor proxies, or purely conducting multivariate analysis such as principal component analysis. However, the former approach may suffer from the bias while the latter can not incorporate additional information. We propose to bridge these two approaches while allowing the number of factor proxies to diverge, and hence make the latent factor model estimation robust, flexible, and statistically more accurate. As a bonus, the number of factors is also allowed to grow. At the heart of our method is a penalized reduced rank regression to combine information. To further deal with heavy-tailed data, a computationally attractive penalized robust reduced rank regression method is proposed. We establish faster rates of convergence compared with the benchmark. Extensive simulations and real examples are used to illustrate the advantages.

0 citations0 reviewsNo code linkedOpen access
preprint2023arXiv

Science with the Ultraviolet Explorer (UVEX)

UVEX is a proposed medium class Explorer mission designed to provide crucial missing capabilities that will address objectives central to a broad range of modern astrophysics. The UVEX design has two co-aligned wide-field imagers operating in the FUV and NUV and a powerful broadband medium resolution spectrometer. In its two-year baseline mission, UVEX will perform a multi-cadence synoptic all-sky survey 50/100 times deeper than GALEX in the NUV/FUV, cadenced surveys of the Large and Small Magellanic Clouds, rapid target of opportunity followup, as well as spectroscopic followup of samples of stars and galaxies. The science program is built around three pillars. First, UVEX will explore the low-mass, low-metallicity galaxy frontier through imaging and spectroscopic surveys that will probe key aspects of the evolution of galaxies by understanding how star formation and stellar evolution at low metallicities affect the growth and evolution of low-metallicity, low-mass galaxies in the local universe. Such galaxies contain half the mass in the local universe, and are analogs for the first galaxies, but observed at distances that make them accessible to detailed study. Second, UVEX will explore the dynamic universe through time-domain surveys and prompt spectroscopic followup capability will probe the environments, energetics, and emission processes in the early aftermaths of gravitational wave-discovered compact object mergers, discover hot, fast UV transients, and diagnose the early stages of stellar explosions. Finally, UVEX will become a key community resource by leaving a large all-sky legacy data set, enabling a wide range of scientific studies and filling a gap in the new generation of wide-field, sensitive optical and infrared surveys provided by the Rubin, Euclid, and Roman observatories. This paper discusses the scientific potential of UVEX, and the broad scientific program.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Adaptive Semi-Supervised Inference for Optimal Treatment Decisions with Electronic Medical Record Data

A treatment regime is a rule that assigns a treatment to patients based on their covariate information. Recently, estimation of the optimal treatment regime that yields the greatest overall expected clinical outcome of interest has attracted a lot of attention. In this work, we consider estimation of the optimal treatment regime with electronic medical record data under a semi-supervised setting. Here, data consist of two parts: a set of `labeled' patients for whom we have the covariate, treatment and outcome information, and a much larger set of `unlabeled' patients for whom we only have the covariate information. We proposes an imputation-based semi-supervised method, utilizing `unlabeled' individuals to obtain a more efficient estimator of the optimal treatment regime. The asymptotic properties of the proposed estimators and their associated inference procedure are provided. Simulation studies are conducted to assess the empirical performance of the proposed method and to compare with a fully supervised method using only the labeled data. An application to an electronic medical record data set on the treatment of hypotensive episodes during intensive care unit (ICU) stays is also given for further illustration.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Frequency-dependent polarization of repeating fast radio bursts-implications for their origin

The polarization of fast radio bursts (FRBs), bright astronomical transients, contains crucial information about their environments. We report polarization measurements of five repeating FRBs, the abundant signals of which enable wide-band observations with two telescopes. A clear trend of lower polarization at lower frequencies was found, which can be well characterized by a single parameter rotation-measure-scatter (σRM) and modeled by multi-path scatter. Sources with higher σRM have higher RM magnitude and scattering timescales. The two sources with the most substantial σRM, FRB 20121102A and FRB 20190520B, are associated with a compact persistent radio source. These properties indicate a complex environment near the repeating FRBs, such as a supernova remnant or a pulsar wind nebula, consistent with their arising from young populations.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

From Pericenter and Back: Full Debris Stream Evolution in Tidal Disruption Events

When a star passes too close to a supermassive black hole, it gets disrupted by strong tidal forces. The stellar debris then evolves into an elongated stream of gas that partly falls back towards the black hole. We present an analytical model describing for the first time the full stream evolution during such a tidal disruption event (TDE). Our framework consists in dividing the stream into different sections of elliptical geometry, whose properties are independently evolved in their co-moving frame under the tidal, pressure, and self-gravity forces. Through an explicit treatment of the tidal force and the inclusion of the gas angular momentum, we can accurately follow the stream evolution near pericenter. Our model evolves the longitudinal stream stretching and both transverse widths simultaneously. For the latter, we identify two regimes depending on whether the dynamics is entirely dominated by the tidal force (ballistic regime) or additionally influenced by pressure and self-gravity (hydrostatic regime). We find that the stream undergoes transverse collapses both shortly after the stellar disruption and upon its return near the black hole, at specific locations determined by the regime of evolution considered. The stream evolution predicted by our model can be used to determine the subsequent interactions experienced by this gas that are at the origin of most of the electromagnetic emission from TDEs. Our results suggest that the accretion disk may be fed at a rate that differs from the standard fallback rate, which would provide novel observational signatures dependent on black hole spin.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Hydrodynamics of Collisions and Close Encounters between Stellar Black Holes and Main-sequence Stars

Recent analyses have shown that close encounters between stars and stellar black holes occur frequently in dense star clusters. Depending upon the distance at closest approach, these interactions can lead to dissipating encounters such as tidal captures and disruptions, or direct physical collisions, all of which may be accompanied by bright electromagnetic transients. In this study, we perform a wide range of hydrodynamic simulations of close encounters between black holes and main-sequence stars that collectively cover the parameter space of interest, and we identify and classify the various possible outcomes. In the case of nearly head-on collisions, the star is completely disrupted with roughly half of the stellar material becoming bound to the black hole. For more distant encounters near the classical tidal-disruption radius, the star is only partially disrupted on the first pericenter passage. Depending upon the interaction details, the partially disrupted stellar remnant may be tidally captured by the black hole or become unbound (in some cases, receiving a sufficiently large impulsive kick from asymmetric mass loss to be ejected from its host cluster). In the former case, the star will undergo additional pericenter passages before ultimately being disrupted fully. Based on the properties of the material bound to the black hole at the end of our simulations (in particular, the total bound mass and angular momentum), we comment upon the expected accretion process and associated electromagnetic signatures that are likely to result.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

On Learning and Testing of Counterfactual Fairness through Data Preprocessing

Machine learning has become more important in real-life decision-making but people are concerned about the ethical problems it may bring when used improperly. Recent work brings the discussion of machine learning fairness into the causal framework and elaborates on the concept of Counterfactual Fairness. In this paper, we develop the Fair Learning through dAta Preprocessing (FLAP) algorithm to learn counterfactually fair decisions from biased training data and formalize the conditions where different data preprocessing procedures should be used to guarantee counterfactual fairness. We also show that Counterfactual Fairness is equivalent to the conditional independence of the decisions and the sensitive attributes given the processed non-sensitive attributes, which enables us to detect discrimination in the original decision using the processed data. The performance of our algorithm is illustrated using simulated data and real-world applications.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Propagation of Alfvén waves in the charge starvation regime

We present numerical simulation results for the propagation of Alfvén waves in the charge starvation regime. This is the regime where the plasma density is below the critical value required to supply the current for the wave. We analyze a conservative scenario where Alfvén waves pick up charges from the region where the charge density exceeds the critical value and advect them along at a high Lorentz factor. The system consisting of the Alfvén wave and charges being carried with it, which we call charge-carrying Alfvén wave (CC-AW), moves through a medium with small, but non-zero, plasma density. We find that the interaction between CC-AW and the stationary medium has a 2-stream like instability which leads to the emergence of a strong electric field along the direction of the unperturbed magnetic field. The growth rate of this instability is of order the plasma frequency of the medium encountered by the CC-AW. Our numerical code follows the system for hundreds of wave periods. The numerical calculations suggest that the final strength of the electric field is of order a few percent of the Alfvén wave amplitude. Little radiation is produced by the sinusoidally oscillating currents associated with the instability during the linear growth phase. However, in the nonlinear phase, the fluctuating current density produces strong EM radiation near the plasma frequency and limits the growth of the instability.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Spectropolarimetry of the tidal disruption event AT 2019qiz: a quasispherical reprocessing layer

We present optical spectropolarimetry of the tidal disruption event (TDE) AT 2019qiz on days $+0$ and $+29$ relative to maximum brightness. Continuum polarization, which informs the shape of the electron-scattering surface, was found to be consistent with 0 per cent at peak brightness. On day $+29$, the continuum polarization rose to $\sim 1$ per cent, making this the first reported spectropolarimetric evolution of a TDE. These findings are incompatible with a naked eccentric disc that lacks significant mass outflow. Instead, the spectropolarimetry paints a picture wherein, at maximum brightness, high-frequency emission from the accretion disc is reprocessed into the optical band by a nearly spherical, optically thick, electron-scattering photosphere located far away from the black hole. We estimate the radius of the scattering photosphere to be $\sim 100\rm\, au$ at maximum brightness -- significantly larger than the tidal radius ($\sim 1\rm\, au$) and the thermalisation radius ($\sim 30\rm\, au$) where the optical continuum is formed. A month later, as the fallback rate drops and the scattering photosphere recedes, the continuum polarization increases, revealing a moderately aspherical interior. We also see evidence for smaller-scale density variations in the scattering photosphere, inferred from the scatter of the data in the Stokes $q-u$ plane. On day $+29$, the H$α$ emission-line peak is depolarized to $\sim 0.3$ per cent (compared to $\sim 1$ per cent continuum polarization), and displays a gradual rise toward the line's redder wavelengths. This observation indicates the H$α$ line formed near the electron-scattering radius.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Temporal Scattering, Depolarization, and Persistent Radio Emission from Magnetized Inhomogeneous Environments Near Repeating Fast Radio Burst Sources

Some repeating fast radio burst (FRB) sources exhibit complex polarization behaviors, including frequency-dependent depolarization, variation of rotation measure (RM), and oscillating spectral structures of polarized components. Very recently, Feng et al. (2022) reported that active repeaters exhibit conspicuous frequency-dependent depolarization and a strong correlation between RM scatter ($σ_{\rm RM}$) and the temporal scattering time ($τ_{\rm s}$), $σ_{\rm RM}\proptoτ_{\rm s}^{1.0\pm0.2}$, both of which can be well described by multi-path propagation through a magnetized inhomogeneous plasma screen. This observation strongly suggests that the temporal scattering and RM scatter originate from the same region. Besides, a particular finding of note in Feng et al. (2022) is that the FRBs with compact persistent radio sources (PRS) tend to have extreme $σ_{\rm RM}$. In this work, we focus on some theoretical predictions on the relations among temporal scattering, depolarization by RM scatter, and PRS contributed by the magnetized plasma environment close to a repeating FRB source. The behaviors of the RM scatter imply that the magnetized plasma environment is consistent with a supernova remnant or a pulsar wind nebula, and the predicted $σ_{\rm RM}$-$τ_{\rm s}$ relation is $σ_{\rm RM}\proptoτ_{\rm s}^{(0.54-0.83)}$ for different astrophysical scenarios. We further make a general discussion on PRS that does not depend on specific astrophysical scenarios. We show that the specific luminosity of a PRS should have a positive correlation with the RM contributed by the plasma screen. This is consistent with the observations of FRB 121102 and FRB 190520B.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The nozzle shock in tidal disruption events

Tidal disruption events (TDEs) occur when a star gets torn apart by the strong tidal forces of a supermassive black hole, which results in the formation of a debris stream that partly falls back towards the compact object. This gas moves along inclined orbital planes that intersect near pericenter, resulting in a so-called "nozzle shock". We perform the first dedicated study of this interaction, making use of a two-dimensional simulation that follows the transverse gas evolution inside a given section of stream. This numerical approach circumvents the lack of resolution encountered near pericenter passage in global three-dimensional simulations using particle-based methods. As it moves inward, we find that the gas motion is purely ballistic, which near pericenter causes strong vertical compression that squeezes the stream into a thin sheet. Dissipation takes place at the resulting nozzle shock, inducing a rise in pressure that causes the collapsing gas to bounce back, although without imparting significant net expansion. As it recedes to larger distances, this matter continues to expand while remaining thin despite the influence of pressure forces. This gas evolution specifies the strength of the subsequent self-crossing shock, which we find to be more affected by black hole spin than previously estimated. We also evaluate the impact of general-relativistic effects, viscous dissipation, magnetic fields and radiative processes on the nozzle shock. This study represents an important step forward in the theoretical understanding of TDEs, bridging the gap between our robust knowledge of the fallback rate and the more complex following stages, during which most of the emission occurs.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The spins of compact objects born from helium stars in binary systems

The angular momentum (AM) content of massive stellar cores helps to determine the natal spin rates of neutron stars and black holes. Asteroseismic measurements of low-mass stars have proven that stellar cores rotate slower than predicted by most prior work, so revised models are necessary. In this work, we apply an updated AM transport model based on the Tayler instability to massive helium stars in close binaries, in which tidal spin-up can greatly increase the star's AM. Consistent with prior work, these stars can produce highly spinning black holes upon core-collapse if the orbital period is less than $P_{\rm orb} \lesssim \! 1 \, {\rm day}$. For neutron stars, we predict a strong correlation between the pre-explosion mass and the neutron star rotation rate, with millisecond periods ($P_{\rm NS} \lesssim 5 \, {\rm ms}$) only achievable for massive ($M \gtrsim 10 \, M_\odot$) helium stars in tight ($P_{\rm orb} \lesssim 1 \, {\rm day}$) binaries. Finally, we discuss our models in relation to type Ib/c supernovae, superluminous supernove, gamma-ray bursts, and LIGO measurements of black hole spins. Our models are roughly consistent with the rates and energetics of these phenomena, with the exception of broad-lined Ic supernovae, whose high rates and ejecta energies are difficult to explain.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

The former companion of the hyper-velocity star S5-HVS1

The hyper-velocity star S5-HVS1, ejected 5 Myr ago from the Galactic Center at 1800 km/s, was most likely produced by tidal break-up of a tight binary by the supermassive black hole SgrA*. Taking a Monte Carlo approach, we show that the former companion of S5-HVS1 was likely a main-sequence star between 1.2 and 6 solar masses and was captured into a highly eccentric orbit with pericenter distance in the range 1-10 AU and semimajor axis about $10^3$ AU. We then explore the fate of the captured star. We find that the heat deposited by tidally excited stellar oscillation modes leads to runaway disruption if the pericenter distance is smaller than about 3 AU. Over the past 5 Myr, its angular momentum has been significantly modified by orbital relaxation, which may stochastically drive the pericenter inwards below 3 AU and cause tidal disruption. We find an overall survival probability in the range 5% to 50%, depending on the local relaxation time in the close environment of the captured star, and the initial pericenter at capture. The pericenter distance of the surviving star has migrated to 10-100 AU, making it potentially the most extreme member of the S-star cluster. From the ejection rate of S5-HVS1-like stars, we estimate that there may currently be a few stars in such highly eccentric orbits. They should be detectable (typically Ks < 18.5 mag) by the GRAVITY instrument and by future Extremely Large Telescopes and hence provide an extraordinary probe of the spin of SgrA*.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

The nascent milliquasar VT J154843.06+220812.6: tidal disruption event or extreme accretion-state change?

We present detailed multiwavelength follow up of a nuclear radio flare, VT J154843.06+220812.6, hereafter VT J1548. VT J1548 was selected as a ${\sim}1$ mJy radio flare in 3 GHz observations from the VLA Sky Survey (VLASS). It is located in the nucleus of a low mass ($\log M_{\rm BH}/M_\odot \sim6$) host galaxy with weak or no past AGN activity. VT J1548 is associated with a slow rising (multiple year), bright mid IR flare in the WISE survey, peaking at ${\sim}10\%L_{\rm edd.}$. No associated optical transient is detected, although we cannot rule out a short, early optical flare given the limited data available. Constant late time (${\sim}3$ years post-flare) X-ray emission is detected at ${\sim}10^{42}$ erg s$^{-1}$. The radio SED is consistent with synchrotron emission from an outflow incident on an asymmetric medium. A follow-up, optical spectrum shows transient, bright, high-ionization coronal line emission ($[{\rm Fe\,X}]\,λ6375,[{\rm Fe\,XI}]\,λ7894,[{\rm S\,XII}]\,λ7612$). Transient broad H$α$ is also detected but without corresponding broad H$β$ emission, suggesting high nuclear extinction. We interpret this event as either a tidal disruption event or an extreme flare of an active galactic nucleus, in both cases obscured by a dusty torus. Although these individual properties have been observed in previous transients, the combination is unprecedented. This event highlights the importance of searches across all wave bands for assembling a sample of nuclear flares that spans the range of observable properties and possible triggers.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

UCX Programming Interface for Remote Function Injection and Invocation

Network library APIs have historically been developed with the emphasis on data movement, placement, and communication semantics. Many communication semantics are available across a large variety of network libraries, such as send-receive, data streaming, put/get/atomic, RPC, active messages, collective communication, etc. In this work we introduce new compute and data movement APIs that overcome the constraints of the single-program, multiple-data (SPMD) programming model by allowing users to send binary executable code between processing elements. Our proof-of-concept implementation of the API is based on the UCX communication framework and leverages the RDMA network for fast compute migration. We envision the API being used to dispatch user functions from a host CPU to a SmartNIC (DPU), computational storage drive (CSD), or remote servers. In addition, the API can be used by large-scale irregular applications (such as semantic graph analysis), composed of many coordinating tasks operating on a data set so big that it has to be stored on many physical devices. In such cases, it may be more efficient to dynamically choose where code runs as the applications progresses.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A New Framework for Online Testing of Heterogeneous Treatment Effect

We propose a new framework for online testing of heterogeneous treatment effects. The proposed test, named sequential score test (SST), is able to control type I error under continuous monitoring and detect multi-dimensional heterogeneous treatment effects. We provide an online p-value calculation for SST, making it convenient for continuous monitoring, and extend our tests to online multiple testing settings by controlling the false discovery rate. We examine the empirical performance of the proposed tests and compare them with a state-of-art online test, named mSPRT using simulations and a real data. The results show that our proposed test controls type I error at any time, has higher detection power and allows quick inference on online A/B testing.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A unified picture of Galactic and cosmological fast radio bursts

The discovery of a fast radio burst (FRB) in our galaxy associated with a magnetar (neutron star with strong magnetic field) has provided a critical piece of information to help us finally understand these enigmatic transients. We show that the volumetric rate of Galactic-FRB like events is consistent with the faint end of the cosmological FRB rate, and hence they most likely belong to the same class of transients. The Galactic FRB had an accompanying X-ray burst but many X-ray bursts from the same object had no radio counterpart. Their relative rates suggest that for every FRB there are roughly 10^2 to 10^3 X-ray bursts. The radio lightcurve of the galactic FRB had two spikes separated by 30 ms in the 400-800 MHz frequency band. This is an important clue and highly constraining of the class of models where the radio emission is produced outside the light-cylinder of the magnetar. We suggest that magnetic disturbances close to the magnetar surface propagate to a distance of a few tens of neutron star radii where they damp and produce radio emission. The coincident hard X-ray spikes associated with the two FRB pulses seen in this burst and the flux ratio between the two frequency bands can be understood in this scenario. This model provides a unified picture for faint bursts like the Galactic FRB as well as the bright events seen at cosmological distances.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Causal Effect Estimation and Optimal Dose Suggestions in Mobile Health

In this article, we propose novel structural nested models to estimate causal effects of continuous treatments based on mobile health data. To find the treatment regime that optimizes the expected short-term outcomes for patients, we define a weighted lag-K advantage as the value function. The optimal treatment regime is then defined to be the one that maximizes the value function. Our method imposes minimal assumptions on the data generating process. Statistical inference is provided for the estimated parameters. Simulation studies and an application to the Ohio type 1 diabetes dataset show that our method could provide meaningful insights for dose suggestions with mobile health data.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Deceleration of relativistic jets with lateral expansion

We present a model for the hydrodynamics of a relativistic jet interacting with the circum-stellar medium (CSM). The shocked CSM and the jet material are assumed to be in an infinitely thin surface, so the original 2D problem is effectively reduced to 1D. From general conservation laws, we derive the equation of motion for each fluid element along this surface, taking into account the deceleration along the surface normal due to newly swept-up mass and lateral expansion due to pressure gradient in the tangential direction. The pressure and energy density of the shocked CSM are given by the jump conditions at the forward shock. The method is implemented with a finite-differencing numerical scheme, along with calculation of synchrotron emission and absorption from shock-accelerated electrons, in a new code $\mathtt{Jedi}$ (for &#34;jet dynamics&#34;). We present a number of test cases, including top-hat jet, power-law structured jet, &#34;boosted fireball&#34; profile, and CSM with density jump at the wind termination shock. Based on the agreement with other analytical and numerical calculations, we conclude that our simplified method provides a good approximation for the hydrodynamics and afterglow emission for a wide variety of jet structures and CSM density profiles. Efficient modeling of the afterglow from e.g., neutron star mergers, will provide important information on the jet energetics, CSM properties, and the viewing angle.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Does the Markov Decision Process Fit the Data: Testing for the Markov Property in Sequential Decision Making

The Markov assumption (MA) is fundamental to the empirical validity of reinforcement learning. In this paper, we propose a novel Forward-Backward Learning procedure to test MA in sequential decision making. The proposed test does not assume any parametric form on the joint distribution of the observed data and plays an important role for identifying the optimal policy in high-order Markov decision processes and partially observable MDPs. We apply our test to both synthetic datasets and a real data example from mobile health studies to illustrate its usefulness.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Implications of CHIME repeating fast radio bursts

CHIME has now detected 18 repeating fast radio bursts (FRBs). We explore what can be learned about the energy distribution and activity level of the repeaters by fitting realistic FRB population models to the data. For a power-law energy distribution dN/dE ~ E^{-γ} for the repeating bursts, there is a critical index γ_c that controls whether the dispersion measure (DM, a proxy for source distance) distribution of repeaters is bottom or top-heavy. We find γ_c = 7/4 for Poisson wait-time distribution of repeaters in Euclidean space and further demonstrate how it is affected by temporal clustering of repetitions and cosmological effects. It is especially interesting that two of the CHIME repeaters (FRB 181017 and 190417) have large DM ~ 1000 pc/cm^3. These can be understood if: (i) the energy distribution is shallow γ= 1.7^{+0.3}_{-0.1} (68% confidence) or (ii) a small fraction of sources are extremely active. In the second scenario, these two high-DM sources should be repeating more than 100 times more frequently than FRB 121102, and the energy index is constrained to be γ= 1.9^{+0.3}_{-0.2} (68% confidence). In either case, this γis consistent with the energy dependence of the non-repeating ASKAP sample, which suggests that they are drawn from the same population. Finally, our model predicts how the CHIME repeating fraction should decrease with redshift, and this can be compared with observations to infer the distribution of activity level in the whole population.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Imprint of local environment on fast radio burst observations

When fast radio burst (FRB) waves propagate through the local (< 1 pc) environment of the FRB source, electrons in the plasma undergo large-amplitude oscillations. The finite-amplitude effects cause the effective plasma frequency and cyclotron frequency to be dependent on the wave strength. The dispersion measure and rotation measure should therefore vary slightly from burst to burst for a repeating source, depending on the luminosity and frequency of the individual burst. Furthermore, free-free absorption of strong waves is suppressed due to the accelerated electrons&#39; reduced energy exchange in Coulomb collisions. This allows bright low-frequency bursts to propagate through an environment that would be optically thick to low-amplitude waves. Given a large sample of bursts from a repeating source, it would be possible to use the deficit of low-frequency and low-luminosity bursts to infer the emission measure of the local intervening plasma and its distance from the source. Information about the local environment will shed light on the nature of FRB sources.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Kernel Assisted Learning for Personalized Dose Finding

An individualized dose rule recommends a dose level within a continuous safe dose range based on patient level information such as physical conditions, genetic factors and medication histories. Traditionally, personalized dose finding process requires repeating clinical visits of the patient and frequent adjustments of the dosage. Thus the patient is constantly exposed to the risk of underdosing and overdosing during the process. Statistical methods for finding an optimal individualized dose rule can lower the costs and risks for patients. In this article, we propose a kernel assisted learning method for estimating the optimal individualized dose rule. The proposed methodology can also be applied to all other continuous decision-making problems. Advantages of the proposed method include robustness to model misspecification and capability of providing statistical inference for the estimated parameters. In the simulation studies, we show that this method is capable of identifying the optimal individualized dose rule and produces favorable expected outcomes in the population. Finally, we illustrate our approach using data from a warfarin dosing study for thrombosis patients.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Radiation Forces Constrain the FRB Mechanism

We provide constraints on Fast Radio Burst (FRB) models by careful considerations of radiation forces associated with these powerful transients. We find that the induced-Compton scatterings of the coherent radiation by electrons/positrons accelerate particles to very large Lorentz factors (LF) in and around the source of this radiation. This severely restricts those models for FRBs that invoke relativistic shocks and maser type instabilities at distances less than about $10^{13}$ cm of the neutron star. Radiation traveling upstream, in these models, forces particles to move away from the shock with a LF larger than the LF of the shock front. This suspends the photon generation process after it has been operating for less than ~0.1 ms (observer frame duration). We show that masers operating in shocks at distances larger than $10^{13}$ cm cannot simultaneously account for the burst duration of 1 ms or more and the observed ~1 GHz frequencies of FRBs without requiring an excessive energy budget ($10^{46}$ erg); the energy is not calculated by imposing any efficiency consideration, or other details, for the maser mechanism, but is entirely the result of ensuring that particle acceleration by induced-Compton forces upstream of the shock front does not choke off the maser process. For the source to operate more or less continuously for a few ms, it should be embedded in a strong magnetic field - cyclotron frequency $\gg$ wave frequency - so that radiation forces do not disperse the plasma and shut-off the engine.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The multiwavelength counterparts of fast radio bursts

The engines that produce extragalactic fast radio bursts (FRBs), and the mechanism by which the emission is generated, remain unknown. Many FRB models predict prompt multi-wavelength counterparts, which can be used to refine our knowledge of these fundamentals of the FRB phenomenon. However, several previous targeted searches for prompt FRB counterparts have yielded no detections, and have additionally not reached sufficient sensitivity with respect to the predictions. In this work, we demonstrate a technique to estimate the ratio, $η$, between the energy outputs of FRB counterparts at various wavelengths and the radio-wavelength emission. Our technique combines the fluence distribution of the FRB population with results from several wide-field blind surveys for fast transients from the optical to the TeV bands. We present constraints on $η$ that improve upon previous observations even in the case that all unclassified transient events in existing surveys are FRB counterparts. In some scenarios for the FRB engine and emission mechanism, we find that FRB counterparts should have already been detected, thus demonstrating that our technique can successfully test predictions for $η$. However, it is possible that FRB counterparts are lurking amongst catalogs of unclassified transient events. Although our technique is robust to the present uncertainty in the FRB fluence distribution, its ultimate application to accurately estimate or bound $η$ will require the careful analysis of all candidate fast-transient events in multi-wavelength survey data sets.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Wind-reprocessed Transients

We consider the situation where the luminosity from a transient event is reprocessed by an optically thick wind. Potential applications are the tidal disruption of stars by black holes, engine-powered supernovae, and unique fast transients found by current and future wide-field surveys. We derive relations between the injected and observed luminosity for steady and time dependent winds, and discuss how the temperature is set for scattering-dominated radiative transport. We apply this framework to specific examples of tidal disruption events and the formation of a black hole by a massive star, as well as discuss other applications such as deriving observables from detailed hydrodynamic simulations. We conclude by exploring what is inferred about the mass loss rate and underlying engine powering AT2018cow if it is explained as a wind-reprocessed transient, demonstrating that its optical emission is consistent with reprocessing of the observed soft X-rays.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Zwicky Transient Facility constraints on the optical emission from the nearby repeating FRB 180916.J0158+65

The discovery rate of fast radio bursts (FRBs) is increasing dramatically thanks to new radio facilities. Meanwhile, wide-field instruments such as the 47 deg$^2$ Zwicky Transient Facility (ZTF) survey the optical sky to study transient and variable sources. We present serendipitous ZTF observations of the CHIME repeating source FRB 180916.J0158+65, that was localized to a spiral galaxy 149 Mpc away and is the first FRB suggesting periodic modulation in its activity. While 147 ZTF exposures corresponded to expected high-activity periods of this FRB, no single ZTF exposure was at the same time as a CHIME detection. No $>3σ$ optical source was found at the FRB location in 683 ZTF exposures, totalling 5.69 hours of integration time. We combined ZTF upper limits and expected repetitions from FRB 180916.J0158+65 in a statistical framework using a Weibull distribution, agnostic of periodic modulation priors. The analysis yielded a constraint on the ratio between the optical and radio fluences of $η\lesssim 200$, corresponding to an optical energy $E_{\rm opt} \lesssim 3 \times 10^{46}$ erg for a fiducial 10 Jy ms FRB (90% confidence). A deeper (but less statistically robust) constraint of $η\lesssim 3$ can be placed assuming a rate of $r(>5$ Jy ms)= hr$^{-1}$ and $1.2\pm 1.1$ FRB occurring during exposures taken in high-activity windows. The constraint can be improved with shorter per-image exposures and longer integration time, or observing FRBs at higher Galactic latitudes. This work demonstrated how current surveys can statistically constrain multi-wavelength counterparts to FRBs even without deliberately scheduled simultaneous radio observation.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Jets from Tidal Disruption Events

The discovery of jets from tidal disruption events (TDEs) rejuvenated the old field of relativistic jets powered by accretion onto supermassive black holes. In this Chapter, we first review the extensive multi-wavelength observations of jetted TDEs. Then, we show that these events provide valuable information on many aspects of jet physics from a new prospective, including the on-and-off switch of jet launching, jet propagation through the ambient medium, $γ/$X-ray radiation mechanism, jet composition, and the multi-messenger picture. Finally, open questions and future prospects in this field are summarized.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Self-intersection of the Fallback Stream in Tidal Disruption Events

We propose a semi-analytical model for the self-intersection of the fallback stream in tidal disruption events (TDEs). When the initial periapsis is less than about 15 gravitational radii, a large fraction of the shocked gas is unbound in the form of a collision-induced outflow (CIO). This is because large apsidal precession causes the stream to self-intersect near the local escape speed at radius much below the apocenter. The rest of the fallback gas is left in more tightly bound orbits and quickly joins the accretion flow. We propose that the CIO is responsible for reprocessing the hard emission from the accretion flow into the optical band. This picture naturally explains the large photospheric radius (or low blackbody temperature) and typical widths of the H and/or He emission lines seen in optical TDEs. We predict the CIO-reprocessed spectrum in the infrared to be L_ν~ ν^{~0.5}, shallower than a blackbody. The partial sky coverage of the CIO also provides a unification of the diverse X-ray behaviors of optical TDEs. According to this picture, optical surveys filter out a large fraction of TDEs with low-mass blackholes due to lack of a reprocessing layer, and the volumetric rate of optical TDEs is nearly flat wrt. the blackhole mass for M < 10^7 solar masses. This filtering causes the optical TDE rate to be lower than the total rate by a factor of ~10 or more. When the CIO is decelerated by the ambient medium, radio emission at the level of that in ASASSN-14li may be produced, but the timescales and peak luminosities can be highly diverse. Finally, our method paves the way for global simulations of the disk formation process by injecting gas at the intersection point according to the prescribed velocity and density profiles.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Simulating realistic disc formation in tidal disruption events

A star coming too close to a supermassive black hole gets disrupted by the tidal force of the compact object in a tidal disruption event, or TDE. Following this encounter, the debris evolves into an elongated stream, half of which coming back to pericenter. Relativistic apsidal precession then leads to a self-crossing shock that initiates the formation of an accretion disc. We perform the first simulation of this process considering a realistic stellar trajectory and black hole mass, which has so far eluded investigations for computational reasons. This numerical issue is alleviated by using as initial conditions the outflow launched by the self-crossing shock according the local simulation of Lu & Bonnerot (2019). We find that the gas leaving the intersection point experiences numerous secondary shocks that result in the rapid formation of a thick and marginally-bound disc. The mass distribution features two overdensities identified as spiral shocks that drive slow gas inflow along the mid-plane. Inward motion primarily takes place along the funnels of the newly-formed torus, from which a fraction of the matter can get accreted. Further out, the gas moves outward forming an extended envelope completely surrounding the accretion flow. Secondary shocks heat the debris at a rate of a few times $10^{44} \, \rm erg \, s^{-1}$ with a large fraction likely participating to the bolometric luminosity. These results pave the way towards a complete understanding of the early radiation from TDEs that progressively becomes accessible from observations.

0 citations0 reviewsNo code linkedOpen access