title: Probing the scalar potential via double Higgs boson production at hadron colliders

abstract: We present a sensitivity study on the cubic and quartic self couplings in double Higgs production via gluon fusion at hadron colliders. Considering the relevant operators in the Standard Model Effective Field Theory up to dimension eight, we calculate the dominant contributions up to two-loop level, where the first dependence on the quartic interaction appears. Our approach allows to study the independent variations of the two self couplings and to clearly identify the terms necessary to satisfy gauge invariance and to obtain UV-finite results order by order in perturbation theory. We focus on the $b \bar b \gamma \gamma$ signature for simplicity and provide the expected bounds for the cubic and quartic self couplings at the 14 TeV LHC with 3000 fb$^{-1}$ (HL-LHC) and for a future 100 TeV collider (FCC-100) with 30 ab$^{-1}$. We find that while the HL-LHC will provide very limited sensitivity on the quartic self coupling, precision measurements of double Higgs production at a FCC-100 will offer the opportunity to set competitive bounds. We show that combining information from double and triple Higgs production leads to significantly improved prospects for the determination of the quartic self coupling.


title: Quantum Corrections to Central Charges and Supersymmetric Casimir Energy in AdS$_3$/CFT$_2$

abstract: We study the Casimir energy of bulk fields in AdS$_3$ and its relation to subleading terms in the central charge of the dual CFT$_2$. Computing both sides of the standard CFT$_2$ relation $E=-c/12$ independently we show that this relation is not necessarily satisfied at the level of individual bulk supergravity states, but in theories with sufficient supersymmetry it is restored at the level of bulk supermultiplets. Assuming only $(0,2)$ supersymmetry (or more), we improve the situation by relating quantum corrections to the central charge and the supersymmetric Casimir energy which in turn is related to an index. These relations adapt recent progress on the AdS$_5$/CFT$_4$ correspondence to AdS$_3$/CFT$_2$ holography. We test our formula successfully in several examples, including the $(0,4)$ MSW theory describing classes of 4D black holes and the large $(4,4)$ theory that is interesting for higher spin holography. We also make predictions for the subleading central charges in several recently proposed $(2,2)$ dualities where the CFT$_2$ is not yet well-understood.


title: There is no isolated interface edge in very supercritical percolation

abstract: We consider the Bernoulli bond percolation model in a box $\Lambda$ (not necessarily parallel to the directions of the lattice) in the regime where the percolation parameter is close to $1$. We condition the configuration on the event that two opposite faces of the box are disconnected. We couple this configuration with an unconstrained percolation configuration. The interface edges are the edges which differ in the two configurations. We prove that, typically, each interface edge is within a distance of order $\ln|\Lambda|$ of another interface edge or of a pivotal edge. We derive an estimate for the law of an edge which is far from the cut and the interface edges.


title: Small Hazard-free Transducers

abstract: Ikenmeyer et al. (JACM'19) proved an unconditional exponential separation between the hazard-free complexity and (standard) circuit complexity of explicit functions. This raises the question: which classes of functions permit efficient hazard-free circuits? In this work, we prove that circuit implementations of transducers with small state space are such a class. A transducer is a finite state machine that transcribes, symbol by symbol, an input string of length $n$ into an output string of length $n$. We present a construction that transforms any function arising from a transducer into an efficient circuit of size $\mathcal{O}(n)$ computing the hazard-free extension of the function. More precisely, given a transducer with $s$ states, receiving $n$ input symbols encoded by $l$ bits, and computing $n$ output symbols encoded by $m$ bits, the transducer has a hazard-free circuit of size $2^{\mathcal{O}(s+\ell)} m n$ and depth $\mathcal{O}(s\log n + \ell)$; in particular, if $s, \ell,m\in \mathcal{O}(1)$, size and depth are asymptotically optimal. In light of the strong hardness results by Ikenmeyer et al. (JACM'19), we consider this a surprising result.


title: On the local H\"older boundary smoothness of an analytic function in the unit ball compared with the smoothness of its modulus

abstract: Local boundary smoothness of an analytic function f on the unit ball of C^n is compared to the smoothness of its modulus. We prove that different conditions imposed on the zeros of f imply different drops of the smoothness. We also show that some of the drops are the best possible.


title: Mixed Non-Gaussianity from Axion-Gauge Field Dynamics

abstract: We study scalar-tensor-tensor cross correlation $\langle \zeta hh \rangle$ generated by the dynamics of interacting axion and SU(2) gauge fields during inflation. We quantize the quadratic action and solve the linear equations by taking into account mixing terms in a non-perturbative manner. Combining that with the in-in formalism, we compute contributions from cubic interactions to the bispectrum $B_{\zeta hh}$. We find that the bispectrum is peaked at the folded configuration, which is a unique feature encoded by the scalar mixing and localized production of tensor modes. With our parameter choice, the amplitude of the bispectrum is $k^6 B_{\zeta hh} \sim 10^{-16}$. The unique shape dependence, together with the parity-violating nature, is thus a distinguishing feature to search for in the CMB observables.


title: Numerical modelling of heat transfer and experimental validation in Powder-Bed Fusion with the Virtual Domain Approximation

abstract: Among metal additive manufacturing technologies, powder-bed fusion features very thin layers and rapid solidification rates, leading to long build jobs and a highly localized process. Many efforts are being devoted to accelerate simulation times for practical industrial applications. The new approach suggested here, the virtual domain approximation, is a physics-based rationale for spatial reduction of the domain in the thermal finite-element analysis at the part scale. Computational experiments address, among others, validation against a large physical experiment of 17.5 $\mathrm{[cm^3]}$ of deposited volume in 647 layers. For fast and automatic parameter estimation at such level of complexity, a high-performance computing framework is employed. It couples FEMPAR-AM, a specialized parallel finite-element software, with Dakota, for the parametric exploration. Compared to previous state-of-the-art, this formulation provides higher accuracy at the same computational cost. This sets the path to a fully virtualized model, considering an upwards-moving domain covering the last printed layers.


title: Diverse Image Synthesis from Semantic Layouts via Conditional IMLE

abstract: Most existing methods for conditional image synthesis are only able to generate a single plausible image for any given input, or at best a fixed number of plausible images. In this paper, we focus on the problem of generating images from semantic segmentation maps and present a simple new method that can generate an arbitrary number of images with diverse appearance for the same semantic layout. Unlike most existing approaches which adopt the GAN framework, our method is based on the recently introduced Implicit Maximum Likelihood Estimation (IMLE) framework. Compared to the leading approach, our method is able to generate more diverse images while producing fewer artifacts despite using the same architecture. The learned latent space also has sensible structure despite the lack of supervision that encourages such behaviour. Videos and code are available at https://people.eecs.berkeley.edu/~ke.li/projects/imle/scene_layouts/.


title: The Ultraviolet Detection of Diffuse Gas in Galaxy Groups

abstract: A small survey of the UV-absorbing gas in 12 low-$z$ galaxy groups has been conducted using the Cosmic Origins Spectrograph (COS) on-board the Hubble Space Telescope (HST). Targets were selected from a large, homogeneously-selected sample of groups found in the Sloan Digital Sky Survey (SDSS). A critical selection criterion excluded sight lines that pass close ($<1.5$ virial radii) to a group galaxy, to ensure absorber association with the group as a whole. Deeper galaxy redshift observations are used both to search for closer galaxies and also to characterize these $10^{13.5}$ to $10^{14.5} M_{\odot}$ groups, the most massive of which are highly-virialized with numerous early-type galaxies (ETGs). This sample also includes two spiral-rich groups, not yet fully-virialized. At group-centric impact parameters of 0.3-2 Mpc, these $\mathrm{S/N}=15$-30 spectra detected HI absorption in 7 of 12 groups; high (OVI) and low (SiIII) ion metal lines are present in 2/3 of the absorption components. None of the three most highly-virialized, ETG-dominated groups are detected in absorption. Covering fractions $\gtrsim50$% are seen at all impact parameters probed, but do not require large filling factors despite an enormous extent. Unlike halo clouds in individual galaxies, group absorbers have radial velocities which are too low to escape the group potential well without doubt. This suggests that these groups are "closed boxes" for galactic evolution in the current epoch. Evidence is presented that the cool and warm group absorbers are not a pervasive intra-group medium (IGrM), requiring a hotter ($T\sim10^6$ to $10^7$ K) IGrM to be present to close the baryon accounting.


title: Holographic entanglement entropy from probe M-theory branes

abstract: We compute the holographic entanglement entropy contribution from planar two-dimensional defects in six-dimensional $\mathcal{N}=(2,0)$ superconformal field theory, holographically dual to probe M2- and M5-branes in $AdS_7 \times S^4$. In particular, we test the viability of the universal contribution of the defect to entanglement entropy as a candidate $C$-function. We find that this coefficient is not monotonic under defect renormalization group flows triggered by the vacuum expectation value of a marginal operator. Another candidate $C$-function, the on-shell action inside the entanglement wedge, monotonically decreases under the flows we study.


title: Tensor network analysis of critical coupling in two dimensional $\phi^{4}$ theory

abstract: We make a detailed analysis of the spontaneous $Z_{2}$-symmetry breaking in the two dimensional real $\phi^{4}$ theory with the tensor renormalization group approach, which allows us to take the thermodynamic limit easily and determine the physical observables without statistical uncertainties. We determine the critical coupling in the continuum limit employing the tensor network formulation for scalar field theories proposed in our previous paper. We obtain $\left[ \lambda / \mu_{\mathrm{c}}^{2} \right]_{\mathrm{cont.}} = 10.913(56)$ with the quartic coupling $\lambda$ and the renormalized critical mass $\mu_{\mathrm{c}}$. The result is compared with previous results obtained by different approaches.


title: Combining Refinement of Parametric Models with Goal-Oriented Reduction of Dynamics

abstract: Parametric models abstract part of the specification of dynamical models by integral parameters. They are for example used in computational systems biology, notably with parametric regulatory networks, which specify the global architecture (interactions) of the networks, while parameterising the precise rules for drawing the possible temporal evolutions of the states of the components. A key challenge is then to identify the discrete parameters corresponding to concrete models with desired dynamical properties. This paper addresses the restriction of the abstract execution of parametric regulatory (discrete) networks by the means of static analysis of reachability properties (goal states). Initially defined at the level of concrete parameterised models, the goal-oriented reduction of dynamics is lifted to parametric networks, and is proven to preserve all the minimal traces to the specified goal states. It results that one can jointly perform the refinement of parametric networks (restriction of domain of parameters) while reducing the necessary transitions to explore and preserving reachability properties of interest.


title: Solving the 3D High-Frequency Helmholtz Equation using Contour Integration and Polynomial Preconditioning

abstract: We propose an iterative solution method for the 3D high-frequency Helmholtz equation that exploits a contour integral formulation of spectral projectors. In this framework, the solution in certain invariant subspaces is approximated by solving complex-shifted linear systems, resulting in faster GMRES iterations due to the restricted spectrum. The shifted systems are solved by exploiting a polynomial fixed-point iteration, which is a robust scheme even if the magnitude of the shift is small. Numerical tests in 3D indicate that $O(n^{1/3})$ matrix-vector products are needed to solve a high-frequency problem with a matrix size $n$ with high accuracy. The method has a small storage requirement, can be applied to both dense and sparse linear systems, and is highly parallelizable.


title: Evolution of magneto-crystalline anisotropies in Mn$_{1-x}$Fe$_x$Si and Mn$_{1-x}$Co$_x$Si as inferred from small-angle neutron scattering and bulk properties

abstract: We report a comprehensive small-angle neutron scattering~(SANS) study of Mn$_{1-x}$Fe$_{x}$Si at zero magnetic field. To delineate changes of magneto-crystalline anisotropies (MCAs) from effects due to defects and disorder, we recorded complementary susceptibility and specific heat data, and investigated selected compositions of Mn$_{1-x}$Co$_{x}$Si. For all systems studied the transition temperature and magnetic phase diagrams evolve monotonically with composition consistent with literature. The SANS patterns of the magnetic order recorded under zero-field cooling display strong changes of the directions of the intensity maxima and smeared out intensity distributions as a function of composition. We show that cubic MCAs account for the complex evolution of the SANS patterns, where for increasing $x$ the character of the MCAs shifts from terms that are fourth-order to terms that are sixth order in spin--orbit coupling. The magnetic field dependence of the susceptibility and SANS establishes that the helix reorientation as a function of magnetic field for Fe- or Co-doped MnSi is dominated by pinning due to defects and disorder. The presence of thermodynamic anomalies of the specific heat at the phase boundaries of the skyrmion lattice phase in the doped samples and properties observed in Mn$_{1-x}$Co$_{x}$Si establishes that the pinning due to defects and disorder remains, however, weak and comparable to the field scale of the helix reorientation. The observation that MCAs, that are sixth order in spin-orbit coupling, play an important role for the spontaneous order in Mn$_{1-x}$Fe$_{x}$Si and Mn$_{1-x}$Co$_{x}$Si, offering a fresh perspective for a wide range of topics in cubic chiral magnets such as the generic magnetic phase diagram, the morphology of topological spin textures, the paramagnetic-to-helical transition, and quantum phase transitions.


title: Mixed field theory

abstract: We consider scalar field theory defined over a direct product of the real and $p$-adic numbers. An adjustable dynamical scaling exponent $z$ enters into the microscopic lagrangian, so that the Gaussian theories provide a line of fixed points. We argue that at $z=1/3$, a branch of Wilson-Fisher fixed points joins onto the line of Gaussian theories. We compute standard critical exponents at the Wilson-Fisher fixed points in the region where they are perturbatively accessible, including a loop correction to the dynamical critical exponent. We show that the classical propagator contains oscillatory behavior in the real direction, though the amplitude of these oscillations can be made exponentially small without fine-tuning parameters of the theory. Similar oscillatory behavior emerges in Fourier space from two-loop corrections, though again it can be highly suppressed. We also briefly consider compact $p$-adic extra dimensions, showing in non-linear, classical, scalar field theories that a form of consistent truncation allows us to retain only finitely many Kaluza-Klein modes in an effective theory formulated on the non-compact directions.


title: Expected performance of the TT-PET scanner

abstract: The TT-PET collaboration is developing a small animal TOF-PET scanner based on silicon detectors featuring 30 ps RMS time resolution and intended to be inserted in an existing MRI scanner. The TT-PET scanner makes use of a stack of layers of high-Z photon-converter and 100 $\mathrm{\mu m}$ thick silicon sensors, to achieve a scanner with 0.5 $\mathrm{\times}$ 0.5 $\mathrm{\times}$ 0.2 $\mathrm{mm^{3}}$ granularity, with precise depth-of-interaction measurement. In this paper we present the results of the Monte Carlo studies for the expected data rate, time resolution on the TOF measurements, spatial resolution and image reconstruction with and without the use of the timing information. Most of the studies have been performed according to the international standards used to assess the performance of small-animal PET system.


title: Network architecture of energy landscapes in mesoscopic quantum systems

abstract: Mesoscopic quantum systems exhibit complex many-body quantum phenomena, where interactions between spins and charges give rise to collective modes and topological states. Even simple, non-interacting theories display a rich landscape of energy states --- distinct many-particle configurations connected by spin- and energy-dependent transition rates. The collective energy landscape is difficult to characterize or predict, especially in regimes of frustration where many-body effects create a multiply degenerate landscape. Here we use network science to characterize the complex interconnection patterns of these energy-state transitions. Using an experimentally verified computational model of electronic transport through quantum antidots, we construct networks where nodes represent accessible energy states and edges represent allowed transitions. We then explore how physical changes in currents and voltages are reflected in the network topology. We find that the networks exhibit Rentian scaling, which is characteristic of efficient transportation systems in computer circuitry, neural circuitry, and human mobility, and can be used to measure the interconnection complexity of a network. Remarkably, networks corresponding to points of frustration in quantum transport (due, for example, to spin-blockade effects) exhibit an enhanced topological complexity relative to networks not experiencing frustration. Our results demonstrate that network characterizations of the abstract topological structure of energy landscapes can capture salient properties of quantum transport. More broadly, our approach motivates future efforts to use network science in understanding the dynamics and control of complex quantum systems.


title: A systematic path to non-Markovian dynamics: New response pdf evolution equations under Gaussian coloured noise excitation

abstract: Determining evolution equations governing the probability density function (pdf) of non-Markovian responses to random differential equations (RDEs) excited by coloured noise, is an important issue arising in various problems of stochastic dynamics, advanced statistical physics and uncertainty quantification of macroscopic systems. In the present work, such equations are derived for a scalar, nonlinear RDE under additive coloured Gaussian noise excitation, through the stochastic Liouville equation. The latter is an exact, yet non-closed equation, involving aver-ages over the time history of the non-Markovian response. This nonlocality is treated by applying an extension of the Novikov-Furutsu theorem and a novel approximation, employing a stochastic Volterra-Taylor functional expansion around instantaneous response moments, leading to efficient, closed, approximate equations for the response pdf. These equations retain a tractable amount of nonlocality and nonlinearity, and they are valid in both the transient and long-time regimes for any correlation function of the excitation. Also, they include as special cases various existing relevant models, and generalize Hanggi's ansatz in a rational way. Numerical results for a bistable nonlinear RDE confirm the accuracy and the efficiency of the new equations. Extension to the multidimensional case (systems of RDEs) is feasible, yet laborious.


title: Dwarf galaxies at the steep and faint end of the Tully Fisher relation

abstract: Nearby isolated galaxies ($z\sim 0$), are considered immerse within a thermal bath at 2.73 K. However the dwarf galaxies orbiting these galaxies are also subject to additional radiation from their hosts, so they are within a thermal bath slightly warmer. We claim that this thermal effect can explain several properties of the dwarf galaxies, such as their rotation curves, their radial acceleration relations (RAR) and why the velocity-mass dependence at the faint end of the Tully Fisher Relation (TFR) is steeper for these galaxies. In the Debye Gravitational Theory (DGT), the galaxies properties, such as the rotation curves, the RARs, and the TFRs are isothermal curves; they depend explicitly only of temperature, of the thermal bath in which are immersed. We show that the steep and faint end of the TFR is composed of a mixture of dwarf galaxies immersed in thermal baths with temperatures from 2.73 K to 3.80 K. A comparison among DGT's predictions for the dwarf galaxies relations with data obtained from the literature show a satisfactory agreement between them.


title: B-brane transport in anomalous (2,2) models and localization

abstract: We study how B-branes in two-dimensional N=(2,2) anomalous models behave as we vary the energy scale and bulk parameters in the quantum K\"ahler moduli space. We focus on (2,2) theories defined by abelian gauged linear sigma models (GLSM). Guided by the hemisphere partition function we find how B-branes split in arbitrary phases into components on the Higgs branch and other branches: this generalizes the band restriction rules of Herbst-Hori-Page to (abelian) anomalous models. Secondly, we address divergences in non-compact models, through the central example of GLSMs for Hirzebruch-Jung resolutions of cyclic surface singularities. For a brane with compact support we explain how to regularize and compute the hemisphere partition function and extract its Higgs branch component, which we match in the zero-instanton sector to the geometric central charge of the brane. To this aim, we clarify the definition of zero-instanton geometric central charge for objects in the derived category of a non-compact toric orbifold.


title: Tree-Structured Recurrent Switching Linear Dynamical Systems for Multi-Scale Modeling

abstract: Many real-world systems studied are governed by complex, nonlinear dynamics. By modeling these dynamics, we can gain insight into how these systems work, make predictions about how they will behave, and develop strategies for controlling them. While there are many methods for modeling nonlinear dynamical systems, existing techniques face a trade off between offering interpretable descriptions and making accurate predictions. Here, we develop a class of models that aims to achieve both simultaneously, smoothly interpolating between simple descriptions and more complex, yet also more accurate models. Our probabilistic model achieves this multi-scale property through a hierarchy of locally linear dynamics that jointly approximate global nonlinear dynamics. We call it the tree-structured recurrent switching linear dynamical system. To fit this model, we present a fully-Bayesian sampling procedure using Polya-Gamma data augmentation to allow for fast and conjugate Gibbs sampling. Through a variety of synthetic and real examples, we show how these models outperform existing methods in both interpretability and predictive capability.


title: 2D Hexagonal Cellular Automata: The Complexity of the Forms

abstract: We created two dimensional hexagonal cellular automata to obtain complexity. Considering the game of life rules, Wolfram's works about life-like structures and John von Neumann's self-replication, self-maintenance, self-reproduction problems, we developed 2-states and 3-states hexagonal growing algorithms that reach large populations through random initial states. Unlike the game of life, we used six neighbourhoods cellular automata instead of eight or four neighbourhoods. First simulations explained that whether we are able to obtain sort of oscillators, blinkers, and gliders. Inspired by Wolfram's 1D cellular automata complexity and life-like structures, we simulated 2D synchronous, discrete, deterministic cellular automata to reach life-like forms with 2-states cells. The life-like formations and the oscillators have been explained how they contribute to initiating self-maintenance together with self-reproduction and self-replication. After comparing the simulation results, we decided to develop the algorithm for another step. Appending a new state to the same algorithm, that we used for reaching life-like structures, led us to experiment new branching and fractal forms. All these studies tried to demonstrate that complex life forms might come from uncomplicated rules.


title: Partially composite supersymmetry

abstract: We consider a supersymmetric model that uses partial compositeness to explain the fermion mass hierarchy and predict the sfermion mass spectrum. The Higgs and third-generation matter superfields are elementary, while the first two matter generations are composite. Linear mixing between elementary superfields and supersymmetric operators with large anomalous dimensions is responsible for simultaneously generating the fermion and sfermion mass hierarchies. After supersymmetry is broken by the strong dynamics, partial compositeness causes the first- and second-generation sfermions to be split from the much lighter gauginos and third-generation sfermions. This occurs even though the tree-level soft masses of the elementary fields are subject to large radiative corrections from the composite sector, which we calculate in the gravitational dual theory using the AdS/CFT correspondence. The sfermion mass scale is constrained by the observed 125 GeV Higgs boson, leading to stop masses and gauginos around 10-100 TeV and the first two generation sfermion masses around 100-1000 TeV. This gives rise to a splitlike supersymmetric model that explains the fermion mass hierarchy while simultaneously predicting an inverted sfermion mass spectrum consistent with LHC and flavor constraints. Finally, the lightest supersymmetric particle is a gravitino in the keV to TeV range, which can play the role of dark matter.


title: Fast compressive Raman bio-imaging via matrix completion

abstract: Raman microscopy is a powerful method combining non-invasiveness with no special sample preparation. Because of this remarkable simplicity, it has been widely exploited in many fields, ranging from life and materials sciences, to engineering. Notoriously, due to the required imaging speeds for bio-imaging, it has remained a challenge how to use this technique for dynamic and large-scale imaging. Recently, compressive Raman has been put forward, allowing for fast imaging, therefore solving the issue of speed. Yet, due to the need of strong a priori information of the species forming the hyperspectrum, it has remained elusive how to apply this technique for microspectroscopy of (dynamic) biological tissues. Combining an original spectral under-sampling measurement technique with matrix completion framework for reconstruction, we demonstrate fast and inexpensive label-free molecular imaging of biological specimens (brain tissues and single cells). Therefore, our results open interesting perspectives for clinical and cell biology applications using the much faster compressive Raman framework.


title: Mass Hierarchy and Vacuum Energy

abstract: A hierarchically small weak scale does not generally coincide with enhanced symmetry, but it may still be exceptional with respect to vacuum energy. By analyzing the classical vacuum energy as a function of parameters such as the Higgs mass, we show how near-criticality, i.e. fine-tuning, corresponds universally to boundaries where the vacuum energy transitions from exactly flat to concave down. In the presence of quantum corrections, these boundary regions can easily be perturbed to become maxima of the vacuum energy. After introducing a dynamical scalar field $\phi$ which scans the Higgs sector parameters, we propose several possible mechanisms by which this field could be localized to the maximum. One possibility is that the $\phi$ potential has many vacua, with those near the maximum vacuum energy expanding faster during a long period of cosmic inflation and hence dominating the volume of the Universe. Alternately, we describe scenarios in which vacua near the maximum could be anthropically favored, due to selection of the late-time cosmological constant or dark matter density. Independent of these specific approaches, the physical value of the weak scale in our proposal is generated naturally and dynamically from loops of heavy states coupled to the Higgs. These states are predicted to be a loop factor heavier than in models without this mechanism, avoiding tension with experimental null results.


title: Photon Fluence and Dose Estimation in Computed Tomography using a Discrete Ordinates Boltzmann Solver

abstract: In this study, cone-beam single projection and axial CT scans are modeled with a software package - DOCTORS, which solves the linear Boltzmann equation using the discrete ordinates method. Phantoms include a uniform 35 cm diameter water cylinder and a non-uniform abdomen phantom. Series simulations were performed with different simulation parameters, including the number of quadrature angles, the order of Legendre polynomial expansions, and coarse and fine mesh grid. Monte Carlo simulations were also performed to benchmark DOCTORS simulations. A quantitative comparison was made between the simulation results obtained using DOCTORS and Monte Carlo methods. The deterministic simulation was in good agreement with the Monte Carlo simulation on dose estimation, with a root-mean-square-deviation (RMSD) difference of around 2.87 percent. It was found that the contribution of uncollided photon fluence directly from the source dominates the local absorbed dose in the diagnostic X-ray energy range. The uncollided photon fluence can be calculated accurately using a ray-tracing algorithm. The accuracy of collided photon fluence estimation is largely affected by the pre-calculated multigroup cross-sections. The primary benefit of DOCTORS lies in its rapid computation speed. Using DOCTORS, parallel computing using GPU enables the cone-beam CT dose estimation nearly in real-time.


title: When Left and Right Turns Inside Out: A Geometric and Categorical Introduction to an Inverse Problem in Persistence

abstract: In this paper we introduce the problem of counting embedded spheres in R^3 whose projection to the z-axis yields a level set barcode of a particular type. Two embedded spheres are considered height equivalent if they are related by a z-level set preserving isotopy. A formula previously used to count functions on the interval with the same sub-level set persistence barcode provides a lower bound on height equivalence classes with the same level set persistence. A conjectured upper bound is provided as well.


title: Continuous-variable quantum repeater based on quantum scissors and mode multiplexing

abstract: Quantum repeaters are indispensable for high-rate, long-distance quantum communications. The vision of a future quantum internet strongly hinges on realizing quantum repeaters in practice. Numerous repeaters have been proposed for discrete-variable (DV) single-photon-based quantum communications. Continuous variable (CV) encodings over the quadrature degrees of freedom of the electromagnetic field mode offer an attractive alternative. For example, CV transmission systems are easier to integrate with existing optical telecom systems compared to their DV counterparts. Yet, repeaters for CV have remained elusive. We present a novel quantum repeater scheme for CV entanglement distribution over a lossy bosonic channel that beats the direct transmission exponential rate-loss tradeoff. The scheme involves repeater nodes consisting of a) two-mode squeezed vacuum (TMSV) CV entanglement sources, b) the quantum scissors operation to perform nondeterministic noiseless linear amplification of lossy TMSV states, c) a layer of switched, mode multiplexing inspired by second-generation DV repeaters, which is the key ingredient apart from probabilistic entanglement purification that makes DV repeaters work, and d) a non-Gaussian entanglement swap operation. We report our exact results on the rate-loss envelope achieved by the scheme.


title: Models of hypercolor based on symplectic gauge group with three heavy vectorlike hyperquarks

abstract: We study possibilities to extend the Standard Model (SM) by three flavors of vectorlike heavy quarks in pseudoreal representation of symplectic hypercolor gauge group. This extension of SM predicts a rich spectra of heavy composite hypermesons and hyperbaryons (all of them carry integer spins) including fourteen pseudo-Nambu--Goldstone states emerging in dynamical breaking of the global symmetry group of the H-quarks, SU(6), to its Sp(6) subgroup. The properties of the lightest states depend strongly on the choice of heavy-quark hypercharges. Our focus is placed on the variants of the model with partially composite Higgs boson, i.e. the experimentally observed boson is comprised of the elementary SM Higgs and a mixture of H-hadrons.


title: CNN-Cert: An Efficient Framework for Certifying Robustness of Convolutional Neural Networks

abstract: Verifying robustness of neural network classifiers has attracted great interests and attention due to the success of deep neural networks and their unexpected vulnerability to adversarial perturbations. Although finding minimum adversarial distortion of neural networks (with ReLU activations) has been shown to be an NP-complete problem, obtaining a non-trivial lower bound of minimum distortion as a provable robustness guarantee is possible. However, most previous works only focused on simple fully-connected layers (multilayer perceptrons) and were limited to ReLU activations. This motivates us to propose a general and efficient framework, CNN-Cert, that is capable of certifying robustness on general convolutional neural networks. Our framework is general -- we can handle various architectures including convolutional layers, max-pooling layers, batch normalization layer, residual blocks, as well as general activation functions; our approach is efficient -- by exploiting the special structure of convolutional layers, we achieve up to 17 and 11 times of speed-up compared to the state-of-the-art certification algorithms (e.g. Fast-Lin, CROWN) and 366 times of speed-up compared to the dual-LP approach while our algorithm obtains similar or even better verification bounds. In addition, CNN-Cert generalizes state-of-the-art algorithms e.g. Fast-Lin and CROWN. We demonstrate by extensive experiments that our method outperforms state-of-the-art lower-bound-based certification algorithms in terms of both bound quality and speed.


title: A Physical Picture of Bispectrum Baryon Acoustic Oscillations in the Interferometric Basis

abstract: We present a picture of the matter bispectrum in a novel "interferometric" basis designed to highlight interference of the baryon acoustic oscillations (BAO) in the power spectra composing it. Triangles where constructive interference amplifies BAO provide stronger cosmic distance constraints than triangles with destructive interference. We show that the amplitude of the BAO feature in the full cyclically summed bispectrum can be decomposed into simpler contributions from single terms or pairs of terms in the perturbation theory bispectrum, and that across large swathes of our parameter space the full BAO amplitude is described well by the amplitude of BAO in a single term. The dominant term is determined largely by the $F^{(2)}$ kernel of Eulerian standard perturbation theory. We present a simple physical picture of the BAO amplitude in each term; the BAO signal is strongest in triangle configurations where two wavenumbers differ by a multiple of the BAO fundamental wavelength.


title: Komondor: a Wireless Network Simulator for Next-Generation High-Density WLANs

abstract: Komondor is a wireless network simulator for next-generation wireless local area networks (WLANs). The simulator has been conceived as an accessible (ready-to-use) open source tool for research on wireless networks and academia. An important advantage of Komondor over other well-known wireless simulators lies in its high event processing rate, which is furnished by the simplification of the core operation. This allows outperforming the execution time of other simulators like ns-3, thus supporting large-scale scenarios with a huge number of nodes. In this paper, we provide insights into the Komondor simulator and overview its main features, development stages and use cases. The operation of Komondor is validated in a variety of scenarios against different tools: the ns-3 simulator and two analytical tools based on Continuous Time Markov Networks (CTMNs) and the Bianchi's DCF model. Results show that Komondor captures the IEEE 802.11 operation very similarly to ns-3. Finally, we discuss the potential of Komondor for simulating complex environments -- even with machine learning support -- in next-generation WLANs by easily developing new user-defined modules of code


title: Cosmological Filaments in the light of Excursion Set of Saddle Points

abstract: The universe in large scales is structured as a network known as cosmic web. Filaments are one of the structural components of this web, which can be introduced as a novel probe to study the formation and evolution of structures and as a probe to study the cosmological models and to address the missing baryon problem. The aim of this work is to introduce an analytical framework to study the statistics of filaments such as number density of them and also to obtain the length-mass relation. For this objective, we model filaments as collapsed objects which have an extension in one direction, accordingly we use the ellipsoidal collapse to study the evolution of an over-dense region via gravitational instability. %In this context, we consider the critical density of filament formation, which will be a crucial parameter for number density count. We find that the nonlinear density of filaments in the epoch of formation is almost mass independent and is in order of $\sim 30$. By introducing {\it{filament's extended condition}} we find a fitting function for length-mass relation. For the statistics of filaments, we propose a novel framework named excursion set of saddle points. In this approach, we count the saddle points of the density field Hessian matrix, and relate it to the count of filaments. In addition, we addressed the filament in filament problem with up-crossing approximation.


title: An extremal problem of regular simplices-The higher dimensional case

abstract: The new result of this paper connected with the following problem: Consider a supporting hyperplane of a regular simplex and its re ected image at this hyperplane. When will be the volume of the convex hull of these two simplices maximal? We prove that in the case when the dimension is less or equal to 4, the maximal volume achieves in that case when the hyperplane goes through on a vertex and orthogonal to the height of the simplex at this vertex. More interesting that in the higher dimensional cases this position is not optimal. We also determine the optimal position of hyperplane in the 5-dimensional case. This corrects an erroneous statement in my paper [3].


title: Phases of 5d SCFTs from M-/F-theory on Non-Flat Fibrations

abstract: We initiate the systematic investigation of non-flat resolutions of non-minimal singularities in elliptically fibered Calabi--Yau threefolds. Compactification of M-theory on these geometries provides an alternative approach to studying phases of five-dimensional superconformal field theories (5d SCFTs). We argue that such resolutions capture non-trivial holonomies in the circle reduction of the 6d conformal matter theory that is the F-theory interpretation of the singular fibration. As these holonomies become mass deformations in the 5d theory, non-flat resolutions furnish a novel method in the attempt to classify 5d SCFTs through 6d SCFTs on a circle. A particularly pleasant aspect of this proposal is the explicit embedding of the 5d SCFT's enhanced flavor group inside that of the parent 6d SCFT, which can be read off from the geometry. We demonstrate these features in toric examples which realize 5d theories up to rank four.


title: A public framework for Feynman calculations and post-Newtonian gravity

abstract: We report here on a line of work that has played a key role in formally establishing and going beyond the state of the art in the effective field theory (EFT) approach and in post-Newtonian (PN) gravity. We also outline here how this comprehensive framework in fact forms the outset of a prospective rich research program, building on the public Feynman and PN technology developed.


title: On the Implicit Assumptions of GANs

abstract: Generative adversarial nets (GANs) have generated a lot of excitement. Despite their popularity, they exhibit a number of well-documented issues in practice, which apparently contradict theoretical guarantees. A number of enlightening papers have pointed out that these issues arise from unjustified assumptions that are commonly made, but the message seems to have been lost amid the optimism of recent years. We believe the identified problems deserve more attention, and highlight the implications on both the properties of GANs and the trajectory of research on probabilistic models. We recently proposed an alternative method that sidesteps these problems.


title: New Convergence Aspects of Stochastic Gradient Algorithms

abstract: The classical convergence analysis of SGD is carried out under the assumption that the norm of the stochastic gradient is uniformly bounded. While this might hold for some loss functions, it is violated for cases where the objective function is strongly convex. In Bottou et al. (2018), a new analysis of convergence of SGD is performed under the assumption that stochastic gradients are bounded with respect to the true gradient norm. We show that for stochastic problems arising in machine learning such bound always holds; and we also propose an alternative convergence analysis of SGD with diminishing learning rate regime. We then move on to the asynchronous parallel setting, and prove convergence of Hogwild! algorithm in the same regime in the case of diminished learning rate. It is well-known that SGD converges if a sequence of learning rates $\{\eta_t\}$ satisfies $\sum_{t=0}^\infty \eta_t \rightarrow \infty$ and $\sum_{t=0}^\infty \eta^2_t < \infty$. We show the convergence of SGD for strongly convex objective function without using bounded gradient assumption when $\{\eta_t\}$ is a diminishing sequence and $\sum_{t=0}^\infty \eta_t \rightarrow \infty$. In other words, we extend the current state-of-the-art class of learning rates satisfying the convergence of SGD.


title: Non-thermal electrons from solar nanoflares

abstract: Context. We introduce a model for including accelerated particles in pure magnetohydrodynamics (MHD) simulations of the solar atmosphere. Aims. We show that the method is viable and produces results that enhance the realism of MHD simulations of the solar atmosphere. Methods. The acceleration of high-energy electrons in solar flares is an accepted fact, but is not included in the most advanced 3D simulations of the solar atmosphere. The effect of the acceleration is not known, and here we introduce a simple method to account for the ability of the accelerated electrons to move energy from the reconnection sites and into the dense transition zone and chromosphere. Results. The method was only run for a short time and with low reconnection energies, but this showed that the reconnection process itself changes, and that there is a clear effect on the observables at the impact sites of the accelerated electrons. Further work will investigate the effect on the reconnection sites and the impact sites in detail.


title: CHILES: HI morphology and galaxy environment at z=0.12 and z=0.17

abstract: We present a study of 16 HI-detected galaxies found in 178 hours of observations from Epoch 1 of the COSMOS HI Large Extragalactic Survey (CHILES). We focus on two redshift ranges between 0.108 <= z <= 0.127 and 0.162 <= z <= 0.183 which are among the worst affected by radio frequency interference (RFI). While this represents only 10% of the total frequency coverage and 18% of the total expected time on source compared to what will be the full CHILES survey, we demonstrate that our data reduction pipeline recovers high quality data even in regions severely impacted by RFI. We report on our in-depth testing of an automated spectral line source finder to produce HI total intensity maps which we present side-by-side with significance maps to evaluate the reliability of the morphology recovered by the source finder. We recommend that this become a common place manner of presenting data from upcoming HI surveys of resolved objects. We use the COSMOS 20k group catalogue, and we extract filamentary structure using the topological DisPerSE algorithm to evaluate the \hi\ morphology in the context of both local and large-scale environments and we discuss the shortcomings of both methods. Many of the detections show disturbed HI morphologies suggesting they have undergone a recent interaction which is not evident from deep optical imaging alone. Overall, the sample showcases the broad range of ways in which galaxies interact with their environment. This is a first look at the population of galaxies and their local and large-scale environments observed in HI by CHILES at redshifts beyond the z=0.1 Universe.


title: Coherent control of electron-nuclear states of rare-earth ions in crystals using radio-frequency and microwave radiation

abstract: We have demonstrated electron-electron and electron-nuclear spin manipulations of Gd3+ ion in CaWO4 crystal. The results suggest that the studied system is perspective for multiqubit implementation in quantum computing.

categories: cond-mat.mes-hall cond-mat.mtrl-sci cond-mat.str-el quant-ph


title: On the properties of spectral effect algebras

abstract: The aim of this paper is to show that there can be either only one or uncountably many contexts in any spectral effect algebra, answering a question posed in [S. Gudder, Convex and Sequential Effect Algebras, (2018), arXiv:1802.01265]. We also provide some results on the structure of spectral effect algebras and their state spaces and investigate the direct products and direct convex sums of spectral effect algebras. In the case of spectral effect algebras with sharply determining state space, stronger properties can be proved: the spectral decompositions are essentially unique, the algebra is sharply dominating and the set of its sharp elements is an orthomodular lattice. The article also contains a list of open questions that might provide interesting future research directions.


title: From Context to Concept: Exploring Semantic Relationships in Music with Word2Vec

abstract: We explore the potential of a popular distributional semantics vector space model, word2vec, for capturing meaningful relationships in ecological (complex polyphonic) music. More precisely, the skip-gram version of word2vec is used to model slices of music from a large corpus spanning eight musical genres. In this newly learned vector space, a metric based on cosine distance is able to distinguish between functional chord relationships, as well as harmonic associations in the music. Evidence, based on cosine distance between chord-pair vectors, suggests that an implicit circle-of-fifths exists in the vector space. In addition, a comparison between pieces in different keys reveals that key relationships are represented in word2vec space. These results suggest that the newly learned embedded vector representation does in fact capture tonal and harmonic characteristics of music, without receiving explicit information about the musical content of the constituent slices. In order to investigate whether proximity in the discovered space of embeddings is indicative of `semantically-related' slices, we explore a music generation task, by automatically replacing existing slices from a given piece of music with new slices. We propose an algorithm to find substitute slices based on spatial proximity and the pitch class distribution inferred in the chosen subspace. The results indicate that the size of the subspace used has a significant effect on whether slices belonging to the same key are selected. In sum, the proposed word2vec model is able to learn music-vector embeddings that capture meaningful tonal and harmonic relationships in music, thereby providing a useful tool for exploring musical properties and comparisons across pieces, as a potential input representation for deep learning models, and as a music generation device.


title: Operational nonlocality

abstract: An operational concept of locality whose quantum violation is indicated independently of any other assumption(s) seems to be lacking in the quantum foundations literature so far. Bell's theorem only shows that quantum correlations violate the conjunction of the ontological assumptions of localism and determinism. Taking a cue from computational complexity theory, here we define such a concept of locality in terms of a class of operational decision problems, and propose that Einstein-Podolsky-Rosen (EPR) steering in its basic asymmetric formulation gives a specific realization of this class. Various nonclassical (convex) operational theories, including quantum mechanics, are shown to be nonlocal according to this operational criterion. We discuss several ramifications arising from this result. It indicates a basic difference between quantum no-signaling and relativistic no-signaling, and suggests an information theoretic derivation of the former from other basic principles in the convex framework. It elucidates the connection between two different bounds on nonlocality, due to Oppenheim and Wehner (2010) and Banik et al. (2015), thereby highlighting the interplay and contrast between uncertainty and measurement incompatibility. Our result provides an argument supporting the realist interpretation of the quantum state, and helps clarify why Spekkens' toy model (2007) features steering but not Bell nonlocality.


title: The Impact of Enhanced Halo Resolution on the Simulated Circumgalactic Medium

abstract: Traditional cosmological hydrodynamics simulations fail to spatially resolve the circumgalatic medium (CGM), the reservoir of tenuous gas surrounding a galaxy and extending to its virial radius. We introduce the technique of Enhanced Halo Resolution (EHR), enabling more realistic physical modeling of the simulated CGM by consistently forcing gas refinement to smaller scales throughout the virial halo of a simulated galaxy. We investigate the effects of EHR in the Tempest simulations, a suite of Enzo-based cosmological zoom simulations following the evolution of an L* galaxy, resolving spatial scales of 500 comoving pc out to 100 comoving kpc in galactocentric radius. Among its many effects, EHR (1) changes the thermal balance of the CGM, increasing its cool gas content and decreasing its warm/hot gas content; (2) preserves cool gas structures for longer periods; and (3) enables these cool clouds to exist at progressively smaller size scales. Observationally, this results in a boost in "low ions" like H I and a drop in "high ions" like O VI throughout the CGM. These effects of EHR do not converge in the Tempest simulations, but extrapolating these trends suggests that the CGM in reality is a mist consisting of ubiquitous, small, long-lived, cool clouds suspended in a hot medium at the virial temperature of the halo. Additionally, we explore the physical mechanisms to explain why EHR produces the above effects, proposing that it works both by (1) better sampling the distribution of CGM phases enabling runaway cooling in the denser, cooler tail of the phase distribution; and (2) preventing cool gas clouds from artificially mixing with the ambient hot halo and evaporating. Evidence is found for both EHR mechanisms occurring in the Tempest simulations.


title: High-resolution Charge Exchange Spectra with L-shell Nickel Show Striking Differences from Models

abstract: We present the first high-resolution laboratory spectra of X-ray emission following L-shell charge exchange between nickel ions and neutral H2 and He. We employ the commonly used charge exchange models found in XSPEC and SPEX, ACX and SPEX-CX, to simulate our experimental results. We show that significant differences between data and models exist in both line energies and strengths. In particular, we find that configuration mixing may play an important role in generating lines from core-excited states, and may be improperly treated in models. Our results indicate that if applied to astrophysical data, these models may lead to incorrect assumptions of the physical and chemical parameters of the region of interest.


title: First detection of scale-dependent linear halo bias in $N$-body simulations with massive neutrinos

abstract: Using $N$-body simulations with massive neutrino density perturbations, we detect the scale-dependent linear halo bias with high significance. This is the first time that this effect is detected in simulations containing neutrino density perturbations on all scales, confirming the same finding from separate universe simulations. The scale dependence is the result of the additional scale in the system, i.e. the massive neutrino free-streaming length, and it persists even if the bias is defined with respect to the cold dark matter plus baryon (instead of total matter) power spectrum. The separate universe approach provides a good model for the scale-dependent linear bias, and the effect is approximately $0.25f_\nu$ and $0.43f_\nu$ for halos with bias of 1.7 and 3.5, respectively. While the size of the effect is small, it is ${\it not}$ insignificant in terms of $f_\nu$ and should therefore be included to accurately constrain neutrino mass from clustering statistics of biased tracers. More importantly, this feature is a distinct signature of free-streaming particles and cannot be mimicked by other components of the standard cosmological model.


title: Phat ELVIS: The inevitable effect of the Milky Way's disk on its dark matter subhaloes

abstract: We introduce an extension of the ELVIS project to account for the effects of the Milky Way galaxy on its subhalo population. Our simulation suite, Phat ELVIS, consists of twelve high-resolution cosmological dark matter-only (DMO) zoom simulations of Milky Way-size $\Lambda$CDM~ haloes ($M_{\rm v} = 0.7-2 \times 10^{12} \,\mathrm{M}_\odot$) along with twelve re-runs with embedded galaxy potentials grown to match the observed Milky Way disk and bulge today. The central galaxy potential destroys subhalos on orbits with small pericenters in every halo, regardless of the ratio of galaxy mass to halo mass. This has several important implications. 1) Most of the $\mathtt{Disk}$ runs have no subhaloes larger than $V_{\rm max} = 4.5$ km s$^{-1}$ within $20$ kpc and a significant lack of substructure going back $\sim 8$ Gyr, suggesting that local stream-heating signals from dark substructure will be rare. 2) The pericenter distributions of Milky Way satellites derived from $\mathit{Gaia}$ data are remarkably similar to the pericenter distributions of subhaloes in the $\mathtt{Disk}$ runs, while the DMO runs drastically over-predict galaxies with pericenters smaller than 20 kpc. 3) The enhanced destruction produces a tension opposite to that of the classic `missing satellites' problem: in order to account for ultra-faint galaxies known within $30$ kpc of the Galaxy, we must populate haloes with $V_\mathrm{peak} \simeq 7$ km s$^{-1}$ ($M \simeq 3 \times 10^{7} \,\mathrm{M}_\odot$ at infall), well below the atomic cooling limit of $V_\mathrm{peak} \simeq 16$ km s$^{-1}$ ($M \simeq 5 \times 10^{8} \,\mathrm{M}_\odot$ at infall). 4) If such tiny haloes do host ultra-faint dwarfs, this implies the existence of $\sim 1000$ satellite galaxies within 300 kpc of the Milky Way.


title: The Little Dippers: Transits of Star-grazing Exocomets?

abstract: We describe EPIC 205718330 and EPIC 235240266, two systems identified in the K2 data whose light curves contain episodic drops in brightness with shapes and durations similar to those of the young "dipper" stars, yet shallower by ~1-2 orders of magnitude. These "little dippers" have diverse profile shapes with durations of ~0.5-1.0 days and depths of ~0.1-1.0% in flux; however, unlike most of the young dipper stars, these do not exhibit any detectable infrared excess indicative of protoplanetary disks, and our ground-based follow-up spectra lack any signatures of youth while indicating these objects as kinematically old. After ruling out instrumental and/or data processing artifacts as sources of the dimming events, we investigate possible astrophysical mechanisms based on the light curve and stellar properties. We argue that the little dippers are consistent with transits of star-grazing exocomets, and speculate that they are signposts of massive non-transiting exoplanets driving the close-approach orbits.


title: Interacting spin-3/2 fermions in a Luttinger semimetal: Competing phases and their selection in the global phase diagram

abstract: We compute the effects of electronic interactions on gapless spin-3/2 excitations that in a noninteracting system emerge at a bi-quadratic touching of Kramers degenerate valence and conduction bands, known as Luttinger semimetal. This model can describe the low-energy physics of HgTe, gray-Sn, 227 pyrochlore iridates and half-Heuslers. For the sake of concreteness we only consider the short-range components of the Coulomb interaction. By combining mean-field analysis with a renormalization group (RG) calculation, we construct multiple cuts of the global phase diagram of interacting spin-3/2 fermions at zero and finite temperature and chemical doping. Such phase diagrams display a rich confluence of competing orders, among which rotational symmetry breaking nematic insulators and time reversal symmetry breaking magnetic orders are the prominent excitonic phases. We also show that even repulsive interactions can be conducive for both s-wave and d-wave pairings. The reconstructed band structure inside the ordered phases allows us to organize them according to the energy (entropy) gain in the following (reverse) order: s-wave pairing, nematic phases, magnetic orders and d-wave pairings, at zero chemical doping. But, the paired states are always energetically superior over the excitonic ones for finite doping. The phase diagrams obtained from the RG analysis show that an ordered phase with higher energy (entropy) gain is realized at low (high) temperature. In addition, we establish a "selection rule" between the interaction channels and the resulting ordered phases, suggesting that repulsive interactions in the magnetic (nematic) channels are conducive for the nucleation of d-wave (s-wave) pairing. The proposed methodology can shed light on the global phase diagram of other strongly interacting multi-band systems, such as doped Dirac semimetal, topological insulators and the like.


title: Evidence of Particle Acceleration in the Superbubble 30 Doradus C with NuSTAR

abstract: We present evidence of diffuse, non-thermal X-ray emission from the superbubble 30 Doradus C (30 Dor C) using hard X-ray images and spectra from NuSTAR observations. For this analysis, we utilize data from a 200 ks targeted observation of 30 Dor C as well as 2.8 Ms of serendipitous off-axis observations from the monitoring of nearby SN 1987A. The complete shell of 30 Dor C is detected up to 20 keV, and the young supernova remnant MCSNR J0536-6913 in the southeast of 30 Dor C is not detected above 8 keV. Additionally, six point sources identified in previous Chandra and XMM-Newton investigations have hard X-ray emission coincident with their locations. Joint spectral fits to the NuSTAR and XMM-Newton spectra across the 30 Dor C shell confirm the non-thermal nature of the diffuse emission. Given the best-fit rolloff frequencies of the X-ray spectra, we find maximum electron energies of 70-110 TeV (assuming a B-field strength of 4$\mu$G), suggesting 30 Dor C is accelerating particles. Particles are either accelerated via diffusive shock acceleration at locations where the shocks have not stalled behind the H$\alpha$ shell, or cosmic-rays are accelerated through repeated acceleration of low-energy particles via turbulence and magnetohydrodynamic waves in the bubble's interior.


title: Star formation and gas in the minor merger UGC 10214

abstract: UGC 10214 is a minor merger in which a dwarf galaxy has interacted with a large spiral galaxy $\sim$250 Myr ago and produced a perturbed disk and a giant tidal tail. We use a multiwavelength dataset in order to study the present and past star formation rate (SFR) and its relation to the gas and stellar mass at a spatial resolution down to 4 kpc. UGC 10214 is a very massive (stellar mass $M_{{\rm \star}}$ = $1.28\times 10^{11}$$M_\odot$) galaxy with a low gas fraction ($M_{{\rm gas}}$/$M_{{\rm \star}}$ = 0.24), a high molecular gas fraction ($M_{{\rm H}_2}$/$M_{\rm HI}$ = 0.4) and a modest SFR (2-5 $M_\odot$ yr$^{-1}$). The comparison of the molecular gas mass and current SFR gives a molecular gas depletion time of about $\sim$ 2 Gyr (based on H$\alpha$), comparable to those of normal spiral galaxies. Both from a comparison of the H$\alpha$ emission, tracing the current SFR, and far-ultraviolet (FUV) emission, tracing the recent SFR during the past tens of Myr, as well as from spectral energy distribution (SED) fitting with CIGALE, we find that the SFR has increased by a factor of about 2-3 during the recent past. This increase is particularly noticeable in the centre of the galaxy. A pixel-to-pixel comparison of the SFR, molecular gas mass and stellar mass shows that the central region has had a depressed FUV-traced SFR, both compared to the molecular gas and the stellar mass, whereas the H$\alpha$-traced SFR shows a normal level. The atomic and molecular gas distribution is asymmetric, but the position-velocity diagram along the major axis shows a pattern of regular rotation. We conclude that the minor merger has most likely caused variations in the SFR in the past resulting in a moderate increase of the SFR, but it has not perturbed the gas significantly so that the molecular depletion time remains normal.


title: Efficient simulation of finite-temperature open quantum systems

abstract: Chain-mapping techniques in combination with the time-dependent density matrix renormalization group are a powerful tool for the simulation of open-system quantum dynamics. For finite-temperature environments, however, this approach suffers from an unfavorable algorithmic scaling with increasing temperature. We prove that the system dynamics under thermal environments can be non-perturbatively described by temperature-dependent system-environmental couplings with the initial environment state being in its pure vacuum state, instead of a mixed thermal state. As a consequence, as long as the initial system state is pure, the global system-environment state remains pure at all times. The resulting speedup and relaxed memory requirements of this approach enable the efficient simulation of open quantum systems interacting with highly structured environments in any temperature range, with applications extending from quantum thermodynamics to quantum effects in mesoscopic systems.


title: Gravity with antisymmetric components

abstract: This work proposes a new gravitational theory formulated in terms of the vierbein field. The vierbein contains components which can be shifted by local Lorentz transformations and therefore do not show up in the spacetime metric. These components are given dynamics and become physical in our setup. They enter the massless theory in the form of an antisymmetric tensor field which makes the action reminiscent of the bosonic sector of supergravity. We then demonstrate that both the metric and the antisymmetric tensor can be made massive by adding a potential term for the vierbein. The form of this mass potential is inspired by ghost-free massive gravity. We confirm the absence of additional and potentially pathological degrees of freedom in an ADM analysis. However, at the linearized level around maximally symmetric solutions, the fluctuation of the antisymmetric tensor has a tachyonic mass pole.


title: Using a Recurrent Neural Network to Reconstruct Quantum Dynamics of a Superconducting Qubit from Physical Observations

abstract: At its core, Quantum Mechanics is a theory developed to describe fundamental observations in the spectroscopy of solids and gases. Despite these practical roots, however, quantum theory is infamous for being highly counterintuitive, largely due to its intrinsically probabilistic nature. Neural networks have recently emerged as a powerful tool that can extract non-trivial correlations in vast datasets. They routinely outperform state-of-the-art techniques in language translation, medical diagnosis and image recognition. It remains to be seen if neural networks can be trained to predict stochastic quantum evolution without a priori specifying the rules of quantum theory. Here, we demonstrate that a recurrent neural network can be trained in real time to infer the individual quantum trajectories associated with the evolution of a superconducting qubit under unitary evolution, decoherence and continuous measurement from raw observations only. The network extracts the system Hamiltonian, measurement operators and physical parameters. It is also able to perform tomography of an unknown initial state without any prior calibration. This method has potential to greatly simplify and enhance tasks in quantum systems such as noise characterization, parameter estimation, feedback and optimization of quantum control.


title: Magneto-immutable turbulence in weakly collisional plasmas

abstract: We propose that pressure anisotropy causes weakly collisional turbulent plasmas to self-organize so as to resist changes in magnetic-field strength. We term this effect "magneto-immutability" by analogy with incompressibility (resistance to changes in pressure). The effect is important when the pressure anisotropy becomes comparable to the magnetic pressure, suggesting that in collisionless, weakly magnetized (high-$\beta$) plasmas its dynamical relevance is similar to that of incompressibility. Simulations of magnetized turbulence using the weakly collisional Braginskii model show that magneto-immutable turbulence is surprisingly similar, in most statistical measures, to critically balanced MHD turbulence. However, in order to minimize magnetic-field variation, the flow direction becomes more constrained than in MHD, and the turbulence is more strongly dominated by magnetic energy (a nonzero "residual energy"). These effects represent key differences between pressure-anisotropic and fluid turbulence, and should be observable in the $\beta\gtrsim1$ turbulent solar wind.


title: Candidate Massive Galaxies at $z \sim 4$ in the Dark Energy Survey

abstract: Using stellar population models, we predicted that the Dark Energy Survey (DES) - due to its special combination of area (5000 deg. sq.) and depth ($i = 24.3$) - would be in the position to detect massive ($\gtrsim 10^{11}$ M$_{\odot}$) galaxies at $z \sim 4$. We confront those theoretical calculations with the first $\sim 150$ deg. sq. of DES data reaching nominal depth. From a catalogue containing $\sim 5$ million sources, $\sim26000$ were found to have observed-frame $g-r$ vs $r-i$ colours within the locus predicted for $z \sim 4$ massive galaxies. We further removed contamination by stars and artefacts, obtaining 606 galaxies lining up by the model selection box. We obtained their photometric redshifts and physical properties by fitting model templates spanning a wide range of star formation histories, reddening and redshift. Key to constrain the models is the addition, to the optical DES bands $g$, $r$, $i$, $z$, and $Y$, of near-IR $J$, $H$, $K_{s}$ data from the Vista Hemisphere Survey. We further applied several quality cuts to the fitting results, including goodness of fit and a unimodal redshift probability distribution. We finally select 233 candidates whose photometric redshift probability distribution function peaks around $z\sim4$, have high stellar masses ($\log($M$^{*}$/M$_{\odot})\sim 11.7$ for a Salpeter IMF) and ages around 0.1 Gyr, i.e. formation redshift around 5. These properties match those of the progenitors of the most massive galaxies in the local universe. This is an ideal sample for spectroscopic follow-up to select the fraction of galaxies which is truly at high redshift. These initial results and those at the survey completion, which we shall push to higher redshifts, will set unprecedented constraints on galaxy formation, evolution, and the re-ionisation epoch.


title: Variational optimization in the AI era: Computational Graph States and Supervised Wave-function Optimization

abstract: Representing a target quantum state by a compact, efficient variational wave-function is an important approach to the quantum many-body problem. In this approach, the main challenges include the design of a suitable variational ansatz and optimization of its parameters. In this work, we address both of these challenges. First, we define the variational class of Computational Graph States (CGS) which gives a uniform framework for describing all computable variational ansatz. Secondly, we develop a novel optimization scheme, supervised wave-function optimization (SWO), which systematically improves the optimized wave-function by drawing on ideas from supervised learning. While SWO can be used independently of CGS, utilizing them together provides a flexible framework for the rapid design, prototyping and optimization of variational wave-functions. We demonstrate CGS and SWO by optimizing for the ground state wave-function of 1D and 2D Heisenberg models on nine different variational architectures including architectures not previously used to represent quantum many-body wave-functions and find they are energetically competitive to other approaches. One interesting application of this architectural exploration is that we show that fully convolution neural network wave-functions can be optimized for one system size and, using identical parameters, produce accurate energies for a range of system sizes. We expect these methods to increase the rate of discovery of novel variational ansatz and bring further insights to the quantum many body problem.


title: Elastic shocks in relativistic rigid rods and balls

abstract: We study the free boundary problem for the "hard phase" material introduced by Christodoulou, both for rods in (1+1)-dimensional Minkowski spacetime and for spherically symmetric balls in (3+1)-dimensional Minkowski spacetime. Unlike Christodoulou, we do not consider a "soft phase", and so we regard this material as an elastic medium, capable of both compression and stretching. We prove that shocks must be null hypersurfaces, and derive the conditions to be satisfied at a free boundary. We solve the equations of motion of the rods explicitly, and we prove existence of solutions to the equations of motion of the spherically symmetric balls for an arbitrarily long (but finite) time, given initial conditions sufficiently close to those for the relaxed ball at rest. In both cases we find that the solutions contain shocks if and only if the pressure or its time derivative do not vanish at the free boundary initially. These shocks interact with the free boundary, causing it to lose regularity.


title: Classifying Snapshots of the Doped Hubbard Model with Machine Learning

abstract: Quantum gas microscopes for ultracold atoms can provide high-resolution real-space snapshots of complex many-body systems. We implement machine learning to analyze and classify such snapshots of ultracold atoms. Specifically, we compare the data from an experimental realization of the two-dimensional Fermi-Hubbard model to two theoretical approaches: a doped quantum spin liquid state of resonating valence bond type, and the geometric string theory, describing a state with hidden spin order. This approach considers all available information without a potential bias towards one particular theory by the choice of an observable and can therefore select the theory which is more predictive in general. Up to intermediate doping values, our algorithm tends to classify experimental snapshots as geometric-string-like, as compared to the doped spin liquid. Our results demonstrate the potential for machine learning in processing the wealth of data obtained through quantum gas microscopy for new physical insights.


title: Floquet time crystals in clock models

abstract: We construct a class of period-$n$-tupling discrete time crystals based on $\mathbb{Z}_n$ clock variables, for all the integers $n$. We consider two classes of systems where this phenomenology occurs, disordered models with short-range interactions and fully connected models. In the case of short-range models we provide a complete classification of time-crystal phases for generic $n$. For the specific cases of $n=3$ and $n=4$ we study in details the dynamics by means of exact diagonalisation. In both cases, through an extensive analysis of the Floquet spectrum, we are able to fully map the phase diagram. In the case of infinite-range models, the mapping onto an effective bosonic Hamiltonian allows us to investigate the scaling to the thermodynamic limit. After a general discussion of the problem, we focus on $n=3$ and $n=4$, representative examples of the generic behaviour. Remarkably, for $n=4$ we find clear evidence of a new crystal-to-crystal transition between period $n$-tupling and period $n/2$-tupling.


title: Black Hole Accretion Discs and Luminous Transients in Failed Supernovae from Non-Rotating Supergiants

abstract: We show that for supergiants, net angular momentum is not a necessary condition for forming accretion discs during core-collapse. Even absent net rotation, convective motions in the outer parts of supergiants generate mean horizontal flows at a given radius with velocities of order 1 km/s; the direction of the mean flow will vary as a function of height through the convection zone. We confirm these analytic estimates using Cartesian Boussinesq convection simulations. These mean horizontal flows lead to a random angular momentum in supergiant convection zones that exceeds that of the last stable circular orbit of a black hole by a factor of order 10. As a result, failed explosions of supergiants - in which the accretion shock onto the neutron star does not revive, leading to black hole formation - may often produce accretion discs that can power day-week (blue supergiants) or week-year (yellow and red supergiants) non-thermal and thermal transients through winds and jets. These transients will be especially time variable because the angular momentum of the accreting material will vary substantially in time. Observed sources such as Swift J1644+57, iPTF14hls, and SN 2018cow, as well as energetic Type II supernovae (OGLE-2014-SN-073) may be produced by this mechanism.


title: Upper bounds on the superfluid stiffness and superconducting $T_c$: Applications to twisted-bilayer graphene and ultra-cold Fermi gases

abstract: Understanding the material parameters that control the superconducting transition temperature $T_c$ is a problem of fundamental importance. In many novel superconductors, phase fluctuations determine $T_c$, rather than the collapse of the pairing amplitude. We derive rigorous upper bounds on the superfluid phase stiffness for multi-band systems, valid in any dimension. This in turn leads to an upper bound on $T_c$ in two dimensions (2D), which holds irrespective of pairing mechanism, interaction strength, or order-parameter symmetry. Our bound is particularly useful for the strongly correlated regime of low-density and narrow-band systems, where mean field theory fails. For a simple parabolic band in 2D with Fermi energy $E_F$, we find that $k_BT_c \leq E_F/8$, an exact result that has direct implications for the 2D BCS-BEC crossover in ultra-cold Fermi gases. Applying our multi-band bound to magic-angle twisted bilayer graphene (MA-TBG), we find that band structure results constrain the maximum $T_c$ to be close to the experimentally observed value. Finally, we discuss the question of deriving rigorous upper bounds on $T_c$ in 3D.


title: Ameso Optimization: a Relaxation of Discrete Midpoint Convexity

abstract: In this paper we introduce the Ameso optimization problem, a special class of discrete optimization problems. We establish its basic properties and investigate the relation between Ameso optimization and the convex optimization. Further, we design an algorithm to solve a multi-dimensional Ameso problem by solving a sequence of one-dimensional Ameso problems. Finally, we demonstrate how the knapsack problem can be solved using the Ameso optimization framework.


title: Dinucleon and Nucleon Decay to Two-Body Final States with no Hadrons in Super-Kamiokande

abstract: Using 0.37 megaton$\cdot$years of exposure from the Super-Kamiokande detector, we search for 10 dinucleon and nucleon decay modes that have a two-body final state with no hadrons. These baryon and lepton number violating modes have the potential to probe theories of unification and baryogenesis. For five of these modes the searches are novel, and for the other five modes we improve the limits by more than one order of magnitude. No significant evidence for dinucleon or nucleon decay is observed and we set lower limits on the partial lifetime of oxygen nuclei and on the nucleon partial lifetime that are above $4\times 10^{33}$ years for oxygen via the dinucleon decay modes and up to about $4 \times 10^{34}$ years for nucleons via the single nucleon decay modes.


title: Extended Higgs sector of 2HDM with real singlet facing LHC data

abstract: We study the phenomenology of the two Higgs doublet model with a real singlet scalar $S$ (2HDMS). The model predicts three CP-even Higgses $h_{1,2,3}$, one CP-odd $A^0$ and a pair of charged Higgs. We discuss the consistency of the 2HDMS with theoretical as well as with all available experimental data. In contrast with previous studies, we focus on the scenario where $h_2$ is the Standard Model (SM) 125 GeV Higgs, while $h_1$ is lighter than $h_2$ which may open a window for Higgs to Higgs decays. We perform an extensive scan into the parameter space of 2HDMS of type I and explore the effect of the singlet-doublet admixture. We found that a large singlet-doublet admixture is still compatible with the recent Higgs data from LHC. Moreover, we show that $h_1$ could be quasi-fermiophobic and would decay dominantly into two photons. We also study in details the consistency of the non-detected decay of $h_2\to h_1 h_1$ with LHC data followed by $h_1\to \gamma \gamma$ which leads to four photons final state at LHC: $pp\to h_2\to h_1 h_1\to 4\gamma$. Using the results of null searches of multi-photons carried by the ATLAS collaboration, we have found that a large area of the parameter space is still allowed. We also demonstrate that various neutral Higgs of the 2HDMS could have several exotic decays.


title: AdaFrame: Adaptive Frame Selection for Fast Video Recognition

abstract: We present AdaFrame, a framework that adaptively selects relevant frames on a per-input basis for fast video recognition. AdaFrame contains a Long Short-Term Memory network augmented with a global memory that provides context information for searching which frames to use over time. Trained with policy gradient methods, AdaFrame generates a prediction, determines which frame to observe next, and computes the utility, i.e., expected future rewards, of seeing more frames at each time step. At testing time, AdaFrame exploits predicted utilities to achieve adaptive lookahead inference such that the overall computational costs are reduced without incurring a decrease in accuracy. Extensive experiments are conducted on two large-scale video benchmarks, FCVID and ActivityNet. AdaFrame matches the performance of using all frames with only 8.21 and 8.65 frames on FCVID and ActivityNet, respectively. We further qualitatively demonstrate learned frame usage can indicate the difficulty of making classification decisions; easier samples need fewer frames while harder ones require more, both at instance-level within the same class and at class-level among different categories.


title: The Plane's The Thing: The Case for Wide-Fast-Deep Coverage of the Galactic Plane and Bulge

abstract: We argue that the exclusion of the Galactic Plane and Bulge from the uniform wide-fast-deep (WFD) LSST survey cadence is fundamentally inconsistent with two of the main science drivers of LSST: Mapping the Milky Way and Exploring the Transient Optical Sky. We outline the philosophical basis for this claim and then describe a number of important science goals that can only be addressed by WFD-like coverage of the Plane and Bulge.


title: Multigrid methods for saddle point problems: Karush-Kuhn-Tucker systems

abstract: We construct multigrid methods for an elliptic distributed optimal control problem that are robust with respect to a regularization parameter. We prove the uniform convergence of the $W$-cycle algorithm and demonstrate the performance of $V$-cycle and $W$-cycle algorithms in two and three dimensions through numerical experiments.


title: Direct Detection Experiments at the Neutrino Dipole Portal Frontier

abstract: Heavy sterile neutrinos are typically invoked to accommodate the observed neutrino masses, by positing a new Yukawa term connecting these new states to the neutrinos in the electroweak doublet. However, given our ignorance of the neutrino sector we should explore additional interactions such sterile neutrinos may have with the SM. In this paper, we study the dimension-5 operator which couples the heavy state to a light neutrino and the photon. We find that the recent XENON1T direct detection data can improve the limits on this "Neutrino Dipole Portal" by up to an order of magnitude over previous bounds. Future direct detection experiments may be able to extend these bounds down to the level probed by SN1987A.


title: Freeform Diffractive Metagrating Design Based on Generative Adversarial Networks

abstract: A key challenge in metasurface design is the development of algorithms that can effectively and efficiently produce high performance devices. Design methods based on iterative optimization can push the performance limits of metasurfaces, but they require extensive computational resources that limit their implementation to small numbers of microscale devices. We show that generative neural networks can train from images of periodic, topology-optimized metagratings to produce high-efficiency, topologically complex devices operating over a broad range of deflection angles and wavelengths. Further iterative optimization of these designs yields devices with enhanced robustness and efficiencies, and these devices can be utilized as additional training data for network refinement. In this manner, generative networks can be trained, with a onetime computation cost, and used as a design tool to facilitate the production of near-optimal, topologically-complex device designs. We envision that such data-driven design methodologies can apply to other physical sciences domains that require the design of functional elements operating across a wide parameter space.


title: Subhalo destruction in the Apostle and Auriga simulations

abstract: N-body simulations make unambiguous predictions for the abundance of substructures within dark matter halos. However, the inclusion of baryons in the simulations changes the picture because processes associated with the presence of a large galaxy in the halo can destroy subhalos and substantially alter the mass function and velocity distribution of subhalos. We compare the effect of galaxy formation on subhalo populations in two state-of-the-art sets of hydrodynamical CDM simulations of Milky Way mass halos, APOSTLE and AURIGA. We introduce a new method for tracking the orbits of subhalos between simulation snapshots that gives accurate results down to a few kiloparsecs from the centre of the halo. Relative to a dark matter-only simulation, the abundance of subhalos in APOSTLE is reduced by 50% near the centre and by 10% within r200. In AURIGA the corresponding numbers are 80% and 40%. The velocity distributions of subhalos are also affected by the presence of the galaxy, much more so in AURIGA than in APOSTLE . The differences on subhalo properties in the two simulations can be traced back to the mass of the central galaxies, which in AURIGA are typically twice as massive as those in APOSTLE . We show that some of the results from previous studies are inaccurate due to systematic errors in the modelling of subhalo orbits near the centre of halos.


title: Infrared Behaviour of SU(2) Gauge Theory with $N_f$ fundamental flavours

abstract: We review our recent results on the infrared behaviour of the SU(2) gauge theory with $N_f$ massless fundamental flavour fermions. We have analyzed the running of the coupling in SU(2) gauge theories with six and eight fermionic flavours using the gradient flow step scaling method. From the running of the coupling, we see a clear indication of an infrared stable fixed point in theories with six and eight flavours. These results are confirmed by our mass spectrum study, where we varied the number of flavours from two to six. We also compute the anomalous dimensions of mass and coupling.


title: Response of the skyrmion lattice in MnSi to cubic magnetocrystalline anisotropies

abstract: We report high-precision small angle neutron scattering of the orientation of the skyrmion lattice in a spherical sample of MnSi under systematic changes of the magnetic field direction. For all field directions the skyrmion lattice may be accurately described as a triple-$\vec{Q}$ state, where the modulus $\vert \vec{Q} \vert$ is constant and the wave vectors enclose rigid angles of $120^{\circ}$. Along a great circle across $\langle 100\rangle$, $\langle 110\rangle$, and $\langle 111\rangle$ the normal to the skyrmion-lattice plane varies systematically by $\pm3^{\circ}$ with respect to the field direction, while the in-plane alignment displays a reorientation by $15^{\circ}$ for magnetic field along $\langle 100\rangle$. Our observations are qualitatively and quantitatively in excellent agreement with an effective potential, that is determined by the symmetries of the tetrahedral point group $T$ and includes contributions up to sixth-order in spin-orbit coupling, providing a full account of the effect of cubic magnetocrystalline anisotropies on the skyrmion lattice in MnSi.


title: Rolled up or crumpled: phases of asymmetric tethered membranes

abstract: We show that inversion-asymmetric tethered membranes exhibit a new double-spiral phase with long range orientational order not present in symmetric membranes. We calculate the universal algebraic spiral shapes of these membranes in this phase. Asymmetry can trigger the crumpling of these membranes as well. In-vitro experiments on lipid, red blood cell membrane extracts, and on graphene coated on one side, could test these predictions.


title: Infinitely Many Primes Using Generating Functions

abstract: In this short paper we present an elementary proof of the infinitude of primes. Our proof is similar in spirit to Euler's proof that the reciprocals of primes diverges and only uses tools from elementary number theory and calculus. In addition, our proof highlights an interesting use of generating functions.


title: Microstructure and formation mechanisms of {\delta}-hydrides in variable grain size Zircaloy-4 studied by electron backscatter diffraction

abstract: Microstructure and crystallography of {\delta} phase hydrides in as-received fine grain and 'blocky' alpha large grain Zircaloy-4 (average grain size ~11 {\mu}m and >200 {\mu}m, respectively) were examined using electron backscatter diffraction. Results suggest that the the matrix-hydride orientation relationship is {0001}{\alpha}||{111}{\delta};<11-20>{\alpha}||<110>{\delta} for all the cases studied. The habit plane of intragranular hydrides and some intergranular hydrides has been found to be {10-17} of the surrounding matrix. The morphology of intergranular hydrides can vary depending upon the angle between the grain boundary and the hydride habit plane. The misfit strain between {\alpha}-Zr and {\delta}-hydride is accommodated by high density of dislocations and twin structures in the hydrides, and a mechanism of twin formation in the hydrides has been proposed. The growth of hydrides across grain boundaries is achieved through an auto-catalytic manner similar to the growth pattern of intragranular hydrides. Easy collective shear along <1-100> makes it possible for hydride nucleation at any grain boundaries, while the process seems to favour grain boundaries with low (<40{\deg}) and high (>80{\deg}) c-axis misorientation angles. Moreover, the angle between the grain boundary and the adjacent basal planes does not influence the propensity for hydride nucleation.


title: Metrological Nonlinear Squeezing Parameter

abstract: The well known metrological linear squeezing parameters (such as quadrature or spin squeezing) efficiently quantify the sensitivity of Gaussian states. Yet, these parameters are insufficient to characterize the much wider class of highly sensitive non-Gaussian states. Here, we introduce a class of metrological nonlinear squeezing parameters obtained by analytical optimization of measurement observables among a given set of accessible (possibly nonlinear) operators. This allows for the metrological characterization of non-Gaussian quantum states of discrete and continuous variables. Our results lead to optimized and experimentally-feasible recipes for high-precision moment-based estimation of a phase parameter and can be used to systematically construct multipartite entanglement and non-classicality witnesses for complex quantum states.


title: Flow Shape Design for Microfluidic Devices Using Deep Reinforcement Learning

abstract: Microfluidic devices are utilized to control and direct flow behavior in a wide variety of applications, particularly in medical diagnostics. A particularly popular form of microfluidics -- called inertial microfluidic flow sculpting -- involves placing a sequence of pillars to controllably deform an initial flow field into a desired one. Inertial flow sculpting can be formally defined as an inverse problem, where one identifies a sequence of pillars (chosen, with replacement, from a finite set of pillars, each of which produce a specific transformation) whose composite transformation results in a user-defined desired transformation. Endemic to most such problems in engineering, inverse problems are usually quite computationally intractable, with most traditional approaches based on search and optimization strategies. In this paper, we pose this inverse problem as a Reinforcement Learning (RL) problem. We train a DoubleDQN agent to learn from this environment. The results suggest that learning is possible using a DoubleDQN model with the success frequency reaching 90% in 200,000 episodes and the rewards converging. While most of the results are obtained by fixing a particular target flow shape to simplify the learning problem, we later demonstrate how to transfer the learning of an agent based on one target shape to another, i.e. from one design to another and thus be useful for a generic design of a flow shape.


title: Optimal Decision Rules for Simple Hypothesis Testing under General Criterion Involving Error Probabilities

abstract: The problem of simple $M-$ary hypothesis testing under a generic performance criterion that depends on arbitrary functions of error probabilities is considered. Using results from convex analysis, it is proved that an optimal decision rule can be characterized as a randomization among at most two deterministic decision rules, of the form reminiscent to Bayes rule, if the boundary points corresponding to each rule have zero probability under each hypothesis. Otherwise, a randomization among at most $M(M-1)+1$ deterministic decision rules is sufficient. The form of the deterministic decision rules are explicitly specified. Likelihood ratios are shown to be sufficient statistics. Classical performance measures including Bayesian, minimax, Neyman-Pearson, generalized Neyman-Pearson, restricted Bayesian, and prospect theory based approaches are all covered under the proposed formulation. A numerical example is presented for prospect theory based binary hypothesis testing.


title: Constraints on $U(1)_{L_\mu-L_\tau}$ from LHC Data

abstract: In this study, we apply LHC data to constrain the extension of the Standard Model by an anomaly-free $U(1)_{L_\mu-L_\tau}$ gauge group; this model contains a new gauge boson ($Z^\prime$) and a scalar dark matter particle ($\phi_{\rm DM}$). We recast a large number of LHC analyses from ATLAS and CMS of multi-lepton final states. We find that for $10$ GeV $< m_{Z^\prime} < 60$ GeV the strongest constraint comes from a dedicated $Z^\prime$ search in the $4\mu$ final state by the CMS collaboration; for larger $Z^\prime$ masses, searches for final states with three leptons plus missing $E_T$ are more sensitive. Searches for final states with two leptons and missing $E_T$, which are sensitive to $Z^\prime$ decays into dark matter particles, can only probe regions of parameter space that are excluded by searches in the $3$ and $4$ lepton channels. The combination of LHC data excludes values of $Z^\prime$ mass and coupling constant that can explain the deficit in $g_\mu-2$ for $4$ GeV $\leq m_{Z^\prime} \leq 500$ GeV. However, for much of this range the LHC bound is weaker than the bound that can be derived from searches for trident events in neutrino-nucleus scattering.


title: On the Performance of Reed-Muller Codes with respect to Random Errors and Erasures

abstract: This work proves new results on the ability of binary Reed-Muller codes to decode from random errors and erasures. We obtain these results by proving improved bounds on the weight distribution of Reed-Muller codes of high degrees. Specifically, given weight $\beta \in (0,1)$ we prove an upper bound on the number of codewords of relative weight at most $\beta$. We obtain new results in two different settings: for weights $\beta < 1/2$ and for weights that are close to $1/2$. Our new bounds on the weight distribution imply that RM codes with $m$ variables and degree $\gamma m$, for some explicit constant $\gamma$, achieve capacity for random erasures (i.e. for the binary erasure channel) and for random errors (for the binary symmetric channel). Earlier, it was known that RM codes achieve capacity for the binary symmetric channel for degrees $r = o(m)$. For the binary erasure channel it was known that RM codes achieve capacity for degree $o(m)$ or $r \in [m/2 \pm O(\sqrt{m})]$. Thus, our result provide a new range of parameters for which RM achieve capacity for these two well studied channels. In addition, our results imply that for every $\epsilon > 0$ (in fact we can get up to $\epsilon = \Omega\left(\sqrt{\frac{\log m}{m}}\right)$) RM codes of degree $r<(1/2-\epsilon)m$ can correct a fraction of $1-o(1)$ random erasures with high probability. We also show that, information theoretically, such codes can handle a fraction of $1/2-o(1)$ random errors with high probability. Thus, for example, given noisy evaluations of a degree $0.499m$ polynomial, it is possible to interpolate it even if a random $0.499$ fraction of the evaluations were corrupted, with high probability. While the $o(1)$ terms are not the correct ones to ensure capacity, these results show that RM codes of such degrees are in some sense close to achieving capacity.


title: On estimates for fully nonlinear elliptic equations with Neumann boundary conditions on Riemannian Manifolds

abstract: We derive gradient and second order {\em a priori} estimates for solutions of the Neumann problem for a general class of fully nonlinear elliptic equations on compact Riemannian manifolds with boundary. These estimates yield regularity and existence results.


title: An Adaptive Finite Element DtN Method for Maxwell's Equations in Biperiodic Structures

abstract: Consider the diffraction of an electromagnetic plane wave by a biperiodic structure where the wave propagation is governed by the three-dimensional Maxwell equations. Based on transparent boundary condition, the grating problem is formulated into a boundary value problem in a bounded domain. Using a duality argument technique, we derive an a posteriori error estimate for the finite element method with the truncation of the nonlocal Dirichlet-to-Neumann (DtN) boundary operator. The a posteriori error consists of both the finite element approximation error and the truncation error of boundary operator which decays exponentially with respect to the truncation parameter. An adaptive finite element algorithm is developed with error controlled by the a posterior error estimate, which determines the truncation parameter through the truncation error and adjusts the mesh through the finite element approximation error. Numerical experiments are presented to demonstrate the competitive behavior of the proposed adaptive method.