Source author record

Liang Dai

Liang Dai appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

43works
19topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

43 published item(s)

preprint2024arXiv

Enhancement of cell membrane poration by the antimicrobial peptide Melp5

Melittin, a natural antimicrobial peptide comprising 26 amino acid residues, can kill bacteria by inducing pores in cell membranes. Clinical applications of melittin as an antibiotic require a thorough understanding of its poration mechanism and mutations that enhance its antimicrobial activity. Previous experiments showed Melp5, a variant of melittin with five mutations, exhibits a higher poration ability. However, the mechanism of the enhanced poration ability is not fully understood. Here, we investigated the mechanism by comparing the poration of melittin and Melp5 using coarse-grained (CG) and all-atom (AA) molecular dynamics (MD) simulations. We observe that Melp5 is likely to form a pore with 5 peptides (pentameric), while melittin is likely to form a pore with 4 peptides (tetrameric). Our atomistic MD simulations show that the pentameric pore of Melp5 has a higher water permeability than the tetrameric pore of melittin. We also analyze the stability of the pores of melittin and Melp5 by calculating the interaction energies of the pores. In particular, we investigate the effects of mutant residues on pore stability by calculating electrostatic and LJ interactions. These results should provide insights on the enhanced poration ability of Melp5 and push it toward applications.

preprint2023arXiv

Mapping the Likelihood of GW190521 with Diverse Mass and Spin Priors

We map the likelihood of GW190521, the heaviest detected binary black hole (BBH) merger, by sampling under different mass and spin priors designed to be uninformative. We find that a source-frame total mass of $\sim$$150 M_{\odot}$ is consistently supported, but posteriors in mass ratio and spin depend critically on the choice of priors. We confirm that the likelihood has a multi-modal structure with peaks in regions of mass ratio representing very different astrophysical scenarios. The unequal-mass region ($m_2 / m_1 < 0.3$) has an average likelihood $\sim$$e^6$ times larger than the equal-mass region and a maximum likelihood $\sim$$e^2$ larger. Using ensembles of samples across priors, we examine the implications of qualitatively different BBH sources that fit the data. We find that the equal-mass solution has poorly constrained spins and at least one black hole mass that is difficult to form via stellar collapse due to (pulsational) pair instability. The unequal-mass solution can avoid this mass gap entirely but requires a negative effective spin and a precessing primary. Both of these scenarios are more easily produced by dynamical formation channels than field binary co-evolution. Drawing representative samples from each region of the likelihood map, we find a sensitive comoving volume-time $\mathcal{O}(10)$ times larger in the mass gap region than the gap-avoiding region. Accounting for this distance effect, the likelihood still reverses the advantage to favor the gap-avoiding scenario by a factor of $\mathcal{O}(100)$ before including mass and spin priors. Posterior samplers can be driven away from this high-likelihood region by common prior choices meant to be uninformative, making GW190521 parameter inference sensitive to the choice of mass and spin priors. This may be a generic issue for similarly heavy events given current detector sensitivity and waveform degeneracies.

preprint2022arXiv

Possible Ongoing Merger Discovered by Photometry and Spectroscopy in the Field of the Galaxy Cluster PLCK G165.7+67.0

We present a detailed study of the Planck-selected binary galaxy cluster PLCK G165.7+67.0 (G165; $z$=0.348). A multiband photometric catalog is generated that incorporates new imaging from the Large Binocular Telescope/Large Binocular Camera and Spitzer/IRAC to existing imaging. To cope with the different image characteristics, robust methods are applied in the extraction of the matched-aperture photometry. Photometric redshifts are estimated for 143 galaxies in the 4 arcmin$^{2}$ field of overlap covered by all these data. We confirm that strong lensing effects yield 30 images of 11 background galaxies, of which we contribute photometric redshift estimates for three image multiplicities. These constraints enable the construction of a revised lens model that confirms the bimodal structure, and from which we measure a mass of M$_{600 kpc}$=(2.36$\pm$0.23)$\times$10$^{14}$M$_{\odot}$. In parallel, new spectroscopy using MMT/Binospec and archival data contributes thirteen galaxies which meet our velocity and transverse radius criteria for cluster membership. The two cluster components have a pair-wise velocity of $\lessapprox$100 kms$^{-1}$, favoring an orientation in the plane of the sky with a transverse velocity of 100-1700 kms$^{-1}$. At the same time, the brightest cluster galaxy is offset in velocity from the systemic mean value. New LOFAR and VLA radio maps uncover the BCG and a large red galaxy in the northeastern side to be head-tail galaxies, suggesting that this component has already traversed southwestern side and is now exiting the cluster to the northeast.

preprint2022arXiv

Signs of Higher Multipoles and Orbital Precession in GW151226

We present a reanalysis of GW151226, the second binary black hole merger discovered by the LIGO--Virgo Collaboration. Previous analysis showed that the best-fit waveform for this event corresponded to the merger of a $\sim 14 \, M_\odot$ black hole with a $\sim 7.5 \, M_\odot$ companion, and the posterior distribution in mass ratio ($q \leq 1$) is rather flat. In this work, we perform parameter estimation using a waveform model that includes the effects of orbital precession and higher-order radiative multipole modes, and we find that the source parameters of GW151226 shift towards the low $q$ and high effective spin ($χ_{\rm eff}$) region and that $q$ is better measured. The new solution has a log likelihood roughly two points higher than when either higher multipoles or orbital precession is neglected and can alter the astrophysical interpretation of GW151226. Additionally, we find it useful to use a flat-in-$χ{\rm eff}$ prior, which does not penalize the large $|χ_{\rm eff}|$ region, in order to uncover the higher likelihood region for GW151226. Our solution has several interesting properties: (a) the secondary black hole mass is close to the upper limit of the hypothesized lower mass gap of astrophysical black hole population; and (b) orbital precession is driven by the primary black hole spin, which has a dimensionless magnitude as large as $\sim 0.85$ and is tilted away from the orbital angular momentum at an angle of $\sim 57^\circ$. Since GW151226 is a relatively weak signal, an unambiguous claim of the detection of these effects in the signal cannot be made.

preprint2021arXiv

Mode-by-mode Relative Binning: Fast Likelihood Estimation for Gravitational Waveforms with Spin-Orbit Precession and Multiple Harmonics

Faster likelihood evaluation enhances the efficiency of gravitational wave signal analysis. We present Mode-by-mode Relative Binning (MRB), a new method designed for obtaining fast and accurate likelihoods for advanced waveform models that include spin-orbit precession effects and multiple radiation harmonics from compact binary coalescence. Leveraging the "twisting-up" procedure of constructing precessing waveform modes from non-precessing ones, the new method mitigates degrade of relative binning accuracy due to interference from superimposed modes. Additionally, we supplement algorithms for optimizing the choice of frequency bins specific to any given strain signal under analysis. Using the new method, we are able to evaluate the likelihood with up to an order of magnitude reduction in the number of waveform model calls per frequency compared to the previously used relative binning scheme, and achieve better likelihood accuracy than is sufficient for obtaining source parameter posterior distributions that are indistinguishable from the exact ones.

preprint2020arXiv

Asymmetric Surface Brightness Structure of Caustic Crossing Arc in SDSS J1226+2152: A Case for Dark Matter Substructure

We study the highly magnified arc SGAS J122651.3+215220 caused by a star-forming galaxy at $z_s=2.93$ crossing the lensing caustic cast by the galaxy cluster SDSS J1226+2152 ($z_l=0.43$), using Hubble Space Telescope observations. We report in the arc several asymmetric surface brightness features whose angular separations are a fraction of an arcsecond from the lensing critical curve and appear to be highly but unequally magnified image pairs of underlying compact sources, with one brightest pair having clear asymmetry consistently across four filters. One explanation of unequal magnification is microlensing by intracluster stars, which induces independent flux variations in the images of individual or groups of source stars in the lensed galaxy. For a second possibility, intracluster dark matter subhalos invisible to telescopes effectively perturb lensing magnifications near the critical curve and give rise to persistently unequal image pairs. Our modeling suggests, at least for the most prominent identified image pair, that the microlensing hypothesis is in tension with the absence of notable asymmetry variation over a six-year baseline, while subhalos of $\sim 10^6$--$10^8\,M_\odot$ anticipated from structure formation with Cold Dark Matter typically produce stationary and sizable asymmetries. We judge that observations at additional times and more precise lens models are necessary to stringently constrain temporal variability and robustly distinguish between the two explanations. The arc under this study is a scheduled target of a Director's Discretionary Early Release Science program of the James Webb Space Telescope, which will provide deep images and a high-resolution view with integral field spectroscopy.

preprint2020arXiv

New binary black hole mergers in the second observing run of Advanced LIGO and Advanced Virgo

We report the detection of new binary black hole merger events in the publicly available data from the second observing run of advanced LIGO and advanced Virgo (O2). The mergers were discovered using the new search pipeline described in Venumadhav et al. [Phys. Rev. D 100, 023011 (2019)], and are above the detection thresholds as defined in Abbott et al. (LIGO Scientific and Virgo Collaborations) [Phys. Rev. X 9, 031040 (2019).] Three of the mergers (GW170121, GW170304, GW170727) have inferred probabilities of being of astrophysical origin $p_{\rm astro} > 0.98$. The remaining three (GW170425, GW170202, GW170403) are less certain, with $p_{\rm astro}$ ranging from 0.5 to 0.8. The newly found mergers largely share the statistical properties of previously reported events, with the exception of GW170403, the least secure event, which has a highly negative effective spin parameter $χ_{\rm eff}$ . The most secure new event, GW170121 ($p_{\rm astro} > 0.99$), is also notable due to its inferred negative value of $χ_{\rm eff}$, which is inconsistent with being positive at the ~95.8% confidence level. The new mergers nearly double the sample of gravitational wave events reported from O2, and present a substantial opportunity to explore the statistics of the binary black hole population in the Universe. The number of detected events is not surprising since we estimate that the detection volume of our pipeline may be larger than that of other pipelines by as much as a factor of two (with significant uncertainties in the estimate). The increase in volume is larger when the constituent detectors of the network have very different sensitivities, as is likely to be the case in current and future runs.

preprint2020arXiv

Search for Lensed Gravitational Waves Including Morse Phase Information: An Intriguing Candidate in O2

We search for strongly lensed and multiply imaged gravitational wave signals in the second observing run of Advanced LIGO and Advanced Virgo (O2). We exploit a new source of information, the so-called Morse phase, which further mitigates the search background and constrains viable lenses. The best candidate we find is consistent with a strongly lensed signal from a massive binary black hole (BBH) merger, with three detected images consisting of the previously catalogued events GW170104 and GW170814, and a subthreshold trigger, GWC170620. Given the number of BBH events detected so far, we estimate an overall false alarm probability $\sim 10^{-4}$ for the observed high degree of parameter coincidence between the three events. On the flip side, we measure the Morse phase differences which suggest a complex and atypical lens system, with at least five images including a magnified image at a local maximum of the Fermat potential. The low prior probability for multiple lensed images and the amount of fine tuning required in the lens model reduce the credibility of the lensing hypothesis. The long time delays between lensed images point toward a galaxy cluster lens with an internal velocity dispersion $σ\sim 650\,{\rm km/s}$, and the observed strain amplitudes imply a likely range $0.4 < z \lesssim 0.7$ for the source redshift. We provide an error ellipse of $\sim 16\,{\rm deg}^2$ for the sky location of the source together with additional specific constraints on the lens-host system, and encourage follow-up efforts to confirm or rule out any viable lens. If this is indeed a lensed event, successfully pinpointing the system would offer a unique opportunity to identify the host galaxy of a BBH merger, and even localize the source within it.

preprint2020arXiv

Statistical Microlensing Toward Magnified High-Redshift Star Clusters

We study light variability of gravitationally magnified high-redshift star clusters induced by a foreground population of microlenses. This arises as the incoherent superposition of light variations from a large number of source stars traversing the random magnification pattern on the source plane. The light curve resembles a scale-invariant, Gaussian process on timescales of years to decades, while exhibits rapid and frequent micro-caustic crossing flares of larger amplitudes on timescales of days to months. For a concrete example, we study a young Lyman-continuum-leaking star cluster recently discovered in the lensed Sunburst Arc at $z=2.37$. We show that one magnified image happens to be intervened by a faint foreground galaxy, and hence should exhibit a variable flux at the $1$--$2\%$ level, which is measurable in space with $\sim 1$--$3\,$ks exposures on the Hubble Space Telescope and more easily with the James Webb Space Telescope, or from the ground using a $\sim$4-meter telescope without adaptive optics. Detailed measurement of this variability will enable us to determine the absolute macro magnification and hence the intrinsic mass and length scales of the star cluster, test synthetic models of stellar population, and probe multiplicity of massive stars. We furthermore suggest that monitoring the other lensed images of the star cluster, which are free from significant intervention by foreground stellar microlenses, will allow us to probe planetary to stellar mass compact objects constituting as little as just a few percent of the dark matter. Given the typical surface density of intracluster stars, we expect this phenomenon to be relevant for many other gravitationally magnified star clusters at Cosmic Noon behind galaxy cluster lenses.

preprint2020arXiv

Testing exchangeability with martingale for change-point detection

This work proposes a new exchangeability test for a random sequence through a martingale based approach. Its main contributions include: 1) an additive martingale which is more amenable for designing exchangeability tests by exploiting the Hoeffding-Azuma lemma; 2) different betting functions for constructing the additive martingale are studied. By choosing the underlying probability density function of p-values as a betting function, it can be shown that, when a change-point appears, a satisfying trade-off between the smoothness and expected one-step increment of the martingale sequence can be obtained. An online algorithm based on Beta distribution parametrization for constructing this betting function is discussed in detail as well.

preprint2019arXiv

Gravitational Lensing Signatures of Axion Dark Matter Minihalos in Highly Magnified Stars

Axions are a viable candidate for Cold Dark Matter (CDM) which should generically form minihalos of sub-planetary masses from white-noise isocurvature density fluctuations if the Peccei-Quinn phase transition occurs after inflation. Despite being denser than the larger halos formed out of adiabatic fluctuations from inflation, axion minihalos have surface densities much smaller than the critical value required for gravitational lensing to produce multiple images or high magnification, and hence are practically undetectable as lenses in isolation. However, their lensing effect can be enhanced when superposed near critical curves of other lenses. We propose a method to detect them through photometric monitoring of recently discovered caustic transiting stars behind cluster lenses, under extreme magnification factors $μ\gtrsim 10^3$--$10^4$ as the lensed stars cross microlensing caustics induced by intracluster stars. For masses of the first gravitationally collapsed minihalos in the range $\sim 10^{-15}$--$10^{-8}\,h^{-1}\,M_\odot$, we show that axion minihalos in galaxy clusters should collectively produce subtle surface density fluctuations of amplitude $\sim 10^{-4}$--$10^{-3}$ on projected length scales of $\sim 10$--$10^4\,$AU, which imprint irregularities in the microlensing light curves of caustic transiting stars. We estimate that, inside a cluster halo and over the age of the Universe, most of these minihalos are likely to avoid dynamic disruption by encounters with stars or other minihalos.

preprint2016arXiv

Antisymmetric galaxy cross-correlations as a cosmological probe

The auto-correlation between two members of a galaxy population is symmetric under the interchange of the two galaxies being correlated. The cross-correlation between two different types of galaxies, separated by a vector $\bf{r}$, is not necessarily the same as that for a pair separated by $-\bf{r}$. Local anisotropies in the two-point cross-correlation function may thus indicate a specific direction which when mapped as a function of position trace out a vector field. This vector field can then be decomposed into longitudinal and transverse components, and those transverse components written as positive- and negative-helicity components. A locally asymmetric cross-correlation of the longitudinal type arises naturally in halo clustering, even with Gaussian initial conditions, and could be enhanced with local-type non-Gaussianity. Early-Universe scenarios that introduce a vector field may also give rise to such effects. These antisymmetric cross-correlations also provide a new possibility to seek a preferred cosmic direction correlated with the hemispherical power asymmetry in the cosmic microwave background and to seek a preferred location associated with the CMB cold spot. New ways to seek cosmic parity breaking are also possible.

preprint2016arXiv

Lensing of Fast Radio Bursts as a Probe of Compact Dark Matter

The possibility that part of the dark matter is made of massive compact halo objects (MACHOs) remains poorly constrained over a wide range of masses, and especially in the $20-100\, M_\odot$ window. We show that strong gravitational lensing of extragalactic fast radio bursts (FRBs) by MACHOs of masses larger than $\sim20\,M_\odot$ would result in repeated FRBs with an observable time delay. Strong lensing of an FRB by a lens of mass $M_L$ induces two images, separated by a typical time delay $\sim$ few $\times(M_L/30\, M_\odot)$ milliseconds. Considering the expected FRB detection rate by upcoming experiments, such as CHIME, of $10^4$ FRBs per year, we should observe from tens to hundreds of repeated bursts yearly, if MACHOs in this window make up all the dark matter. A null search for echoes with just $10^4$ FRBs, would constrain the fraction $f_{\rm DM}$ of dark matter in MACHOs to $f_{\rm DM}\lesssim 0.08$ for $M_L\gtrsim 20\,M_\odot$.

preprint2016arXiv

Search for Compensated Isocurvature Perturbations with Planck Power Spectra

In the standard inflationary scenario, primordial perturbations are adiabatic. The amplitudes of most types of isocurvature perturbations are generally constrained by current data to be small. If, however, there is a baryon-density perturbation that is compensated by a dark-matter perturbation in such a way that the total matter density is unperturbed, then this compensated isocurvature perturbation (CIP) has no observable consequence in the cosmic microwave background (CMB) at linear order in the CIP amplitude. Here we search for the effects of CIPs on CMB power spectra to quadratic order in the CIP amplitude. An analysis of the Planck temperature data leads to an upper bound $Δ_{\rm rms}^2 \leq 7.1\times 10^{-3}$, at the 68\% confidence level, to the variance $Δ_{\rm rms}^2$ of the CIP amplitude. This is then strengthened to $Δ_{\rm rms}^2\leq 5.0\times 10^{-3}$ if Planck small-angle polarization data are included. A cosmic-variance-limited CMB experiment could improve the $1σ$ sensitivity to CIPs to $Δ^2_{\rm rms} \lesssim 9\times 10^{-4}$. It is also found that adding CIPs to the standard $Λ$CDM model can improve the fit of the observed smoothing of CMB acoustic peaks just as much as adding a non-standard lensing amplitude.

preprint2016arXiv

Sequential Monte Carlo Methods for System Identification

One of the key challenges in identifying nonlinear and possibly non-Gaussian state space models (SSMs) is the intractability of estimating the system state. Sequential Monte Carlo (SMC) methods, such as the particle filter (introduced more than two decades ago), provide numerical solutions to the nonlinear state estimation problems arising in SSMs. When combined with additional identification techniques, these algorithms provide solid solutions to the nonlinear system identification problem. We describe two general strategies for creating such combinations and discuss why SMC is a natural tool for implementing these strategies.

preprint2015arXiv

Conformal Fermi Coordinates

Fermi Normal Coordinates (FNC) are a useful frame for isolating the locally observable, physical effects of a long-wavelength spacetime perturbation. Their cosmological application, however, is hampered by the fact that they are only valid on scales much smaller than the horizon. We introduce a generalization that we call Conformal Fermi Coordinates (CFC). CFC preserve all the advantages of FNC, but in addition are valid outside the horizon. They allow us to calculate the coupling of long- and short-wavelength modes on all scales larger than the sound horizon of the cosmological fluid, starting from the epoch of inflation until today, by removing the complications of the second order Einstein equations to a large extent, and eliminating all gauge ambiguities. As an application, we present a calculation of the effect of long-wavelength tensor modes on small scale density fluctuations. We recover previous results, but clarify the physical content of the individual contributions in terms of locally measurable effects and "projection" terms.

preprint2015arXiv

Detecting the integrated Sachs-Wolfe effect with high-redshift 21-cm surveys

We investigate the possibility to detect the integrated Sachs-Wolfe (ISW) effect by cross-correlating 21-cm surveys at high redshifts with galaxies, in a way similar to the usual CMB-galaxy cross-correlation. The high-redshift 21-cm signal is dominated by CMB photons that travel freely without interacting with the intervening matter, and hence its late-time ISW signature should correlate extremely well with that of the CMB at its peak frequencies. Using the 21-cm temperature brightness instead of the CMB would thus be a further check of the detection of the ISW effect, measured by different instruments at different frequencies and suffering from different systematics. We also study the ISW effect on the photons that are scattered by HI clouds. We show that a detection of the unscattered photons is achievable with planned radio arrays, while one using scattered photons will require advanced radio interferometers, either an extended version of the planned Square Kilometre Array or futuristic experiments such as a lunar radio array.

preprint2015arXiv

On Separate Universes

(abridged version) The separate universe conjecture states that in General Relativity a density perturbation behaves locally (i.e. on scales much smaller than the wavelength of the mode) as a separate universe with different background density and curvature. We prove this conjecture for a spherical compensated tophat density perturbation of arbitrary amplitude and radius in $Λ$CDM. We then use Conformal Fermi Coordinates to generalize this result to scalar perturbations of arbitrary configuration and scale. In this case, the separate universe conjecture holds for the isotropic part of the perturbations. The anisotropic part on the other hand is exactly captured by a tidal field in the Newtonian form. We show that the separate universe picture is restricted to scales larger than the sound horizons of all fluid components. We then derive an expression for the locally measured matter bispectrum induced by a long-wavelength mode of arbitrary wavelength. We show that nonlinear gravitational dynamics does not generate observable contributions that scale like local-type non-Gaussianity $f_{\rm NL}^{\rm loc}$, and hence does not contribute to a scale-dependent galaxy bias $Δb \propto k^{-2}$ on large scales; rather, the locally measurable long-short mode coupling assumes a form essentially identical to subhorizon perturbation theory results, once the long-mode density perturbation is replaced by the synchronous-comoving gauge density perturbation. Apparent $f_{\rm NL}^{\rm loc}$-type contributions arise through projection effects on photon propagation, which depend on the specific large-scale structure tracer and observable considered, and are in principle distinguishable from the local mode coupling induced by gravity. We conclude that any observation of $f_{\rm NL}^{\rm loc}$ beyond these projection effects signals a departure from standard single-clock inflation.

preprint2015arXiv

On the exponential convergence of the Kaczmarz algorithm

The Kaczmarz algorithm (KA) is a popular method for solving a system of linear equations. In this note we derive a new exponential convergence result for the KA. The key allowing us to establish the new result is to rewrite the KA in such a way that its solution path can be interpreted as the output from a particular dynamical system. The asymptotic stability results of the corresponding dynamical system can then be leveraged to prove exponential convergence of the KA. The new bound is also compared to existing bounds.

preprint2015arXiv

Revisiting the anomalous bending elasticity of sharply bent DNA

Several recent experiments suggest that sharply bent DNA has a surprisingly high bending flexibility, but the cause of this flexibility is poorly understood. Although excitation of flexible defects can explain these results, whether such excitation can occur with the level of DNA bending in these experiments remains unclear. Intriguingly, the DNA contained preexisting nicks in most of these experiments but whether nicks might play a role in flexibility has never been considered in the interpretation of experimental results. Here, using full-atom molecular dynamics simulations, we show that nicks promote DNA basepair disruption at the nicked sites, which drastically reduces DNA bending energy. In addition, lower temperatures suppress the nick-dependent basepair disruption. In the absence of nicks, basepair disruption can also occur but requires a higher level of DNA bending. Therefore, basepair disruption inside B-form DNA can be suppressed if the DNA contains preexisting nicks. Overall, our results suggest that the reported mechanical anomaly of sharply bent DNA is likely dependent on preexisting nicks, therefore the intrinsic mechanisms of sharply bent nick-free DNA remain an open question.

preprint2014arXiv

Linking the BICEP2 result and the hemispherical power asymmetry through spatial variation of $r$

For the simplest inflation models, the recent detection of a large primordial B-mode polarization signal by the BICEP2 experiment indicates a slight tension with the upper limit on the tensor-to-scalar ratio, $r$, from the Planck satellite. Here, we discuss spatially varying $r$ as a possible explanation for this discrepancy. This idea seems attractive since it may also explain part of the hemispherical temperature power asymmetry seen by WMAP and Planck at large angular scales. If these two aspects are indeed connected, the model suggests that in the Northern hemisphere $r$ should be much smaller, a hypothesis that could be confirmed with future B-mode experiments, providing a test for the stationarity of primordial tensor contributions across the sky. The BICEP2 measurement furthermore rules out that a simple dipolar modulation of $r$ alone can be responsible for the full hemispherical power asymmetry.

preprint2014arXiv

New operator approach to the CMB aberration kernels in harmonic space

Aberration kernels describe how harmonic-space multipole coefficients of cosmic microwave background (CMB) observables transform under Lorentz boosts of the reference frame. For spin-weighted CMB observables, transforming like the CMB temperature (i.e. Doppler weight $d = 1$), we show that the aberration kernels are the matrix elements of a unitary boost operator in harmonic space. Algebraic properties of the rotation and boost generators then give simple, exact recursion relations that allow us to raise or lower the multipole quantum numbers $\ell$ and $m$, and the spin weight $s$. Further recursion relations express kernels of other Doppler weights $d \neq 1$ in terms of the $d = 1$ kernels. From those we show that on the full sky, to all orders in $β=v/c$, $E$- and $B$-mode polarization observables do not mix under aberration if and only if $d = 1$. The new relations, fully non-linear in the boost velocity $β$, form the basis of a practical recursive algorithm to accurately compute any aberration kernel. In addition, we develop a second, fast algorithm in which aberration kernels are obtained using a set of ordinary differential equations. This system can also be approximately solved at small scales, providing simple asymptotic formulae for the aberration kernels. The results of this work will be useful for further studying the effect of aberration on future CMB temperature and polarization analysis, and might provide a basis for relativistic radiative transfer schemes.

preprint2014arXiv

On the Nuclear Norm heuristic for a Hankel matrix Recovery Problem

This note addresses the question if and why the nuclear norm heuristic can recover an impulse response generated by a stable single-real-pole system, if elements of the upper-triangle of the associated Hankel matrix were given. Since the setting is deterministic, theories based on stochastic assumptions for low-rank matrix recovery do not apply here. A 'certificate' which guarantees the completion is constructed by exploring the structural information of the hidden matrix. Experimental results and discussions regarding the nuclear norm heuristic applied to a more general setting are also given.

preprint2014arXiv

Reheating constraints to inflationary models

Evidence from the BICEP2 experiment for a significant gravitational-wave background has focussed attention on inflaton potentials $V(ϕ) \propto ϕ^α$ with $α=2$ ("chaotic" or "$m^2ϕ^2$" inflation) or with smaller values of $α$, as may arise in axion-monodromy models. Here we show that reheating considerations may provide additional constraints to these models. The reheating phase preceding the radiation era is modeled by an effective equation-of-state parameter $w_{\rm re}$. The canonical reheating scenario is then described by $w_{\rm re}=0$. The simplest $α=2$ models are consistent with $w_{\rm re} = 0$ for values of $n_s$ well within the current $1σ$ range. Models with $α=1$ or $α=2/3$ require a more exotic reheating phase, with $-1/3<w_{\rm re}<0$, unless $n_s$ falls above the current $1σ$ range. Likewise, models with $α=4$ require a physically implausible $w_{\rm re}>1/3$, unless $n_s$ is close to the lower limit of the $2σ$ range. For $m^2ϕ^2$ inflation and canonical reheating as a benchmark, we derive a relation $\log_{10}\left(T_{\rm re}/10^6\,{\rm GeV} \right) \simeq 2000\,(n_s-0.96)$ between the reheat temperature $T_{\rm re}$ and the scalar spectral index $n_s$. Thus, if $n_s$ is close to its central value, then $T_{\rm re}\lesssim 10^6$~GeV, just above the electroweak scale. If the reheat temperature is higher, as many theorists may prefer, then the scalar spectral index should be closer to $n_s\simeq0.965$ (at the pivot scale $k=0.05\,{\rm Mpc}^{-1}$), near the upper limit of the $1σ$ error range. Improved precision in the measurement of $n_s$ should allow $m^2ϕ^2$, axion-monodromy, and $ϕ^4$ models to be distinguished, even without precise measurement of $r$, and to test the $m^2ϕ^2$ expectation of $n_s\simeq0.965$.

preprint2014arXiv

Rotation of the cosmic microwave background polarization from weak gravitational lensing

When a cosmic microwave background (CMB) photon travels from the surface of last scatter through spacetime metric perturbations, the polarization vector may rotate about its direction of propagation. This gravitational rotation is distinct from, and occurs in addition to, the lensing deflection of the photon trajectory. This rotation can be sourced by linear vector or tensor metric perturbations and is fully coherent with the curl deflection field. Therefore, lensing corrections to the CMB polarization power spectra as well as the temperature-polarization cross-correlations due to non-scalar perturbations are modified. The rotation does not affect lensing by linear scalar perturbations, but needs to be included when calculations go to higher orders. We present complete results for weak lensing of the full-sky CMB power spectra by general linear metric perturbations, taking into account both deflection of the photon trajectory and rotation of the polarization. For the case of lensing by gravitational waves, we show that the B modes induced by the rotation largely cancel those induced by the curl component of deflection.

preprint2014arXiv

Sparse Estimation From Noisy Observations of an Overdetermined Linear System

This note studies a method for the efficient estimation of a finite number of unknown parameters from linear equations, which are perturbed by Gaussian noise. In case the unknown parameters have only few nonzero entries, the proposed estimator performs more efficiently than a traditional approach. The method consists of three steps: (1) a classical Least Squares Estimate (LSE), (2) the support is recovered through a Linear Programming (LP) optimization problem which can be computed using a soft-thresholding step, (3) a de-biasing step using a LSE on the estimated support set. The main contribution of this note is a formal derivation of an associated ORACLE property of the final estimate. That is, when the number of samples is large enough, the estimate is shown to equal the LSE based on the support of the {\em true} parameters.

preprint2014arXiv

Spectral distortions from the dissipation of tensor perturbations

Spectral distortions of the cosmic microwave background (CMB) may become a powerful probe of primordial perturbations at small scales. Existing studies of spectral distortions focus almost exclusively on primordial scalar metric perturbations. Similarly, vector and tensor perturbations should source CMB spectral distortions. In this paper, we give general expressions for the effective heating rate caused by these types of perturbations, including previously neglected contributions from polarization states and higher multipoles. We then focus our discussion on the dissipation of tensors, showing that for nearly scale invariant tensor power spectra, the overall distortion is some six orders of magnitudes smaller than from the damping of adiabatic scalar modes. We find simple analytic expressions describing the effective heating rate from tensors using a quasi-tight coupling approximation. In contrast to adiabatic modes, tensors cause heating without additional photon diffusion and thus over a wider range of scales, as recently pointed out by Ota et. al 2014. Our results are in broad agreement with their conclusions, but we find that small-scale modes beyond k< 2x10^4 Mpc^{-1} cannot be neglected, leading to a larger distortion, especially for very blue tensor power spectra. At small scales, also the effect of neutrino damping on the tensor amplitude needs to be included.

preprint2014arXiv

Tensor-induced B modes with no temperature fluctuations

The recent indications for a tensor-to-scalar ratio $ r \simeq 0.2$ from BICEP2 measurements of the cosmic microwave background (CMB) B-mode polarization present some tension with upper limits $r \lesssim 0.1$ from measurements of CMB temperature fluctuations. Here we point out that tensor perturbations can induce B modes in the CMB polarization without inducing any temperature fluctuations nor E-mode polarization whatsoever, but only, at the expense of violating the Copernican principle. We present this mathematical possibility as a new ingredient for the model-builder's toolkit in case the tension between B modes and temperature fluctuations cannot be resolved with more conventional ideas.

preprint2014arXiv

The redshift-space galaxy two-point correlation function and baryon acoustic oscillations

Future galaxy surveys will measure baryon acoustic oscillations (BAOs) with high significance, and a complete understanding of the anisotropies of BAOs in redshift space will be important to exploit the cosmological information in BAOs. Here we describe the anisotropies that arise in the redshift-space galaxy two-point correlation function (2PCF) and elucidate the origin of features that arise in the dependence of the BAOs on the angle between the orientation of the galaxy pair and the line of sight. We do so with a derivation of the configuration-space 2PCF using streaming model. We find that, contrary to common belief, the locations of BAO peaks in the redshift-space 2PCF are anisotropic even in the linear theory. Anisotropies in BAO depend strongly on the method of extracting the peak, showing maximum 3 % angular variation. We also find that extracting the BAO peak of $r^2ξ(r,μ)$ significantly reduces the anisotropy to sub-percent level angular variation. When subtracting the tilt due to the broadband behavior of the 2PCF, the BAO bump is enhanced along the line of sight because of local infall velocities toward the BAO bump. Precise measurement of the angular dependence of the redshift-space 2PCF will allow new geometrical tests of dark energy beyond the BAO.

preprint2013arXiv

Anisotropic imprint of long-wavelength tensor perturbations on cosmic structure

Inflationary models predict a correlation between primordial density perturbations (scalar metric perturbations) and gravitational waves (tensor metric perturbations) in the form of a scalar-scalar-tensor three-point correlation, or bispectrum in Fourier space. The squeezed limit of this bispectrum implies a quadrupolar asymmetry in the observed local power spectrum for matter and galaxies. Here we show (like others before) that an infrared divergence in the amplitude of this power asymmetry predicted in single-field slow-roll models is canceled by projection effects when considering the observed power spectrum. We then further evaluate the nonzero, but finite, residual quadrupolar power asymmetry that remains after the divergences are canceled. While the quadrupolar power asymmetry is small, it is conceptually important. Our calculation moreover clarifies how the predictions for this power asymmetry may change in models with different scalar-scalar-tensor bispectra, and shows that convincing detection of the quadrupolar power asymmetry would rule out the single-field slow-roll models of inflation.

preprint2013arXiv

Multiple scattering Sunyaev-Zeldovich signal I: lowest order effect

Future high-resolution, high-sensitivity Sunyaev-Zeldovich (SZ) observations of individual clusters will allow us to study the dynamical state of the intra cluster medium (ICM). To exploit the full potential of this new observational window, it is crucial to understand the different contributions to the total SZ signal. Here, we investigate the signature caused by multiple scatterings at lowest order of the electron temperature. Previous analytic discussions of this problem used the isotropic scattering approximation (ISA), which even for the simplest cluster geometries is rather rough. We take a step forward and consistently treat the anisotropy of the ambient radiation field caused by the first scattering. We show that the multiple scattering SZ signal directly probes line of sight anisotropies of the ICM, thereby delivering a new set of SZ observables which could be used for 3D cluster-profile reconstruction. The multiple scattering signal should furthermore correlate spatially with the cluster's X-ray and SZ polarization signals, an effect that could allow enhancing the detectability of this contribution.

preprint2013arXiv

Multiple scattering Sunyaev-Zeldovich signal II: relativistic effects

We study the multiple scattering Sunyaev-Zeldovich (SZ) signature, extending our previous analysis to high-temperature clusters. We consistently treat the anisotropy of the ambient radiation field caused by the first scattering and also consider lowest order kinematic terms. We show that due to temperature corrections monopole through octupole anisotropy of the singly scattered SZ signal attain different spectra in the second scattering. The difference becomes more pronounced at high temperature, and thus could be used to constrain individual line of sight moments of the electron density and temperature profiles. While very challenging from the observational point of view, this further extends the list of possible SZ observables that will be important for 3D cluster-profile reconstruction, possibly helping to break geometric degeneracies caused by projection effects. We also briefly discuss the scattering of primordial CMB anisotropies by SZ clusters.

preprint2013arXiv

Seeking Inflation Fossils in the Cosmic Microwave Background

If during inflation the inflaton couples to a "fossil" field, some new scalar, vector, or tensor field, it typically induces a scalar-scalar-fossil bispectrum. Even if the fossil field leaves no direct physical trace after inflation, it gives rise to correlations between different Fourier modes of the curvature or, equivalently, a nonzero curvature trispectrum, but without a curvature bispectrum. Here we quantify the effects of a fossil field on the cosmic microwave background (CMB) temperature fluctuations in terms of bipolar spherical harmonics (BiPoSHs). The effects of vector and tensor fossils can be distinguished geometrically from those of scalars through the parity of the BiPoSHs they induce. However, the two-dimensional nature of the CMB sky does not allow vectors to be distinguished geometrically from tensors. We estimate the detectability of a signal in terms of the scalar-scalar-fossil coupling for scalar, vector, and tensor fossils, assuming a local-type coupling. We comment on a divergence that arises in the quadrupolar BiPoSH from the scalar-scalar-tensor correlation in single-field slow-roll inflation.

preprint2013arXiv

The effect of aberration on partial-sky measurements of the cosmic microwave background temperature power spectrum

Our motion relative to the cosmic-microwave-background (CMB) rest frame deflects light rays giving rise to shifts as large as L -> L(1+beta), where beta=0.00123 is our velocity (in units of the speed of light) on measurements CMB fluctuations. Here we present a novel harmonic-space approach to this CMB aberration that improves upon prior work by allowing us to (i) go to higher orders in beta, thus extending the validity of the analysis to measurements at L > 1/beta ~ 800; and (ii) treat the effects of window functions and pixelization in a more accurate and computationally efficient manner. We calculate precisely the magnitude of the systematic bias in the power spectrum inferred from the partial sky, and show that aberration shifts the multipole moment by Delta L/L ~ beta<cos(theta)>, with <cos(theta)> averaged over the survey footprint. Such a shift, if ignored, would bias the measurement of the sound-horizon size theta_* at the 0.01%-level, which is comparable to the measurement uncertainties of Planck. The bias can then propagate into cosmological parameters such as the angular-diameter distance, Hubble parameter and dark-energy equation of state. We study the effect of aberration for current Planck, South Pole Telescope (SPT) and Atacama Cosmology Telescope (ACT) data and show that the bias cannot be neglected. We suggest that the small tension between Planck, ACT, and SPT may be due partially to aberration. An Appendix shows how the near constancy of the full-sky power spectrum under aberration follows from unitarity of the aberration kernel.

preprint2013arXiv

The Pesky Power Asymmetry

Physical models for the hemispherical power asymmetry in the cosmic microwave background (CMB) reported by the Planck Collaboration must satisfy CMB constraints to the homogeneity of the Universe and quasar constraints to power asymmetries. We survey a variety of models for the power asymmetry and show that consistent models include a modulated scale-dependent isocurvature contribution to the matter power spectrum or a modulation of the reionization optical depth, gravitational-wave amplitude, or scalar spectral index. We propose further tests to distinguish between the different scenarios.

preprint2013arXiv

Wigner-Eckart theorem in cosmology: Bispectra for total-angular-momentum waves

Total-angular-momentum (TAM) waves provide a set of basis functions for scalar, vector, and tensor fields that can be used in place of plane waves and that reflect the rotational symmetry of the spherical sky. Here we discuss three-point correlation functions, or bispectra in harmonic space, for scalar, vector, and tensor fields in terms of TAM waves. The Wigner-Eckart theorem dictates that the expectation value, assuming statistical isotropy, of the product of three TAM waves is the product of a Clebsch-Gordan coefficient (or Wigner-3j symbol) times a function only of the total-angular-momentum quantum numbers. Here we show how this works, and we provide explicit expressions relating the bispectra for TAM waves in terms of the more commonly used Fourier-space bispectra. This formalism will be useful to simplify calculations of projections of three-dimensional bispectra onto the spherical sky.

preprint2012arXiv

Total Angular Momentum Waves for Scalar, Vector, and Tensor Fields

Most calculations in cosmological perturbation theorydecompose those perturbations into plane waves (Fourier modes). However, for some calculations, particularly those involving observations performed on a spherical sky, a decomposition into waves of fixed total angular momentum (TAM) may be more appropriate. Here we introduce TAM waves, solutions of fixed total angular momentum to the Helmholtz equation, for three-dimensional scalar, vector, and tensor fields. The vector TAM waves of given total angular momentum can be decomposed further into a set of three basis functions of fixed orbital angular momentum (OAM), a set of fixed helicity, or a basis consisting of a longitudinal (L) and two transverse (E and B) TAM waves. The symmetric traceless rank-2 tensor TAM waves can be similarly decomposed into a basis of fixed OAM or fixed helicity, or a basis that consists of a longitudinal (L), two vector (VE and VB, of opposite parity), and two tensor (TE and TB, of opposite parity) waves. We show how all of the vector and tensor TAM waves can be obtained by applying derivative operators to scalar TAM waves. This operator approach then allows one to decompose a vector field into three covariant scalar fields for the L, E, and B components and symmetric-traceless-tensor fields into five covariant scalar fields for the L, VE, VB, TE, and TB components. We provide projections of the vector and tensor TAM waves onto vector and tensor spherical harmonics. We provide calculational detail to facilitate the assimilation of this formalism into cosmological calculations. As an example, we calculate the power spectra of the deflection angle for gravitational lensing by density perturbations and by gravitational waves. We comment on an alternative approach to CMB fluctuations based on TAM waves. Our work may have applications elsewhere in field theory and in general relativity.

preprint2012arXiv

Tree amplitudes and color decomposition in broken SU(2)

We propose a color decomposition for general tree amplitudes in a SU(2) gauge theory which is spontaneously broken via the Higgs mechanism. Working in the unitary gauge, we construct color-ordered amplitudes by explicitly presenting a set of color-ordered Feynman rules. Those primitive amplitudes are gauge-invariant, and they preserve perturbative unitarity in the high-energy limit. Serving as building blocks of color-dressed tree amplitudes, they allow for efficient evaluation of tree-level scattering amplitudes involving gauge bosons and the Higgs boson via the Berends-Giele recursion relations for color-ordered currents. We demonstrate the efficiency of this computational scheme by calculating on-shell amplitudes for scattering of five, six and nine W-bosons in the limit of vanishing Weinberg angle.

preprint2011arXiv

Next-to-leading order QCD predictions for $A^{0}γ$ associated production at the CERN Large Hadron Collider

We calculate the complete next-to-leading-order (NLO) QCD corrections (including SUSY QCD corrections) to the inclusive total cross sections of the associated production processes $pp\rightarrow A^{0}γ+X$ in the minimal supersymmetric standard model (MSSM) at the CERN Large Hadron Collider (LHC). Our results show that the enhancement of the total cross sections from the NLO QCD corrections can reach $15%\sim20%$ for 200 GeV$<m_{A}<500$ GeV and $\tanβ=50$. The scale dependence of the total cross section is improved by the NLO corrections, which is less than 5%. We also show the Monte Carlo simulation results for the $τ^{+}τ^{-}+γ$ signature including the complete NLO QCD effects, and find an observable signature above the standard model (SM) background for a normal luminosity of 100 fb$^{-1}$ at the LHC.

preprint2011arXiv

One-loop Helicity Amplitudes for Top Quark Pair Production in Randall-Sundrum Model

In this paper, we show how to calculate analytically the one-loop helicity amplitudes for the process $q\bar{q} rightarrow t\bar{t}$ induced by KK gluon, using the spinor-helicity formalism. A minimal set of Feynman rules which are uniquely fixed by gauge invariance and the color representation of the KK gluon are derived and used in the calculation. Our results can be applied to a variety of models containing a massive color octet vector boson.