title: Spectral multiplexing using stacked VPHGs - Part I

abstract: Many focal-reducer spectrographs, currently available at state-of-the art telescopes facilities, would benefit from a simple refurbishing that could increase both the resolution and spectral range in order to cope with the progressively challenging scientific requirements but, in order to make this update appealing, it should minimize the changes in the existing structure of the instrument. In the past, many authors proposed solutions based on stacking subsequently layers of dispersive elements and record multiple spectra in one shot (multiplexing). Although this idea is promising, it brings several drawbacks and complexities that prevent the straightforward integration of a such device in a spectrograph. Fortunately nowadays, the situation has changed dramatically thanks to the successful experience achieved through photopolymeric holographic films, used to fabricate common Volume Phase Holographic Gratings (VPHGs). Thanks to the various advantages made available by these materials in this context, we propose an innovative solution to design a stacked multiplexed VPHGs that is able to secure efficiently different spectra in a single shot. This allows to increase resolution and spectral range enabling astronomers to greatly economize their awarded time at the telescope. In this paper, we demonstrate the applicability of our solution, both in terms of expected performance and feasibility, supposing the upgrade of the Gran Telescopio CANARIAS (GTC) Optical System for Imaging and low-Intermediate-Resolution Integrated Spectroscopy (OSIRIS).


title: Virtual Resonant Emission and Oscillatory Long-Range Tails in van der Waals Interactions of Excited States: QED Treatment and Applications

abstract: We report on a quantum electrodynamic (QED) investigation of the interaction between a ground state atom with another atom in an excited state. General expressions, applicable to any atom, are indicated for the long-range tails which are due to virtual resonant emission and absorption into and from vacuum modes whose frequency equals the transition frequency to available lower-lying atomic states. For identical atoms, one of which is in an excited state, we also discuss the mixing term which depends on the symmetry of the two-atom wave function (these evolve into either the gerade or the ungerade state for close approach), and we include all nonresonant states in our rigorous QED treatment. In order to illustrate the findings, we analyze the fine-structure resolved van der Waals interaction for nD-1S hydrogen interactions with n=8,10,12 and find surprisingly large numerical coefficients.


title: Coverage and Rate Analysis of Super Wi-Fi Networks Using Stochastic Geometry

abstract: Recent regulatory changes proposed by the Federal Communications Commission (FCC) permitting unlicensed use of television white space (TVWS) channels present new opportunities for designing wireless networks that make efficient use of this spectrum. The favorable propagation characteristics of these channels and their widespread availability, especially in rural areas, make them well-suited for providing broadband services in sparsely populated regions where economic factors hinder deployment of such services on licensed spectrum. In this context, this paper explores the deployment of an outdoor Wi-Fi-like network operating in TVWS channels, referred to commonly as a Super Wi-Fi network. Since regulations governing unlicensed use of these channels allow (a) mounting fixed devices up to a height of 30 m and operation at transmit powers of up to 4 W EIRP, and (b) operation at transmit powers of up to 100 mW EIRP for portable devices, such networks can provide extended coverage and higher rates than traditional Wi-Fi networks. However, these gains are subject to the viability of the uplink from the portable devices (clients) to the fixed devices (access points (AP)) because of tighter restrictions on transmit power of clients compared to APs. This paper leverages concepts from stochastic geometry to study the performance of such networks with specific focus on the effect of (a) transmit power asymmetry between APs and clients and its impact on uplink viability and coverage, and (b) the interplay between height and transmit power of APs in determining the network throughput. Such an analysis reveals that (a) maximum coverage of no more than 700 m is obtained even when APs are deployed at 30 m height, and (b) operating APs at transmit power of more than 1 W is beneficial only at sparse deployment densities when rate is prioritized over coverage.


title: Two Mott Insulator Theory of Superconductivity in K$_3$X (X: picene, .. p-terphenyl, .. C$_{60}$)

abstract: We look for unifying aspects behind superconductivity in aromatic hydrocarbon and fullerene family K$_3$X (X: picene, .. p-terphenyl, .. C$_{60}$). Aromatic hydrocarbon molecules support RVB states. Consequent stability (aromaticity) makes them reluctant electron acceptors. We argue that X accepts only two (not all three) electrons from K$_3$ and creates charged RVB's in X$^{2-}$, and becomes a (molecular) Cooper pair box. A weak Josephson coupling between X$^{2-}$ molecules creates a Bose Mott insulator, a potential high Tc superconductor. Remaining lone electron in the complex (K$_3)^{2+}$ occupies a suitable metal orbital hybrid. They hybridize weakly through X$^{2-}$ molecular bridges, to form a half filled band of renormalized K atom orbitals, a Fermionic Mott insulator. An interplay of RVB physics and charge transfer (mutual doping) or external doping leads to superconductivity in one or both Mott insulators. In our theory there is room for room temperature superconductivity.


title: Structure of $^{20}$Ne states in the resonance $^{16}$O+$\alpha$ elastic scattering

abstract: Background The nuclear structure of the cluster bands in $^{20}$Ne presents a challenge for different theoretical approaches. It is especially difficult to explain the broad 0$^+$, 2$^+$ states at 9 MeV excitation energy. Simultaneously, it is important to obtain more reliable experimental data for these levels in order to quantitatively assess the theoretical framework. Purpose To obtain new data on $^{20}$Ne $\alpha$ cluster structure. Method Thick target inverse kinematics technique was used to study the $^{16}$O+$\alpha$ resonance elastic scattering and the data were analyzed using an \textit{R} matrix approach. The $^{20}$Ne spectrum, the cluster and nucleon spectroscopic factors were calculated using cluster-nucleon configuration interaction model (CNCIM). Results We determined the parameters of the broad resonances in \textsuperscript{20}Ne: 0$^+$ level at 8.77 $\pm$ 0.150 MeV with a width of 750 (+500/-220) keV; 2$^+$ level at 8.75 $\pm$ 0.100 MeV with the width of 695 $\pm$ 120 keV; the width of 9.48 MeV level of 65 $\pm$ 20 keV and showed that 9.19 MeV, 2$^+$ level (if exists) should have width $\leq$ 10 keV. The detailed comparison of the theoretical CNCIM predictions with the experimental data on cluster states was made. Conclusions Our experimental results by the TTIK method generally confirm the adopted data on $\alpha$ cluster levels in $^{20}$Ne. The CNCIM gives a good description of the $^{20}$Ne positive parity states up to an excitation energy of $\sim$ 7 MeV, predicting reasonably well the excitation energy of the states and their cluster and single particle properties. At higher excitations, the qualitative disagreement with the experimentally observed structure is evident, especially for broad resonances.


title: On sequents of $\Sigma$ formulas

abstract: We investigate the position that foundational theories should be modelled on ordinary computability. In this context, we investigate the metamathematics of $\Sigma$ formulas. We consider theories whose axioms are implications between $\Sigma$ formulas, and we show that arbitrarily strong such theories prove their own correctness. We also show that a natural extension of such a theory proves the validity of intuitionistic reasoning for that theory. Finally, we show the equivalence of two completeness principles appropriate to a potentialist conception of the universe of sets.


title: Locally optimal 2-periodic sphere packings

abstract: The sphere packing problem is an old puzzle. We consider packings with m spheres in the unit cell (m-periodic packings). For the case m = 1 (lattice packings), Voronoi proved there are finitely many inequivalent local optima and presented an algorithm to enumerate them, and this computation has been implemented in up to d = 8 dimensions. We generalize Voronoi's method to m > 1 and present a procedure to enumerate all locally optimal 2-periodic sphere packings in any dimension, provided there are finitely many. We implement this computation in d = 3, 4, and 5 and show that no 2-periodic packing surpasses the density of the optimal lattices in these dimensions. A partial enumeration is performed in d = 6.


title: R&D On Beam Injection and Bunching Schemes In The Fermilab Booster

abstract: Fermilab is committed to upgrade its accelerator complex to support HEP experiments at the intensity frontier. The ongoing Proton Improvement Plan (PIP) enables us to reach 700 kW beam power on the NuMI neutrino targets. By the end of the next decade, the current 400 MeV normal conducting LINAC will be replaced by an 800 MeV superconducting LINAC (PIP-II) with an increased beam power >50% of the PIP design goal. Both in PIP and PIP-II era, the existing Booster is going to play a very significant role, at least for next two decades. In the meanwhile, we have recently developed an innovative beam injection and bunching scheme for the Booster called "early injection scheme" that continues to use the existing 400 MeV LINAC and implemented into operation. This scheme has the potential to increase the Booster beam intensity by >40% from the PIP design goal. Some benefits from the scheme have already been seen. In this paper, I will describe the basic principle of the scheme, results from recent beam experiments, our experience with the new scheme in operation, current status, issues and future plans. This scheme fits well with the current and future intensity upgrade programs at Fermilab.


title: Anomaly-free local horizontal symmetry and anomaly-full rare B-decays

abstract: The largest global symmetry that can be made local in the Standard Model + 3$\nu_R$ while being compatible with Pati-Salam unification is $SU(3)_H\times U(1)_{B-L}$. The gauge bosons of this theory would induce flavour effects involving both quarks and leptons, and are a potential candidate to explain the recent reports of lepton universality violation in rare B meson decays. In this letter we characterise this type of models and show how they can accommodate the data and naturally be within reach of direct searches.


title: Proximity Eliashberg theory of electrostatic field-effect-doping in superconducting films

abstract: We calculate the effect of a static electric field on the critical temperature of a s-wave one band superconductor in the framework of proximity effect Eliashberg theory. In the weak electrostatic field limit the theory has no free parameters while, in general, the only free parameter is the thickness of the surface layer where the electric field acts. We conclude that the best situation for increasing the critical temperature is to have a very thin film of a superconducting material with a strong increase of electron-phonon (boson) constant upon charging.


title: From Fixed-X to Random-X Regression: Bias-Variance Decompositions, Covariance Penalties, and Prediction Error Estimation

abstract: In statistical prediction, classical approaches for model selection and model evaluation based on covariance penalties are still widely used. Most of the literature on this topic is based on what we call the "Fixed-X" assumption, where covariate values are assumed to be nonrandom. By contrast, it is often more reasonable to take a "Random-X" view, where the covariate values are independently drawn for both training and prediction. To study the applicability of covariance penalties in this setting, we propose a decomposition of Random-X prediction error in which the randomness in the covariates contributes to both the bias and variance components. This decomposition is general, but we concentrate on the fundamental case of least squares regression. We prove that in this setting the move from Fixed-X to Random-X prediction results in an increase in both bias and variance. When the covariates are normally distributed and the linear model is unbiased, all terms in this decomposition are explicitly computable, which yields an extension of Mallows' Cp that we call $RCp$. $RCp$ also holds asymptotically for certain classes of nonnormal covariates. When the noise variance is unknown, plugging in the usual unbiased estimate leads to an approach that we call $\hat{RCp}$, which is closely related to Sp (Tukey 1967), and GCV (Craven and Wahba 1978). For excess bias, we propose an estimate based on the "shortcut-formula" for ordinary cross-validation (OCV), resulting in an approach we call $RCp^+$. Theoretical arguments and numerical simulations suggest that $RCP^+$ is typically superior to OCV, though the difference is small. We further examine the Random-X error of other popular estimators. The surprising result we get for ridge regression is that, in the heavily-regularized regime, Random-X variance is smaller than Fixed-X variance, which can lead to smaller overall Random-X error.


title: Stability of zero-growth economics analysed with a Minskyan model

abstract: As humanity is becoming increasingly confronted by Earth's finite biophysical limits, there is increasing interest in questions about the stability and equitability of a zero-growth capitalist economy, most notably: if one maintains a positive interest rate for loans, can a zero-growth economy be stable? This question has been explored on a few different macroeconomic models, and both `yes' and `no' answers have been obtained. However, economies can become unstable whether or not there is ongoing underlying growth in productivity with which to sustain growth in output. Here we attempt, for the first time, to assess via a model the relative stability of growth versus no-growth scenarios. The model employed draws from Keen's model of the Minsky financial instability hypothesis. The analysis focuses on dynamics as opposed to equilibrium, and scenarios of growth and no-growth of output (GDP) are obtained by tweaking a productivity growth input parameter. We confirm that, with or without growth, there can be both stable and unstable scenarios. To maintain stability, firms must not change their debt levels or target debt levels too quickly. Further, according to the model, the wages share is higher for zero-growth scenarios, although there are more frequent substantial drops in employment.


title: Tripartite entanglement detection through tripartite quantum steering in one-sided and two-sided device-independent scenarios

abstract: In the present work, we study tripartite quantum steering of quantum correlations arising from two local dichotomic measurements on each side in the two types of partially device-independent scenarios: $1$-sided device-independent scenario where one of the parties performs untrusted measurements while the other two parties perform trusted measurements and $2$-sided device-independent scenario where one of the parties performs trusted measurements while the other two parties perform untrusted measurements. We demonstrate that tripartite steering in the $2$-sided device-independent scenario is weaker than tripartite steering in the $1$-sided device-independent scenario by using two families of quantum correlations. That is these two families of quantum correlations in the $2$-sided device-independent framework detect tripartite entanglement through tripartite steering for a larger region than that in the $1$-sided device-independent framework. It is shown that tripartite steering in the $2$-sided device-independent scenario implies the presence of genuine tripartite entanglement of $2\times 2 \times 2$ quantum system, even if the correlation does not exhibit genuine nonlocality or genuine steering.


title: Classical and Quantum Aspects of Particle Propagation in External Gravitational Fields

abstract: In the study of covariant wave equations, linear gravity manifests itself through the metric deviation $\gamma_{\mu\nu}$ and a two-point vector potential $K_{\lambda}$ itself constructed from $\gamma_{\mu\nu}$ and its derivatives. The simultaneous presence of the two gravitational potentials is non contradictory. Particles also assume the character of quasiparticles and $K_{\lambda}$ carries information about the matter with which it interacts. We consider the influence of $K_{\lambda}$ on the dispersion relations of the particles involved, the particles' motion, quantum tunneling through a horizon, radiation, energy-momentum dissipation and flux quantization. % No {\it REVTEX} limit to number of lines.


title: Spin Torque Efficiency and Analytic Error Rate Estimates of Skyrmion Racetrack Memory

abstract: In this paper the thermal stability of skyrmion bubbles and the critical currents to move them over pinning sites is investigated. For the used pinning geometries and the used parameters, the unexpected behavior is reported that the energy barrier to overcome the pinning site is larger than the energy barrier of the annihilation of a skyrmion. The annihilation takes place at boundaries by current driven motion as well as due to the excitation over energy barriers, in the absence of currents, without forming Bloch points. It is reported that the pinning sites, which are required to allow thermally stable bits, significantly increase the critical current densities to move the bits in skyrmion like structures to about $j_{crit}$ = 0.62 TA/m$^2$. These currents are similar to those obtained experimentally to move stable skyrmions at room temperature. By calculating the thermal stability as well as the critical current, we can derive the spin torque efficiency $\eta$ = $\Delta/I_c = 0.19 k_BT_{300}~/\mu$A, which is in a similar range to the simulated spin torque efficiency of MRAM structures. Finally, it is shown that the stochastic depinning process of any racetrack like device requires extremely narrow depinning time distribution smaller than ~6% of the current pulse length to reach bit error rates smaller than ${10^{ - 9}}$.


title: A Generalization of Convolutional Neural Networks to Graph-Structured Data

abstract: This paper introduces a generalization of Convolutional Neural Networks (CNNs) from low-dimensional grid data, such as images, to graph-structured data. We propose a novel spatial convolution utilizing a random walk to uncover the relations within the input, analogous to the way the standard convolution uses the spatial neighborhood of a pixel on the grid. The convolution has an intuitive interpretation, is efficient and scalable and can also be used on data with varying graph structure. Furthermore, this generalization can be applied to many standard regression or classification problems, by learning the the underlying graph. We empirically demonstrate the performance of the proposed CNN on MNIST, and challenge the state-of-the-art on Merck molecular activity data set.


title: The Spectrum of SS 433 in the H and K Bands

abstract: SS~433 is an X-ray binary and the source of sub-relativistic, precessing, baryonic jets. We present high-resolution spectrograms of SS 433 in the infrared H and K bands. The spectrum is dominated by hydrogen and helium emission lines. The precession phase of the emission lines from the jet continues to be described by a constant period, P_jet= 162.375 d. The limit on any secularly changing period is $|\dot P| \lesssim 10^{-5}$. The He I 2.0587 micron line has complex and variable P Cygni absorption features produced by an inhomogeneous wind with a maximum outflow velocity near 900 km/s. The He II emission lines in the spectrum also arise in this wind. The higher members of the hydrogen Brackett lines show a double-peaked profile with symmetric wings extending more than +/-1500 km/s from the line center. The lines display radial velocity variations in phase with the radial velocity variation expected of the compact star, and they show a distortion during disk eclipse that we interpret as a rotational distortion. We fit the line profiles with a model in which the emission comes from the surface of a symmetric, Keplerian accretion disk around the compact object. The outer edge of the disk has velocities that vary from 110 to 190 km/s. These comparatively low velocities place an important constraint on the mass of the compact star: Its mass must be less than 2.2 M_solar and is probably less than 1.6 M_solar.


title: Colombeau algebras without asymptotics

abstract: We present a construction of algebras of generalized functions of Colombeau-type which, instead of asymptotic estimates with respect to a regularization parameter, employs only topological estimates on certain spaces of kernels for its definition.


title: Implications of the $R_K$ and $R_{K^*}$ anomalies

abstract: We discuss the implications of the recently reported $R_K$ and $R_{K^*}$ anomalies, the lepton flavor non-universality in the $B\to K\ell^+\ell^-$ and $B\to K^*\ell^+\ell^-$. Using two sets of hadronic inputs of form factors, we perform a fit of the new physics to the $R_K$ and $R_{K^*}$ data, and significant new physics contributions are found. We suggest to study the lepton flavor universality in a number of related rare $B, B_s, B_c$ and $\Lambda_b$ decay channels, and in particular we give the predictions for the $\mu$-to-$e$ ratios of decay widths with different polarizations of the final state particles, and of the $b\to d\ell^+\ell^-$ processes which are presumably more sensitive to the structure of the underlying new physics. With the new physics contributions embedded in Wilson coefficients, we present theoretical predictions for lepton flavor non-universality in these processes.


title: Noisy entanglement-assisted classical capacity as a security framework for two-way quantum key distribution protocols

abstract: We use the noisy entanglement-assisted classical capacity formula [arXiv:1609.08592] to create a coherent-attack security framework for Gaussian two-way quantum key distribution protocols in the asymptotic region.


title: Performance Characterization of a Solenoid-type Gas Valve for the $H^{-}$ Magnetron Source at FNAL

abstract: The magnetron-style $H^{-}$ ion sources currently in operation at Fermilab use piezoelectric gas valves to function. This kind of gas valve is sensitive to small changes in ambient temperature, which affect the stability and performance of the ion source. This motivates the need to find an alternative way of feeding H2 gas into the source. A solenoid-type gas valve has been characterized in a dedicated off-line test stand to assess the feasibility of its use in the operational ion sources. $H^{-}$ ion beams have been extracted at 35 keV using this valve. In this study, the performance of the solenoid gas valve has been characterized measuring the beam current output of the magnetron source with respect to the voltage and pulse width of the signal applied to the gas valve.


title: Asymptotical properties of social network dynamics on time scales

abstract: In this paper we develop conditions for various types of stability in social networks governed by Imitation of Success principle. Considering so-called Prisoner's Dilemma as the base of node-to-node game in the network we obtain well-known Hopfield neural network model. Asymptotic behavior of the original model and dynamic Hopfield model has a certain correspondence. To obtain more general results, we consider Hopfield model dynamic system on time scales. Developed stability conditions combine main parameters of network structure such as network size and maximum relative nodes' degree with the main characteristics of time scale, nodes' inertia and resistance, rate of input-output response.


title: Comment on "Ultrathin Interface Regime of Core-Shell Magnetic Nanoparticles for Effective Magnetism Tailoring"

abstract: There has been renewed interest in the use of magnetic nanoparticles to convert high frequency electromagnetic energy into heat. During the last decade, numerous examples in the field of catalysis, lightweight thermoplastic composites for aeronautical and automotive engineering, and biomedical applications as well, have flourished. In a recent Letter "Ultrathin Interface Regime of Core-Shell Magnetic Nanoparticles for Effective Magnetism Tailoring", Moon et al. DOI: 10.1021/acs.nanolett.6b04016 claim to have observed a new interface spin canting phenomenon in core-shell ferrrites that pushes heating efficiencies above 10,000 W/g. If this were true, it represents new opportunities for further tailoring nanoscale agents to be used in drug delivery and selective destruction of tumours by hyperthermia. But here I argue that these claims are inconsistent with the facts of physics.


title: Scalar products of Bethe vectors in the models with $\mathfrak{gl}(m|n)$ symmetry

abstract: We study scalar products of Bethe vectors in the models solvable by the nested algebraic Bethe ansatz and described by $\mathfrak{gl}(m|n)$ superalgebra. Using coproduct properties of the Bethe vectors we obtain a sum formula for their scalar products. This formula describes the scalar product in terms of a sum over partitions of Bethe parameters. We also obtain recursions for the Bethe vectors. This allows us to find recursions for the highest coefficient of the scalar product.


title: Structures far below sub-Planck scale in quantum phase-space through superoscillations

abstract: In 2001, Zurek derived the generic minimum scale $a_{Z}$ for the area of structures of Wigner's quantum phase distribution. Here we show by construction, using superoscillatory functions, that the Wigner distribution can locally show regular spotty structures on scales much below Zurek's scale $a_{Z}$. The price to pay for the presence of such structures is their exponential smallness. For the case we construct there is no increased interferometric sensitivity from the presence of patches with superoscillatory structure in phase-space.


title: High-Frequency Jump Analysis of the Bitcoin Market

abstract: We use the database leak of Mt. Gox exchange to analyze the dynamics of the price of bitcoin from June 2011 to November 2013. This gives us a rare opportunity to study an emerging retail-focused, highly speculative and unregulated market with trader identifiers at a tick transaction level. Jumps are frequent events and they cluster in time. The order flow imbalance and the preponderance of aggressive traders, as well as a widening of the bid-ask spread predict them. Jumps have short-term positive impact on market activity and illiquidity and see a persistent change in the price.


title: Ring Oscillators for Clocking Reversible Superconducting Circuits and Fabrication Process Benchmarking

abstract: Existing concepts of reversible superconducting circuits as well as demonstrated adiabatic circuits require three-phase bias/clock signals generated by room temperature sources. A while ago, we suggested that a multi-phase bias/clock could be provided by a local Josephson junction-based generator. The generator looks like a long annular Josephson junction, only composed of discreet elements - junctions and inductors, and closed into a ring via a flux pump to inject any required number of vortices into the ring. A steady motion of the vortices forced by a uniformly distributed dc bias current applied to the ring is accompanied by a nearly harmonic ac currents flowing via the Josephson junctions (JJs) connected in series with small inductors. These ac currents serve as multi-phase bias/clock for nSQUID-based circuitry. To verify this concept and compare the dissipated energy with kBTln2 threshold, we developed a ring composed of 256 unshunted JJs with 20 {\mu}A target critical current, Ic. We investigated the behavior of the ring oscillator at each vortex count from 0 to 256. The measured critical current of the ring with vortices was about 0.1 {\mu}A per one JJ, which can be explained by unavoidable nonuniformity of the ring components and the influence of fluxes frozen near the ring. The corresponding energy dissipation, about 10kBT per passage of one vortex through one JJ, should be reduced further for prospective experiments with reversible circuits. However, obtained I-V characteristics could be of interest for scientists working with long Josephson junctions. Superiority of the fabrication process used in this work is demonstrated by the obtained about 200 times reduction of Ic of the ring with vortices with respect to a single comprising JJ, much larger than in any previously described case.


title: Recent Operation of the FNAL Magnetron $H^{-}$ Ion Source

abstract: This paper will detail changes in the operational paradigm of the Fermi National Accelerator Laboratory (FNAL) magnetron $H^{-}$ ion source due to upgrades in the accelerator system. Prior to November of 2012 the $H^{-}$ ions for High Energy Physics (HEP) experiments were extracted at ~18 keV vertically downward into a 90 degree bending magnet and accelerated through a Cockcroft-Walton accelerating column to 750 keV. Following the upgrade in the fall of 2012 the $H^{-}$ ions are now directly extracted from a magnetron at 35 keV and accelerated to 750 keV by a Radio Frequency Quadrupole (RFQ). This change in extraction energy as well as the orientation of the ion source required not only a redesign of the ion source, but an updated understanding of its operation at these new values. Discussed in detail are the changes to the ion source timing, arc discharge current, hydrogen gas pressure, and cesium delivery system that were needed to maintain consistent operation at >99% uptime for HEP, with an increased ion source lifetime of over 9 months.


title: MASTER OT J004207.99+405501.1/M31LRN 2015 Luminous Red Nova in M31: Discovery, Light Curve, Hydrodynamics, Evolution

abstract: We report the discovery and multicolor (VRIW) photometry of a rare explosive star MASTER OT J004207.99+405501.1 - a luminous red nova - in the Andromeda galaxy M31N2015-01a. We use our original light curve acquired with identical MASTER Global Robotic Net telescopes in one photometric system: VRI during first 30 days and W (unfiltered) during 70 days. Also we added publishied multicolor photometry data to estimate the mass and energy of the ejected shell, and discuss the likely formation scenarios of outbursts of this type. We propose the interpretation of the explosion, that is consistent with the evolutionary scenario where star merger is a natural stage of the evolution of close-mass stars and may serve as an extra channel for the formation of nova outbursts.


title: GRB 081024B and GRB 140402A: two additional short GRBs from binary neutron star mergers

abstract: Theoretical and observational evidences have been recently gained for a two-fold classification of short bursts: 1) short gamma-ray flashes (S-GRFs), with isotropic energy $E_{iso}<10^{52}$~erg and no BH formation, and 2) the authentic short gamma-ray bursts (S-GRBs), with isotropic energy $E_{iso}>10^{52}$~erg evidencing a BH formation in the binary neutron star merging process. The signature for the BH formation consists in the on-set of the high energy ($0.1$--$100$~GeV) emission, coeval to the prompt emission, in all S-GRBs. No GeV emission is expected nor observed in the S-GRFs. In this paper we present two additional S-GRBs, GRB 081024B and GRB 140402A, following the already identified S-GRBs, i.e., GRB 090227B, GRB 090510 and GRB 140619B. We also return on the absence of the GeV emission of the S-GRB 090227B, at an angle of $71^{\rm{o}}$ from the \textit{Fermi}-LAT boresight. All the correctly identified S-GRBs correlate to the high energy emission, implying no significant presence of beaming in the GeV emission. The existence of a common power-law behavior in the GeV luminosities, following the BH formation, when measured in the source rest-frame, points to a commonality in the mass and spin of the newly-formed BH in all S-GRBs.


title: Entanglement generation between a charge qubit and its bosonic environment during pure dephasing - dependence on environment size

abstract: We study entanglement generated between a charge qubit and a bosonic bath due to their joint evolution which leads to pure dephasing of the qubit. We tune the parameters of the interaction, so that the decoherence is quantitatively independent of the number of bosonic modes taken into account and investigate, how the entanglement generated depends on the size of the environment. A second parameter of interest is the mixedness of the initial state of the environment which is controlled by temperature. We show analytically that for a pure initial state of the environment, entanglement does not depend on environment size. For mixed initial states of the environment, the generated entanglement decreases with the increase of environment size. This effect is stronger for larger temperatures, when the environment is initially more mixed, but in the limit of an infinitely large environment, no entanglement is created at any finite temperature.


title: Roadmap for the international, accelerator-based neutrino programme

abstract: In line with its terms of reference the ICFA Neutrino Panel has developed a roadmapfor the international, accelerator-based neutrino programme. A "roadmap discussion document" was presented in May 2016 taking into account the peer-group-consultation described in the Panel's initial report. The "roadmap discussion document" was used to solicit feedback from the neutrino community---and more broadly, the particle- and astroparticle-physics communities---and the various stakeholders in the programme. The roadmap, the conclusions and recommendations presented in this document take into account the comments received following the publication of the roadmap discussion document. With its roadmap the Panel documents the approved objectives and milestones of the experiments that are presently in operation or under construction. Approval, construction and exploitation milestones are presented for experiments that are being considered for approval. The timetable proposed by the proponents is presented for experiments that are not yet being considered formally for approval. Based on this information, the evolution of the precision with which the critical parameters governinger the neutrino are known has been evaluated. Branch or decision points have been identified based on the anticipated evolution in precision. The branch or decision points have in turn been used to identify desirable timelines for the neutrino-nucleus cross section and hadro-production measurements that are required to maximise the integrated scientific output of the programme. The branch points have also been used to identify the timeline for the R&D required to take the programme beyond the horizon of the next generation of experiments. The theory and phenomenology programme, including nuclear theory, required to ensure that maximum benefit is derived from the experimental programme is also discussed.


title: Provably secure key establishment against quantum adversaries

abstract: At Crypto 2011, some of us had proposed a family of cryptographic protocols for key establishment capable of protecting quantum and classical legitimate parties unconditionally against a quantum eavesdropper in the query complexity model. Unfortunately, our security proofs were unsatisfactory from a cryptographically meaningful perspective because they were sound only in a worst-case scenario. Here, we extend our results and prove that for any e > 0, there is a classical protocol that allows the legitimate parties to establish a common key after O(N) expected queries to a random oracle, yet any quantum eavesdropper will have a vanishing probability of learning their key after O(N^{1.5-e}) queries to the same oracle. The vanishing probability applies to a typical run of the protocol. If we allow the legitimate parties to use a quantum computer as well, their advantage over the quantum eavesdropper becomes arbitrarily close to the quadratic advantage that classical legitimate parties enjoyed over classical eavesdroppers in the seminal 1974 work of Ralph Merkle. Along the way, we develop new tools to give lower bounds on the number of quantum queries required to distinguish two probability distributions. This method in itself could have multiple applications in cryptography. We use it here to study average-case quantum query complexity, for which we develop a new composition theorem of independent interest.


title: A Dunkl Analogue of Operators Including Two-variable Hermite polynomials

abstract: The aim of this paper is to introduce a Dunkl generalization of the operators including two variable Hermite polynomials which are defined by Krech [14](Krech, G. A note on some positive linear operators associated with the Hermite polynomials, Carpathian J. Math., 32 (1) (2016), 71--77) and to investigate approximating properties for these operators by means of the classical modulus of continuity, second modulus of continuity and Peetre's K-functional.


title: Elementary considerations for classes of meromorphic univalent functions

abstract: In this article we consider functions $f$ meromorphic in the unit disk. We give an elementary proof for a condition that is sufficient for the univalence of such functions. This condition simplifies and generalizes known conditions. We present some typical problems of geometrical function theory and give elementary solutions in the case of the above functions.


title: s+is superconductivity with incipient bands: doping dependence and STM signatures

abstract: Motivated by the recent observations of small Fermi energies and comparatively large superconducting gaps, present also on bands not crossing the Fermi energy (incipient bands) in iron-based superconductors, we analyse the doping evolution of superconductivity in a four-band model across the Lifshitz transition including BCS-BEC crossover effects on the shallow bands. Similar to the BCS case we find that with hole doping the phase difference between superconducting order parameters of the hole bands changes from $0$ to $\pi$ through an intermediate $s+is$ state breaking time-reversal symmetry. The transition however occurs in the region where electron bands are incipient and chemical potential renormalization in the superconducting state leads to a significant broadening of the $s+is$ region. We further present the qualitative features of the $s+is$ state that can be observed in scanning tunnelling microscopy (STM) experiments also taking incipient bands into account.


title: On the local $M$-derivative

abstract: We introduce a new fractional derivative that generalizes the so-called alternative fractional derivative recently proposed by Katugampola. We denote this new differential operator by $\mathscr{D}_{M}^{\alpha,\beta }$, where the parameter $\alpha$, associated with the order, is such that $0<\alpha<1$, $\beta>0$ and $M$ is used to denote that the function to be derived involves a Mittag-Leffler function with one parameter. This new derivative satisfies some properties of integer-order calculus, e.g.\ linearity, product rule, quotient rule, function composition and the chain rule. Besides as in the case of the Caputo derivative, the derivative of a constant is zero. Because Mittag-Leffler function is a natural generalization of the exponential function, we can extend some of the classical results of integer-order calculus, namely: Rolle's theorem, the mean value theorem and its extension. Further, when the order of the derivative is $\alpha=1$ and the parameter of the Mittag-Leffler function is also unitary, our definition is equivalent to the definition of the ordinary derivative of order one. Finally, we present the corresponding fractional integral from which, as a natural consequence, new results emerge which can be interpreted as applications. Specifically, we generalize the inversion property of the fundamental theorem of calculus and prove a theorem associated with the classical integration by parts.


title: A new truncated $M$-fractional derivative type unifying some fractional derivative types with classical properties

abstract: We introduce a truncated $M$-fractional derivative type for $\alpha$-differentiable functions that generalizes four other fractional derivatives types recently introduced by Khalil et al., Katugampola and Sousa et al., the so-called conformable fractional derivative, alternative fractional derivative, generalized alternative fractional derivative and $M$-fractional derivative, respectively. We denote this new differential operator by $_{i}\mathscr{D}_{M}^{\alpha,\beta }$, where the parameter $\alpha$, associated with the order of the derivative is such that $ 0 <\alpha<1 $, $\beta>0$ and $ M $ is the notation to designate that the function to be derived involves the truncated Mittag-Leffler function with one parameter. The definition of this truncated $M$-fractional derivative type satisfies the properties of the integer-order calculus. We also present, the respective fractional integral from which emerges, as a natural consequence, the result, which can be interpreted as an inverse property. Finally, we obtain the analytical solution of the $M$-fractional heat equation and present a graphical analysis.


title: $f(R)$ constant-roll inflation

abstract: The previously introduced class of two-parametric phenomenological inflationary models in General Relativity in which the slow-roll assumption is replaced by the more general, constant-roll condition is generalized to the case of $f(R)$ gravity. A simple constant-roll condition is defined in the original Jordan frame, and exact expressions for a scalaron potential in the Einstein frame, for a function $f(R)$ (in the parametric form) and for inflationary dynamics are obtained. The region of the model parameters permitted by the latest observational constraints on the scalar spectral index and the tensor-to-scalar ratio of primordial metric perturbations generated during inflation is determined.


title: Properties of Ultra Gamma Function

abstract: In this paper we study the integral of type \[_{\delta,a}\Gamma_{\rho,b}(x) =\Gamma(\delta,a;\rho,b)(x)=\int_{0}^{\infty}t^{x-1}e^{-\frac{t^{\delta}}{a}-\frac{t^{-\rho}}{b}}dt.\] Different authors called this integral by different names like ultra gamma function, generalized gamma function, Kratzel integral, inverse Gaussian integral, reaction-rate probability integral, Bessel integral etc. We prove several identities and recurrence relation of above said integral, we called this integral as Four Parameter Gamma Function. Also we evaluate relation between Four Parameter Gamma Function, p-k Gamma Function and Classical Gamma Function. With some conditions we can evaluate Four Parameter Gamma Function in term of Hypergeometric function.


title: Generalized Convex Functions and Their Applications

abstract: This study focuses on convex functions and their generalized. Thus, we start this study by giving the definition of convex functions and some of their properties and discussing a simple geometric property. Then we generalize E-convex functions and establish some their properties. Moreover, we give generalized $ s $-convex functions in the second sense and present some new inequalities of generalized Hermite-Hadamard type for the class of functions whose second local fractional derivatives of order $ \alpha $ in absolute value at certain powers are generalized $ s $-convex functions in the second sense. At the end, some examples that these inequalities are able to be applied to some special means are showed.


title: Note on extensions of the beta function

abstract: The classical beta function B(x; y) is one of the most fundamental special functions, due to its important role in various fields in the mathematical, physical, engineering and statistical sciences. Useful extensions of the classical Beta function has been considered by many authors. In the present paper, our main objective is to study the convergence of extensions of classical beta function and introduce modified extension of classical beta function. It is interpreted numerically and geometrically in the view of convergence, further properties and integral presentations are established.


title: Missing Data and Prediction

abstract: Missing data are a common problem for both the construction and implementation of a prediction algorithm. Pattern mixture kernel submodels (PMKS) - a series of submodels for every missing data pattern that are fit using only data from that pattern - are a computationally efficient remedy for both stages. Here we show that PMKS yield the most predictive algorithm among all standard missing data strategies. Specifically, we show that the expected loss of a forecasting algorithm is minimized when each pattern-specific loss is minimized. Simulations and a re-analysis of the SUPPORT study confirms that PMKS generally outperforms zero-imputation, mean-imputation, complete-case analysis, complete-case submodels, and even multiple imputation (MI). The degree of improvement is highly dependent on the missingness mechanism and the effect size of missing predictors. When the data are Missing at Random (MAR) MI can yield comparable forecasting performance but generally requires a larger computational cost. We see that predictions from the PMKS are equivalent to the limiting predictions for a MI procedure that uses a mean model dependent on missingness indicators (the MIMI model). Consequently, the MIMI model can be used to assess the MAR assumption in practice. The focus of this paper is on out-of-sample prediction behavior, implications for model inference are only briefly explored.


title: Evaluating openEHR for storing computable representations of electronic health record phenotyping algorithms

abstract: Electronic Health Records (EHR) are data generated during routine clinical care. EHR offer researchers unprecedented phenotypic breadth and depth and have the potential to accelerate the pace of precision medicine at scale. A main EHR use-case is creating phenotyping algorithms to define disease status, onset and severity. Currently, no common machine-readable standard exists for defining phenotyping algorithms which often are stored in human-readable formats. As a result, the translation of algorithms to implementation code is challenging and sharing across the scientific community is problematic. In this paper, we evaluate openEHR, a formal EHR data specification, for computable representations of EHR phenotyping algorithms.


title: A formula for the nonsymmetric Opdam's hypergeometric function of type $A_2$

abstract: The aim of this paper is to give an explicit formula for the nonsymmetric Heckman-Opdam's hypergeometric function of type $A_2$. This is obtained by differentiating the corresponding symmetric hypergeometric function.


title: Moving-centre monotonicity formulae for minimal submanifolds and related equations

abstract: Monotonicity formulae play a crucial role for many geometric PDEs, especially for their regularity theories. For minimal submanifolds in a Euclidean ball, the classical monotonicity formula implies that if such a submanifold passes through the centre of the ball, then its area is at least that of the equatorial disk. Recently Brendle and Hung proved a sharp area bound for minimal submanifolds when the prescribed point is not the centre of the ball, which resolved a conjecture of Alexander, Hoffman and Osserman. Their proof involves asymptotic analysis of an ingeniously chosen vector field, and the divergence theorem. In this article we prove a sharp `moving-centre' monotonicity formula for minimal submanifolds, which implies the aforementioned area bound. We also describe similar moving-centre monotonicity formulae for stationary $p$-harmonic maps, mean curvature flow and the harmonic map heat flow.


title: Diversity of knot solitons in liquid crystals manifested by linking of preimages in torons and hopfions

abstract: Topological solitons are knots in continuous physical fields classified by non-zero Hopf index values. Despite arising in theories that span many branches of physics, from elementary particles to condensed matter and cosmology, they remain experimentally elusive and poorly understood. We introduce a method of experimental and numerical analysis of such localized structures in liquid crystals that, similar to the mathematical Hopf maps, relates all points of the medium's order parameter space to their closed-loop preimages within the three-dimensional solitons. We uncover a surprisingly large diversity of naturally occurring and laser-generated topologically nontrivial solitons with differently knotted nematic fields, which previously have not been realized in theories and experiments alike. We discuss the implications of the liquid crystal's non-polar nature on the knot soliton topology and how the medium's chirality, confinement and elastic anisotropy help to overcome the constrains of the Hobart-Derrick theorem, yielding static three-dimensional solitons without or with additional defects. Our findings will establish chiral nematics as a model system for experimental exploration of topological solitons and may impinge on understanding of such nonsingular field configurations in other branches of physics, as well as may lead to technological applications


title: Character Networks and Book Genre Classification

abstract: We compare the social character networks of biographical, legendary and fictional texts, in search for marks of genre differentiation. We examine the degree distribution of character appearance and find a power law that does not depend on the literary genre or historical content. We also analyze local and global complex networks measures, in particular, correlation plots between the recently introduced Lobby (or Hirsh $H(1)$) index and Degree, Betweenness and Closeness centralities. Assortativity plots, which previous literature claims to separate fictional from real social networks, were also studied. We've found no relevant differences in the books for these network measures and we give a plausible explanation why the previous assortativity result is not correct.


title: Planck scale Corrections to the Harmonic Oscillator, Coherent and Squeezed States

abstract: The Generalized Uncertainty Principle (GUP) is a modification of Heisenberg's Principle predicted by several theories of Quantum Gravity. It consists of a modified commutator between position and momentum. In this work we compute potentially observable effects that GUP implies for the harmonic oscillator, coherent and squeezed states in Quantum Mechanics. In particular, we rigorously analyze the GUP-perturbed harmonic oscillator Hamiltonian, defining new operators that act as ladder operators on the perturbed states. We use these operators to define the new coherent and squeezed states. We comment on potential applications.


title: Perpetual integrals convergence and extinctions in population dynamics

abstract: In this article we use a criterion for the integrability of paths of one-dimensional diffusion processes from which we derive new insights on allelic fixation in several situations. This well known criterion involves a simple necessary and sufficient condition based on scale function and speed measure. We provide a new simple proof for this result and also obtain explicit bounds for the moments of such integrals. We also extend this criterion to non-homogeneous processes by use of Girsanov's transform. We apply our results to multi-type population dynamics: using the criterion with appropriate time changes, we characterize the behavior of proportions of each type before population extinction in different situations.


title: Quadratically-Regularized Optimal Transport on Graphs

abstract: Optimal transportation provides a means of lifting distances between points on a geometric domain to distances between signals over the domain, expressed as probability distributions. On a graph, transportation problems can be used to express challenging tasks involving matching supply to demand with minimal shipment expense; in discrete language, these become minimum-cost network flow problems. Regularization typically is needed to ensure uniqueness for the linear ground distance case and to improve optimization convergence; state-of-the-art techniques employ entropic regularization on the transportation matrix. In this paper, we explore a quadratic alternative to entropic regularization for transport over a graph. We theoretically analyze the behavior of quadratically-regularized graph transport, characterizing how regularization affects the structure of flows in the regime of small but nonzero regularization. We further exploit elegant second-order structure in the dual of this problem to derive an easily-implemented Newton-type optimization algorithm.


title: Transient instabilities in swelling dynamics

abstract: We investigate the swelling dynamics driven by solvent absorption in a hydrogel sphere immersed in a solvent bath, through an accurate computational model and numerical study. We extensively describe the transient process from dry to wet and discuss the onset of surface instabilities through a measure of the lack of smoothness of the outer surface and a morphological pattern of that surface with respect to the two material parameters driving the swelling dynamics.


title: Compressed sensing in the quaternion algebra

abstract: The article concerns compressed sensing methods in the quaternion algebra. We prove that it is possible to uniquely reconstruct - by $\ell_1$-norm minimization - a sparse quaternion signal from a limited number of its linear measurements, provided the quaternion measurement matrix satisfies so-called restricted isometry property with a sufficiently small constant. We also provide error estimates for the reconstruction of a non-sparse quaternion signal in the noisy and noiseless cases.


title: Uniqueness of minimal energy solutions for a semilinear problem involving the fractional laplacian

abstract: In this paper we study a semilinear problem for the fractional laplacian that are the counterpart of the Neumann problems in the classical setting. We show uniqueness of minimal energy solutions for small domains.


title: Deterministic creation of a four-qubit W state using one- and two-qubit gates

abstract: We propose an optical scheme to build an entangled network composed of $W$ state based on polarization encoded qubits (photons). This new setup consists of 2 $cNOT$ gates, 4 $V$ gates, 2 $Hadamard$ gates and basic optical tools such as polarizing beamsplitters ($PBS$s) and path couplers ($PC$s). $V$ gate is a specially-designed tool acting as a two-qubit gate which is composed of a $cNOT$ gate, 3 $PBS$s and a $PC$. By using this gate, one benefits from the temporarily generated optical degree of freedom, which is the spatial mode of a photon in the proposed scheme. Using an extra degree of freedom allows us to perform more capable processing for $W$-state creation protocols. We use four photons as input, which means we do not need entanglement as a resource. Also, we show that the proposed scheme can be implemented by operating the quantum optical gates which can be realized via current photonics technology.


title: On 2-Verma modules for quantum $\mathfrak{sl}_2$

abstract: In this paper we study the superalgebra $A_n$, introduced by the authors in previous work on categorification of Verma modules for quantum $\mathfrak{sl}_2$. The superalgebra $A_n$ is akin to the nilHecke algebra, and shares similar properties. In particular, we prove a uniqueness result about 2-Verma modules on $\Bbbk$-linear 2-categories.


title: Renormalization group procedure for potential $-g/r^2$

abstract: Schr\"odinger equation with potential $-g/r^2$ exhibits a limit cycle, described in the literature in a broad range of contexts using various regularizations of the singularity at $r=0$. Instead, we use the renormalization group transformation based on Gaussian elimination, from the Hamiltonian eigenvalue problem, of high momentum modes above a finite, floating cutoff scale. The procedure identifies a richer structure than the one we found in the literature. Namely, it directly yields an equation that determines the renormalized Hamiltonians as functions of the floating cutoff: solutions to this equation exhibit, in addition to the limit-cycle, also the asymptotic-freedom, triviality, and fixed-point behaviors, the latter in vicinity of infinitely many separate pairs of fixed points in different partial waves for different values of $g$.


title: LHC as an Axion Factory: Probing an Axion Explanation for $(g-2)_\mu$ with Exotic Higgs Decays

abstract: We argue that a large region of so far unconstrained parameter space for axion-like particles (ALPs), where their couplings to the Standard Model are of order $(0.01\!-\!1)\,\mbox{TeV}^{-1}$, can be explored by searches for the exotic Higgs decays $h\to Za$ and $h\to aa$ in Run-2 of the LHC. Almost the complete region in which ALPs can explain the anomalous magnetic moment of the muon can be probed by searches for these decays with subsequent decay $a\to\gamma\gamma$, even if the relevant couplings are loop suppressed and the $a\to\gamma\gamma$ branching ratio is less than~1.


title: Existence and uniqueness of global classical solutions to a two species cancer invasion haptotaxis model

abstract: We consider a haptotaxis cancer invasion model that includes two families of cancer cells. Both families, migrate on the extracellular matrix and proliferate. Moreover the model describes an epithelial-to-mesenchymal-like transition between the two families, as well as a degradation and a self-reconstruction process of the extracellular matrix. We prove positivity and conditional global existence and uniqueness of the classical solutions of the problem for large initial data.


title: Cooling-Rate Effects in Sodium Silicate Glasses: Bridging the Gap between Molecular Dynamics Simulations and Experiments

abstract: Although molecular dynamics (MD) simulations are commonly used to predict the structure and properties of glasses, they are intrinsically limited to short time scales, necessitating the use of fast cooling rates. It is therefore challenging to compare results from MD simulations to experimental results for glasses cooled on typical laboratory time scales. Based on MD simulations of a sodium silicate glass with varying cooling rate (from 0.01 to 100 K/ps), here we show that thermal history primarily affects the medium-range order structure, while the short-range order is largely unaffected over the range of cooling rates simulated. This results in a decoupling between the enthalpy and volume relaxation functions, where the enthalpy quickly plateaus as the cooling rate decreases, whereas density exhibits a slower relaxation. Finally, we demonstrate that the outcomes of MD simulations can be meaningfully compared to experimental values if properly extrapolated to slower cooling rates.


title: Discontinuous classical ground state magnetic response as an even-odd effect in higher order rotationally invariant exchange interactions

abstract: The classical ground state magnetic response of the Heisenberg model when rotationally invariant exchange interactions of integer order q>1 are added is found to be discontinuous, even though the interactions lack magnetic anisotropy. This holds even in the case of bipartite lattices which are not frustrated, as well as for the frustrated triangular lattice. The total number of discontinuities is associated with even-odd effects as it depends on the parity of q via the relative strength of the bilinear and higher order exchange interactions, and increases with q. These results demonstrate that the precise form of the microscopic interactions is important for the ground state magnetization response.


title: Vibrational-ground-state zero-width resonances for laser filtration: An extended semiclassical analysis

abstract: A semiclassical model supporting the destructive interference interpretation of zero-width resonances (ZWR) is extended to wavelengths inducing c_minus-type curve crossing situations in Na2 strong field dissociation. This opens the possibility to get critical couples of wavelengths lambda and field intensities I to reach ZWRs associated with the ground vibrationless level v = 0, that, contrary to other vibrational states (v > 0), is not attainable for the commonly referred c+-type crossings. The morphology of such ZWRs in the laser (I; lambda) parameter plane and their usefulness in filtration strategies aiming at molecular cooling down to the ground v = 0 state are examined within the frame of an adiabatic transport scheme.


title: Localization and dynamics of sulfur-oxidizing microbes in natural sediment

abstract: Organic material in anoxic sediment represents a globally significant carbon reservoir that acts to stabilize Earth's atmospheric composition. The dynamics by which microbes organize to consume this material remain poorly understood. Here we observe the collective dynamics of a microbial community, collected from a salt marsh, as it comes to steady state in a two-dimensional ecosystem, covered by flowing water and under constant illumination. Microbes form a very thin front at the oxic-anoxic interface that moves towards the surface with constant velocity and comes to rest at a fixed depth. Fronts are stable to all perturbations while in the sediment, but develop bioconvective plumes in water. We observe the transient formation of parallel fronts. We model these dynamics to understand how they arise from the coupling between metabolism, aerotaxis, and diffusion. These results identify the typical timescale for the oxygen flux and penetration depth to reach steady state.


title: High-Dimensional Function Approximation: Breaking the Curse with Monte Carlo Methods

abstract: In this dissertation we study the tractability of the information-based complexity $n(\varepsilon,d)$ for $d$-variate function approximation problems. In the deterministic setting for many unweighted problems the curse of dimensionality holds, that means, for some fixed error tolerance $\varepsilon>0$ the complexity $n(\varepsilon,d)$ grows exponentially in $d$. For integration problems one can usually break the curse with the standard Monte Carlo method. For function approximation problems, however, similar effects of randomization have been unknown so far. The thesis contains results on three more or less stand-alone topics. For an extended five page abstract, see the section "Introduction and Results". Chapter 2 is concerned with lower bounds for the Monte Carlo error for general linear problems via Bernstein numbers. This technique is applied to the $L_{\infty}$-approximation of certain classes of $C^{\infty}$-functions, where it turns out that randomization does not affect the tractability classification of the problem. Chapter 3 studies the $L_{\infty}$-approximation of functions from Hilbert spaces with methods that may use arbitrary linear functionals as information. For certain classes of periodic functions from unweighted periodic tensor product spaces, in particular Korobov spaces, we observe the curse of dimensionality in the deterministic setting, while with randomized methods we achieve polynomial tractability. Chapter 4 deals with the $L_1$-approximation of monotone functions via function values. It is known that this problem suffers from the curse in the deterministic setting. An improved lower bound shows that the problem is still intractable in the randomized setting. However, Monte Carlo breaks the curse, in detail, for any fixed error tolerance $\varepsilon>0$ the complexity $n(\varepsilon,d)$ grows exponentially in $\sqrt{d}$ only.


title: Nilpotency and the number of word maps of a finite group

abstract: For a finite group $G$ and a non-negative integer $d$, denote by $\Omega_d(G)$ the number of functions $G^d\rightarrow G$ that are induced by substitution into a word with variables among $X_1,\ldots,X_d$. In this note, we show that nilpotency of $G$ can be characterized through the asymptotic growth rate of $\Omega_d(G)$ as $d\to\infty$.


title: A Rate Adaptation Algorithm for Tile-based 360-degree Video Streaming

abstract: In the 360-degree immersive video, a user only views a part of the entire raw video frame based on her viewing direction. However, today's 360-degree video players always fetch the entire panoramic view regardless of users' head movement, leading to significant bandwidth waste that can be potentially avoided. In this paper, we propose a novel adaptive streaming scheme for 360-degree videos. The basic idea is to fetch the invisible portion of a video at the lowest quality based on users' head movement prediction and to adaptively decide the video playback quality for the visible portion based on bandwidth prediction. Doing both in a robust manner requires overcome a series of challenges, such as jointly considering the spatial and temporal domains, tolerating prediction errors, and achieving low complexity. To overcome these challenges, we first define quality of experience (QoE) metrics for adaptive 360-degree video streaming. We then formulate an optimization problem and solve it at a low complexity. The algorithm strategically leverages both future bandwidth and the distribution of users' head positions to determine the quality level of each tile (i.e., a sub-area of a raw frame). We further provide theoretical proof showing that our algorithm achieves optimality under practical assumptions. Numerical results show that our proposed algorithms significantly boost the user QoE by at least 20\% compared to baseline algorithms.


title: Massive pre-main sequence stars in M17

abstract: The formation process of massive stars is still poorly understood. Massive young stellar objects (mYSOs) are deeply embedded in their parental clouds, they are rare and thus typically distant, and their reddened spectra usually preclude the determination of their photospheric parameters. M17 is one of the best studied HII regions in the sky, is relatively nearby, and hosts a young stellar population. With X-shooter on the ESO Very Large Telescope we have obtained optical to near-infrared spectra of candidate mYSOs, identified by Hanson et al. (1997), and a few OB stars in this region. The large wavelength coverage enables a detailed spectroscopic analysis of their photospheres and circumstellar disks. We confirm the pre-main sequence (PMS) nature of six of the stars and characterise the O stars. The PMS stars have radii consistent with being contracting towards the main sequence and are surrounded by a remnant accretion disk. The observed infrared excess and the (double-peaked) emission lines provide the opportunity to measure structured velocity profiles in the disks. We compare the observed properties of this unique sample of young massive stars with evolutionary tracks of massive protostars by Hosokawa & Omukai (2009), and propose that these mYSOs near the western edge of the HII region are on their way to become main-sequence stars ($\sim 6 - 20$ $M_{\odot}$) after having undergone high mass-accretion rates (${\dot{M}_{\rm acc}} \sim 10^{-4} - 10^{-3}$ $M_{\odot}$ $\rm yr^{-1}$). Their spin distribution upon arrival at the zero age main sequence (ZAMS) is consistent with that observed for young B stars, assuming conservation of angular momentum and homologous contraction.


title: Resonances and $CP$ violation in $B_s^0$ and $\overline{B}_s^0 \to J/\psi K^+K^-$ decays in the mass region above the $\phi(1020)$

abstract: The decays of $B_s^0$ and $\overline{B}_s^0$ mesons into the $J/\psi K^+K^-$ final state are studied in the $K^+K^-$ mass region above the $\phi(1020)$ meson in order to determine the resonant substructure and measure the $CP$-violating phase, $\phi_s$, the decay width, $\Gamma_s$, and the width difference between light and heavy mass eigenstates, $\Delta\Gamma_s$. A decay-time dependent amplitude analysis is employed. The data sample corresponds to an integrated luminosity of $3{\rm fb}^{-1}$ produced in 7 and 8 Tev $pp$ collisions at the LHC, collected by the LHCb experiment. The measurement determines $\phi_s = 119\pm107\pm34 {\rm \, mrad}$. A combination with previous LHCb measurements using similar decays into the $J/\psi \pi^+\pi^-$ and $J/\psi\phi(1020)$ final states gives $\phi_s=1\pm37 {\rm \,mrad}$, consistent with the Standard Model prediction.


title: New region force for variational models in image segmentation and high dimensional data clustering

abstract: We propose an effective framework for multi-phase image segmentation and semi-supervised data clustering by introducing a novel region force term into the Potts model. Assume the probability that a pixel or a data point belongs to each class is known a priori. We show that the corresponding indicator function obeys the Bernoulli distribution and the new region force function can be computed as the negative log-likelihood function under the Bernoulli distribution. We solve the Potts model by the primal-dual hybrid gradient method and the augmented Lagrangian method, which are based on two different dual problems of the same primal problem. Empirical evaluations of the Potts model with the new region force function on benchmark problems show that it is competitive with existing variational methods in both image segmentation and semi-supervised data clustering.


title: Dynamic algorithms for visibility polygons in simple polygons

abstract: We devise the following dynamic algorithms for both maintaining as well as querying for the visibility and weak visibility polygons amid vertex insertions and/or deletions to the simple polygon. * A fully-dynamic algorithm for maintaining the visibility polygon of a fixed point located interior to the simple polygon amid vertex insertions and deletions to the simple polygon. The time complexity to update the visibility polygon of a point $q$ due to the insertion (resp. deletion) of vertex $v$ to (resp. from) the current simple polygon is expressed in terms of the number of combinatorial changes needed to the visibility polygon of $q$ due to the insertion (resp. deletion) of $v$. * An output-sensitive query algorithm to answer the visibility polygon query corresponding to any point $p$ in $\mathbb{R}^2$ amid vertex insertions and deletions to the simple polygon. If $p$ is not exterior to the current simple polygon, then the visibility polygon of $p$ is computed. Otherwise, our algorithm outputs the visibility polygon corresponding to the exterior visibility of $p$. * An incremental algorithm to maintain the weak visibility polygon of a fixed-line segment located interior to the simple polygon amid vertex insertions to the simple polygon. The time complexity to update the weak visibility polygon of a line segment $pq$ due to the insertion of vertex $v$ to the current simple polygon is expressed in terms of the sum of the number of combinatorial updates needed to the geodesic shortest path trees rooted at $p$ and $q$ due to the insertion of $v$. * An output-sensitive algorithm to compute the weak visibility polygon corresponding to any query line segment located interior to the simple polygon amid both the vertex insertions and deletions to the simple polygon. Each of these algorithms requires preprocessing the initial simple polygon.


title: Higgs-differential cross section at NNLO in dimensional regularisation

abstract: We present an analytic computation of the Higgs production cross section in the gluon fusion channel, which is differential in the components of the Higgs momentum and inclusive in the associated partonic radiation through NNLO in perturbative QCD. Our computation includes the necessary higher order terms in the dimensional regulator beyond the finite part that are required for renormalisation and collinear factorisation at N$^3$LO. We outline in detail the computational methods which we employ. We present numerical predictions for realistic final state observables, specifically distributions for the decay products of the Higgs boson in the $\gamma\gamma$ decay channel.


title: Numerical studies on the link between radioisotopic signatures on Earth and the formation of the Local Bubble. I. 60Fe transport to the solar system by turbulent mixing of ejecta from nearby supernovae into a locally homogeneous ISM

abstract: The discovery of radionuclides like 60Fe with half-lives of million years in deep-sea crusts and sediments offers the unique possibility to date and locate nearby supernovae. We want to quantitatively establish that the 60Fe enhancement is the result of several supernovae which are also responsible for the formation of the Local Bubble, our Galactic habitat. We performed three-dimensional hydrodynamic adaptive mesh refinement simulations (with resolutions down to subparsec scale) of the Local Bubble and the neighbouring Loop I superbubble in different homogeneous, self-gravitating environments. For setting up the Local and Loop I superbubble, we took into account the time sequence and locations of the generating core-collapse supernova explosions, which were derived from the mass spectrum of the perished members of certain stellar moving groups. The release of 60Fe and its subsequent turbulent mixing process inside the superbubble cavities was followed via passive scalars, where the yields of the decaying radioisotope were adjusted according to recent stellar evolution calculations. The models are able to reproduce both the timing and the intensity of the 60Fe excess observed with rather high precision, provided that the external density does not exceed 0.3 cm-3 on average. Thus the two best-fit models presented here were obtained with background media mimicking the classical warm ionised and warm neutral medium. We also found that 60Fe (which is condensed onto dust grains) can be delivered to Earth via two physical mechanisms: either through individual fast-paced supernova blast waves, which cross the Earth's orbit sometimes even twice as a result of reflection from the Local Bubble's outer shell, or, alternatively, through the supershell of the Local Bubble itself, injecting the 60Fe content of all previous supernovae at once, but over a longer time range.


title: Influence of fracture criteria on dynamic fracture propagation in a discrete chain

abstract: The extent to which time-dependent fracture criteria affect the dynamic behavior of fracture in a discrete structure is discussed in this work. The simplest case of a semi-infinite isotropic chain of oscillators has been studied. Two history-dependent criteria are compared to the classical one of threshold elongation for linear bonds. The results show that steady-state regimes can be reached in the low subsonic crack speed range where it is impossible according to the classical criterion. Repercussions in terms of load and crack opening versus velocity are explained in detail. A strong qualitative influence of history-dependent criteria is observed at low subsonic crack velocities, especially in relation to achievable steady-state propagation regimes.


title: Communication games reveal preparation contextuality

abstract: A communication game consists of distributed parties attempting to jointly complete a task with restricted communication. Such games are useful tools for studying limitations of physical theories. A theory exhibits preparation contextuality whenever its predictions cannot be explained by a preparation noncontextual model. Here, we show that communication games performed in operational theories reveal the preparation contextuality of that theory. For statistics obtained in a particular family of communication games, we show a direct correspondance with correlations in space-like separated events obeying the no-signaling principle. Using this, we prove that all mixed quantum states of any finite dimension are preparation contextual. We report on an experimental realization of a communication game involving three-level quantum systems from which we observe a strong violation of the constraints of preparation noncontextuality.


title: Punny Captions: Witty Wordplay in Image Descriptions

abstract: Wit is a form of rich interaction that is often grounded in a specific situation (e.g., a comment in response to an event). In this work, we attempt to build computational models that can produce witty descriptions for a given image. Inspired by a cognitive account of humor appreciation, we employ linguistic wordplay, specifically puns, in image descriptions. We develop two approaches which involve retrieving witty descriptions for a given image from a large corpus of sentences, or generating them via an encoder-decoder neural network architecture. We compare our approach against meaningful baseline approaches via human studies and show substantial improvements. We find that when a human is subject to similar constraints as the model regarding word usage and style, people vote the image descriptions generated by our model to be slightly wittier than human-written witty descriptions. Unsurprisingly, humans are almost always wittier than the model when they are free to choose the vocabulary, style, etc.


title: Crystal and Magnetic Structures in Layered, Transition Metal Dihalides and Trihalides

abstract: Materials composed of two dimensional layers bonded to one another through weak van der Waals interactions often exhibit strongly anisotropic behaviors and can be cleaved into very thin specimens and sometimes into monolayer crystals. Interest in such materials is driven by the study of low dimensional physics and the design of functional heterostructures. Binary compounds with the compositions MX2 and MX3 where M is a metal cation and X is a halogen anion often form such structures. Magnetism can be incorporated by choosing a transition metal with a partially filled d-shell for M, enabling ferroic responses for enhanced functionality. Here a brief overview of binary transition metal dihalides and trihalides is given, summarizing their crystallographic properties and long-range-ordered magnetic structures, focusing on those materials with layered crystal structures and partially filled d-shells required for combining low dimensionality and cleavability with magnetism.


title: From minimal Lagrangian to J-minimal submanifolds: persistence and uniqueness

abstract: Given a minimal Lagrangian submanifold L in a negative Kaehler--Einstein manifold M, we show that any small Kaehler--Einstein perturbation of M induces a deformation of L which is minimal Lagrangian with respect to the new structure. This provides a new source of examples of minimal Lagrangians. More generally, the same is true for the larger class of totally real J-minimal submanifolds in Kaehler manifolds with negative definite Ricci curvature.


title: Accelerating Stochastic Gradient Descent For Least Squares Regression

abstract: There is widespread sentiment that it is not possible to effectively utilize fast gradient methods (e.g. Nesterov's acceleration, conjugate gradient, heavy ball) for the purposes of stochastic optimization due to their instability and error accumulation, a notion made precise in d'Aspremont 2008 and Devolder, Glineur, and Nesterov 2014. This work considers these issues for the special case of stochastic approximation for the least squares regression problem, and our main result refutes the conventional wisdom by showing that acceleration can be made robust to statistical errors. In particular, this work introduces an accelerated stochastic gradient method that provably achieves the minimax optimal statistical risk faster than stochastic gradient descent. Critical to the analysis is a sharp characterization of accelerated stochastic gradient descent as a stochastic process. We hope this characterization gives insights towards the broader question of designing simple and effective accelerated stochastic methods for more general convex and non-convex optimization problems.


title: Density solutions to a class of integro-differential equations

abstract: We consider the integro-differential equation ${\rm I}^{\alpha}_{0+}f= x^m f$ on the half-line. We show that there exists a density solution, which is then unique and can be expressed in terms of the Beta distribution, if and only if $m> \alpha.$ These density solutions extend the class of generalized one-sided stable distributions introduced in Schneider (1987) and more recently investigated in Pakes (2014). We study various analytical aspects of these densities, and we solve the open problems about infinite divisibility formulated in Pakes (2014).


title: Generalized G-estimation and Model Selection

abstract: Dynamic treatment regimes (DTRs) aim to formalize personalized medicine by tailoring treatment decisions to individual patient characteristics. G-estimation for DTR identification targets the parameters of a structural nested mean model known as the blip function from which the optimal DTR is derived. Despite considerable work deriving such estimation methods, there has been little focus on extending G-estimation to the case of non-additive effects, non-continuous outcomes or on model selection. We demonstrate how G-estimation can be more widely applied through the use of iteratively-reweighted least squares procedures, and illustrate this for log-linear models. We then derive a quasi-likelihood function for G-estimation within the DTR framework, and show how it can be used to form an information criterion for blip model selection. These developments are demonstrated through application to a variety of simulation studies as well as data from the Sequenced Treatment Alternatives to Relieve Depression study.


title: Maximum a posteriori CMB lensing reconstruction

abstract: Gravitational lensing of the CMB is a valuable cosmological signal that correlates to tracers of large-scale structure and acts as a important source of confusion for primordial $B$-mode polarization. State-of-the-art lensing reconstruction analyses use quadratic estimators, which are easily applicable to data. However, these estimators are known to be suboptimal, in particular for polarization, and large improvements are expected to be possible for high signal-to-noise polarization experiments. We develop a method and numerical code, $\rm{LensIt}$, that is able to find efficiently the most probable lensing map, introducing no significant approximations to the lensed CMB likelihood, and applicable to beamed and masked data with inhomogeneous noise. It works by iteratively reconstructing the primordial unlensed CMB using a deflection estimate and its inverse, and removing residual lensing from these maps with quadratic estimator techniques. Roughly linear computational cost is maintained due to fast convergence of iterative searches, combined with the local nature of lensing. The method achieves the maximal improvement in signal to noise expected from analytical considerations on the unmasked parts of the sky. Delensing with this optimal map leads to forecast tensor-to-scalar ratio parameter errors improved by a factor $\simeq 2 $ compared to the quadratic estimator in a CMB stage IV configuration.


title: Estimating the Coefficients of a Mixture of Two Linear Regressions by Expectation Maximization

abstract: We give convergence guarantees for estimating the coefficients of a symmetric mixture of two linear regressions by expectation maximization (EM). In particular, we show that the empirical EM iterates converge to the target parameter vector at the parametric rate, provided the algorithm is initialized in an unbounded cone. In particular, if the initial guess has a sufficiently large cosine angle with the target parameter vector, a sample-splitting version of the EM algorithm converges to the true coefficient vector with high probability. Interestingly, our analysis borrows from tools used in the problem of estimating the centers of a symmetric mixture of two Gaussians by EM. We also show that the population EM operator for mixtures of two regressions is anti-contractive from the target parameter vector if the cosine angle between the input vector and the target parameter vector is too small, thereby establishing the necessity of our conic condition. Finally, we give empirical evidence supporting this theoretical observation, which suggests that the sample based EM algorithm performs poorly when initial guesses are drawn accordingly. Our simulation study also suggests that the EM algorithm performs well even under model misspecification (i.e., when the covariate and error distributions violate the model assumptions).


title: High-precision stellar limb-darkening in exoplanetary transits

abstract: Characterization of the atmospheres of transiting exoplanets relies on accurate measurements of the extent of the optically thick area of the planet at multiple wavelengths with a precision $\lesssim$100 parts per million (ppm). Next-generation instruments onboard the James Webb Space Telescope (JWST) are expected to achieve $\sim$10 ppm precision for several tens of targets. A similar precision can be obtained in modelling only if other astrophysical effects, including the stellar limb-darkening, are accounted for properly. In this paper, we explore the limits on precision due to the mathematical formulas currently adopted to approximate the stellar limb-darkening, and to the use of limb-darkening coefficients obtained either from stellar-atmosphere models or empirically. We propose a new limb-darkening law with two coefficients, `power-2', which outperforms other two-coefficient laws adopted in the literature in most cases, and particularly for cool stars. Empirical limb-darkening based on two-coefficient formulas can be significantly biased, even if the light-curve residuals are nearly photon-noise limited. We demonstrate an optimal strategy to fitting for the four-coefficients limb-darkening in the visible, using prior information on the exoplanet orbital parameters to break some of the degeneracies that otherwise would prevent the convergence of the fit. Infrared observations taken with the James Webb Space Telescope (JWST) will provide accurate measurements of the exoplanet orbital parameters with unprecedented precision, which can be used as priors to improve the stellar limb-darkening characterization, and therefore the inferred exoplanet parameters, from observations in the visible, such as those taken with Kepler/K2, JWST, other past and future instruments.


title: Preimage problems for deterministic finite automata

abstract: Given a subset of states $S$ of a deterministic finite automaton and a word $w$, the preimage is the subset of all states mapped to a state in $S$ by the action of $w$. We study three natural problems concerning words giving certain preimages. The first problem is whether, for a given subset, there exists a word \emph{extending} the subset (giving a larger preimage). The second problem is whether there exists a \emph{totally extending} word (giving the whole set of states as a preimage)---equivalently, whether there exists an \emph{avoiding} word for the complementary subset. The third problem is whether there exists a \emph{resizing} word. We also consider variants where the length of the word is upper bounded, where the size of the given subset is restricted, and where the automaton is strongly connected, synchronizing, or binary. We conclude with a summary of the complexities in all combinations of the cases.


title: Optimal excess-of-loss reinsurance and investment problem for an insurer with default risk under a stochastic volatility model

abstract: In this paper, we study an optimal excess-of-loss reinsurance and investment problem for an insurer in defaultable market. The insurer can buy reinsurance and invest in the following securities: a bank account, a risky asset with stochastic volatility and a defaultable corporate bond. We discuss the optimal investment strategy into two subproblems: a pre-default case and a post-default case. We show the existence of a classical solution to a pre-default case via super-sub solution techniques and give an explicit characterization of the optimal reinsurance and investment policies that maximize the expected CARA utility of the terminal wealth. We prove a verification theorem establishing the uniqueness of the solution. Numerical results are presented in the case of the Scott model and we discuss economic insights obtained from these results.


title: Decremental Data Structures for Connectivity and Dominators in Directed Graphs

abstract: We introduce a new dynamic data structure for maintaining the strongly connected components (SCCs) of a directed graph (digraph) under edge deletions, so as to answer a rich repertoire of connectivity queries. Our main technical contribution is a decremental data structure that supports sensitivity queries of the form "are $ u $ and $ v $ strongly connected in the graph $ G \setminus w $?", for any triple of vertices $ u, v, w $, while $ G $ undergoes deletions of edges. Our data structure processes a sequence of edge deletions in a digraph with $n$ vertices in $O(m n \log{n})$ total time and $O(n^2 \log{n})$ space, where $m$ is the number of edges before any deletion, and answers the above queries in constant time. We can leverage our data structure to obtain decremental data structures for many more types of queries within the same time and space complexity. For instance for edge-related queries, such as testing whether two query vertices $u$ and $v$ are strongly connected in $G \setminus e$, for some query edge $e$. As another important application of our decremental data structure, we provide the first nontrivial algorithm for maintaining the dominator tree of a flow graph under edge deletions. We present an algorithm that processes a sequence of edge deletions in a flow graph in $O(m n \log{n})$ total time and $O(n^2 \log{n})$ space. For reducible flow graphs we provide an $O(mn)$-time and $O(m + n)$-space algorithm. We give a conditional lower bound that provides evidence that these running times may be tight up to subpolynomial factors.


title: Dirac $\delta$-function potential in quasiposition representation of a minimal-length scenario

abstract: A minimal-length scenario can be considered as an effective description of quantum gravity effects. In quantum mechanics the introduction of a minimal length can be accomplished through a generalization of Heisenberg's uncertainty principle. In this scenario, state eigenvectors of the position operator are no longer physical states and the representation in momentum space or a representation in a quasiposition space must be used. In this work, we solve the Schroedinger equation with Dirac $\delta$-function potential in quasiposition space. We calculate the bound state energy and the coefficients of reflection and transmission for scattering states. We show that leading corrections are of order of the minimal length $({\sl O}(\sqrt{\beta}))$ and the coefficients of reflection and transmission are no longer the same for the Dirac delta well and barrier as in ordinary quantum mechanics. Furthermore, assuming that the equivalence of the 1s state energy of the hydrogen atom and the bound state energy of the Dirac $\delta$-function potential in 1-dim is kept in a minimal-length scenario, we also find that the leading correction term for the ground state energy of the hydrogen atom is of order of the minimal length and $\Delta x_{min} \le 10^{-25}$ m.


title: Multi-spatial-mode effects in squeezed-light-enhanced interferometric gravitational wave detectors

abstract: Proposed near-future upgrades of the current advanced interferometric gravitational wave detectors include the usage of frequency dependent squeezed light to reduce the current sensitivity-limiting quantum noise. We quantify and describe the degradation effects that spatial mode-mismatches between optical resonators have on the squeezed field. These mode-mismatches can to first order be described by scattering of light into second-order Gaussian modes. As a demonstration of principle, we also show that squeezing the second-order Hermite-Gaussian modes $\mathrm{HG}_{02}$ and $\mathrm{HG}_{20}$, in addition to the fundamental mode, has the potential to increase the robustness to spatial mode-mismatches. This scheme, however, requires independently optimized squeeze angles for each squeezed spatial mode, which would be challenging to realise in practise.


title: Gravitational allocation for uniform points on the sphere

abstract: Given a collection $\mathcal L$ of $n$ points on a sphere $\mathbf{S}^2_n$ of surface area $n$, a fair allocation is a partition of the sphere into $n$ parts each of area $1$, and each associated with a distinct point of $\mathcal L$. We show that if the $n$ points are chosen uniformly at random and the partition is defined by considering the gravitational field defined by the $n$ points, then the expected distance between a point on the sphere and the associated point of $\mathcal L$ is $O(\sqrt{\log n})$. We use our result to define a matching between two collections of $n$ independent and uniform points on the sphere, and prove that the expected distance between a pair of matched points is $O(\sqrt{\log n})$, which is optimal by a result of Ajtai, Koml\'os, and Tusn\'ady.


title: Exploring Application Performance on Emerging Hybrid-Memory Supercomputers

abstract: Next-generation supercomputers will feature more hierarchical and heterogeneous memory systems with different memory technologies working side-by-side. A critical question is whether at large scale existing HPC applications and emerging data-analytics workloads will have performance improvement or degradation on these systems. We propose a systematic and fair methodology to identify the trend of application performance on emerging hybrid-memory systems. We model the memory system of next-generation supercomputers as a combination of "fast" and "slow" memories. We then analyze performance and dynamic execution characteristics of a variety of workloads, from traditional scientific applications to emerging data analytics to compare traditional and hybrid-memory systems. Our results show that data analytics applications can clearly benefit from the new system design, especially at large scale. Moreover, hybrid-memory systems do not penalize traditional scientific applications, which may also show performance improvement.


title: Towards Printable Natural Dielectric Cloaks via Inverse Scattering Techniques

abstract: The synthesis of non-magnetic 2D dielectric cloaks as proper solutions of an inverse scattering problem is addressed in this paper. Adopting the relevant integral formulation governing the scattering phenomena, analytic and numerical approaches are exploited to provide new insights on how frequency and direction of arrival of the incoming wave may influence the cloaking mechanism in terms of permittivity distribution within the cover region. In quasi-static (subwavelength) regime a solution is analytically derived in terms of homogeneous artificial dielectric cover with $\varepsilon<\varepsilon_0$ which is found to be a necessary and sufficient condition for achieving omnidirectional cloaking. On the other hand, beyond quasi-static regime, the cloaking problem is addressed as an optimization task looking for only natural dielectric coatings with $\varepsilon>\varepsilon_0$ able to hide the object for a given number of directions of the incident field. Simulated results confirm the validity of both analytic and numerical methodologies and allow to estimate effective bandwidths both in terms of frequency range and direction of arrival of the impinging field.


title: The Complexity of Computing a Robust Flow

abstract: Robust network flows are a concept for dealing with uncertainty and unforeseen failures in the network infrastructure. They and their dual counterpart, network flow interdiction, have received steady attention within the operations research community over the past years. One of the most basic models is the Maximum Robust Flow Problem: Given a network and an integer k, the task is to find a path flow of maximum robust value, i.e., the guaranteed value of surviving flow after removal of any k arcs in the network. The complexity of this problem appeared to have been settled almost a decade ago: Aneja et al. (2001) showed that the problem can be solved efficiently when k = 1, while an article by Du and Chandrasekaran (2007) established that the problem is NP-hard for any constant value of k larger than 1. We point to a flaw in the proof of the latter result, leaving the complexity for constant k open once again. For the case that k is not bounded by a constant, we present a new hardness proof, establishing that the problem is strongly NP-hard, even when only two different capacity values occur and the number of paths is polynomial in the size of the input. We further show that computing optimal integral solutions is already NP-hard for k = 2 (whereas for k=1, an efficient algorithm is known) and give a positive result for the case that capacities are in {1, 2}.


title: Experimental Two-dimensional Quantum Walk on a Photonic Chip

abstract: Quantum walks, in virtue of the coherent superposition and quantum interference, possess exponential superiority over its classical counterpart in applications of quantum searching and quantum simulation. The quantum enhanced power is highly related to the state space of quantum walks, which can be expanded by enlarging the photon number and/or the dimensions of the evolution network, but the former is considerably challenging due to probabilistic generation of single photons and multiplicative loss. Here we demonstrate a two-dimensional continuous-time quantum walk by using the external geometry of photonic waveguide arrays, rather than the inner degree of freedoms of photons. Using femtosecond laser direct writing, we construct a large-scale three-dimensional structure which forms a two-dimensional lattice with up to 49X49 nodes on a photonic chip. We demonstrate spatial two-dimensional quantum walks using heralded single photons and single-photon-level imaging. We analyze the quantum transport properties via observing the ballistic evolution pattern and the variance profile, which agree well with simulation results. We further reveal the transient nature that is the unique feature for quantum walks of beyond one dimension. An architecture that allows a walk to freely evolve in all directions and a large scale, combining with defect and disorder control, may bring up powerful and versatile quantum walk machines for classically intractable problems.