title: Exclusive Backward-Angle Omega Meson Electroproduction

abstract: This work is a pioneering study of backward-angle omega cross sections through the exclusive 1H(e,e'p)omega reaction using the missing mass reconstruction technique. The extracted cross sections are separated into the transverse (T), longitudinal (L), and LT, TT interference terms. The analyzed data were part of experiment E01-004 (Fpi-2), which used 2.6-5.2GeV electron beams and HMS+SOS spectrometers in Jefferson Lab Hall C. The primary objective was to detect coincidence pion in the forward-angle, where the backward-angle omega events were fortuitously detected. The experiment has central Q2 values of 1.60 and 2.45GeV2, at W=2.21GeV. There was significant coverage in phi and epsilon, which allowed separation of sigma_T,L,LT,TT. The data set has a unique u coverage of -u~0, which corresponds to -t>4GeV2. The separated sigma_T result suggest a flat ~1/Q^(1.33+/-1.21) dependence, whereas sigma_L seems to hold a stronger 1/Q^(9.43+/-6.28) dependence. The sigma_L/sigma_T ratio indicate sigma_T dominance at Q2=2.45 GeV2 at the ~90% confidence level. After translating the results into the -t space of the published CLAS data, our data show evidence of a backward-angle omega electroproduction peak at both Q2 settings. Previously, this phenomenon showing both forward and backward-angle peaks was only observed in the meson photoproduction data. Through comparison of our sigma_T data with the prediction of the Transition Distribution Amplitude (TDA) model, and signs of sigma_T dominance, promising indications of the applicability of the TDA factorization are demonstrated at a much lower Q2 value than its preferred range of Q2>10GeV2. These studies have opened a new means to study the transition of the nucleon wavefunction through backward-angle experimental observables.


title: Superconducting and normal-state anisotropy of the doped topological insulator Sr$_{0.1}$Bi$_2$Se$_3$

abstract: Sr$_x$Bi$_2$Se$_3$ and the related compounds Cu$_x$Bi$_2$Se$_3$ and Nb$_x$Bi$_2$Se$_3$ have attracted considerable interest, as these materials may be realizations of unconventional topological superconductors. Superconductivity with T$_c$ ~ 3 K in Sr$_x$Bi$_2$Se$_3$ arises upon intercalation of Sr into the layered topological insulator Bi$_2$Se$_3$. Here we elucidate the anisotropy of the normal and superconducting state of Sr$_{0.1}$Bi$_2$Se$_3$ with angular dependent magnetotransport and thermodynamic measurements. High resolution x-ray diffraction studies underline the high crystalline quality of the samples. We demonstrate that the normal state electronic and magnetic properties of Sr$_{0.1}$Bi$_2$Se$_3$ are isotropic in the basal plane while we observe a large two-fold in-plane anisotropy of the upper critical field in the superconducting state. Our results support the recently proposed odd-parity nematic state characterized by a nodal gap of $E_u$ symmetry in Sr$_x$Bi$_2$Se$_3$.


title: Qatar Exoplanet Survey: Qatar-6b -- a grazing transiting hot Jupiter

abstract: We report the discovery of Qatar-6b, a new transiting planet identified by the Qatar Exoplanet Survey (QES). The planet orbits a relatively bright (V=11.44), early-K main-sequence star at an orbital period of P~3.506 days. An SED fit to available multi-band photometry, ranging from the near-UV to the mid-IR, yields a distance of d = 101 +/- 6 pc to the system. From a global fit to follow-up photometric and spectroscopic observations, we calculate the mass and radius of the planet to be Mp = 0.67 +/- 0.07 Mjup and Rp = 1.06 +/- 0.07 Rjup, respectively. We use multi-color photometric light curves to show that the transit is grazing, making Qatar-6b one of the few exoplanets known in a grazing transit configuration. It adds to the short list of targets that offer the best opportunity to look for additional bodies in the host planetary system through variations in the transit impact factor and duration.


title: Basic Thresholding Classification

abstract: In this thesis, we propose a light-weight sparsity-based algorithm, basic thresholding classifier (BTC), for classification applications (such as face identification, hyper-spectral image classification, etc.) which is capable of identifying test samples extremely rapidly and performing high classification accuracy. Originally BTC is a linear classifier which works based on the assumption that the samples of the classes of a given dataset are linearly separable. However, in practice those samples may not be linearly separable. In this context, we also propose another algorithm namely kernel basic thresholding classifier (KBTC) which is a non-linear kernel version of the BTC algorithm. KBTC can achieve promising results especially when the given samples are linearly non-separable. For both proposals, we introduce sufficient identification conditions (SICs) under which BTC and KBTC can identify any test sample in the range space of a given dictionary. By using SICs, we develop parameter estimation procedures which do not require any cross validation. Both BTC and KBTC algorithms provide efficient classifier fusion schemes in which individual classifier outputs are combined to produce better classification results. For instance, for the application of face identification, this is done by combining the residuals having different random projectors. For spatial applications such as hyper-spectral image classification, the fusion is carried out by incorporating the spatial information, in which the output residual maps are filtered using a smoothing filter. Numerical results on publicly available face and hyper-spectral datasets show that our proposal outperforms well-known support vector machines (SVM)-based techniques, multinomial logistic regression (MLR)-based methods, and sparsity-based approaches like $l_1$-minimization and simultaneous orthogonal matching pursuit (SOMP).


title: Realizing a Circuit Analog of an Optomechanical System with Longitudinally Coupled Superconducting Resonators

abstract: We realize a superconducting circuit analog of the generic cavity-optomechanical Hamiltonian by longitudinally coupling two superconducting resonators, which are an order of magnitude different in frequency. We achieve longitudinal coupling by embedding a superconducting quantum interference device (SQUID) into a high frequency resonator, making its resonance frequency depend on the zero point current fluctuations of a nearby low frequency LC-resonator. By employing sideband drive fields we enhance the intrinsic coupling strength of about 15 kHz up to 280 kHz by controlling the amplitude of the drive field. Our results pave the way towards the exploration of optomechanical effects in a fully superconducting platform and could enable quantum optics experiments with photons in the yet unexplored radio frequency band.


title: Strong convergence of quantum channels: continuity of the Stinespring dilation and discontinuity of the unitary dilation

abstract: We show that a sequence $\{\Phi_n\}$ of quantum channels strongly converges to a quantum channel $\Phi_0$ if and only if there exist a common environment for all the channels and a corresponding sequence $\{V_n\}$ of Stinespring isometries strongly converging to a Stinespring isometry $V_0$ of the channel $\Phi_0$. We also give a quantitative description of the above characterization of the strong convergence in terms of the appropriate metrics on the sets of quantum channels and Stinespring isometries. As a result, the uniform selective continuity of the complementary operation with respect to the strong convergence is established. We show discontinuity of the unitary dilation by constructing a strongly converging sequence of channels which can not be represented as a reduction of a strongly converging sequence of unitary channels. The Stinespring representation of strongly converging sequences of quantum channels allows to prove the lower semicontinuity of the entropic disturbance as a function of a pair (channel, input ensemble). Some corollaries of this property are considered.

categories: quant-ph math-ph math.FA math.MP math.OA


title: Updated determination of $D^0$-$\overline{D}{}^0$ mixing and CP violation parameters with $D^0\to K^+\pi^-$ decays

abstract: We report measurements of charm-mixing parameters based on the decay-time-dependent ratio of $D^0\to K^+\pi^-$ to $D^0\to K^-\pi^+$ rates. The analysis uses a data sample of proton-proton collisions corresponding to an integrated luminosity of $5.0$ fb$^{-1}$ recorded by the LHCb experiment from 2011 through 2016. Assuming charge-parity (CP) symmetry, the mixing parameters are determined to be $x'^2=(3.9 \pm 2.7) \times10^{-5}$, $y'=(5.28 \pm 0.52) \times 10^{-3}$, and $R_D=(3.454 \pm 0.031)\times10^{-3}$. Without this assumption, the measurement is performed separately for $D^0$ and $\overline{D}{}^0$ mesons, yielding a direct CP-violating asymmetry $A_D =(-0.1\pm9.1)\times10^{-3}$, and magnitude of the ratio of mixing parameters $1.00< |q/p| <1.35$ at the $68.3\%$ confidence level. All results include statistical and systematic uncertainties and improve significantly upon previous single-measurement determinations. No evidence for CP violation in charm mixing is observed.


title: Nuclear modification of scalar, axial and tensor charges from lattice QCD

abstract: Complete flavour decompositions of the scalar, axial and tensor charges of the proton, deuteron, diproton and $^3$He at SU(3)-symmetric values of the quark masses corresponding to a pion mass $m_\pi\sim806$ MeV are determined using lattice QCD. At the physical quark masses, the scalar charges constrain mean-field models of nuclei and the low-energy interactions of nuclei with potential dark matter candidates. The axial and tensor charges of nuclei constrain their spin content, integrated transversity and the quark contributions to their electric dipole moments. External fields are used to directly access the quark-line connected matrix elements of quark bilinear operators, and a combination of stochastic estimation techniques is used to determine the disconnected sea-quark contributions. Significant nuclear modifications are found, with particularly large, O(10%), effects in the scalar charges. Typically, these nuclear effects reduce the effective charge of the nucleon (quenching), although in some cases an enhancement is not excluded. Given the size of the nuclear modifications of the scalar charges resolved here, contributions from correlated multi-nucleon effects should be quantified in the analysis of dark matter direct-detection experiments using nuclear targets.


title: Nanophotonic Particle Simulation and Inverse Design Using Artificial Neural Networks

abstract: We propose a method to use artificial neural networks to approximate light scattering by multilayer nanoparticles. We find the network needs to be trained on only a small sampling of the data in order to approximate the simulation to high precision. Once the neural network is trained, it can simulate such optical processes orders of magnitude faster than conventional simulations. Furthermore, the trained neural network can be used solve nanophotonic inverse design problems by using back- propogation - where the gradient is analytical, not numerical.


title: S-Shaped vs. V-Shaped Transfer Functions for Antlion Optimization Algorithm in Feature Selection Problems

abstract: Feature selection is an important preprocessing step for classification problems. It deals with selecting near optimal features in the original dataset. Feature selection is an NP-hard problem, so meta-heuristics can be more efficient than exact methods. In this work, Ant Lion Optimizer (ALO), which is a recent metaheuristic algorithm, is employed as a wrapper feature selection method. Six variants of ALO are proposed where each employ a transfer function to map a continuous search space to a discrete search space. The performance of the proposed approaches is tested on eighteen UCI datasets and compared to a number of existing approaches in the literature: Particle Swarm Optimization, Gravitational Search Algorithm, and two existing ALO-based approaches. Computational experiments show that the proposed approaches efficiently explore the feature space and select the most informative features, which help to improve the classification accuracy.


title: Boltzmann games in heterogeneous consensus dynamics

abstract: We consider a constrained hierarchical opinion dynamics in the case of leaders' competition and with complete information among leaders. Each leaders' group tries to drive the followers' opinion towards a desired state accordingly to a specific strategy. By using the Boltzmann-type control approach we analyze the best-reply strategy for each leaders' population. Derivation of the corresponding Fokker-Planck model permits to investigate the asymptotic behaviour of the solution. Heterogeneous followers populations are then considered where the effect of knowledge impacts the leaders' credibility and modifies the outcome of the leaders' competition.


title: Positivity-preserving, energy stable numerical schemes for the Cahn-Hilliard equation with logarithmic potential

abstract: We present and analyze finite difference numerical schemes for the Allen Cahn/Cahn-Hilliard equation with a logarithmic Flory Huggins energy potential. Both the first order and second order accurate temporal algorithms are considered. In the first order scheme, we treat the nonlinear logarithmic terms and the surface diffusion term implicitly, and update the linear expansive term and the mobility explicitly. We provide a theoretical justification that, this numerical algorithm has a unique solution such that the positivity is always preserved for the logarithmic arguments. In particular, our analysis reveals a subtle fact: the singular nature of the logarithmic term around the values of $-1$ and 1 prevents the numerical solution reaching these singular values, so that the numerical scheme is always well-defined as long as the numerical solution stays similarly bounded at the previous time step. Furthermore, an unconditional energy stability of the numerical scheme is derived, without any restriction for the time step size. The unique solvability and the positivity-preserving property for the second order scheme are proved using similar ideas, in which the singular nature of the logarithmic term plays an essential role. For both the first and second order accurate schemes, we are able to derive an optimal rate convergence analysis, which gives the full order error estimate. The case with a non-constant mobility is analyzed as well. We also describe a practical and efficient multigrid solver for the proposed numerical schemes, and present some numerical results, which demonstrate the robustness of the numerical schemes.


title: Redundant edges in Ramsey graphs

abstract: For graphs $G$, $F$ and $H$, let $G\rightarrow (F,H)$ signify that any edge coloring of $G$ in red and blue contains a red $F$ or a blue $H$. The Ramsey number $R(F,H)=\min\{r|\; K_r\rightarrow (F,H)\}$. In this note, we consider redundant edges in Ramsey graphs, which are associate with critical Ramsey numbers. For an integer $k\ge 1$, let ${\mathbb G}=\{G_k,G_{k+1},\dots \}$ be a class of graphs with $\delta(G_n)\ge 1$. We define the critical Ramsey number $R_{\mathbb G}(F,H)$ with respect to $\mathbb G$ to be $\max\{n|\; K_r\setminus G_n \rightarrow(F,H),\,G_n\in{\mathbb G}\big\}$, where $r=R(F,H)$. We shall determine some $R_{\mathbb G}(F,H)$, where ${\mathbb G}$ consists of stars, matchings and complete graphs, respectively.


title: A Local-Realistic Model of Quantum Mechanics Based on a Discrete Spacetime (Extended version)

abstract: This paper presents a realistic, stochastic, and local model that reproduces nonrelativistic quantum mechanics (QM) results without using its mathematical formulation. The proposed model only uses integer-valued quantities and operations on probabilities, in particular assuming a discrete spacetime under the form of a Euclidean lattice. Individual (spinless) particle trajectories are described as random walks. Transition probabilities are simple functions of a few quantities that are either randomly associated to the particles during their preparation, or stored in the lattice nodes they visit during the walk. QM predictions are retrieved as probability distributions of similarly-prepared ensembles of particles. The scenarios considered to assess the model comprise of free particle, constant external force, harmonic oscillator, particle in a box, the Delta potential, particle on a ring, particle on a sphere and include quantization of energy levels and angular momentum.


title: Music Transcription by Deep Learning with Data and "Artificial Semantic" Augmentation

abstract: In this progress paper the previous results of the single note recognition by deep learning are presented. The several ways for data augmentation and "artificial semantic" augmentation are proposed to enhance efficiency of deep learning approaches for monophonic and polyphonic note recognition by increase of dimensions of training data, their lossless and lossy transformations.


title: A Theoretical Study of Aqueous Humor Secretion Based on a Continuum Model Coupling Electrochemical and Fluid-Dynamical Transmembrane Mechanisms

abstract: Intraocular pressure, resulting from the balance of aqueous humor (AH) production and drainage, is the only approved treatable risk factor in glaucoma. AH production is determined by the concurrent function of ionic pumps and aquaporins in the ciliary processes but their individual contribution is difficult to characterize experimentally. In this work, we propose a novel unified modeling and computational framework for the finite element simulation of the role of the main ionic pumps involved in AH secretion, namely, the sodium potassium pump, the calcium-sodium pump, the anion channel and the hydrogenate-sodium pump. The theoretical model is developed at the cellular scale and is based on the coupling between electrochemical and fluid-dynamical transmembrane mechanisms characterized by a novel description of the electric pressure exerted by the ions on the intrachannel fluid that includes electrochemical and osmotic corrections. Considering a realistic geometry of the ionic pumps, the proposed model is demonstrated to correctly predict their functionality as a function of (1) the permanent electric charge density over the channel pump surface; (2) the osmotic gradient coefficient; (3) the stoichiometric ratio between the ionic pump currents enforced at the inlet and outlet sections of the channel. In particular, theoretical predictions of the transepithelial membrane potential for each simulated pump/channel allow us to perform a first significant model comparison with experimental data for monkeys. This is a significant step for future multidisciplinary studies on the action of molecules on AH production.


title: Ghosts of Jupiter's past: is 2017 UV43 a relative of comet Shoemaker-Levy 9?

abstract: The present-day orbit of 2017 UV43 resembles the average one of the fragments of comet Shoemaker-Levy 9. Here, we present numerical evidence suggesting that 2017 UV43 might be related to Shoemaker-Levy 9. With our exploratory results we intend to encourage a search for precovery images of 2017 UV43 and perhaps even follow-up observations that may help in improving its poorly determined orbit so its true dynamical nature can be understood.


title: What can we learn from the stochastic gravitational wave background produced by oscillons?

abstract: The stochastic gravitational wave (GW) background provides a fascinating window to the physics of the very early universe. Beyond the nearly scale-invariant primordial GW spectrum produced during inflation, a spectrum with a much richer structure is typically generated during the preheating phase after inflation (or after some other phase transition at lower energies). This raises the question of what one can learn from a future observation of the stochastic gravitational wave background spectrum about the underlying physics during preheating. Recently, it has been shown that during preheating non-perturbative quasi-stable objects like oscillons can act as strong sources for GW, leading to characteristic features such as distinct peaks in the spectrum. In this paper, we study the GW production from oscillons using semi-analytical techniques. In particular, we discuss how the GW spectrum is affected by the parameters that characterise a given oscillon system, e.g. by the background cosmology, the asymmetry of the oscillons and the evolution of the number density of the oscillons. We compare our semi-analytic results with numerical lattice simulations for a hilltop inflation model and a KKLT scenario, which differ strongly in some of these characteristics, and find very good agreement.


title: All-optical signatures of Strong-Field QED in the vacuum emission picture

abstract: We study all-optical signatures of the effective nonlinear couplings among electromagnetic fields in the quantum vacuum, using the collision of two focused high-intensity laser pulses as an example. The experimental signatures of quantum vacuum nonlinearities are encoded in signal photons, whose kinematic and polarization properties differ from the photons constituting the macroscopic laser fields. We implement an efficient numerical algorithm allowing for the theoretical investigation of such signatures in realistic field configurations accessible in experiment. This algorithm is based on a vacuum emission scheme and can readily be adapted to the collision of more laser beams or further involved field configurations. We solve the case of two colliding pulses in full 3+1 dimensional spacetime, and identify experimental geometries and parameter regimes with improved signal-to-noise ratios.


title: On the Largest Part Size and Its Multiplicity of a Random Integer Partition

abstract: Let $\lambda$ be a partition of the positive integer $n$ chosen umiformly at random among all such partitions. Let $L_n=L_n(\lambda)$ and $M_n=M_n(\lambda)$ be the largest part size and its multiplicity, respectively. For large $n$, we focus on a comparison between the partition statistics $L_n$ and $L_n M_n$. In terms of convergence in distribution, we show that they behave in the same way. However, it turns out that the expectation of $L_n M_n -L_n$ grows as fast as $\frac{1}{2}\log{n}$ We obtain a precise asymptotic expansion for this expectation and conclude with an open problem arising from this study.


title: Primitive ideal space of Higher-rank graph $C^*$-algebras and decomposability

abstract: In this paper, we describe primitive ideal space of the $C^*$-algebra $C^*(\Lambda)$ associated to any locally convex row-finite $k$-graph $\Lambda$. To do this, we will apply the Farthing's desourcifying method on a recent result of Carlsen, Kang, Shotwell, and Sims. We also characterize certain maximal ideals of $C^*(\Lambda)$. Furthermore, we study the decomposability of $C^*(\Lambda)$. We apply the description of primitive ideals to show that if $I$ is a direct summand of $C^*(\Lambda)$, then it is gauge-invariant and isomorphic to a certain $k$-graph $C^*$-algebra. So, we may characterize decomposable higher-rank $C^*$-algebras by giving necessary and sufficient conditions for the underlying $k$-graphs. Moreover, we determine all such $C^*$-algebras which can be decomposed into a direct sum of finitely many indecomposable $C^*$-algebras.


title: AdS$_7$/CFT$_6$ with orientifolds

abstract: AdS$_7$ solutions of massive type IIA have been classified, and are dual to a large class of six-dimensional $(1,0)$ SCFT's whose tensor branch deformations are described by linear quivers of SU groups. Quivers and AdS vacua depend solely on the group theory data of the NS5-D6-D8 brane configurations engineering the field theories. This has allowed for a direct holographic match of their $a$ conformal anomaly. In this paper we extend the match to cases where O6 and O8-planes are present, thereby introducing SO and USp groups in the quivers. In all of them we show that the $a$ anomaly computed in supergravity agrees with the holographic limit of the exact field theory result, which we extract from the anomaly polynomial. As a byproduct we construct special AdS$_7$ vacua dual to nonperturbative F-theory configurations. Finally, we propose a holographic $a$-theorem for six-dimensional Higgs branch RG flows.


title: Black hole mass formula in the membrane paradigm

abstract: The membrane paradigm approach adopts a timelike surface, stretched out off the null event horizon, to study several important black hole properties. We use this powerful tool to give a direct derivation of the black hole mass formula in the static and stationary cases without and with electric field. Since here the membrane is a self-gravitating material system we go beyond the usual applicability on test particles and test fields of the paradigm.


title: Late time afterglow observations reveal a collimated relativistic jet in the ejecta of the binary neutron star merger GW170817

abstract: The binary neutron star (BNS) merger GW170817 was the first astrophysical source detected in gravitational waves and multi-wavelength electromagnetic radiation. The almost simultaneous observation of a pulse of gamma-rays proved that BNS mergers are associated with at least some short gamma-ray bursts (GRBs). However, the gamma-ray pulse was faint, casting doubts on the association of BNS mergers with the luminous, highly relativistic outflows of canonical short GRBs. Here we show that structured jets with a relativistic, energetic core surrounded by slower and less energetic wings produce afterglow emission that brightens characteristically with time, as recently seen in the afterglow of GW170817. Initially, we only see the relatively slow material moving towards us. As time passes, larger and larger sections of the outflow become visible, increasing the luminosity of the afterglow. The late appearance and increasing brightness of the multi-wavelength afterglow of GW170817 allow us to constrain the geometry of its ejecta and thus reveal the presence of an off-axis jet pointing about 30 degrees away from Earth. Our results confirm a single origin for BNS mergers and short GRBs: GW170817 produced a structured outflow with a highly relativistic core and a canonical short GRB. We did not see the bright burst because it was beamed away from Earth. However, approximately one in 20 mergers detected in gravitational waves will be accompanied by a bright, canonical short GRB.


title: Fibonacci Topological Superconductor

abstract: We introduce a model of interacting Majorana fermions that describes a superconducting phase with a topological order characterized by the Fibonacci topological field theory. Our theory, which is based on a $SO(7)_1/(G_2)_1$ coset factorization, leads to a solvable one dimensional model that is extended to two dimensions using a network construction. In addition to providing a description of the Fibonacci phase without parafermions, our theory predicts a closely related "anti-Fibonacci" phase, whose topological order is characterized by the tricritical Ising model. We show that Majorana fermions can split into a pair of Fibonacci anyons, and propose an interferometer that generalizes the $Z_2$ Majorana interferometer and directly probes the Fibonacci non-Abelian statistics.


title: Electron Heating in Low Mach Number Perpendicular shocks. II. Dependence on the Pre-Shock Conditions

abstract: Recent X-ray observations of merger shocks in galaxy clusters have shown that the post-shock plasma is two-temperature, with the protons being hotter than the electrons. In this work, the second of a series, we investigate by means of two-dimensional particle-in-cell simulations the efficiency of electron irreversible heating in perpendicular low Mach number shocks. We consider values of plasma beta (ratio of thermal and magnetic pressures) in the range $4\lesssim \beta_{p0}\lesssim 32$ and sonic Mach number (ratio of shock speed to pre-shock sound speed) in the range $2\lesssim M_{s}\lesssim 5$, as appropriate for galaxy cluster shocks. As shown in Paper I, magnetic field amplification - induced by shock compression of the pre-shock field, or by strong proton cyclotron and mirror modes accompanying the relaxation of proton temperature anisotropy - can drive the electron temperature anisotropy beyond the threshold of the electron whistler instability. The growth of whistler waves breaks the electron adiabatic invariance, and allows for efficient entropy production. We find that the post-shock electron temperature $T_{e2}$ exceeds the adiabatic expectation $T_{e2,\rm ad}$ by an amount $(T_{e2}-T_{e2,\rm ad})/T_{e0}\simeq 0.044 \,M_s (M_s-1)$ (here, $T_{e0}$ is the pre-shock temperature), which depends only weakly on the plasma beta, over the range $4\lesssim \beta_{p0}\lesssim 32$ which we have explored, and on the proton-to-electron mass ratio (the coefficient of $\simeq 0.044$ is measured for our fiducial $m_i/m_e=49$, and we estimate that it will decrease to $\simeq 0.03$ for the realistic mass ratio). Our results have important implications for current and future observations of galaxy cluster shocks in the radio band (synchrotron emission and Sunyaev-Zel'dovich effect) and at X-ray frequencies.


title: GW170817 Most Likely Made a Black Hole

abstract: There are two outstanding issues regarding the neutron-star merger event GW170817: the nature of the compact remnant and the interstellar shock. The mass of the remnant of GW170817, $\sim$2.7 $M_\odot$, implies the remnant could be either a massive, rotating, neutron star, or a black hole. We report Chandra Director's Discretionary Time observations made in 2017 December and 2018 January, and we reanalyze earlier observations from 2017 August and 2017 September, in order to address these unresolved issues. We estimate the X-ray flux from a neutron star remnant and compare that to the measured X-ray flux. If we assume that the spin-down luminosity of any putative neutron star is converted to pulsar wind nebula X-ray emission in the 0.5-8 keV band with an efficiency of $10^{-3}$, for a dipole magnetic field with $3 \times 10^{11}$ G < $B$ < $10^{14}$ G, a rising X-ray signal would result and would be brighter than that observed by day 107, we therefore conclude that the remnant of GW170817 is most likely a black hole. Independent of any assumptions of X-ray efficiency, however, if the remnant is a rapidly-rotating, magnetized, neutron star, the total energy in the external shock should rise by a factor $\sim$$10^2$ (to $\sim$$10^{52}$ erg) after a few years, therefore, Chandra observations over the next year or two that do not show substantial brightening will rule out such a remnant. The same observations can distinguish between two different models for the relativistic outflow, either an angular or radially varying structure.


title: KELT-21b: A Hot Jupiter Transiting the Rapidly-Rotating Metal-Poor Late-A Primary of a Likely Hierarchical Triple System

abstract: We present the discovery of KELT-21b, a hot Jupiter transiting the $V=10.5$ A8V star HD 332124. The planet has an orbital period of $P=3.6127647\pm0.0000033$ days and a radius of $1.586_{-0.040}^{+0.039}$ $R_J$. We set an upper limit on the planetary mass of $M_P<3.91$ $M_J$ at $3\sigma$ confidence. We confirmed the planetary nature of the transiting companion using this mass limit and Doppler tomographic observations to verify that the companion transits HD 332124. These data also demonstrate that the planetary orbit is well-aligned with the stellar spin, with a sky-projected spin-orbit misalignment of $\lambda=-5.6_{-1.9}^{+1.7 \circ}$. The star has $T_{\mathrm{eff}}=7598_{-84}^{+81}$ K, $M_*=1.458_{-0.028}^{+0.029}$ $M_{\odot}$, $R_*=1.638\pm0.034$ $R_{\odot}$, and $v\sin I_*=146$ km s$^{-1}$, the highest projected rotation velocity of any star known to host a transiting hot Jupiter. The star also appears to be somewhat metal-poor and $\alpha$-enhanced, with [Fe/H]$=-0.405_{-0.033}^{+0.032}$ and [$\alpha$/Fe]$=0.145 \pm 0.053$; these abundances are unusual, but not extraordinary, for a young star with thin-disk kinematics like KELT-21. High-resolution imaging observations revealed the presence of a pair of stellar companions to KELT-21, located at a separation of 1.2" and with a combined contrast of $\Delta K_S=6.39 \pm 0.06$ with respect to the primary. Although these companions are most likely physically associated with KELT-21, we cannot confirm this with our current data. If associated, the candidate companions KELT-21 B and C would each have masses of $\sim0.12$ $M_{\odot}$, a projected mutual separation of $\sim20$ AU, and a projected separation of $\sim500$ AU from KELT-21. KELT-21b may be one of only a handful of known transiting planets in hierarchical triple stellar systems.


title: Steeper stellar cusps in galactic centers from binary disruption

abstract: The relaxed distribution of stars around a massive black hole is known to follow a cusp profile $\rho(r)\propto r^{-\alpha}$ with characteristic slope $\alpha=7/4$. This follows from energy conservation and a scattering rate as given by two body encounters. However, we show that injection of stars close to the black hole, i.e. a source term in the standard cusp picture, modifies this profile. In the steady-state configuration, the cusp develops a central region with typical slope $\alpha=9/4$ in which stars diffuse outward. Binary disruption by the intense tidal field of the massive black hole is among the phenomena that take place in the Galactic Center. In such disruption, one of the binary members remains bound to the black hole, thus providing a source term of stars close to the black hole. Assuming a binary fraction of $0.1$ and an orbital circularization efficiency of $0.35$, we show that this source is strong enough to modify the cusp profile within $\approx 0.07$ pc in the Galactic Center. If the binary fraction at the influence radius is of order unity and the orbits of all the captured stars are efficiently circularized, the steeper cusp extends almost as far as the radius of influence of the black hole.


title: Non-Perturbative Dynamical Effects in Nearly Scale Invariant Systems: The Action of Breaking Scale Invariance

abstract: In this work we develop a general formalism that categorizes the action of broken scale invariance on the non-equilibrium dynamics of non-relativistic quantum systems. This approach is equally applicable to both strongly and weakly interacting systems. We show that any small deviation from the strongly interacting fixed point, in three spatial dimensions, leads to non-pertubative effects in the long time dynamics, dramatically altering the dynamics observed at the scale invariant fixed point. As a concrete example, we apply this approach to the non-equilibrium dynamics for the interacting two-body problem, and for a non-interacting quantum gas in the presence of an impurity, both in three spatial dimensions. Slightly away from the resonantly-interacting scale invariant fixed point, we show that the dynamics are altered by a non-perturbative log-periodic beat. The presence of the beat depends only on deviating from the resonant fixed point, while the frequency depends on the microscopic parameters of the system.


abstract: We list 4d interacting $\mathcal{N}=2$ SCFT with minimal flavor central charge from the theory space constructed using 6d $(2,0)$ theory. For $ADE$ and $C_N$ flavor groups, our theory saturates the bound found using bootstrap method, but other cases have higher values. We find interesting rank one SCFTs with $B_3, G_2, F_4, C_4\times U(1), C_1\times U(1)$ flavor symmetry. Many physical properties of these theories are also studied.


title: Identical synchronization of nonidentical oscillators: when only birds of different feathers flock together

abstract: An outstanding problem in the study of networks of heterogeneous dynamical units concerns the development of rigorous methods to probe the stability of synchronous states when the differences between the units are not small. Here, we address this problem by presenting a generalization of the master stability formalism that can be applied to heterogeneous oscillators with large mismatches. Our approach is based on the simultaneous block diagonalization of the matrix terms in the variational equation, and it leads to dimension reduction that simplifies the original equation significantly. This new formalism allows the systematic investigation of scenarios in which the oscillators need to be nonidentical in order to reach an identical state, where all oscillators are completely synchronized. In the case of networks of identically coupled oscillators, this corresponds to breaking the symmetry of the system as a means to preserve the symmetry of the dynamical state---a recently discovered effect termed asymmetry-induced synchronization (AISync). Our framework enables us to identify communication delay as a new and potentially common mechanism giving rise to AISync, which we demonstrate using networks of delay-coupled Stuart-Landau oscillators. The results also have potential implications for control, as they reveal oscillator heterogeneity as an attribute that may be manipulated to enhance the stability of synchronous states.


title: Priority-Aware Near-Optimal Scheduling for Heterogeneous Multi-Core Systems with Specialized Accelerators

abstract: To deliver high performance in power limited systems, architects have turned to using heterogeneous systems, either CPU+GPU or mixed CPU-hardware systems. However, in systems with different processor types and task affinities, scheduling tasks becomes more challenging than in homogeneous multi-core systems or systems without task affinities. The problem is even more complex when specialized accelerators and task priorities are included. In this paper, we provide a formal proof for the optimal scheduling policy for heterogeneous systems with arbitrary number of resource types, including specialized accelerators, independent of the task arrival rate, task size distribution, and resource processing order. We transform the optimal scheduling policy to a nonlinear integer optimization problem and propose a fast, near-optimal algorithm. An additional heuristic is proposed for the case of priority-aware scheduling. Our experimental results demonstrate that the proposed algorithm is only 0.3% from the optimal and superior to conventional scheduling policies.


title: On the size-Ramsey number of tight paths

abstract: For any $r\geq 2$ and $k\geq 3$, the $r$-color size-Ramsey number $\hat R(\mathcal{G},r)$ of a $k$-uniform hypergraph $\mathcal{G}$ is the smallest integer $m$ such that there exists a $k$-uniform hypergraph $\mathcal{H}$ on $m$ edges such that any coloring of the edges of $\mathcal{H}$ with $r$ colors yields a monochromatic copy of $\mathcal{G}$. Let $\mathcal{P}_{n,k-1}^{(k)}$ denote the $k$-uniform tight path on $n$ vertices. Dudek, Fleur, Mubayi and R\H{o}dl showed that the size-Ramsey number of tight paths $\hat R(\mathcal{P}_{n,k-1}^{(k)}, 2) = O(n^{k-1-\alpha} (\log n)^{1+\alpha})$ where $\alpha = \frac{k-2}{\binom{k-1}{2}+1}$. In this paper, we improve their bound by showing that $\hat R(\mathcal{P}_{n,k-1}^{(k)}, r) = O(r^k (n\log n)^{k/2})$ for all $k\geq 3$ and $r\geq 2$.


title: Anomalous dynamics triggered by a non-convex equation of state in relativistic flows

abstract: The non-monotonicity of the local speed of sound in dense matter at baryon number densities much higher than the nuclear saturation density ($n_0 \approx 0.16\,$fm$^{-3}$) suggests the possible existence of a non-convex thermodynamics which will lead to a non-convex dynamics. Here, we explore the rich and complex dynamics that an equation of state (EoS) with non-convex regions in the pressure-density plane may develop as a result of genuinely relativistic effects, without a classical counterpart. To this end, we have introduced a phenomenological EoS, whose parameters can be restricted heeding to causality and thermodynamic stability constraints. This EoS shall be regarded as a toy-model with which we may mimic realistic (and far more complex) EoS of practical use in the realm of Relativistic Hydrodynamics.


title: Social Emotion Mining Techniques for Facebook Posts Reaction Prediction

abstract: As of February 2016 Facebook allows users to express their experienced emotions about a post by using five so-called `reactions'. This research paper proposes and evaluates alternative methods for predicting these reactions to user posts on public pages of firms/companies (like supermarket chains). For this purpose, we collected posts (and their reactions) from Facebook pages of large supermarket chains and constructed a dataset which is available for other researches. In order to predict the distribution of reactions of a new post, neural network architectures (convolutional and recurrent neural networks) were tested using pretrained word embeddings. Results of the neural networks were improved by introducing a bootstrapping approach for sentiment and emotion mining on the comments for each post. The final model (a combination of neural network and a baseline emotion miner) is able to predict the reaction distribution on Facebook posts with a mean squared error (or misclassification rate) of 0.135.


title: Merger-Induced Metallicity Dilution in Cosmological Galaxy Formation Simulations

abstract: Observational studies have revealed that galaxy pairs tend to have lower gas-phase metallicity than isolated galaxies. This metallicity deficiency can be caused by inflows of low-metallicity gas due to the tidal forces and gravitational torques associated with galaxy mergers, diluting the metal content of the central region. In this work we demonstrate that such metallicity dilution occurs in state-of-the-art cosmological simulations of galaxy formation. We find that the dilution is typically 0.1 dex for major mergers, and is noticeable at projected separations smaller than $40$ kpc. For minor mergers the metallicity dilution is still present, even though the amplitude is significantly smaller. Consistent with previous analysis of observed galaxies we find that mergers are outliers from the \emph{fundamental metallicity relation}, with deviations being larger than expected for a Gaussian distribution of residuals. Our large sample of mergers within full cosmological simulations also makes it possible to estimate how the star formation rate enhancement and gas consumption timescale behave as a function of the merger mass ratio. We confirm that strong starbursts are likely to occur in major mergers, but they can also arise in minor mergers if more than two galaxies are participating in the interaction, a scenario that has largely been ignored in previous work based on idealised isolated merger simulations.


title: Short Proofs for Slow Consistency

abstract: Let $\operatorname{Con}(\mathbf T)\!\restriction\!x$ denote the finite consistency statement "there are no proofs of contradiction in $\mathbf T$ with $\leq x$ symbols". For a large class of natural theories $\mathbf T$, Pudl\'ak has shown that the lengths of the shortest proofs of $\operatorname{Con}(\mathbf T)\!\restriction\!n$ in the theory $\mathbf T$ itself are bounded by a polynomial in $n$. At the same time he conjectures that $\mathbf T$ does not have polynomial proofs of the finite consistency statements $\operatorname{Con}(\mathbf T+\operatorname{Con}(\mathbf T))\!\restriction\!n$. In contrast we show that Peano arithmetic ($\mathbf{PA}$) has polynomial proofs of $\operatorname{Con}(\mathbf{PA}+\operatorname{Con}^*(\mathbf{PA}))\!\restriction\!n$, where $\operatorname{Con}^*(\mathbf{PA})$ is the slow consistency statement for Peano arithmetic, introduced by S.-D. Friedman, Rathjen and Weiermann. We also obtain a new proof of the result that the usual consistency statement $\operatorname{Con}(\mathbf{PA})$ is equivalent to $\varepsilon_0$ iterations of slow consistency. Our argument is proof-theoretic, while previous investigations of slow consistency relied on non-standard models of arithmetic.


title: Higher Integrability for Constrained Minimizers of Integral Functionals with (p,q)-Growth in low dimension

abstract: We prove higher summability for the gradient of minimizers of strongly convex integral functionals of the Calculus of Variations with (p,q)-Growth conditions in low dimension. Our procedure is set in the framework of Fractional Sobolev Spaces and renders the desired regularity as the result of an approximation technique relying on estimates obtained through a careful use of difference quotients.


title: Estimation of the mechanical properties of the eye through the study of its vibrational modes

abstract: Measuring the eye's mechanical properties in vivo and with minimally invasive techniques can be the key for individualized solutions to a number of eye pathologies. The development of such techniques largely relies on a computational modelling of the eyeball and, it optimally requires the synergic interplay between experimentation and numerical simulation. In Astrophysics and Geophysics the remote measurement of structural properties of the systems of their realm is performed on the basis of (helio-)seismic techniques. As a biomechanical system, the eyeball possesses normal vibrational modes encompassing rich information about its structure and mechanical properties. However, the integral analysis of the eyeball vibrational modes has not been performed yet. Here we develop a new finite difference method to compute both the spheroidal and, specially, the toroidal eigenfrequencies of the human eye. Using this numerical model, we show that the vibrational eigenfrequencies of the human eye fall in the interval 100 Hz - 10 MHz. We find that compressible vibrational modes may release a trace on high frequency changes of the intraocular pressure, while incompressible normal modes could be registered analyzing the scattering pattern that the motions of the vitreous humour leave on the retina. Existing contact lenses with embebed devices operating at high sampling frequency could be used to register the microfluctuations of the eyeball shape we obtain. We advance that an inverse problem to obtain the mechanical properties of a given eye (e.g., Young's modulus, Poisson ratio) measuring its normal frequencies is doable. These measurements can be done using non-invasive techniques, opening very interesting perspectives to estimate the mechanical properties of eyes in vivo. Future research might relate various ocular pathologies with anomalies in measured vibrational frequencies of the eye.


title: Accelerating K-mer Frequency Counting with GPU and Non-Volatile Memory

abstract: The emergence of Next Generation Sequencing (NGS) platforms has increased the throughput of genomic sequencing and in turn the amount of data that needs to be processed, requiring highly efficient computation for its analysis. In this context, modern architectures including accelerators and non-volatile memory are essential to enable the mass exploitation of these bioinformatics workloads. This paper presents a redesign of the main component of a state-of-the-art reference-free method for variant calling, SMUFIN, which has been adapted to make the most of GPUs and NVM devices. SMUFIN relies on counting the frequency of \textit{k-mers} (substrings of length $k$) in DNA sequences, which also constitutes a well-known problem for many bioinformatics workloads, such as genome assembly. We propose techniques to improve the efficiency of k-mer counting and to scale-up workloads like \sm that used to require 16 nodes of \mn to a single machine with a GPU and NVM drives. Results show that although the single machine is not able to improve the time to solution of 16 nodes, its CPU time is 7.5x shorter than the aggregate CPU time of the 16 nodes, with a reduction in energy consumption of 5.5x.


title: Painting galaxies into dark matter halos using machine learning

abstract: We develop a machine learning (ML) framework to populate large dark matter-only simulations with baryonic galaxies. Our ML framework takes input halo properties including halo mass, environment, spin, and recent growth history, and outputs central galaxy and halo baryonic properties including stellar mass ($M_*$), star formation rate (SFR), metallicity ($Z$), neutral ($\rm HI$) and molecular ($\rm H_2$) hydrogen mass. We apply this to the MUFASA cosmological hydrodynamic simulation, and show that it recovers the mean trends of output quantities with halo mass highly accurately, including following the sharp drop in SFR and gas in quenched massive galaxies. However, the scatter around the mean relations is under-predicted. Examining galaxies individually, at $z=0$ the stellar mass and metallicity are accurately recovered ($\sigma\lesssim 0.2$~dex), but SFR and $\rm HI$ show larger scatter ($\sigma\gtrsim 0.3$~dex); these values improve somewhat at $z=1,2$. Remarkably, ML quantitatively recovers second parameter trends in galaxy properties, e.g. that galaxies with higher gas content and lower metallicity have higher SFR at a given $M_*$. Testing various ML algorithms, we find that none perform significantly better than the others, nor does ensembling improve performance, likely because none of the algorithms reproduce the large observed scatter around the mean properties. For the random forest algorithm, we find that halo mass and nearby ($\sim 200$~kpc) environment are the most important predictive variables followed by growth history, while halo spin and $\sim$Mpc scale environment are not important. Finally we study the impact of additionally inputting key baryonic properties $M_*$, SFR and $Z$, as would be available e.g. from an equilibrium model, and show that particularly providing the SFR enables $\rm HI$ to be recovered substantially more accurately.


title: ALMA Detection of Interstellar Methoxymethanol (CH$_3$OCH$_2$OH)

abstract: We report the detection of interstellar methoxymethanol (CH$_3$OCH$_2$OH) in ALMA Bands 6 and 7 toward the MM1 core in the high-mass star-forming region NGC 6334I at ~0.1" - 1" spatial resolution. A column density of 4(2) x $10^{18}$ cm$^{-2}$ at $T_{ex}$ = 200 K is derived toward MM1, ~34 times less abundant than methanol (CH$_3$OH), and significantly higher than predicted by astrochemical models. Probable formation and destruction pathways are discussed, primarily through the reaction of the CH$_3$OH photodissociation products, the methoxy (CH$_3$O) and hydroxymethyl (CH$_2$OH) radicals. Finally, we comment on the implications of these mechanisms on gas-phase vs grain-surface routes operative in the region, and the possibility of electron-induced dissociation of CH$_3$OH rather than photodissociation.


title: Transformational Sparse Coding

abstract: A fundamental problem faced by object recognition systems is that objects and their features can appear in different locations, scales and orientations. Current deep learning methods attempt to achieve invariance to local translations via pooling, discarding the locations of features in the process. Other approaches explicitly learn transformed versions of the same feature, leading to representations that quickly explode in size. Instead of discarding the rich and useful information about feature transformations to achieve invariance, we argue that models should learn object features conjointly with their transformations to achieve equivariance. We propose a new model of unsupervised learning based on sparse coding that can learn object features jointly with their affine transformations directly from images. Results based on learning from natural images indicate that our approach matches the reconstruction quality of traditional sparse coding but with significantly fewer degrees of freedom while simultaneously learning transformations from data. These results open the door to scaling up unsupervised learning to allow deep feature+transformation learning in a manner consistent with the ventral+dorsal stream architecture of the primate visual cortex.


title: Equidistribution of Farey sequences on horospheres in covers of SL(n+1,Z)\SL(n+1,R) and applications

abstract: We establish the limiting distribution of certain subsets of Farey sequences, i.e., sequences of primitive rational points, on expanding horospheres in covers $\Delta\backslash\mathrm{SL}(n+1,\mathbb{R})$ of $\mathrm{SL}(n+1,\mathbb{Z})\backslash\mathrm{SL}(n+1,\mathbb{R})$, where $\Delta$ is a finite index subgroup of $\mathrm{SL}(n+1,\mathbb{Z})$. These subsets can be obtained by projecting to the hyperplane $\{(x_1,\ldots,x_{n+1})\in\mathbb{R}^{n+1}:x_{n+1}=1\}$ sets of the form $\mathbf{A}=\bigcup_{j=1}^J\boldsymbol{a}_j\Delta$, where for all $j$, $\boldsymbol{a}_j$ is a primitive lattice point in $\mathbb{Z}^{n+1}$. Our method involves applying the equidistribution of expanding horospheres in quotients of $\mathrm{SL}(n+1,\mathbb{R})$ developed by Marklof and Str\"{o}mbergsson, and more precisely understanding how the full Farey sequence distributes in $\Delta\backslash\mathrm{SL}(n+1,\mathbb{R})$ when embedded on expanding horospheres as done in previous work by Marklof. For each of the Farey sequence subsets, we extend the statistical results by Marklof regarding the full multidimensional Farey sequences, and solutions by Athreya and Ghosh to Diophantine approximation problems of Erd\H{o}s-Sz\"{u}sz-Tur\'an and Kesten. We also prove that Marklof's result on the asymptotic distribution of Frobenius numbers holds for sets of primitive lattice points of the form $\mathbf{A}$.


title: Development and cytotoxic response of two proliferative MDA- MB-231 and non-proliferative SUM1315 three-dimensional cell culture models of triple-negative basal-like breast cancer cell lines

abstract: Triple-Negative Basal-Like tumors, representing 15 to 20% of breast cancers, are very aggressive and with poor prognosis. Targeted therapies have been developed extensively in preclinical and clinical studies to open the way for new treatment strategies. The present study has focused on developing 3D cell cultures from SUM1315 and MDA-MB-231, two triple-negative basal-like (TNBL) breast cancer cell lines, using the liquid overlay technique. Extracellular matrix concentration, cell density, proliferation, cell viability, topology and ultrastructure parameters were determined. The results showed that for both cell lines, the best conditioning regimen for compact and homogeneous spheroid formation was to use 1000 cells per well and 2% Geltrex. This conditioning regimen highlighted two 3D cell models: non-proliferative SUM1315 spheroids and proliferative MDA-MB-231 spheroids. In both cell lines, the comparison of 2D vs 3D cell culture viability in the presence of increasing concentrations of chemotherapeutic agents i.e. cisplatin, docetaxel and epirubicin, showed that spheroids were clearly less sensitive than monolayer cell cultures. Moreover, a proliferative or non-proliferative 3D cell line property would enable determination of cytotoxic and/or cytostatic drug activity. 3D cell culture could be an excellent tool in addition to the arsenal of techniques currently used in preclinical studies. http://www.impactjournals.com/oncotarget/ Oncotarget, Advance Publications 2017


title: Unconditional stability of semi-implicit discretizations of singular flows

abstract: A popular and efficient discretization of evolutions involving the singular $p$-Laplace operator is based on a factorization of the differential operator into a linear part which is treated implicitly and a regularized singular factor which is treated explicitly. It is shown that an unconditional energy stability property for this semi-implicit time stepping strategy holds. Related error estimates depend critically on a required regularization parameter. Numerical experiments reveal reduced experimental convergence rates for smaller regularization parameters and thereby confirm that this dependence cannot be avoided in general.


title: Detection of Orbital Angular Momentum with Metasurface at Microwave Band

abstract: An orbital angular momentum (OAM) detection approach at microwave band is proposed. A transmittance function is exploited to model a transmissive metasurface. Then the metasurface is designed to convert an OAM wave to multiple waves, only one of which is gaussian. The radiation direction of the gaussian wave is distinguishable according to the order of incident OAM. Consequently, by locating the gaussian wave, the incident OAM can be conveniently determined. We use a simple field source to simulate the incident OAM wave in fullwave simulation. It largely simplifies the simulation process when an incident wave carrying OAM is needed. Both numerical and full-wave simulation results are provided to validate our design and they show a good agreement with each other. Then, the metasurface is optimized for high directivity. Our work can provide an efficient and effective way for OAM detection in radio communications.


title: Correlation between Foam-Bubble Size and Drag Coefficient in Hurricane Conditions

abstract: Recently proposed model of foam impact on the air sea drag coefficient Cd has been employed for the estimation of the efficient foam-bubble radius Rb variation with wind speed U10 in hurricane conditions. The model relates Cd (U10) with the efficient roughness length Zeff (U10) represented as a sum of aerodynamic roughness lengths of the foam free and foam covered sea surfaces Zw (U10 ), and Zf (U10) weighted with the foam coverage coefficient. This relation is treated for known phenomenological distributions Cd (U10), Zw (U10) at strong wind speeds as an inverse problem for the efficient roughness parameter of foam-covered sea surface Zf (U10).


title: The optimal multilinear Bohnenblust-Hille constants: a computational solution for the real case

abstract: The Bohnenblust-Hille inequality for $m$-linear forms was proven in 1931 as a generalization of the famous 4/3-Littlewood inequality. The optimal constants (or at least their asymptotic behavior as $m$ grows) is unknown, but significant for applications. A recent result, obtained by Cavalcante, Teixeira and Pellegrino, provides a kind of algorithm, composed by finitely many elementary steps, giving as the final outcome the optimal truncated Bohnenblust-Hille constants of any order. But the procedure of Cavalcante \textit{et al.} has a fairly large number of calculations and computer assistance cannot be avoided. In this paper we present a computational solution to the problem, using the Wolfram Language. We also use this approach to investigate a conjecture raised by Pellegrino and Teixeira, asserting that $C_m=2^{1-1/m}$ for all $m\in\mathbb{N}$ and to reveal interesting unknown facts about the geometry of $B_{\mathcal{L}(^3\mathbb{R}^3)}$.


title: SU(3) Polyakov linear-sigma model: Magnetic properties of QCD matter in thermal and dense medium

abstract: The linear-sigma model, in which information about confining gluons is included through the Polyakov-loop potential (PLSM), is considered in order to perform a systematic study for various magnetic properties of QCD matter under extreme conditions of high temperatures and densities and finite magnetic field strengths. The introduction of magnetic field to the PLSM Lagrangian requires suitable utilization of Landau quantization, modification in the dispersion relations, and momentum-space dimension-reduction. We observed that increasing the magnetic field leads to filling-up lower Landau levels first and decreasing the number of occupied levels. We conclude that the population of Landau levels is most sensitive to the magnetic field and to the quark charges. The influences of finite magnetic field on the temperature dependence of chiral and deconfinement order-parameter(s) are studied. We present estimations for the magnetization, the magnetic susceptibility, the permeability and the catalysis properties of QCD matter as functions of temperature. The dependences of the resulting freezeout parameters, temperatures and baryon chemical potentials on the corresponding magnetic field strengths have been analyzed, as well. These calculations are compared with recent lattice QCD simulations, whenever available. We conclude that the QCD matter seems to have paramagnetic property at temperatures greater than the critical one. There is an evidence for weak diamagnetic property at low temperatures. Last but not least, we observe that the magnetic catalysis is inverse, namely the critical temperatures decrease with increasing the magnetic field.


title: Heat kernel estimates for Dirichlet fractional Laplacian with gradient perturbation

abstract: We give a direct proof of the sharp two-sided estimates, recently established in [4,9], for the Dirichlet heat kernel of the fractional Laplacian with gradient perturbation in $C^{1, 1}$ open sets by using Duhamel formula. We also obtain a gradient estimate for the Dirichlet heat kernel. Our assumption on the open set is slightly weaker in that we only require $D$ to be $C^{1,\theta}$ for some $\theta\in (\alpha/2, 1]$.


title: Reconstructing the invisible with matrix elements

abstract: We propose a fully flexible method to perform an hypothesis test between signal and background based on the Matrix Element Method in the presence of multiple invisible particles. The proposed method performs a mapping of the measured final state onto its minimal hypersurface of degrees of freedom for a given process and then maximises the matrix element on this hypersurface separately for signal and background. To show how performant the method is in separating signal from background, we apply it to the prominent partly invisible decay of a Higgs boson into a muon-antimuon pair and two muon-neutrinos via two W bosons.


title: Performance of Analog Nonlinear Filtering for Impulsive Noise Mitigation in OFDM-based PLC Systems

abstract: Asynchronous and cyclostationary impulsive noise can severely impact the bit-error-rate (BER) of OFDM-based powerline communication systems. In this paper, we analyze an adaptive nonlinear analog front end filter that mitigates various types of impulsive noise without detrimental effects such as self-interference and out-of-band power leakage caused by other nonlinear approaches like clipping and blanking. Our proposed Adaptive Nonlinear Differential Limiter (ANDL) is constructed from a linear analog filter by applying a feedback-based nonlinearity, controlled by a single resolution parameter. We present a simple practical method to find the value of this resolution parameter that ensures the mitigation of impulsive without impacting the desired OFDM signal. Unlike many prior approaches for impulsive noise mitigation that assume a statistical noise model, ANDL is blind to the exact nature of the noise distribution, and is designed to be fully compatible with existing linear front end filters. We demonstrate the potency of ANDL by simulating the OFDM-based narrowband PLC compliant with the IEEE standards. We show that the proposed ANDL outperforms other approaches in reducing the BER in impulsive noise environments.


title: On Freud's superpotential in General Relativity and in Einstein-Cartan theory

abstract: The identification of a suitable gravitational energy in theories of gravity has a long history, and it is well known that a unique answer cannot be given. In the first part of this paper we present a streamlined version of the derivation of Freud's superpotential in general relativity. It is found if we once integrate the gravitational field equation by parts. This allows us to extend these results directly to the Einstein-Cartan theory. Interestingly, Freud's original expression, first stated in 1939, remains valid even when considering gravitational theories in Riemann-Cartan or, more generally, in metric-affine spacetimes.


title: Numerical methods for thermally stressed shallow shell equations

abstract: We develop efficient and accurate numerical methods to solve a class of shallow shell problems of the von Karman type. The governing equations form a fourth-order coupled system of nonlinear biharnomic equations for the transverse deflection and Airy's stress function. A second-order finite difference discretization with three iterative methods (Picard, Newton and Trust-Region Dogleg) are proposed for the numerical solution of the nonlinear PDE system. Three simple boundary conditions and two application-motivated mixed boundary conditions are considered. Along with the nonlinearity of the system, boundary singularities that appear when mixed boundary conditions are specified are the main numerical challenges. Two approaches that use either a transition function or local corrections are developed to deal with these boundary singularities. All the proposed numerical methods are validated using carefully designed numerical tests, where expected orders of accuracy and rates of convergence are observed. A rough run-time performance comparison is also conducted to illustrate the efficiency of our methods. As an application of the methods, a snap-through thermal buckling problem is considered. The critical thermal loads of shell buckling with various boundary conditions are numerically calculated, and snap-through bifurcation curves are also obtained using our numerical methods together with a pseudo-arclength continuation method. Our results are consistent with previous studies.


title: Spacetime Singularities vs. Topologies of Zeeman-G\"obel Class

abstract: In this article we first observe that the Path topology of Hawking, King and MacCarthy is an analogue, in curved spacetimes, of a topology that was suggested by Zeeman as an alternative topology to his so-called Fine topology in Minkowski spacetime. We then review a result of a recent paper on spaces of paths and the Path topology, and see that there are at least five more topologies in the class $\mathfrak{Z}-\mathfrak{G}$ of Zeeman-G\"obel topologies which admit a countable basis, incorporate the causal and conformal structures, but the Limit Curve Theorem fails to hold. The "problem" that L.C.T. does not hold can be resolved by "adding back" the light-cones in the basic-open sets of these topologies, and create new basic open sets for new topologies. But, the main question is: do we really need the L.C.T. to hold, and why? Why is the manifold topology, under which the group of homeomorphisms of a spacetime is vast and of no physical significance (Zeeman), more preferable from an appropriate topology in the class $\mathfrak{Z}-\mathfrak{G}$ under which a homeomorphism is an isometry (G\"obel)? Since topological conditions that come as a result of a causality requirement are key in the existence of singularities in general relativity, the global topological conditions that one will supply the spacetime manifold might play an important role in describing the transition from the quantum non-local theory to a classical local theory.


title: Nontrivial temperature behavior of the carrier concentration in the nanostructure "graphene channel on ferroelectric substrate with domain walls"

abstract: This work explores a nontrivial temperature behavior of the carriers concentration, which governs graphene channel conductance in the nano-structure "graphene channel on ferroelectric substrate" that is a basic element for FETs in non-volatile memory units of new generation. We revealed the transition from a single to double antiferroelectric-like hysteresis loop of the concentration voltage dependence that happens with the temperature increase and then exist in a wide temperature range (from 350 to 500 K). Unexpectedly we revealed the double loops of polarization and concentration can have irregular shape that remains irregular as long as the computation takes place, and the voltage position of the different features (jumps, secondary maxima, etc.) changes from one period to another, leading to the impression of quasi-chaotic behavior. It appeared that these effects originate from the nonlinear screening of ferroelectric polarization by graphene carriers, as well as it is conditioned by the temperature evolution of the domain structure kinetics in ferroelectric substrate. The nonlinearity rules the voltage behavior of polarization screening by graphene 2D-layer and at the same time induces the motion of separated domain walls accompanied by the motion of p-n junction along the graphene channel (2D-analog of Hann effect). Since the domain walls structure, period and kinetics can be controlled by varying the temperature, we concluded that the considered nano-structures based on graphene-on-ferroelectric are promising for the fabrication of new generation of modulators based on the graphene p-n junctions.


title: Skyrme N2LO pseudo-potential for calculations of properties of atomic nuclei

abstract: We present recent developments obtained in the so-called N2LO extension of the usual Skyrme pseudo-potential. In particular, we discuss the isovector splitting mass in infinite nuclear matter and the pairing gaps of selected semi-magic even-even nuclei.


title: Near-field refrigeration and tunable heat exchange through four-wave mixing

abstract: We modify and extend a recently proposed four-wave mixing scheme [Opt. Express 25 (19),23164 (2017)] for achieving near-field thermal upconversion and energy transfer, to demonstrate efficient thermal refrigeration at low intensities $\sim 10^{-9}$W/m$^2$ over a wide range of gap sizes (from tens to hundreds of nanometers) and operational temperatures (from tens to hundreds of Kelvins). We further exploit the scheme to achieve magnitude and directional tunability of near-field heat exchange between bodies held at different temperatures.


title: Multidimensional Models of Type Ia Supernova Nebular Spectra: Strong Emission Lines from Stripped Companion Gas Rule Out Classic Single Degenerate Systems

abstract: The classic single-degenerate model for the progenitors of Type Ia Supernova (SN Ia) predicts that the supernova ejecta should be enriched with solar-like abundance material stripped from the companion star. Spectroscopic observations of normal SNe Ia at late times, however, have not resulted in definite detection of hydrogen. In this Letter, we study line formation in SNe Ia at nebular times using non-LTE spectral modeling. We present, for the first time, multidimensional radiative transfer calculations of SNe Ia with stripped material mixed in the ejecta core, based on hydrodynamical simulations of ejecta-companion interaction. We find that interaction models with main sequence companions produce significant H$\alpha$ emission at late times, ruling out this type of binaries being viable progenitors of SNe Ia. We also predict significant He I line emission at optical and near-infrared wavelengths for both hydrogen-rich or helium-rich material, providing an additional observational probe of stripped ejecta. We produce models with reduced stripped masses and find a more stringent mass limit of $M_{st} \lesssim 1\times 10^{-4} M_\odot$ of stripped companion material for SN 2011fe.


title: Improved Search for Heavy Neutrinos in the Decay $\pi\rightarrow e\nu$

abstract: A search for massive neutrinos has been made in the decay $\pi\rightarrow e^+ \nu$. No evidence was found for extra peaks in the positron energy spectrum indicative of pion decays involving massive neutrinos ($\pi\rightarrow e^+ \nu_h$). Upper limits (90 \% C.L.) on the neutrino mixing matrix element $|U_{ei}|^2$ in the neutrino mass region 60--135 MeV/$c^2$ were set, which are %representing an order of magnitude improvement over previous results.


title: Linear stability analysis of detonations via numerical computation and dynamic mode decomposition

abstract: We introduce a new method to investigate linear stability of gaseous detonations that is based on an accurate shock-fitting numerical integration of the linearized reactive Euler equations with a subsequent analysis of the computed solution via the dynamic mode decomposition. The method is applied to the detonation models based on both the standard one-step Arrhenius kinetics and two-step exothermic-endothermic reaction kinetics. Stability spectra for all cases are computed and analyzed. The new approach is shown to be a viable alternative to the traditional normal-mode analysis used in detonation theory.


title: On Bipartite Operators Defined by Completely Different Permutations

abstract: We introduce a class of bipartite operators acting on $\mathcal{H} \otimes \mathcal{H}$ ($\mathcal{H}$ being an $n$-dimensional Hilbert space) defined by a set of $n$ Completely Different Permutations CDP. Bipartite operators are of particular importance in quantum information theory to represent states and observables of composite quantum systems. It turns out that any set of CDPs gives rise to a certain direct sum decomposition of the total Hilbert space which enables one simple construction of the corresponding bipartite operator. Interestingly, if set of CDPs defines an abelian group then the corresponding bipartite operator displays an additional property -- the partially transposed operator again corresponds to (in general different) set of CDPs. Therefore, our technique may be used to construct new classes of so called PPT states which are of great importance for quantum information. Using well known relation between bipartite operators and linear maps one use also construct linear maps related to CDPs.


title: Using Black-box Compression Algorithms for Phase Retrieval

abstract: Compressive phase retrieval refers to the problem of recovering a structured $n$-dimensional complex-valued vector from its phase-less under-determined linear measurements. The non-linearity of measurements makes designing theoretically-analyzable efficient phase retrieval algorithms challenging. As a result, to a great extent, algorithms designed in this area are developed to take advantage of simple structures such as sparsity and its convex generalizations. The goal of this paper is to move beyond simple models through employing compression codes. Such codes are typically developed to take advantage of complex signal models to represent the signals as efficiently as possible. In this work, it is shown how an existing compression code can be treated as a black box and integrated into an efficient solution for phase retrieval. First, COmpressive PhasE Retrieval (COPER) optimization, a computationally-intensive compression-based phase retrieval method, is proposed. COPER provides a theoretical framework for studying compression-based phase retrieval. The number of measurements required by COPER is connected to $\kappa$, the $\alpha$-dimension (closely related to the rate-distortion dimension) of the given family of compression codes. To finds the solution of COPER, an efficient iterative algorithm called gradient descent for COPER (GD-COPER) is proposed. It is proven that under some mild conditions on the initialization, if the number of measurements is larger than $ C \kappa^2 \log^2 n$, where $C$ is a constant, GD-COPER obtains an accurate estimate of the input vector in polynomial time. In the simulation results, JPEG2000 is integrated in GD-COPER to confirm the superb performance of the resulting algorithm on real-world images.


title: Pseudochemotaxis in inhomogeneous active Brownian systems

abstract: We study dynamical properties of confined, self-propelled Brownian particles in an inhomogeneous activity profile. Using Brownian dynamics simulations, we calculate the probability to reach a fixed target and the mean first passage time to the target of an active particle. We show that both these quantities are strongly influenced by the inhomogeneous activity. When the activity is distributed such that high-activity zone is located between the target and the starting location, the target finding probability is increased and the passage time is decreased in comparison to a uniformly active system. Moreover, for a continuously distributed profile, the activity gradient results in a drift of active particle up the gradient bearing resemblance to chemotaxis. Integrating out the orientational degrees of freedom, we derive an approximate Fokker-Planck equation and show that the theoretical predictions are in very good agreement with the Brownian dynamics simulations.


title: Nintendo Super Smash Bros. Melee: An "Untouchable" Agent

abstract: Nintendo's Super Smash Bros. Melee fighting game can be emulated on modern hardware allowing us to inspect internal memory states, such as character positions. We created an AI that avoids being hit by training using these internal memory states and outputting controller button presses. After training on a month's worth of Melee matches, our best agent learned to avoid the toughest AI built into the game for a full minute 74.6% of the time.


title: Fast Low-Rank Matrix Estimation without the Condition Number

abstract: In this paper, we study the general problem of optimizing a convex function $F(L)$ over the set of $p \times p$ matrices, subject to rank constraints on $L$. However, existing first-order methods for solving such problems either are too slow to converge, or require multiple invocations of singular value decompositions. On the other hand, factorization-based non-convex algorithms, while being much faster, require stringent assumptions on the \emph{condition number} of the optimum. In this paper, we provide a novel algorithmic framework that achieves the best of both worlds: asymptotically as fast as factorization methods, while requiring no dependency on the condition number. We instantiate our general framework for three important matrix estimation problems that impact several practical applications; (i) a \emph{nonlinear} variant of affine rank minimization, (ii) logistic PCA, and (iii) precision matrix estimation in probabilistic graphical model learning. We then derive explicit bounds on the sample complexity as well as the running time of our approach, and show that it achieves the best possible bounds for both cases. We also provide an extensive range of experimental results, and demonstrate that our algorithm provides a very attractive tradeoff between estimation accuracy and running time.


title: Relativistic anisotropic hydrodynamics

abstract: In this paper we review recent progress in relativistic anisotropic hydrodynamics. We begin with a pedagogical introduction to the topic which takes into account the advances in our understanding of this topic since its inception. We consider both conformal and non-conformal systems and demonstrate how one can implement a realistic equation of state using a quasiparticle approach. We then consider the inclusion of non-spheroidal (non-ellipsoidal) corrections to leading-order anisotropic hydrodynamics and present the findings of the resulting second-order viscous anisotropic hydrodynamics framework. We compare the results obtained in both the conformal and non-conformal cases with exact solutions to the Boltzmann equation and demonstrate that, in all known cases, anisotropic hydrodynamics best reproduces the exact solutions. Based on this success, we then discuss the phenomenological application of anisotropic hydrodynamics. Along these lines, we review techniques which can be used to convert a momentum-space anisotropic fluid into hadronic degrees of freedom by generalizing the original idea of Cooper-Frye freeze-out to momentum-space anisotropic systems. And, finally, we present phenomenological results of 3+1d quasiparticle anisotropic hydrodynamic simulations and compare them to experimental data produced in 2.76 TeV Pb-Pb collisions at the LHC. Our results indicate that anisotropic hydrodynamics provides a promising framework for describing the dynamics of the momentum-space anisotropic QGP created in heavy-ion collisions.


title: Spatio-temporal complexity in dual delay nonlinear laser dynamics: chimeras and dissipative solitons

abstract: We demonstrate for a nonlinear photonic system that two highly asymmetric feedback delays can induce a variety of emergent patterns which are highly robust during the system's global evolution. Explicitly, two-dimensional chimeras and dissipative solitons become visible upon a space-time transformation. Switching between chimeras and dissipative solitons requires only adjusting two system parameters, demonstrating self-organization exclusively based on the system's dynamical properties. Experiments were performed using a tunable semiconductor laser's transmission through a Fabry-Perot resonator resulting in an Airy function as nonlinearity. Resulting dynamics were band-pass filtered and propagated along two feedback paths whose time delays differ by two orders of magnitude. An excellent agreement between experimental results and theoretical model given by modified Ikeda equations was achieved.


title: Followup Observations of SDSS and CRTS Candidate Cataclysmic Variables II

abstract: Spectra of 38 candidate or known cataclysmic variables are presented. Most are candidate dwarf novae or systems containing possible highly magnetic white dwarfs, while a few (KR Aur, LS Peg, V380 Oph and V694 Mon) are previously known objects caught in unusual states. Individual spectra are used to confirm a dwarf nova nature or other classification while radial velocities of 15 systems provide orbital periods and velocity amplitudes that aid in determining the nature of the objects. Our results substantiate a polar nature for four objects, find an eclipsing SW Sex star below the period gap, another as a likely intermediate polar, as well as two dwarf novae with periods in the middle of the gap.


title: Non-commutative massive unquenched ABJM

abstract: In this paper we study non-commutative massive unquenched Chern-Simons matter theory using its gravity dual. We construct this novel background by applying a TsT-transformation on the known parent commutative solution. We discuss several aspects of this solution to the Type IIA supergravity equations of motion and, amongst others, check that it preserves ${\cal N}=1$ supersymmetry. We then turn our attention to applications and investigate how dynamical flavor degrees of freedom affect numerous observables of interest. Our framework can be regarded as a key step towards the construction of holographic quantum Hall states on a non-commutative plane.


title: Effects of energy conservation on equilibrium properties of hot asymmetric nuclear matter

abstract: Based on the relativistic Vlasov-Uehling-Uhlenbeck transport model, which includes relativistic scalar and vector potentials on baryons, we consider a $N-\Delta-\pi$ system in a box with periodic boundary conditions to study the effects of energy conservation in particle production and absorption processes on the equilibrium properties of the system. The density and temperature of the matter in the box are taken to be similar to the hot dense matter formed in heavy ion collisions at intermediate energies. We find that to maintain the equilibrium numbers of $N$, $\Delta$ and $\pi$, which depend on the mean-field potentials of $N$ and $\Delta$, requires the inclusion of these potentials in the energy conservation condition that determines the momenta of outgoing particles after a scattering or decay process. We further find that the baryon scalar potentials mainly affect the $\Delta$ and pion equilibrium numbers, while the baryon vector potentials have considerable effects on the effective charged pion ratio at equilibrium. Our results thus indicate that it is essential to include in the transport model the effect of potentials in the energy conservation of a scattering or decay process, which is ignored in most transport models, for studying pion production in heavy ion collisions.


title: Shock electromagnetic wave of large amplitude in the nonlinear transmission line based on distributed semiconductor diode

abstract: An analytical solution to the problem of formation and propagation of large amplitude shock electromagnetic waves in strip transmission line (TL) based on the distributed semiconductor diode was obtained. This solution correctly takes into account non-linearity, time dispersion and dissipation of TL. The results of the theory has been used to estimate the parameters of TL as peaker of the high voltage pulse applied to transmission line input.


title: The Evens and Odds of CMB Anomalies

abstract: The lack of power of large--angle CMB anisotropies is known to increase its statistical significance at higher Galactic latitudes, where a string--inspired pre--inflationary scale $\Delta$ can also be detected. Considering the Planck 2015 data, and relying largely on a Bayesian approach, we show that the effect is mostly driven by the \emph{even}--$\ell$ harmonic multipoles with $\ell \lesssim 20$, which appear sizably suppressed in a way that is robust with respect to Galactic masking, along with the corresponding detections of $\Delta$. On the other hand, the first \emph{odd}--$\ell$ multipoles are only suppressed at high Galactic latitudes. We investigate this behavior in different sky masks, constraining $\Delta$ through even and odd multipoles, and we elaborate on possible implications. We include low--$\ell$ polarization data which, despite being noise--limited, help in attaining confidence levels of about 3 $\sigma$ in the detection of $\Delta$. We also show by direct forecasts that a future all--sky $E$--mode cosmic--variance--limited polarization survey may push the constraining power for $\Delta$ beyond 5 $\sigma$.


title: Chaotic bursting in semiconductor lasers

abstract: We investigate the dynamic mechanisms for low frequency fluctuations in semiconductor lasers subject to delayed optical feedback, using the Lang-Kobayashi model. This system of delay differential equations displays pronounced envelope dynamics, ranging from erratic, so called low frequency fluctuations to regular pulse packages, if the time scales of fast oscillations and envelope dynamics are well separated. We investigate the parameter regions where low frequency fluctuations occur and compute their Lyapunov spectrum. Using geometric singular perturbation theory, we study this intermittent chaotic behavior and characterize these solutions as bursting slow-fast oscillations.


title: Device-to-Device Networking Meets Cellular via Network Coding

abstract: Utilizing device-to-device (D2D) connections among mobile devices is promising to meet the increasing throughput demand over cellular links. In particular, when mobile devices are in close proximity of each other and are interested in the same content, D2D connections such as Wi-Fi Direct can be opportunistically used to construct a cooperative (and jointly operating) cellular and D2D networking system. However, it is crucial to understand, quantify, and exploit the potential of network coding for cooperating mobile devices in the joint cellular and D2D setup. In this paper, we consider this problem, and (i) develop a network coding framework, namely NCMI, for cooperative mobile devices in the joint cellular and D2D setup, where cellular and D2D link capacities are the same, and (ii) characterize the performance of the proposed network coding framework, where we use packet completion time, which is the number of transmission slots to recover all packets, as a performance metric. We demonstrate the benefits of our network coding framework through simulations.


title: Input-to-State Stability of Periodic Orbits of Systems with Impulse Effects via Poincar\'e Analysis

abstract: In this paper we investigate the relation between robustness of periodic orbits exhibited by systems with impulse effects and robustness of their corresponding Poincar\'e maps. In particular, we prove that input-to-state stability (ISS) of a periodic orbit under external excitation in both continuous and discrete time is equivalent to ISS of the corresponding 0-input fixed point of the associated \emph{forced} Poincar\'e map. This result extends the classical Poincar\'e analysis for asymptotic stability of periodic solutions to establish orbital input-to-state stability of such solutions under external excitation. In our proof, we define the forced Poincar\'e map, and use it to construct ISS estimates for the periodic orbit in terms of ISS estimates of this map under mild assumptions on the input signals. As a consequence of the availability of these estimates, the equivalence between exponential stability (ES) of the fixed point of the 0-input (unforced) Poincar\'e map and ES of the corresponding orbit is recovered. The results can be applied naturally to study the robustness of periodic orbits of continuous-time systems as well. Although our motivation for extending classical Poincar\'e analysis to address ISS stems from the need to design robust controllers for limit-cycle walking and running robots, the results are applicable to a much broader class of systems that exhibit periodic solutions.


title: Evaluating the Robustness of Rogue Waves Under Perturbations

abstract: Rogue waves, and their periodic counterparts, have been shown to exist in a number of integrable models. However, relatively little is known about the existence of these objects in models where an exact formula is unattainable. In this work, we develop a novel numerical perspective towards identifying such states as localized solutions in space-time. Importantly, we illustrate that this methodology in addition to benchmarking known solutions (and confirming their numerical propagation under controllable error) enables the continuation of such solutions over parametric variations to non-integrable models. As a result, we can answer in the positive the question about the parametric robustness of Peregrine-like waveforms and even of generalizations thereof on a cnoidal wave background.


title: Experimental and Numerical Investigation of Turbulent Jets Issuing Through a Realistic Pipeline Geometry: Asymmetry Effects for Air, Helium, and Hydrogen

abstract: Experiments and numerical simulations were conducted to investigate the dispersion of turbulent jets issuing from realistic pipe geometries. The effect of jet densities and Reynolds numbers on vertical buoyant jets were investigated, as they emerged from the side wall of a circular pipe, through a round orifice. Particle image velocimetry (PIV) and planar laser-induced fluorescence (PLIF) techniques were employed simultaneously to provide time-averaged flow velocity and concentrations fields. Large eddy simulation (LES) was applied to provide further detail with regards to the three-dimensionality of air, helium, and hydrogen jets. These realistic jets were always asymmetric and found to deflect about the vertical axis. This deflection was influenced by buoyancy, where heavier gases deflected more than lighter gases. Significant turbulent mixing was also observed in the near field. The realistic jets, therefore, experienced faster velocity decay, and asymmetric jet spreading compared to round jets. These findings indicate that conventional round jet assumptions are, to some extent, inadequate to predict gas concentration, entrainment rates and, consequently, the extent of the flammability envelope of realistic gas leaks.


title: Toward computing gravitational initial data without elliptic solvers

abstract: Two new methods have been proposed for solving the gravitational constraints without using elliptic solvers by formulating them as either an algebraic-hyperbolic or parabolic-hyperbolic system. Here, we compare these two methods and present a unified computational infrastructure for their implementation as numerical evolution codes. An important potential application of these methods is the prescription of initial data for the simulation of black holes. This paper is meant to support progress and activity in that direction.


title: Momentum-Space Cluster Dual Fermion Method

abstract: Recent years have seen the development of two types of non-local extensions to the single-site dynamical mean field theory. On one hand, cluster approximations, such as the dynamical cluster approximation, recover short-range momentum-dependent correlations non-perturbatively. On the other hand, diagrammatic extensions, such as the dual fermion theory, recover long ranged corrections perturbatively. The correct treatment of both strong short-ranged and weak long-ranged correlations within the same framework is therefore expected to lead to a quick convergence of results, and offers the potential of obtaining smooth self-energies in non-perturbative regimes of phase space. In this paper, we present an exact cluster dual fermion method based on an expansion around the dynamical cluster approximation. Unlike previous formulations, our method does not employ a coarse graining approximation to the interaction, which we show to be the leading source of error at high temperature, and converges to the exact result independent of the size of the underlying cluster. We illustrate the power of the method with results for the second-order cluster dual fermion approximation to the single-particle self-energies and double occupancies.


title: Nonparametric Composite Hypothesis Testing in an Asymptotic Regime

abstract: We investigate the nonparametric, composite hypothesis testing problem for arbitrary unknown distributions in the asymptotic regime where both the sample size and the number of hypotheses grow exponentially large. Such asymptotic analysis is important in many practical problems, where the number of variations that can exist within a family of distributions can be countably infinite. We introduce the notion of \emph{discrimination capacity}, which captures the largest exponential growth rate of the number of hypotheses relative to the sample size so that there exists a test with asymptotically vanishing probability of error. Our approach is based on various distributional distance metrics in order to incorporate the generative model of the data. We provide analyses of the error exponent using the maximum mean discrepancy (MMD) and Kolmogorov-Smirnov (KS) distance and characterize the corresponding discrimination rates, i.e., lower bounds on the discrimination capacity, for these tests. Finally, an upper bound on the discrimination capacity based on Fano's inequality is developed. Numerical results are presented to validate the theoretical results.


title: On the Longest Spanning Tree with Neighborhoods

abstract: We study a maximization problem for geometric network design. Given a set of $n$ compact neighborhoods in $\mathbb{R}^d$, select a point in each neighborhood, so that the longest spanning tree on these points (as vertices) has maximum length. Here we give an approximation algorithm with ratio $0.511$, which represents the first, albeit small, improvement beyond $1/2$. While we suspect that the problem is NP-hard already in the plane, this issue remains open.


title: Neumann Optimizer: A Practical Optimization Algorithm for Deep Neural Networks

abstract: Progress in deep learning is slowed by the days or weeks it takes to train large models. The natural solution of using more hardware is limited by diminishing returns, and leads to inefficient use of additional resources. In this paper, we present a large batch, stochastic optimization algorithm that is both faster than widely used algorithms for fixed amounts of computation, and also scales up substantially better as more computational resources become available. Our algorithm implicitly computes the inverse Hessian of each mini-batch to produce descent directions; we do so without either an explicit approximation to the Hessian or Hessian-vector products. We demonstrate the effectiveness of our algorithm by successfully training large ImageNet models (Inception-V3, Resnet-50, Resnet-101 and Inception-Resnet-V2) with mini-batch sizes of up to 32000 with no loss in validation error relative to current baselines, and no increase in the total number of steps. At smaller mini-batch sizes, our optimizer improves the validation error in these models by 0.8-0.9%. Alternatively, we can trade off this accuracy to reduce the number of training steps needed by roughly 10-30%. Our work is practical and easily usable by others -- only one hyperparameter (learning rate) needs tuning, and furthermore, the algorithm is as computationally cheap as the commonly used Adam optimizer.


title: Posterior distribution existence and error control in Banach spaces in the Bayesian approach to UQ in inverse problems

abstract: We generalize the results of \cite{Capistran2016} on expected Bayes factors (BF) to control the numerical error in the posterior distribution to an infinite dimensional setting when considering Banach functional spaces and now in a prior setting. The main result is a bound on the absolute global error to be tolerated by the Forward Map numerical solver, to keep the BF of the numerical vs. the theoretical model near to 1, now in this more general setting, possibly including a truncated, finite dimensional approximate prior measure. In so doing we found a far more general setting to define and prove existence of the infinite dimensional posterior distribution than that depicted in, for example, \cite{Stuart2010}. Discretization consistency and rates of convergence are also investigated in this general setting for the Bayesian inverse problem.


title: Weak Fraisse categories

abstract: We develop the theory of weak Fraisse categories, where the crucial concept is the weak amalgamation property, discovered relatively recently in model theory. We show that, in a suitable framework, every weak Fraisse category has its unique limit, a special object in a bigger category, characterized by certain variant of injectivity. This significantly extends the present theory of Fraisse limits.