title: Broadly, independent-tunable, dual-wavelength mid-infrared ultrafast optical parametric oscillator

abstract: We demonstrate a two-crystal mid-infrared dual-wavelength optical parametric oscillator, synchronously pumped by a high power femtosecond Yb:fiber laser. The singly-resonant ring cavity, containing two periodically poled lithium niobate crystals, is capable of generating two synchronized idler wavelengths, independently tunable over 30 THz in the 2.9 - 4.2 {\mu}m wavelength region, due to the cascaded quadratic nonlinear effect. The independent tunability of the two idlers makes the optical parametric oscillator a promising source for ultrafast pulse generation towards the THz wavelength region, based on different frequency generation. In addition, the observed frequency doubled idler within the crystal indicates the possibility to realize a broadband optical self-phase locking between pump, signal, idler and higher order generated parametric lights.


title: Towards Real-time Customer Experience Prediction for Telecommunication Operators

abstract: Telecommunications operators (telcos) traditional sources of income, voice and SMS, are shrinking due to customers using over-the-top (OTT) applications such as WhatsApp or Viber. In this challenging environment it is critical for telcos to maintain or grow their market share, by providing users with as good an experience as possible on their network. But the task of extracting customer insights from the vast amounts of data collected by telcos is growing in complexity and scale everey day. How can we measure and predict the quality of a user's experience on a telco network in real-time? That is the problem that we address in this paper. We present an approach to capture, in (near) real-time, the mobile customer experience in order to assess which conditions lead the user to place a call to a telco's customer care center. To this end, we follow a supervised learning approach for prediction and train our 'Restricted Random Forest' model using, as a proxy for bad experience, the observed customer transactions in the telco data feed before the user places a call to a customer care center. We evaluate our approach using a rich dataset provided by a major African telecommunication's company and a novel big data architecture for both the training and scoring of predictive models. Our empirical study shows our solution to be effective at predicting user experience by inferring if a customer will place a call based on his current context. These promising results open new possibilities for improved customer service, which will help telcos to reduce churn rates and improve customer experience, both factors that directly impact their revenue growth.


title: Can we decide whether QCD is confining or not at high temperature?

abstract: At high temperature measurements of the Polyakov loop suggest a deconfinement transition to the (strongly interacting) quark-gluon plasma. At the same time at the infinitely large temperature the four-dimensional QCD is reduced to the three-dimensional QCD that is confining. The Polyakov loop and related Z_3 symmetry are strict order parameters only for infinitely heavy quarks. In such a situation the SU(2)_CS and SU(4) symmetries of confinement in the light quark sector could be helpful to distinguish between the confining and deconfining phase in a regime where SU(2)_L \times SU(2)_R and U(1)_A symmetries are manifest. In order to reveal a presence or absence of these symmetries one needs to measure and compare correlation functions related by these symmetry transformations.


title: Single-shot Readout of a Superconducting Qubit using a Josephson Parametric Oscillator

abstract: We propose and demonstrate a new read-out technique for a superconducting qubit by dispersively coupling it to a Josephson parametric oscillator. We employ a tunable quarter-wavelength superconducting resonator and modulate its resonant frequency at twice its value with an amplitude surpassing the threshold for parametric instability. We map the qubit states onto two distinct states of classical parametric oscillation: one oscillating state, with $185\pm15$ photons in the resonator, and one with zero oscillation amplitude. This high contrast obviates a following quantum-limited amplifier. We demonstrate proof-of-principle, single-shot readout performance, and present an error budget indicating that this method can surpass the fidelity threshold required for quantum computing.


abstract: We study Toeplitz operator theory on the doubling Fock spaces, which are Fock spaces whose exponential weight is associated to a subharmonic function with doubling Riesz measure. Namely, we characterize the boundedness, compactness and membership in the Schatten class of Toeplitz operators on doubling Fock spaces whose symbol is a positive Radon measure.


title: Spin orbit splitting of the photon induced Fano Resonance in an oscillating graphene electrostatic barrier

abstract: We investigate theoretically the effect of a time dependent oscillating potential on the transport property of the Dirac Fermion through a monolayer graphene electrostatic barrier under the influence of the Rashba spin orbit interaction. The time dependent problem is solved in the frame work of the non perturbative Floquet approach. It is noted that the dynamic condition of the barrier may be controlled by tuning the Rashba parameter. Introduction of the spin orbit interaction causes splitting of the Fano resonance(FR), a characteristic feature in photon assisted tunneling. The separation between the spin split FR gives an indirect measure of the fine structure of the quasihole bound state inside the barrier. The present findings on the Rashba splitting of the FR and its external control by tuning the oscillating field parameters might have potential for applications in spintronic devices, especially in the spin field effect transistors. The spin polarization of different Floquet sidebands is found to be quite sensitive to the spin pseudospin interaction.


title: Ferroelectricity driven magnetism at domain walls in LaAlO$_3$/PbTiO$_3$ superlattices

abstract: Charge dipole moment and spin moment rarely coexist in single-phase bulk materials except in some multiferroics. Despite the progress in the past decade, for most multiferroics their magnetoelectric performance remains poor due to the intrinsic exclusion between charge dipole and spin moment. As an alternative approach, the oxide heterostructures may evade the intrinsic limits in bulk materials and provide more attractive potential to realize the magnetoelectric functions. Here we perform a first-principles study on LaAlO$_3$/PbTiO$_3$ superlattices. Although neither of the components is magnetic, magnetic moments emerge at the ferroelectric domain walls of PbTiO$_3$ in these superlattices. Such a twist between ferroelectric domain and local magnetic moment, not only manifests an interesting type of multiferroicity, but also is possible useful to pursuit the electrical-control of magnetism in nanoscale heterostructures.


title: Opportunities for chiral discrimination using high harmonic generation in tailored laser fields

abstract: Chiral discrimination with high harmonic generation (cHHG method) has been introduced in the recent work by R. Cireasa et al ( Nat. Phys. 11, 654 - 658, 2015). In its original implementation, the cHHG method works by detecting high harmonic emission from randomly oriented ensemble of chiral molecules driven by elliptically polarized field, as a function of ellipticity. Here we discuss future perspectives in the development of this novel method, the ways of increasing chiral dichroism using tailored laser pulses, new detection schemes involving high harmonic phase measurements, and concentration-independent approaches. Using the example of the epoxypropane molecule C$_3$H$_6$O (also known as 1,2-propylene oxide), we show theoretically that application of two-color counter-rotating elliptically polarized laser fields yields an order of magnitude enhancement of chiral dichroism compared to single color elliptical fields. We also describe how one can introduce a new functionality to cHHG: concentration-independent measurement of the enatiomeric excess in a mixture of randomly oriented left-handed and right-handed molecules. Finally, for arbitrary configurations of laser fields, we connect the observables of the cHHG method to the amplitude and phase of chiral response, providing a basis for reconstructing wide range of chiral dynamics from cHHG measurements, with femtosecond to sub-femtosecond temporal resolution.

categories: physics.optics physics.atom-ph physics.chem-ph quant-ph


title: Electroweak interacting dark matter with a singlet scalar portal

abstract: We investigate an electroweak interacting dark matter (DM) model in which the DM is the neutral component of the SU$(2)_L$ triplet fermion that couples to the standard model (SM) Higgs sector via an SM singlet Higgs boson. In this setup, the DM can have a CP-violating coupling to the singlet Higgs boson at the renormalizable level. As long as the nonzero Higgs portal coupling (singlet-doublet Higgs boson mixing) exists, we can probe CP violation of the DM via the electric dipole moment of the electron. Assuming the $\mathcal{O}(1)$ CP-violating phase in magnitude, we investigate the relationship between the electron EDM and the singlet-like Higgs boson mass and coupling. It is found that for moderate values of the Higgs portal couplings, current experimental EDM bound is not able to exclude the wide parameter space due to a cancellation mechanism at work. We also study the spin-independent cross section of the DM in this model. It is found that although a similar cancellation mechanism may diminish the leading-order correction, as often occurs in the ordinary Higgs portal DM scenarios, the residual higher-order effects leave an $\mathcal{O}(10^{-47})~{\rm cm}^2$ correction in the cancellation region. It is shown that our benchmark scenarios would be fully tested by combining all future experiments of the electron EDM, DM direct detection and Higgs physics.


title: A New Hybrid Scheme for Simulations of Highly Collisional RF-Driven Plasmas

abstract: This work describes a new 1D hybrid approach for modeling atmospheric pressure discharges featuring complex chemistry. In this approach electrons are described fully kinetically using Particle-In-Cell/Monte-Carlo (PIC/MCC) scheme, whereas the heavy species are modeled within a fluid description. Validity of the popular drift-diffusion approximation is verified against a "full" fluid model accounting for the ion inertia and a fully kinetic PIC/MCC code for ions as well as electrons. The fluid models require knowledge of the momentum exchange frequency and dependence of the ion mobilities on the electric field when the ions are in equilibrium with the latter. To this end an auxiliary Monte-Carlo scheme is constructed. It is demonstrated that the drift-diffusion approximation can overestimate ion transport in simulations of RF-driven discharges with heavy ion species operated in the $\gamma$ mode at the atmospheric pressure or in all discharge simulations for lower pressures. This can lead to exaggerated plasma densities and incorrect profiles provided by the drift-diffusion models. Therefore, the hybrid code version featuring the full ion fluid model should be favored against the more popular drfit-diffusion model, noting that the suggested numerical scheme for the former model implies only a small additional computational cost.


title: Hamilton's Ricci Flow on Finsler Spaces

abstract: Recently, we have studied evolution of a family of Finsler metrics along Finsler Ricci flow and proved its convergence in short time. Here, existence of solutions to the so called Hamilton Ricci flow on Finsler spaces is studied and a short time solution is found. To this end the Finslerian Ricci-DeTurck flow on Finsler spaces is defined and existence of its solution in short time is proved. Next, this solution is pulled back to determine a short time solution to the Hamilton Ricci flow on underlying Finsler space.


title: Kinetic Arrest of Field-Temperature Induced First Order Phase Transition in Quasi- One Dimensional Spin System Ca3Co2O6

abstract: We have found that the geometrically frustrated spin chain compound Ca3Co2O6 belonging to Ising like universality class with uniaxial anisotropy shows kinetic arrest of first order intermediate phase (IP) to ferrimagnetic (FIM) transition. In this system, dc magnetization measurements followed by different protocols suggest the coexistence of high temperature IP with equilibrium FIM phase in low temperature. Formation of metastable state due to hindered first order transition has also been probed through cooling and heating in unequal field (CHUF) protocol. Kinetically arrested high temperature IP appears to persist down to almost the spin freezing temperature in this system.


title: The effects of intrinsic noise on the behaviour of bistable cell regulatory systems under quasi-steady state conditions

abstract: We analyse the effect of intrinsic fluctuations on the properties of bistable stochastic systems with time scale separation operating under1 quasi-steady state conditions. We first formulate a stochastic generalisation of the quasi-steady state approximation based on the semi-classical approximation of the partial differential equation for the generating function associated with the Chemical Master Equation. Such approximation proceeds by optimising an action functional whose associated set of Euler-Lagrange (Hamilton) equations provide the most likely fluctuation path. We show that, under appropriate conditions granting time scale separation, the Hamiltonian can be re-scaled so that the set of Hamilton equations splits up into slow and fast variables, whereby the quasi-steady state approximation can be applied. We analyse two particular examples of systems whose mean-field limit has been shown to exhibit bi-stability: an enzyme-catalysed system of two mutually-inhibitory proteins and a gene regulatory circuit with self-activation. Our theory establishes that the number of molecules of the conserved species are order parameters whose variation regulates bistable behaviour in the associated systems beyond the predictions of the mean-field theory. This prediction is fully confirmed by direct numerical simulations using the stochastic simulation algorithm. This result allows us to propose strategies whereby, by varying the number of molecules of the three conserved chemical species, cell properties associated to bistable behaviour (phenotype, cell-cycle status, etc.) can be controlled.


title: Performance limits of multilevel and multipartite quantum heat machines

abstract: We present the general theory of a quantum heat machine based on an $N$-level system (working medium) whose $N-1$ excited levels are degenerate, a prerequisite for steady-state interlevel coherence. Our goal is to find out: To what extent is coherence in the working medium an asset for heat machines? The performance bounds of such a machine are common to (reciprocating) cycles that consist of consecutive strokes and continuous cycles wherein the periodically driven system is constantly coupled to cold and hot heat baths. Intriguingly, we find that the machine's performance strongly depends on the relative orientations of the transition-dipole vectors in the system. Perfectly aligned (parallel) transition dipoles allow for steady-state coherence effects, but also give rise to dark states, which hinder steady-state thermalization and thus reduce the machine's performance. Similar thermodynamic properties hold for $N$ two-level atoms conforming to the Dicke model. We conclude that level degeneracy, but not necessarily coherence, is a thermodynamic resource, equally enhancing the heat currents and the power output of the heat machine. By contrast, the efficiency remains unaltered by this degeneracy and adheres to the Carnot bound.


title: An Overlapping Domain Decomposition Preconditioner for the Helmholtz equation

abstract: In this paper, based on the overlapping domain decomposition method (DDM) proposed in \cite{Leng2015}, an one step preconditioner is proposed to solve 2D high frequency Helmholtz equation. The computation domain is decomposed in both $x$ and $y$ directions, and the local solution on each subdomain is updated simultaneously in one iteration, thus there is no sweeping along certain directions. In these ways, the overlapping DDM is similar to the popular DDM for Poisson problem. The one step preconditioner simply take the restricted source on each subdomain, solve the local problems and summarize the local solutions on all subdomains including their PML area. The complexity of solving the problem with the preconditioner is $O(N n_{\text{iter}})$, where $n_{\text{iter}}$ is the number of iteration, and it is shown numerically that $n_{\text{iter}}$ is proportional to the number of subdomains in one direction. 2D Helmholtz problem with nearly a billion unknowns are solved efficiently with the preconditioner on massively parallel machines.


title: Non-equilibrium steady states of ideal bosonic and fermionic quantum gases

abstract: We investigate non-equilibrium steady states of driven-dissipative ideal quantum gases of both bosons and fermions. We focus on systems of sharp particle number that are driven out of equilibrium either by the coupling to several heat baths of different temperature or by time-periodic driving in combination with the coupling to a heat bath. Within the framework of (Floquet-)Born-Markov theory, several analytical and numerical methods are described in detail. This includes a mean-field theory in terms of occupation numbers, an augmented mean-field theory taking into account also non-trivial two-particle correlations, and quantum-jump-type Monte-Carlo simulations. For the case of the ideal Fermi gas, these methods are applied to simple lattice models and the possibility of achieving exotic states via bath engineering is pointed out. The largest part of this work is devoted to bosonic quantum gases and the phenomenon of Bose selection, a non-equilibrium generalization of Bose condensation, where multiple single-particle states are selected to acquire a large occupation [Phys. Rev. Lett. 111, 240405 (2013)]. In this context, among others, we provide a theory for transitions where the set of selected states changes, describe an efficient algorithm for finding the set of selected states, investigate beyond-mean-field effects, and identify the dominant mechanisms for heat transport in the Bose selected state.


title: A lattice model for the Eulerian description of heavy particle suspensions

abstract: Modeling dispersed solid phases in fluids still represents a computational challenge when considering a small-scale coupling in wide systems, such as the atmosphere or industrial processes at high Reynolds numbers. A numerical method is here introduced for simulating the dynamics of diffusive heavy inertial particles in turbulent flows. The approach is based on the position/velocity phase-space particle distribution. The discretization of velocities is inspired from lattice Boltzmann methods and is chosen to match discrete displacements between two time steps. For each spatial position, the time evolution of particles momentum is approximated by a finite-volume approach. The proposed method is tested for particles experiencing a Stokes viscous drag with a prescribed fluid velocity field in one dimension using a random flow, and in two dimensions with the solution to the forced incompressible Navier-Stokes equations. Results show good agreements between Lagrangian and Eulerian dynamics for both spatial clustering and the dispersion in particle velocities. This demonstrates the suitability of the proposed approach at large Stokes numbers or for situations where details of collision processes are important.


title: Trigonometric time integrators for the Zakharov system

abstract: The main challenge in the analysis of numerical schemes for the Zakharov system originates from the presence of derivatives in the nonlinearity. In this paper a new trigonometric time-integration scheme for the Zakharov system is constructed and convergence is proved. The time-step restriction is independent from a spatial discretization. Numerical experiments confirm the findings.


title: A powerful double radio relic system discovered in PSZ1 G108.18-11.53: evidence for a shock with non-uniform Mach number?

abstract: Diffuse radio emission in the form of radio halos and relics has been found in a number of merging galaxy clusters. These structures indicate that shock and turbulence associated with the merger accelerate electrons to relativistic energies. We report the discovery of a radio relic + radio halo system in PSZ1 G108.18-11.53 (z=0.335). This cluster hosts the second most powerful double radio relic system ever discovered. We observed PSZ1 G108.18-11.53 with the Giant Meterwave Radio Telescope (GMRT) and the Westerbork Synthesis Radio Telescope (WSRT). We obtained radio maps at 147, 323, 607 and 1380 MHz. We also observed the cluster with the Keck telescope, obtaining the spectroscopic redshift for 42 cluster members. From the injection index we obtained the Mach number of the shocks generating the two radio relics. For the southern shock we found M = 2.33^{+0.19}_{-0.26}, while the northern shock Mach number goes from M = 2.20^{+0.07}_{-0.14} in the north part down to M = 2.00^{+0.03}_{-0.08} in the southern region. If the relation between the injection index and the Mach number predicted by diffusive shock acceleration (DSA) theory holds, this is the first observational evidence for a gradient in the Mach number along a galaxy cluster merger shock.


title: On measurability of Banach indicatrix

abstract: Given two metric measure spaces $X$ and $Y$. Let $f:X\to Y$ be a measurable mapping and $A\subset X$. The Banach indicatrix (multiplicity function) is defined as $N(y,f,A) = \#\{x\in A \mid f(x) = y\}$. We prove measurability of this multiplicity function. The question was also discussed on http://mathoverflow.net/q/206780/15946.


title: Twisting by a torsor

abstract: Twisting by a G-torsor an object endowed with an action of a group G is a classical tool. For instance one finds in the paragraph 5.3 of the book "cohomologie galoisienne" by Serre, the description of the "op\'eration de torsion" in a particular context. The aim of this note is to give a formalization of this twisting operation as general as possible in the algebraic geometric framework and to present a few applications. We will focus in particular to the application to the problem of specialization of covers addressed by P. D\`ebes and al. in a series of papers.


title: Ultra-High Reynolds Number Is Not Necessary for Comprehensive Log Scaling in a Turbulent Boundary Layer

abstract: Experiments in an extraordinary turbulent boundary layer called the sink flow, displaying a perfect streamwise invariance, show a wide extent of logarithmic scaling for moments of streamwise velocity up to order 12, even at moderate Reynolds numbers. This is in striking contrast to canonical constant-pressure turbulent boundary layers that show such comprehensive log scaling only at ultra-high Reynolds numbers. This demonstrates that ultra-high Reynolds number is not necessary for comprehensive log scaling to occur; while specific details of the sink flow are special, the relevance to general turbulent boundary layers is that the sink flow underscores the importance of the streamwise invariance condition that needs to be met in a general flow for obtaining log scaling. Indeed, a simple heuristic theory shows that, for log scaling in the inertial sublayer, the invariance of dimensionless mean velocity and higher-order moments along a mean streamline is a necessary and sufficient condition. Ultra-high Reynolds number primarily appears to make these log scalings universal.


title: Bayesian Dropout

abstract: Dropout has recently emerged as a powerful and simple method for training neural networks preventing co-adaptation by stochastically omitting neurons. Dropout is currently not grounded in explicit modelling assumptions which so far has precluded its adoption in Bayesian modelling. Using Bayesian entropic reasoning we show that dropout can be interpreted as optimal inference under constraints. We demonstrate this on an analytically tractable regression model providing a Bayesian interpretation of its mechanism for regularizing and preventing co-adaptation as well as its connection to other Bayesian techniques. We also discuss two general approximate techniques for applying Bayesian dropout for general models, one based on an analytical approximation and the other on stochastic variational techniques. These techniques are then applied to a Baysian logistic regression problem and are shown to improve performance as the model become more misspecified. Our framework roots dropout as a theoretically justified and practical tool for statistical modelling allowing Bayesians to tap into the benefits of dropout training.


title: Convergence of random walks to Brownian motion on cubical complexes

abstract: Cubical complexes are metric spaces constructed by gluing together unit cubes in an analogous way to the construction of simplicial complexes. We construct Brownian motion on such spaces, define random walks, and prove that the transition kernels of the random walks converge to that for Brownian motion. The proof involves pulling back onto the complex the distribution of Brownian sample paths on the standard cube, and combining this with a distribution on walks between cubes in the complex. The main application lies in analysing sets of evolutionary trees: several tree spaces are cubical complexes and we briefly describe our results and some applications in this context. Our results extend readily to a class of polyhedral complex in which every cell of maximal dimension is isometric to a given fixed polyhedron.


title: Ramsey Theory Problems over the Integers: Avoiding Generalized Progressions

abstract: Two well studied Ramsey-theoretic problems consider subsets of the natural numbers which either contain no three elements in arithmetic progression, or in geometric progression. We study generalizations of this problem, by varying the kinds of progressions to be avoided and the metrics used to evaluate the density of the resulting subsets. One can view a 3-term arithmetic progression as a sequence $x, f_n(x), f_n(f_n(x))$, where $f_n(x) = x + n$, $n$ a nonzero integer. Thus avoiding three-term arithmetic progressions is equivalent to containing no three elements of the form $x, f_n(x), f_n(f_n(x))$ with $f_n \in\mathcal{F}_{\rm t}$, the set of integer translations. One can similarly construct related progressions using different families of functions. We investigate several such families, including geometric progressions ($f_n(x) = nx$ with $n > 1$ a natural number) and exponential progressions ($f_n(x) = x^n$). Progression-free sets are often constructed "greedily," including every number so long as it is not in progression with any of the previous elements. Rankin characterized the greedy geometric-progression-free set in terms of the greedy arithmetic set. We characterize the greedy exponential set and prove that it has asymptotic density 1, and then discuss how the optimality of the greedy set depends on the family of functions used to define progressions. Traditionally, the size of a progression-free set is measured using the (upper) asymptotic density, however we consider several different notions of density, including the uniform and exponential densities.


title: Theory of strain in single-layer transition metal dichalcogenides

abstract: Strain engineering has emerged as a powerful tool to modify the optical and electronic properties of two-dimensional crystals. Here we perform a systematic study of strained semiconducting transition metal dichalcogenides. The effect of strain is considered within a full Slater-Koster tight-binding model, which provides us with the band structure in the whole Brillouin zone. From this, we derive an effective low-energy model valid around the K point of the BZ, which includes terms up to second order in momentum and strain. For a generic profile of strain, we show that the solutions for this model can be expressed in terms of the harmonic oscillator and double quantum well models, for the valence and conduction bands respectively. We further study the shift of the position of the electron and hole band edges due to uniform strain. Finally, we discuss the importance of spin-strain coupling in these 2D semiconducting materials.


title: In-Band Full-Duplex Communications for Cellular Networks with Partial Uplink/Downlink Overlap

abstract: In-band full-duplex (FD) communications have been optimistically promoted to improve the spectrum utilization in cellular networks. However, the explicit impact of spatial interference, imposed by FD communications, on uplink and downlink transmissions has been overlooked in the literature. This paper presents an extensive study of the explicit effect of FD communications on the uplink and downlink performances. For the sake of rigorous analysis, we develop a tractable framework based on stochastic geometry toolset. The developed model accounts for uplink truncated channel inversion power control in FD cellular networks. The study shows that FD communications improve the downlink throughput at the expense of significant degradation in the uplink throughput. Therefore, we propose a novel fine-grained duplexing scheme, denoted as $\alpha$-duplex scheme, which allows a partial overlap between uplink and downlink frequency bands. To this end, we show that the amount of the overlap can be optimized via adjusting $\alpha$ to achieve a certain design objective.


title: Composites with invisible inclusions: eigenvalues of R-linear problem

abstract: An new eigenvalue $\mathbb R$-linear problem arisen in the theory of metamaterials is stated and constructively investigated for circular non-overlapping inclusions. An asymptotic formula for eigenvalues is deduced when the radii of inclusions tend to zero. The nodal domains conjecture related to univalent eigenfunctions is posed. Demonstration of the conjecture allows to justify that a set of inclusions can be made neutral by surrounding it with an appropriate coating.


title: Exoplanet Transmission Spectroscopy using KMOS

abstract: KMOS (K-Band Multi Object Spectrograph) is a novel integral field spectrograph installed in the VLT's ANTU unit. The instrument offers an ability to observe 24 2.8"$\times$2.8" sub-fields positionable within a 7.2' patrol field, each sub-field producing a spectrum with a 14$\times$14-pixel spatial resolution. The main science drivers for KMOS are the study of galaxies, star formation, and molecular clouds, but its ability to simultaneously measure spectra of multiple stars makes KMOS an interesting instrument for exoplanet atmosphere characterization via transmission spectroscopy. We set to test whether transmission spectroscopy is practical with KMOS, and what are the conditions required to achieve the photometric precision needed, based on observations of a partial transit of WASP-19b, and full transits of GJ 1214b and HD 209458b. Our analysis uses the simultaneously observed comparison stars to reduce the effects from instrumental and atmospheric sources, and Gaussian processes to model the residual systematics. We show that KMOS can, in theory, deliver the photometric precision required for transmission spectroscopy. However, this is shown to require a) pre-imaging to ensure accurate centering and b) a very stable night with optimal observing conditions (seeing $\sim$0.8"). Combining these two factors with the need to observe several transits, each with a sufficient out-of-transit baseline (and with the fact that similar or better precision can be reached with telescopes and instruments with smaller pressure,) we conclude that transmission spectroscopy is not the optimal science case to take advantage of the abilities offered by KMOS and VLT.


title: Strong interface-induced spin-orbit coupling in graphene on WS2

abstract: Interfacial interactions allow the electronic properties of graphene to be modified, as recently demonstrated by the appearance of satellite Dirac cones in the band structure of graphene on hexagonal boron nitride (hBN) substrates. Ongoing research strives to explore interfacial interactions in a broader class of materials in order to engineer targeted electronic properties. Here we show that at an interface with a tungsten disulfide (WS2) substrate, the strength of the spin-orbit interaction (SOI) in graphene is very strongly enhanced. The induced SOI leads to a pronounced low-temperature weak anti-localization (WAL) effect, from which we determine the spin-relaxation time. We find that spin-relaxation time in graphene is two-to-three orders of magnitude smaller on WS2 than on SiO2 or hBN, and that it is comparable to the intervalley scattering time. To interpret our findings we have performed first-principle electronic structure calculations, which both confirm that carriers in graphene-on-WS2 experience a strong SOI and allow us to extract a spin-dependent low-energy effective Hamiltonian. Our analysis further shows that the use of WS2 substrates opens a possible new route to access topological states of matter in graphene-based systems.


title: The Electronic Behavior of Zinc-Finger Protein Binding Sites in the Context of the DNA Extended Ladder Model

abstract: The eukaryotic Cys2His2 zinc finger proteins bind to DNA ubiquitously at highly conserved domains, responsible for gene regulation and the spatial organization of DNA. To study and understand the zinc finger DNA-protein interaction, we use the extended ladder in the DNA model proposed by Zhu, Rasmussen, Balatsky \& Bishop (2007) \cite{Zhu-2007}. Considering one single spinless electron in each nucleotide $\pi$-orbital along a double DNA chain (dDNA), we find a typical pattern for the {\color{black} bottom of the occupied molecular orbital (BOMO)}, highest occupied molecular orbital (HOMO) and lowest unoccupied orbital (LUMO) along the binding sites. We specifically looked at two members of zinc finger protein family: specificity protein 1 (SP1) and early grown response 1 transcription factors (EGR1). When the valence band is filled, we find electrons in the purines along the nucleotide sequence, compatible with the electric charges of the binding amino acids in SP1 and EGR1 zinc finger.


title: Unconventional Bose-Einstein Condensation in a System of Two-species Bosons in the $p$-orbital Bands of a Bipartite Lattice

abstract: In the context of Gross-Pitaevskii theory, we investigate the unconventional Bose-Einstein condensations in the two-species mixture with $p$-wave symmetry in the second band of a bipartite optical lattice. A new imaginary-time propagation method is developed to numerically determine the $p$-orbital condensation. Different from the single-species case, the two-species boson mixture exhibits two non-equivalent complex condensates in the intraspecies-interaction-dominating regime, exhibiting the vortex-antivortex lattice configuration in the charge and spin channels, respectively. When the interspecies interaction is tuned across the SU(2) invariant point, the system undergoes a quantum phase transition toward a checkerboard-like spin density wave state with a real-valued condensate wavefunction. The influence of lattice asymmetry on the quantum phase transition is addressed. Finally, we present a phase-sensitive measurement scheme for experimentally detecting the UBEC in our model.


title: Xe Irradiation of Graphene on Ir(111): From Trapping to Blistering

abstract: Using X-ray photoelectron spectroscopy, thermal desorption spectroscopy, and scanning tunneling microscopy we show that upon keV Xe + irradiation of graphene on Ir(111), Xe atoms are trapped under the graphene. Upon annealing, aggregation of Xe leads to graphene bulges and blisters. The efficient trapping is an unexpected and remarkable phenomenon, given the absence of chemical binding of Xe to Ir and to graphene, the weak interaction of a perfect graphene layer with Ir(111), as well as the substantial damage to graphene due to irradiation. By combining molecular dynamics simulations and density functional theory calculations with our experiments, we uncover the mechanism of trapping. We describe ways to avoid blister formation during graphene growth, and also demonstrate how ion implantation can be used to intentionally create blisters without introducing damage to the graphene layer. Our approach may provide a pathway to synthesize new materials at a substrate - 2D material interface or to enable confined reactions at high pressures and temperatures.


title: Quaternions in Hamiltonian dynamics of a rigid body -- Part I

abstract: This paper showed that Poisson brackets in quaternion variables can be obtained directly from canonical Poisson brackets on cotangent bundle of $SE(3)$ (or $SO(3)$) endowed by canonical symplectic geometry. Quaternion parameters in our case are just dynamic variables in canonical Hamiltonian mechanics of a rigid body on $T^*SE(3)$ The obtained results based on quaternions representation as explicit functions of rotation matrix elements of $SO(3)$ group. The relation of obtained Poisson structure to the canonical Poisson and symplectic structures on $T^*S^3$ were investigated. To derive the motion equations of Hamiltonian dynamics in quaternionic variables it is proposed to use the mixed frame of reference where translational degrees of freedom describes in the inertial frame of reference and degree of rotational freedom in the body frame. It turns out that motion equations for system with Hamiltonian of the rather general form can be written in algebraic operations on quaternions.


title: Chimera states in networks of phase oscillators: the case of two small populations

abstract: Chimera states are dynamical patterns in networks of coupled oscillators in which regions of synchronous and asynchronous oscillation coexist. Although these states are typically observed in large ensembles of oscillators and analyzed in the continuum limit, chimeras may also occur in systems with finite (and small) numbers of oscillators. Focusing on networks of $2N$ phase oscillators that are organized in two groups, we find that chimera states, corresponding to attracting periodic orbits, appear with as few as two oscillators per group and demonstrate that for $N>2$ the bifurcations that create them are analogous to those observed in the continuum limit. These findings suggest that chimeras, which bear striking similarities to dynamical patterns in nature, are observable and robust in small networks that are relevant to a variety of real-world systems.


title: On the minimal modules for exceptional Lie algebras: Jordan blocks and stabilisers

abstract: Let G be a simple simple-connected exceptional algebraic group of type G_2, F_4, E_6 or E_7 over an algebraically closed field k of characteristic p>0 with \g=Lie(G). For each nilpotent orbit G.e of \g, we list the Jordan blocks of the action of e on the minimal induced module V_min of \g. We also establish when the centralisers G_v of vectors v\in V_min and stabilisers \Stab_G<v> of 1-spaces <v>\subset V_min are smooth; that is, when \dim G_v=\dim\g_v or \dim \Stab_G<v>=\dim\Stab_\g<v>.


title: Identification of Insurance Models with Multidimensional Screening

abstract: This paper addresses the identification of insurance models with multidimensional screening where insurees have private information about their risk and risk aversion. The model includes a random damage and the possibility of several claims. Screening of insurees relies on their certainty equivalence. The paper then investigates how data availability on the number of offered coverages and reported claims affects the identification of the model primitives under four different scenarios. We show that the model structure is identified despite bunching due to multidimensional screening and/or a finite number of offered coverages. The observed number of claims plays a key role in the identification of the joint distribution of risk and risk aversion. In addition, the paper derives all the restrictions imposed by the model on observables. Our results are constructive with explicit equations for estimation and model testing.


title: On melting and freezing for the 2d radial Stefan problem

abstract: We consider the two dimensional free boundary Stefan problem describing the evolution of a spherically symmetric ice ball $\{r\leq \lambda(t)\}$. We revisit the pioneering analysis of [20] and prove the existence in the radial class of finite time melting regimes $$ \lambda(t)=\left\{\begin{array}{ll} (T-t)^{1/2}e^{-\frac{\sqrt{2}}{2}\sqrt{|\ln(T-t)|}+O(1)}\\ (c+o(1))\frac{(T-t)^{\frac{k+1}{2}}}{|\ln (T-t)|^{\frac{k+1}{2k}}}, \ \ k\in \Bbb N^*\end{array}\right. \quad\text{ as } t\to T $$ which respectively correspond to the fundamental stable melting rate, and a sequence of codimension $k\in \Bbb N^*$ excited regimes. Our analysis fully revisits a related construction for the harmonic heat flow in [42] by introducing a new and canonical functional framework for the study of type II (i.e. non self similar) blow up. We also show a deep duality between the construction of the melting regimes and the derivation of a discrete sequence of global-in-time freezing regimes $$ \lambda_\infty - \lambda(t)\sim\left\{\begin{array}{ll} \frac{1}{\log t}\\ \frac{1}{t^{k}(\log t)^{2}}, \ \ k\in \Bbb N^*\end{array}\right. \quad\text{ as } t\to +\infty $$ which correspond respectively to the fundamental stable freezing rate, and excited regimes which are codimension $k$ stable.


title: Simultaneous electronic and the magnetic excitation of a ferromagnet by intense THz pulses

abstract: The speed of magnetization reversal is a key feature in magnetic data storage. Magnetic fields from intense THz pulses have been recently shown to induce small magnetization dynamics in Cobalt thin film on the sub-picosecond time scale. Here, we show that at higher field intensities, the THz electric field starts playing a role, strongly changing the dielectric properties of the cobalt thin film. Both the electronic and magnetic responses are found to occur simultaneously, with the electric field response persistent on a time scale orders of magnitude longer than the THz stimulus


title: Tunable Fermi surface topology and Lifshitz transition in bilayer graphene

abstract: Bilayer graphene is a highly tunable material: not only can one tune the Fermi energy using standard gates, as in single-layer graphene, but the band structure can also be modified by external perturbations such as transverse electric fields or strain. We review the theoretical basics of the band structure of bilayer graphene and study the evolution of the band structure under the influence of these two external parameters. We highlight their key role concerning the ease to experimentally probe the presence of a Lifshitz transition, which consists in a change of Fermi contour topology as a function of energy close to the edges of the conduction and valence bands. Using a device geometry that allows the application of exceptionally high displacement fields, we then illustrate in detail the way to probe the topology changes experimentally using quantum Hall effect measurements in a gapped bilayer graphene system.


title: Gauge Freedom in Path Integrals in Abelian Gauge Theory

abstract: We extend gauge symmetry of Abelian gauge field to incorporate quantum gauge degrees of freedom. We twice apply the Harada--Tsutsui gauge recovery procedure to gauge-fixed theories. First, starting from the Faddeev--Popov path integral in the Landau gauge, we recover the gauge symmetry by introducing an additional field as an extended gauge degree of freedom. Fixing the extended gauge symmetry by the usual Faddeev--Popov procedure, we obtain the theory of Type I gaugeon formalism. Next, applying the same procedure to the resulting gauge-fixed theory, we obtain a theory equivalent to the extended Type I gaugeon formalism.


title: Point Group Symmetry and Deformation Induced Symmetry Breaking of Superlattice Materials

abstract: The point group symmetry of materials is closely related to their physical properties and quite important for material modeling. However, superlattice materials have more complex symmetry conditions than crystals due to their multilevel structural feature. Thus, a theoretical framework is proposed to characterize and determine the point group symmetry of nonmagnetic superlattice materials systematically. A variety of examples are presented to show the symmetry features of superlattice materials in different dimensions and scales. In addition, the deformation induced symmetry breaking phenomenon is also studied for superlattice materials, which has potential application in tuning physical properties by imposing a strain field.


title: RCR: Robust Compound Regression for Robust Estimation of Errors-in-Variables Model

abstract: The errors-in-variables (EIV) regression model, being more realistic by accounting for measurement errors in both the dependent and the independent variables, is widely adopted in applied sciences. The traditional EIV model estimators, however, can be highly biased by outliers and other departures from the underlying assumptions. In this paper, we develop a novel nonparametric regression approach - the robust compound regression (RCR) analysis method for the robust estimation of EIV models. We first introduce a robust and efficient estimator called least sine squares (LSS). Taking full advantage of both the new LSS method and the compound regression analysis method developed in our own group, we subsequently propose the RCR approach as a generalization of those two, which provides a robust counterpart of the entire class of the maximum likelihood estimation (MLE) solutions of the EIV model, in a 1-1 mapping. Technically, our approach gives users the flexibility to select from a class of RCR estimates the optimal one with a predefined regression efficiency criterion satisfied. Simulation studies and real-life examples are provided to illustrate the effectiveness of the RCR approach.


title: Two-dimensional quantum transport of multivalley (111) surface state in topological crystalline insulator SnTe thin films

abstract: Magneto-transport properties of (111)-oriented single-crystal thin films of SnTe were investigated. SnTe (111) thin films were epitaxially grown on a BaF2 substrate by molecular beam epitaxy. By optimizing the growth conditions and the thickness of the films, the bulk carrier density could be reduced, making it possible to detect the surface transport. In the magneto-conductance (MC) measurement, a cusp-like feature around zero magnetic field was observed, which is attributed to the weak-antilocalization effect of the transport in the topological surface state. Detailed analysis of this negative MC reveals a reduced number of transport channels contributing to the surface transport, suggesting a strong coupling between Dirac valleys on the SnTe (111) surface, as a characteristic feature of the transport in the multivalley structure of TCI.


title: Sylvester's Theorem and the Non-integrality of a certain Binomial Sum

abstract: In this note, we show that $S(n,r):=\sum_{k=0}^{n} \binom{n}{k}\frac{k}{k+r}$ is not an integer for any positive integer $n$ and $r\in \{1,2,3,4,5,6\}$ and for $n\le r-1$. This gives a partial answer to a conjecture of [3].


title: Karl Weierstrass Bicentenary

abstract: Academic biography of Karl Weierstrass, his basic works, influence of his doctrine on the development of mathematics.


title: On the codimension of subalgebras of the algebra of matrices over a field

abstract: In this paper we provide an elementary and easy proof that a proper subalgebra of the matrix algebra $ \mathbb{K}^{n,n }$, with $n \geq3$ and $\mathbb{K}$ an arbitrary field, has dimension strictly less than $n^2-1$.


title: End-to-End Optimization of High Throughput DNA Sequencing

abstract: At the core of high throughput DNA sequencing platforms lies a bio-physical surface process that results in a random geometry of clusters of homogenous short DNA fragments typically hundreds of base pairs long - bridge amplification. The statistical properties of this random process and length of the fragments are critical as they affect the information that can be subsequently extracted, i.e., density of successfully inferred DNA fragment reads. The ensemble of overlapping DNA fragment reads are then used to computationally reconstruct the much longer target genome sequence, e.g, ranging from hundreds of thousands to billions of base pairs. The success of the reconstruction in turn depends on having a sufficiently large ensemble of DNA fragments that are sufficiently long. In this paper using stochastic geometry we model and optimize the end-to-end process linking and partially controlling the statistics of the physical processes to the success of the computational step. This provides, for the first time, a framework capturing salient features of such sequencing platforms that can be used to study cost, performance or sensitivity of the sequencing process.


title: Constraints on the distribution of supernova remnants with Galactocentric radius

abstract: Supernova remnants (SNRs) in the Galaxy are an important source of energy injection into the interstellar medium, and also of cosmic rays. Currently there are 294 known SNRs in the Galaxy, and their distribution with Galactocentric radius is of interest for various studies. Here I discuss some of the statistics of Galactic SNRs, including the observational selection effects that apply, and difficulties in obtaining distances for individual remnants from the `Sigma-D' relation. Comparison of the observed Galactic longitude distribution of a sample of bright Galactic SNRs -- which are not strongly affected by selection effects -- with those expected from models is used to constrain the Galactic distribution of SNRs. The best-fitting power-law/exponential model is more concentrated towards the Galactic centre than the widely used distribution obtained by Case & Bhattacharya (1998).


title: All of zeros of Riemann's Zeta-Function are on $\sigma$=1/2

abstract: The research shows that Riemann proved that all of zeros of Riemann's zeta function are on $\sigma=1/2$ based on the functional equation \begin{align*} \pi^{-\frac{s}{2}}\Gamma \left( \frac{s}{2} \right) \zeta(s)&={\frac{1}{s(s-1)} + \int\limits_1^\infty \psi(x) \left( x^{\frac{s}{2} - 1} + x^{-\frac{1+s}{2}} \right) \,dx,}\quad\qquad{s}=\sigma+it, \end{align*} which is in Riemann's ``\"{U}ber die Anzahl der Primzahlen unter einer gegebenen Grosse". According to the geometric meaning of the functional equation and the argument principle, we obtain the number of zeros $N_0(T)$ of the Riemann zeta function on the critical segment $\sigma=1/2,0\leq{t}\leq{T}$ and the number of zeros $N(T)$ of the Riemann zeta function in the rectangular region $-1\leq\sigma\leq{2},0\leq{t}\leq{T}$, respectively. The result is \begin{align*} N(T)&=N_0(T)=\frac{\arg{\left[\pi^{-\frac{s}{2}}\Gamma\left(\frac{s}{2} \right)\zeta(s)\right]}}{\pi}+1\\ &=\frac{T}{2\pi}\log\frac{T}{2\pi}-\frac{T}{2\pi}+O(\log{T}),\qquad{s=1/2+iT}. \end{align*}


title: No Regret Bound for Extreme Bandits

abstract: Algorithms for hyperparameter optimization abound, all of which work well under different and often unverifiable assumptions. Motivated by the general challenge of sequentially choosing which algorithm to use, we study the more specific task of choosing among distributions to use for random hyperparameter optimization. This work is naturally framed in the extreme bandit setting, which deals with sequentially choosing which distribution from a collection to sample in order to minimize (maximize) the single best cost (reward). Whereas the distributions in the standard bandit setting are primarily characterized by their means, a number of subtleties arise when we care about the minimal cost as opposed to the average cost. For example, there may not be a well-defined "best" distribution as there is in the standard bandit setting. The best distribution depends on the rewards that have been obtained and on the remaining time horizon. Whereas in the standard bandit setting, it is sensible to compare policies with an oracle which plays the single best arm, in the extreme bandit setting, there are multiple sensible oracle models. We define a sensible notion of "extreme regret" in the extreme bandit setting, which parallels the concept of regret in the standard bandit setting. We then prove that no policy can asymptotically achieve no extreme regret.


title: Permutations $r_j$ such that $\sum_{i=1}^n \prod_{j=1}^k r_j(i)$ is maximized or minimized

abstract: We consider the problem of finding the set of permutations $r_j$ of $\{1,\cdots , n\}$ such that $\sum_{i=1}^n \prod_{j=1}^k r_j(i)$ is maximized or minimized. While the set of permutations maximizing this value are easily determined, finding the set of permutations minimizing this value appears to be an open problem. We show values of $k$ and $n$ for which an explicit solution exists and comment on computational issues in determining the general problem. We also look at the dual problem of finding the permutations such that $\prod_{i=1}^n \sum_{j=1}^k r_j(i)$ is maximized or minimized. As part of this study we also look at a variant of a rearrangement inequality.


title: Newton's method in practice: finding all roots of polynomials of degree one million efficiently

abstract: We use Newton's method to find all roots of several polynomials in one complex variable of degree up to and exceeding one million and show that the method, applied to appropriately chosen starting points, can be turned into an algorithm that can be applied routinely to find all roots without deflation and with the inherent numerical stability of Newton's method. We specify an algorithm that provably terminates and finds all roots of any polynomial of arbitrary degree, provided all roots are distinct and exact computation is available. It is known that Newton's method is inherently stable, so computing errors do not accumulate; we provide an exact bound on how much numerical precision is sufficient.


abstract: We study $L^\infty$-variational problems associated to measurable Finsler structures in Euclidean spaces. We obtain existence and uniqueness results for the absolute minimizers.


title: An overview of some recent results on the Euler system of isentropic gas dynamics

abstract: In this overview we discuss some recent results of non--uniqueness for the isentropic Euler equations of gas dynamics with particular attention to the role of some admissibility criteria proposed in the literature.


title: Uniqueness of Solution of the Unsteady Filtration Problem in Heterogeneous Porous Media

abstract: We establish uniqueness of the solution of the unsteady state dam problem in the heterogeneous and rectangular case assuming the dam wet at the bottom and dry near to the top.


title: Scale invariant Volkov-Akulov Supergravity

abstract: A scale invariant Goldstino theory coupled to Supergravity is obtained as a standard supergravity dual of a rigidly scale invariant higher--curvature Supergravity with a nilpotent chiral scalar curvature. The bosonic part of this theory describes a massless scalaron and a massive axion in a de Sitter Universe.


title: Mixed Integer Reformulations of Integer Programs and the Affine TU-dimension of a Matrix

abstract: We study the reformulation of integer linear programs by means of a mixed integer linear program with fewer integer variables. Such reformulations can be solved efficiently with mixed integer linear programming techniques. We exhibit examples that demonstrate how integer programs can be reformulated using far fewer integer variables. To this end, we introduce a generalization of total unimodularity called the \emph{affine TU-dimension} of a matrix and study related theory and algorithms for determining the affine TU-dimension of a matrix. We also present bounds on the number of integer variables needed to represent certain integer hulls.


title: Atomically flat carbon monolayer as an extremely unstable quasi-2D mesoscopic quantum system

abstract: The carbon monolayer band structure calculated in the approximation of weakly interacting {\pi} electrons corresponds to massless electron excitations known as Dirac fermions not previously observed in any other material. However, if strong Coulomb and exchange interactions between {\pi} electrons are taken into account, another picture of the {\pi} electron state emerges. These interactions result in {\pi} electron localization and electron crystal formation. The atomically flat layer can be regarded as a simplest quasi-two-dimensional mesoscopic quantum system consisting of a carbon ion plane and two {\pi} electron crystals on opposite sides of the plane. Such a system must have dielectric and pronounced diamagnetic properties and a high sensitivity to external factors distorting its electron crystals. The instability manifests itself in a tendency of the monolayer to be transformed into a more stable carbon modification with a rolled-up or wrapped-up carbon skeleton which is observed as a monolayer corrugation (or ripples). The corrugated monolayer is characterized by the presence of excited {\pi} electrons which are responsible for its physical and chemical properties. The approach can prove useful for investigation of carbon nanotubes, fullerenes, graphite, and also topological insulators and complicated quantum systems with a layered structure. Calculations of the quasi-two-dimensional quantum system faces a many-body problem and are beyond the band-structure description, which forces us to confine ourselves to purely qualitative analysis


title: Galton's Family Heights Data Revisited

abstract: Galton's family heights data has been a preeminent historical dataset in regression analysis, on which the original model and basic results have survived the close scrutiny of statisticians for 125 years. However by revisiting Galton's family data, we challenge whether Galton's classic model and his regression towards mean interpretation are proper. Using Galton's data as a benchmark for different regression methods, such as least squares, orthogonal regression, geometric mean regression, and least sine squares regression - a newly developed nonparametric robust regression approach, we elucidate that his regression model has fundamental drawbacks not only in variable and model selection by "transmuting" women into men thus the simple linear model, but also a strong bias in least squares regression leading to otherwise alternative conclusions on the true relationships between the heights of the child and his or her parents.


title: The canonical central exact sequence for locally compact quantum groups

abstract: For a locally compact quantum group $\mathbb{G}$ we define its center, $\mathscr{Z}(\mathbb{G})$, and its quantum group of inner automorphisms, $\mathrm{Inn}(\mathbb{G})$. We show that one obtains a natural isomorphism between $\mathrm{Inn}(\mathbb{G})$ and $\mathbb{G}/\!\mathscr{Z}(\mathbb{G})$, we characterize normal quantum subgroups of a compact quantum group as those left invariant by the action of the quantum group of inner automorphisms and discuss several examples.


title: Topological phases of shaken quantum Ising lattices

abstract: The quantum compass model consists of a two-dimensional square spin lattice where the orientation of the spin-spin interactions depends on the spatial direction of the bonds. It has remarkable symmetry properties and the ground state shows topological degeneracy. The implementation of the quantum compass model in quantum simulation setups like ultracold atoms and trapped ions is far from trivial, since spin interactions in those sytems typically are independent of the spatial direction. Ising spin interactions, on the contrary, can be induced and controlled in atomic setups with state-of-the art experimental techniques. In this work, we show how the quantum compass model on a rectangular lattice can be simulated by the use of the photon-assisted tunneling induced by periodic drivings on a quantum Ising spin model. We describe a procedure to adiabatically prepare one of the doubly-degenerate ground states of this model by adiabatically ramping down a transverse magnetic field, with surprising differences depending on the parity of the lattice size. Exact diagonalizations confirm the validity of this approach for small lattices. Specific implementations of this scheme are presented with ultracold atoms in optical lattices in the Mott insulator regime, as well as with Rydberg atoms.


title: Gene expression dynamics with stochastic bursts: exact results for a coarse-grained model

abstract: We present a theoretical framework to analyze the dynamics of gene expression with stochastic bursts. Beginning with an individual-based model which fully accounts for the messenger RNA (mRNA) and protein populations, we propose a novel expansion of the master equation for the joint process. The resulting coarse-grained model reduces the dimensionality of the system, describing only the protein population while fully accounting for the effects of discrete and fluctuating mRNA population. Closed form expressions for the stationary distribution of the protein population and mean first-passage times of the coarse-grained model are derived and large-scale Monte Carlo simulations show that the analysis accurately describes the individual-based process accounting for mRNA population, in contrast to the failure of commonly proposed diffusion-type models.


title: A Hausdorff dimension for finite sets

abstract: The classical Hausdorff dimension of finite or countable sets is zero. We define an analog for finite sets, called finite Hausdorff dimension which is non-trivial. It turns out that a finite bound for the finite Hausdorff dimension guarantees that every point of the set has "nearby" neighbors. This property is important for many computer algorithms of great practical value, that obtain solutions by finding nearest neighbors. We also define an analog for finite sets of the classical box-counting dimension, and compute examples. The main result of the paper is a Convergence Theorem. It gives conditions under which, if a sequence of finite sets converges to a compact set (convergence of compact subsets of Euclidean space under the Hausdorff metric), then the finite Hausdorff dimension of the finite sets will converge to the classical Hausdorff dimension of the compact set.


title: Ab initio transport results for strongly correlated fermions

abstract: Quantum transport of strongly correlated fermions is of central interest in condensed matter physics. Here, we present first-principle nonequilibrium Green functions results using $T$-matrix selfenergies for finite Hubbard clusters of dimension $1,2,3$. We compute the expansion dynamics following a potential quench and predict its dependence on the interaction strength and particle number. We discover a universal scaling, allowing an extrapolation to infinite-size systems, which shows excellent agreement with recent cold atom diffusion experiments [Schneider et al., Nat. Phys. 8, 213 (2012)].


title: A LOFAR census of millisecond pulsars

abstract: We report the detection of 48 millisecond pulsars (MSPs) out of 75 observed thus far using the LOFAR in the frequency range 110-188 MHz. We have also detected three MSPs out of nine observed in the frequency range 38-77 MHz. This is the largest sample of MSPs ever observed at these low frequencies, and half of the detected MSPs were observed for the first time at frequencies below 200 MHz. We present the average pulse profiles of the detected MSPs, their effective pulse widths, and flux densities and compare these with higher observing frequencies. The flux-calibrated, multifrequency LOFAR pulse profiles are publicly available via the EPN Database of Pulsar Profiles. We also present average values of dispersion measures (DM) and discuss DM and profile variations. About 35% of the MSPs show strong narrow profiles, another 25% exhibit scattered profiles, and the rest are only weakly detected. A qualitative comparison of LOFAR profiles with those at higher radio frequencies shows constant separation between profile components. Similarly, the profile widths are consistent with those observed at higher frequencies, unless scattering dominates at the lowest frequencies. This is very different from what is observed for normal pulsars and suggests a compact emission region in the MSP magnetosphere. The amplitude ratio of the profile components, on the other hand, can dramatically change towards low frequencies, often with the trailing component becoming dominant. As previously demonstrated this can be caused by aberration and retardation. This data set enables high-precision studies of pulse profile evolution with frequency, dispersion, Faraday rotation, and scattering in the interstellar medium. Characterising and correcting these systematic effects may improve pulsar-timing precision at higher observing frequencies, where pulsar timing array projects aim to directly detect gravitational waves.


title: Conformal weights of charged Renyi entropy twist operators for free scalar fields

abstract: I compute the conformal weights of the twist operators of free scalar fields for charged R\'enyi entropy in both odd and even dimensions. Explicit expressions can be found, in odd dimensions as a function of the chemical potential in the absence of a conical singularity and thence by images for all integer coverings. This method, developed some time ago, is equivalent, in results, to the replica technique. A review is given. The same method applies for even dimensions but a general form is more immediately available. For no chemical potential, the the closed form in the covering order is written in an alternative way related to old trigonometric sums. Some derivatives are obtained. An analytical proof is given of a conjecture made by Bueno, Myers and Witczac--Krempa regarding the relation between the conformal weights and a corner coefficient (a universal quantity) in the R\'enyi entropy.


title: XMM-Newton Large Program on SN1006 - I: Methods and Initial Results of Spatially-Resolved Spectroscopy

abstract: Based on our newly developed methods and the XMM-Newton large program of SN1006, we extract and analyze the spectra from 3596 tessellated regions of this SNR each with 0.3-8 keV counts $>10^4$. For the first time, we map out multiple physical parameters, such as the temperature ($kT$), electron density ($n_e$), ionization parameter ($n_et$), ionization age ($t_{ion}$), metal abundances, as well as the radio-to-X-ray slope ($\alpha$) and cutoff frequency ($\nu_{cutoff}$) of the synchrotron emission. We construct probability distribution functions of $kT$ and $n_et$, and model them with several Gaussians, in order to characterize the average thermal and ionization states of such an extended source. We construct equivalent width (EW) maps based on continuum interpolation with the spectral model of each regions. We then compare the EW maps of OVII, OVIII, OVII K$\delta-\zeta$, Ne, Mg, SiXIII, SiXIV, and S lines constructed with this method to those constructed with linear interpolation. We further extract spectra from larger regions to confirm the features revealed by parameter and EW maps, which are often not directly detectable on X-ray intensity images. For example, O abundance is consistent with solar across the SNR, except for a low-abundance hole in the center. This "O Hole" has enhanced OVII K$\delta-\zeta$ and Fe emissions, indicating recently reverse shocked ejecta, but also has the highest $n_et$, indicating forward shocked ISM. Therefore, a multi-temperature model is needed to decompose these components. The asymmetric metal distributions suggest there is either an asymmetric explosion of the SN or an asymmetric distribution of the ISM.


title: Light Pseudoscalar Mesons in Bethe-Salpeter Equation with Instantaneous Interaction

abstract: The light pseudoscalar mesons play a twofold role: they may or have to be regarded both as low-lying bound states of the fundamental degrees of freedom of quantum chromodynamics as well as the (pseudo-) Goldstone bosons of the spontaneously broken chiral symmetries of quantum chromodynamics. We interrelate these aspects in a single quantum-field-theoretic approach relying on the Bethe-Salpeter formalism in instantaneous approximation by very simple means: the shape of the pseudoscalar-meson Bethe-Salpeter wave function dictated by chiral symmetry is used in Bethe-Salpeter equations for bound states of vanishing mass, in order to deduce analytically the interactions which govern the bound states under study. In this way, we obtain exact Bethe-Salpeter solutions for pseudoscalar mesons, in the sense of establishing the rigorous relationship between, on the one hand, the relevant interactions and, on the other hand, the Bethe-Salpeter amplitudes that characterize the bound states.


title: An inverse theorem for compact Lipschitz regions in $R^d$ using localized kernel bases

abstract: While inverse estimates in the context of radial basis function approximation on boundary-free domains have been known for at least ten years, such theorems for the more important and difficult setting of bounded domains have been notably absent. This article develops inverse estimates for finite dimensional spaces arising in radial basis function approximation and meshless methods. The inverse estimates we consider control Sobolev norms of linear combinations of a localized basis by the $L_p$ norm over a bounded domain. The localized basis is generated by forming local Lagrange functions for certain types of RBFs (namely Mat\'ern and surface spline RBFs). In this way it extends the boundary-free construction of Fuselier, Hangelbroek, Narcowich, Ward and Wright.


title: Transport Through a Correlated Interface: Auxiliary Master Equation Approach

abstract: We present improvements of a recently introduced numerical method [Arrigoni etal, Phys. Rev. Lett. 110, 086403 (2013)] to compute steady state properties of strongly correlated electronic systems out of equilibrium. The method can be considered as a non-equilibrium generalization of exact diagonalization based dynamical mean-field theory (DMFT). The key modification for the non-equilibrium situation consists in addressing the DMFT impurity problem within an auxiliary system consisting of the correlated impurity, $N_b$ uncorrelated bath sites and two Markovian environments (sink and reservoir). Algorithmic improvements in the impurity solver allow to treat efficiently larger values of $N_b$ than previously in DMFT. This increases the accuracy of the results and is crucial for a correct description of the physical behavior of the system in the relevant parameter range including a semi-quantitative description of the Kondo regime. To illustrate the approach we consider a monoatomic layer of correlated orbitals, described by the single-band Hubbard model, attached to two metallic leads. The non-equilibrium situation is driven by a bias-voltage applied to the leads. For this system, we investigate the spectral function and the steady state current-voltage characteristics in the weakly as well as in the strongly interacting limit. In particular we investigate the non-equilibrium behavior of quasi-particle excitations within the Mott gap of the correlated layer. We find for low bias voltage Kondo like behavior in the vicinity of the insulating phase. In particular we observe a splitting of the Kondo resonance as a function of the bias voltage.


title: Minimal Length Maximal Green Sequences and Triangulations of Polygons

abstract: We use combinatorics of quivers and the corresponding surfaces to study maximal green sequences of minimal length for quivers of type $\mathbb{A}$. We prove that such sequences have length $n+t$, where $n$ is the number of vertices and $t$ is the number of 3-cycles in the quiver. Moreover, we develop a procedure that yields these minimal length maximal green sequences.


title: Effect of scattering on the transonic solution topology and intrinsic variability of line-driven stellar winds

abstract: For line-driven winds from hot, luminous OB stars, we examine the subtle but important role of diffuse, scattered radiation in determining both the topology of steady-state solutions and intrinsic variability in the transonic wind base. We use a smooth source function formalism to obtain nonlocal, integral expressions for the direct and diffuse components of the line-force that account for deviations from the usual localized, Sobolev forms. As the scattering source function is reduced, we find the solution topology in the transonic region transitions from X-type, with a unique wind solution, to a nodal type, characterized by a degenerate family of solutions. Specifically, in the idealized case of an optically thin source function and a uniformly bright stellar disk, the unique X-type solution proves to be a stable attractor to which time-dependent numerical radiation-hydrodynamical simulations relax. But in models where the scattering strength is only modestly reduced, the topology instead turns nodal, with the associated solution degeneracy now manifest by intrinsic structure and variability that persist down to the photospheric wind base. This highlights the potentially crucial role of diffuse radiation for the dynamics and variability of line-driven winds, and seriously challenges the use of Sobolev theory in the transonic wind region. Since such Sobolev-based models are commonly used in broad applications like stellar evolution and feedback, this prompts development of new wind models, not only for further quantifying the intrinsic variability found here, but also for computing new theoretical predictions of global properties like velocity laws and mass-loss rates.


title: Design of a low noise, wide band, active dipole antenna for a cosmic ray radiodetection experiment

abstract: An active dipole antenna has been designed to measure transient electric field induced by ultra high energy cosmic rays for the CODALEMA experiment. The main requirements for this detector, composed of a low noise preamplifier placed close to a dipole antenna, are a wide bandwidth ranging from 100 kHz to 100 MHz and a good sensitivity on the whole spectrum.


title: Statistical framework for estimating GNSS bias

abstract: We present a statistical framework for estimating global navigation satellite system (GNSS) non-ionospheric differential time delay bias. The biases are estimated by examining differences of measured line integrated electron densities (TEC) that are scaled to equivalent vertical integrated densities. The spatio-temporal variability, instrumentation dependent errors, and errors due to inaccurate ionospheric altitude profile assumptions are modeled as structure functions. These structure functions determine how the TEC differences are weighted in the linear least-squares minimization procedure, which is used to produce the bias estimates. A method for automatic detection and removal of outlier measurements that do not fit into a model of receiver bias is also described. The same statistical framework can be used for a single receiver station, but it also scales to a large global network of receivers. In addition to the Global Positioning System (GPS), the method is also applicable to other dual frequency GNSS systems, such as GLONASS (Globalnaya Navigazionnaya Sputnikovaya Sistema). The use of the framework is demonstrated in practice through several examples. A specific implementation of the methods presented here are used to compute GPS receiver biases for measurements in the MIT Haystack Madrigal distributed database system. Results of the new algorithm are compared with the current MIT Haystack Observatory MAPGPS bias determination algorithm. The new method is found to produce estimates of receiver bias that have reduced day-to-day variability and more consistent coincident vertical TEC values.


title: Algorithmic Design of Majorizers for Large-Scale Inverse Problems

abstract: Iterative majorize-minimize (MM) (also called optimization transfer) algorithms solve challenging numerical optimization problems by solving a series of "easier" optimization problems that are constructed to guarantee monotonic descent of the cost function. Many MM algorithms replace a computationally expensive Hessian matrix with another more computationally convenient majorizing matrix. These majorizing matrices are often generated using various matrix inequalities, and consequently the set of available majorizers is limited to structures for which these matrix inequalities can be efficiently applied. In this paper, we present a technique to algorithmically design matrix majorizers with wide varieties of structures. We use a novel duality-based approach to avoid the high computational and memory costs of standard semidefinite programming techniques. We present some preliminary results for 2D X-ray CT reconstruction that indicate these more exotic regularizers may significantly accelerate MM algorithms.


title: Mountain Peak Detection in Online Social Media

abstract: We present a system for the classification of mountain panoramas from user-generated photographs followed by identification and extraction of mountain peaks from those panoramas. We have developed an automatic technique that, given as input a geo-tagged photograph, estimates its FOV (Field Of View) and the direction of the camera using a matching algorithm on the photograph edge maps and a rendered view of the mountain silhouettes that should be seen from the observer's point of view. The extraction algorithm then identifies the mountain peaks present in the photograph and their profiles. We discuss possible applications in social fields such as photograph peak tagging on social portals, augmented reality on mobile devices when viewing a mountain panorama, and generation of collective intelligence systems (such as environmental models) from massive social media collections (e.g. snow water availability maps based on mountain peak states extracted from photograph hosting services).


title: Pore-scale lattice Boltzmann simulation of laminar and turbulent flow through a sphere pack

abstract: The lattice Boltzmann method can be used to simulate flow through porous media with full geometrical resolution. With such a direct numerical simulation, it becomes possible to study fundamental effects which are difficult to assess either by developing macroscopic mathematical models or experiments. We first evaluate the lattice Boltzmann method with various boundary handling of the solid-wall and various collision operators to assess their suitability for large scale direct numerical simulation of porous media flow. A periodic pressure drop boundary condition is used to mimic the pressure driven flow through the simple sphere pack in a periodic domain. The evaluation of the method is done in the Darcy regime and the results are compared to a semi-analytic solution. Taking into account computational cost and accuracy, we choose the most efficient combination of the solid boundary condition and collision operator. We apply this method to perform simulations for a wide range of Reynolds numbers from Stokes flow over seven orders of magnitude to turbulent flow. Contours and streamlines of the flow field are presented to show the flow behavior in different flow regimes. Moreover, unknown parameters of the Forchheimer, the Barree--Conway and friction factor models are evaluated numerically for the considered flow regimes.


title: A Dzyaloshinskii-Moriya Anisotropy in Nanomagnets with in-plane Magnetization

abstract: We report on a new source of in-plane anisotropy in nanomagnets due to the presence of Dzyaloshinskii-Moriya interaction (DMI). This anisotropy depends on the shape of the magnet, and is orthogonal to the demagnetization shape anisotropy. This effect originates from the DMI energy reduction due to an out-of-plane tilt of the spins at edges oriented perpendicular to the magnetization. Our investigation combining experimental, numerical and analytical results demonstrate that this energy reduction can compensate the demagnetization shape anisotropy energy in magnets of elongated shape, provided that their volumes is small enough and thus that their magnetization is quasi-uniform.


title: The Dirichlet Problem with Prescribed Asymptotic Singularities

abstract: We solve the nonlinear Dirichlet problem (uniquely) for functions with prescribed asymptotic singularities at a finite number of points, and with arbitrary continuous boundary data, on a domain in euclidean space. The main results apply, in particular, to subequations with a Riesz characteristic $p \geq 2$. In this case it is shown that, without requiring uniform ellipticity, the Dirichlet problem can be solved uniquely for arbitrary continuous boundary data with singularities asymptotic to the Riesz kernel: $\Theta_j K_p(x - x_j)$, where $K_p(x) = - {1\over|x|^{p-2}}$ for $p>2$ and $K_2(x) = \log |x|$, at any prescribed finite set of points $x_1,...,x_k$ in the domain and any finite set of positive real numbers $\Theta_1,..., \Theta_k$. This sharpens a previous result of the authors concerning the discreteness of high-density sets of subsolutions. Uniqueness and existence results are also established for finite-type singularities such as $\Theta_j |x - x_j|^{2-p}$ for $1\leq p<2$. The main results apply similarly with prescribed singularities asymptotic to the fundamental solutions of Armstrong-Sirakov-Smart (in the uniformly elliptic case).


title: The varieties of Heisenberg vertex operator algebras

abstract: For a vertex operator algebra $V$ with conformal vector $\omega$, we consider a class of vertex operator subalgebras and their conformal vectors. They are called semi-conformal vertex operator subalgebras and semi-conformal vectors of $(V,\omega)$, respectively, and were used to study duality theory of vertex operator algebras via coset constructions. Using these objects attached to $(V,\omega)$, we shall understand the structure of the vertex operator algebra $(V,\omega)$. At first, we define the set $\on{Sc}(V,\omega)$ of semi-conformal vectors of $V$, then we prove that $\on{Sc}(V,\omega)$ is an affine algebraic variety with a partial ordering and an involution map. Corresponding to each semi-conformal vector, there is a unique maximal semi-conformal vertex operator subalgebra containing it. The properties of these subalgebras are invariants of vertex operator algebras. As an example, we describe the corresponding varieties of semi-conformal vectors for Heisenberg vertex operator algebras. As an application, we give two characterizations of Heisenberg vertex operator algebras using the properties of these varieties.


title: Binary words avoiding xx^Rx and strongly unimodal sequences

abstract: In previous work, Currie and Rampersad showed that the growth of the number of binary words avoiding the pattern xxx^R was intermediate between polynomial and exponential. We now show that the same holds for the growth of the number of binary words avoiding the pattern xx^Rx. Curiously, the analysis for xx^Rx is much simpler than that for xxx^R. We derive our results by giving a bijection between the set of binary words avoiding xx^Rx and a class of sequences closely related to the class of "strongly unimodal sequences."


title: Linearly constrained evolutions of critical points and an application to cohesive fractures

abstract: We introduce a novel constructive approach to define time evolution of critical points of an energy functional. Our procedure, which is different from other more established approaches based on viscosity approximations in infinite dimension, is prone to efficient and consistent numerical implementations, and allows for an existence proof under very general assumptions. We consider in particular rather nonsmooth and nonconvex energy functionals, provided the domain of the energy is finite dimensional. Nevertheless, in the infinite dimensional case study of a cohesive fracture model, we prove a consistency theorem of a discrete-to-continuum limit. We show that a quasistatic evolution can be indeed recovered as a limit of evolutions of critical points of finite dimensional discretizations of the energy, constructed according to our scheme. To illustrate the results, we provide several numerical experiments both in one and two dimensions. These agree with the crack initiation criterion, which states that a fracture appears only when the stress overcomes a certain threshold, depending on the material.


title: Quantifying Double McCormick

abstract: When using the standard McCormick inequalities twice to convexify trilinear monomials, as is often the practice in modeling and software, there is a choice of which variables to group first. For the important case in which the domain is a nonnegative box, we calculate the volume of the resulting relaxation, as a function of the bounds defining the box. In this manner, we precisely quantify the strength of the different possible relaxations defined by all three groupings, in addition to the trilinear hull itself. As a by product, we characterize the best double McCormick relaxation. We wish to emphasize that, in the context of spatial branch-and-bound for factorable formulations, our results do not only apply to variables in the input formulation. Our results apply to monomials that involve auxiliary variables as well. So, our results apply to the product of any three (possibly complicated) expressions in a formulation.


title: Quantum anomalous Hall effect and related topological electronic states

abstract: Over a long period of exploration, the successful observation of quantized version of anomalous Hall effect (AHE) in thin film of magnetically-doped topological insulator completed a quantum Hall trio---quantum Hall effect (QHE), quantum spin Hall effect (QSHE), and quantum anomalous Hall effect (QAHE). On the theoretical front, it was understood that intrinsic AHE is related to Berry curvature and U(1) gauge field in momentum space. This understanding established connection between the QAHE and the topological properties of electronic structures characterized by the Chern number. With the time reversal symmetry broken by magnetization, a QAHE system carries dissipationless charge current at edges, similar to the QHE where an external magnetic field is necessary. The QAHE and corresponding Chern insulators are also closely related to other topological electronic states, such as topological insulators and topological semimetals, which have been extensively studied recently and have been known to exist in various compounds. First-principles electronic structure calculations play important roles not only for the understanding of fundamental physics in this field, but also towards the prediction and realization of realistic compounds. In this article, a theoretical review on the Berry phase mechanism and related topological electronic states in terms of various topological invariants will be given with focus on the QAHE and Chern insulators. We will introduce the Wilson loop method and the band inversion mechanism for the selection and design of topological materials, and discuss the predictive power of first-principles calculations. Finally, remaining issues, challenges and possible applications for future investigations in the field will be addressed.


title: Space-efficient detection of unusual words

abstract: Detecting all the strings that occur in a text more frequently or less frequently than expected according to an IID or a Markov model is a basic problem in string mining, yet current algorithms are based on data structures that are either space-inefficient or incur large slowdowns, and current implementations cannot scale to genomes or metagenomes in practice. In this paper we engineer an algorithm based on the suffix tree of a string to use just a small data structure built on the Burrows-Wheeler transform, and a stack of $O(\sigma^2\log^2 n)$ bits, where $n$ is the length of the string and $\sigma$ is the size of the alphabet. The size of the stack is $o(n)$ except for very large values of $\sigma$. We further improve the algorithm by removing its time dependency on $\sigma$, by reporting only a subset of the maximal repeats and of the minimal rare words of the string, and by detecting and scoring candidate under-represented strings that $\textit{do not occur}$ in the string. Our algorithms are practical and work directly on the BWT, thus they can be immediately applied to a number of existing datasets that are available in this form, returning this string mining problem to a manageable scale.


title: Anharmonic nuclear motion and the relative stability of hexagonal and cubic ice

abstract: We use extensive first-principles quantum mechanical calculations to show that, although the static lattice and harmonic vibrational energies are almost identical, the anharmonic vibrational energy of hexagonal ice is significantly lower than that of cubic ice. This difference in anharmonicity is crucial, stabilising hexagonal ice compared with cubic ice by at least 1.4 meV/H2O, in agreement with experimental estimates. The difference in anharmonicity arises predominantly from molecular O-H bond stretching vibrational modes and is related to the different stacking of atomic layers.


title: Realizing orbit categories as stable module categories - a complete classification

abstract: We classify all triangulated orbit categories of path-algebras of Dynkin diagrams that are triangle equivalent to a stable module category of a representation-finite self-injective standard algebra. For each triangulated orbit category T we give an explicit description of a representation-finite self-injective standard algebra with stable module category triangle equivalent to T.


title: Regularized Non-Gaussian Image Denoising

abstract: In image denoising problems, one widely-adopted approach is to minimize a regularized data-fit objective function, where the data-fit term is derived from a physical image acquisition model. Typically the regularizer is selected with two goals in mind: (a) to accurately reflect image structure, such as smoothness or sparsity, and (b) to ensure that the resulting optimization problem is convex and can be solved efficiently. The space of such regularizers in Gaussian noise settings is well studied; however, non-Gaussian noise models lead to data-fit expressions for which entirely different families of regularizers may be effective. These regularizers have received less attention in the literature because they yield non-convex optimization problems in Gaussian noise settings. This paper describes such regularizers and a simple reparameterization approach that allows image reconstruction to be accomplished using efficient convex optimization tools. The proposed approach essentially modifies the objective function to facilitate taking advantage of tools such as proximal denoising routines. We present examples of imaging denoising under exponential family (Bernoulli and Poisson) and multiplicative noise.


title: Harmonic maps and para-Sasakian geometry

abstract: The purpose of this paper is to study the harmonicity of maps to or from para-Sasakian manifolds. We derive the condition for the tension field of paraholomorphic map between almost para-Hermitian manifold and para-Sasakian manifold. The necessary and sufficient condition for a paraholomorphic map between para-Sasakian manifolds to be parapluriharmonic are shown and a non-trivial example is presented for its illustrations.


title: On the Effect of Bias Estimation on Coverage Accuracy in Nonparametric Inference

abstract: Nonparametric methods play a central role in modern empirical work. While they provide inference procedures that are more robust to parametric misspecification bias, they may be quite sensitive to tuning parameter choices. We study the effects of bias correction on confidence interval coverage in the context of kernel density and local polynomial regression estimation, and prove that bias correction can be preferred to undersmoothing for minimizing coverage error and increasing robustness to tuning parameter choice. This is achieved using a novel, yet simple, Studentization, which leads to a new way of constructing kernel-based bias-corrected confidence intervals. In addition, for practical cases, we derive coverage error optimal bandwidths and discuss easy-to-implement bandwidth selectors. For interior points, we show that the MSE-optimal bandwidth for the original point estimator (before bias correction) delivers the fastest coverage error decay rate after bias correction when second-order (equivalent) kernels are employed, but is otherwise suboptimal because it is too "large". Finally, for odd-degree local polynomial regression, we show that, as with point estimation, coverage error adapts to boundary points automatically when appropriate Studentization is used; however, the MSE-optimal bandwidth for the original point estimator is suboptimal. All the results are established using valid Edgeworth expansions and illustrated with simulated data. Our findings have important consequences for empirical work as they indicate that bias-corrected confidence intervals, coupled with appropriate standard errors, have smaller coverage error and are less sensitive to tuning parameter choices in practically relevant cases where additional smoothness is available.


title: Birationally rigid Pfaffian Fano 3-folds

abstract: We classify birationally rigid orbifold Fano 3-folds of index one defined by $5 \times 5$ Pfaffians. We give a sharp criterion for birational rigidity of these families based on the type of singularities that the varieties admit. Various conjectures are born out of our study, highlighting a possible approach to the classification of terminal Fano 3-folds. The birationally rigid cases are the first known rigid examples of Fanos that are not (weighted) complete intersection.


title: Permutation totally symmetric self-complementary plane partitions

abstract: Alternating sign matrices and totally symmetric self-complementary plane partitions are equinumerous sets of objects for which no explicit bijection is known. In this paper, we identify a subset of totally symmetric self-complementary plane partitions corresponding to permutations by giving a statistic-preserving bijection to permutation matrices, which are a subset of alternating sign matrices. We use this bijection to define a new partial order on permutations, and prove this new poset contains both the Tamari lattice and the Catalan distributive lattice as subposets. We also study a new partial order on totally symmetric self-complementary plane partitions arising from this perspective and show that this is a distributive lattice related to Bruhat order when restricted to permutations.


title: Electronic structure and correlations of vitamin B12 studied within the Haldane-Anderson impurity model

abstract: We study the electronic structure and correlations of vitamin B12 (cyanocobalamine) by using the framework of the multi-orbital single-impurity Haldane-Anderson model of a transition-metal impurity in a semiconductor host. The parameters of the effective Haldane-Anderson model are obtained within the Hartree-Fock (HF) approximation. The quantum Monte Carlo (QMC) technique is then used to calculate the one-electron and magnetic correlation functions of this effective model. We observe that new states form inside the semiconductor gap found by HF due to the intra-orbital Coulomb interaction at the impurity 3d orbitals. In particular, the lowest unoccupied states correspond to an impurity bound state, which consists of states from mainly the CN axial ligand and the corring ring as well as the Co e_g-like orbitals. We also observe that the Co(3d) orbitals can develop antiferromagnetic correlations with the surrounding atoms depending on the filling of the impurity bound states. In addition, we make comparisons of the HF+QMC data with the density functional theory calculations. We also discuss the photoabsorption spectrum of cyanocobalamine.


title: A massively parallel multi-level approach to a domain decomposition method for the optical flow estimation with varying illumination

abstract: We consider a variational method to solve the optical flow problem with varying illumination. We apply an adaptive control of the regularization parameter which allows us to preserve the edges and fine features of the computed flow. To reduce the complexity of the estimation for high resolution images and the time of computations, we implement a multi-level parallel approach based on the domain decomposition with the Schwarz overlapping method. The second level of parallelism uses the massively parallel solver MUMPS. We perform some numerical simulations to show the efficiency of our approach and to validate it on classical and real-world image sequences.


title: Scarred eigenstates for arithmetic toral point scatterers

abstract: We investigate eigenfunctions of the Laplacian perturbed by a delta potential on the standard tori $\mathbb{R}^d/2 \pi\mathbb{Z}^d$ in dimensions $d=2,3$. Despite quantum ergodicity holding for the set of "new" eigenfunctions we show that there is scarring in the momentum representation for $d=2,3$, as well as in the position representation for $d=2$ (i.e., the eigenfunctions fail to equidistribute in phase space along an infinite subsequence of new eigenvalues.) For $d=3$, scarred eigenstates are quite rare, but for $d=2$ scarring in the momentum representation is very common --- with $N_{2}(x) \sim x/\sqrt{\log x}$ denoting the counting function for the new eigenvalues below $x$, there are $\gg N_{2}(x)/\log^A x$ eigenvalues corresponding to momentum scarred eigenfunctions.


title: Holomorphic families of $[\lambda]-$primitive themes

abstract: This article is the continuation of [B. 13-b] where we show how the isomorphism class of a $[\lambda]-$primitive theme with a given Bernstein polynomial may be characterized by a (small) finite number of complex parameters. We construct here a corresponding locally versal holomorphic deformation of $ [\lambda]-$primitive themes for each given Bernstein polynomial. Then we prove the universality of the corresponding "canonical family" in many cases. We also give some examples where no local universal family exists.