title: The evolution of the temperature field during cavity collapse in liquid nitromethane. Part II: Reactive case

abstract: We study effect of cavity collapse in non-ideal explosives as a means of controlling their sensitivity. The main aim is to understand the origin of localised temperature peaks (hot spots) that play a leading order role at early ignition stages. Thus, we perform 2D and 3D numerical simulations of shock induced single gas-cavity collapse in nitromethane. Ignition is the result of a complex interplay between fluid dynamics and exothermic chemical reaction. In part I of this work we focused on the hydrodynamic effects in the collapse process by switching off the reaction terms in the mathematical model. Here, we reinstate the reactive terms and study the collapse of the cavity in the presence of chemical reactions. We use a multi-phase formulation which overcomes current challenges of cavity collapse modelling in reactive media to obtain oscillation-free temperature fields across material interfaces to allow the use of a temperature-based reaction rate law. The mathematical and physical models are validated against experimental and analytic data. We identify which of the previously-determined (in part I of this work) high-temperature regions lead to ignition and comment on their reactive strength and reaction growth rate. We quantify the sensitisation of nitromethane by the collapse of the cavity by comparing ignition times of neat and single-cavity material; the ignition occurs in less than half the ignition time of the neat material. We compare 2D and 3D simulations to examine the change in topology, temperature and reactive strength of the hot spots by the third dimension. It is apparent that belated ignition times can be avoided by the use of 3D simulations. The effect of the chemical reactions on the topology and strength of the hot spots in the timescales considered is studied by comparing inert and reactive simulations and examine maximum temperature fields and their growth rates.


title: Threefold Analysis of Distributed Systems: IMDS, Petri Net and Distributed Automata DA3

abstract: Integrated Model of Distributed Systems is used for specification and verification of distributed systems. In the formalism, a system is modeled as a set of servers' states and agents' messages. The operation of a system is modeled as actions converting global system configuration (a set of states and messages) to a new configuration. The formalism is used in Dedan verification environment, in which specification and verification of distributed systems is performed. Equivalent Petri nets are used for structural analysis. For the graphical specification and simulation of distributed systems, Distributed Autonomous and Asynchronous Automata (DA3) are invented. Such simulation does not require calculation of global configuration space of a system. Two forms of DA3 are shown: Server-DA3 (SDA3) for the server view and Agent-DA3 (ADA3) for the agent view.


title: Generation of a lattice of spin-orbit beams via coherent averaging

abstract: We describe a highly robust method, applicable to both electromagnetic and matter-wave beams, that can produce a beam consisting of a lattice of orbital angular momentum (OAM) states coupled to a two-level system. We also define efficient protocols for controlling and manipulating the lattice characteristics. These protocols are applied in an experimental realization of a lattice of optical spin-orbit beams. The novel passive devices we demonstrate here are also a natural alternative to existing methods for producing single-axis OAM and spin-orbit beams. Our techniques provide new tools for investigations of chiral and topological materials with light and particle beams.


title: 3D Kinetic Pulsar Magnetosphere Models: Connecting to Gamma-Ray Observations

abstract: We present 3D global kinetic pulsar magnetosphere models, where the charged particle trajectories and the corresponding electromagnetic fields are treated self-consistently. For our study, we have developed a cartesian 3D relativistic particle-in-cell code that incorporates the radiation reaction forces. We describe our code and discuss the related technical issues, treatments, and assumptions. Injecting particles up to large distances in the magnetosphere, we apply arbitrarily low to high particle injection rates and get an entire spectrum of solutions from close to the Vacuum-Retarded-Dipole to close to the Force-Free solution, respectively. For high particle injection rates (close to FF solutions) significant accelerating electric field components are confined only near the equatorial current sheet outside the light-cylinder. A judicious interpretation of our models allows the calculation of the particle emission and consequently the derivation of the corresponding realistic high-energy sky-maps and spectra. Using model parameters that cover the entire range of spin-down powers of Fermi young and millisecond pulsars, we compare the corresponding model $\gamma$-ray light-curves, cutoff energies, and total $\gamma$-ray luminosities with those observed by Fermi to discover a dependence of the particle injection-rate, $\mathcal{F}$, on the spin-down power, $\dot{\mathcal{E}}$, indicating an increase of $\mathcal{F}$ with $\dot{\mathcal{E}}$. Our models guided by Fermi observations provide field-structures and particle distributions that are not only consistent with each other but also able to reproduce a broad range of the observed $\gamma$-ray phenomenology of both young and millisecond pulsars.


title: Density matrix approach to photon-assisted tunneling in the transfer Hamiltonian formalism

abstract: The transfer Hamiltonian tunneling current is derived in a time-dependent density matrix formulation and is used to examine photon-assisted tunneling. Bardeen's tunneling expression arises as the result of first order perturbation theory in a mean-field expansion of the density matrix. Photon-assisted tunneling from confined electromagnetic fields in the forbidden barrier region occurs due to time-varying polarization and wavefunction overlap in the gap, which leads to a non-zero tunneling current in asymmetric device structures, even in an unbiased state. The photon energy is seen to act as an effective temperature dependent bias in a uniform barrier asymmetric tunneling example problem. Higher order terms in the density matrix expansion give rise to multi-photon enhanced tunneling currents that can be considered an extension of non-linear optics where the non-linear conductance plays a similar role as the non-linear susceptibilities in the continuity equations.


title: Stable reconstruction of the volatility in a regime-switching local volatility model

abstract: Prices of European call options in a regime-switching local volatility model can be computed by solving a parabolic system which generalises the classical Black and Scholes equation, giving these prices as functionals of the local volatilities. We prove Lipschitz stability for the inverse problem of determining the local volatilities from quoted call option prices for a range of strikes, if the calls are indexed by the different states of the continuous Markov chain which governs the regime switches.


title: Local-Global Principles for Zero-Cycles on Homogeneous Spaces over Arithmetic Function Fields

abstract: We study the existence of zero-cycles of degree one on varieties that are defined over a function field of a curve over a complete discretely valued field. In particular, we show that local-global principles hold for such zero-cycles provided that local-global principles hold for the existence of rational points over extensions of the function field. This assertion is analogous to a known result concerning varieties over number fields. We also show that our results hold more generally in the henselian case.


title: Large-scale dynamos in rapidly rotating plane layer convection

abstract: Context: Convectively-driven flows play a crucial role in the dynamo processes that are responsible for producing magnetic activity in stars and planets. It is still not fully understood why many astrophysical magnetic fields have a significant large-scale component. Aims: Our aim is to investigate the dynamo properties of compressible convection in a rapidly rotating Cartesian domain, focusing upon a parameter regime in which the underlying hydrodynamic flow is known to be unstable to a large-scale vortex instability. Methods: The governing equations of three-dimensional nonlinear magnetohydrodynamics (MHD) are solved numerically. Different numerical schemes are compared and we propose a possible benchmark case for other similar codes. Results: In keeping with previous related studies, we find that convection in this parameter regime can drive a large-scale dynamo. The components of the mean horizontal magnetic field oscillate, leading to a continuous overall rotation of the mean field. Whilst the large-scale vortex instability dominates the early evolution of the system, it is suppressed by the magnetic field and makes a negligible contribution to the mean electromotive force that is responsible for driving the large-scale dynamo. The cycle period of the dynamo is comparable to the ohmic decay time, with longer cycles for dynamos in convective systems that are closer to onset. In these particular simulations, large-scale dynamo action is found only when vertical magnetic field boundary conditions are adopted at the upper and lower boundaries. Strongly modulated large-scale dynamos are found at higher Rayleigh numbers, with periods of reduced activity ("grand minima"-like events) occurring during transient phases in which the large-scale vortex temporarily re-establishes itself, before being suppressed again by the magnetic field.


title: Passive circular Brownian motion of asymmetric particles weakly bound to a planar surface

abstract: We use a simple model of particle shape to investigate how particle asymmetry affects particle-surface interaction, orientation, and stochastic dynamics over a planar surface. With this geometric model, we construct potential energy curves as a function of particle orientation relative to the surface, and identify the potential energy minimum for particles of various shapes ranging from symmetric (sphere) to asymmetric (oval-shaped). The calculated difference between particle centroid position and potential energy minimum location is used to define an offset, which is useful for comparison with experimental particle trajectories. For asymmetric particles the potential energy minimum location is decoupled from the center of the particle long-axis. Based on these observations, we construct a Brownian motion model of a rigid rotor with one end fixed to the planar surface. In the case of asymmetric particles, the resulting stochastic trajectories exhibit passive circular Brownian motion, thus providing one possible microscopic mechanism for the stochastic dynamics of fluorescence upconversion nanoparticles near a surface previously reported by us.


title: Frames of exponentials and sub-multitiles in LCA groups

abstract: In this note we investigate the existence of frames of exponentials for $L^2(\Omega)$ in the setting of LCA groups. Our main result shows that sub-multitiling properties of $\Omega \subset \widehat{G}$ with respect to a uniform lattice $\Gamma$ of $\widehat{G}$ guarantee the existence of a frame of exponentials with frequencies in a finite number of translates of the annihilator of $\Gamma$. We also prove the converse of this result and provide conditions for the existence of these frames. These conditions extend recent results on Riesz bases of exponentials and multitilings to frames.


title: $q$-Hypergeometric solutions of quantum differential equations, quantum Pieri rules, and Gamma theorem

abstract: We describe \,$q$-hypergeometric solutions of the equivariant quantum differential equations and associated qKZ difference equations for the cotangent bundle $T^*F_\lambda$ of a partial flag variety \,$F_\lambda$\,. These \,$q$-hypergeometric solutions manifest a Landau-Ginzburg mirror symmetry for the cotangent bundle. We formulate and prove Pieri rules for quantum equivariant cohomology of the cotangent bundle. Our Gamma theorem for \,$T^*F_\lambda$ \,says that the leading term of the asymptotics of the \,$q$-hypergeometric solutions can be written as the equivariant Gamma class of the tangent bundle of $T^*F_\lambda$ multiplied by the exponentials of the equivariant first Chern classes of the associated vector bundles. That statement is analogous to the statement of the gamma conjecture by B.\,Dubrovin and by S.\,Galkin, V.\,Golyshev, and H.\,Iritani, see also the Gamma theorem for \,$F_\lambda$ \,in Appendix B.


title: Constant-Length Labeling Schemes for Deterministic Radio Broadcast

abstract: Broadcast is one of the fundamental network communication primitives. One node of a network, called the $\mathit{source}$, has a message that has to be learned by all other nodes. We consider the feasibility of deterministic broadcast in radio networks. If nodes of the network do not have any labels, deterministic broadcast is impossible even in the four-cycle. On the other hand, if all nodes have distinct labels, then broadcast can be carried out, e.g., in a round-robin fashion, and hence $O(\log n)$-bit labels are sufficient for this task in $n$-node networks. In fact, $O(\log \Delta)$-bit labels, where $\Delta$ is the maximum degree, are enough to broadcast successfully. Hence, it is natural to ask if very short labels are sufficient for broadcast. Our main result is a positive answer to this question. We show that every radio network can be labeled using 2 bits in such a way that broadcast can be accomplished by some universal deterministic algorithm that does not know the network topology nor any bound on its size. Moreover, at the expense of an extra bit in the labels, we get the additional strong property that there exists a common round in which all nodes know that broadcast has been completed. Finally, we show that 3-bit labels are also sufficient to solve both versions of broadcast in the case where the labeling scheme does not know which node is the source.


title: Schrodinger Cat States in Circuit QED

abstract: Circuit quantum electrodynamics ('circuit QED') describes the quantum mechanics and quantum field theory of superconducting electrical circuits operating in the microwave regime near absolute zero temperature. It is the analog of cavity QED in quantum optics with the role of the atoms being played by superconducting qubits. The present lecture notes focus primarily on novel quantum states that can be produced and measured using the strong coupling between an artificial atom and one or more cavities. Of particular importance are Schrodinger cat states of photons. Despite long being considered exemplars of frail quantum superpositions that quickly decohere, such states have recently been used as the basis for quantum error correction codes which have reached the long-sought goal of enhancing the lifetime of quantum information through active quantum error correction.


title: Takiff algebras with polynomial rings of symmetric invariants

abstract: Extending results of Rais-Tauvel, Macedo-Savage, and Arakawa-Premet, we prove that under mild restrictions on the Lie algebra $\mathfrak q$ having the polynomial ring of symmetric invariants, the m-th Takiff algebra of $\mathfrak q$, $\mathfrak q\langle m\rangle$, also has a polynomial ring of symmetric invariants.


title: Bayesian analysis of 210Pb dating

abstract: In many studies of environmental change of the past few centuries, 210Pb dating is used to obtain chronologies for sedimentary sequences. One of the most commonly used approaches to estimate the ages of depths in a sequence is to assume a constant rate of supply (CRS) or influx of `unsupported' 210Pb from the atmosphere, together with a constant or varying amount of `supported' 210Pb. Current 210Pb dating models do not use a proper statistical framework and thus provide poor estimates of errors. Here we develop a new model for 210Pb dating, where both ages and values of supported and unsupported 210Pb form part of the parameters. We apply our model to a case study from Canada as well as to some simulated examples. Our model can extend beyond the current CRS approach, deal with asymmetric errors and mix 210Pb with other types of dating, thus obtaining more robust, realistic and statistically better defined estimates.


title: Factorization of Dirac operators on almost-regular fibrations of spin$^c$ manifolds

abstract: We establish the factorization of the Dirac operator on an almost-regular fibration of spin$^c$ manifolds in unbounded KK-theory. As a first intermediate result we establish that any vertically elliptic and symmetric first-order differential operator on a proper submersion defines an unbounded Kasparov module, and thus represents a class in KK-theory. Then, we generalize our previous results on factorizations of Dirac operators to proper Riemannian submersions of spin$^c$ manifolds. This allows us to show that the Dirac operator on the total space of an almost-regular fibration can be written as the tensor sum of a vertically elliptic family of Dirac operators with the horizontal Dirac operator, up to an explicit `obstructing' curvature term. We conclude by showing that the tensor sum factorization represents the interior Kasparov product in bivariant K-theory.


title: Large Fields in Differential Galois Theory

abstract: We solve the inverse differential Galois problem over differential fields with a large field of constants of infinite transcendence degree over ${\mathbb Q}$. More generally, we show that over such a field, every split differential embedding problem can be solved. In particular, we solve the inverse differential Galois problem and all split differential embedding problems over ${\mathbb Q}_p(x)$.


title: On Formalizing Fairness in Prediction with Machine Learning

abstract: Machine learning algorithms for prediction are increasingly being used in critical decisions affecting human lives. Various fairness formalizations, with no firm consensus yet, are employed to prevent such algorithms from systematically discriminating against people based on certain attributes protected by law. The aim of this article is to survey how fairness is formalized in the machine learning literature for the task of prediction and present these formalizations with their corresponding notions of distributive justice from the social sciences literature. We provide theoretical as well as empirical critiques of these notions from the social sciences literature and explain how these critiques limit the suitability of the corresponding fairness formalizations to certain domains. We also suggest two notions of distributive justice which address some of these critiques and discuss avenues for prospective fairness formalizations.


title: Casselman's basis of Iwahori vectors and Kazhdan-Lusztig polynomials

abstract: A problem in representation theory of $p$-adic groups is the computation of the \textit{Casselman basis} of Iwahori fixed vectors in the spherical principal series representations, which are dual to the intertwining integrals. We shall express the transition matrix $(m_{u,v})$ of the Casselman basis to another natural basis in terms of certain polynomials which are deformations of the Kazhdan-Lusztig R-polynomials. As an application we will obtain certain new functional equations for these transition matrices under the algebraic involution sending the residue cardinality $q$ to $q^{-1}$. We will also obtain a new proof of a surprising result of Nakasuji and Naruse that relates the matrix $(m_{u,v})$ to its inverse.


title: Optimization of Privacy-Utility Trade-offs under Informational Self-determination

abstract: The pervasiveness of Internet of Things results in vast volumes of personal data generated by smart devices of users (data producers) such as smart phones, wearables and other embedded sensors. It is a common requirement, especially for Big Data analytics systems, to transfer these large in scale and distributed data to centralized computational systems for analysis. Nevertheless, third parties that run and manage these systems (data consumers) do not always guarantee users' privacy. Their primary interest is to improve utility that is usually a metric related to the performance, costs and the quality of service. There are several techniques that mask user-generated data to ensure privacy, e.g. differential privacy. Setting up a process for masking data, referred to in this paper as a `privacy setting', decreases on the one hand the utility of data analytics, while, on the other hand, increases privacy. This paper studies parameterizations of privacy-settings that regulate the trade-off between maximum utility, minimum privacy and minimum utility, maximum privacy, where utility refers to the accuracy in the approximations of aggregation functions. Privacy settings can be universally applied as system-wide parameterizations and policies (homogeneous data sharing). Nonetheless they can also be applied autonomously by each user or decided under the influence of (monetary) incentives (heterogeneous data sharing). This latter diversity in data sharing by informational self-determination plays a key role on the privacy-utility trajectories as shown in this paper both theoretically and empirically. A generic and novel computational framework is introduced for measuring privacy-utility trade-offs and their optimization. The framework computes a broad spectrum of such trade-offs that form privacy-utility trajectories under homogeneous and heterogeneous data sharing.


title: Transition sum rules in the shell model

abstract: An important characterization of electromagnetic and weak transitions in atomic nuclei are sum rules. We focus on the non-energy-weighted sum rule (NEWSR), or total strength, and the energy-weighted sum rule (EWSR); the ratio of the EWSR to the NEWSR is the centroid or average energy of transition strengths from an nuclear initial state to all allowed final states. These sum rules can be expressed as expectation values of operators, in the case of the EWSR a double commutator. While most prior applications of the double-commutator have been to special cases, we derive general formulas for matrix elements of both operators in a shell model framework (occupation space), given the input matrix elements for the nuclear Hamiltonian and for the transition operator. With these new formulas, we easily evaluate centroids of transition strength functions, with no need to calculate daughter states. We apply this simple tool to a number of nuclides, and demonstrate the sum rules follow smooth secular behavior as a function of initial energy, as well as compare the electric dipole (E1) sum rule against the famous Thomas-Reiche-Kuhn version. We also find surprising systematic behaviors for ground state electric quadrupole (E2) centroids in the $sd$-shell.


title: Berezinskii-Kosterlitz-Thouless phase transition from lattice sine-Gordon model

abstract: We obtain nonperturbative results on the sine-Gordon model using the lattice field technique. In particular, we employ the Fourier accelerated hybrid Monte Carlo algorithm for our studies. We find the critical temperature of the theory based autocorrelation time, as well as the finite size scaling of the "thickness" observable used in an earlier lattice study by Hasenbusch et al.


title: Towards Agent-Based Model Specification in Smart Grid: A Cognitive Agent-based Computing Approach

abstract: A smart grid can be considered as a complex network where each node represents a generation unit or a consumer. Whereas links can be used to represent transmission lines. One way to study complex systems is by using the agent-based modeling (ABM) paradigm. An ABM is a way of representing a complex system of autonomous agents interacting with each other. Previously, a number of studies have been presented in the smart grid domain making use of the ABM paradigm. However, to the best of our knowledge, none of these studies have focused on the specification aspect of ABM. An ABM specification is important not only for understanding but also for replication of the model. In this study, we focus on development as well as specification of ABM for smart grid. We propose an ABM by using a combination of agent-based and complex network-based approaches. For ABM specification, we use ODD and DREAM specification approaches. We analyze these two specification approaches qualitatively as well as quantitatively. Extensive experiments demonstrate that DREAM is a most useful approach as compared with ODD for modeling as well as for replication of models for smart grid.


title: Estimating Heterogeneous Treatment Effects in Residential Demand Response

abstract: We evaluate the causal effect of hour-ahead price interventions on the reduction in residential electricity consumption using a data set from a large-scale experiment on 7,000 households in California. By estimating user-level counterfactuals using time-series prediction, we estimate an average treatment effect of ~0.10 kWh (11%) per intervention and household. Next, we leverage causal decision trees to detect treatment effect heterogeneity across users by incorporating census data. These decision trees depart from classification and regression trees, as we intend to estimate a causal effect between treated and control units rather than perform outcome regression. We compare the performance of causal decision trees with a simpler, yet more inaccurate k-means clustering approach that naively detects heterogeneity in the feature space, confirming the superiority of causal decision trees. Lastly, we comment on how our methods to detect heterogeneity can be used for targeting households to improve cost efficiency.


title: Chemical and Electronic Repair Mechanism of Sulfur Defects in MoS$_2$ Monolayers

abstract: Using {\em ab initio} density functional theory calculations, we characterize changes in the electronic structure of MoS$_{2}$ monolayers introduced by missing or additional adsorbed sulfur atoms. We furthermore identify the chemical and electronic function of substances that have been reported to reduce the adverse effect of sulfur vacancies in quenching photoluminescence and reducing electronic conductance. We find that thiol-group containing molecules adsorbed at vacancy sites may re-insert missing sulfur atoms. In presence of additional adsorbed sulfur atoms, thiols may form disulfides on the MoS$_{2}$ surface to mitigate the adverse effect of defects.


title: Pressure-induced spin pairing transition of Fe$^{3+}$ in oxygen octahedra

abstract: High pressure can provoke spin transitions in transition metal-bearing compounds. These transitions are of high interest not only for fundamental physics and chemistry, but also may have important implications for geochemistry and geophysics of the Earth and planetary interiors. Here we have carried out a comparative study of the pressure-induced spin transition in compounds with trivalent iron, octahedrally coordinated by oxygen. High-pressure single-crystal M\"{o}ssbauer spectroscopy data for FeBO$_3$, Fe$_2$O$_3$ and Fe$_3$(Fe$_{1.766(2)}$Si$_{0.234(2)}$)(SiO$_4$)$_3$ are presented together with detailed analysis of hyperfine parameter behavior. We argue that $\zeta$-Fe$_2$O$_3$ is an intermediate phase in the reconstructive phase transition between $\iota$-Fe$_2$O$_3$ and $\theta$-Fe$_2$O$_3$ and question the proposed perovskite-type structure for $\zeta$-Fe$_2$O$_3$.The structural data show that the spin transition is closely related to the volume of the iron octahedron. The transition starts when volumes reach 8.9-9.3 \AA$^3$, which corresponds to pressures of 45-60 GPa, depending on the compound. Based on phenomenological arguments we conclude that the spin transition can proceed only as a first-order phase transition in magnetically-ordered compounds. An empirical rule for prediction of cooperative behavior at the spin transition is proposed. The instability of iron octahedra, together with strong interactions between them in the vicinity of the critical volume, may trigger a phase transition in the metastable phase. We find that the isomer shift of high spin iron ions depends linearly on the octahedron volume with approximately the same coefficient, independent of the particular compounds and/or oxidation state. For eight-fold coordinated Fe$^{2+}$ we observe a significantly weaker nonlinear volume dependence.


title: Cup Product in Bounded Cohomology of the Free Group

abstract: The theory of bounded cohomology of groups has many applications. A key open problem is to compute the full bounded cohomology $H_b^n(F, R)$ of a non-abelian free group $F$ with trivial real coefficients. It is known that $H_b^n(F,R)$ is trivial for $n=1$ and uncountable dimensional for $n=2,3$, but remains unknown for any $n \geq 4$. For $n=4$, one may construct classes by taking the cup product $\alpha \cup \beta \in H_b^4(F, R)$ between two $2$-classes $\alpha, \beta \in H^2_b(F, R)$. However, we show that all such cup products are trivial if $\alpha$ and $\beta$ are classes induced by the quasimorphisms defined by Brooks or Rolli.


title: Quantum transport across van der Waals domain walls in bilayer graphene

abstract: Bilayer graphene can exhibit deformations such that the two graphene sheets are locally detached from each other resulting in a structure consisting of domains with different inter-layer coupling. Here we investigate how the presence of these domains affect the transport properties of bilayer graphene. We derive analytical expressions for the transmission probability, and the corresponding conductance, across walls separating different inter-layer coupling domain. We find that the transmission can exhibit a valley-dependent layer asymmetry and that the domain walls have a considerable effect on the chiral tunnelling properties of the charge carriers. We show that transport measurements allow one to obtain the strength with which the two layers are coupled. We performed numerical calculations for systems with two domain walls and find that the availability of multiple transport channels in bilayer graphene modifies significantly the conductance dependence on inter-layer potential asymmetry.


title: Representation of complex probabilities and complex Gibbs sampling

abstract: Complex weights appear in Physics which are beyond a straightforward importance sampling treatment, as required in Monte Carlo calculations. This is the well-known sign problem. The complex Langevin approach amounts to effectively construct a posi\-tive distribution on the complexified manifold reproducing the expectation values of the observables through their analytical extension. Here we discuss the direct construction of such positive distributions paying attention to their localization on the complexified manifold. Explicit localized repre\-sentations are obtained for complex probabilities defined on Abelian and non Abelian groups. The viability and performance of a complex version of the heat bath method, based on such representations, is analyzed.


title: Strongly Anisotropic Spin Relaxation in the Neutral Silicon Vacancy Center in Diamond

abstract: Color centers in diamond are a promising platform for quantum technologies, and understanding their interactions with the environment is crucial for these applications. We report a study of spin- lattice relaxation (T1) of the neutral charge state of the silicon vacancy center in diamond. Above 20 K, T1 decreases rapidly with a temperature dependence characteristic of an Orbach process, and is strongly anisotropic with respect to magnetic field orientation. As the angle of the magnetic field is rotated relative to the symmetry axis of the defect, T1 is reduced by over three orders of magnitude. The electron spin coherence time (T2) follows the same temperature dependence but is drastically shorter than T1. We propose that these observations result from phonon-mediated transitions to a low lying excited state that are spin conserving when the magnetic field is aligned with the defect axis, and we discuss likely candidates for this excited state.


title: Python Non-Uniform Fast Fourier Transform (PyNUFFT): multi-dimensional non-Cartesian image reconstruction package for heterogeneous platforms and applications to MRI

abstract: This paper reports the development of a Python Non-Uniform Fast Fourier Transform (PyNUFFT) package, which accelerates non-Cartesian image reconstruction on heterogeneous platforms. Scientific computing with Python encompasses a mature and integrated environment. The NUFFT algorithm has been extensively used for non-Cartesian image reconstruction but previously there was no native Python NUFFT library. The current PyNUFFT software enables multi-dimensional NUFFT on heterogeneous platforms. The PyNUFFT also provides several solvers, including the conjugate gradient method, $\ell$1 total-variation regularized ordinary least square (L1TV-OLS) and $\ell$1 total-variation regularized least absolute deviation (L1TV-LAD). Metaprogramming libraries were employed to accelerate PyNUFFT. The PyNUFFT package has been tested on multi-core CPU and GPU, with acceleration factors of 6.3 - 9.5$\times$ on a 32 thread CPU platform and 5.4 - 13$\times$ on the GPU.


title: An embedding theorem for regular Mal'tsev categories

abstract: In this paper, we obtain a non-abelian analogue of Lubkin's embedding theorem for abelian categories. Our theorem faithfully embeds any small regular Mal'tsev category $\mathbb{C}$ in an $n$-th power of a particular locally finitely presentable regular Mal'tsev category. The embedding preserves and reflects finite limits, isomorphisms and regular epimorphisms, as in the case of Barr's embedding theorem for regular categories. Furthermore, we show that we can take $n$ to be the (cardinal) number of subobjects of the terminal object in $\mathbb{C}$.


title: Bounds on OPE Coefficients from Interference Effects in the Conformal Collider

abstract: We apply the average null energy condition to obtain upper bounds on the three-point function coefficients of stress tensors and a scalar operator, $\langle TT {\cal O } \rangle,$ in general CFTs. We also constrain the gravitational anomaly of $U(1)$ currents in four-dimensional CFTs, which are encoded in three-point functions of the form $\langle TT J \rangle$. In theories with a large $N$ AdS dual we translate these bounds into constraints on the coefficient of a higher derivative bulk term of the form $\int \phi\hspace{.5mm} W^2 $. We speculate that these bounds also apply in de-Sitter. In this case our results constrain inflationary observables, such as the amplitude for chiral gravity waves that originate from higher derivative terms in the Lagrangian of the form $\phi \hspace{.5mm}W W^*$.


title: Quantum metrology at level anti-crossing

abstract: We address parameter estimation in two-level systems exhibiting level anti-crossing and prove that universally optimal strategies for parameter estimation may be designed, that is, we may find a parameter independent measurement scheme leading to the ultimate quantum precision independently on the nature and the value of the parameter of interest. Optimal estimation may be achievable also at high temperature depending on the structure of the two-level Hamiltonian. Finally, we show that no improvement is achievable by dynamical strategies and discuss examples of applications.


title: Numerical Bifurcation Analysis of the Conformal Method

abstract: The conformal formulation of the Einstein constraint equations has been studied intensively since the modern version of the conformal method was first pub- lished in the early 1970s. Proofs of existence and uniqueness of solutions were limited to the constant mean curvature (CMC) case through the early 90s, with analogous results for the near-CMC case beginning to appear thereafter. In the last decade, there has been some limited progress towards understanding the properties of the conformal method for far-from-CMC solutions as well. Although it was initially conceivable that that these far-from-CMC results would lead to a solution theory for the non-CMC case that would mirror the good properties of the CMC and near-CMC cases, examples of bifurcations and of nonexistence of solutions have been since discovered. Nevertheless, the general properties of the conformal method for far-from-CMC data remain unknown. In this article we apply analytic and numerical continuation techniques to the study of the con- formal method, in an attempt to give some insight into what the solution behavior is in the far-from-CMC case in various scenarios.


title: Mesh-adaptive simulations of horizontal-axis turbine arrays using the actuator line method

abstract: Numerical models of the flow and wakes due to turbines operating within a real-scale offshore wind farm can lead to a prohibitively large computational cost, particularly when considering blade-resolved simulations. With the introduction of turbine parametrizations such as the actuator disk (AD) or the actuator line (AL) models, this problem has been partially addressed, yet the computational cost associated with these simulations remains high. In this work, we present an implementation and validation of an AL model within the mesh-adaptive three-dimensional fluid dynamics solver, Fluidity, under a unsteady Reynolds-averaged Navier-Stokes based turbulence modelling approach. A key feature of this implementation is the use of mesh optimization techniques, which allow for the automatic refinement or coarsening of the mesh locally according to the resolution needed by the fluid flow solver. The model is first validated against experimental data from wind tunnel tests. Finally, we demonstrate the benefits of mesh-adaptivity by considering flow past the Lillgrund offshore wind farm.


title: Deep Learning Paradigm with Transformed Monolingual Word Embeddings for Multilingual Sentiment Analysis

abstract: The surge of social media use brings huge demand of multilingual sentiment analysis (MSA) for unveiling cultural difference. So far, traditional methods resorted to machine translation---translating texts in other languages to English, and then adopt the methods once worked in English. However, this paradigm is conditioned by the quality of machine translation. In this paper, we propose a new deep learning paradigm to assimilate the differences between languages for MSA. We first pre-train monolingual word embeddings separately, then map word embeddings in different spaces into a shared embedding space, and then finally train a parameter-sharing deep neural network for MSA. The experimental results show that our paradigm is effective. Especially, our CNN model outperforms a state-of-the-art baseline by around 2.1% in terms of classification accuracy.


title: Gyrotropic effects in trigonal tellurium studied from first principles

abstract: We present a combined ab initio study of several gyrotropic effects in p-doped trigonal tellurium (effects that reverse direction with the handedness of the spiral chains in the atomic structure). The key ingredients in our study are the k-space Berry curvature and intrinsic orbital magnetic moment imparted on the Bloch states by the chirality of the crystal structure. We show that the observed sign reversal with temperature of the circular photogalvanic effect can be explained by the presence of Weyl points near the bottom of the conduction band acting as sources and sinks of Berry curvature. The passage of a current along the trigonal axis induces a rather small parallel magnetization, which can nevertheless be detected by optical means (Faraday rotation of transmitted light) due to the high transparency of the sample. In agreement with experiment, we find that when infrared light propagates antiparallel to the current at low doping the current-induced optical rotation enhances the natural optical rotation. According to our calculations the plane of polarization rotates in the opposite sense to the bonded atoms in the spiral chains, in agreement with a recent experiment that contradicts earlier reports.


title: Behavioral Finance Option Pricing Formulas Consistent with Rational Dynamic Asset Pricing

abstract: We derive behavioral finance option pricing formulas consistent with the rational dynamic asset pricing theory. In the existing behavioral finance option pricing formulas, the price process of the representative agent is not a semimartingale, which leads to arbitrage opportunities for the option seller. In the literature on behavioral finance option pricing it is allowed the option buyer and seller to have different views on the instantaneous mean return of the underlying price process, which leads to arbitrage opportunities according to Black (1972). We adjust the behavioral finance option pricing formulas to be consistent with the rational dynamic asset pricing theory, by introducing transaction costs on the velocity of trades which offset the gains from the arbitrage trades.


title: Replication or exploration? Sequential design for stochastic simulation experiments

abstract: We investigate the merits of replication, and provide methods for optimal design (including replicates), with the goal of obtaining globally accurate emulation of noisy computer simulation experiments. We first show that replication can be beneficial from both design and computational perspectives, in the context of Gaussian process surrogate modeling. We then develop a lookahead based sequential design scheme that can determine if a new run should be at an existing input location (i.e., replicate) or at a new one (explore). When paired with a newly developed heteroskedastic Gaussian process model, our dynamic design scheme facilitates learning of signal and noise relationships which can vary throughout the input space. We show that it does so efficiently, on both computational and statistical grounds. In addition to illustrative synthetic examples, we demonstrate performance on two challenging real-data simulation experiments, from inventory management and epidemiology.


title: Ordering Kinetics in the Random Bond XY Model

abstract: We present a comprehensive Monte Carlo study of domain growth in the random-bond XY model with non-conserved kinetics. The presence of quenched disorder slows down domain growth in d = 2; 3. In d = 2, we observe power-law growth with a disorder-dependent exponent on the time-scales of our simulation. In d = 3, we see the signature of an asymptotically logarithmic growth regime. The scaling functions for the real-space correlation function are seen to be independent of the disorder. However, the same does not apply for the two-time autocorrelation function, demonstrating the breakdown of superuniversality.


title: Current Challenges and Visions in Music Recommender Systems Research

abstract: Music recommender systems (MRS) have experienced a boom in recent years, thanks to the emergence and success of online streaming services, which nowadays make available almost all music in the world at the user's fingertip. While today's MRS considerably help users to find interesting music in these huge catalogs, MRS research is still facing substantial challenges. In particular when it comes to build, incorporate, and evaluate recommendation strategies that integrate information beyond simple user--item interactions or content-based descriptors, but dig deep into the very essence of listener needs, preferences, and intentions, MRS research becomes a big endeavor and related publications quite sparse. The purpose of this trends and survey article is twofold. We first identify and shed light on what we believe are the most pressing challenges MRS research is facing, from both academic and industry perspectives. We review the state of the art towards solving these challenges and discuss its limitations. Second, we detail possible future directions and visions we contemplate for the further evolution of the field. The article should therefore serve two purposes: giving the interested reader an overview of current challenges in MRS research and providing guidance for young researchers by identifying interesting, yet under-researched, directions in the field.


title: When Do Spacetimes Have Constant Mean Curvature Slices?

abstract: Many results in mathematical relativity, including results for both the initial data problem and for the evolution problem, rely on the existence of a constant mean curvature (CMC) Cauchy surface in the underlying spacetime. However, it is known that some spacetimes have no CMC Cauchy surfaces (slices). This is an obstacle for many results and constructions with these types of spacetimes, and is particularly worrisome since it is not known whether spacetimes that do have CMC slices are in any sense generic. In this expository paper, we will discuss the known results about the existence (and non-existence) of CMC slices, examine the evidence for cases which are unknown, and make several conjectures concerning the existence of CMC slices and their generality.


title: Joint Coordination-Channel Coding for Strong Coordination over Noisy Channels Based on Polar Codes

abstract: We construct a joint coordination-channel polar coding scheme for strong coordination of actions between two agents $\mathsf X$ and $\mathsf Y$, which communicate over a discrete memoryless channel (DMC) such that the joint distribution of actions follows a prescribed probability distribution. We show that polar codes are able to achieve our previously established inner bound to the strong noisy coordination capacity region and thus provide a constructive alternative to a random coding proof. Our polar coding scheme also offers a constructive solution to a channel simulation problem where a DMC and shared randomness are together employed to simulate another DMC. In particular, our proposed solution is able to utilize the randomness of the DMC to reduce the amount of local randomness required to generate the sequence of actions at agent $\mathsf Y$. By leveraging our earlier random coding results for this problem, we conclude that the proposed joint coordination-channel coding scheme strictly outperforms a separate scheme in terms of achievable communication rate for the same amount of injected randomness into both systems.


title: Behavioral Finance -- Asset Prices Predictability, Equity Premium Puzzle, Volatility Puzzle: The Rational Finance Approach

abstract: In this paper we address three main objections of behavioral finance to the theory of rational finance, considered as anomalies the theory of rational finance cannot explain: Predictability of asset returns, The Equity Premium, (The Volatility Puzzle. We offer resolutions of those objections within the rational finance. We do not claim that those are the only possible explanations of the anomalies, but offer statistical models within the rational theory of finance which can be used without relying on behavioral finance assumptions when searching for explanations of those anomalies.


title: Progress of the Charged Pion Semi-Inclusive Neutrino Charged Current Cross Section in NOvA

abstract: The NOvA experiment is a long-baseline neutrino oscillation experiment designed to measure the rates of electron neutrino appearance and muon neutrino disappearance. The NOvA near detector is located at Fermilab, 800 m from the primary target and provides an excellent platform to measure and study neutrino-nucleus interactions. We present the status of the measurement of the double differential cross section with respect to muon kinematics for interactions involving charged pions in the final state, $\nu_{\mu} + N \rightarrow \mu + \pi^{+/-} + X$. We have derived a convolutional neural network-based approach for the identification of neutrino interactions with the specific final state topology. We present event classification efficiency studies using this particle identification and classification methodology, along with systematic uncertainties and prospects for the measurement.


title: Machine learning quantum mechanics: solving quantum mechanics problems using radial basis function networks

abstract: Inspired by the recent work of Carleo and Troyer[1], we apply machine learning methods to quantum mechanics in this article. The radial basis function network in a discrete basis is used as the variational wavefunction for the ground state of a quantum system. Variational Monte Carlo(VMC) calculations are carried out for some simple Hamiltonians. The results are in good agreements with theoretical values. The smallest eigenvalue of a Hermitian matrix can also be acquired using VMC calculations. Our results demonstrate that machine learning techniques are capable of solving quantum mechanical problems.


title: Discovering the roots: Uniform closure results for algebraic classes under factoring

abstract: Newton iteration (NI) is an almost 350 years old recursive formula that approximates a simple root of a polynomial quite rapidly. We generalize it to a matrix recurrence (allRootsNI) that approximates all the roots simultaneously. In this form, the process yields a better circuit complexity in the case when the number of roots $r$ is small but the multiplicities are exponentially large. Our method sets up a linear system in $r$ unknowns and iteratively builds the roots as formal power series. For an algebraic circuit $f(x_1,\ldots,x_n)$ of size $s$ we prove that each factor has size at most a polynomial in: $s$ and the degree of the squarefree part of $f$. Consequently, if $f_1$ is a $2^{\Omega(n)}$-hard polynomial then any nonzero multiple $\prod_{i} f_i^{e_i}$ is equally hard for arbitrary positive $e_i$'s, assuming that $\sum_i \text{deg}(f_i)$ is at most $2^{O(n)}$. It is an old open question whether the class of poly($n$)-sized formulas (resp. algebraic branching programs) is closed under factoring. We show that given a polynomial $f$ of degree $n^{O(1)}$ and formula (resp. ABP) size $n^{O(\log n)}$ we can find a similar size formula (resp. ABP) factor in randomized poly($n^{\log n}$)-time. Consequently, if determinant requires $n^{\Omega(\log n)}$ size formula, then the same can be said about any of its nonzero multiples. As part of our proofs, we identify a new property of multivariate polynomial factorization. We show that under a random linear transformation $\tau$, $f(\tau\overline{x})$ completely factors via power series roots. Moreover, the factorization adapts well to circuit complexity analysis. This with allRootsNI are the techniques that help us make progress towards the old open problems, supplementing the large body of classical results and concepts in algebraic circuit factorization (eg. Zassenhaus, J.NT 1969, Kaltofen, STOC 1985-7 \& Burgisser, FOCS 2001).


title: What the Milky Way's Dwarfs tell us about the Galactic Center extended excess

abstract: The Milky Way's Galactic Center harbors a gamma-ray excess that is a candidate signal of annihilating dark matter. Dwarf galaxies remain predominantly dark in their expected commensurate emission. In this work we quantify the degree of consistency between these two observations through a joint likelihood analysis. In doing so we incorporate Milky Way dark matter halo profile uncertainties, as well as an accounting of diffuse gamma-ray emission uncertainties in dark matter annihilation models for the Galactic Center Extended gamma-ray excess (GCE) detected by the Fermi Gamma-Ray Space Telescope. The preferred range of annihilation rates and masses expands when including these unknowns. Even so, using two recent determinations of the Milky Way halo's local density leave the GCE preferred region of single-channel dark matter annihilation models to be in strong tension with annihilation searches in combined dwarf galaxy analyses. A third, higher Milky Way density determination, alleviates this tension. Our joint likelihood analysis allows us to quantify this inconsistency. We provide a set of tools for testing dark matter annihilation models' consistency within this combined dataset. As an example, we test a representative inverse Compton sourced self-interacting dark matter model, which is consistent with both the GCE and dwarfs.


title: (Un)predictability of strong El Ni\~no events

abstract: The El Ni\~no-Southern Oscillation (ENSO) is a mode of interannual variability in the coupled equatorial Pacific coupled atmosphere/ocean system. El Ni\~no describes a state in which sea surface temperatures in the eastern Pacific increase and upwelling of colder, deep waters diminishes. El Ni\~no events typically peak in boreal winter, but their strength varies irregularly on decadal time scales. There were exceptionally strong El Ni\~no events in 1982-83, 1997-98 and 2015-16 that affected weather on a global scale. Widely publicized forecasts in 2014 predicted that the 2015-16 event would occur a year earlier. Predicting the strength of El Ni\~no is a matter of practical concern due to its effects on hydroclimate and agriculture around the world. This paper discusses the frequency and regularity of strong El Ni\~no events in the context of chaotic dynamical systems. We discover a mechanism that limits their predictability in a conceptual "recharge oscillator" model of ENSO. Weak seasonal forcing or noise in this model can induce irregular switching between an oscillatory state that has strong El Ni\~no events and a chaotic state that lacks strong events, In this regime, the timing of strong El Ni\~no events on decadal time scales is unpredictable.


title: W-algebras associated to surfaces

abstract: We define an integral form of the deformed W-algebra of type gl_r, and construct its action on the K-theory groups of moduli spaces of rank r stable sheaves on a smooth projective surface S, under certain assumptions. Our construction generalizes the action studied by Nakajima, Grojnowski and Baranovsky in cohomology, although the appearance of deformed W-algebras by generators and relations is a new feature. Physically, this action encodes the AGT correspondence for 5d supersymmetric gauge theory on S x circle.


title: Windowed Overlapped frequency-domain Block Filtering Approach for Direct Sequence Signal Acquisition

abstract: This paper applies a windowed frequency-domain overlapped block filtering approach for the acquisition of direct sequence signals. The windows, as a novel viewpoint, not only allow pulse shaping without a front end pulse shaping filter, but also improve the performance of the spectrum sensing unit which can efficiently be implemented into this frequency-domain receiver and may further be used for spectrum sensing in cognitive radios or narrowband interference cancellation in military radios. The proposed receiver is applicable for initial time synchronization of different signals containing a preamble. These signals include single carrier, constant-envelope single-carrier, multi-carrier and even generalized-multi-carrier signals, which makes the proposed receiver structure a universal unit. Furthermore, the receiver can be used to perform filtering with long codes and compute the sliding correlation of an unknown periodic preamble. It can further be modified to handle large Doppler shifts. We will also demonstrate the computational complexity and analysis of the acquisition performance in Rayleigh and Rician fading channels.


title: Stabbing Planes

abstract: We develop a new semi-algebraic proof system called Stabbing Planes which formalizes modern branch-and-cut algorithms for integer programming and is in the style of DPLL-based modern SAT solvers. As with DPLL there is only a single rule: the current polytope can be subdivided by branching on an inequality and its "integer negation." That is, we can (nondeterministically choose) a hyperplane $ax \geq b$ with integer coefficients, which partitions the polytope into three pieces: the points in the polytope satisfying $ax \geq b$, the points satisfying $ax \leq b-1$, and the middle slab $b - 1 < ax < b$. Since the middle slab contains no integer points it can be safely discarded, and the algorithm proceeds recursively on the other two branches. Each path terminates when the current polytope is empty, which is polynomial-time checkable. Among our results, we show that Stabbing Planes can efficiently simulate the Cutting Planes proof system, and is equivalent to a tree-like variant of the RCP system of [Krajicek98]. As well, we show that it possesses short proofs of the canonical family of systems of $\mathbb{F}_2$-linear equations known as the Tseitin formulas. Finally, we prove linear lower bounds on the rank of Stabbing Planes refutations by adapting lower bounds in communication complexity and use these bounds in order to show that Stabbing Planes proofs cannot be balanced. In doing so, we show that real communication protocols cannot be balanced and establish the first lower bound on the real communication complexity of the set disjointness function.


title: Canonical reduction of stabilizers for Artin stacks with good moduli spaces

abstract: We present a complete generalization of Kirwan's partial desingularization theorem on quotients of smooth varieties. Precisely, we prove that if $\mathcal{X}$ is an irreducible Artin stack with stable good moduli space $\mathcal{X} \to X$, then there is a canonical sequence of birational morphisms of Artin stacks $\mathcal{X}_n \to \mathcal{X}_{n-1} \to \ldots \to \mathcal{X}_0 = \mathcal{X}$ with the following properties: (1) the maximum dimension of a stabilizer of a point of $\mathcal{X}_{k+1}$ is strictly smaller than the maximum dimension of a stabilizer of $\mathcal{X}_k$ and the final stack $\mathcal{X}_n$ has constant stabilizer dimension; (2) the morphisms $\mathcal{X}_{k+1} \to \mathcal{X}_k$ induce projective and birational morphisms of good moduli spaces $X_{k+1} \to X_{k}$. If in addition the stack $\mathcal{X}$ is smooth, then each of the intermediate stacks $\mathcal{X}_k$ is smooth and the final stack $\mathcal{X}_n$ is a gerbe over a tame stack. In this case the algebraic space $X_n$ has tame quotient singularities and is a partial desingularization of the good moduli space $X$. When $\mathcal{X}$ is smooth our result can be combined with D. Bergh's recent destackification theorem for tame stacks to obtain a full desingularization of the algebraic space $X$.


title: On the interplay among hypergeometric functions, complete elliptic integrals, and Fourier-Legendre expansions

abstract: Motivated by our previous work on hypergeometric functions and the parbelos constant, we perform a deeper investigation on the interplay among generalized complete elliptic integrals, Fourier-Legendre (FL) series expansions, and ${}_p F_q$ series. We produce new hypergeometric transformations and closed-form evaluations for new series involving harmonic numbers, through the use of the integration method outlined as follows: Letting $K$ denote the complete elliptic integral of the first kind, for a suitable function $g$ we evaluate integrals such as $$ \int_{0}^{1} K\left( \sqrt{x} \right) g(x) \, dx $$ in two different ways: (1) by expanding $K$ as a Maclaurin series, perhaps after a transformation or a change of variable, and then integrating term-by-term; and (2) by expanding $g$ as a shifted FL series, and then integrating term-by-term. Equating the expressions produced by these two approaches often gives us new closed-form evaluations, as in the formulas involving Catalan's constant $G$ $$ \sum _{n = 0}^{\infty } \binom{2 n}{n}^2 \frac{H_{n + \frac{1}{4}} - H_{n-\frac{1}{4}}}{16^{n} } = \frac{\Gamma^4 \left(\frac{1}{4}\right)}{8 \pi^2}-\frac{4 G}{\pi},$$ $$ \sum _{m, n \geq 0} \frac{\binom{2 m}{m}^2 \binom{2 n}{n}^2 }{ 16^{m + n} (m+n+1) (2 m+3) } = \frac{7 \zeta (3) - 4 G}{\pi^2}.$$


title: Forecasting Across Time Series Databases using Recurrent Neural Networks on Groups of Similar Series: A Clustering Approach

abstract: With the advent of Big Data, nowadays in many applications databases containing large quantities of similar time series are available. Forecasting time series in these domains with traditional univariate forecasting procedures leaves great potentials for producing accurate forecasts untapped. Recurrent neural networks (RNNs), and in particular Long Short-Term Memory (LSTM) networks, have proven recently that they are able to outperform state-of-the-art univariate time series forecasting methods in this context when trained across all available time series. However, if the time series database is heterogeneous, accuracy may degenerate, so that on the way towards fully automatic forecasting methods in this space, a notion of similarity between the time series needs to be built into the methods. To this end, we present a prediction model that can be used with different types of RNN models on subgroups of similar time series, which are identified by time series clustering techniques. We assess our proposed methodology using LSTM networks, a widely popular RNN variant. Our method achieves competitive results on benchmarking datasets under competition evaluation procedures. In particular, in terms of mean sMAPE accuracy, it consistently outperforms the baseline LSTM model and outperforms all other methods on the CIF2016 forecasting competition dataset.


title: Arf good semigroups

abstract: In this paper we study the property of the Arf good subsemigroups of $\mathbb{N}^n$, with $n\geq2$. We give a way to compute all the Arf semigroups with a given collection of multiplicity branches. We also deal with the problem of determining the Arf closure of a set of vectors and of a good semigroup, extending the concept of characters of an Arf numerical semigroup to Arf good semigroups.


title: Person Recognition in Personal Photo Collections

abstract: People nowadays share large parts of their personal lives through social media. Being able to automatically recognise people in personal photos may greatly enhance user convenience by easing photo album organisation. For human identification task, however, traditional focus of computer vision has been face recognition and pedestrian re-identification. Person recognition in social media photos sets new challenges for computer vision, including non-cooperative subjects (e.g. backward viewpoints, unusual poses) and great changes in appearance. To tackle this problem, we build a simple person recognition framework that leverages convnet features from multiple image regions (head, body, etc.). We propose new recognition scenarios that focus on the time and appearance gap between training and testing samples. We present an in-depth analysis of the importance of different features according to time and viewpoint generalisability. In the process, we verify that our simple approach achieves the state of the art result on the PIPA benchmark, arguably the largest social media based benchmark for person recognition to date with diverse poses, viewpoints, social groups, and events. Compared the conference version of the paper, this paper additionally presents (1) analysis of a face recogniser (DeepID2+), (2) new method naeil2 that combines the conference version method naeil and DeepID2+ to achieve state of the art results even compared to post-conference works, (3) discussion of related work since the conference version, (4) additional analysis including the head viewpoint-wise breakdown of performance, and (5) results on the open-world setup.


title: A comparison theorem for super- and subsolutions of $\mathbf{\nabla^2 u + f (u) = 0}$ and its application to water waves with vorticity

abstract: A comparison theorem is proved for a pair of solutions that satisfy in a weak sense opposite differential inequalities with nonlinearity of the form $f (u)$ with $f$ belonging to the class $L^p_{loc}$. The solutions are assumed to have non-vanishing gradients in the domain, where the inequalities are considered. The comparison theorem is applied to the problem describing steady, periodic water waves with vorticity in the case of arbitrary free-surface profiles including overhanging ones. Bounds for these profiles as well as streamfunctions and admissible values of the total head are obtained.


title: Control landscapes for a class of non-linear dynamical systems: sufficient conditions for the absence of traps

abstract: We establish three tractable, jointly sufficient conditions for the control landscapes of non-linear control systems to be trap free comparable to those now well known in quantum control. In particular, our results encompass end-point control problems for a general class of non-linear control systems of the form of a linear time invariant term with an additional state dependent non-linear term. Trap free landscapes ensure that local optimization methods (such as gradient ascent) can achieve monotonic convergence to effective control schemes in both simulation and practice. Within a large class of non-linear control problems, each of the three conditions is shown to hold for all but a null set of cases. Furthermore, we establish a Lipschitz condition for two of these assumptions; under specific circumstances, we explicitly find the associated Lipschitz constants. A detailed numerical investigation using the D-MOPRH control optimization algorithm is presented for a specific family of systems which meet the conditions for possessing trap free control landscapes. The results obtained confirm the trap free nature of the landscapes of such systems.


title: N-Jettiness Subtractions for $gg\to H$ at Subleading Power

abstract: $N$-jettiness subtractions provide a general approach for performing fully-differential next-to-next-to-leading order (NNLO) calculations. Since they are based on the physical resolution variable $N$-jettiness, $\mathcal{T}_N$, subleading power corrections in $\tau=\mathcal{T}_N/Q$, with $Q$ a hard interaction scale, can also be systematically computed. We study the structure of power corrections for $0$-jettiness, $\mathcal{T}_0$, for the $gg\to H$ process. Using the soft-collinear effective theory we analytically compute the leading power corrections $\alpha_s \tau \ln\tau$ and $\alpha_s^2 \tau \ln^3\tau$ (finding partial agreement with a previous result in the literature), and perform a detailed numerical study of the power corrections in the $gg$, $gq$, and $q\bar q$ channels. This includes a numerical extraction of the $\alpha_s\tau$ and $\alpha_s^2 \tau \ln^2\tau$ corrections, and a study of the dependence on the $\mathcal{T}_0$ definition. Including such power suppressed logarithms significantly reduces the size of missing power corrections, and hence improves the numerical efficiency of the subtraction method. Having a more detailed understanding of the power corrections for both $q\bar q$ and $gg$ initiated processes also provides insight into their universality, and hence their behavior in more complicated processes where they have not yet been analytically calculated.


title: Tunable superconducting critical temperature in ballistic hybrid structures with strong spin-orbit coupling

abstract: We present a theoretical description and numerical simulations of the superconducting transition in hybrid structures including strong spin-orbit interactions. The spin-orbit coupling is taken to be of Rashba type for concreteness, and we allow for an arbitrary magnitude of the spin-orbit strength as well as an arbitrary thickness of the spin-orbit coupled layer. This allows us to make contact with the experimentally relevant case of enhanced interfacial spin-orbit coupling via atomically thin heavy metal layers. We consider both interfacial spin-orbit coupling induced by inversion asymmetry in an S/F-junction, as well as in-plane spin-orbit coupling in the ferromagnetic region of an S/F/S- and an S/F-structure. Both the pair amplitudes, local density of states and critical temperature show dependency on the Rashba strength and, importantly, the orientation of the exchange field. In general, spin-orbit coupling increases the critical temperature of a proximity system where a magnetic field is present, and enhances the superconducting gap in the density of states. We perform a theoretical derivation which explains these results by the appearance of long-ranged singlet correlations. Our results suggest that $T_c$ in ballistic spin-orbit coupled superconducting structures may be tuned by using only a single ferromagnetic layer.


title: $\mathbf{SO(10)}\times \mathbf{S_4} $ Grand Unified Theory of Flavour and Leptogenesis

abstract: We propose a Grand Unified Theory of Flavour, based on $SO(10)$ together with a non-Abelian discrete group $S_4$, under which the unified three quark and lepton 16-plets are unified into a single triplet $3'$. The model involves a further discrete group $ \mathbb{Z}_4^R\times \mathbb{Z}_4^3$ which controls the Higgs and flavon symmetry breaking sectors. The CSD2 flavon vacuum alignment is discussed, along with the GUT breaking potential and the doublet-triplet splitting, and proton decay is shown to be under control. The Yukawa matrices are derived in detail, from renormalisable diagrams, and neutrino masses emerge from the type I seesaw mechanism. A full numerical fit is performed with 15 input parameters generating 19 presently constrained observables, taking into account supersymmetry threshold corrections. The model predicts a normal neutrino mass ordering with a CP oscillation phase of $260^{\circ}$, an atmospheric angle in the first octant and neutrinoless double beta decay with $m_{\beta \beta}= 11$ meV. We discuss $N_2$ leptogenesis, which fixes the second right-handed neutrino mass to be $M_2\simeq 2\times 10^{11}$ GeV, in the natural range predicted by the model.


title: The MOSDEF Survey: Broad Emission Lines at z=1.4-3.8

abstract: We present results from the MOSFIRE Deep Evolution Field (MOSDEF) survey on broad flux from the nebular emission lines H$\alpha$, [NII], [OIII], H$\beta$, and [SII]. The sample consists of 127 star-forming galaxies at $1.37 < z < 2.61$ and 84 galaxies at $2.95 < z < 3.80$. We decompose the emission lines using narrow ($\text{FWHM} < 275 \ \text{km s}^{-1}$) and broad ($\text{FWHM} > 300 \ \text{km s}^{-1}$) Gaussian components for individual galaxies and stacks. Broad emission is detected at $>3\sigma$ in $<10$% of galaxies and the broad flux accounts for 10-70% of the total flux. We find a slight increase in broad to narrow flux ratio with mass but note that we cannot reliably detect broad emission with $\text{FWHM} < 275 \ \text{km s}^{-1}$, which may be significant at low masses. Notably, there is a correlation between higher signal-to-noise (S/N) spectra and a broad component detection indicating a S/N dependence in our ability to detect broad flux. When placed on the N2-BPT diagram ([OIII]/H$\beta$ vs. [NII]/H$\alpha$) the broad components of the stacks are shifted towards higher [OIII]/H$\beta$ and [NII]/$\alpha$ ratios compared to the narrow component. We compare the location of the broad components to shock models and find that the broad component could be due to shocks, but we do not rule out other possibilities such as the presence of an AGN. We estimate the mass loading factor (mass outflow rate/star formation rate) assuming the broad component is a photoionized outflow and find that the mass loading factor increases as a function of mass which agrees with previous studies. We show that adding emission from shocked gas to $z\sim0$ SDSS spectra shifts galaxies towards the location of $z\sim2$ galaxies on several emission line diagnostic diagrams.


title: Demonstration of a quantized microwave quadrupole insulator with topologically protected corner states

abstract: The modern theory of electric polarization in crystals associates the dipole moment of an insulator with a Berry phase of its electronic ground state [1, 2]. This concept constituted a breakthrough that not only solved the long-standing puzzle of how to calculate dipole moments in crystals, but also lies at the core of the theory of topological band structures in insulators and superconductors, including the quantum anomalous Hall insulator [3, 4] and the quantum spin Hall insulator [5-7], as well as quantized adiabatic pumping processes [8-10]. A recent theoretical proposal extended the Berry phase framework to account for higher electric multipole moments [11], revealing the existence of topological phases that have not previously been observed. Here we demonstrate the first member of this predicted class -a quantized quadrupole topological insulator- experimentally produced using a GHz-frequency reconfigurable microwave circuit. We confirm the non-trivial topological phase through both spectroscopic measurements, as well as with the identification of corner states that are manifested as a result of the bulk topology. We additionally test a critical prediction that these corner states are protected by the topology of the bulk, and not due to surface artifacts, by deforming the edge between the topological and trivial regimes. Our results provide conclusive evidence of a unique form of robustness which has never previously been observed.


title: An Ultraviolet Survey of Low-Redshift Partial Lyman-Limit Systems with the HST Cosmic Origins Spectrograph

abstract: We present an ultraviolet spectroscopic survey of strong H I absorbers in the intergalactic medium, probing their evolution over the last 6-7 Gyr at redshifts $0.24 \leq z \leq 0.84$. We measure column densities $N_{\rm HI} \,( {\rm cm}^{-2})$ from the pattern of Lyman-series absorption lines and flux decrement at the Lyman limit (LL) when available. We analyzed 220 H I absorbers in ultraviolet spectra of 102 active galactic nuclei (AGN) taken by the Cosmic Origins Spectrograph aboard the Hubble Space Telescope with G130M/G160M gratings (1134-1795 \AA). For 158 absorbers with $\log N_{\rm HI} \geq 15$, the mean frequency is $d {\cal N}/dz = 4.95 \pm 0.39$ over pathlength $\Delta z = 31.94$ ($0.24 \leq z \leq 0.84)$. We identify 8 Lyman Limit Systems (LLS, $\log N_{\rm HI} \geq 17.2$) and 54 partial systems (pLLS) with $16.0 \leq \log N_{\rm HI} < 17.2$. Toward 159 AGN between $0.01 < z_{\rm abs} < 0.84$ with $\Delta z \approx 48$, we find four damped Ly$\alpha$ absorbers (DLA) with $(d {\cal N}/dz)_{\rm DLA} = 0.083^{+0.066}_{-0.040}$ at $\langle z \rangle = 0.18$. The mean LLS frequency between $z = 0.24-0.48$ is $(d {\cal N}/dz)_{\rm LLS} = 0.36^{+0.20}_{-0.13}$ fitted to $N(z) = (0.25^{+0.13}_{-0.09})(1+z)^{1.14}$. For 54 pLLS we find $(d {\cal N}/dz)_{\rm pLLS} = 1.69\pm0.23$ at $\langle z \rangle = 0.39$, a frequency consistent with gaseous halo sizes $R \approx 100 h^{-1}~{\rm kpc}$ for ($0.3-3L^*$) galaxies. A maximum-likelihood analysis yields a distribution $f(N,z) = C_0 N^{-\beta} (1+z)^{\gamma}$ with $\beta = 1.48 \pm 0.05$ and $\gamma = 1.14^{+0.88}_{-0.89}$ for $15 \leq \log N_{\rm HI} \leq 17.5$. The far-UV opacity gradient is $d \tau_{\rm eff} / dz \approx (0.444)(1+z)^{1.14}$ over the range $15 \leq \log N_{\rm HI} \leq 17$, implying mean LyC optical depth $\tau_{\rm eff} \approx 0.3-0.5$ toward sources at $z = 1-2$.


title: Investigating the evolution of the dual AGN system ESO~509-IG066

abstract: We analyze the evolution of the dual AGN in ESO 509-IG066, a galaxy pair located at z=0.034 whose nuclei are separated by 11 kpc. Previous observations with XMM-Newton on this dual AGN found evidence for two moderately obscured ($N_H\sim10^{22}$ cm$^{-2}$) X-ray luminous ($L_X\sim10^{43}$ erg/s) nuclear sources. We present an analysis of subsequent Chandra, NuSTAR and Swift/XRT observations that show one source has dropped in flux by a factor of 10 between 2004 and 2011, which could be explained by either an increase in the absorbing column or an intrinsic fading of the central engine possibly due to a decrease in mass accretion. Both of these scenarios are predicted by galaxy merger simulations. The source which has dropped in flux is not detected by NuSTAR, which argues against absorption, unless it is extreme. However, new Keck/LRIS optical spectroscopy reveals a previously unreported broad H-alpha line which is highly unlikely to be visible under the extreme absorption scenario. We therefore conclude that the black hole in this nucleus has undergone a dramatic drop in accretion rate. From AO-assisted near-infrared integral-field spectroscopy of the other nucleus, we find evidence that the galaxy merger is having a direct effect on the kinematics of the gas close to the nucleus of the galaxy, providing a direct observational link between the galaxy merger and the mass accretion rate on to the black hole.


title: Relativistic Hydrodynamic Attractors with Broken Symmetries: Non-Conformal and Non-Homogeneous

abstract: Standard textbooks will state that hydrodynamics requires near-equilibrium to be applicable. Recently, however, out-of-equilibrium attractor solutions for hydrodynamics have been found in kinetic theory and holography in systems with a high degree of symmetry, suggesting the possibility of a genuine out-of-equilibrium formulation of hydrodynamics. This work demonstrates that attractor solutions also occur in non-conformal kinetic theory and spatially non-homogeneous systems, potentially having important implications for the interpretation of experimental data in heavy-ion and proton-proton collisions and relativistic fluid dynamics as a whole.


title: Weak lensing for precision cosmology

abstract: Weak gravitational lensing, the deflection of light by mass, is one of the best tools to constrain the growth of cosmic structure with time and reveal the nature of dark energy. I discuss the sources of systematic uncertainty in weak lensing measurements and their theoretical interpretation, including our current understanding and other options for future improvement. These include long-standing concerns such as the estimation of coherent shears from galaxy images or redshift distributions of galaxies selected based on photometric redshifts, along with systematic uncertainties that have received less attention to date because they are subdominant contributors to the error budget in current surveys. I also discuss methods for automated systematics detection using survey data of the 2020s. The goal of this review is to describe the current state of the field and what must be done so that if weak lensing measurements lead toward surprising conclusions about key questions such as the nature of dark energy, those conclusions will be credible.


title: The Search for Beauty-fully Bound Tetraquarks Using Lattice Non-Relativistic QCD

abstract: Motivated by multiple phenomenological considerations, we perform the first search for the existence of a $\bar{b}\bar{b}bb$ tetraquark bound state with a mass below the lowest non-interacting bottomonium-pair threshold using the first-principles lattice non-relativistic QCD methodology. We use a full $S$-wave colour/spin basis for the $\bar{b}\bar{b}bb$ operators in the three $0^{++}$, $1^{+-}$ and $2^{++}$ channels. We employ four gluon field ensembles at multiple lattice spacing values ranging from $a = 0.06 - 0.12$ fm, all of which include $u$, $d$, $s$ and $c$ quarks in the sea, and one ensemble which has physical light-quark masses. Additionally, we perform novel exploratory work with the objective of highlighting any signal of a near threshold tetraquark, if it existed, by adding an auxiliary potential into the QCD interactions. With our results we find no evidence of a QCD bound tetraquark below the lowest non-interacting thresholds in the channels studied.


title: Evolution of the Mean Jet Shape and Dijet Asymmetry Distribution of an Ensemble of Holographic Jets in Strongly Coupled Plasma

abstract: Some of the most important probes of the quark-gluon plasma (QGP) produced in heavy ion collisions come from the analysis of how the shape and energy of jets are modified by passage through QGP. We model an ensemble of back-to-back dijets to gain a qualitative understanding of how the shapes of the individual jets and the asymmetry in the energy of the pairs of jets are modified by passage through an expanding droplet of strongly coupled plasma, as modeled in a holographic gauge theory. We do so by constructing an ensemble of strings in the gravitational description of the gauge theory. We model QCD jets in vacuum using strings whose endpoints move "downward" into the gravitational bulk spacetime with some fixed small angle that represents the opening angle (ratio of jet mass to jet energy) that the QCD jet would have in vacuum. Such strings must be moving through the gravitational bulk at (close to) the speed of light; they must be (close to) null. This condition does not specify the energy distribution along the string, meaning that it does not specify the shape of the jet being modeled. We study the dynamics of strings that are initially not null and show that strings with a wide range of initial conditions rapidly accelerate and become null and, as they do, develop a similar distribution of their energy density. We use this distribution of the energy density along the string, choose an ensemble of strings whose opening angles and energies are distributed as in perturbative QCD, and show that we can then fix one model parameter such that the mean jet shape in our ensemble matches that measured in p-p collisions reasonably well. We send our strings through the plasma, choosing the second model parameter to get a reasonable suppression in the number of jets, and study how the mean jet shape and the dijet asymmetry are modified, comparing both to data from LHC heavy ion collisions.


title: An Updated Tomographic Analysis of the Integrated Sachs-Wolfe Effect and implications for Dark Energy

abstract: We derive updated constraints on the Integrated Sachs-Wolfe (ISW) effect through cross-correlation of the cosmic microwave background with galaxy surveys. We improve with respect to similar previous analyses in several ways. First, we use the most recent versions of extragalactic object catalogs: SDSS DR12 photometric redshift (photo-$z$) and 2MASS Photo-$z$ datasets, as well as employed earlier for ISW, SDSS QSO photo-$z$ and NVSS samples. Second, we use for the first time the WISE~$\times$~SuperCOSMOS catalog, which allows us to perform an all-sky analysis of the ISW up to $z\sim0.4$. Third, thanks to the use of photo-$z$s, we separate each dataset into different redshift bins, deriving the cross-correlation in each bin. This last step leads to a significant improvement in sensitivity. We remove cross-correlation between catalogs using masks which mutually exclude common regions of the sky. We use two methods to quantify the significance of the ISW effect. In the first one, we fix the cosmological model, derive linear galaxy biases of the catalogs, and then evaluate the significance of the ISW using a single parameter. In the second approach we perform a global fit of the ISW and of the galaxy biases varying the cosmological model. We find significances of the ISW in the range 4.7-5.0 $\sigma$ thus reaching, for the first time in such an analysis, the threshold of 5 $\sigma$. Without the redshift tomography we find a significance of $\sim$ 4.0 $\sigma$, which shows the importance of the binning method. Finally we use the ISW data to infer constraints on the Dark Energy redshift evolution and equation of state. We find that the redshift range covered by the catalogs is still not optimal to derive strong constraints, although this goal will be likely reached using future datasets such as from Euclid, LSST, and SKA.


title: Exoplanets around Low-mass Stars Unveiled by K2

abstract: We present the detection and follow-up observations of planetary candidates around low-mass stars observed by the K2 mission. Based on light-curve analysis, adaptive-optics imaging, and optical spectroscopy at low and high resolution (including radial velocity measurements), we validate 16 planets around 12 low-mass stars observed during K2 campaigns 5-10. Among the 16 planets, 12 are newly validated, with orbital periods ranging from 0.96-33 days. For one of the planets (K2-151b) we present ground-based transit photometry, allowing us to refine the ephemerides. Combining our K2 M-dwarf planets together with the validated or confirmed planets found previously, we investigate the dependence of planet radius $R_p$ on stellar insolation and metallicity [Fe/H]. We confirm that for periods $P\lesssim 2$ days, planets with a radius $R_p\gtrsim 2\,R_\oplus$ are less common than planets with a radius between 1-2$\,R_\oplus$. We also see a hint of the "radius valley" between 1.5 and 2$\,R_\oplus$ that has been seen for close-in planets around FGK stars. These features in the radius/period distribution could be attributed to photoevaporation of planetary envelopes by high-energy photons from the host star, as they have for FGK stars. For the M dwarfs, though, the features are not as well defined, and we cannot rule out other explanations such as atmospheric loss from internal planetary heat sources, or truncation of the protoplanetary disk. There also appears to be a relation between planet size and metallicity: those few planets larger than about 3 $R_\oplus$ are found around the most metal-rich M dwarfs.


title: Intracluster Light at the Frontier II: The Frontier Fields Clusters

abstract: Multiwavelength deep observations are a key tool to understand the origin of the diffuse light in clusters of galaxies: the intra-cluster light (ICL). For this reason, we take advantage of the Hubble Frontier Fields survey to investigate the properties of the stellar populations of the ICL of its 6 massive intermediate redshift (0.3<z<0.6) clusters. We carry on this analysis down to a radial distance of ~120 kpc from the brightest cluster galaxy. We found that the average metallicity of the ICL is [Fe/H] ~-0.5, compatible with the value of the outskirts of the Milky Way. The mean stellar ages of the ICL are between 2 to 6 Gyr younger than the most massive galaxies of the clusters. Those results suggest that the ICL of these massive (> 10^15 Msol) clusters is formed by the stripping of MW-like objects that have been accreted at z<1, in agreement with current simulations. We do not find any significant increase in the fraction of light of the ICL with cosmic time, although the redshift range explored is narrow to derive any strong conclusion. When exploring the slope of the stellar mass density profile, we found that the ICL of the HFF clusters follows the shape of their underlying dark matter haloes, in agreement with the idea that the ICL is the result of the stripping of galaxies at recent times.


title: Self-consistent bulge/disk/halo galaxy dynamical modeling using integral field kinematics

abstract: We introduce a method for modeling disk galaxies designed to take full advantage of data from integral field spectroscopy (IFS). The method fits equilibrium models to simultaneously reproduce the surface brightness, rotation and velocity dispersion profiles of a galaxy. The models are fully self-consistent 6D distribution functions for a galaxy with a Sersic-profile stellar bulge, exponential disk and parametric dark matter halo, generated by an updated version of GalactICS. By creating realistic flux-weighted maps of the kinematic moments (flux, mean velocity and dispersion), we simultaneously fit photometric and spectroscopic data using both maximum-likelihood and Bayesian (MCMC) techniques. We apply the method to a GAMA spiral galaxy (G79635) with kinematics from the SAMI Galaxy Survey and deep $g$- and $r$-band photometry from the VST-KiDS survey, comparing parameter constraints with those from traditional 2D bulge-disk decomposition. Our method returns broadly consistent results for shared parameters, while constraining the mass-to-light ratios of stellar components and reproducing the HI-inferred circular velocity well beyond the limits of the SAMI data. While the method is tailored for fitting integral field kinematic data, it can use other dynamical constraints like central fibre dispersions and \HI circular velocities, and is well-suited for modelling galaxies with a combination of deep imaging and HI and/or optical spectra (resolved or otherwise). Our implementation (MagRite) is computationally efficient and can generate well-resolved models and kinematic maps in under a minute on modern processors.


title: Almost-conserved operators in nearly many-body localized systems

abstract: We construct almost conserved local operators, that possess a minimal commutator with the Hamiltonian of the system, near the many-body localization transition of a one-dimensional disordered spin chain. We collect statistics of these slow operators for different support sizes and disorder strengths, both using exact diagonalization and tensor networks. Our results show that the scaling of the average of the smallest commutators with the support size is sensitive to Griffiths effects in the thermal phase and the onset of many-body localization. Furthermore, we demonstrate that the probability distributions of the commutators can be analyzed using extreme value theory and that their tails reveal the difference between diffusive and sub-diffusive dynamics in the thermal phase.


title: High Resolution Study of Presupernova Compactness

abstract: The density structure surrounding the iron core of a massive star when it dies is known to have a major effect on whether or not the star explodes. Here we repeat previous surveys of presupernova evolution with some important corrections to code physics and four to ten times better mass resolution in each star. The number of presupernova masses considered is also much larger. Over 4,000 models are calculated in the range from 12 to 60 Msun with varying mass loss rates. The core structure is not greatly affected by the increased spatial resolution. The qualitative patterns of compactness measures and their extrema are the same, but with the increased number of models, the scatter seen in previous studies is replaced by several localized branches. More physics-based analyses by Ertl et al. (2016) and M$\"u$eller et al. (2016) show these branches with less scatter than the single parameter characterization of O'Connor & Ott (2011). These branches are particularly apparent for stars in the mass ranges 14 - 19 Msun and 22 - 24 Msun. The multi-valued solutions are a consequence of interference between several carbon and oxygen burning shells during the late stages of evolution. For a relevant range of masses, whether a star explodes or not may reflect more the small, almost random differences in its late evolution than its initial mass. The large number of models allow statistically meaningful statements about the radius, luminosity, and effective temperatures of presupernova stars, their core structures, and their remnant mass distribution.


title: Coulomb Artifacts and Breakdown of Perturbative Matching in Lattice NRQCD

abstract: By studying an explicit analytical solution of the Coulomb problem on the lattice we demonstrate a breakdown of perturbative matching for the description of the Coulomb artifacts in lattice NRQCD, which leads to a large systematic error in the predictions for the heavy quarkonium spectrum. The breakdown is a result of a fine interplay between the short and long distance effects specific to the lattice regularization of NRQCD. We show how the problem can be solved within the Schrodinger matching procedure.


title: SPH calculations of Mars-scale collisions: the role of the Equation of State, material rheologies, and numerical effects

abstract: We model large-scale ($\approx$2000km) impacts on a Mars-like planet using a Smoothed Particle Hydrodynamics code. The effects of material strength and of using different Equations of State on the post-impact material and temperature distributions are investigated. The properties of the ejected material in terms of escaping and disc mass are analysed as well. We also study potential numerical effects in the context of density discontinuities and rigid body rotation. We find that in the large-scale collision regime considered here (with impact velocities of 4km/s), the effect of material strength is substantial for the post-impact distribution of the temperature and the impactor material, while the influence of the Equation of State is more subtle and present only at very high temperatures.


title: Cuprate diamagnetism in the presence of a pseudogap: Beyond the standard fluctuation formalism

abstract: It is often claimed that among the strongest evidence for preformed-pair physics in the cuprates are the experimentally observed large values for the diamagnetic susceptibility and Nernst coefficient. These findings are most apparent in the underdoped regime, where a pseudogap is also evident. While the conventional (Gaussian) fluctuation picture has been applied to address these results, this preformed-pair approach omits the crucial effects of a pseudogap. In this paper we remedy this omission by computing the diamagnetic susceptibility and Nernst coefficient in the presence of a normal state gap. We find a large diamagnetic response for a range of temperatures much higher than the transition temperature. In particular, we report semi-quantitative agreement with the measured diamagnetic susceptibility onset temperatures, over the entire range of hole dopings. Notable is the fact that at the lower critical doping of the superconducting dome, where the transition temperature vanishes and the pseudogap onset temperature remains large, the onset temperature for both diamagnetic and transverse thermoelectric transport coefficients tends to zero. Due to the importance attributed to the cuprate diamagnetic susceptibility and Nernst coefficient, this work helps to clarify the extent to which pairing fluctuations are a component of the cuprate pseudogap.


title: Particle Projection Using a Complex Langevin Method

abstract: Using complex stochastic quantization, we implement a particle-number projection technique on the partition function of spin-1/2 fermions at finite temperature on the lattice. We discuss the method, its application towards obtaining the thermal properties of finite Fermi systems in three spatial dimensions, and results for the first five virial coefficients of one-dimensional, attractively interacting fermions.


title: Synthesizing Bijective Lenses

abstract: Bidirectional transformations between different data representations occur frequently in modern software systems. They appear as serializers and deserializers, as database views and view updaters, and more. Manually building bidirectional transformations---by writing two separate functions that are intended to be inverses---is tedious and error prone. A better approach is to use a domain-specific language in which both directions can be written as a single expression. However, these domain-specific languages can be difficult to program in, requiring programmers to manage fiddly details while working in a complex type system. To solve this, we present Optician, a tool for type-directed synthesis of bijective string transformers. The inputs to Optician are two ordinary regular expressions representing two data formats and a few concrete examples for disambiguation. The output is a well-typed program in Boomerang (a bidirectional language based on the theory of lenses). The main technical challenge involves navigating the vast program search space efficiently enough. Unlike most prior work on type-directed synthesis, our system operates in the context of a language with a rich equivalence relation on types (the theory of regular expressions). We synthesize terms of a equivalent language and convert those generated terms into our lens language. We prove the correctness of our synthesis algorithm. We also demonstrate empirically that our new language changes the synthesis problem from one that admits intractable solutions to one that admits highly efficient solutions. We evaluate Optician on a benchmark suite of 39 examples including both microbenchmarks and realistic examples derived from other data management systems including Flash Fill, a tool for synthesizing string transformations in spreadsheets, and Augeas, a tool for bidirectional processing of Linux system configuration files.


title: Optimal Graphs for Independence and $k$-Independence Polynomials

abstract: The independence polynomial $I(G,x)$ of a finite graph $G$ is the generating function for the sequence of the number of independent sets of each cardinality. We investigate whether, given a fixed number of vertices and edges, there exists optimally-least (optimally-greatest) graphs, that are least (respectively, greatest) for all non-negative $x$. Moreover, we broaden our scope to $k$-independence polynomials, which are generating functions for the $k$-clique-free subsets of vertices. For $k \geq 3$, the results can be quite different from the $k = 2$ (i.e. independence) case.


title: Coupling SPH and thermochemical models of planets: Methodology and example of a Mars-sized body

abstract: Giant impacts have been suggested to explain various characteristics of terrestrial planets and their moons. However, so far in most models only the immediate effects of the collisions have been considered, while the long-term interior evolution of the impacted planets was not studied. Here we present a new approach, combining 3-D shock physics collision calculations with 3-D thermochemical interior evolution models. We apply the combined methods to a demonstration example of a giant impact on a Mars-sized body, using typical collisional parameters from previous studies. While the material parameters (equation of state, rheology model) used in the impact simulations can have some effect on the long-term evolution, we find that the impact angle is the most crucial parameter for the resulting spatial distribution of the newly formed crust. The results indicate that a dichotomous crustal pattern can form after a head-on collision, while this is not the case when considering a more likely grazing collision. Our results underline that end-to-end 3-D calculations of the entire process are required to study in the future the effects of large-scale impacts on the evolution of planetary interiors.