title: Gap solitons for the repulsive Gross-Pitaevskii equation with periodic potential: coding and method for computation

abstract: The paper is devoted to nonlinear localized modes ("gap solitons") for the spatially one-dimensional Gross-Pitaevskii equation (1D GPE) with a periodic potential and repulsive interparticle interactions. It has been recently shown (G. L. Alfimov, A. I. Avramenko, Physica D, 254, 29 (2013)) that under certain conditions all the stationary modes for the 1D GPE can be coded by bi-infinite sequences of symbols of some finite alphabet (called "codes" of the solutions). We present and justify a numerical method which allows to reconstruct the profile of a localized mode by its code. As an example, the method is applied to compute the profiles of gap solitons for 1D GPE with a cosine potential.

categories: nlin.PS cond-mat.quant-gas math-ph math.MP nlin.CD


title: On the quantum principles of cognitive learning

abstract: Pothos & Busemeyer's (P&B's) query about whether quantum probability can provide a foundation for the cognitive modeling embodies so many underlying implications that the subject is far from exhausted. In this brief commentary, however, I suggest that the conceptual thresholds of the meaningful learning give rise to a typical Boltzmann's weighting measure, which indicates a statistical verisimilitude of quantum behavior in the human cognitive ensemble.


title: Rotation sensing with trapped ions

abstract: We present a protocol for using trapped ions to measure rotations via matter-wave Sagnac interferometry. The trap allows the interferometer to enclose a large area in a compact apparatus through repeated round-trips in a Sagnac geometry. We show how a uniform magnetic field can be used to close the interferometer over a large dynamic range in rotation speed and measurement bandwidth without losing contrast. Since this technique does not require the ions to be confined in the Lamb-Dicke regime, thermal states with many phonons should be sufficient for operation.


title: Gibbs states defined by biorthogonal sequences

abstract: Motivated by the growing interest on PT-quantum mechanics, in this paper we discuss some facts on generalized Gibbs states and on their related KMS-like conditions. To achieve this, we first consider some useful connections between similar (Hamiltonian) operators and we propose some extended version of the Heisenberg algebraic dynamics, deducing some of their properties, useful for our purposes.


title: Localization by Fusing a Group of Fingerprints via Multiple Antennas in Indoor Environment

abstract: Most existing fingerprints-based indoor localization approaches are based on some single fingerprints, such as received signal strength (RSS), channel impulse response (CIR), and signal subspace. However, the localization accuracy obtained by the single fingerprint approach is rather susceptible to the changing environment, multi-path, and non-line-of-sight (NLOS) propagation. Furthermore, building the fingerprints is a very time consuming process. In this paper, we propose a novel localization framework by Fusing A Group Of fingerprinTs (FAGOT) via multiple antennas for the indoor environment. We first build a GrOup Of Fingerprints (GOOF), which includes five different fingerprints, namely, RSS, covariance matrix, signal subspace, fractional low order moment, and fourth-order cumulant, which are obtained by different transformations of the received signals from multiple antennas in the offline stage. Then, we design a parallel GOOF multiple classifiers based on AdaBoost (GOOF-AdaBoost) to train each of these fingerprints in parallel as five strong multiple classifiers. In the online stage, we input the corresponding transformations of the real measurements into these strong classifiers to obtain independent decisions. Finally, we propose an efficient combination fusion algorithm, namely, MUltiple Classifiers mUltiple Samples (MUCUS) fusion algorithm to improve the accuracy of localization by combining the predictions of multiple classifiers with different samples. As compared with the single fingerprint approaches, the prediction probability of our proposed approach is improved significantly. The process for building fingerprints can also be reduced drastically. We demonstrate the feasibility and performance of the proposed algorithm through extensive simulations as well as via real experimental data using a Universal Software Radio Peripheral (USRP) platform with four antennas.


title: Escape rate and diffusion of a random walker

abstract: We determine the rate of escape from a potential well, and the diffusion coefficient in a periodic potential, of a random walker that moves under the influence of the potential in between successive collisions with the heat bath. In the overdamped limit, both the escape rate and the diffusion coefficient coincide with those of a Langevin particle. Conversely, in the underdamped limit the two dynamics have a different temperature dependence. In particular, at low temperature the random walk has a smaller escape rate, but a larger diffusion coefficient.


title: Semiclassical theory of anisotropic transport at LaAlO3/SrTiO3 interfaces under in-plane magnetic field

abstract: The unconventional magnetotransport at the interface between transition-metal oxides $LaAlO_3$ (LAO) and $SrTiO_3$ (STO) is frequently related to mobile electrons interacting with localized magnetic moments. However nature and properties of magnetism at this interface are not well understood so far. In this paper, we focus on transport effects driven by spin-orbit coupling and intentionally neglect possible strong correlations. The electrical resistivity tensor is calculated as a function of the magnitude and orientation of an external magnetic field parallel to the interface. The semiclassical Boltzmann equation is solved numerically for the two-dimensional system of spin-orbit coupled electrons accelerated by an electric field and scattered by spatially-correlated impurities. At temperatures of a few Kelvin and densities such that the chemical potential crosses the second pair of spin-orbit split bands, we find a strongly anisotropic modulation of the (negative) magnetoresistance above 10 T, characterized by multiple maxima and minima away from the crystalline axes. Along with the drop of the magnetoresistance, an abrupt enhancement of the transverse resistivity occurs. The angular modulation of the latter considerably deviates from a (low-field) sinusoidal dependence to a (high-field) step-like behaviour. These peculiar features are the consequences of the anisotropy of both (intra-band and inter-band ) scattering-amplitudes in the Brillouin zone when the relevant energy scales in the system - chemical potential, spin-orbit interaction and Zeeman energy - are all comparable to each other. The theory provides good qualitative agreement with experimental data in the literature.


title: Transparent Clouds: An Enhancement to Abstraction

abstract: With the introduction of various hardware/software technologies such as Cloud Technologies or Virtualization technologies, there has been a great potential to reuse ICT artifacts thanks to Abstraction and also Exchangeability features achieved via these technologies. These technologies also provide various advantages with respect to sustainability including resource consumption reduction (in the use phase only or in the whole life cycle). However, there is an additional but untapped potential associated with the anonymization of resources introduced by both abstraction and exchangeability features. By realizing on this potential, we can improve cloud solutions and reduce their by-product opacity, which usually prevents leveraging on the specialized but tweakable (i.e., nonessential modifications without changing the main function) features of components that are captured in the component models. This is especially a challenge in the case heterogeneous/disaggregated infrastructure where developing models to cover everything is practically impossible. In this work, by leveraging on the concept of pathways, we develop a few mechanisms that enable transparency and therefore tweakability of features even in the presence of abstraction and heterogeneity. In particular, the layered-stack approach to system decomposition is considered because of its role in both software defined networking (SDN) and Network Function Virtualization (NFV) system decompositions. For a concrete example, the case of dynamic frequency scaling of processors is considered and it is shown that the associated consumption could be considerably reduced without requiring additional changes to the middle components.


title: $PT$-symmetric graphene under a magnetic field

abstract: We propose a $PT$-symmetrically deformed version of the graphene tight-binding model under a magnetic field. We analyze the structure of the spectra and the eigenvectors of the Hamiltonians around the $K$ and $K'$ points, both in the $PT$-symmetric and $PT$-broken regions. In particular we show that the presence of the deformation parameter $V$ produces several interesting consequences, including the asymmetry of the zero-energy states of the Hamiltonians and the breakdown of the completeness of the eigenvector sets. We also discuss the biorthogonality of the eigenvectors, which {turns out to be} different in the $PT$-symmetric and $PT$-broken regions.


title: A Note on Moments of Limit Log Infinitely Divisible Stochastic Measures of Bacry and Muzy

abstract: A multiple integral representation of single and joint moments of the total mass of the limit log-infinitely divisible stochastic measure of Bacry and Muzy [$\textit{Comm. Math. Phys.}$ ${\bf 236}$: 449-475, 2003] is derived. The covariance structure of the total mass of the measure is shown to be logarithmic. A generalization of the Selberg integral corresponding to single moments of the limit measure is proposed and shown to satisfy a recurrence relation. The joint moments of the limit lognormal measure, classical Selberg integral with $\lambda_1=\lambda_2=0,$ and Morris integral are represented in the form of multiple binomial sums. For application, low moments of the limit log-Poisson measure are computed exactly and low joint moments of the limit lognormal measure are considered in detail.


title: (Almost) Closing the $\nu$MSM Sterile Neutrino Dark Matter Window with NuSTAR

abstract: We use NuSTAR observations of the Galactic Center to search for X-ray lines from the radiative decay of sterile neutrino dark matter. Finding no evidence of unknown lines, we set limits on the sterile neutrino mass and mixing angle. In most of the mass range 10-50 keV, these are now the strongest limits, at some masses improving upon previous limits by a factor of ~10. In the neutrino minimal standard model framework, where additional constraints from dark matter production and structure formation apply, the allowed parameter space is reduced by more than half. Future NuSTAR observations may be able to cover much of the remaining parameter space.


title: Structural Simplicity as a Restraint on the Structure of Amorphous Silicon

abstract: Understanding the structural origins of the properties of amorphous materials remains one of the most important challenges in structural science. In this study we demonstrate that local 'structural simplicity', embodied by the degree to which atomic environments within a material are similar to each other, is powerful concept for rationalising the structure of canonical amorphous material amorphous silicon (a-Si). We show, by restraining a reverse Monte Carlo refinement against pair distribution function (PDF) data to be simpler, that the simplest model consistent with the PDF is a continuous random network (CRN). A further effect of producing a simple model of a-Si is the generation of a (pseudo)gap in the electronic density of states, suggesting that structural homogeneity drives electronic homogeneity. That this method produces models of a-Si that approach the state-of-the-art without the need for chemically specific restraints (beyond the assumption of homogeneity) suggests that simplicity-based refinement approaches may allow experiment-driven structural modelling techniques to be developed for the wide variety of amorphous semiconductors with strong local order.


title: Bootstrapping a Five-Loop Amplitude Using Steinmann Relations

abstract: The analytic structure of scattering amplitudes is restricted by Steinmann relations, which enforce the vanishing of certain discontinuities of discontinuities. We show that these relations dramatically simplify the function space for the hexagon function bootstrap in planar maximally supersymmetric Yang-Mills theory. Armed with this simplification, along with the constraints of dual conformal symmetry and Regge exponentiation, we obtain the complete five-loop six-particle amplitude.


title: An EM based Iterative Method for Solving Large Sparse Linear Systems

abstract: We propose a novel iterative algorithm for solving a large sparse linear system. The method is based on the EM algorithm. If the system has a unique solution, the algorithm guarantees convergence with a geometric rate. Otherwise, convergence to a minimal Kullback--Leibler divergence point is guaranteed. The algorithm is easy to code and competitive with other iterative algorithms.


title: Structural Convergence Results for Approximation of Dominant Subspaces from Block Krylov Spaces

abstract: This paper is concerned with approximating the dominant left singular vector space of a real matrix $A$ of arbitrary dimension, from block Krylov spaces generated by the matrix $AA^T$ and the block vector $AX$. Two classes of results are presented. First are bounds on the distance, in the two and Frobenius norms, between the Krylov space and the target space. The distance is expressed in terms of principal angles. Second are quality of approximation bounds, relative to the best approximation in the Frobenius norm. For starting guesses $X$ of full column-rank, the bounds depend on the tangent of the principal angles between $X$ and the dominant right singular vector space of $A$. The results presented here form the structural foundation for the analysis of randomized Krylov space methods. The innovative feature is a combination of traditional Lanczos convergence analysis with optimal approximations via least squares problems.


title: The Inflation Technique for Causal Inference with Latent Variables

abstract: The problem of causal inference is to determine if a given probability distribution on observed variables is compatible with some causal structure. The difficult case is when the causal structure includes latent variables. We here introduce the $\textit{inflation technique}$ for tackling this problem. An inflation of a causal structure is a new causal structure that can contain multiple copies of each of the original variables, but where the ancestry of each copy mirrors that of the original. To every distribution of the observed variables that is compatible with the original causal structure, we assign a family of marginal distributions on certain subsets of the copies that are compatible with the inflated causal structure. It follows that compatibility constraints for the inflation can be translated into compatibility constraints for the original causal structure. Even if the constraints at the level of inflation are weak, such as observable statistical independences implied by disjoint causal ancestry, the translated constraints can be strong. We apply this method to derive new inequalities whose violation by a distribution witnesses that distribution's incompatibility with the causal structure (of which Bell inequalities and Pearl's instrumental inequality are prominent examples). We describe an algorithm for deriving all such inequalities for the original causal structure that follow from ancestral independences in the inflation. For three observed binary variables with pairwise common causes, it yields inequalities that are stronger in at least some aspects than those obtainable by existing methods. We also describe an algorithm that derives a weaker set of inequalities but is more efficient. Finally, we discuss which inflations are such that the inequalities one obtains from them remain valid even for quantum (and post-quantum) generalizations of the notion of a causal model.

categories: quant-ph math.ST stat.ME stat.ML stat.TH


title: Helical Kink Instability in a Confined Solar Eruption

abstract: A model for strongly writhing confined solar eruptions suggests an origin in the helical kink instability of a coronal flux rope which remains stable against the torus instability. This model is tested against the well observed filament eruption on 2002 May 27 in a parametric MHD simulation study which comprises all phases of the event. Good agreement with the essential observed properties is obtained. These include the confinement, terminal height, writhing, distortion, and dissolution of the filament, and the flare loops. The agreement is robust against variations in a representative range of parameter space. Careful comparisons with the observation data constrain the ratio of the external toroidal and poloidal field components to $B_\mathrm{et}/B_\mathrm{ep}\approx1$ and the initial flux rope twist to $\Phi\approx4\pi$. Different from ejective eruptions, two distinct phases of strong magnetic reconnection can occur. First, the erupting flux is cut by reconnection with overlying flux in the helical current sheet formed by the instability. If the resulting flux bundles are linked as a consequence of the erupting rope's strong writhing, they subsequently reconnect in the vertical current sheet between them. This reforms the overlying flux and a far less twisted flux rope, offering a pathway to homologous eruptions.


title: On the representation number of a crown graph

abstract: A graph $G=(V,E)$ is word-representable if there exists a word $w$ over the alphabet $V$ such that letters $x$ and $y$ alternate in $w$ if and only if $xy$ is an edge in $E$. It is known that any word-representable graph $G$ is $k$-word-representable for some $k$, that is, there exists a word $w$ representing $G$ such that each letter occurs exactly $k$ times in $w$. The minimum such $k$ is called $G$'s representation number. A crown graph $H_{n,n}$ is a graph obtained from the complete bipartite graph $K_{n,n}$ by removing a perfect matching. In this paper we show that for $n\geq 5$, $H_{n,n}$'s representation number is $\lceil n/2 \rceil$. This result not only provides a complete solution to the open Problem 7.4.2 in \cite{KL}, but also gives a negative answer to the question raised in Problem 7.2.7 in \cite{KL} on 3-word-representability of bipartite graphs. As a byproduct we obtain a new example of a graph class with a high representation number.


title: Limiting empirical distribution of zeros and critical points of random polynomials agree in general

abstract: In this article, we study critical points (zeros of derivative) of random polynomials. Take two deterministic sequences $\{a_n\}_{n\geq1}$ and $\{b_n\}_{n\geq1}$ of complex numbers whose limiting empirical measures are same. By choosing $\xi_n = a_n$ or $b_n$ with equal probability, define the sequence of polynomials by $P_n(z)=(z-\xi_1)\dots(z-\xi_n)$. We show that the limiting measure of zeros and critical points agree for this sequence of random polynomials under some assumption. We also prove a similar result for triangular array of numbers. A similar result for zeros of generalized derivative (can be thought as random rational function) is also proved. Pemantle and Rivin initiated the study of critical points of random polynomials. Kabluchko proved the result considering the zeros to be i.i.d. random variables.


title: A Search for the LHCb Charmed 'Pentaquark' using Photo-Production of $J/{\psi}$ at Threshold in Hall C at Jefferson Lab

abstract: We propose to measure the photo-production cross section of $J/{\psi}$ near threshold, in search of the recently observed LHCb hidden-charm resonances $P_c$(4380) and $P_c$(4450) consistent with 'pentaquarks'. The observation of these resonances in photo-production will provide strong evidence of the true resonance nature of the LHCb states, distinguishing them from kinematic enhancements. A bremsstrahlung photon beam produced with an 11 GeV electron beam at CEBAF covers the energy range of $J/{\psi}$ production from the threshold photo-production energy of 8.2 GeV, to an energy beyond the presumed $P_c$(4450) resonance. The experiment will be carried out in Hall C at Jefferson Lab using a 50{\mu}A electron beam incident on a 9% copper radiator. The resulting photon beam passes through a 15 cm liquid hydrogen target, producing $J/{\psi}$ mesons through a diffractive process in the $t$-channel, or through a resonant process in the $s$- and $u$-channel. The decay $e^+e^-$ pair of the $J/{\psi}$ will be detected in coincidence using the two high-momentum spectrometers of Hall C. The spectrometer settings have been optimized to distinguish the resonant $s$- and $u$-channel production from the diffractive $t$-channel $J/{\psi}$ production. The $s$- and $u$-channel production of the charmed 5-quark resonance dominates the $t$-distribution at large $t$. The momentum and angular resolution of the spectrometers is sufficient to observe a clear resonance enhancement in the total cross section and $t$-distribution. We request a total of 11 days of beam time with 9 days to carry the main experiment and 2 days to acquire the needed $t$-channel elastic $J/{\psi}$ production data for a calibration measurement. This calibration measurement in itself will greatly enhance our knowledge of $t$-channel elastic $J/{\psi}$ production near threshold.


title: Detection and Estimation of Multiple DoA Targets with Single Snapshot Measurements

abstract: In this paper, we explore the problems of detecting the number of narrow-band, far-field targets and estimating their corresponding directions of arrivals (DoAs) from single snapshot measurements. We use the principles of sparse signal recovery (SSR) for detection and estimation of multiple targets. In the SSR framework, the DoA estimation problem is grid based and can be posed as the lasso optimization problem. The corresponding DoA detection problem reduces to estimating the optimal regularization parameter ($\tau$) of the lasso problem for achieving the required probability of correct detection ($P_c$). We propose finite sample and asymptotic test statistics for detecting the number of sources with the required $P_c$ at moderate to high signal to noise ratios. Once the number of sources are detected, or equivalently the optimal $\hat{\tau}$ is estimated, the corresponding DoAs can be estimated by solving the lasso with regularization parameter set to $\hat{\tau}$.


title: Particle dynamics and Stochastic Resonance in Periodic potentials

abstract: We have studied the dynamics of a particle in a periodically driven underdamped periodic potential. Recent studies have reported the occurrence of Stochastic Resonance (SR) in such systems in the high frequency regime, using input energy per period of external drive as a quantifier. The particle trajectories in these systems can be in two dynamical states characterised by their definite energy and phase relation with the external drive. SR is due to the noise assisted transition of the particles between these two states. We study the role of damping on the occurrence of SR. We show that a driven underdamped periodic system exhibits SR only if the damping is below a particular limit. To explain this we study the syatem in the deterministic regime. The existence of the two dynamical states in the deterministic regime is dependent on the amount of damping and the amplitude od external drive. We also study the input energy distributions and phase difference of the response amplitude with the external drive as afunction of the friction parameter.


title: G64-12 and G64-37 are Carbon-Enhanced Metal-Poor Stars

abstract: We present new high-resolution chemical-abundance analyses for the well-known high proper-motion subdwarfs G64-12 and G64-37, based on very high signal-to-noise spectra (S/N ~ 700/1) with resolving power R ~ 95,000. These high-quality data enable the first reliable determination of the carbon abundances for these two stars; we classify them as carbon-enhanced metal-poor (CEMP) stars based on their carbonicities, which both exceed [C/Fe] = +1.0. They are sub-classified as CEMP- no Group-II stars, based on their location in the Yoon-Beers diagram of absolute carbon abundance, A(C) vs. [Fe/H], as well as on the conventional diagnostic [Ba/Fe]. The relatively low absolute carbon abundances of CEMP-no stars, in combination with the high effective temperatures of these two stars (Teff ~ 6500 K) weakens their CH molecular features to the point that accurate carbon abundances can only be estimated from spectra with very high S/N. A comparison of the observed abundance patterns with the predicted yields from massive, metal-free supernova models reduces the inferred progenitor masses by factors of ~ 2-3, and explosion energies by factors of ~ 10-15, compared to those derived using previously claimed carbon abundance estimates. There are certainly many more warm CEMP-no stars near the halo main-sequence turnoff that have been overlooked in past studies, directly impacting the derived frequencies of CEMP-no stars as a function of metallicity, a probe that provides important constraints on Galactic chemical evolution models, the initial mass function in the early Universe, and first-star nucleosynthesis.


title: Accurate De Novo Prediction of Protein Contact Map by Ultra-Deep Learning Model

abstract: Recently exciting progress has been made on protein contact prediction, but the predicted contacts for proteins without many sequence homologs is still of low quality and not very useful for de novo structure prediction. This paper presents a new deep learning method that predicts contacts by integrating both evolutionary coupling (EC) and sequence conservation information through an ultra-deep neural network formed by two deep residual networks. This deep neural network allows us to model very complex sequence-contact relationship as well as long-range inter-contact correlation. Our method greatly outperforms existing contact prediction methods and leads to much more accurate contact-assisted protein folding. Tested on three datasets of 579 proteins, the average top L long-range prediction accuracy obtained our method, the representative EC method CCMpred and the CASP11 winner MetaPSICOV is 0.47, 0.21 and 0.30, respectively; the average top L/10 long-range accuracy of our method, CCMpred and MetaPSICOV is 0.77, 0.47 and 0.59, respectively. Ab initio folding using our predicted contacts as restraints can yield correct folds (i.e., TMscore>0.6) for 203 test proteins, while that using MetaPSICOV- and CCMpred-predicted contacts can do so for only 79 and 62 proteins, respectively. Further, our contact-assisted models have much better quality than template-based models. Using our predicted contacts as restraints, we can (ab initio) fold 208 of the 398 membrane proteins with TMscore>0.5. By contrast, when the training proteins of our method are used as templates, homology modeling can only do so for 10 of them. One interesting finding is that even if we do not train our prediction models with any membrane proteins, our method works very well on membrane protein prediction. Finally, in recent blind CAMEO benchmark our method successfully folded 5 test proteins with a novel fold.


title: Fractal basin boundaries and escape dynamics in a multiwell potential

abstract: The escape dynamics in a two-dimensional multiwell potential is explored. A thorough numerical investigation is conducted in several types of two-dimensional planes and also in a three-dimensional subspace of the entire four-dimensional phase space in order to distinguish between non-escaping (ordered and chaotic) and escaping orbits. The determination of the location of the basins of escape towards the different escape channels and their correlations with the corresponding escape time of the orbits is undoubtedly an issue of paramount importance. It was found that in all examined cases regions of non-escaping motion coexist with several basins of escape. Furthermore, we monitor how the percentages of all types of orbits evolve when the total orbital energy varies. The larger escape periods have been measured for orbits with initial conditions in the fractal basin boundaries, while the lowest escape rates belong to orbits with initial conditions inside the basins of escape. The Newton-Raphson basins of attraction of the equilibrium points of the system have also been determined. We hope that our numerical analysis will be useful for a further understanding of the escape mechanism of orbits in open Hamiltonian systems with two degrees of freedom.


title: Unifying Markov Chain Approach for Disease and Rumor Spreading in Complex Networks

abstract: Spreading processes are ubiquitous in natural and artificial systems. They can be studied via a plethora of models, depending on the specific details of the phenomena under study. Disease contagion and rumor spreading are among the most important of these processes due to their practical relevance. However, despite the similarities between them, current models address both spreading dynamics separately. In this paper, we propose a general information spreading model that is based on discrete time Markov chains. The model includes all the transitions that are plausible for both a disease contagion process and rumor propagation. We show that our model not only covers the traditional spreading schemes, but that it also contains some features relevant in social dynamics, such as apathy, forgetting, and lost/recovering of interest. The model is evaluated analytically to obtain the spreading thresholds and the early time dynamical behavior for the contact and reactive processes in several scenarios. Comparison with Monte Carlo simulations shows that the Markov chain formalism is highly accurate while it excels in computational efficiency. We round off our work by showing how the proposed framework can be applied to the study of spreading processes occurring on social networks.


title: Pairwise, Magnitude, or Stars: What's the Best Way for Crowds to Rate?

abstract: We compare three popular techniques of rating content: the ubiquitous five star rating, the less used pairwise comparison, and the recently introduced (in crowdsourcing) magnitude estimation approach. Each system has specific advantages and disadvantages, in terms of required user effort, achievable user preference prediction accuracy and number of ratings required. We design an experiment where the three techniques are compared in an unbiased way. We collected 39'000 ratings on a popular crowdsourcing platform, allowing us to release a dataset that will be useful for many related studies on user rating techniques.


title: Attaining the quantum limit of passive imaging

abstract: We consider the problem, where a camera is tasked with determining one of two hypotheses: first with an incoherently-radiating quasi-monochromatic point source and the second with two identical closely spaced point sources. We are given that the total number of photons collected over an integration time is assumed to be the same under either hypothesis. For the one-source hypothesis, the source is taken to be on-axis along the line of sight and for the two-source hypothesis, we give ourselves the prior knowledge of the angular separation of the sources, and they are assumed to be identical and located symmetrically off-axis. This problem was studied by Helstrom in 1973, who evaluated the probability of error achievable using a sub-optimal optical measurement, with an unspecified structured realization. In this paper, we evaluate the quantum Chernoff bound, a lower bound on the minimum probability of error achievable by any physically-realizable receiver, which is exponentially tight in the regime that the integration time is high. We give an explicit structured receiver that separates three orthogonal spatial modes of the aperture field followed by quantum-noise-limited time-resolved photon measurement and show that this achieves the quantum Chernoff bound. In other words, the classical Chernoff bound of our mode-resolved detector exactly matches the quantum Chernoff bound for this problem. Finally, we evaluate the classical Chernoff bound on the error probability achievable using an ideal focal plane array---a signal shot-noise limited continuum photon-detection receiver with infinitely many infinitesimally-tiny pixels---and quantify its performance gap with the quantum limit.


title: Majorana quasiparticles of inhomogeneous Rashba chain

abstract: We investigate the inhomogeneous Rashba chain coupled to a superconducting substrate, hosting the Majorana quasiparticles near its edges. We discuss its subgap spectrum and study how robust are the zero-energy quasiparticles against the diagonal and off-diagonal disorder. Studying the $\mathbb{Z}_2$ topological invariant we show that disorder induced transition from the topologically non-trivial to trivial phases is manifested by characteristic features in the spatially-resolved quasiparticle spectrum at zero energy. We provide evidence for the non-local nature of the zero-energy Majorana quasiparticles, that are well preserved upon partitioning the chain into separate pieces. Even though the Majorana quasiparticles are not completely immune to inhomogeneity we show that they can spread onto other (normal) nanoscopic objects via the proximity effect.


title: Single photon in hierarchical architecture for physical reinforcement learning: Photon intelligence

abstract: Understanding and using natural processes for intelligent functionalities, referred to as natural intelligence, has recently attracted interest from a variety of fields, including post-silicon computing for artificial intelligence and decision making in the behavioural sciences. In a past study, we successfully used the wave-particle duality of single photons to solve the two-armed bandit problem, which constitutes the foundation of reinforcement learning and decision making. In this study, we propose and confirm a hierarchical architecture for single-photon-based reinforcement learning and decision making that verifies the scalability of the principle. Specifically, the four-armed bandit problem is solved given zero prior knowledge in a two-layer hierarchical architecture, where polarization is autonomously adapted in order to effect adequate decision making using single-photon measurements. In the hierarchical structure, the notion of layer-dependent decisions emerges. The optimal solutions in the coarse layer and in the fine layer, however, conflict with each other in some contradictive problems. We show that while what we call a tournament strategy resolves such contradictions, the probabilistic nature of single photons allows for the direct location of the optimal solution even for contradictive problems, hence manifesting the exploration ability of single photons. This study provides insights into photon intelligence in hierarchical architectures for future artificial intelligence as well as the potential of natural processes for intelligent functionalities.


title: An invariance principle for sums and record times of regularly varying stationary sequences

abstract: We prove a sequence of limiting results about weakly dependent stationary and regularly varying stochastic processes in discrete time. After deducing the limiting distribution for individual clusters of extremes, we present a new type of point process convergence theorem. It is designed to preserve the entire information about the temporal ordering of observations which is typically lost in the limit after time scaling. By going beyond the existing asymptotic theory, we are able to prove a new functional limit theorem. Its assumptions are satisfied by a wide class of applied time series models, for which standard limiting theory in the space $D$ of \cadlag\ functions does not apply. To describe the limit of partial sums in this more general setting, we use the space~$E$ of so--called decorated \cadlag\ functions. We also study the running maximum of partial sums for which a corresponding functional theorem can be still expressed in the familiar setting of space $D$. We further apply our method to analyze record times in a sequence of dependent stationary observations, even when their marginal distribution is not necessarily regularly varying. Under certain restrictions on dependence among the observations, we show that the record times after scaling converge to a relatively simple compound scale invariant Poisson process.


title: Boltzmann-Gibbs states in topological quantum walks and associated many-body systems: Fidelity and Uhlmann parallel transport analysis of Phase Transitions

abstract: We perform the fidelity analysis for Boltzmann-Gibbs-like states in order to investigate whether the topological order of 1D fermionic systems at zero temperature is maintained at finite temperatures. We use quantum walk protocols that are known to simulate topological phases and the respective quantum phase transitions for chiral symmetric Hamiltonians. Using the standard approaches of the fidelity analysis and the study of edge states, we conclude that no thermal-like phase transitions occur as temperature increases, i.e., the topological behaviour is washed out gradually. We also show that the behaviour of the Uhlmann geometric factor associated to the considered fidelity exhibits the same behaviour as the latter, thus confirming the results obtained using the previously established approaches.


title: Ensemble Learned Vaccination Uptake Prediction using Web Search Queries

abstract: We present a method that uses ensemble learning to combine clinical and web-mined time-series data in order to predict future vaccination uptake. The clinical data is official vaccination registries, and the web data is query frequencies collected from Google Trends. Experiments with official vaccine records show that our method predicts vaccination uptake effectively (4.7 Root Mean Squared Error). Whereas performance is best when combining clinical and web data, using solely web data yields comparative performance. To our knowledge, this is the first study to predict vaccination uptake using web data (with and without clinical data).


title: Integrability of the hyperbolic reduced Maxwell-Bloch equations for strongly correlated Bose-Einstein condensates

abstract: We derive and study the hyperbolic reduced Maxwell-Bloch equations (HRMB), a simplified model for the dynamics of strongly correlated Bose-Einstein condensates (BECs), and in particular for the interaction between the BEC atoms and its evaporated atoms under the strong interactions. This equation is one among four which are proven to be integrable via the existence of a Lax pair, and thus the method of inverse scattering transform. Another equation is the reduced Maxwell-Bloch equation of quantum optics and the two others do not have physical applications yet. By studying the linear stability of the constant solutions of these four equations we observe various regimes, from stable, to modulational unstable and unstable at all frequencies. The finite dimensional reduction of the RMB equations is also used to give more insight into the constant solutions of these equations. From this study, we find that the HRMB equation arising from strongly correlated BECS is stable under the particular condition that the transition rate of evaporation is not too large compared to the number of evaporated atoms. We then derive explicit soliton solutions of the RMB equations and use numerical simulations to show collisions of solitons and kink solitons.

categories: math-ph cond-mat.quant-gas math.MP nlin.SI


title: Multilevel Monte Carlo for Reliability Theory

abstract: As the size of engineered systems grows, problems in reliability theory can become computationally challenging, often due to the combinatorial growth in the cut sets. In this paper we demonstrate how Multilevel Monte Carlo (MLMC) - a simulation approach which is typically used for stochastic differential equation models - can be applied in reliability problems by carefully controlling the bias-variance tradeoff in approximating large system behaviour. In this first exposition of MLMC methods in reliability problems we address the canonical problem of estimating the expectation of a functional of system lifetime and show the computational advantages compared to classical Monte Carlo methods. The difference in computational complexity can be orders of magnitude for very large or complicated system structures.


title: Review of recent results on heavy-ion physics and astroparticle physics in ALICE at the LHC

abstract: In this work we present a summary of the most relevant results on heavy-ion and astroparticle physics in ALICE. The summary includes a brief overview of the current status on the characterization of the hot and dense QCD medium created in the heavy-ion collisions produced at the LHC, as well as the intriguing finding of collective-like phenomena in small collision systems.


title: Tunneling Conductance in Normal-Insulator-Superconductor junctions of Silicene

abstract: We theoretically investigate the transport properties of a normal-insulator-superconductor (NIS) junction of silicene in the thin barrier limit. Similar to graphene the tunneling conductance in such NIS structure exhibits an oscillatory behavior as a function of the strength of the barrier in the insulating region. However, unlike in graphene, the tunneling conductance in silicene can be controlled by an external electric field owing to its buckled structure. We also demonstrate the change in behavior of the tunneling conductance across the NIS junction as we change the chemical potential in the normal silicene region. In addition, at high doping levels in the normal region, the period of oscillation of the tunneling conductance as a function of the barrier strength changes from $\pi/2$ to $\pi$ with the variation of doping in the superconducting region of silicene.


title: Two computationally efficient polynomial-iteration infeasible interior-point algorithms for linear programming

abstract: Since the beginning of the development of interior-point methods, there exists a puzzling gap between the results in theory and the observations in numerical experience, i.e., algorithms with good polynomial bound are not computationally efficient and algorithms demonstrated efficiency in computation do not have a good or any polynomial bound. Todd raised a question in 2002: "Can we find a theoretically and practically efficient way to reoptimize?" This paper is an effort to close the gap. We propose two arc-search infeasible interior-point algorithms with infeasible central path neighborhood wider than all existing infeasible interior-point algorithms that are proved to be convergent. We show that the first algorithm is polynomial and its simplified version, if it terminates in finite iterations, has a complexity bound equal to the best known complexity bound for all (feasible or infeasible) interior-point algorithms. We demonstrate the computational efficiency of the proposed algorithms by testing all Netlib linear programming problems in standard form and comparing the numerical results to those obtained by Mehrotra's predictor-corrector algorithm and a recently developed more efficient arc-search algorithm (the convergence of these two algorithms is unknown). We conclude that the newly proposed algorithms are not only polynomial but also computationally competitive comparing to both Mehrotra's predictor-corrector algorithm and the efficient arc-search algorithm.


title: MRT-SS Calculator: An R Shiny Application for Sample Size Calculation in Micro-Randomized Trials

abstract: The micro-randomized trial (MRT) is a new experimental design which allows for the investigation of the proximal effects of a "just-in-time" treatment, often provided via a mobile device as part of a mobile health intervention. As with a traditional randomized controlled trial, computing the minimum required sample size to achieve a desired power is a crucial step in designing an MRT. We present MRT-SS Calculator, an online sample-size calculator for micro-randomized trials, built with R Shiny. MRT-SS Calculator requires specification of time-varying patterns for the proximal treatment effect and expected treatment availability. We illustrate the implementation of MRT-SS Calculator using a mobile health trial, HeartSteps. The application can be accessed from https://pengliao.shinyapps.io/mrt-calculator.


title: Adaptive Bayesian Spectral Analysis of Nonstationary Biomedical Time Series

abstract: Many studies of biomedical time series signals aim to measure the association between frequency-domain properties of time series and clinical and behavioral covariates. However, the time-varying dynamics of these associations are largely ignored due to a lack of methods that can assess the changing nature of the relationship through time. This article introduces a method for the simultaneous and automatic analysis of the association between the time-varying power spectrum and covariates. The procedure adaptively partitions the grid of time and covariate values into an unknown number of approximately stationary blocks and nonparametrically estimates local spectra within blocks through penalized splines. The approach is formulated in a fully Bayesian framework, in which the number and locations of partition points are random, and fit using reversible jump Markov chain Monte Carlo techniques. Estimation and inference averaged over the distribution of partitions allows for the accurate analysis of spectra with both smooth and abrupt changes. The proposed methodology is used to analyze the association between the time-varying spectrum of heart rate variability and self-reported sleep quality in a study of older adults serving as the primary caregiver for their ill spouse.


title: Deterministic and Stochastic Becker-D\"oring equations: Past and Recent Mathematical Developments

abstract: We present a survey on the results on a particular coagulation-fragmentation model given by the Becker-D\"oring equations. For both the deterministic and stochastic versions, we include well-posedness, long-time behavior, convergence rate towards equilibrium, coarsening and relation to transport equations, time-dependent properties, metastability and classical nucleation theory. All along this survey, we highlight recent results and open questions.


title: Photothermal effects in ultra-precisely stabilized tunable microcavities

abstract: We study the mechanical stability of a tunable high-finesse microcavity under ambient conditions and investigate light-induced effects that can both suppress and excite mechanical fluctuations. As an enabling step, we demonstrate the ultra-precise electronic stabilization of a microcavity. We then show that photothermal mirror expansion can provide high-bandwidth feedback and improve cavity stability by almost two orders of magnitude. At high intracavity power, we observe self-oscillations of mechanical resonances of the cavity. We explain the observations by a dynamic photothermal instability, leading to parametric driving of mechanical motion. For an optimized combination of electronic and photothermal stabilization, we achieve a feedback bandwidth of $500\,$kHz and a noise level of $1.1 \times 10^{-13}\,$m rms.


title: Approximate orthogonality of powers for ergodic affine unipotent diffeomorphisms on nilmanifolds

abstract: Let $ G $ be a connected, simply connected nilpotent Lie group and $ \Gamma < G $ a lattice. We prove that each ergodic diffeomorphism $ \phi(x\Gamma)=uA(x)\Gamma $ on the nilmanifold $ G/\Gamma $, where $ u\in G $ and $ A:G\to G $ is a unipotent automorphism satisfying $ A(\Gamma)=\Gamma $, enjoys the property of asymptotically orthogonal powers (AOP). Two consequences follow: (i) Sarnak's conjecture on M\"obius orthogonality holds in every uniquely ergodic model of an ergodic affine unipotent diffeomorphism; (ii) For ergodic affine unipotent diffeomorphisms themselves, the M\"obius orthogonality holds on so called typical short interval: $ \frac1 M\sum_{M\leq m<2M}\left|\frac1H\sum_{m\leq n<m+H} f(\phi^n(x\Gamma))\mu (n)\right|\to 0$ as $ H\to\infty $ and $ H/M\to0 $ for each $ x\Gamma\in G/\Gamma $ and each $ f\in C(G/\Gamma) $. In particular, the results in (i) and (ii) hold for ergodic nil-translations. Moreover, we prove that each nilsequence is orthogonal to the M\"obius function $\mu$ on a typical short interval. We also study the problem of lifting of the AOP property to induced actions and derive some applications on uniform distribution.


title: Towards the Quantization of Eddington-inspired-Born-Infeld Theory

abstract: The quantum effects close to the classical big rip singularity within the Eddington-inspired-Born-Infeld theory (EiBI) are investigated through quantum geometrodynamics. It is the first time that this approach is applied to a modified theory constructed upon Palatini formalism. The Wheeler-DeWitt (WDW) equation is obtained and solved based on an alternative action proposed in Ref.[1], under two different factor ordering choices. This action is dynamically equivalent to the original EiBI action while it is free of square root of the spacetime curvature. We consider a homogeneous, isotropic and spatially flat universe, which is assumed to be dominated by a phantom perfect fluid whose equation of state is a constant. We obtain exact solutions of the WDW equation based on some specific conditions. In more general cases, we propose a qualitative argument with the help of a Wentzel-Kramers-Brillouin (WKB) approximation to get further solutions. Besides, we also construct an effective WDW equation by simply promoting the classical Friedmann equations. We find that for all the approaches considered, the DeWitt condition hinting singularity avoidance is satisfied. Therefore the big rip singularity can be avoided through the quantum approach within the EiBI theory.


title: High-frequency, scaled MoS2 transistors

abstract: The interest in MoS2 for radio-frequency (RF) application has recently increased. However, little is known on the scaling behavior of transistors made from MoS2 for RF applications, which is important for establishing performance limits for electronic circuits based on 2D semiconductors on flexible and rigid substrates. Here, we present a systematic study of top-gated trilayer MoS2 RF transistors with gate lengths scaled down to 70 and 40 nm. In addition, by introducing edge-contacted injection of electrons in trilayer MoS2 devices, we decrease the contact resistance and as a result obtain the highest cutoff frequency of 6 GHz before the de-embedding procedure and 25 GHz after the de-embedding procedure.


title: Numerical solution of a semilinear parabolic degenerate Hamilton-Jacobi-Bellman equation with singularity

abstract: We consider a semilinear parabolic degenerated Hamilton-Jacobi-Bellman (HJB) equation with singularity which is related to a stochastic control problem with fuel constraint. The fuel constraint translates into a singular initial condition for the HJB equation. We first propose a transformation based on a change of variables that gives rise to an equivalent HJB equation with nonsingular initial condition but irregular coefficients. We then construct explicit and implicit numerical schemes for solving the transformed HJB equation and prove their convergences by establishing an extension to the result of Barles and Souganidis (1991).


title: Slanted snaking of localized Faraday waves

abstract: We report on an experimental, theoretical and numerical study of slanted snaking of spatially localized parametrically excited waves on the surface of a water-surfactant mixture in a Hele-Shaw cell. We demonstrate experimentally the presence of a hysteretic transition to spatially extended parametrically excited surface waves when the acceleration amplitude is varied, as well as the presence of spatially localized waves exhibiting slanted snaking. The latter extend outside the hysteresis loop. We attribute this behavior to the presence of a conserved quantity, the liquid volume, and introduce a universal model based on symmetry arguments, which couples the wave amplitude with such a conserved quantity. The model captures both the observed slanted snaking and the presence of localized waves outside the hysteresis loop, as demonstrated by numerical integration of the model equations.


title: A comparison of Redlich-Kister polynomial and cubic spline representations of the chemical potential in phase field computations

abstract: Free energies play a central role in many descriptions of equilibrium and non-equilibrium properties of solids. Continuum partial differential equations (PDEs) of atomic transport, phase transformations and mechanics often rely on first and second derivatives of a free energy function. The stability, accuracy and robustness of numerical methods to solve these PDEs are sensitive to the particular functional representations of the free energy. In this communication we investigate the influence of different representations of thermodynamic data on phase field computations of diffusion and two-phase reactions in the solid state. First-principles statistical mechanics methods were used to generate realistic free energy data for HCP titanium with interstitially dissolved oxygen. While Redlich-Kister polynomials have formed the mainstay of thermodynamic descriptions of multi-component solids, they require high order terms to fit oscillations in chemical potentials around phase transitions. Here we demonstrate that high fidelity fits to rapidly fluctuating free energy functions are obtained with spline functions. Spline functions that are many degrees lower than Redlich-Kister polynomials provide equal or superior fits to chemical potential data and, when used in phase field computations, result in solution times approaching an order of magnitude speed up relative to the use of Redlich-Kister polynomials.


title: Classification of the stable solution to the fractional $2<s<3$ Lane-Emden equation

abstract: We classify the stable solutions (positive or sign-changing, radial or not) to the following nonlocal Lane-Emden equation: $(-\Delta)^s u=|u|^{p-1}u$ in $\mathbb{R}^n$ for $2<s<3$.


title: Families from Supergroups and Predictions for Leptonic CP Violation

abstract: As was shown in 1984 by Caneschi, Farrar, and Schwimmer, decomposing representations of the supergroup SU(M|N), can give interesting anomaly-free sets of fermion representations of SU(M) x SU(N) x U(1). It is shown here that such groups can be used to construct realistic grand unified models with non-abelian gauged family symmetries. A particularly simple three-family example based on SU(5) x SU(2) x U(1) is studied. The forms of the mass matrices, including that of the right-handed neutrinos, are determined in terms of SU(2) Clebsch coefficients; and the model is able to fit the lepton sector and predict the Dirac CP-violating phase of the neutrinos.. Models of this type would have a rich phenomenology if part of the family symmetry is broken near the electroweak scale.


title: Magnetic properties of Proxima Centauri b analogues

abstract: The discovery of a planet around the closest star to our Sun, Proxima Centauri, represents a quantum leap in the testability of exoplanetary models. Unlike any other discovered exoplanet, models of Proxima b could be contrasted against near future telescopic observations and far future in-situ measurements. In this paper we aim at predicting the planetary radius and the magnetic properties (dynamo lifetime and magnetic dipole moment) of Proxima b analogues (solid planets with masses of $\sim 1-3\,M_\oplus$, rotation periods of several days and habitable conditions). For this purpose we build a grid of planetary models with a wide range of compositions and masses. For each point in the grid we run the planetary evolution model developed in \citet{Zuluaga2013}. Our model assumes small orbital eccentricity, negligible tidal heating and earth-like radiogenic mantle elements abundances. We devise a statistical methodology to estimate the posterior distribution of the desired planetary properties assuming simple prior distributions for the orbital inclination and bulk composition. Our model predicts that Proxima b would have a mass $1.3\leq M_{\rm p}\leq 2.3\,M_{\oplus}$ and a radius $R_{\rm p}=1.4^{+0.3}_{-0.2}\,R_{\oplus}$. In our simulations, most Proxima b analogues develop intrinsic dynamos that last for $\geq$4 Gyr (the estimated age of the host star). If alive, the dynamo of Proxima b have a dipole moment ${\cal M}_{\rm dip}>0.32^{\times 2.3}_{\div 2.9}{\cal M}_{\rm dip,\oplus}$. These results are not restricted to Proxima b but they also apply to earth-like planets having similar observed properties.


title: Operators and higher genus mirror curves

abstract: We perform further tests of the correspondence between spectral theory and topological strings, focusing on mirror curves of genus greater than one with nontrivial mass parameters. In particular, we analyze the geometry relevant to the SU(3) relativistic Toda lattice, and the resolved C^3/Z_6 orbifold. Furthermore, we give evidence that the correspondence holds for arbitrary values of the mass parameters, where the quantization problem leads to resonant states. We also explore the relation between this correspondence and cluster integrable systems.


title: Scientific Return of a Lunar Elevator

abstract: The concept of a space elevator dates back to Tsilokovsky, but they are not commonly considered in near-term plans for space exploration, perhaps because a terrestrial elevator would not be possible without considerable improvements in tether material. A Lunar Space Elevator (LSE), however, can be built with current technology using commercially available tether polymers. This paper considers missions leading to infrastructure capable of shortening the time, lowering the cost and enhancing the capabilities of robotic and human explorers. These missions use planetary scale tethers, strings many thousands of kilometers long stabilized either by rotation or by gravitational gradients. These systems promise major reduction in transport costs versus chemical rockets, in a rapid timeframe, for a modest investment. Science will thus benefit as well as commercial activities.


title: Flexible Bayesian Quantile Regression in Ordinal Models

abstract: The paper introduces an estimation method for flexible Bayesian quantile regression in ordinal (FBQROR) models i.e., an ordinal quantile regression where the error follows a generalized asymmetric Laplace (GAL) distribution. The GAL distribution, unlike the asymmetric Laplace (AL) distribution, allows to fix specific quantiles while simultaneously letting the mode, skewness and tails to vary. We also introduce the cumulative distribution function (necessary for constructing the likelihood) and the moment generating function of the GAL distribution. The algorithm is illustrated in multiple simulation studies and implemented to analyze public opinion on homeownership as the best long-term investment in the United States.


title: Generalized Spatial and Spatiotemporal Autoregressive Conditional Heteroscedasticity

abstract: In this paper, we introduce a new spatial model that incorporates heteroscedastic variance depending on neighboring locations. The proposed process is regarded as the spatial equivalent to the temporal autoregressive conditional heteroscedasticity (ARCH) model. We show additionally how the introduced spatial ARCH model can be used in spatiotemporal settings. In contrast to the temporal ARCH model, in which the distribution is known given the full information set of the prior periods, the distribution is not straightforward in the spatial and spatiotemporal setting. However, it is possible to estimate the parameters of the model using the maximum-likelihood approach. Via Monte Carlo simulations, we demonstrate the performance of the estimator for a specific spatial weighting matrix. Moreover, we combine the known spatial autoregressive model with the spatial ARCH model assuming heteroscedastic errors. Eventually, the proposed autoregressive process is illustrated using an empirical example. Specifically, we model lung cancer mortality in 3108 U.S. counties and compare the introduced model with two benchmark approaches.


title: Model structures on the category of complexes of quiver representations

abstract: In this paper, we study the category $C(Rep(\mathcal{Q}, \mathcal{A}))$ of complexes of representations of quiver $\mathcal{Q}$ with values in an abelian category $\mathcal{A}$. We develop a method for constructing some model structures on $C(Rep(\mathcal{Q}, \mathcal{A}))$ based on componentwise notion. Moreover, we also show that these model structures are monoidal. As an application of these model structures, we introduce some descriptions of the derived category of complexes of representations of $\mathcal{Q}$ in $\Mod R$.


title: Superdense coding over optical fiber links with complete Bell-state measurements

abstract: Adopting quantum communication to modern networking requires transmitting quantum information through fiber-based infrastructure. We report the first demonstration of superdense coding over optical fiber links, taking advantage of a complete Bell-state measurement enabled by time-polarization hyperentanglement, linear optics, and common single-photon detectors. We demonstrate the highest single-qubit channel capacity to date utilizing linear optics, $1.665 \pm 0.018$, and we provide a full experimental implementation of a hybrid, quantum-classical communication protocol for image transfer.


title: Analytical On-shell Calculation of Low Energy Higher Order Scattering

abstract: We demonstrate that the use of analytical on-shell methods involving calculation of the discontinuity across the t-channel cut associated with the exchange of a pair of massless particles (photons or gravitons) can be used to evaluate one-loop contributions to electromagnetic and gravitational scattering, with and without polarizability, reproducing via simple algebraic manipulations, results obtained previously, generally using Feynman diagram techniques. In the gravitational case the use of factorization permits a straightforward and algebraic calculation of higher order scattering without consideration of ghost contributions or of triple-graviton couplings, which made previous evaluations considerably more arduous.


title: Rarefied elliptic hypergeometric functions

abstract: Two exact evaluation formulae for multiple rarefied elliptic beta integrals related to the simplest lens space are proved. They generalize evaluations of the type I and II elliptic beta integrals attached to the root system $C_n$. In a special $n=1$ case, the simplest $p\to 0$ limit is shown to lead to a new class of $q$-hypergeometric identities. Symmetries of a rarefied elliptic analogue of the Euler-Gauss hypergeometric function are described and the respective generalization of the hypergeometric equation is constructed. Some extensions of the latter function to $C_n$ and $A_n$ root systems and corresponding symmetry transformations are considered. An application of the rarefied type II $C_n$ elliptic hypergeometric function to some eigenvalue problems is briefly discussed.


title: Lectures on Inflation

abstract: Planning to explore the beginning of the Universe? A lightweight introductory guide to the theory of Inflation.


title: Very Massive Tracers and Higher Derivative Biases

abstract: Most of the upcoming cosmological information will come from analyzing the clustering of the Large Scale Structures (LSS) of the universe through LSS or CMB observations. It is therefore essential to be able to understand their behavior with exquisite precision. The Effective Field Theory of Large Scale Structures (EFTofLSS) provides a consistent framework to make predictions for LSS observables in the mildly non-linear regime. In this paper we focus on biased tracers. We argue that in calculations at a given order in the dark matter perturbations, highly biased tracers will underperform because of their larger higher derivative biases. A natural prediction of the EFTofLSS is therefore that by simply adding higher derivative biases, all tracers should perform comparably well. We implement this prediction for the halo-halo and the halo-matter power spectra at one loop, and the halo-halo-halo, halo-halo-matter, and halo-matter-matter bispectra at tree-level, and compare with simulations. We find good agreement with the prediction: for all tracers, we are able to match the bispectra up to $k\simeq0.17\,h/$Mpc at $z=0$ and the power spectra to a higher wavenumber.


title: Convolutional Neural Networks for Text Categorization: Shallow Word-level vs. Deep Character-level

abstract: This paper reports the performances of shallow word-level convolutional neural networks (CNN), our earlier work (2015), on the eight datasets with relatively large training data that were used for testing the very deep character-level CNN in Conneau et al. (2016). Our findings are as follows. The shallow word-level CNNs achieve better error rates than the error rates reported in Conneau et al., though the results should be interpreted with some consideration due to the unique pre-processing of Conneau et al. The shallow word-level CNN uses more parameters and therefore requires more storage than the deep character-level CNN; however, the shallow word-level CNN computes much faster.


title: Peacock Bundles: Bundle Coloring for Graphs with Globality-Locality Trade-off

abstract: Bundling of graph edges (node-to-node connections) is a common technique to enhance visibility of overall trends in the edge structure of a large graph layout, and a large variety of bundling algorithms have been proposed. However, with strong bundling, it becomes hard to identify origins and destinations of individual edges. We propose a solution: we optimize edge coloring to differentiate bundled edges. We quantify strength of bundling in a flexible pairwise fashion between edges, and among bundled edges, we quantify how dissimilar their colors should be by dissimilarity of their origins and destinations. We solve the resulting nonlinear optimization, which is also interpretable as a novel dimensionality reduction task. In large graphs the necessary compromise is whether to differentiate colors sharply between locally occurring strongly bundled edges ("local bundles"), or also between the weakly bundled edges occurring globally over the graph ("global bundles"); we allow a user-set global-local tradeoff. We call the technique "peacock bundles". Experiments show the coloring clearly enhances comprehensibility of graph layouts with edge bundling.


title: Characterization of single layer anti-reflective coatings for bolometer-based rare event searches

abstract: Combining analysis from phonon signals and photon signals is a powerful technique for reducing backgrounds in bolometer-based rare event searches. Anti-reflective coatings can significantly increase the performance of the secondary light-sensing bolometer in these experiments. As a first step toward these improvements, coatings of SiO$_2$, HfO$_2$, and TiO$_2$ on Ge and Si wafers were fabricated and characterized at room temperature and multiple angles of incidence.


title: Motivic splitting principle

abstract: A variation on the splitting principle


title: Quantum walks and gravitational waves

abstract: A new family of discrete-time quantum walks (DTQWs) propagating on a regular $(1+2)$D spacetime lattice is introduced. The continuous limit of these DTQWs is shown to coincide with the dynamics of a Dirac fermion interacting with an arbitrary relativistic gravitational field. This family is used to model the influence of arbitrary linear gravitational waves (GWs) on DTQWs. Pure shear GWs are studied in detail. On large spatial scales, pure shear GWs do not modify the polarizations of the fermion eigen-modes, but rescale all energies by a common factor. On smaller scales typically comparable to two or three lattice steps, both polarizations and energies are modified in a non-trivial way. This effect is particularly salient on interference patterns between the fermion eigen-modes.


title: Superluminality, Black Holes and Effective Field Theory

abstract: Under the assumption that a UV theory does not display superluminal behavior, we ask what constraints on superluminality are satisfied in the effective field theory (EFT). We study two examples of effective theories: quantum electrodynamics (QED) coupled to gravity after the electron is integrated out, and the flat-space galileon. The first is realized in nature, the second is more speculative, but they both exhibit apparent superluminality around non-trivial backgrounds. In the QED case, we attempt, and fail, to find backgrounds for which the superluminal signal advance can be made larger than the putative resolving power of the EFT. In contrast, in the galileon case it is easy to find such backgrounds, indicating that if the UV completion of the galileon is (sub)luminal, quantum corrections must become important at distance scales of order the Vainshtein radius of the background configuration, much larger than the naive EFT strong coupling distance scale. Such corrections would be reminiscent of the non-perturbative Schwarzschild scale quantum effects that are expected to resolve the black hole information problem. Finally, a byproduct of our analysis is a calculation of how perturbative quantum effects alter charged Reissner-Nordstrom black holes.


title: QSL Squasher: A Fast Quasi-Separatrix Layer Map Calculator

abstract: Quasi-Separatrix Layers (QSLs) are a useful proxy for the locations where current sheets can develop in the solar corona, and give valuable information about the connectivity in complicated magnetic field configurations. However, calculating QSL maps even for 2-dimensional slices through 3-dimensional models of coronal magnetic fields is a non-trivial task as it usually involves tracing out millions of magnetic field lines with immense precision. Thus, extending QSL calculations to three dimensions has rarely been done until now. In order to address this challenge, we present QSL Squasher -- a public, open-source code, which is optimized for calculating QSL maps in both two and three dimensions on GPUs. The code achieves large processing speeds for three reasons, each of which results in an order-of-magnitude speed-up. 1) The code is parallelized using OpenCL. 2) The precision requirements for the QSL calculation are drastically reduced by using perturbation theory. 3) A new boundary detection criterion between quasi-connectivity domains is used, which quickly identifies possible QSL locations which need to be finely sampled by the code. That boundary detection criterion relies on finding the locations of abrupt field-line length changes, which we do by introducing a new Field-line Length Edge (FLEDGE) map. We find FLEDGE maps useful on their own as a quick-and-dirty substitute for QSL maps. QSL Squasher allows constructing high-resolution 3D FLEDGE maps in a matter of minutes, which is two orders of magnitude faster than calculating the corresponding 3D QSL maps. We include a sample of calculations done using QSL Squasher to demonstrate its capabilities as a QSL calculator, as well as to compare QSL and FLEDGE maps.


title: Equation of state effects and one-arm spiral instability in hypermassive neutron stars formed in eccentric neutron star mergers

abstract: We continue our investigations of the development and importance of the one-arm spiral instability in long-lived hypermassive neutron stars (HMNSs) formed in dynamical capture binary neutron star mergers. Employing hydrodynamic simulations in full general relativity, we find that the one-arm instability is generic in that it can develop in HMNSs within a few tens of milliseconds after merger for all equations of state in our survey. We find that mergers with stiffer equations of state tend to produce HMNSs with stronger $m=2$ azimuthal mode density deformations, and weaker $m=1$ components, relative to softer equations of state. We also find that for equations of state that can give rise to double-core HMNSs, large $m=1$ density modes can already be present due to asymmetries in the two cores. This results in the generation of $l=2$, $m=1$ gravitational wave modes even before the dominance of a one-arm mode that ultimately arises following merger of the two cores. Our results further suggest that stiffer equations of state give rise to HMNSs generating lower $m=1$ gravitational wave frequencies. Thus, if gravitational waves from the one-arm instability are detected, they could in principle constrain the neutron star equation of state. We estimate that, depending on the equation of state, the one-arm mode could potentially be detectable by second generation gravitational wave detectors at $\sim 10$ Mpc and by third generation ones at $\sim 100$ Mpc. Finally, we provide estimates of the properties of dynamical ejecta, as well as the accompanying kilonovae signatures.


title: Zodiacal Exoplanets in Time (ZEIT) IV: seven transiting planets in the Praesepe cluster

abstract: Open clusters and young stellar associations are attractive sites to search for planets and to test theories of planet formation, migration, and evolution. We present our search for, and characterization of, transiting planets in the ~800 Myr old Praesepe (Beehive, M44) Cluster from K2 light curves. We identify seven planet candidates, six of which we statistically validate to be real planets, the last of which requires more data. For each host star we obtain high-resolution NIR spectra to measure its projected rotational broadening and radial velocity, the latter of which we use to confirm cluster membership. We combine low-resolution spectra with the known cluster distance and metallicity to provide precise temperatures, masses, radii, and luminosities for the host stars. Combining our measurements of rotational broadening, rotation periods, and our derived stellar radii, we show that all planetary orbits are consistent with alignment to their host star's rotation. We fit the K2 light curves, including priors on stellar density to put constraints on the planetary eccentricities, all of which are consistent with zero. The difference between the number of planets found in Praesepe and Hyades (8 planets) and a similar dataset for Pleiades (0 planets, ~125 Myr) suggests a trend with age, but may be due to incompleteness of current search pipelines for younger, faster-rotating stars. We see increasing evidence that some planets continue to lose atmosphere past 800 Myr, as now two planets at this age have radii significantly larger than their older counterparts from Kepler.


title: No evidence for Population III stars or a Direct Collapse Black Hole in the z = 6.6 Lyman-$\alpha$ emitter 'CR7'

abstract: The z = 6.6 Lyman-$\alpha$ emitter 'CR7' has been claimed to have a Population III-like stellar population, or alternatively, be a candidate Direct Collapse Black Hole (DCBH). In this paper we investigate the evidence for these exotic scenarios using recently available, deeper, optical, near-infrared and mid-infrared imaging. We find strong Spitzer/IRAC detections for the main component of CR7 at 3.6 and 4.5 microns, and show that it has a blue colour ([3.6] - [4.5] $= -1.2\pm 0.3$). This colour cannot be reproduced by current Pop. III or pristine DCBH models. Instead, the results suggest that the [3.6] band is contaminated by the [OIII]4959,5007 emission line with an implied rest-frame equivalent width of EW_0 (H$\beta$ + [OIII]) $\gtrsim 2000$\AA. Furthermore, we find that new near-infrared data from the UltraVISTA survey supports a weaker He II 1640 emission line than previously measured, with EW_0 $= 40 \pm 30$\AA. For the fainter components of CR7 visible in Hubble Space Telescope imaging, we find no evidence that they are particularly red as previously claimed, and show that the derived masses and ages are considerably uncertain. In light of the likely detection of strong [OIII] emission in CR7 we discuss other more standard interpretations of the system that are consistent with the data. We find that a low-mass, narrow-line AGN can reproduce the observed features of CR7, including the lack of radio and X-ray detections. Alternatively, a young, low-metallicity (~1/200 solar) starburst, modelled including binary stellar pathways, can reproduce the inferred strength of the He II line and simultaneously the strength of the observed [OIII] emission, but only if the gas shows super-solar $\alpha$-element abundances (O/Fe ~ 5 O/Fe solar).


title: Gaia Parallax Zero Point From RR Lyrae Stars

abstract: Like Hipparcos, Gaia is designed to give absolute parallaxes, independent of any astrophysical reference system. And indeed, Gaia's internal zero-point error for parallaxes is likely to be smaller than any individual parallax error. Nevertheless, due in part to mechanical issues of unknown origin, there are many astrophysical questions for which the parallax zero-point error $\sigma(\pi_0)$ will be the fundamentally limiting constraint. These include the distance to the Large Magellanic Cloud and the Galactic Center. We show that by using the photometric parallax estimates for RR Lyrae stars (RRL) within 8kpc, via the ultra-precise infrared period-luminosity relation, one can independently determine a hyper-precise value for $\pi_{0}$. Despite their paucity relative to bright quasars, we show that RRL are competitive due to their order-of-magnitude improved parallax precision for each individual object relative to bright quasars. We show that this method is mathematically robust and well-approximated by analytic formulae over a wide range of relevant distances.


title: Minkowski flux vacua of type II supergravities

abstract: We study flux compactifications of 10d type II supergravities to 4d Minkowski space-time, supported by parallel orientifold Op-planes with 3 $\leq$ p $\leq$ 8. With some restrictions, the 4d Ricci scalar can be written as a negative sum of squares involving BPS-like conditions. Setting all squares to zero provides automatically a solution to 10d equations of motion. This way, we characterize a broad class, if not the complete set, of Minkowski flux vacua. We also conjecture an extension to include non-geometric fluxes. None of our results rely on supersymmetry.


title: Constrained simulations and excursion sets: understanding the risks and benefits of `genetically modified' haloes

abstract: Constrained realisations of Gaussian random fields are used in cosmology to design special initial conditions for numerical simulations. We review this approach and its application to density peaks providing several worked-out examples. We then critically discuss the recent proposal to use constrained realisations to modify the linear density field within and around the Lagrangian patches that form dark-matter haloes. The ambitious concept is to forge `genetically modified' haloes with some desired properties after the non-linear evolution. We demonstrate that the original implementation of this method is not exact but approximate because it tacitly assumes that protohaloes sample a set of random points with a fixed mean overdensity. We show that carrying out a full genetic modification is a formidable and daunting task requiring a mathematical understanding of what determines the biased locations of protohaloes in the linear density field. We discuss approximate solutions based on educated guesses regarding the nature of protohaloes. We illustrate how the excursion-set method can be adapted to predict the non-linear evolution of the modified patches and thus fine tune the constraints that are necessary to obtain preselected halo properties. This technique allows us to explore the freedom around the original algorithm for genetic modification. We find that the quantity which is most sensitive to changes is the halo mass-accretion rate at the mass scale on which the constraints are set. Finally we discuss constraints based on the protohalo angular momenta.


title: Photometry and spectroscopy of SN 2015bh in the galaxy NGC 2770

abstract: We present medium-resolution spectroscopy and multicolor photometry for the optical transient PSN J09093496+3307204 (SN 2015bh) in the galaxy NGC 2770, which has transferred into the supernova phase. The observations were carried out between February 2015 to May 2016. Both at the phase of the SN impostor (2015a) and at the supernova phase (2015b), besides Balmer emissions, the strong Fe II emissions are seen in the spectrum; so, these spectra resemble those of Williams Fe II type classical novae. The star is located near the edge of a dark nebula and notably absorbed (Av = 1.14 +/-0.15 mag). Taking into account this absorption, we determined maximum absolute magnitudes of Mv = -15.0 +/-0.3 mag at the 2015a phase and of Mv = -18.14 +/-0.30 mag at the 2015b phase. The light curve at the 2015b phase is similar to those of SN IIL. The supernova progenitor is a luminous blue variable (LBV) star with the powerful H{\alpha} emission. We considered several hypotheses of supernovae explosions following optical transients related with LBV. The hypothesis of core collapse of an evolved massive star interrupting the process of its merging with massive companion in a binary system (a failed luminous red nova) was chosen as the preferable one for this event.


title: 4D Scattering Amplitudes and Asymptotic Symmetries from 2D CFT

abstract: We reformulate the scattering amplitudes of 4D flat space gauge theory and gravity in the language of a 2D CFT on the celestial sphere. The resulting CFT structure exhibits an OPE constructed from 4D collinear singularities, as well as infinite-dimensional Kac-Moody and Virasoro algebras encoding the asymptotic symmetries of 4D flat space. We derive these results by recasting 4D dynamics in terms of a convenient foliation of flat space into 3D Euclidean AdS and Lorentzian dS geometries. Tree-level scattering amplitudes take the form of Witten diagrams for a continuum of (A)dS modes, which are in turn equivalent to CFT correlators via the (A)dS/CFT dictionary. The Ward identities for the 2D conserved currents are dual to 4D soft theorems, while the bulk-boundary propagators of massless (A)dS modes are superpositions of the leading and subleading Weinberg soft factors of gauge theory and gravity. In general, the massless (A)dS modes are 3D Chern-Simons gauge fields describing the soft, single helicity sectors of 4D gauge theory and gravity. Consistent with the topological nature of Chern-Simons theory, Aharonov-Bohm effects record the "tracks" of hard particles in the soft radiation, leading to a simple characterization of gauge and gravitational memories. Soft particle exchanges between hard processes define the Kac-Moody level and Virasoro central charge, which are thereby related to the 4D gauge coupling and gravitational strength in units of an infrared cutoff. Finally, we discuss a toy model for black hole horizons via a restriction to the Rindler region.


title: A Fixed-point Scheme for the Numerical Construction of Magnetohydrostatic Atmospheres in Three Dimensions

abstract: Magnetohydrostatic models of the solar atmosphere are often based on idealized analytic solutions because the underlying equations are too difficult to solve in full generality. Numerical approaches, too, are often limited in scope and have tended to focus on the two-dimensional problem. In this article we develop a numerical method for solving the nonlinear magnetohydrostatic equations in three dimensions. Our method is a fixed-point iteration scheme that extends the method of Grad and Rubin (Proc. 2nd Int. Conf. on Peaceful Uses of Atomic Energy 31, 190, 1958) to include a finite gravity force. We apply the method to a test case to demonstrate the method in general and our implementation in code in particular.


title: Role of gravity waves in vertical coupling during sudden stratospheric warmings

abstract: Gravity waves are primarily generated in the lower atmosphere, and can reach thermospheric heights in the course of their propagation. This paper reviews the recent progress in understanding the role of gravity waves in vertical coupling during sudden stratospheric warmings. Modeling of gravity wave effects is briefly reviewed, and the recent developments in the field are presented. Then, the impact of these waves on the general circulation of the upper atmosphere is outlined. Finally, the role of gravity waves in vertical coupling between the lower and the upper atmosphere is discussed in the context of sudden stratospheric warmings.


title: Complexity of quantum impurity problems

abstract: We give a quasi-polynomial time classical algorithm for estimating the ground state energy and for computing low energy states of quantum impurity models. Such models describe a bath of free fermions coupled to a small interacting subsystem called an impurity. The full system consists of $n$ fermionic modes and has a Hamiltonian $H=H_0+H_{imp}$, where $H_0$ is quadratic in creation-annihilation operators and $H_{imp}$ is an arbitrary Hamiltonian acting on a subset of $O(1)$ modes. We show that the ground energy of $H$ can be approximated with an additive error $2^{-b}$ in time $n^3 \exp{[O(b^3)]}$. Our algorithm also finds a low energy state that achieves this approximation. The low energy state is represented as a superposition of $\exp{[O(b^3)]}$ fermionic Gaussian states. To arrive at this result we prove several theorems concerning exact ground states of impurity models. In particular, we show that eigenvalues of the ground state covariance matrix decay exponentially with the exponent depending very mildly on the spectral gap of $H_0$. A key ingredient of our proof is Zolotarev's rational approximation to the $\sqrt{x}$ function. We anticipate that our algorithms may be used in hybrid quantum-classical simulations of strongly correlated materials based on dynamical mean field theory. We implemented a simplified practical version of our algorithm and benchmarked it using the single impurity Anderson model.

categories: quant-ph cond-mat.str-el math-ph math.MP


title: A robust covariance testing approach for high-throughput data

abstract: The problem of testing changes in covariance has received increasing attention in recent years, especially in the context of high-dimensional testing. A number of approaches have been proposed, all limited to the two-sample problem and involving varying statistics and assumptions on the number of features $p$ vs. the sample size $n$. There are no general approaches to test association of covariances with a continuous outcome. We propose a uniform framework for testing association of covariances with an experimental variable, whether discrete or continuous. The approach is not limited by the data dimensions. Our test procedure (i) does not rely on parametric assumptions, (ii) works well for a range of $p$ and $n$ (e.g., does not require $n > p$), (iii) provides correct type I error control, and (iv) includes four different statistics, to ensure power and flexibility under various settings, including a new "connectivity" statistic that is sensitive to changes in overall covariance magnitude. We demonstrate that, for the two-sample special case, the proposed statistics are permutationally equivalent or similar to existing proposed statistics. We demonstrate the power and utility of our approaches via simulation and analysis of real data. The approach is implemented in an $R$ package.


title: Space Development and Space Science Together, an Historic Opportunity

abstract: The national space programs have an historic opportunity to help solve the global-scale economic and environmental problems of Earth while becoming more effective at science through the use of space resources. Space programs will be more cost-effective when they work to establish a supply chain in space, mining and manufacturing then replicating the assets of the supply chain so it grows to larger capacity. This has become achievable because of advances in robotics and artificial intelligence. It is roughly estimated that developing a lunar outpost that relies upon and also develops the supply chain will cost about 1/3 or less of the existing annual budgets of the national space programs. It will require a sustained commitment of several decades to complete, during which time science and exploration become increasingly effective. At the end, this space industry will capable of addressing global-scale challenges including limited resources, clean energy, economic development, and preservation of the environment. Other potential solutions, including nuclear fusion and terrestrial renewable energy sources, do not address the root problem of our limited globe and there are real questions whether they will be inadequate or too late. While industry in space likewise cannot provide perfect assurance, it is uniquely able to solve the root problem, and it gives us an important chance that we should grasp. What makes this such an historic opportunity is that the space-based solution is obtainable as a side-benefit of doing space science and exploration within their existing budgets. Thinking pragmatically, it may take some time for policymakers to agree that setting up a complete supply chain is an achievable goal, so this paper describes a strategy of incremental progress.

categories: astro-ph.IM astro-ph.EP physics.pop-ph physics.soc-ph


title: Harder-Narasimhan theory for linear codes

abstract: In this text we develop some aspects of Harder-Narasimhan theory, slopes, semistability and canonical filtration, in the setting of combinatorial lattices. Of noticeable importance is the Harder-Narasimhan structure associated to a Galois connection between two lattices. It applies, in particular, to matroids. We then specialize this to linear codes. This could be done from at least three different approaches: using the sphere-packing analogy, or the geometric view, or the Galois connection construction just introduced. A remarkable fact is that these all lead to the same notion of semistability and canonical filtration. Relations to previous propositions towards a classification of codes, and to Wei's generalized Hamming weight hierarchy, are also discussed. Last, we study the important question of the preservation of semistability (or more generally the behaviour of slopes) under duality, and under tensor product. The former essentially follows from Wei's duality theorem for higher weights---and its matroid version---which we revisit in an appendix, developing analogues of the Riemann-Roch, Serre duality, Clifford, and gap and gonality sequence theorems for codes. Likewise the latter is closely related to the bound on higher weights of a tensor product, conjectured by Wei and Yang, and proved by Schaathun in the geometric language, which we reformulate directly in terms of codes. From this material we then derive semistability of tensor product.


title: Bypassing damaged nervous tissue

abstract: We show the principal ability of bypassing damaged demyelinated portions of nervous tissue, thereby restoring its normal function for the passage of action potentials. We carry out a theoretical analysis on the basis of the synchronization mechanism of action potential propagation along a bundle of neurons, proposed recently in [1]. And we discuss the feasibility of implement a bypass to restore damaged nervous tissue and creating an artificial neuron network.


title: Fully Automated Treatment Planning for Head and Neck Radiotherapy using a Voxel-Based Dose Prediction and Dose Mimicking Method

abstract: Recent works in automated radiotherapy treatment planning have used machine learning based on historical treatment plans to infer the spatial dose distribution for a novel patient directly from the planning image. We present an atlas-based approach which learns a dose prediction model for each patient (atlas) in a training database, and then learns to match novel patients to the most relevant atlases. The method creates a spatial dose objective, which specifies the desired dose-per-voxel, and therefore replaces any requirement for specifying dose-volume objectives for conveying the goals of treatment planning. A probabilistic dose distribution is inferred from the most relevant atlases, and is scalarized using a conditional random field to determine the most likely spatial distribution of dose to yield a specific dose prior (histogram) for relevant regions of interest. Voxel-based dose mimicking then converts the predicted dose distribution to a deliverable treatment plan dose distribution. In this study, we investigated automated planning for right-sided oropharaynx head and neck patients treated with IMRT and VMAT. We compare four versions of our dose prediction pipeline using a database of 54 training and 12 independent testing patients. Our preliminary results are promising; automated planning achieved a higher number of dose evaluation criteria in 7 patients and an equal number in 4 patients compared with clinical. Overall, the relative number of criteria achieved was higher for automated planning versus clinical (17 vs 8) and automated planning demonstrated increased sparing for organs at risk (52 vs 44) and better target coverage/uniformity (41 vs 31).


title: A Conic Section Approach to the Relativistic Reflection Law

abstract: We consider the reflection of light, from a stationary source, off of a uniformly moving flat mirror, and derive the relativistic reflection law using well-known properties of conic sections. The effective surface of reflection (ESR) is defined as the loci of intersection of all beams, emanating from the source at a given time, with the moving mirror. Fermat principle of least time is then applied to ESR and it is shown that, assuming the independence of speed of light, the result is identical with the relativistic reflection law. For a uniformly moving mirror ESR is a conic and the reflection law becomes a case of bi-angular equation of the conic, with the incident and reflected beams coinciding with the focal rays of the conic. A short calculus-based proof for accelerating mirrors is also given.


title: Influence of coupling on thermal forces and dynamic friction in plasmas with multiple ion species

abstract: The recently proposed effective potential theory [Phys. Rev. Lett. 110, 235001 (2013)] is used to investigate the influence of coupling on inter-ion-species diffusion and momentum exchange in multi-component plasmas. Thermo-diffusion and the thermal force are found to diminish rapidly as strong coupling onsets. For the same coupling parameters, the dynamic friction coefficient is found to tend to unity. These results provide an impetus for addressing the role of coupling on diffusive processes in inertial confinement fusion experiments.


title: Extended self-similarity in moment-generating-functions in wall-bounded turbulence at high Reynolds number

abstract: In wall-bounded turbulence, the moment generating functions (MGFs) of the streamwise velocity fluctuations $\left<\exp(qu_z^+)\right>$ develop power-law scaling as a function of the wall normal distance $z/\delta$. Here $u$ is the streamwise velocity fluctuation, $+$ indicates normalization in wall units (averaged friction velocity), $z$ is the distance from the wall, $q$ is an independent variable and $\delta$ is the boundary layer thickness. Previous work has shown that this power-law scaling exists in the log-region {\small $3Re_\tau^{0.5}\lesssim z^+$, $z\lesssim 0.15\delta$}, where $Re_\tau$ is the friction velocity-based Reynolds numbers. Here we present empirical evidence that this self-similar scaling can be extended, including bulk and viscosity-affected regions $30<z^+$, $z<\delta$, provided the data are interpreted with the Extended-Self-Similarity (ESS), i.e. self-scaling of the MGFs as a function of one reference value, $q_o$. ESS also improves the scaling properties, leading to more precise measurements of the scaling exponents. The analysis is based on hot-wire measurements from boundary layers at $Re_\tau$ ranging from $2700$ to $13000$ from the Melbourne High-Reynolds-Number-Turbulent-Boundary-Layer-Wind-Tunnel. Furthermore, we investigate the scalings of the filtered, large-scale velocity fluctuations $u^L_z$ and of the remaining small-scale component, $u^S_z=u_z-u^L_z$. The scaling of $u^L_z$ falls within the conventionally defined log region and depends on a scale that is proportional to {\small $l^+\sim Re_\tau^{1/2}$}; the scaling of $u^{S}_z$ extends over a much wider range from $z^+\approx 30$ to $z\approx 0.5\delta$. Last, we present a theoretical construction of two multiplicative processes for $u^L_z$ and $u^S_z$ that reproduce the empirical findings concerning the scalings properties as functions of $z^+$ and in the ESS sense.


title: Which subsets of an infinite random graph look random?

abstract: Given a countable graph, we say a set $A$ of its vertices is \emph{universal} if it contains every countable graph as an induced subgraph, and $A$ is \emph{weakly universal} if it contains every finite graph as an induced subgraph. We show that, for almost every graph on $\mathbb N$, $(1)$ every set of positive upper density is universal, and $(2)$ every set with divergent reciprocal sums is weakly universal. We show that the second result is sharp (i.e., a random graph on $\mathbb N$ will almost surely contain non-universal sets with divergent reciprocal sums) and, more generally, that neither of these two results holds for a large class of partition regular families.