title: Newtonian self-gravity in trapped quantum systems and experimental tests

abstract: No experimental evidence exists, to date, whether or not the gravitational field must be quantised. Theoretical arguments in favour of quantisation are inconclusive. The most straightforward alternative to quantum gravity, a coupling between classical gravity and quantum matter according to the semi-classical Einstein equations, yields a nonlinear modification of the Schr\"odinger equation. Here, effects of this so-called Schr\"odinger-Newton equation are discussed, which allow for technologically feasible experimental tests.


title: Optimal Control for Nonlinear Hybrid Systems via Convex Relaxations

abstract: This paper considers the optimal control for hybrid systems whose trajectories transition between distinct subsystems when state-dependent constraints are satisfied. Though this class of systems is useful while modeling a variety of physical systems undergoing contact, the construction of a numerical method for their optimal control has proven challenging due to the combinatorial nature of the state-dependent switching and the potential discontinuities that arise during switches. This paper constructs a convex relaxation-based approach to solve this optimal control problem. Our approach begins by formulating the problem in the space of relaxed controls, which gives rise to a linear program whose solution is proven to compute the globally optimal controller. This conceptual program is solved by constructing a sequence of semidefinite programs whose solutions are proven to converge from below to the true solution of the original optimal control problem. Finally, a method to synthesize the optimal controller is developed. Using an array of examples, the performance of the proposed method is validated on problems with known solutions and also compared to a commercial solver.


title: A storm is Coming: A Modern Probabilistic Model Checker

abstract: We launch the new probabilistic model checker storm. It features the analysis of discrete- and continuous-time variants of both Markov chains and MDPs. It supports the PRISM and JANI modeling languages, probabilistic programs, dynamic fault trees and generalized stochastic Petri nets. It has a modular set-up in which solvers and symbolic engines can easily be exchanged. It offers a Python API for rapid prototyping by encapsulating storm's fast and scalable algorithms. Experiments on a variety of benchmarks show its competitive performance.


title: Magnetism of monomer MnO and heterodimer FePt@MnO nanoparticles

abstract: We report about the magnetic properties of antiferromagnetic (AF) MnO nanoparticles (NPs) with different sizes (6-19nm). Using a combination of polarized neutron scattering and magnetometry we were able to resolve previously observed peculiarities. Magnetometry, on the one hand, reveals a peak in the zero field cooled (ZFC) magnetization curves at low temperatures (~25K) but no feature around the N\'eel temperature at 118K. On the other hand, polarized neutron scattering shows the expected behavior of the AF order parameter vanishing around 118K. Moreover, hysteresis curves measured at various temperatures reveal an exchange bias effect indicating a coupling of an AF core to a ferromagnetic (FM)-like shell. ZFC data measured at various fields exclude a purely superparamagnetic (SPM) scenario. We conclude that the magnetic behavior of MnO particles can be explained by a superposition of SPM-like thermal fluctuations of the AF-N\'eel vector inside the AF core and a strong magnetic coupling to a ferrimagnetic Mn$_2$O$_3$ or Mn$_3$O$_4$ shell. In addition, we have studied heterodimer ('Janus') particles, where a FM FePt particle is attached to the AF MnO particle. Via the exchange bias effect, the magnetic moment of the FePt subunit is stabilized by the MnO.


title: Terminal-Pairability in Complete Bipartite Graphs

abstract: We investigate the terminal-pairibility problem in the case when the base graph is a complete bipartite graph, and the demand graph is also bipartite with the same color classes. We improve the lower bound on maximum value of $\Delta(D)$ which still guarantees that the demand graph $D$ is terminal-pairable in this setting. We also prove a sharp theorem on the maximum number of edges such a demand graph can have.


title: Cosmology with XMM galaxy clusters: the X-CLASS/GROND catalogue and photometric redshifts

abstract: The XMM Cluster Archive Super Survey (X-CLASS) is a serendipitously-detected X-ray-selected sample of 845 galaxy clusters based on 2774 XMM archival observations and covering approximately 90 deg$^2$ spread across the high-Galactic latitude ($|b|>20$ deg) sky. The primary goal of this survey is to produce a well-selected sample of galaxy clusters on which cosmological analyses can be performed. This article presents the photometric redshift followup of a high signal-to-noise subset of 266 of these clusters with declination $\delta<+20$ deg with GROND, a seven channel ($grizJHK$) simultaneous imager on the MPG 2.2m telescope at the ESO La Silla Observatory. We use a newly developed technique based on the red sequence colour-redshift relation, enhanced with information coming from the X-ray detection to provide photometric redshifts for this sample. We determine photometric redshifts for 236 clusters, finding a median redshift of $z=0.39$ with an accuracy of $\Delta z = 0.02 (1+z)$ when compared to a sample of 76 spectroscopically confirmed clusters. We also compute X-ray luminosities for the entire sample and find a median bolometric luminosity of $7.2\times10^{43} \mathrm{erg\ s^{-1}}$ and a median temperature 2.9 keV. We compare our results to the XMM-XCS and XMM-XXL surveys, finding good agreement in both samples. The X-CLASS catalogue is available online at http://xmm-lss.in2p3.fr:8080/l4sdb/.


title: Optimal design problems for the first $p-$fractional eigenvalue with mixed boundary conditions

abstract: In this paper we study an optimal shape design problem for the first eigenvalue of the fractional $p-$laplacian with mixed boundary conditions. The optimization variable is the set where the Dirichlet condition is imposed (that is restricted to have measure equal than a prescribed quantity, $\alpha$). We show existence of an optimal design and analyze the asymptotic behavior when the fractional parameter $s\uparrow 1$ obtaining asymptotic bounds that are independent of $\alpha$.


title: Acceleration of the Implicit-Explicit Non-hydrostatic Unified Model of the Atmosphere (NUMA) on Manycore Processors

abstract: We present the acceleration of an IMplicit-EXplicit (IMEX) non-hydrostatic atmospheric model on manycore processors such as GPUs and Intel's MIC architecture. IMEX time integration methods sidestep the constraint imposed by the Courant-Friedrichs-Lewy condition on explicit methods through corrective implicit solves within each time step. In this work, we implement and evaluate the performance of IMEX on manycore processors relative to explicit methods. Using 3D-IMEX at Courant number C=15 , we obtained a speedup of about 4X relative to an explicit time stepping method run with the maximum allowable C=1. In addition, we demonstrate a much larger speedup of 100X at C=150 using 1D-IMEX due to the unconditional stability of the method in the vertical direction. Several improvements on the IMEX procedure were necessary in order to outperform our results with explicit methods: a) reducing the number of degrees of freedom of the IMEX formulation by forming the Schur complement; b) formulating a horizontally-explicit vertically-implicit (HEVI) 1D-IMEX scheme that has a lower workload and potentially better scalability than 3D-IMEX; c) using high-order polynomial preconditioners to reduce the condition number of the resulting system; d) using a direct solver for the 1D-IMEX method by performing and storing LU factorizations once to obtain a constant cost for any Courant number. Without all of these improvements, explicit time integration methods turned out to be difficult to beat. We discuss in detail the IMEX infrastructure required for formulating and implementing efficient methods on manycore processors. Finally, we validate our results with standard benchmark problems in NWP and evaluate the performance and scalability of the IMEX method using up to 4192 GPUs and 16 Knights Landing processors.


title: Analysis of the Extended Coupled-Cluster Method in Quantum Chemistry

abstract: The mathematical foundation of the so-called extended coupled-cluster method for the solution of the many-fermion Schr\"odinger equation is here developed. We prove an existence and uniqueness result, both in the full infinite-dimensional amplitude space as well as for discretized versions of it. The extended coupled-cluster method is formulated as a critical point of an energy function using a generalization of the Rayleigh-Ritz principle: the bivariational principle. This gives a quadratic bound for the energy error in the discretized case. The existence and uniqueness results are proved using a type of monotonicity property for the flipped gradient of the energy function. Comparisons to the analysis of the standard coupled-cluster method is made, and it is argued that the bivariational principle is a useful tool, both for studying coupled-cluster type methods, and for developing new computational schemes in general.

categories: math.NA physics.atom-ph physics.chem-ph physics.comp-ph quant-ph


title: Non-Hamiltonian isotopic Lagrangians on the one-point blow-up of CP^2

abstract: We show that two Hamiltonian isotopic Lagrangians in (CP^2,\omega_\textup{FS}) induce two Lagrangian submanifolds in the one-point blow-up (\widetilde{CP}^2,\widetilde{\omega}_\rho) that are not Hamiltonian isotopic. Furthermore, we show that for any integer k>1 there are k Hamiltonian isotopic Lagrangians in (CP^2,\omega_\textup{FS}) that induce k Lagrangian submanifolds in the one-point blow-up such that no two of them are Hamiltonian isotopic.


title: Superresolution via Structured Illumination Quantum Correlation Microscopy (SIQCM)

abstract: We propose to use intensity correlation microscopy in combination with structured illumination to image quantum emitters that exhibit antibunching with a spatial resolution reaching far beyond the Rayleigh limit. Combining intensity measurements and intensity auto correlations up to order $m$ creates an effective PSF with FWHM shrunk by the factor $\sqrt{m}$. Structured Illumination microscopy on the other hand introduces a resolution improvement of factor 2 by use of the principle of moir\'e fringes. Here, we show that for linear low-intensity excitation and linear optical detection the simultaneous use of both techniques leads to an in theory unlimited resolution power with the improvement scaling favorably as $m + \sqrt{m}$ in dependence of the correlation order $m$. Hence, yielding this technique to be of interest in microscopy for imaging a variety of samples including biological ones. We present the underlying theory and simulations demonstrating the highly increased spatial superresolution, and point out requirements for an experimental implementation.


title: Arbitrarily Tight Bounds on a Singularly Perturbed Linear-Quadratic Optimal Control Problem

abstract: We calculate arbitrarily tight upper and lower bounds on an unconstrained control, linear-quadratic, singularly perturbed optimal control problem whose exact solution is computationally intractable. It is well known that for the aforementioned problem, an approximate solution $\bar{V}^N(\epsilon)$ can be constructed such that it is asymptotically equivalent in $\epsilon$ to the solution $V(\epsilon)$ of the singularly perturbed problem in the sense that $|V(\epsilon)-\bar{V}^N(\epsilon)| =O(\epsilon^{N+1})$ for any integer $N\geq0$ as $\epsilon \rightarrow 0$. For this approximation to be considered useful, the parameter $\epsilon$ is typically restricted to be in some sufficiently small set; however, for values of $\epsilon$ outside this set, a poor approximation can result. We improve on this approximation by incorporating a duality theory into the singularly perturbed optimal control problem and derive an upper bound $\chi^N_u(\epsilon)$ and a lower bound $\chi^N_l(\epsilon)$ of $V(\epsilon)$ that hold for arbitrary $\epsilon$ and, furthermore, satisfy the inequality $|\chi^N_u(\epsilon)-\chi^N_l(\epsilon)|=O(\epsilon^{N+1})$ for any integer $N \geq 0$ as $\epsilon \rightarrow 0$.


title: Interaction Between Two Ring Shaped YBaCuO Superconductors with currents ia and ib

abstract: Two rings with currents ia and ib become magnets that can either attract or repel each other. If mechanical work is extracted from this ring system, either by bringing closer the rings during an attraction or by separating them as the result of a repulsion, the energy of the rings, the only existing source of energy in the system, decreases: $\frac{1}{2}LI_{a}^{2}+\frac{1}{2}LI_{b}^{2}$ If an energy either of attraction or repulsion is communicated to the system, (working against the forces) its energy increases, i.e. the currents in the rings increase. This happens because the currents existing in the rings are free, i.e. not linked to any device which imposes a current (power supply, constant current). Similarly, the effect of electromagnetic induction can increase or decrease the currents with the corresponding mechanical work linked to these developments. Not only does the current in the rings decrease as they get closer during the attraction, but it also can be cancelled and, if the effect continues, it may even change its direction. An application that already exists is the levitation train. It is attracted by the rail and held on it, but if it is very close to this due to the weight of the goods a strong repulsion (levitation) occurs. Another application would be the accumulation of energy in these rings, where they can spend large amounts of current, not only because of its high critical currents but also by the fact that the magnetic flux remains constant, cannot increase. Another advantage in this storage system is its reversibility in the exchange of electromagnetic energy and mechanical work. In the energy storage system, the greater the current flowing through the ring, the greater its energy is. In order to attain this, there has been made a study of the ideal conditions to store as much energy as possible.


title: Parameterized Algorithms for Recognizing Monopolar and 2-Subcolorable Graphs

abstract: A graph $G$ is a $(\Pi_A,\Pi_B)$-graph if $V(G)$ can be bipartitioned into $A$ and $B$ such that $G[A]$ satisfies property $\Pi_A$ and $G[B]$ satisfies property $\Pi_B$. The $(\Pi_{A},\Pi_{B})$-Recognition problem is to recognize whether a given graph is a $(\Pi_A,\Pi_B)$-graph. There are many $(\Pi_{A},\Pi_{B})$-Recognition problems, including the recognition problems for bipartite, split, and unipolar graphs. We present efficient algorithms for many cases of $(\Pi_A,\Pi_B)$-Recognition based on a technique which we dub inductive recognition. In particular, we give fixed-parameter algorithms for two NP-hard $(\Pi_{A},\Pi_{B})$-Recognition problems, Monopolar Recognition and 2-Subcoloring. We complement our algorithmic results with several hardness results for $(\Pi_{A},\Pi_{B})$-Recognition.


title: Don't cry over spilled records: Memory elasticity of data-parallel applications and its application to cluster scheduling

abstract: Understanding the performance of data-parallel workloads when resource-constrained has significant practical importance but unfortunately has received only limited attention. This paper identifies, quantifies and demonstrates memory elasticity, an intrinsic property of data-parallel tasks. Memory elasticity allows tasks to run with significantly less memory that they would ideally want while only paying a moderate performance penalty. For example, we find that given as little as 10% of ideal memory, PageRank and NutchIndexing Hadoop reducers become only 1.2x/1.75x and 1.08x slower. We show that memory elasticity is prevalent in the Hadoop, Spark, Tez and Flink frameworks. We also show that memory elasticity is predictable in nature by building simple models for Hadoop and extending them to Tez and Spark. To demonstrate the potential benefits of leveraging memory elasticity, this paper further explores its application to cluster scheduling. In this setting, we observe that the resource vs. time trade-off enabled by memory elasticity becomes a task queuing time vs task runtime trade-off. Tasks may complete faster when scheduled with less memory because their waiting time is reduced. We show that a scheduler can turn this task-level trade-off into improved job completion time and cluster-wide memory utilization. We have integrated memory elasticity into Apache YARN. We show gains of up to 60% in average job completion time on a 50-node Hadoop cluster. Extensive simulations show similar improvements over a large number of scenarios.


title: Dilution of Ferromagnets via a Random Graph-based Strategy

abstract: The dynamics and behavior of ferromagnets have a great relevance even beyond the domain of statistical physics. In this work, we propose a Monte Carlo method, based on random graphs, for modeling their dilution. In particular, we focus on ferromagnets with dimension $D \ge 4$, which can be approximated by the Curie-Weiss model. Since the latter has as graphic counterpart a complete graph, a dilution can be in this case viewed as a pruning process. Hence, in order to exploit this mapping, the proposed strategy uses a modified version of the Erd\H{o}s-Renyi graph model. In doing so, we are able both to simulate a continuous dilution, and to realize diluted ferromagnets in one step. The proposed strategy is studied by means of numerical simulations, aimed to analyze main properties and equilibria of the resulting diluted ferromagnets. To conclude, we also provide a brief description of further applications of our strategy in the field of complex networks.


title: Quantitative stratification for some free-boundary problems

abstract: In this paper we prove the rectifiability of and measure bounds on the singular set of the free boundary for minimizers of a functional first considered by Alt-Caffarelli. Our main tools are the Quantitative Stratification and Rectifiable-Reifenberg framework of Naber-Valtorta, which allow us to do a type of "effective dimension-reduction." The arguments are sufficiently robust that they apply to a broad class of related free boundary problems as well.


title: Pentazole and Ammonium Pentazolate: Crystalline Hydro-Nitrogens at High Pressure

abstract: Two new crystalline compounds, pentazole (N_{5}H) and ammonium pentazolate (NH_{4})(N_{5}), both featuring cyclo-{\rm N_{5}^{-}} are discovered using first principles evolutionary search of the nitrogen-rich portion of the hydro-nitrogen binary phase diagram (N_{x}H_{y}, x\geqy) at high pressures. Both crystals consist of the pentazolate N_{5}^{-} anion and ammonium NH_{4}^{+} or hydrogen H^{+} cations. These two crystals are predicted to be thermodynamically stable at pressures above 30 GPa for (NH_{4})(N_{5}) and 50 GPa for pentazole N_{5}H. The chemical transformation of ammonium azide (NH_{4})(N_{3}) mixed with di-nitrogen (N_{2}) to ammonium pentazolate (NH_{4})(N_{5}) is predicted to become energetically favorable above 12.5 GPa. To assist in identification of newly synthesized compounds in future experiments, the Raman spectra of both crystals are calculated and mode assignments are made as a function of pressure up to 75 GPa.


title: The Biot-Savart operator of a bounded domain

abstract: We construct the analog of the Biot-Savart integral for bounded domains. Specifically, we show that the velocity field of an incompressible fluid with tangency boundary conditions on a bounded domain can be written in terms of its vorticity using an integral kernel $K_\Omega(x,y)$ that has an inverse-square singularity on the diagonal.


title: Measurements of CKM angles at Belle

abstract: In this review recent studies on $CP$ violation and related hadronic $B$ decays by the Belle experiment, in particular measurements of CKM angles $\phi_1$ and $\phi_2$ are reported.


title: Analysis of extreme values with random location

abstract: Analysis of the rare and extreme values through statistical modeling is an important issue in economical crises, climate forecasting, and risk management of financial portfolios. Extreme value theory provides the probability models needed for statistical modeling of the extreme values. There are generally two ways to identifying the extreme values in a data set, the block-maxima and the peak-over threshold method. The block-maxima method uses the Generalized Extreme Value distribution and the peak-over threshold method uses the Generalized Pareto distribution. It is common that the location of these distributions kept fixed. It is possible that some unobserved variables produce heterogeneity in the location of the assumed distribution. In this article we focus on modeling this unobserved heterogeneity in block-maxima method.


title: A Nonparametric Bayesian Approach for Sparse Sequence Estimation

abstract: A nonparametric Bayes approach is proposed for the problem of estimating a sparse sequence based on Gaussian random variables. We adopt the popular two-group prior with one component being a point mass at zero, and the other component being a mixture of Gaussian distributions. Although the Gaussian family has been shown to be suboptimal for this problem, we find that Gaussian mixtures, with a proper choice on the means and mixing weights, have the desired asymptotic behavior, e.g., the corresponding posterior concentrates on balls with the desired minimax rate. To achieve computation efficiency, we propose to obtain the posterior distribution using a deterministic variational algorithm. Empirical studies on several benchmark data sets demonstrate the superior performance of the proposed algorithm compared to other alternatives.


title: Low-energy phonons in Bi$_2$Sr$_2$CaCu$_2$O$_{8+\delta}$ and their possible interaction with electrons measured by inelastic neutron scattering

abstract: Electron-phonon interaction in copper oxide superconductors is still enigmatic. Strong coupling for certain optic phonons is now well established experimentally, but theoretical understanding is challenging. Scattering of electrons near the Fermi surface by the longitudinal acoustic (LA) phonons is expected from basic theory because these phonons modulate electron density. We used inelastic neutron scattering on a large single crystal sample of optimally-doped Bi$_2$Sr$_2$CaCu$_2$O$_{8+\delta}$ to show that low-energy LA phonons could couple to electronic density fluctuations only at small phonon wavevectors, which naturally limits any interaction to forward scattering. Such scattering should not be pairbreaking in the case of the d-wave gap. We also found that previously the reported low energy phonon spectral weight half-way to the zone boundary is consistent with conventional lattice dynamics and does not reflect an incipient charge density wave.


title: Adaptive Reduced-Order Modeling for Non-Linear Fluid-Structure Interaction

abstract: We present an adaptive reduced-order model for the efficient time-resolved simulation of fluid-structure interaction problems with complex and non-linear deformations. The model is based on repeated linearizations of the structural balance equations. Upon each linearization step, the number of unknowns is strongly decreased by using modal reduction, which leads to a substantial gain in computational efficiency. Through adaptive re-calibration and truncation augmentation whenever a non-dimensional deformation threshold is exceeded, we ensure that the reduced modal basis maintains arbitrary accuracy for small and large deformations. Our novel model is embedded into a partitioned, loosely coupled finite volume - finite element framework, in which the structural interface motion within the Eulerian fluid solver is accounted for by a conservative cut-element immersed-boundary method. Applications to the aeroelastic instability of a flat plate at supersonic speeds, to an elastic panel placed within a shock tube, and to the shock induced buckling of an inflated thin semi-sphere demonstrate the efficiency and accuracy of the method.


title: On the Relevance of Auditory-Based Gabor Features for Deep Learning in Automatic Speech Recognition

abstract: Previous studies support the idea of merging auditory-based Gabor features with deep learning architectures to achieve robust automatic speech recognition, however, the cause behind the gain of such combination is still unknown. We believe these representations provide the deep learning decoder with more discriminable cues. Our aim with this paper is to validate this hypothesis by performing experiments with three different recognition tasks (Aurora 4, CHiME 2 and CHiME 3) and assess the discriminability of the information encoded by Gabor filterbank features. Additionally, to identify the contribution of low, medium and high temporal modulation frequencies subsets of the Gabor filterbank were used as features (dubbed LTM, MTM and HTM respectively). With temporal modulation frequencies between 16 and 25 Hz, HTM consistently outperformed the remaining ones in every condition, highlighting the robustness of these representations against channel distortions, low signal-to-noise ratios and acoustically challenging real-life scenarios with relative improvements from 11 to 56% against a Mel-filterbank-DNN baseline. To explain the results, a measure of similarity between phoneme classes from DNN activations is proposed and linked to their acoustic properties. We find this measure to be consistent with the observed error rates and highlight specific differences on phoneme level to pinpoint the benefit of the proposed features.


title: Galaxy chemical evolution models: The role of molecular gas formation

abstract: In our grid of multiphase chemical evolution models (Moll\'a & D\'iaz, 2005), star formation in the disk occurs in two steps: first, molecular gas forms, and then stars are created by cloud-cloud collisions or interactions of massive stars with the surrounding molecular clouds. The formation of both molecular clouds and stars are treated through the use of free parameters we refer to as efficiencies. In this work we modify the formation of molecular clouds based on several new prescriptions existing in the literature, and we compare the results obtained for a chemical evolution model of the Milky Way Galaxy regarding the evolution of the Solar region, the radial structure of the Galactic disk, and the ratio between the diffuse and molecular components, HI/H$_2$. Our results show that the six prescriptions we have tested reproduce fairly consistent most of the observed trends, differing mostly in their predictions for the (poorly-constrained) outskirts of the Milky Way and the evolution in time of its radial structure. Among them, the model proposed by Ascasibar et al. (2017), where the conversion of diffuse gas into molecular clouds depends on the local stellar and gas densities as well as on the gas metallicity, seems to provide the best overall match to the observed data.


title: Anomalies in (semi)-leptonic $B$ decays $B^{\pm} \to \tau^{\pm} \nu$, $B^{\pm} \to D \tau^{\pm} \nu$ and $B^{\pm} \to D^* \tau^{\pm} \nu$, and possible resolution with sterile neutrino

abstract: The universality of the weak interactions can be tested in semileptonic $b \to c$ transitions, and in particular in the ratios $R(D^{(*)}) \equiv \Gamma(B \rightarrow D^{(*)} \tau \nu )/\Gamma(B \rightarrow D^{(*)} \ell \nu )$ (where $\ell = \mu$ or $e$). Due to the recent differences between the experimental measurements of these observables by BaBar, Belle and LHCb on the one hand and the Standard Model predicted values on the other hand, we study the predicted ratios $R(D^{(*)}) = \Gamma(B \rightarrow D^{(*)} \tau +$ "missing" $)/\Gamma(B \rightarrow D^{(*)} \ell \nu )$ in scenarios with an additional sterile heavy neutrino of mass $\sim 1$ GeV. Further, we evaluate the newly defined ratio $R(0) \equiv \Gamma(B \to \tau +$ "missing" $)/\Gamma(B \to \mu \nu )$ in such scenarios, in view of the future possibilities of measuring the quantity at Belle-II.


title: The Standard Model Algebra - Leptons, Quarks, and Gauge from the Complex Clifford Algebra Cl6

abstract: A simple geometric algebra is shown to contain automatically the leptons and quarks of a family of the Standard Model, and the electroweak and color gauge symmetries, without predicting extra particles and symmetries. The algebra is already naturally present in the Standard Model, in two instances of the Clifford algebra $\mathbb{C}\ell_6$, one being algebraically generated by the Dirac algebra and the weak symmetry generators, and the other by a complex three-dimensional representation of the color symmetry, which generates a Witt decomposition which leads to the decomposition of the algebra into ideals representing leptons and quarks. The two instances being isomorphic, the minimal approach is to identify them, resulting in the model proposed here. The Dirac and Lorentz algebras appear naturally as subalgebras acting on the ideals representing leptons and quarks. The resulting representations on the ideals are invariant to the electromagnetic and color symmetries, which are generated by the bivectors of the algebra. The electroweak symmetry is also present, and it is already broken by the geometry of the algebra. The model predicts a bare Weinberg angle $\theta_W$ given by $\sin^2\theta_W=0.25$. The model shares common ideas with previously known models, particularly with Chisholm and Farwell, 1996, Trayling and Baylis, 2004, and Furey, 2016.


title: Spectral Properties of Continuum Fibonacci Schr\"odinger Operators

abstract: We study continuum Schr\"odinger operators on the real line whose potentials are comprised of two compactly supported square-integrable functions concatenated according to an element of the Fibonacci substitution subshift over two letters. We show that the Hausdorff dimension of the spectrum tends to one in the small-coupling and high-energy regimes, regardless of the shape of the potential pieces.


title: Beyond the perturbative description of the nonlinear optical response of low-index materials

abstract: We show that standard approximations in nonlinear optics are violated for situations involving a small value of the linear refractive index. Consequently, the conventional equation for the intensity-dependent refractive index, $n(I) = n_0 + n_2 I$, becomes inapplicable in epsilon-near-zero and low-index media, even in the presence of only third-order effects. For the particular case of indium tin oxide, we find that the $\chi^{(3)}$, $\chi^{(5)}$ and $\chi^{(7)}$ contributions to refraction eclipse the linear term; thus, the nonlinear response can no longer be interpreted as a perturbation in these materials. Although the response is non-perturbative, we find no evidence that the power series expansion of the material polarization diverges.


title: GRB 161219B-SN 2016jca: a powerful stellar collapse

abstract: We report observations and analysis of the nearby gamma-ray burst GRB\,161219B (redshift $z=0.1475$) and the associated Type Ic supernova (SN) 2016jca. GRB\,161219B had an isotropic gamma-ray energy of $\sim 1.6 \times 10^{50}$\,erg. Its afterglow is likely refreshed at an epoch preceding the first photometric points (0.6\,d), which slows down the decay rates. Combined analysis of the SN light curve and multiwavelength observations of the afterglow suggest that the GRB jet was broad during the afterglow phase (full opening angle $\sim 42^\circ \pm 3^\circ$). Our spectral series shows broad absorption lines typical of GRB supernovae (SNe), which testify to the presence of material with velocities up to $\sim 0.25$c. The spectrum at 3.73\,d allows for the very early identification of a SN associated with a GRB. Reproducing it requires a large photospheric velocity ($35,000 \pm 7000$\,\kms). The kinetic energy of the SN is estimated through models to be \KE $\approx 4 \times 10^{52}$\,erg in spherical symmetry. The ejected mass in the explosion was \Mej $\approx 6.5 \pm 1.5$\,\Msun, much less than that of other GRB-SNe, demonstrating diversity among these events. The total amount of \Nifs\ in the explosion was $0.27 \pm 0.05$\,\Msun. The observed spectra require the presence of freshly synthesised \Nifs\ at the highest velocities, at least 3 times more than a standard GRB-SN. We also find evidence for a decreasing \Nifs\ abundance as a function of decreasing velocity. This suggests that SN\,2016jca was a highly aspherical explosion viewed close to on-axis, powered by a compact remnant. Applying a typical correction for asymmetry, the energy of SN\,2016jca was $\sim$ (1--3) $\times 10^{52}$\,erg, confirming that most of the energy produced by GRB-SNe goes into the kinetic energy of the SN ejecta.


title: Emergent universal critical behavior of the 2D $N$-color Ashkin-Teller model in the presence of correlated disorder

abstract: We study the critical behavior of the 2D $N$-color Ashkin-Teller model in the presence of random bond disorder whose correlations decays with the distance $r$ as a power-law $r^{-a}$. We consider the case when the spins of different colors sitting at the same site are coupled by the same bond and map this problem onto the 2D system of $N/2$ flavors of interacting Dirac fermions in the presence of correlated disorder. Using renormalization group we show that for $N=2$, a "weakly universal" scaling behavior at the continuous transition becomes universal with new critical exponents. For $N>2$, the first-order phase transition is rounded by the correlated disorder and turns into a continuous one.


title: The Greater Taurus-Auriga Ecosystem I: There Is A Distributed Older Population

abstract: The Taurus-Auriga association and its associated molecular cloud are a benchmark population for studies of star and planet formation. The census of Taurus-Auriga has been assembled over seven decades and has inherited the biases, incompleteness, and systematic uncertainties of the input studies. The notably unusual shape of the inferred IMF and the existence of several isolated disk-bearing stars suggest that additional (likely disk-free) members might remain to be discovered. We therefore have begun a global reassessment of the membership of Taurus-Auriga that exploits new data and better definitions of youth and kinematic membership. As a first step, we reconsider the membership of all disk-free candidate members from the literature with spectral type $\ge$F0, $3^h50^m<\alpha<5^h40^m$, and $14^{\circ}<\delta<34^{\circ}$. We combine data from the literature with Keck/HIRES and UH88/SNIFS spectra to test the membership of these candidates using HR diagram positions, proper motions, RVs, H$\alpha$, lithium, and surface gravity. We find 218 confirmed or likely Taurus members, 160 confirmed or likely interlopers, and only 18 that still lack sufficient evidence to draw firm conclusions. A significant fraction of these stars (81/218=37%) are not included in the most recent canonical member lists. Intriguingly, there are few additional members in the immediate vicinity of the molecular clouds, preserving the IMFs that have been deemed anomalous in past work. Many of the likely Taurus members are distributed broadly across the search area. When combined with known disk hosts, our updated census reveals two regimes: a high-density population with a high disk fraction (indicative of youth) that broadly traces the molecular clouds, and a low-density population with low disk fraction (hence likely older) that most likely represents previous generations of star formation.


title: BranchHull: Convex bilinear inversion from the entrywise product of signals with known signs

abstract: We consider the bilinear inverse problem of recovering two vectors, $x$ and $w$, in $\mathbb{R}^L$ from their entrywise product. For the case where the vectors have known signs and belong to known subspaces, we introduce the convex program BranchHull, which is posed in the natural parameter space that does not require an approximate solution or initialization in order to be stated or solved. Under the structural assumptions that $x$ and $w$ are members of known $K$ and $N$ dimensional random subspaces, we present a recovery guarantee for the noiseless case and a noisy case. In the noiseless case, we prove that the BranchHull recovers $x$ and $w$ up to the inherent scaling ambiguity with high probability when $L\ \gg\ 2(K+N)$. The analysis provides a precise upper bound on the coefficient for the sample complexity. In a noisy case, we show that with high probability the BranchHull is robust to small dense noise when $L = \Omega(K+N)$. BranchHull is motivated by the sweep distortion removal task in dielectric imaging, where one of the signals is a nonnegative reflectivity, and the other signal lives in a known wavelet subspace. Additional potential applications are blind deconvolution and self-calibration.


title: 3DNA Printer: A Tool for Automated DNA Origami

abstract: In the last two decades, DNA self-assembly has grown into a major area of research attracting people from diverse background. It has numerous potential applications such as targeted drug delivery, artificial photosynthesis etc. In the last decade, another area received wide attention known as DNA origami, where using M13 virus and carefully designed staple strands one can fold the DNA into desired 2-D and 3-D shapes. In 2016, a group of researchers at MIT have developed an automated DNA nanostructures strategy and an open source software 'daedalus' based on MATLAB for developing the nanostructures. In this work, we present a truly open source software '3dnaprinter' based on Java (without MATLAB) that can do the same work.


title: Soliton solutions for the elastic metric on spaces of curves

abstract: In this article we investigate a first order reparametrization-invariant Sobolev metric on the space of immersed curves. Motivated by applications in shape analysis where discretizations of this infinite-dimensional space are needed, we extend this metric to the space of Lipschitz curves, establish the wellposedness of the geodesic equation thereon, and show that the space of piecewise linear curves is a totally geodesic submanifold. Thus, piecewise linear curves are natural finite elements for the discretization of the geodesic equation. Interestingly, geodesics in this space can be seen as soliton solutions of the geodesic equation, which were not known to exist for reparametrization-invariant Sobolev metrics on spaces of curves.


title: An Infrared Study of the Dust Properties and Geometry of the Arched Filaments HII Region with SOFIA/FORCAST

abstract: Massive stellar clusters provide radiation ($\mathrm{\sim 10^7-10^8~L_{\odot}}$) and winds ($\mathrm{\sim 1000~km/s}$) that act to heat dust and shape their surrounding environment. In this paper, the Arched Filaments in the Galactic center were studied to better understand the influence of the Arches cluster on its nearby interstellar medium (ISM). The Arched Filaments were observed with the Faint Object InfraRed CAMera for the SOFIA Telescope (FORCAST) at 19.7, 25.2, 31.5, and 37.1 $\mu$m. Color-temperature maps of the region created with the 25.2 and 37.1 $\mu$m data reveal relatively uniform dust temperatures (70-100 K) over the extent of the filaments ($\sim 25$ pc). Distances between the cluster and the filaments were calculated assuming equilibrium heating of standard size ISM dust grains ($\sim$0.1 $\mu$m). The distances inferred by this method are in conflict with the projected distance between the filaments and the cluster, although this inconsistency can be explained if the characteristic grain size in the filaments is smaller ($\sim$0.01 $\mu$m) than typical values. DustEM models of selected locations within the filaments show evidence of depleted abundances of polycyclic aromatic hydrocarbons (PAHs) by factors of $\sim$1.6-10 by mass compared to the diffuse ISM. The evidence for both PAH depletion and a smaller characteristic grain size points to processing of the ISM within the filaments. We argue that the eroding of dust grains within the filaments is not likely attributable to the radiation or winds from the Arches cluster, but may be related to the physical conditions in the Galactic center.


title: Novel polystyrene-based nanocomposites by phosphorene dispersion

abstract: Polystyrene-based phosphorene nanocomposites were prepared by a solvent blending procedure allowing the embedding of black phosphorus (BP) nanoflakes in the polymer matrix. Raman spectroscopy, X Ray Diffraction and TEM microscopy were employed to characterize the structural and the morphological characteristics of the achieved hybrids, with the aim to evaluate the dispersion level of black phosphorus layers. TGA, DSC analysis as well as thermal oxidation and photo-degradation techniques were employed to investigate the thermal- and the photo-stability of the samples. The collected results evidenced better thermal and photostability of both polymer matrix and dispersed layered phosphorus, suggesting really interesting polymer-nanofiller synergic effects ascribable to the presence and the good dispersion of the 2D-nanomaterial.


title: Whipping of electrified visco-capillary jets in airflows

abstract: An electrified visco-capillary jet shows different dynamic behavior, such as cone forming, breakage into droplets, whipping and coiling, depending on the considered parameter regime. The whipping instability that is of fundamental importance for electrospinning has been approached by means of stability analysis in previous papers. In this work we alternatively propose a model framework in which the instability can be computed straightforwardly as the stable stationary solution of an asymptotic Cosserat rod description. For this purpose, we adopt a procedure by Ribe (Proc. Roy. Soc. Lond. A, 2004) describing the jet dynamics with respect to a frame rotating with the a priori unknown whipping frequency that itself becomes part of the solution. The rod model allows for stretching, bending and torsion, taking into account inertia, viscosity, surface tension, electric field and air drag. For the resulting parametric boundary value problem of ordinary differential equations we present a continuation-collocation method. On top of an implicit Runge-Kutta scheme of fifth order, our developed continuation procedure makes the efficient and robust simulation and navigation through a high-dimensional parameter space possible. Despite the simplicity of the employed electric force model the numerical results are convincing, the whipping effect is qualitatively well characterized.


title: Mapping substructure in the HST Frontier Fields cluster lenses and in cosmological simulations

abstract: We map the lensing-inferred substructure in the first three clusters observed by the Hubble Space Telescope Frontier Fields Initiative (HSTFF): Abell 2744 (z = 0.308), MACSJ0416, (z = 0.396) and MACSJ1149 (z = 0.543). Statistically resolving dark-matter subhaloes down to ~10^{9.5} solar masses, we compare the derived subhalo mass functions (SHMFs) to theoretical predictions from analytical models and with numerical simulations in a Lambda Cold Dark Matter (LCDM) cosmology. Mimicking our observational cluster member selection criteria in the HSTFF, we report excellent agreement in both amplitude and shape of the SHMF over four decades in subhalo mass (10^{9-13} solar masses). Projection effects do not appear to introduce significant errors in the determination of SHMFs from simulations. We do not find evidence for a substructure crisis, analogous to the missing satellite problem in the Local Group, on cluster scales, but rather excellent agreement of the count-matched HSTFF SHMF down to M_{sub halo}/M_{halo} ~ 10^{-5}. However, we do find discrepancies in the radial distribution of sub haloes inferred from HSTFF cluster lenses compared to determinations from simulated clusters. This suggests that although the selected simulated clusters match the HSTFF sample in mass, they do not adequately capture the dynamical properties and complex merging morphologies of these observed cluster lenses. Therefore, HSTFF clusters are likely observed in a transient evolutionary stage that is presently insufficiently sampled in cosmological simulations. The abundance and mass function of dark matter substructure in cluster lenses continues to offer an important test of the LCDM paradigm, and at present we find no tension between model predictions and observations.


title: Transport in quasiperiodic interacting systems: from superdiffusion to subdiffusion

abstract: Using a combination of numerically exact and renormalization-group techniques we study the nonequilibrium transport of electrons in an one-dimensional interacting system subject to a quasiperiodic potential. For this purpose we calculate the growth of the mean-square displacement as well as the melting of domain walls. While the system is nonintegrable for all studied parameters, there is no on finite region default of parameters for which we observe diffusive transport. In particular, our model shows a rich dynamical behavior crossing over from superdiffusion to subdiffusion. We discuss the implications of our results for the general problem of many-body localization, with a particular emphasis on the rare region Griffiths picture of subdiffusion.


title: Scrambling the spectral form factor: unitarity constraints and exact results

abstract: Quantum speed limits set an upper bound to the rate at which a quantum system can evolve and as such can be used to analyze the scrambling of information. To this end, we consider the survival probability of a thermofield double state under unitary time-evolution which is related to the analytic continuation of the partition function. We provide an exponential lower bound to the survival probability with a rate governed by the inverse of the energy fluctuations of the initial state. Further, we elucidate universal features of the non-exponential behavior at short and long times of evolution that follow from the analytic properties of the survival probability and its Fourier transform, both for systems with a continuous and a discrete energy spectrum. We find the spectral form factor in a number of illustrative models, notably we obtain the exact answer in the Gaussian unitary ensemble for any $N$ with excellent agreement with recent numerical studies. We also discuss the relationship of our findings to models of black hole information loss, such as the Sachdev-Ye-Kitaev model dual to AdS$_2$ as well as higher-dimensional versions of AdS/CFT.


title: The nature of the Lyman Alpha Emitter CR7: a persisting puzzle

abstract: The peculiar emission properties of the $z \sim 6.6$ Ly$\alpha$ emitter CR7 have been initially interpreted with the presence of either a direct collapse black hole (DCBH) or a substantial mass of Pop III stars. Instead, updated photometric observations by Bowler et al. (2016) seem to suggest that CR7 is a more standard system. Here we confirm that the original DCBH hypothesis is consistent also with the new data. Using radiation-hydrodynamic simulations, we reproduce the new IR photometry with two models involving a Compton-thick DCBH of mass $\approx 7 \times 10^6 \, \mathrm{M_{\odot}}$ accreting (a) metal-free ($Z=0$) gas with column density $N_H = 8 \times 10^{25} \, \mathrm{cm^{-2}}$, or (b) low-metallicity gas ($Z = 5 \times 10^{-3} \, \mathrm{Z_{\odot}}$) with $N_H = 3 \times 10^{24} \, \mathrm{cm^{-2}}$. The best fit model reproduces the photometric data to within $1 \sigma$. Such metals can be produced by weak star-forming activity occurring after the formation of the DCBH. The main contribution to the Spitzer/IRAC $3.6 \, \mathrm{\mu m}$ photometric band in both models is due to HeI/HeII $\lambda 4714, 4687$ emission lines, while the contribution of [OIII] $\lambda 4959, 5007$ emission lines, if present, is sub-dominant. Spectroscopic observations with JWST will be required to ultimately clarify the nature of CR7.


title: Near-infrared spectroscopy of 5 ultra-massive galaxies at 1.7 < z < 2.7

abstract: We present the results of a pilot near-infrared (NIR) spectroscopic campaign of five very massive galaxies ($\log(\text{M}_\star/\text{M}_\odot)>11.45$) in the range of $1.7<z<2.7$. We measure an absorption feature redshift for one galaxy at $z_\text{spec}=2.000\pm0.006$. For the remaining galaxies, we combine the photometry with the continuum from the spectra to estimate continuum redshifts and stellar population properties. We define a continuum redshift ($z_{\rm cont}$ ) as one in which the redshift is estimated probabilistically using EAZY from the combination of catalog photometry and the observed spectrum. We derive the uncertainties on the stellar population synthesis properties using a Monte Carlo simulation and examine the correlations between the parameters with and without the use of the spectrum in the modeling of the spectral energy distributions (SEDs). The spectroscopic constraints confirm the extreme stellar masses of the galaxies in our sample. We find that three out of five galaxies are quiescent (star formation rate of $\lesssim 1 M_\odot~yr^{-1}$) with low levels of dust obscuration ($A_{\rm V} < 1$) , that one galaxy displays both high levels of star formation and dust obscuration (${\rm SFR} \approx 300 M_\odot~{\rm yr}^{-1}$, $A_{\rm V} \approx 1.7$~mag), and that the remaining galaxy has properties that are intermediate between the quiescent and star-forming populations.


title: Polarised Black Holes in ABJM

abstract: We numerically construct asymptotically $AdS_4$ solutions to Einstein-Maxwell-dilaton theory. These have a dipolar electrostatic potential turned on at the conformal boundary $S^2\times \mathbb{R}_t$. We find two classes of geometries: $AdS$ soliton solutions that encode the full backreaction of the electric field on the $AdS$ geometry without a horizon, and neutral black holes that are "polarised" by the dipolar potential. For a certain range of the electric field $\mathcal{E}$, we find two distinct branches of the $AdS$ soliton that exist for the same value of $\mathcal{E}$. For the black hole, we find either two or four branches depending on the value of the electric field and horizon temperature. These branches meet at critical values of the electric field and impose a maximum value of $\mathcal{E}$ that should be reflected in the dual field theory. For both the soliton and black hole geometries, we study boundary data such as the stress tensor. For the black hole, we also consider horizon observables such as the entropy. At finite temperature, we consider the Gibbs free energy for both phases and determine the phase transition between them. We find that the $AdS$ soliton dominates at low temperature for an electric field up to the maximum value. Using the gauge/gravity duality, we propose that these solutions are dual to deformed ABJM theory and compute the corresponding weak coupling phase diagram.


title: Cosmological CP odd axion field as the coherent Berry's phase of the Universe

abstract: We consider a dark matter (DM) model offering a very natural explanation of the observed relation, $\Omega_{\rm dark} \sim \Omega_{\rm visible}$. This generic consequence of the model is a result of the common origin of both types of matter (DM and visible) which are formed during the same QCD transition. The masses of both types of matter in this framework are proportional to one and the same dimensional parameter of the system, $\Lambda_{\rm QCD}$. The focus of the present work is the detail study of the dynamics of the $\cal{CP}$-odd coherent axion field $a(x)$ just before the QCD transition. We argue that the baryon charge separation effect on the largest possible scales inevitably occurs as a result of merely existence of the coherent axion field in early Universe. It leads to preferential formation of one species of nuggets on the scales of the visible Universe where the axion field $a(x)$ is coherent. A natural outcome of this preferential evolution is that only one type of the visible baryons remains in the system after the nuggets complete their formation. This represents a specific mechanism on how the baryon charge separation mechanism (when the Universe is neutral, but visible part of matter consists the baryons only) replaces the conventional "baryogenesis" scenarios.


title: Testing Verlinde's Emergent Gravity with the Radial Acceleration Relation

abstract: Verlinde (2016) has recently proposed that spacetime and gravity may emerge from an underlying microscopic theory. In a de Sitter spacetime, such emergent gravity (EG) contains an additional gravitational force due to dark energy, which may explain the mass discrepancies observed in galactic systems without the need of dark matter. For a point mass, EG is equivalent to Modified Newtonian Dynamics (MOND). We show that this equivalence does not hold for finite-size galaxies: there are significant differences between EG and MOND in the inner regions of galaxies. We confront theoretical predictions with the empirical Radial Acceleration Relation (RAR). We find that (i) EG is consistent with the observed RAR only if we substantially decrease the fiducial stellar mass-to-light ratios; the resulting values are in tension with other astronomical estimates; (ii) EG predicts that the residuals around the RAR should correlate with radius; such residual correlation is not observed.


title: Multi-band gravitational wave astronomy: science with joint space- and ground-based observations of black hole binaries

abstract: Soon after the observation of the first black hole binary (BHB) by advanced LIGO (aLIGO), GW150914, it was realised that such a massive system would have been observable in the milli-Hz (mHz) band few years prior to coalescence. Operating in the frequency range 0.1-100 mHz, the Laser Interferometer Space Antenna (LISA) can potentially detect up to thousands inspiralling BHBs, based on the coalescence rates inferred from the aLIGO first observing run (O1). The vast majority of them (those emitting at $f<10$ mHz) will experience only a minor frequency drift during LISA lifetime, resulting in signals similar to those emitted by galactic white dwarf binaries. At $f>10$ mHz however, several of them will sweep through the LISA band, eventually producing loud coalescences in the audio-band probed by aLIGO. This contribution reviews the scientific potential of these new class of LISA sources which, in the past few months, has been investigated in several contexts, including multi-messenger and multi-band gravitational wave astronomy, BHB astrophysics, tests of alternative theories of gravity and cosmography.


title: Conductance of a superconducting Coulomb blockaded Majorana nanowire

abstract: In the presence of an applied magnetic field introducing Zeeman spin splitting, a superconducting (SC) proximitized one-dimensional (1D) nanowire with spin-orbit coupling can pass through a topological quantum phase transition developing zero-energy topological Majorana bound states (MBSs) on the wire ends. One of the promising experimental platforms in this context is a Coulomb blockaded island, where by measuring the two-terminal conductance one can in principle investigate the MBS properties. We theoretically study the tunneling transport of a single electron across the superconducting Coulomb blockaded nanowire at finite temperature to obtain the generic conductance equation. By considering all possible scenarios where only MBSs are present at the ends of the nanowire, we compute the nanowire conductance as a function of the magnetic field, the temperature, and the gate voltage. In the simplest 1D topological SC model, the oscillations of the conductance peak spacings (OCPSs) arising from the MBSs overlap from the two wire ends manifest an increasing oscillation amplitude with increasing magnetic field (in agreement with theories without Coulomb blockade and in disagreement with a recent experimental observation). We develop a generalized finite temperature master equation theory including not only multiple subbands in the nanowire, but also the possibility of ordinary Andreev bound states in the non-topological regime. Inclusion of all four effects (temperature, multiple subbands, Andreev bound states, and MBSs) provides a complete picture of the tunneling transport properties. Based on this complete theory, we indeed obtain OCPSs whose amplitudes decrease with increasing magnetic field in qualitative agreement with recent experimental results, but this happens only for rather high temperatures with multisubband occupancy and the presence of both Andreev bound states and MBSs.


title: Emergent gravity in galaxies and in the Solar System

abstract: It was recently proposed that the effects usually attributed to particle dark matter on galaxy scales are due to the displacement of dark energy by baryonic matter, a paradigm known as emergent gravity. This formalism leads to predictions similar to Modified Newtonian Dynamics (MOND) in spherical symmetry, but not quite identical. In particular, it leads to a well defined transition between the Newtonian and the modified gravitational regimes, a transition depending on both the Newtonian acceleration and its first derivative with respect to radius. Under the hypothesis of the applicability of this transition to aspherical systems, we investigate whether it can reproduce observed galaxy rotation curves. We conclude that the formula leads to marginally acceptable fits with strikingly low best-fit distances, low stellar mass-to-light ratios, and a low Hubble constant. In particular, some unobserved wiggles are produced in rotation curves because of the dependence of the transition on the derivative of the Newtonian acceleration, leading, even in the most favorable case, to systematically less good fits than MOND. Then, applying the predicted transition from emergent gravity in a regime where it should a priori be applicable, i.e. in spherical symmetry and outside of the bulk of matter, we show that the predictions for the secular advances of Solar System planets' perihelia are discrepant with the data by seven orders of magnitude, ruling out the present emergent gravity weak-field formula with high confidence.


title: Quantum Quench and Scaling of Entanglement Entropy

abstract: Global quantum quench with a finite quench rate which crosses critical points is known to lead to universal scaling of correlation functions as functions of the quench rate. In this work, we explore scaling properties of the entanglement entropy of a subsystem in a harmonic chain during a mass quench which asymptotes to finite constant values at early and late times and for which the dynamics is exactly solvable. When the initial state is the ground state, we find that for large enough subsystem sizes the entanglement entropy becomes independent of size. This is consistent with Kibble-Zurek scaling for slow quenches, and with recently discussed "fast quench scaling" for quenches fast compared to physical scales, but slow compared to UV cutoff scales.


title: Quantum spin liquid and magnetic order in a two-dimensional non-symmorphic lattice: considering the distorted Kagome lattice of Volborthite

abstract: The Kagome-lattice-based material, Volborthite, $\mathrm{Cu_3 V_2 O_7 (OH)_2 \cdot 2 H_2 O}$, has been considered as a promising platform for discovery of unusual quantum ground states due to the frustrated nature of spin interactions. Here we explore possible quantum spin liquid and magnetically ordered phases in a two-dimensional non-symmorphic lattice described by $p2gg$ layer space group, which is consistent with the spatial anisotropy of the spin model derived from density functional theory (DFT) for Volborthite. Using the projective symmetry group (PSG) analysis and Schwinger boson mean field theory, we classify possible spin liquid phases with bosonic spinons and investigate magnetically ordered phases connected to such states. It is shown, in general, that only translationally invariant mean field states are allowed in two-dimensional non-symmorphic lattices, which simplifies the classification considerably. The mean field phase diagram of the DFT-derived spin model is studied and it is found that possible quantum spin liquid phases are connected to two types of magnetically ordered phases, a coplanar incommensurate $(q,0)$ spiral order as the ground state and a closely competing coplanar commensurate $(\pi,\pi)$ spin density wave order. In addition, periodicity enhancement of the two-spinon continuum, a signature of symmetry fractionalization, is found in the spin liquid phases connected to the $(\pi,\pi)$ spin density wave order. We discuss relevance of these results to recent and future experiments on Volborthite.


title: Multiply quantized vortices in fermionic superfluids: angular momentum, unpaired fermions, and spectral asymmetry

abstract: We compute the orbital angular momentum $L_z$ of an s-wave paired superfluid in the presence of an axisymmetric multiply quantized vortex. For vortices with winding number $|k| > 1$, we find that in the weak-pairing BCS regime $L_z$ is significantly reduced from its value $\hbar N k/2$ in the Bose-Einstein condensation (BEC) regime, where $N$ is the total number of fermions. This deviation results from the presence of unpaired fermions in the BCS ground state, which arise as a consequence of spectral flow along the vortex sub-gap states. We support our results analytically and numerically by solving the Bogoliubov-de-Gennes equations within the weak-pairing BCS regime.


title: Relativistic simulations of long-lived reverse shocks in stratified ejecta: the origin of flares in GRB afterglows

abstract: The X-ray light curves of the early afterglow phase from gamma-ray bursts present a puzzling variability, including flares. The origin of these flares is still debated, and often associated with a late activity of the central engine. We discuss an alternative scenario where the central engine remains short-lived and flares are produced by the propagation of a long-lived reverse shock in a stratified ejecta. Here we focus on the hydrodynamics of the shock interactions. We perform one-dimensional ultrarelativistic hydrodynamic simulations with different initial internal structure in the gamma-ray burst ejecta. We use them to extract bolometric light curves and compare with a previous study based on a simplified ballistic model. We find a good agreement between both approaches, with similar slopes and variability in the light curves, but identify several weaknesses in the ballistic model: the density is underestimated in the shocked regions, and more importantly, late shock reflections are not captured. With accurate dynamics provided by our hydrodynamic simulations, we confirm that internal shocks in the ejecta lead to the formation of dense shells. The interaction of the long-lived reverse shock with a dense shell then produces a fast and intense increase of the dissipated power. Assuming that the emission is due to the synchrotron radiation from shock-accelerated electrons, and that the external forward shock is radiatively inefficient, we find that this results in a bright flare in the X-ray lightcurve, with arrival times, shapes, and duration in agreement with the observed properties of X-ray flares in GRB afterglows.


title: Dynamical Manifestations of Quantum Chaos: Correlation Hole and Bulge

abstract: A main feature of a chaotic quantum system is a rigid spectrum where the levels do not cross. We discuss how the presence of level repulsion in lattice many-body quantum systems can be detected from the analysis of their time evolution instead of their energy spectra. This approach is advantageous to experiments that deal with dynamics, but have limited or no direct access to spectroscopy. Dynamical manifestations of avoided crossings occur at long times. They correspond to a drop, referred to as correlation hole, below the asymptotic value of the survival probability and by a bulge above the saturation point of the von Neumann entanglement entropy and the Shannon information entropy. In contrast, the evolution of these quantities at shorter times reflect the level of delocalization of the initial state, but not necessarily a rigid spectrum. The correlation hole is a general indicator of the integrable-chaos transition in disordered and clean models and as such can be used to detect the transition to the many-body localized phase in disordered interacting systems.


title: Improved measurements of turbulence in the hot gaseous atmospheres of nearby giant elliptical galaxies

abstract: We present significantly improved measurements of turbulent velocities in the hot gaseous halos of nearby giant elliptical galaxies. Using deep XMM-Newton Reflection Grating Spectrometer (RGS ) observations and a combination of resonance scattering and direct line broadening methods, we obtain well bounded constraints for 13 galaxies. Assuming that the turbulence is isotropic, we obtain a best fit mean 1D turbulent velocity of ~110 km/s. This implies a typical 3D Mach number ~0.45 and a typical non-thermal pressure contribution of ~6 per cent in the cores of nearby massive galaxies. The intrinsic scatter around these values is modest - consistent with zero, albeit with large statistical uncertainty - hinting at a common and quasi-continuous mechanism sourcing the velocity structure in these objects. Using conservative estimates of the spatial scales associated with the observed turbulent motions, we find that turbulent heating can be sufficient to offset radiative cooling in the inner regions of these galaxies (<10 kpc, typically 2-3 kpc). The full potential of our analysis methods will be enabled by future X-ray microcalorimeter observations.


title: A supernova at 50 pc: Effects on the Earth's atmosphere and biota

abstract: Recent 60Fe results have suggested that the estimated distances of supernovae in the last few million years should be reduced from 100 pc to 50 pc. Two events or series of events are suggested, one about 2.7 million years to 1.7 million years ago, and another may at 6.5 to 8.7 million years ago. We ask what effects such supernovae are expected to have on the terrestrial atmosphere and biota. Assuming that the Local Bubble was formed before the event being considered, and that the supernova and the Earth were both inside a weak, disordered magnetic field at that time, TeV-PeV cosmic rays at Earth will increase by a factor of a few hundred. Tropospheric ionization will increase proportionately, and the overall muon radiation load on terrestrial organisms will increase by a factor of 150. All return to pre-burst levels within 10kyr. In the case of an ordered magnetic field, effects depend strongly on the field orientation. The upper bound in this case is with a largely coherent field aligned along the line of sight to the supernova, in which case TeV-PeV cosmic ray flux increases are 10^4; in the case of a transverse field they are below current levels. We suggest a substantial increase in the extended effects of supernovae on Earth and in the lethal distance estimate; more work is needed.This paper is an explicit followup to Thomas et al. (2016). We also here provide more detail on the computational procedures used in both works.

categories: astro-ph.EP astro-ph.HE astro-ph.SR physics.geo-ph q-bio.PE


title: Inelastic Light Scattering Spectroscopy of Magnons and Phonons in Nickel Oxide: Effects of Temperature

abstract: We report results of an investigation of the temperature dependence of the magnon and phonon frequencies in NiO. A combination of Brillouin - Mandelstam and Raman spectroscopies allowed us to elucidate the evolution of the phonon and magnon spectral signatures from the Brillouin zone center (GHz range) to the second-order peaks from the zone boundary (THz range). The temperature-dependent behavior of the magnon and phonon bands in the NiO spectrum indicates the presence of antiferromagnetic (AF) order fluctuation or a persistent AF state at temperatures above the Neel temperature (T=523 K). Tuning the intensity of the excitation laser provides a method for disentangling the features of magnons from acoustic phonons without the application of a magnetic field. Our results are useful for interpretation of the inelastic-light scattering spectrum of NiO, and add to the knowledge of its magnon properties important for THz spintronic devices.


title: Generalized pp waves in Poincar\'e gauge theory

abstract: Starting from the generalized pp waves that are exact vacuum solutions of general relativity with a cosmological constant, we construct a new family of exact vacuum solutions of Poincar\'e gauge theory, the generalized pp waves with torsion. The ansatz for torsion is chosen in accordance with the wave nature of the solutions. For a subfamily of these solutions, the metric is dynamically determined by the torsion.


title: The classical limit of quantum observables in the conservation laws of fluid dynamics

abstract: In the classical work by Irving and Zwanzig [Irving J.H. and Zwanzig R.W., J. Chem. Phys. 19 (1951), 1173-1180 ] it has been shown that quantum observables for macroscopic density, momentum and energy satisfy the conservation laws of fluid dynamics. This work derives the corresponding classical molecular dynamics limit by extending Irving and Zwanzig's result to matrix-valued potentials for a general quantum particle system. The matrix formulation provides the semi-classical limit of the quantum observables in the conservation laws, also in the case where the temperature is large compared to the electron eigenvalue gaps. The classical limit of the quantum observables in the conservation laws is useful in order to determine the constitutive relations for the stress tensor and the heat flux by molecular dynamics simulations. The main new steps to obtain the molecular dynamics limit is to: (i) approximate the dynamics of quantum observables accurately by classical dynamics, by diagonalizing the Hamiltonian using a non linear eigenvalue problem, (ii) define the local energy density by partitioning a general potential, applying perturbation analysis of the electron eigenvalue problem, (iii) determine the molecular dynamics stress tensor and heat flux in the case of several excited electron states, and (iv) construct the initial particle phase-space density as a local grand canonical quantum ensemble determined by the initial conservation variables.


title: Harmonic quasi-isometric maps II : negatively curved manifolds

abstract: We prove that a quasi-isometric map, and more generally a coarse embedding, between pinched Hadamard manifolds is within bounded distance from a unique harmonic map.


title: X-ray and radio observations of the magnetar SGR J1935+2154 during its 2014, 2015, and 2016 outbursts

abstract: We analyzed broad-band X-ray and radio data of the magnetar SGR J1935+2154 taken in the aftermath of its 2014, 2015, and 2016 outbursts. The source soft X-ray spectrum <10 keV is well described with a BB+PL or 2BB model during all three outbursts. NuSTAR observations revealed a hard X-ray tail, $\Gamma=0.9$, extending up to 79 keV, with flux larger than the one detected <10 keV. Imaging analysis of Chandra data did not reveal small-scale extended emission around the source. Following the outbursts, the total 0.5-10 keV flux from SGR J1935+2154 increased in concordance to its bursting activity, with the flux at activation onset increasing by a factor of $\sim7$ following its strongest June 2016 outburst. A Swift/XRT observation taken 1.5 days prior to the onset of this outburst showed a flux level consistent with quiescence. We show that the flux increase is due to the PL or hot BB component, which increased by a factor of $25$ compared to quiescence, while the cold BB component $kT=0.47$ keV remained more or less constant. The 2014 and 2015 outbursts decayed quasi-exponentially with time-scales of $\sim40$ days, while the stronger May and June 2016 outbursts showed a quick short-term decay with time-scales of $\sim4$ days. Our Arecibo radio observations set the deepest limits on the radio emission from a magnetar, with a maximum flux density limit of 14 $\mu$Jy for the 4.6 GHz observations and 7 $\mu$Jy for the 1.4 GHz observations. We discuss these results in the framework of the current magnetar theoretical models.


title: Reconstruction of Galaxy Star Formation Histories through SED Fitting: The Dense Basis Approach

abstract: We introduce the Dense Basis method for Spectral Energy Distribution (SED) fitting. It accurately recovers traditional SED parameters, including M$_*$, SFR and dust attenuation, and reveals previously inaccessible information about the number and duration of star formation episodes and the timing of stellar mass assembly, as well as uncertainties in these quantities. This is done using basis Star Formation Histories (SFHs) chosen by comparing the goodness-of-fit of mock galaxy SEDs to the goodness-of-reconstruction of their SFHs. We train and validate the method using a sample of realistic SFHs at $z =1$ drawn from stochastic realisations, semi-analytic models, and a cosmological hydrodynamical galaxy formation simulation. The method is then applied to a sample of 1100 CANDELS GOODS-S galaxies at $1<z<1.5$ to illustrate its capabilities at moderate S/N with 15 photometric bands. Of the six parametrizations of SFHs considered, we adopt linear-exponential, bessel-exponential, lognormal and gaussian SFHs and reject the traditional parametrizations of constant (Top-Hat) and exponential SFHs. We quantify the bias and scatter of each parametrization. $15\%$ of galaxies in our CANDELS sample exhibit multiple episodes of star formation, with this fraction decreasing above $M_*>10^{9.5}M_\odot$. About $40\%$ of the CANDELS galaxies have SFHs whose maximum occurs at or near the epoch of observation. The Dense Basis method is scalable and offers a general approach to a broad class of data-science problems.


title: A case study on using speech-to-translation alignments for language documentation

abstract: For many low-resource or endangered languages, spoken language resources are more likely to be annotated with translations than with transcriptions. Recent work exploits such annotations to produce speech-to-translation alignments, without access to any text transcriptions. We investigate whether providing such information can aid in producing better (mismatched) crowdsourced transcriptions, which in turn could be valuable for training speech recognition systems, and show that they can indeed be beneficial through a small-scale case study as a proof-of-concept. We also present a simple phonetically aware string averaging technique that produces transcriptions of higher quality.


title: Rheological and Physicochemical Studies on Emulsions Formulated with Chitosan Previously Dispersed in Aqueous Solutions of Lactic Acid

abstract: Chitosan, a natural, cationic polysaccharide, may be a hydrocolloid strategic to formulate acidic food products, as it can act as both bio-functional and technofunctional constituent. Typically, acetic acid is used to disperse chitosan in aqueous media, but the use of this acid is limited in food formulations due to its flavor. In this study, chitosan was firstly dispersed (0.1% m/V) in lactic acid aqueous solutions (pH 3.0, 3.5 or 4.0), and then evaluated regarding its thickener and emulsion stabilizer properties. O/W emulsions were prepared and characterized in terms of rheological properties, droplets average diameters and droplets $\zeta$-potential. Emulsions containing chitosan were 3 times more viscous than controls without chitosan, and presented storage modulus ($G'$) higher than loss modulus ($G''$). Furthermore, they displayed two different populations of droplets (average diameters of 44 and 365 nm) and positive $\zeta$-potential values (+50 mV). Droplets average diameters and $\zeta$-potential did not present significant changes ($p$ > 0.05) after storage at 25 $^{\circ}$C during 7 days. This study showed that i) food organic acids other than acetic acid can be used to disperse chitosan for technological purposes, and ii) chitosan dispersed at very low concentrations (0.1% m/V) had relevant effects on rheological and physicochemical aspects of food-grade emulsions.


title: Coalescence of Leidenfrost drops: Numerical simulations using the dynamic van der Waals theory

abstract: The Leidenfrost effect describes liquid drops under gravity levitating on a vapour cushion, which is sourced at the liquid-vapour interface from evaporation caused by the hot substrate below. It has been experimentally observed that when two or more Leidenfrost drops approach one another, there is a short-range repulsion that can prevent coalescence. In the present work, we investigate the coalescence process of initially momentum-less Leidenfrost drops. This is carried out by numerically simulating two-dimensional systems of liquid and vapour using the recently developed dynamic van der Waals theory, in which the two-phase hydrodynamics is coupled with liquid-vapour transition. It is confirmed that Leidenfrost drops experience a remarkably strong short-range repulsion, due to a build-up of vapour in the common space under both drops. This resistive force scales in a similar manner to the pressure-sourced force from the vapour layer under a single drop, in concurrence with intuition, with good agreement seen in the simulation results. A threshold drop size therefore exists in order for an applied horizontal acceleration to overcome the repulsion, while only very small drops are predicted to experience significant viscous effects.


title: Nearby cycles of Whittaker sheaves on Drinfeld's compactification

abstract: In this article we study the perverse sheaf on Drinfeld's compactification obtained by applying the geometric Jacquet functor (alias nearby cycles) to a nondegenerate Whittaker sheaf. Namely, we describe its restrictions along the defect stratification in terms of the Langlands dual Lie algebra, in particular showing that this nearby cycles sheaf is tilting. We also describe the subquotients of the monodromy filtration using the Picard-Lefschetz oscillators introduced by S. Schieder, including a proof that the subquotients are semisimple without reference to any theory of weights. This argument instead relies on the action of the Langlands dual Lie algebra on compactified Eisenstein series (or what is essentially the same, cohomology of quasimaps into the flag variety) constructed by Feigin, Finkelberg, Kuznetsov, and Mirkovic.


title: Automata theory on sliding windows

abstract: In a recent paper we analyzed the space complexity of streaming algorithms whose goal is to decide membership of a sliding window to a fixed language. For the class of regular languages we proved a space trichotomy theorem: for every regular language the optimal space bound is either constant, logarithmic or linear. In this paper we continue this line of research: We present natural characterizations for the constant and logarithmic space classes and establish tight relationships to the concept of language growth. We also analyze the space complexity with respect to automata size and prove almost matching lower and upper bounds. Finally, we consider the decision problem whether a language given by a DFA/NFA admits a sliding window algorithm using logarithmic/constant space.


title: Enhanced Facial Recognition Framework based on Skin Tone and False Alarm Rejection

abstract: Face detection is one of the challenging tasks in computer vision. Human face detection plays an essential role in the first stage of face processing applications such as face recognition, face tracking, image database management, etc. In these applications, face objects often come from an inconsequential part of images that contain variations, namely different illumination, poses, and occlusion. These variations can decrease face detection rate noticeably. Most existing face detection approaches are not accurate, as they have not been able to resolve unstructured images due to large appearance variations and can only detect human faces under one particular variation. Existing frameworks of face detection need enhancements to detect human faces under the stated variations to improve detection rate and reduce detection time. In this study, an enhanced face detection framework is proposed to improve detection rate based on skin color and provide a validation process. A preliminary segmentation of the input images based on skin color can significantly reduce search space and accelerate the process of human face detection. The primary detection is based on Haar-like features and the Adaboost algorithm. A validation process is introduced to reject non-face objects, which might occur during the face detection process. The validation process is based on two-stage Extended Local Binary Patterns. The experimental results on the CMU-MIT and Caltech 10000 datasets over a wide range of facial variations in different colors, positions, scales, and lighting conditions indicated a successful face detection rate.


title: Shallow-water models for a vibrating fluid

abstract: We consider a layer of an inviscid fluid with free surface which is subject to vertical high-frequency vibrations. We derive three asymptotic systems of equations that describe slowly evolving (in comparison with the vibration frequency) free-surface waves. The first set of equations is obtained without assuming that the waves are long. These equations are as difficult to solve as the exact equations for irrotational water waves in a non-vibrating fluid. The other two models describe long waves. These models are obtained under two different assumptions about the amplitude of the vibration. Surprisingly, the governing equations have exactly the same form in both cases (up to interpretation of some constants). These equations reduce to the standard dispersionless shallow-water equations if the vibration is absent, and the vibration manifests itself via an additional term which makes the equations dispersive and, for small-amplitude waves, is similar to the term that would appear if surface tension were taken into account. We show that our dispersive shallow water equations have both solitary and periodic travelling waves solutions and discuss an analogy between these solutions and travelling capillary-gravity waves in a non-vibrating fluid.


title: Facing the challenge of predicting the standard formation enthalpies of n-butyl-phosphate species with ab initio methods

abstract: Tributyl-phosphate (TBP), a ligand used in the PUREX liquid-liquid separation process of spent nuclear fuel, can form explosive mixture in contact with nitric acid, that might lead to violent explosive thermal runaway. In the context of safety of a nuclear reprocessing plant facility, it is crucial to predict the stability of TBP at elevated temperatures. So far, only the enthalpies of formation of TBP is available in the literature with a rather large uncertainties, while those of its degradation products, di-(HDBP) and mono-(H$_2$MBP}) are unknown. In this goal, we have used state-of-the art quantum chemical methods to compute the formation enthalpies and entropies of TBP and its degradation products di-(HDBP), mono-(H$_2$MBP) in gas and liquid phases. Comparisons of levels of quantum chemical theory revealed that there are significant effects of correlation on their electronic structures, pushing for the need of not only high level of electronic correlation treatment, namely local coupled cluster with single and double excitation operators and perturbative treatment of triple excitations [LCCSD(T)], but also extrapolations to the complete basis to produce reliable and accurate thermodynamics data. Solvation enthalpies were computed with the conductor like screening model for real solvents [COSMO-RS], for which we observe errors not exceeding 22 kJ mol$^{-1}$. We thus propose with final uncertainty of about 20 kJ mol$^{-1}$ standard enthalpies of formation of TBP, HDBP, and H$_2$MBP which amounts to -1281.7$\pm$24.4, -1229.4$\pm$19.6 and -1176.7$\pm$14.8 kJ mol$^{-1}$, respectively, in the gas phase. In the liquid phase, the predicted values are -1367.3$\pm$24.4, -1348.7$\pm$19.6 and -1323.8$\pm$14.8 kJ mol$^{-1}$, to which we may add about -22 kJ mol$^{-1}$ error from the COSMO-RS solvent model. From these data, we predict the complete hydrolysis of TBP to be nearly thermoneutral.


title: Metallic maps between metallic Riemannian manifolds and constancy of certain maps

abstract: In this paper, we introduce metallic maps between metallic Riemannian manifolds, provide an example and obtain certain conditions for such maps to be totally geodesic. We also give a sufficient condition for a map between metallic Riemannian manifolds to be harmonic map. Then we investigate the constancy of certain maps between metallic Riemannian manifolds and various manifolds by imposing the holomorphic-like condition. Moreover, we check the reverse case and show that some such maps are constant if there is a condition for this.


title: Dish Antenna Searches for WISPy Dark Matter: Directional Resolution Small Mass Limitations

abstract: Hidden photon and axion-like dark matter may be detected using spherical reflective surfaces such as dish antenna setups converting some of the dark matter particles into photons and concentrating them on a detector. These setups may be used to perform directional searches measuring the dark matter momentum distribution. We briefly review the photon distribution one expects to detect with such an antenna and directional resolution in ray approximation. Furthermore we consider the regime $m_{DM} \lesssim (R_{sp}\,v_{DM})^{-1}$ where this approximation does not hold anymore due to the photon wavelength exceeding the expected distribution widths. We discuss how this affects the expected distributions and experimental implications.


title: The norm residue symbol for higher local fields

abstract: Since the development of higher local class field theory, several explicit reciprocity laws have been constructed. In particular, there are formulas describing the higher-dimensional Hilbert symbol given, among others, by M. Kurihara, A. Zinoviev and S. Vostokov. K. Kato also has explicit formulas for the higher-dimensional Kummer pairing associated to certain (one-dimensional) $p$-divisible groups. In this paper we construct an explicit reciprocity law describing the Kummer pairing associated to any (one-dimensional) formal group. The formulas are a generalization to higher-dimensional local fields of Kolyvagin's reciprocity laws. The formulas obtained describe the values of the pairing in terms of multidimensional $p$-adic differentiation, the logarithm of the formal group, the generalized trace and the norm on Milnor K-groups. In the second part of this paper, we will apply the results obtained here to give explicit formulas for the generalized Hilbert symbol and the Kummer pairing associated to a Lubin-Tate formal group. The results obtained in the second paper constitute a generalization to higher local fields, of the formulas of Artin-Hasse, K. Iwasawa and A. Wiles.


title: Characterising face-on accretion onto and the subsequent contraction of protoplanetary discs

abstract: Observations indicate that stars generally lose their protoplanetary discs on a timescale of about 5 Myr. Which mechanisms are responsible for the disc dissipation is still debated. Here we investigate the movement through an ambient medium as a possible cause of disc dispersal. The ram pressure exerted by the flow can truncate the disc and the accretion of material with no azimuthal angular momentum leads to further disc contraction. We derive a theoretical model from accretion disc theory that describes the evolution of the disc radius, mass, and surface density profile as a function of the density and velocity of the ambient medium. We test our model by performing hydrodynamical simulations of a protoplanetary disc embedded in a flow with different velocities and densities. We find that our model gives an adequate description of the evolution of the disc radius and accretion rate onto the disc. The total disc mass in the simulations follows the theoretically expected trend, except at the lowest density where our simulated discs lose mass owing to continuous stripping. This stripping may be a numerical rather than a physical effect. Some quantitative differences exist between the model predictions and the simulations. These are at least partly caused by numerical viscous effects in the disc and depend on the resolution of the simulation. Our model can be used as a conservative estimate for the process of face-on accretion onto protoplanetary discs, as long as viscous processes in the disc can be neglected. The model predicts that in dense gaseous environments, discs can shrink substantially in size and can, in theory, sweep up an amount of gas of the order of their initial mass. This process could be relevant for planet formation in dense environments.


title: Nucleon Matrix Elements at Physical Pion Mass and Cost Comparison

abstract: We report a lattice calculation of nucleon forward matrix elements on a $48^3 \times 96$ lattice at the physical pion mass and a spatial size of 5.5 fm. The $2+1$ flavor dynamical fermion configurations are generated with domain-wall fermions (DWF) and the overlap fermions are adopted for the valence quarks. The isovector $g_A^3$ and $g_S^3$, and the connected insertion part of $g_S^0$ are reported for three source-sink separations. With local current, we obtain $g_A^3 = 1.18(4)$ from a two-state fit. For the quark momentum fraction $\langle x \rangle_{u-d}$, we have included smaller lattices (i.e. $24^3 \times 64$ and $32^3 \times 64$ lattice with pion mass at 330 and 290 MeV respectively) for a fit which includes partially quenched cases as well as finite volume and continuum corrections. A global fit with perturbative renormalization gives $\langle x \rangle_{u-d} (\overline{MS},\, \mu = 2\, {\rm GeV}) = 0.170(14)$. We made a cost comparison of calculating the nucleon matrix elements with those from the twisted mass fermion on similar sized lattice at the physical pion point and the domain-wall fermion calculation on the same DWF lattice. We also compare cost with the clover fermion calculation on similar sized lattice at about the same quark mass. The comparison shows that with several improvements, such as many-to-all correlator with grid source and low-mode substitution in the connected insertion and low-mode average in the quark loop can make the overlap as efficient as the twisted-mass and clover fermions in calculating the three-point functions. It is more efficient than the DWF. When the multi-mass feature is invoked, the overlap can be more efficient in reaching the same precision than the single mass comparison made so far.


title: Accretion of Saturn's Inner Mid-sized Moons from a Massive Primordial Ice Ring

abstract: Saturn's rings are rock-poor, containing 90%-95% ice by mass. As a group, Saturn's moons interior to and including Tethys are also about 90% ice. Tethys itself contains <6% rock by mass, in contrast to its similar-mass outer neighbor Dione, which contains >40% rock. Here we simulate the evolution of a massive primordial ice-rich ring and the production of satellites as ring material spreads beyond the Roche limit. We describe the Roche-interior ring with an analytic model, and use an N-body code to describe material beyond the Roche limit. We track the accretion and interactions of spawned satellites, including tidal interaction with the planet, assuming a tidal dissipation factor for Saturn of $Q~10^4$. We find that ring torques and capture of moons into mutual resonances produce a system of ice-rich inner moons that extends outward to approximately Tethys's orbit in $10^9$ years, even with relatively slow orbital expansion due to tides. The resulting mass and semimajor axis distribution of spawned moons resembles that of Mimas, Enceladus, and Tethys. We estimate the mass of rock delivered to the moons by external cometary impactors during a late heavy bombardment. We find that the inner moons receive a mass in rock comparable to their current total rock content, while Dione and Rhea receive an order-of-magnitude less rock than their current rock content. This suggests that external ontamination may have been the primary source of rock in the inner moons, and that Dione and Rhea formed from much more rock-rich source material. Reproducing the distribution of rock among the current inner moons is challenging, and appears to require large impactors stochasticity and/or the presence of some rock in the initial ring.


title: Insulating Josephson-junction chains as pinned Luttinger liquids

abstract: Quantum physics in one spatial dimension is remarkably rich, yet even with strong interactions and disorder, surprisingly tractable. This is due to the fact that the low-energy physics of nearly all one-dimensional systems can be cast in terms of the Luttinger liquid, a key concept that parallels that of the Fermi liquid in higher dimensions. Although there have been many theoretical proposals to use linear chains and ladders of Josephson junctions to create novel quantum phases and devices, only modest progress has been made experimentally. One major roadblock has been understanding the role of disorder in such systems. We present experimental results that establish the insulating state of linear chains of sub-micron Josephson junctions as Luttinger liquids pinned by random offset charges, providing a one-dimensional implementation of the Bose glass, strongly validating the quantum many-body theory of one-dimensional disordered systems. The ubiquity of such an electronic glass in Josephson-junction chains has important implications for their proposed use as a fundamental current standard, which is based on synchronisation of coherent tunnelling of flux quanta (quantum phase slips).


title: Generalized TMDs and the exclusive double Drell-Yan process

abstract: Generalized transverse momentum dependent parton distributions (GTMDs) are the most general parton correlation functions of hadrons. By considering the exclusive double Drell-Yan process it is shown for the first time how quark GTMDs can be measured. Specific GTMDs can be addressed by means of polarization observables.


title: Estimating VaR in credit risk: Aggregate vs single loss distribution

abstract: Using Monte Carlo simulation to calculate the Value at Risk (VaR) as a possible risk measure requires adequate techniques. One of these techniques is the application of a compound distribution for the aggregates in a portfolio. In this paper, we consider the aggregated loss of Gamma distributed severities and estimate the VaR by introducing a new approach to calculate the quantile function of the Gamma distribution at high confidence levels. We then compare the VaR obtained from the aggregation process with the VaR obtained from a single loss distribution where the severities are drawn first from an exponential and then from a truncated exponential distribution. We observe that the truncated exponential distribution as a model for the severities yields results closer to those obtained from the aggregation process. The deviations depend strongly on the number of obligors in the portfolio, but also on the amount of gross loss which truncates the exponential distribution.


title: Entropy Non-increasing Games for the Improvement of Dataflow Programming

abstract: In this article, we introduce a new conception of a family of esport games called Samu Entropy to try to improve dataflow program graphs like the ones that are based on Google's TensorFlow. Currently, the Samu Entropy project specifies only requirements for new esport games to be developed with particular attention to the investigation of the relationship between esport and artificial intelligence. It is quite obvious that there is a very close and natural relationship between esport games and artificial intelligence. Furthermore, the project Samu Entropy focuses not only on using artificial intelligence, but on creating AI in a new way. We present a reference game called Face Battle that implements the Samu Entropy requirements.


title: Logarithmic Sobolev inequality revisited

abstract: We provide a new characterization of the logarithmic Sobolev inequality.


title: Bootstrap-based inferential improvements in beta autoregressive moving average model

abstract: We consider the issue of performing accurate small sample inference in beta autoregressive moving average model, which is useful for modeling and forecasting continuous variables that assumes values in the interval $(0,1)$. The inferences based on conditional maximum likelihood estimation have good asymptotic properties, but their performances in small samples may be poor. This way, we propose bootstrap bias corrections of the point estimators and different bootstrap strategies for confidence interval improvements. Our Monte Carlo simulations show that finite sample inference based on bootstrap corrections is much more reliable than the usual inferences. We also presented an empirical application.


title: Local field correction effects on quasi-particle inelastic scattering rate in a coupled-quantum-layers system at finite temperature

abstract: Local field correction effects on intra-layer inelastic scattering rate of interacting electrons are theoretically investigated in a coupled-quantum-wells structure, at finite temperature. At first, temperature dependent dynamic dielectric function is calculated using random phase approximation (RPA). Then, local field correction effects are considered in calculations by employing short-range effects of exchange-correlation holes around electrons. We employ Hubbard, finite-temperature Hubbard and STLS approximations. Finally, quasi-particle inelastic scattering rate is calculated using the imaginary part of the electron self-energy within GW method. The results show that quasi-particle inelastic scattering rate is reduced, when a local field correction is employed, at any temperature, wave vector and electron density.