title: A probability theory for non-equilibrium gravitational systems

abstract: This paper uses dynamical invariants to describe the evolution of collisionless systems subject to time-dependent gravitational forces without resorting to maximum-entropy probabilities. We show that collisionless relaxation can be viewed as a special type of diffusion process in the integral-of-motion space. In time-varying potentials with a fixed spatial symmetry the diffusion coefficients are closely related to virial quantities, such as the specific moment of inertia, the virial factor and the mean kinetic and potential energy of microcanonical particle ensembles. The non-equilibrium distribution function (DF) is found by convolving the initial DF with the Green function that solves Einstein's equation for freely diffusing particles. Such a convolution also yields a natural solution to the Fokker-Planck equations in the energy space. Our mathematical formalism can be generalized to potentials with a time-varying symmetry, where diffusion extends over multiple dimensions of the integral-of-motion space. The new probability theory is in many ways analogous to stochastic calculus, with two significant differences: (i) the equations of motion that govern the trajectories of particles are fully deterministic, and (ii) the diffusion coefficients can be derived self-consistently from microcanonical phase-space averages without relying on ergodicity assumptions. For illustration we follow the cold collapse of $N$-body models in a time-dependent logarithmic potential. Comparison between the analytical and numerical results shows excellent agreement in regions where the potential evolution does not depart too strongly from the adiabatic regime.


title: $1/f^\alpha$ noise and generalized diffusion in random Heisenberg spin systems

abstract: We study the `flux noise' spectrum of random-bond quantum Heisenberg spin systems using a real-space renormalization group (RSRG) procedure that accounts for both the renormalization of the system Hamiltonian and of a generic probe that measures the noise. For spin chains, we find that the dynamical structure factor $S_q(f)$, at finite wave-vector $q$, exhibits a power-law behavior both at high and low frequencies $f$, with exponents that are connected to one another and to an anomalous dynamical exponent through relations that differ at $T = 0$ and $T = \infty$. The low-frequency power-law behavior of the structure factor is inherited by any generic probe with a finite band-width and is of the form $1/f^\alpha$ with $0.5 < \alpha < 1$. An analytical calculation of the structure factor, assuming a limiting distribution of the RG flow parameters (spin size, length, bond strength) confirms numerical findings. More generally, we demonstrate that this form of the structure factor, at high temperatures, is a manifestation of anomalous diffusion which directly follows from a generalized spin-diffusion propagator. We also argue that $1/f$-noise is intimately connected to many-body-localization at finite temperatures. In two dimensions, the RG procedure is less reliable; however, it becomes convergent for quasi-one-dimensional geometries where we find that one-dimensional $1/f^\alpha$ behavior is recovered at low frequencies; the latter configurations are likely representative of paramagnetic spin networks that produce $1/f^\alpha$ noise in SQUIDs.


title: A possible geometrical origin of the accelerated expansion of the universe

abstract: The modified geodetic brane cosmology (MGBC) is tested with observational data. The MGBC is derived from the geodetic brane gravity action corrected by the extrinsic curvature of the braneworld. The density parameter coming from this additional term produces an accelerated expansion of geometrical origin. Subject to the Supernovae Ia, Observable Hubble parameter, Baryon Acoustic Oscillations and Cosmic Microwave Background probes, the obtained fit provides enough evidence in the sense that the extrinsic curvature effect is able to reproduce the accelerated expansion of the universe without need of invoking dark energy, exotic matter or cosmological constant. Moreover the MGBC is free of the problems present in other braneworld models.


title: Real time evolution of non-Gaussian cumulants in the QCD critical regime

abstract: We derive a coupled set of equations that describe the non-equilibrium evolution of cumulants of critical fluctuations for space-time trajectories on the cross-over side of the QCD phase diagram. In particular, novel expressions are obtained for the non-equilibrium evolution of non-Gaussian Skewness and Kurtosis cumulants. Utilizing a simple model of the space-time evolution of a heavy- ion collision, we demonstrate that, depending on the relaxation rate of critical fluctuations, Skewness and Kurtosis can differ significantly in magnitude as well as in sign from equilibrium expectations. Memory effects are important and shown to persist even for trajectories that skirt the edge of the critical regime. We use phenomenologically motivated parameterizations of freeze-out curves, and of the beam energy dependence of the net baryon chemical potential, to explore the implications of our model study for the critical point search in heavy-ion collisions.


title: On the spatial distributions of stars and gas in numerical simulations of molecular clouds

abstract: We compare the spatial distribution of stars which form in hydrodynamical simulations to the spatial distribution of the gas, using the $\mathcal{Q}$-parameter. The $\mathcal{Q}$-parameter enables a self-consistent comparison between the stars and gas because it uses a pixelated image of the gas as a distribution of points, in the same way that the stars (sink particles in the simulations) are a distribution of points. We find that, whereas the stars have a substructured, or hierarchical spatial distribution ($\mathcal{Q} \sim 0.4 - 0.7$), the gas is dominated by a smooth, concentrated component and typically has $\mathcal{Q} \sim 0.9$. We also find no statistical difference between the structure of the gas in simulations that form with feedback, and those that form without, despite these two processes producing visually different distributions. These results suggest that the link between the spatial distributions of gas, and the stars which form from them, is non-trivial.


title: Perturbative Renormalization of Neutron-Antineutron Operators

abstract: Two-loop anomalous dimensions and one-loop renormalization scheme matching factors are calculated for six-quark operators responsible for neutron-antineutron transitions. When combined with lattice QCD determinations of the matrix elements of these operators, our results can be used to reliably predict the neutron-antineutron vacuum transition time, $\tau_{n\bar{n}}$, in terms of basic parameters of baryon-number violating beyond-the-Standard-Model theories. The operators are classified by their chiral transformation properties, and a basis in which there is no operator mixing due to strong interactions is identified. Operator projectors that are required for non-perturbative renormalization of the corresponding lattice QCD six-quark operator matrix elements are constructed. A complete calculation of $\delta m = 1/\tau_{n\bar{n}}$ in a particular beyond-the-Standard-Model theory is presented as an example to demonstrate how operator renormalization and results from lattice QCD are combined with experimental bounds on $\delta m$ to constrain the scale of new baryon-number violating physics. At the present computationally accessible lattice QCD matching scale of $\sim$ 2 GeV, the next-to-next-to-leading-order effects calculated in this work correct the leading-order plus next-to-leading-order $\delta m$ predictions of beyond-the-Standard-Model theories by $< 26\%$. Next-to-next-to-next-to-leading-order effects provide additional unknown corrections to predictions of $\delta m$ that are estimated to be $< 7\%$.


title: Stellar Masses and SFRs for 1M Galaxies from SDSS and WISE

abstract: We combine SDSS and WISE photometry for the full SDSS spectroscopic galaxy sample, creating SEDs that cover lambda=0.4-22 micron for an unprecedented large and comprehensive sample of 858,365 present-epoch galaxies. Using MAGPHYS we then model simultaneously and consistently both the attenuated stellar SED and the dust emission at 12 micron and 22 micron, producing robust new calibrations for monochromatic mid-IR star formation rate proxies. These modeling results provide the first mid-IR-based view of the bi-modality in star formation activity among galaxies, exhibiting the sequence of star-forming galaxies (main sequence) with a slope of dlogSFR/dlogM*=0.80 and a scatter of 0.39 dex. We find that these new star-formation rates along the SF main sequence are systematically lower by a factor of 1.4 than those derived from optical spectroscopy. We show that for most present-day galaxies the 0.4-22 micron SED fits can exquisitely predict the fluxes measured by Herschel at much longer wavelengths. Our analysis also illustrates that the majorities of stars in the present-day universe is formed in luminous galaxies (~L*) in and around the green valley of the color-luminosity plane. We make the matched photometry catalog and SED modeling results publicly available.


title: Power-law Optical Conductivity from Unparticles: Application to the Cuprates

abstract: We calculate the optical conductivity using several models for unparticle or scale-invariant matter. Within a Gaussian action for unparticles that is gauged with Wilson lines, we find that the conductivity computed from the Kubo formalism with vertex corrections yields no non-trivial deviation from the free-theory result. This result obtains because at the Gaussian level, unparticles are just a superposition of particle fields and hence any transport property must be consistent with free theory. Beyond the Gaussian approach, we adopt the continuous mass formulation of unparticles and calculate the Drude conductivity directly. We show that unparticles in this context can be tailored to yield an algebraic conductivity that scales as $\omega^{-2/3}$ with the associated phase angle between the imaginary and real parts of $\arctan\frac{\sigma_2}{\sigma_1}=60^\circ$ as is seen in the cuprates. Given the recent results\cite{Donos2014,Rangamani2015,Langley} that gravitational crystals lack a power-law optical conductivity, this constitutes the first consistent account of the $\omega^{-2/3}$ conductivity and the phase angle seen in optimally doped cuprates. Our results indicate that at each frequency in the scaling regime, excitations on all energy scales contribute. Hence, incoherence is at the heart of the power-law in the optical conductivity in strongly correlated systems such as the cuprates.


title: The anomalous Floquet-Anderson insulator as a non-adiabatic quantized charge pump

abstract: Periodically driven quantum systems provide a novel and versatile platform for realizing topological phenomena. Among these are analogs of topological insulators and superconductors, attainable in static systems; however, some of these phenomena are unique to the periodically driven case. Here, we show that disordered, periodically driven systems admit an "anomalous" two dimensional phase, whose quasi-energy spectrum consists of chiral edge modes that coexist with a fully localized bulk - an impossibility for static Hamiltonians. This unique situation serves as the basis for a new topologically-protected non-equilibrium transport phenomenon: quantized non-adiabatic charge pumping. We identify the bulk topological invariant that characterizes the new phase (which we call the "anomalous Floquet Anderson Insulator", or AFAI). We provide explicit models which constitute a proof of principle for the existence of the new phase. Finally, we present evidence that the disorder-driven transition from the AFAI to a trivial, fully localized phase is in the same universality class as the quantum Hall plateau transition.


title: Flavor Changing Heavy Higgs Interactions at the LHC

abstract: A general two Higgs doublet model (2HDM) is adopted to study the signature of flavor changing neutral Higgs (FCNH) decay $\phi^0 \to t\bar{c}+\bar{t}c$, where $\phi^0$ could be a CP-even scalar ($H^0$) or a CP-odd pseudoscalar ($A^0$). Measurement of the light 125 GeV neutral Higgs boson ($h^0$) couplings at the Large Hadron Collider (LHC) favor the decoupling limit or the alignment limit of a 2HDM, in which gauge boson and diagonal fermion couplings of $h^0$ approach Standard Model values. In such limit, FCNH couplings of $h^0$ are naturally suppressed by a small mixing parameter $\cos(\beta-\alpha)$, while the off-diagonal couplings of heavier neutral scalars $\phi^0$ are sustained by $\sin(\beta-\alpha) \sim 1$. We study physics background from dominant processes with realistic acceptance cuts and tagging efficiencies. Promising results are found for the LHC running at 13 or 14 TeV collision energies.


title: The effect of stellar feedback on a Milky Way-like galaxy and its gaseous halo

abstract: We present the study of a set of N-body+SPH simulations of a Milky Way-like system produced by the radiative cooling of hot gas embedded in a dark matter halo. The galaxy and its gaseous halo evolve for 10 Gyr in isolation, which allows us to study how internal processes affect the evolution of the system. We show how the morphology, the kinematics and the evolution of the galaxy are affected by the input supernova feedback energy E$_{\rm SN}$, and we compare its properties with those of the Milky Way. Different values of E$_{\rm SN}$ do not significantly affect the star formation history of the system, but the disc of cold gas gets thicker and more turbulent as feedback increases. Our main result is that, for the highest value of E$_{\rm SN}$ considered, the galaxy shows a prominent layer of extra-planar cold (log(T)<4.3) gas extended up to a few kpc above the disc at column densities of $10^{19}$ cm$^{-2}$. The kinematics of this material is in agreement with that inferred for the HI halos of our Galaxy and NGC 891, although its mass is lower. Also, the location, the kinematics and the typical column densities of the hot (5.3<log(T)<5.7) gas are in good agreement with those determined from the O$_{\rm VI}$ absorption systems in the halo of the Milky Way and external galaxies. In contrast with the observations, however, gas at log(T)<5.3 is lacking in the circumgalactic region of our systems.


title: Opening up the compressed region of stop searches at 13 TeV LHC

abstract: Light top superpartners play a key role in stabilizing the electroweak scale in supersymmetric theories. For $R$-parity conserved supersymmetric models, traditional searches are not sensitive to the compressed regions. In this paper, we propose a new method targeting this region, with stop and neutralino mass splitting ranging from $m_{\tilde t} - m_\chi \gtrsim m_t$ to about 20 GeV. In particular, we focus on the signal process in which a pair of stops are produced in association with a hard jet, and define a new observable $R_M$ whose distribution has a peak in this compressed region. The position of the peak is closely correlated with $m_{\tilde t}$. We show that for 13 TeV LHC with a luminosity of 3000 fb$^{-1}$, this analysis can extend the reach of stop in the compressed region to $m_{\tilde t}$ around 800 GeV.


title: Jordan Form and Quantum Tomography

abstract: In this brief article we indicate a connection between Jordan normal form of a square matrix and the stroboscopic approach to quantum tomography. We show that the index of cyclicy of a generator of evolution, which receives much attention in the stroboscopic tomography, can be defined and computed by referring to the Jordan decomposition of a square matrix. The result presented in this work shows a relation between terminology from quantum tomography and linear algebra.


title: A matter of measurement: rotation velocities and the velocity function of dwarf galaxies

abstract: The velocity function derived from large scale surveys can be compared with the predictions of LCDM cosmology, by matching the measured rotation velocities Vrot of galaxies to the maximum circular velocity of dark matter (DM) halos Vmax. For Vrot<50km/s, a major discrepancy arises between the observed and LCDM velocity functions. However, the manner in which different observational measures of Vrot are associated with Vmax is not straight forward in dwarf galaxies. We instead relate galaxies to DM halos using the empirical baryon- mass to halo-mass relation, and show that different observational measures of Vrot result in very different velocity functions. We show how the W50 velocity function, i.e. using the HI profile line width at 50% of peak HI flux to measure Vrot, can be reconciled with a LCDM cosmology. Our semi-empirical methodology allows us to determine the region of rotation curves that are probed by HI measurements (RHI), and shows that the Vrot of dwarfs are generally measured at a fraction of Rmax, explaining their tendency to have rising rotation curves. We provide fitting formulae for relating RHI and Reff (the effective radius) to the virial radius of DM halos. To continue to use velocity functions as a probe of LCDM cosmology, it is necessary to be precise about how the different measures of rotation velocity are probing the mass of the DM halos, dropping the assumption that any measure of rotational velocity can be equally used as a proxy for Vmax.


title: Starbursts are preferentially interacting: confirmation from the nearest galaxies

abstract: We complement a recent ApJ Letter by Luo et al. by comparing the fraction of starburst galaxies which are interacting with the overall fraction of interacting galaxies in the nearby galaxy population (within 40 Mpc). We confirm that in starburst galaxies the fraction of interacting galaxies is enhanced, by a factor of around 2, but crucially we do so by studying a sample of almost 1500 of the nearest galaxies, including many dwarfs and irregulars. We discuss how adjusting the starburst definition influences the final result and conclude that our result is stable. We find significantly lower fractions of interacting galaxies than Luo et al. did from their larger but more distant sample of galaxies, and argue that the difference is most likely due to various biases in the sample selection, with a representative sample of the nearest galaxies, such as the one used here, being the best possible representation of a general picture. Our overall conclusion is that interactions can and do increase the number of starburst galaxies, and that this result is extremely robust. By far most starburst galaxies, however, show no evidence of a present interaction.


title: Supergravity divergences, supersymmetry and automorphic forms

abstract: We investigate the occurrence of divergences in maximal supergravity in various dimensions from the point of view of supersymmetry constraints on the U-duality invariant threshold functions defining the higher derivative couplings in the effective action. Our method makes use of tensorial differential equations that constrain these couplings. We study in detail the Fourier modes and wave-front sets of the associated automorphic functions and find that they are always associated to special nilpotent orbits.


title: An online parameter identification method for time dependent partial differential equations

abstract: Online parameter identification is of importance, e.g., for model predictive control. Since the parameters have to be identified simultaneously to the process of the modeled system, dynamical update laws are used for state and parameter estimates. Most of the existing methods for infinite dimensional systems either impose strong assumptions on the model or cannot handle partial observations. Therefore we propose and analyze an online parameter identification method that is less restrictive concerning the underlying model and allows for partial observations and noisy data. The performance of our approach is illustrated by some numerical experiments.


title: Calculating the 3D magnetic field of ITER for European TBM studies

abstract: The magnetic perturbation due to the ferromagnetic test blanket modules (TBMs) may deteriorate fast ion confinement in ITER. This effect must be quantified by numerical studies in 3D. We have implemented a combined finite element method (FEM) -- Biot-Savart law integrator method (BSLIM) to calculate the ITER 3D magnetic field and vector potential in detail. Unavoidable geometry simplifications changed the mass of the TBMs and ferritic inserts (FIs) up to 26%. This has been compensated for by modifying the nonlinear ferromagnetic material properties accordingly. Despite the simplifications, the computation geometry and the calculated fields are highly detailed. The combination of careful FEM mesh design and using BSLIM enables the use of the fields unsmoothed for particle orbit-following simulations. The magnetic field was found to agree with earlier calculations and revealed finer details. The vector potential is intended to serve as input for plasma shielding calculations.


title: Moving Beyond Effective Field Theory in Dark Matter Searches at Colliders

abstract: In the past few years, the interest to collider searches for direct dark matter (DM) production has been growing exponentially. A variety of "Mono-X" signatures have been considered, where X stands for a probe particle recoiling against DM particles, which allows for the event to be triggerable. So far, the analysis of these signatures has been largely carried out in the framework of effective field theory (EFT), which allows for a comparison of the collider searches with searches in direct detection experiments. Unfortunately, as it has been recently pointed out by a number of authors, the EFT approach has severe limitations and may result in drastically underestimated or overestimated reach. I'll discuss these limitations and the new ideas in interpreting the collider searches for DM.


title: Classification of Spatio-Temporal Data via Asynchronous Sparse Sampling: Application to Flow Around a Cylinder

abstract: We present a novel method for the classification and reconstruction of time dependent, high-dimensional data using sparse measurements, and apply it to the flow around a cylinder. Assuming the data lies near a low dimensional manifold (low-rank dynamics) in space and has periodic time dependency with a sparse number of Fourier modes, we employ compressive sensing for accurately classifying the dynamical regime. We further show that we can reconstruct the full spatio-temporal behavior with these limited measurements, extending previous results of compressive sensing that apply for only a single snapshot of data. The method can be used for building improved reduced-order models and designing sampling/measurement strategies that leverage time asynchrony.


title: A mutation-selection model for evolution of random dispersal

abstract: We consider a mutation-selection model of a population structured by the spatial variables and a trait variable which is the diffusion rate. Competition for resource is local in spatial variables, but nonlocal in the trait variable. We focus on the asymptotic profile of positive steady state solutions. Our result shows that in the limit of small mutation rate, the solution remains regular in the spatial variables and yet concentrates in the trait variable and forms a Dirac mass supported at the lowest diffusion rate. [Hastings, Theor. Pop. Biol. 24, 244-251, 1983] and [Dockery et al., J. Math. Biol. 37, 61-83, 1998] showed that for two competing species in spatially heterogeneous but temporally constant environment, the slower diffuser always prevails, if all other things are held equal. Our result suggests that their findings may hold for arbitrarily many traits.


title: Local Langlands Duality and a Duality of Conformal Field Theories

abstract: We show that the numerical local Langlands duality for GL_n and the T - duality of two-dimensional quantum gravity arise from one and the same symmetry principle. The unifying theme is that the local Fourier transform in both its l-adic and complex incarnation gives rise to symmetries of arithmetic and geometric local Langlands parameters.


title: Saturn's aurora observed by the Cassini camera at visible wavelengths

abstract: The first observations of Saturn's visible-wavelength aurora were made by the Cassini camera. The aurora was observed between 2006 and 2013 in the northern and southern hemispheres. The color of the aurora changes from pink at a few hundred km above the horizon to purple at 1000-1500 km above the horizon. The spectrum observed in 9 filters spanning wavelengths from 250 nm to 1000 nm has a prominent H-alpha line and roughly agrees with laboratory simulated auroras. Auroras in both hemispheres vary dramatically with longitude. Auroras form bright arcs between 70 and 80 degree latitude north and between 65 and 80 degree latitude south, which sometimes spiral around the pole, and sometimes form double arcs. A large 10,000-km-scale longitudinal brightness structure persists for more than 100 hours. This structure rotates approximately together with Saturn. On top of the large steady structure, the auroras brighten suddenly on the timescales of a few minutes. These brightenings repeat with a period of about 1 hour. Smaller, 1000-km-scale structures may move faster or lag behind Saturn's rotation on timescales of tens of minutes. The persistence of nearly-corotating large bright longitudinal structure in the auroral oval seen in two movies spanning 8 and 11 rotations gives an estimate on the period of 10.65 $\pm$0.15 h for 2009 in the northern oval and 10.8$\pm$ 0.1 h for 2012 in the southern oval. The 2009 north aurora period is close to the north branch of Saturn Kilometric Radiation (SKR) detected at that time.


title: Light stop in the MSSM after LHC Run 1

abstract: The discovery of a Higgs boson with a mass of 126 GeV at the LHC when combined with the non-observation of new physics both in direct and indirect searches imposes strong constraints on supersymmetric models and in particular on the top squark sector. The experiments for direct detection of dark matter have provided with yet more constraints on the neutralino LSP mass and its interactions. After imposing limits from the Higgs, flavour and dark matter sectors, we examine the feasibility for a light stop in the context of the pMSSM, in light of current results for stop and other SUSY searches at the LHC. We only require that the neutralino dark matter explains a fraction of the cosmologically measured dark matter abundance. We find that a stop with mass below $\sim$ 500 GeV is still allowed. We further study various probes of the light stop scenario that could be performed at the LHC Run - II either through direct searches for the light and heavy stop, or SUSY searches not currently available in simplified model results. Moreover we study the characteristics of heavy Higgs for the points in the parameter space allowed by all the available constraints and illustrate the region with large cross sections to fermionic or electroweakino channels. Finally we show that nearly all scenarios with a small stop$-$LSP mass difference will be tested by Xenon1T provided the NLSP is a chargino, thus probing a region hard to access at the LHC.


title: Cutoff for Extensions of Massive Gravity and Bi-Gravity

abstract: Recently there has been interest in extending ghost-free massive gravity, bi-gravity, and multi-gravity by including non-standard kinetic terms and matter couplings. We first review recent proposals for this class of extensions, emphasizing how modifications of the kinetic and potential structure of the graviton and modifications of the coupling to matter are related. We then generalize existing no-go arguments in the metric language to the vielbein language in second-order form. We give an ADM argument to show that the most promising extensions to the kinetic term and matter coupling contain a Boulware-Deser ghost. However, as recently emphasized, we may still be able to view these extensions as effective field theories below some cutoff scale. To address this possibility, we show that there is a decoupling limit where a ghost appears for a wide class of matter couplings and kinetic terms. In particular, we show that there is a decoupling limit where the linear effective vielbein matter coupling contains a ghost. Using the insight we gain from this decoupling limit analysis, we place an upper bound on the cutoff for the linear effective vielbein coupling. This result can be generalized to new kinetic interactions in the vielbein language in second-order form. Combined with recent results, this provides a strong uniqueness argument on the form of ghost-free massive gravity, bi-gravity, and multi-gravity.


title: QCD corrections to $B \to \pi$ form factors from light-cone sum rules

abstract: We compute perturbative corrections to $B \to \pi$ form factors from QCD light-cone sum rules with $B$-meson distribution amplitudes. Applying the method of regions we demonstrate factorization of the vacuum-to-$B$-meson correlation function defined with an interpolating current for pion, at one-loop level, explicitly in the heavy quark limit. The short-distance functions in the factorization formulae of the correlation function involves both hard and hard-collinear scales; and these functions can be further factorized into hard coefficients by integrating out the hard fluctuations and jet functions encoding the hard-collinear information. Resummation of large logarithms in the short-distance functions is then achieved via the standard renormalization-group approach. We further show that structures of the factorization formulae for $f_{B \pi}^{+}(q^2)$ and $f_{B \pi}^{0}(q^2)$ at large hadronic recoil from QCD light-cone sum rules match that derived in QCD factorization. In particular, we perform an exploratory phenomenological analysis of $B \to \pi$ form factors, paying attention to various sources of perturbative and systematic uncertainties, and extract $|V_{ub}|= \left(3.05^{+0.54}_{-0.38} |_{\rm th.} \pm 0.09 |_{\rm exp.}\right) \times 10^{-3}$ with the inverse moment of the $B$-meson distribution amplitude $\phi_B^{+}(\omega)$ determined by reproducing $f_{B \pi}^{+}(q^2=0)$ obtained from the light-cone sum rules with $\pi$ distribution amplitudes. Furthermore, we present the invariant-mass distributions of the lepton pair for $B \to \pi \ell \nu_{\ell}$ ($\ell= \mu \,, \tau$) in the whole kinematic region. Finally, we discuss non-valence Fock state contributions to the $B \to \pi$ form factors $f_{B \pi}^{+}(q^2)$ and $f_{B \pi}^{0}(q^2)$ in brief.


title: A generalized enhanced Fourier law and underlying connections to major frameworks for quasi-ballistic phonon transport

abstract: An enhanced Fourier law (EFL) that accounts for quasi-ballistic phonon transport effects in a formulation entirely in terms of physical observables, is derived from the Boltzmann transport equation, assuming a gray population of quasi-ballistic phonon modes. This equation is generalized to an arbitrary phonon population. Other phonon transport models are analyzed in the context of the generalized EFL and connections are made between the generalized EFL and other models, revealing the essential unity of seemingly disparate models reported in the literature.


title: Concentration and regularization of random graphs

abstract: This paper studies how close random graphs are typically to their expectations. We interpret this question through the concentration of the adjacency and Laplacian matrices in the spectral norm. We study inhomogeneous Erd\"os-R\'enyi random graphs on $n$ vertices, where edges form independently and possibly with different probabilities $p_{ij}$. Sparse random graphs whose expected degrees are $o(\log n)$ fail to concentrate; the obstruction is caused by vertices with abnormally high and low degrees. We show that concentration can be restored if we regularize the degrees of such vertices, and one can do this in various ways. As an example, let us reweight or remove enough edges to make all degrees bounded above by $O(d)$ where $d=\max np_{ij}$. Then we show that the resulting adjacency matrix $A'$ concentrates with the optimal rate: $\|A' - \mathbb{E} A\| = O(\sqrt{d})$. Similarly, if we make all degrees bounded below by $d$ by adding weight $d/n$ to all edges, then the resulting Laplacian concentrates with the optimal rate: $\|L(A') - L(\mathbb{E} A')\| = O(1/\sqrt{d})$. Our approach is based on Grothendieck-Pietsch factorization, using which we construct a new decomposition of random graphs. We illustrate the concentration results with an application to the community detection problem in the analysis of networks.


title: On the time-dependent grade-two model for the magnetohydrodynamic flow: 2D case

abstract: In this paper we discuss the MHD flow of a second grade fluid, in particular we prove the existence and uniqueness of a weak solution of a time-dependent grade two fluid model in a two-dimensional Lipschitz domain. We follow the methodology of Girault and Saadouni (Ref.[3]), i.e , we use a constructive method which can be adapted to the numerical analysis of finite-element schemes for solving this problem numerically.


title: Sample-Optimal Density Estimation in Nearly-Linear Time

abstract: We design a new, fast algorithm for agnostically learning univariate probability distributions whose densities are well approximated by piecewise polynomial functions. Let $f$ be the density function of an arbitrary univariate distribution, and suppose that $f$ is $\mathrm{OPT}$-close in $L_1$-distance to an unknown piecewise polynomial function with $t$ interval pieces and degree $d$. Our algorithm draws $n = O(t(d+1)/\epsilon^2)$ samples from $f$, runs in time $\tilde{O}(n \cdot \mathrm{poly}(d))$, and with probability at least $9/10$ outputs an $O(t)$-piecewise degree-$d$ hypothesis $h$ that is $4 \cdot \mathrm{OPT} +\epsilon$ close to $f$. Our general algorithm yields (nearly) sample-optimal and nearly-linear time estimators for a wide range of structured distribution families over both continuous and discrete domains in a unified way. For most of our applications, these are the first sample-optimal and nearly-linear time estimators in the literature. As a consequence, our work resolves the sample and computational complexities of a broad class of inference tasks via a single "meta-algorithm". Moreover, we experimentally demonstrate that our algorithm performs very well in practice. Our algorithm consists of three "levels": (i) At the top level, we employ an iterative greedy algorithm for finding a good partition of the real line into the pieces of a piecewise polynomial. (ii) For each piece, we show that the sub-problem of finding a good polynomial fit on the current interval can be solved efficiently with a separation oracle method. (iii) We reduce the task of finding a separating hyperplane to a combinatorial problem and give an efficient algorithm for this problem. Combining these three procedures gives a density estimation algorithm with the claimed guarantees.

categories: cs.DS cs.IT cs.LG math.IT math.ST stat.TH


title: Perturbative Vacuum Wavefunctional for Gauge Theories in the Milne Space

abstract: The spectrum of vacuum fluctuations in the Milne space (i.e. the tau-eta coordinate system) is an important ingredient in the thermalization studies in relativistic heavy ion collisions. In this paper, the Schrodinger functional for the gauge theory perturbative vacuum is derived for the Milne space. The Wigner-transform of the corresponding vacuum density functional is also found together with the propagators. We finally identify the fluctuation spectrum in vacuum, and show the equivalence between the present approach and the symplectic product based method.


title: Mutual Dependence: A Novel Method for Computing Dependencies Between Random Variables

abstract: In data science, it is often required to estimate dependencies between different data sources. These dependencies are typically calculated using Pearson's correlation, distance correlation, and/or mutual information. However, none of these measures satisfy all the Granger's axioms for an "ideal measure". One such ideal measure, proposed by Granger himself, calculates the Bhattacharyya distance between the joint probability density function (pdf) and the product of marginal pdfs. We call this measure the mutual dependence. However, to date this measure has not been directly computable from data. In this paper, we use our recently introduced maximum likelihood non-parametric estimator for band-limited pdfs, to compute the mutual dependence directly from the data. We construct the estimator of mutual dependence and compare its performance to standard measures (Pearson's and distance correlation) for different known pdfs by computing convergence rates, computational complexity, and the ability to capture nonlinear dependencies. Our mutual dependence estimator requires fewer samples to converge to theoretical values, is faster to compute, and captures more complex dependencies than standard measures.


title: Projections and Phase retrieval

abstract: We characterize collections of orthogonal projections for which it is possible to reconstruct a vector from the magnitudes of the corresponding projections. As a result we are able to show that in an $M$-dimensional real vector space a vector can be reconstructed from the magnitudes of its projections onto a generic collection of $N \geq 2M-1$ subspaces. We also show that this bound is sharp when $N = 2^k +1$. The results of this paper answer a number of questions raised in \cite{CCPW:13}.


title: The Quantum Hilbert Hotel

abstract: In 1924 David Hilbert conceived a paradoxical tale involving a hotel with an infinite number of rooms to illustrate some aspects of the mathematical notion of "infinity". In continuous-variable quantum mechanics we routinely make use of infinite state spaces: here we show that such a theoretical apparatus can accommodate an analog of Hilbert's hotel paradox. We devise a protocol that, mimicking what happens to the guests of the hotel, maps the amplitudes of an infinite eigenbasis to twice their original quantum number in a coherent and deterministic manner, producing infinitely many unoccupied levels in the process. We demonstrate the feasibility of the protocol by experimentally realising it on the orbital angular momentum of a paraxial field. This new non-Gaussian operation may be exploited for example for enhancing the sensitivity of N00N states, for increasing the capacity of a channel or for multiplexing multiple channels into a single one.


title: Exceptional sets for nonuniformly expanding maps

abstract: Given a rational map of the Riemann sphere and a subset $A$ of its Julia set, we study the $A$-exceptional set, that is, the set of points whose orbit does not accumulate at $A$. We prove that if the topological entropy of $A$ is less than the topological entropy of the full system then the $A$-exceptional set has full topological entropy. Furthermore, if the Hausdorff dimension of $A$ is smaller than the dynamical dimension of the system then the Hausdorff dimension of the $A$-exceptional set is larger than or equal to the dynamical dimension, with equality in the particular case when the dynamical dimension and the Hausdorff dimension coincide. We discuss also the case of a general conformal $C^{1+\alpha}$ dynamical system and, in particular, certain multimodal interval maps on their Julia set.


title: Characterisation of Strongly Stable Matchings

abstract: An instance of a strongly stable matching problem (SSMP) is an undirected bipartite graph $G=(A \cup B, E)$, with an adjacency list of each vertex being a linearly ordered list of ties, which are subsets of vertices equally good for a given vertex. Ties are disjoint and may contain one vertex. A matching $M$ is a set of vertex-disjoint edges. An edge $(x,y) \in E \setminus M$ is a {\em blocking edge} for $M$ if $x$ is either unmatched or strictly prefers $y$ to its current partner in $M$, and $y$ is either unmatched or strictly prefers $x$ to its current partner in $M$ or is indifferent between them. A matching is {\em strongly stable} if there is no blocking edge with respect to it. We present an algorithm for the generation of all strongly stable matchings, thus solving an open problem already stated in the book by Gusfield and Irving \cite{GI}. It has previously been shown that strongly stable matchings form a distributive lattice and although the number of strongly stable matchings can be exponential in the number of vertices, we show that there exists a partial order with $O(m)$ elements representing all strongly stable matchings, where $m$ denotes the number of edges in the graph. We give two algorithms that construct two such representations: one in $O(nm^2)$ time and the other in $O(nm)$ time, where $n$ denotes the number of vertices in the graph. Note that the construction of the second representation has the same time complexity as that of computing a single strongly stable matching.


title: Persistent Phylogeny: A Galled-Tree and Integer Linear Programming Approach

abstract: The Persistent-Phylogeny Model is an extension of the widely studied Perfect-Phylogeny Model, encompassing a broader range of evolutionary phenomena. Biological and algorithmic questions concerning persistent phylogeny have been intensely investigated in recent years. In this paper, we explore two alternative approaches to the persistent-phylogeny problem that grow out of our previous work on perfect phylogeny, and on galled trees. We develop an integer programming solution to the Persistent-Phylogeny Problem; empirically explore its efficiency; and empirically explore the utility of using fast algorithms that recognize galled trees, to recognize persistent phylogeny. The empirical results identify parameter ranges where persistent phylogeny are galled trees with high frequency, and show that the integer programming approach can efficiently identify persistent phylogeny of much larger size than has been previously reported.


title: Photocurrent generation with two-dimensional van der Waals semiconductors

abstract: Two-dimensional (2D) materials have attracted a great deal of interest in recent years. This family of materials allows for the realization of versatile electronic devices and holds promise for next-generation (opto)electronics. Their electronic properties strongly depend on the number of layers, making them interesting from a fundamental standpoint. For electronic applications, semiconducting 2D materials benefit from sizable mobilities and large on/off ratios, due to the large modulation achievable via the gate field-effect. Moreover, being mechanically strong and flexible, these materials can withstand large strain (>10\%) before rupture, making them interesting for strain engineering and flexible devices. Even in their single layer form, semiconducting 2D materials have demonstrated efficient light absorption, enabling large responsivity in photodetectors. Therefore, semiconducting layered 2D materials are strong candidates for optoelectronic applications, especially for photodetection. Here, we review the state-of-the-art in photodetectors based on semiconducting 2D materials, focusing on the transition metal dichalcogenides, novel van der Waals materials, black phosphorus, and heterostructures.


title: Concentration Independent Random Number Generation in Tile Self-Assembly

abstract: In this paper we introduce the \emph{robust random number generation} problem where the goal is to design an abstract tile assembly system (aTAM system) whose terminal assemblies can be split into $n$ partitions such that a resulting assembly of the system lies within each partition with probability 1/$n$, regardless of the relative concentration assignment of the tile types in the system. First, we show this is possible for $n=2$ (a \emph{robust fair coin flip}) within the aTAM, and that such systems guarantee a worst case $\mathcal{O}(1)$ space usage. We accompany our primary construction with variants that show trade-offs in space complexity, initial seed size, temperature, tile complexity, bias, and extensibility, and also prove some negative results. As an application, we combine our coin-flip system with a result of Chandran, Gopalkrishnan, and Reif to show that for any positive integer $n$, there exists a $\mathcal{O}(\log n)$ tile system that assembles a constant-width linear assembly of expected length $n$ for any concentration assignment. We then extend our robust fair coin flip result to solve the problem of robust random number generation in the aTAM for all $n$. Two variants of robust random bit generation solutions are presented: an unbounded space solution and a bounded space solution which incurs a small bias. Further, we consider the harder scenario where tile concentrations change arbitrarily at each assembly step and show that while this is not possible in the aTAM, the problem can be solved by exotic tile assembly models from the literature.


title: Mechanism for collective cell alignment in Myxococcus xanthus bacteria

abstract: Myxococcus xanthus cells self-organize into aligned groups, clusters, at various stages of their lifecycle. Formation of these clusters is crucial for the complex dynamic multi-cellular behavior of these bacteria. However, the mechanism underlying the cell alignment and clustering is not fully understood. Motivated by studies of clustering in self-propelled rods, we hypothesized that M. xanthus cells can align and form clusters through pure mechanical interactions among cells and between cells and substrate. We test this hypothesis using an agent-based simulation framework in which each agent is based on the biophysical model of an individual M. xanthus cell. We show that model agents, under realistic cell flexibility values, can align and form cell clusters but only when periodic reversals of cell directions are suppressed. However, by extending our model to introduce the observed ability of cells to deposit and follow slime trails, we show that effective trail-following leads to clusters in reversing cells. Furthermore, we conclude that mechanical cell alignment combined with slime-trail-following is sufficient to explain the distinct clustering behaviors observed for wild-type and non-reversing M. xanthus mutants in recent experiments. Our results are robust to variation in model parameters, match the experimentally observed trends and can be applied to understand surface motility patterns of other bacterial species.


title: Group buying with bundle discounts: computing efficient, stable and fair solutions

abstract: We model a market in which nonstrategic vendors sell items of different types and offer bundles at discounted prices triggered by demand volumes. Each buyer acts strategically in order to maximize her utility, given by the difference between product valuation and price paid. Buyers report their valuations in terms of reserve prices on sets of items, and might be willing to pay prices different than the market price in order to subsidize other buyers and to trigger discounts. The resulting price discrimination can be interpreted as a redistribution of the total discount. We consider a notion of stability that looks at unilateral deviations, and show that efficient allocations - the ones maximizing the social welfare - can be stabilized by prices that enjoy desirable properties of rationality and fairness. These dictate that buyers pay higher prices only to subsidize others who contribute to the activation of the desired discounts, and that they pay premiums over the discounted price proportionally to their surplus - the difference between their current utility and the utility of their best alternative. Therefore, the resulting price discrimination appears to be desirable to buyers. Building on this existence result, and letting N, M and c be the numbers of buyers, vendors and product types, we propose a O(N^2+NM^c) algorithm that, given an efficient allocation, computes prices that are rational and fair and that stabilize the market. The algorithm first determines the redistribution of the discount between groups of buyers with an equal product choice, and then computes single buyers' prices. Our results show that if a desirable form of price discrimination is implemented then social efficiency and stability can coexists in a market presenting subtle externalities, and computing individual prices from market prices is tractable.


title: Billiards with Bombs

abstract: In this paper, we define a variant of billiards in which the ball bounces around a square grid erasing walls as it goes. We prove that there exist periodic tunnels with arbitrarily large period from any possible starting point, that there exist non-periodic tunnels from any possible starting point, and that there are versions of the problem for which the same starting point and initial direction result in periodic tunnels of arbitrarily large period. We conjecture that there exist starting conditions which do not lead to tunnels, justify the conjecture with simulation evidence, and discuss the difficulty of proving it.


title: Multi-Version Coding - An Information Theoretic Perspective of Consistent Distributed Storage

abstract: In applications of distributed storage systems to distributed computing and implementation of key- value stores, the following property, usually referred to as consistency in computer science and engineering, is an important requirement: as the data stored changes, the latest version of the data must be accessible to a client that connects to the storage system. An information theoretic formulation called multi-version coding is introduced in the paper, in order to study storage costs of consistent distributed storage systems. Multi-version coding is characterized by {\nu} totally ordered versions of a message, and a storage system with n servers. At each server, values corresponding to an arbitrary subset of the {\nu} versions are received and encoded. For any subset of c servers in the storage system, the value corresponding to the latest common version, or a later version as per the total ordering, among the c servers is required to be decodable. An achievable multi-version code construction via linear coding and a converse result that shows that the construction is approximately tight, are provided. An implication of the converse is that there is an inevitable price, in terms of storage cost, to ensure consistency in distributed storage systems.


title: Model-based reinforcement learning for infinite-horizon approximate optimal tracking

abstract: This paper provides an approximate online adaptive solution to the infinite-horizon optimal tracking problem for control-affine continuous-time nonlinear systems with unknown drift dynamics. Model-based reinforcement learning is used to relax the persistence of excitation condition. Model-based reinforcement learning is implemented using a concurrent learning-based system identifier to simulate experience by evaluating the Bellman error over unexplored areas of the state space. Tracking of the desired trajectory and convergence of the developed policy to a neighborhood of the optimal policy are established via Lyapunov-based stability analysis. Simulation results demonstrate the effectiveness of the developed technique.


title: Invariance, existence and uniqueness of solutions of nonlinear valuation PDEs and FBSDEs inclusive of credit risk, collateral and funding costs

abstract: We study conditions for existence, uniqueness and invariance of the comprehensive nonlinear valuation equations first introduced in Pallavicini et al (2011). These equations take the form of semilinear PDEs and Forward-Backward Stochastic Differential Equations (FBSDEs). After summarizing the cash flows definitions allowing us to extend valuation to credit risk and default closeout, including collateral margining with possible re-hypothecation, and treasury funding costs, we show how such cash flows, when present-valued in an arbitrage free setting, lead to semi-linear PDEs or more generally to FBSDEs. We provide conditions for existence and uniqueness of such solutions in a viscosity and classical sense, discussing the role of the hedging strategy. We show an invariance theorem stating that even though we start from a risk-neutral valuation approach based on a locally risk-free bank account growing at a risk-free rate, our final valuation equations do not depend on the risk free rate. Indeed, our final semilinear PDE or FBSDEs and their classical or viscosity solutions depend only on contractual, market or treasury rates and we do not need to proxy the risk free rate with a real market rate, since it acts as an instrumental variable. The equations derivations, their numerical solutions, the related XVA valuation adjustments with their overlap, and the invariance result had been analyzed numerically and extended to central clearing and multiple discount curves in a number of previous works, including Pallavicini et al (2011), Pallavicini et al (2012), Brigo et al (2013), Brigo and Pallavicini (2014), and Brigo et al (2014).


title: The Parabolic variance (PVAR), a wavelet variance based on least-square fit

abstract: This article introduces the Parabolic Variance (PVAR), a wavelet variance similar to the Allan variance, based on the Linear Regression (LR) of phase data. The companion article arXiv:1506.05009 [physics.ins-det] details the $\Omega$ frequency counter, which implements the LR estimate. The PVAR combines the advantages of AVAR and MVAR. PVAR is good for long-term analysis because the wavelet spans over $2 \tau$, the same of the AVAR wavelet; and good for short-term analysis because the response to white and flicker PM is $1/\tau^3$ and $1/\tau^2$, same as the MVAR. After setting the theoretical framework, we study the degrees of freedom and the confidence interval for the most common noise types. Then, we focus on the detection of a weak noise process at the transition - or corner - where a faster process rolls off. This new perspective raises the question of which variance detects the weak process with the shortest data record. Our simulations show that PVAR is a fortunate tradeoff. PVAR is superior to MVAR in all cases, exhibits the best ability to divide between fast noise phenomena (up to flicker FM), and is almost as good as AVAR for the detection of random walk and drift.


title: A non-conforming domain decomposition approximation for the Helmholtz screen problem with hypersingular operator

abstract: We present and analyze a non-conforming domain decomposition approximation for a hypersingular operator governed by the Helmholtz equation in three dimensions. This operator appears when considering the corresponding Neumann problem in unbounded domains exterior to open surfaces. We consider small wave numbers and low-order approximations with Nitsche coupling across interfaces. Under appropriate assumptions on mapping properties of the weakly singular and hypersingular operators with Helmholtz kernel, we prove that this method converges almost quasi-optimally. Numerical experiments confirm our error estimate.


title: Slow light transmission in one-dimensional periodic structures

abstract: We have analyzed the transmission properties of pulses through one-dimensional periodic structures in order to systematically explore the best conditions to achieve the maximum delay with the minimum possible distortion. In the absence of absorption and no layer variation, the transmission coefficient $t_N$ can be well approximated by a sum of Lorentzian resonances. The ratio between their width $\Gamma_{\rm r}$ and their separation $\Delta\omega_{\rm r}$ is a crucial parameter to characterize the distortion of the transmitted pulse. For typical values of the parameters used in telecommunications and high index of refraction contrasts $n_2 / n_1$, the distortion of the transmitted pulse is unacceptably large for frequencies near the edge of the transmission window. We estimate fractional delays achievables in terms of the central frequency used and the pulse bandwidth.


title: Space-time gradient metasurfaces

abstract: Metasurfaces characterized by a transverse gradient of local impedance have recently opened exciting directions for light manipulation at the nanoscale. Here we add a temporal gradient to the picture, showing that spatio-temporal variations over a surface may largely extend the degree of light manipulation in metasurfaces, and break several of their constraints associated to symmetries. As an example, we synthesize a non-reciprocal classical analogue to electromagnetic induced transparency, opening a narrow window of one-way transmission in an otherwise opaque surface. These properties pave the way to magnetic free, planarized non-reciprocal ultrathin surfaces for free-space isolation.


title: Robust Covariance and Scatter Matrix Estimation under Huber's Contamination Model

abstract: Covariance matrix estimation is one of the most important problems in statistics. To accommodate the complexity of modern datasets, it is desired to have estimation procedures that not only can incorporate the structural assumptions of covariance matrices, but are also robust to outliers from arbitrary sources. In this paper, we define a new concept called matrix depth and then propose a robust covariance matrix estimator by maximizing the empirical depth function. The proposed estimator is shown to achieve minimax optimal rate under Huber's $\epsilon$-contamination model for estimating covariance/scatter matrices with various structures including bandedness and sparsity.


title: Universal Central Extension of the Lie Algebra of Hamiltonian Vector Fields

abstract: We determine the universal central extension of the Lie algebra of hamiltonian vector fields, thereby classifying its central extensions. Furthermore, we classify the central extensions of the Lie algebra of symplectic vector fields, of the Poisson Lie algebra, and of its compactly supported version.


title: Route to hyperchaos in Rayleigh-Benard convection

abstract: Transition to hyperchaotic regimes in Rayleigh-Benard convection in a square periodicity cell is studied by three-dimensional numerical simulations. By fixing the Prandtl number at P=0.3 and varying the Rayleigh number as a control parameter, a bifurcation diagram is constructed where a route to hyperchaos involving quasiperiodic regimes with two and three incommensurate frequencies, multistability, chaotic intermittent attractors and a sequence of boundary and interior crises is shown. The three largest Lyapunov exponents exhibit a linear scaling with the Rayleigh number and are positive in the final hyperchaotic attractor. Thus, a route to weak turbulence is found.


title: An Iterative On-Line Mechanism for Demand-Side Aggregation

abstract: This paper considers a demand-side aggregation scheme specifically for large numbers of small loads, such as households and small and medium-sized businesses. We introduce a novel auction format, called a staggered clock-proxy auction (SCPA), for on-line scheduling of these loads. This is a two phase format, consisting of: a sequence of overlapping iterative ascending-price clock auctions, one for each time-slot over a finite decision horizon, and; a set of proxy auctions that begin at the termination of each individual clock auction, and which determine the final price and allocation for each time-slot. The overlapping design of the clock phases grant bidders the ability to effectively bid on inter-temporal bundles of electricity use, thereby focusing on the most relevant parts of the price-quantity space. Since electricity is a divisible good, the proxy auction uses demand-schedule bids, which the aggregator uses to compute a uniform-price partial competitive equilibrium for each time slot. We show that, under mild assumptions on the bidders' utilities functions, the proxy phase implements the Vickrey-Clarke-Groves outcome, which makes straightforward bidding in the proxy phase a Bayes-Nash equilibrium. Furthermore, we demonstrate the SCPA in a scenario comprised of household agents with three different utility function types, and show how the mechanism enables efficient on-line energy use scheduling.


title: On the Skolem Problem for Continuous Linear Dynamical Systems

abstract: The Continuous Skolem Problem asks whether a real-valued function satisfying a linear differential equation has a zero in a given interval of real numbers. This is a fundamental reachability problem for continuous linear dynamical systems, such as linear hybrid automata and continuous-time Markov chains. Decidability of the problem is currently open---indeed decidability is open even for the sub-problem in which a zero is sought in a bounded interval. In this paper we show decidability of the bounded problem subject to Schanuel's Conjecture, a unifying conjecture in transcendental number theory. We furthermore analyse the unbounded problem in terms of the frequencies of the differential equation, that is, the imaginary parts of the characteristic roots. We show that the unbounded problem can be reduced to the bounded problem if there is at most one rationally linearly independent frequency, or if there are two rationally linearly independent frequencies and all characteristic roots are simple. We complete the picture by showing that decidability of the unbounded problem in the case of two (or more) rationally linearly independent frequencies would entail a major new effectiveness result in Diophantine approximation, namely computability of the Diophantine-approximation types of all real algebraic numbers.


title: Stein-Tomas Restriction Theorem via Spectral Measure on Metric Measure Spaces

abstract: The Stein-Tomas restriction theorem on Euclidean space says one can meaningfully restrict $\hat{f}$ to the unit sphere of $\mathbb{R}^n$ provided $f \in L^p(\mathbb{R}^n)$ with $1 < p < 2$. This result can be rewritten in terms of the estimates for the spectral measure of Laplacian. Guillarmou, Hassell and Sikora formulated a sufficient condition of the restriction theorem, via spectral measure, on abstract metric measure spaces. But they only proved the result in a special case. The present note aims to give a complete proof. In the end, it will be applied to the restriction theorem on asymptotically conic manifolds.


title: Approximations of Bond and Swaption Prices in a Black-Karasi\'{n}ski Model

abstract: We derive semi-analytic approximation formulae for bond and swaption prices in a Black-Karasi\'{n}ski interest rate model. Approximations are obtained using a novel technique based on the Karhunen-Lo\`{e}ve expansion. Formulas are easily computable and prove to be very accurate in numerical tests. This makes them useful for numerically efficient calibration of the model.


title: Statistical Machine Translation Features with Multitask Tensor Networks

abstract: We present a three-pronged approach to improving Statistical Machine Translation (SMT), building on recent success in the application of neural networks to SMT. First, we propose new features based on neural networks to model various non-local translation phenomena. Second, we augment the architecture of the neural network with tensor layers that capture important higher-order interaction among the network units. Third, we apply multitask learning to estimate the neural network parameters jointly. Each of our proposed methods results in significant improvements that are complementary. The overall improvement is +2.7 and +1.8 BLEU points for Arabic-English and Chinese-English translation over a state-of-the-art system that already includes neural network features.


title: What Moser Could Have Asked: Counting Hamilton Cycles in Tournaments

abstract: Moser asked for a construction of explicit tournaments on $n$ vertices having at least $(\frac{n}{3e})^n$ Hamilton cycles. We show that he could have asked for rather more.


title: Lattice Boltzmann modeling of self-propelled Leidenfrost droplets on ratchet surfaces

abstract: In this paper, the self-propelled motion of Leidenfrost droplets on ratchet surfaces is numerically investigated with a thermal multiphase lattice Boltzmann model with liquid-vapor phase change. The capability of the model for simulating evaporation is validated via the D2 law. Using the model, we first study the performances of Leidenfrost droplets on horizontal ratchet surfaces. It is numerically shown that the motion of self-propelled Leidenfrost droplets on ratchet surfaces is owing to the asymmetry of the ratchets and the vapor flows beneath the droplets. It is found that the Leidenfrost droplets move in the direction toward the slowly inclined side from the ratchet peaks, which agrees with the direction of droplet motion in experiments [Linke et al., Phys. Rev. Lett., 2006, 96, 154502]. Moreover, the influences of the ratchet aspect ratio are investigated. For the considered ratchet surfaces, a critical value of the ratchet aspect ratio is approximately found, which corresponds to the maximum droplet moving velocity. Furthermore, the processes that the Leidenfrost droplets climb uphill on inclined ratchet surfaces are also studied. Numerical results show that the maximum inclination angle at which a Leidenfrost droplet can still climb uphill successfully is affected by the initial radius of the droplet.


title: A wirebasket preconditioner for the mortar boundary element method

abstract: We present and analyze a preconditioner of the additive Schwarz type for the mortar boundary element method. As a basic splitting, on each subdomain we separate the degrees of freedom related to its boundary from the inner degrees of freedom. The corresponding wirebasket-type space decomposition is stable up to logarithmic terms. For the blocks that correspond to the inner degrees of freedom standard preconditioners for the hypersingular integral operator on open boundaries can be used. For the boundary and interface parts as well as the Lagrangian multiplier space, simple diagonal preconditioners are optimal. Our technique applies to quasi-uniform and non-uniform meshes of shape-regular elements. Numerical experiments on triangular and quadrilateral meshes confirm theoretical bounds for condition and MINRES iteration numbers.


title: Luminosity--time and luminosity--luminosity correlations for GRB prompt and afterglow plateau emissions

abstract: We present an analysis of 123 Gamma-ray bursts (GRBs) with known redshifts possessing an afterglow plateau phase. We reveal that $L_a-T^{*}_a$ correlation between the X-ray luminosity $L_a$ at the end of the plateau phase and the plateau duration, $T^*_a$, in the GRB rest frame has a power law slope different, within more than 2 $\sigma$, from the slope of the prompt $L_{f}-T^{*}_{f}$ correlation between the isotropic pulse peak luminosity, $L_{f}$, and the pulse duration, $T^{*}_{f}$, from the time since the GRB ejection. Analogously, we show differences between the prompt and plateau phases in the energy-duration distributions with the afterglow emitted energy being on average $10\%$ of the prompt emission. Moreover, the distribution of prompt pulse versus afterglow spectral indexes do not show any correlation. In the further analysis we demonstrate that the $L_{peak}-L_a$ distribution, where $L_{peak}$ is the peak luminosity from the start of the burst, is characterized with a considerably higher Spearman correlation coefficient, $\rho=0.79$, than the one involving the averaged prompt luminosity, $L_{prompt}-L_a$, for the same GRB sample, yielding $\rho=0.60$. Since some of this correlation could result from the redshift dependences of the luminosities, namely from their cosmological evolution we use the Efron-Petrosian method to reveal the intrinsic nature of this correlation. We find that a substantial part of the correlation is intrinsic. We apply a partial correlation coefficient to the new de-evolved luminosities showing that the intrinsic correlation exists.


title: Variations in the Star Formation Efficiency of the Dense Molecular Gas across the Disks of Star-Forming Galaxies

abstract: We present a new survey of HCN(1-0) emission, a tracer of dense molecular gas, focused on the little-explored regime of normal star-forming galaxy disks. Combining HCN, CO, and infrared (IR) emission, we investigate the role of dense gas in Star Formation (SF), finding systematic variations in both the apparent dense gas fraction and the apparent SF efficiency (SFE) of dense gas. The latter may be unexpected, given the popularity of gas density threshold models to explain SF scaling relations. We used the IRAM 30-m telescope to observe HCN(1-0) across 29 nearby disk galaxies whose CO(2-1) emission has previously been mapped by the HERACLES survey. Because our observations span a range of galactocentric radii, we are able to investigate the properties of the dense gas as a function of local conditions. We focus on how the IR/CO, HCN/CO, and IR/HCN ratios (observational cognates of the SFE, dense gas fraction, and dense gas SFE) depend on the stellar surface density and the molecular/atomic ratio. The HCN/CO ratio correlates tightly with these two parameters across a range of 2.1 dex and increases in the high surface density parts of galaxies. Simultaneously, the IR/HCN ratio decreases systematically with these same parameters and is ~6-8 times lower near galaxy centers than in the outer regions. For fixed line-mass conversion factors, these results are incompatible with a simple model in which SF depends only on the gas mass above some density threshold. Only a specific set of environment-dependent conversion factors can render our observations compatible with such a model. Whole cloud models, such as the theory of turbulence regulated SF, do a better job of matching our data. We explore one such model in which variations in the Mach number and in the mean density would respectively drive the trends within galaxy disks and the differences between disk and merging galaxies (abridged).


title: Reinvigorating Excited State Coherences in $V$ Level Systems with the use of Trap States

abstract: We investigate potential schemes to reinvigorate coherences between excited states in a $V$ level system under various excitation conditions. Coherent pulsed, cw and noisy pulsed excitation are considered and the existence of a trap state is shown to allow for the resurgence of excited state coherence in the two pulsed scenarios. With a short trap time (strong decoherence), it is possible to restart excited state coherence in the pulsed scenarios. In the case of cw excitation with relaxation back to the ground state, the coherences between the excited states do not vanish in the long time limit but instead approach a non-zero value.


title: New Geometric Representations of the CMB 2pcf

abstract: When searching for deviations of statistical isotropy in CMB, a popular strategy is to write the two-point correlation function (2pcf) as the most general function of four spherical angles (i.e., two unit vectors) in the celestial sphere. Then, using a basis of bipolar spherical harmonics, statistical anisotropy will show up if and only if any coefficient of the expansion with non-trivial bipolar momentum is detected -- although this detection will not in general elucidate the origin of the anisotropy. In this work we show that two new sets of four angles and basis functions exist which completely specifies the 2pcf, while, at the same time, offering a clearer geometrical interpretation of the mechanisms generating the signal. Since the coefficients of these expansions are zero if and only if isotropy holds, they act as a simple and geometrically motivated null test of statistical isotropy, with the advantage of allowing cosmic variance to be controlled in a systematic way. We report the results of the application of these null tests to the latest temperature data released by the Planck collaboration.


title: A general lower bound for the asymptotic convergence factor

abstract: We provide a rather general and very simple to compute lower bound for the asymptotic convergence factor of compact subsets of the set of complex numbers with connected complement and finitely many connected components .


title: Polarization oscillations of near-field thermal emission

abstract: We consider the polarization of thermal emission in the near-field of various materials including dielectrics and metallic systems with resonant surface modes. We find that at thermal equilibrium, the degree of polarization exhibits spatial oscillations with a period of approximately half the optical wavelength, independent of material composition. This result contrasts with that of Setala, Kaivola and Friberg [Phys. Rev. Lett. 88, 123902 (2002)], who find monotonic decay of the degree of polarization for systems in local thermal equilibrium.


title: Geometric Bogomolov conjecture for nowhere degenerate abelian varieties of dimension $5$ with trivial trace

abstract: We prove that the geometric Bogomolov conjecture holds for nowhere degenerate abelian varieties of dimension $5$ with trivial trace. By this result together with our previous work, we see that the conjecture holds for an abelian variety such that the difference between the dimension of its maximal nowhere degenerate abelian subvariety and that of its trace equals $5$.


title: Nature of short high amplitude pulses in a periodic dissipative laminate metamaterial

abstract: We study the evolution of high amplitude stress pulses in periodic dissipative laminates taking into account the nonlinear constitutive equations of the components and their dissipative behavior. Aluminum and Tungsten laminate was taken as an example due to the large difference in acoustic impedances of aluminum and tungsten, the significant nonlinearity of aluminum constitutive equation at the investigated range of stresses, and its possible practical applications. Laminates with different cell size, which controls internal time scale, impacted by the plates with different thicknesses, determining the incoming pulse duration, were investigated. It has been observed that the ratio of the duration of the incoming pulse to the internal characteristic time determines the nature of the high amplitude dissipative propagating waves, the oscillatory shock like profile, the train of localized pulses or a single localized pulse. These localized quasistationary waves resemble solitary waves even in the presence of dissipation: the similar pulses emerged from different initial conditions, indicating that they are inherent properties of the corresponding laminates, their characteristic length scale is determined by the mesostructural scale and the stress amplitude, and they mostly recover their shapes after collision with phase shift and a linear relationship exists between their speed and amplitude. A theoretical description approximating the shape, length scale and speed of these high amplitude dissipative pulses was proposed based on the Korteweg de Vries type equation with a dispersive term dictated by the mesostructure. The nonlinear term is derived from nonlinear constitutive equations of Aluminum and Tungsten, which were extracted from their Hugoniot curves.


title: Joseph ideals and lisse minimal W-algebras

abstract: We consider a lifting of Joseph ideals for the minimal nilpotent orbit closure to the setting of affine Kac-Moody algebras and find new examples of affine vertex algebras whose associated varieties are minimal nilpotent orbit closures. As an application we obtain a new family of lisse ($C_2$-cofinite) W-algebras that are not coming from admissible representations of affine Kac-Moody algebras.


title: Quantifying Creativity in Art Networks

abstract: Can we develop a computer algorithm that assesses the creativity of a painting given its context within art history? This paper proposes a novel computational framework for assessing the creativity of creative products, such as paintings, sculptures, poetry, etc. We use the most common definition of creativity, which emphasizes the originality of the product and its influential value. The proposed computational framework is based on constructing a network between creative products and using this network to infer about the originality and influence of its nodes. Through a series of transformations, we construct a Creativity Implication Network. We show that inference about creativity in this network reduces to a variant of network centrality problems which can be solved efficiently. We apply the proposed framework to the task of quantifying creativity of paintings (and sculptures). We experimented on two datasets with over 62K paintings to illustrate the behavior of the proposed framework. We also propose a methodology for quantitatively validating the results of the proposed algorithm, which we call the "time machine experiment".


title: Reidemeister torsion of a 3-manifold obtained by a Dehn-surgery along the figure-eight knot

abstract: Let M be a closed 3-manifold obtained by Dehn surgery along the figure-eight knot. We give a formula of the Reisdemeisiter torsion of M for any irreducible SL(2;C)-representation. It is described as a rational expression of the trace of the image of the meridian.


title: Aggregate Particles in the Plumes of Enceladus

abstract: Estimates of the total particulate mass of the plumes of Enceladus are important to constrain theories of particle formation and transport at the surface and interior of the satellite. We revisit the calculations of Ingersoll & Ewald (2011), who estimated the particulate mass of the Enceladus plumes from strongly forward scattered light in Cassini ISS images. We model the plume as a combination of spherical particles and irregular aggregates resulting from the coagulation of spherical monomers, the latter of which allows for plumes of lower particulate mass. Though a continuum of solutions are permitted by the model, the best fits to the ISS data consist either of low mass plumes composed entirely of small aggregates or high mass plumes composed of mostly spheres. The high particulate mass plumes have total particulate masses of (166 $\pm$ 42) $\times$ 10$^3$ kg, consistent with the results of Ingersoll & Ewald (2011). The low particulate mass plumes have masses of (25 $\pm$ 4) $\times$ 10$^3$ kg, leading to a solid to vapor mass ratio of 0.07 $\pm$ 0.01 for the plume. If indeed the plumes are made of such aggregates, then a vapor-based origin for the plume particles cannot be ruled out. Finally, we show that the residence time of the monomers inside the plume vents is sufficiently long for Brownian coagulation to form the aggregates before they are ejected to space.


title: Null lifts and projective dynamics

abstract: We describe natural Hamiltonian systems using projective geometry. The null lift procedure endows the tangent bundle with a projective structure where the null Hamiltonian is identified with a projective conic and induces a Weyl geometry. Projective transformations generate a set of known and new dualities between Hamiltonian systems, as for example the phenomenon of coupling-constant metamorphosis. We conclude outlining how this construction can be extended to the quantum case for Eisenhart-Duval lifts.


title: Cell structures for the Yokonuma-Hecke algebra and the algebra of braids and ties

abstract: We construct a faithful tensor representation for the Yokonuma-Hecke algebra Y, and use it to give a concrete isomorphism between Y and Shoji's modified Ariki-Koike algebra. We give a cellular basis for Y and show that the Jucys-Murphy elements for Y are JM-elements in the abstract sense. Finally, we construct a cellular basis for the Aicardi-Juyumaya algebra of braids and ties.


title: Tackling Exascale Software Challenges in Molecular Dynamics Simulations with GROMACS

abstract: GROMACS is a widely used package for biomolecular simulation, and over the last two decades it has evolved from small-scale efficiency to advanced heterogeneous acceleration and multi-level parallelism targeting some of the largest supercomputers in the world. Here, we describe some of the ways we have been able to realize this through the use of parallelization on all levels, combined with a constant focus on absolute performance. Release 4.6 of GROMACS uses SIMD acceleration on a wide range of architectures, GPU offloading acceleration, and both OpenMP and MPI parallelism within and between nodes, respectively. The recent work on acceleration made it necessary to revisit the fundamental algorithms of molecular simulation, including the concept of neighborsearching, and we discuss the present and future challenges we see for exascale simulation - in particular a very fine-grained task parallelism. We also discuss the software management, code peer review and continuous integration testing required for a project of this complexity.


title: Assessing Efficiency-Effectiveness Tradeoffs in Multi-Stage Retrieval Systems Without Using Relevance Judgments

abstract: Large-scale retrieval systems are often implemented as a cascading sequence of phases -- a first filtering step, in which a large set of candidate documents are extracted using a simple technique such as Boolean matching and/or static document scores; and then one or more ranking steps, in which the pool of documents retrieved by the filter is scored more precisely using dozens or perhaps hundreds of different features. The documents returned to the user are then taken from the head of the final ranked list. Here we examine methods for measuring the quality of filtering and preliminary ranking stages, and show how to use these measurements to tune the overall performance of the system. Standard top-weighted metrics used for overall system evaluation are not appropriate for assessing filtering stages, since the output is a set of documents, rather than an ordered sequence of documents. Instead, we use an approach in which a quality score is computed based on the discrepancy between filtered and full evaluation. Unlike previous approaches, our methods do not require relevance judgments, and thus can be used with virtually any query set. We show that this quality score directly correlates with actual differences in measured effectiveness when relevance judgments are available. Since the quality score does not require relevance judgments, it can be used to identify queries that perform particularly poorly for a given filter. Using these methods, we explore a wide range of filtering options using thousands of queries, categorize the relative merits of the different approaches, and identify useful parameter combinations.


title: Peeling Algorithm on Random Hypergraphs with Superlinear Number of Hyperedges

abstract: When we try to solve a system of linear equations, we can consider a simple iterative algorithm in which an equation including only one variable is chosen at each step, and the variable is fixed to the value satisfying the equation. The dynamics of this algorithm is captured by the peeling algorithm. Analyses of the peeling algorithm on random hypergraphs are required for many problems, e.g., the decoding threshold of low-density parity check codes, the inverting threshold of Goldreich's pseudorandom generator, the load threshold of cuckoo hashing, etc. In this work, we deal with random hypergraphs including superlinear number of hyperedges, and derive the tight threshold for the succeeding of the peeling algorithm. For the analysis, Wormald's method of differential equations, which is commonly used for analyses of the peeling algorithm on random hypergraph with linear number of hyperedges, cannot be used due to the superlinear number of hyperedges. A new method called the evolution of the moment generating function is proposed in this work.


title: Potentially crystalline deformation rings in the ordinary case

abstract: We study potentially crystalline deformation rings for a residual, ordinary Galois representation $\overline{\rho}: G_{\mathbb{Q}_p}\rightarrow \mathrm{GL}_3(\mathbb{F}_p)$. We consider deformations with Hodge-Tate weights $(0,1,2)$ and inertial type chosen to contain exactly one Fontaine-Laffaille modular weight for $\overline{\rho}$. We show that, in this setting, the potentially crystalline deformation space is formally smooth over $\mathbb{Z}_p$ and any potentially crystalline lift is ordinary. The proof requires an understanding of the condition imposed by the monodromy operator on Breuil modules with descent datum, in particular, that this locus mod p is formally smooth.


title: 68 Gbps quantum random number generation by measuring laser phase fluctuations

abstract: The speed of a quantum random number generator is essential for practical applications, such as high-speed quantum key distribution systems. Here, we push the speed of a quantum random number generator to 68 Gbps by operating a laser around its threshold level. To achieve the rate, not only high-speed photodetector and high sampling rate are needed, but also a very stable interferometer is required. A practical interferometer with active feedback instead of common temperature control is developed to meet requirement of stability. Phase fluctuations of the laser are measured by the interferometer with a photodetector, and then digitalized to raw random numbers with a rate of 80 Gbps. The min-entropy of the raw data is evaluated by modeling the system and is used to quantify the quantum randomness of the raw data. The bias of the raw data caused by other signals, such as classical and detection noises, can be removed by Toeplitz-matrix hashing randomness extraction. The final random numbers can pass through the standard randomness tests. Our demonstration shows that high-speed quantum random number generators are ready for practical usage.


title: Statistical transmutation in Floquet driven optical lattices

abstract: We show that interacting bosons in a periodically-driven two dimensional (2D) optical lattice may effectively exhibit fermionic statistics. The phenomenon is similar to the celebrated Tonks-Girardeau regime in 1D. The Floquet band of a driven lattice develops the moat shape, i.e. a minimum along a closed contour in the Brillouin zone. Such degeneracy of the kinetic energy favors fermionic quasiparticles. The statistical transmutation is achieved by the Chern-Simons flux attachment similar to the fractional quantum Hall case. We show that the velocity distribution of the released bosons is a sensitive probe of the fermionic nature of their stationary Floquet state.


title: A Distributed Algorithm for Demand Response with Mixed-Integer Variables

abstract: This letter presents a fast distributed algorithm for aggregating a large number of households with mixed-integer variables and intricate couplings between devices. The proposed fast distributed gradient algorithm is applied to the double smoothed dual function of the adopted DR model. The results also show that, with minimal parameter adjustments, the convergence of the dual objective exhibits the same behavior irrespective of the system size.


title: Current-Driven Motion of Magnetic Domain Wall with Many Bloch Lines

abstract: The current-driven motion of a domain wall (DW) in a ferromagnet with many Bloch lines (BLs) via the spin transfer torque is studied theoretically. It is found that the motion of BLs changes the current-velocity ($j$-$v$) characteristic dramatically. Especially, the critical current density to overcome the pinning force is reduced by the factor of the Gilbert damping coefficient $\alpha$ even compared with that of a skyrmion. This is in sharp contrast to the case of magnetic field driven motion, where the existence of BLs reduces the mobility of the DW.


title: The Brownian web, the Brownian net, and their universality

abstract: The Brownian web is a collection of one-dimensional coalescing Brownian motions starting from everywhere in space and time, and the Brownian net is a generalization that also allows branching. They appear in the diffusive scaling limits of many one-dimensional interacting particle systems with branching and coalescence. This article gives an introduction to the Brownian web and net, and how they arise in the scaling limits of various one-dimensional models, focusing mainly on coalescing random walks and random walks in i.i.d. space-time random environments. We will also briefly survey models and results connected to the Brownian web and net, including alternative topologies, population genetic models, true self-repelling motion, planar aggregation, drainage networks, oriented percolation, black noise and critical percolation. Some open questions are discussed at the end.


title: 3D gravity with dust: classical and quantum theory

abstract: We study the Einstein gravity and dust system in three spacetime dimensions as an example of a non-perturbative quantum gravity model with local degrees of freedom. We derive the Hamiltonian theory in the dust time gauge and show that it has a rich class of exact solutions. These include the Ba\~nados-Teitelboim-Zanelli black hole, static solutions with naked singularities and travelling wave solutions with dynamical horizons. We give a complete quantization of the wave sector of the theory, including a definition of a self-adjoint spacetime metric operator. This operator is used to demonstrate the quantization of deficit angle and the fluctuation of dynamical horizons.


title: Introduction to adic tropicalization

abstract: This is an expository article on the adic tropicalization of algebraic varieties. We outline joint work with Sam Payne in which we put a topology and structure sheaf of local topological rings on the exploded tropicalization. The resulting object, which blends polyhedral data of the tropicalization with algebraic data of the associated initial degenerations, is called the adic tropicalization. It satisfies a theorem of the form "Huber analytification is the limit of all adic tropicalizations." We explain this limit theorem in the present article, and illustrate connections between adic tropicalization and the curve complexes of O. Amini and M. Baker.


title: Nonlocal effective medium approximation for metallic nanorod metamaterials

abstract: We present an analytical nonlocal effective medium approximation to describe the optical nonlocal effects in metallic nanorod metamaterials based on Mie scattering theory. It is shown that the developed nonlocal effective medium theory can predict a coexistence state of two modes around the epsilon-near-zero region, where strong optical nonlocal effects lead to the behavior of both positive refraction and negative refraction in the nanorod metamaterials. Outside of the coexistence region, only one mode can be excited and its behavior can be well described using the local effective medium theory.


title: Network assisted analysis to reveal the genetic basis of autism

abstract: While studies show that autism is highly heritable, the nature of the genetic basis of this disorder remains illusive. Based on the idea that highly correlated genes are functionally interrelated and more likely to affect risk, we develop a novel statistical tool to find more potentially autism risk genes by combining the genetic association scores with gene co-expression in specific brain regions and periods of development. The gene dependence network is estimated using a novel partial neighborhood selection (PNS) algorithm, where node specific properties are incorporated into network estimation for improved statistical and computational efficiency. Then we adopt a hidden Markov random field (HMRF) model to combine the estimated network and the genetic association scores in a systematic manner. The proposed modeling framework can be naturally extended to incorporate additional structural information concerning the dependence between genes. Using currently available genetic association data from whole exome sequencing studies and brain gene expression levels, the proposed algorithm successfully identified 333 genes that plausibly affect autism risk.


title: On the variational interpretation of the discrete KP equation

abstract: We study the variational structure of the discrete Kadomtsev-Petviashvili (dKP) equation by means of its pluri-Lagrangian formulation. We consider the dKP equation and its variational formulation on the cubic lattice ${\mathbb Z}^{N}$ as well as on the root lattice $Q(A_{N})$. We prove that, on a lattice of dimension at least four, the corresponding Euler-Lagrange equations are equivalent to the dKP equation.


title: Milgrom's Law and Lambda's Shadow: How Massive Gravity Connects Galactic and Cosmic Dynamics

abstract: Massive gravity provides a natural solution for the dark energy problem of cosmology and is also a candidate for resolving the dark matter problem. I demonstrate that, assuming reasonable scaling relations, massive gravity can provide for Milgrom's law of gravity (or "modified Newtonian dynamics") which is known to remove the need for particle dark matter from galactic dynamics. Milgrom's law comes with a characteristic acceleration, Milgrom's constant, which is observationally constrained to $a_0\approx1.1\times10^{-10}$ m/s$^2$. In the derivation presented here, this constant arises naturally from the cosmologically required mass of gravitons like $a_0\propto c\sqrt{\Lambda}\propto c H_0\sqrt{3\Omega_{\Lambda}}$, with $\Lambda$, $H_0$, and $\Omega_{\Lambda}$ being the cosmological constant, the Hubble constant, and the third cosmological parameter, respectively. My derivation suggests that massive gravity could be the mechanism behind both, dark matter and dark energy.


title: Trajectory Optimization Algorithm Studies

abstract: In complex engineered systems, completing an objective is sometimes not enough. The system must be able to reach a set performance characteristic, such as an unmanned aerial vehicle flying from point A to point B, \textit{under 10 seconds}. This introduces the notion of optimality, what is the most efficient, the fastest, the cheapest way to complete a task. This report explores the two pillars of optimal control, Bellman's Dynamic Programming and Pontryagin's Maximum Principle, and compares implementations of both theories onto simulated systems. Dynamic Programming is realized through a Differential Dynamic Programming Algorithm, where utilizes a forward-backward pass to iteratively optimize a control sequence and trajectory. The Maximum Principle is realized via Gauss Pseudospectral Optimal Control, where the optimal control problem is first approximated through polynomial basis functions, then solved, with optimality being achieved through the costate equations of the Maximum Principle. The results of the report show that, for short time Horizons, DDP can optimize quickly and can generate a trajectory that utilizes less control effort for the same problem formulation. On the other hand Pseudospectral methods can optimize faster for longer time horizons, but require a key understanding of the problem structure. Future work involves completing an implementation of DDP in a C++ code, and testing the speed of convergence for both methods, as well as extended the Pseudospectral Optimal Control framework in to the world of stochastic optimal control.


title: A characterization of nilpotent nonassociative algebras by invertible Leibniz-derivations

abstract: Moens proved that a finite-dimensional Lie algebra over field of characteristic zero is nilpotent if and only if it has an invertible Leibniz-derivation. In this article we prove the analogous results for finite-dimensional Malcev, Jordan, (-1,1)-, quasiassociative, quasialternative, right alternative and Malcev-admissible noncommutative Jordan algebras over the field of characteristic zero. Also, we describe all Leibniz-derivations of semisimple Jordan, right alternative and Malcev algebras.


title: A Remark on Sieving in Biased Coin Convolutions

abstract: In this work, we establish a nontrivial level of distribution for densities on $\{1,\ldots, N\}$ obtained by a biased coin convolution. As a consequence of sieving theory, one then derives the expected lower bound for the weight of such densities on sets of pseudo-primes.


title: Generalized derivations of (color) n-ary algebras

abstract: We generalize the results of Leger and Luks about generalized derivations of Lie algebras to the case of color $n$-ary $\Omega$-algebras. Particularly, we prove some properties of generalized derivations of color $n$-ary algebras; prove that a quasiderivation algebra of a color $n$-ary $\Omega$-algebra can be embedded into the derivation algebra of a larger color $n$-ary $\Omega$-algebra, and describe (anti)commutative $n$-ary algebras satisfying the condition $QDer = End.$


title: Constraints on neutrino masses coming from magnetic dipole moments in a two Higgs doublet model type I and II

abstract: In the framework of a two Higgs doublet model type I and type II, we calculate limits on neutrino masses for the different types of neutrinos, by using the experimental bounds on their magnetic dipole moments. This is carried out by analyzing diagrams of Cherenkov neutrino decays with a charged Higgs into the loop, coming from the two Higgs doublet model (2HDM). Such constraints are translated into allowed regions in the free parameters of the models, for each neutrino flavor. The analysis was performed by sweeping the charged Higgs mass between (100-900)GeV and taking into account the experimental constraints for $\tan\beta $ in the 2HDM type I and II, obtaining contributions close to the experimental thresholds for muon and tau neutrinos, while for electron neutrino the relevant contribution comes from standard model and keeps out of the reach of forthcoming experiments.


title: On The Kantor Product

abstract: We study the algebra of bilinear multiplications of an $n$-dimensional vector space. In particular, we study the Kantor product of some well-known (associative, Lie, alternative, Novikov and some other) multiplications.


title: Remarks on an operator Wielandt inequality

abstract: Let $A$ be a positive operator on a Hilbert space $\mathcal{H}$ with $0<m\leq A\leq M$ and $X$ and $Y$ are two isometries on $\mathcal{H}$ such that $X^{*}Y=0$. For every 2-positive linear map $\Phi$, define $$\Gamma=\left(\Phi(X^{*}AY)\Phi(Y^{*}AY)^{-1}\Phi(Y^{*}AX)\right)^{p}\Phi(X^{*}AX)^{-p}, \, \, \, p>0.$$ We consider several upper bounds for $\frac{1}{2}|\Gamma+\Gamma^{*}|$. These bounds complement a recent result on operator Wielandt inequality.