title: Possible Ergodic-nonergodic regions in the quantum Sherrington-Kirkpatrick spin glass model and quantum annealing

abstract: We explore the behavior of order parameter distribution of quantum Sherrington-Kirkpatrick model in the spin glass phase using Monte Carlo technique for the effective Suzuki-Trotter Hamil- tonian at finite temperatures and that at zero temperature obtained using exact diagonalization method. Our numerical results indicate the existence of low but finite temperature quantum fluc- tuation dominated ergodic region along with the classical fluctuation dominated high temperature nonergodic region in the spin glass phase of the model. In the ergodic region, the order parameter distribution gets narrower around the most probable value of the order parameter as the system size increases. In the other region, the Parisi order distribution function has non-vanishing value every- where in thermodynamic limit, indicating nonergodicity. We also show, that the average annealing time for convergence (to a low energy level of the model; within a small error range) becomes system size independent for annealing down through the (quantum fluctuation dominated) ergodic region. It becomes strongly system size dependent for annealing through the nonergodic region. Possible finite size scaling type behavior for the extent of the ergodic region is also addressed.


title: A simple test for stability of black hole by $S$-deformation

abstract: We study a sufficient condition to prove the stability of a black hole when the master equation for linear perturbation takes the form of the Schr\"odinger equation. If the potential contains a small negative region, usually, the $S$-deformation method was used to show the non-existence of unstable mode. However, in some cases, it is hard to find an appropriate deformation function analytically because the only way known so far to find it is a try-and-error approach. In this paper, we show that it is easy to find a regular deformation function by numerically solving the differential equation such that the deformed potential vanishes everywhere, when the spacetime is stable. Even if the spacetime is almost marginally stable, our method still works. We also discuss a simple toy model which can be solved analytically, and show the condition for the non-existence of a bound state is the same as that for the existence of a regular solution for the differential equation in our method. From these results, we conjecture that our criteria is also a necessary condition.


title: Generalized Entanglement Measure for Continuous Variable Systems

abstract: Concurrence introduced by Hill and Wootters [Phys. Rev. Lett. 78, 5022 (1997)], provides an important measure of entanglement for a general pair of qubits that is strictly positive entangled states and vanishes for all separable states. We present an extension of entanglement measure to general continuous variable states of multiple degrees of freedom by generalizing the Lagrange's identity and wedge product framework proposed by Bhaskara et al. [Quantum Inf. Process. 16, 118 (2017)] for pure discrete variable systems in arbitrary dimensions. A family of faithful entanglement measures is constructed that admit necessary and sufficient conditions for separability across arbitrary bipartitions presented by Vedral et al. [Phys. Rev. Lett. 16 2275 (1997)]. The computed entanglement measure in the present approach for general Gaussian states, pair coherent states and non Gaussian continuous variable Bell states, matches with known results. We also quantify entanglement of phase randomized squeezed states and superposition of squeezed states. Our results also simplify several results in quantum entanglement theory.


title: To Index or Not to Index: Optimizing Exact Maximum Inner Product Search

abstract: Exact Maximum Inner Product Search (MIPS) is an important task that is widely pertinent to recommender systems and high-dimensional similarity search. The brute-force approach to solving exact MIPS is computationally expensive, thus spurring recent development of novel indexes and pruning techniques for this task. In this paper, we show that a hardware-efficient brute-force approach, blocked matrix multiply (BMM), can outperform the state-of-the-art MIPS solvers by over an order of magnitude, for some -- but not all -- inputs. In this paper, we also present a novel MIPS solution, MAXIMUS, that takes advantage of hardware efficiency and pruning of the search space. Like BMM, MAXIMUS is faster than other solvers by up to an order of magnitude, but again only for some inputs. Since no single solution offers the best runtime performance for all inputs, we introduce a new data-dependent optimizer, OPTIMUS, that selects online with minimal overhead the best MIPS solver for a given input. Together, OPTIMUS and MAXIMUS outperform state-of-the-art MIPS solvers by 3.2$\times$ on average, and up to 10.9$\times$, on widely studied MIPS datasets.


title: A Joint Model for Question Answering and Question Generation

abstract: We propose a generative machine comprehension model that learns jointly to ask and answer questions based on documents. The proposed model uses a sequence-to-sequence framework that encodes the document and generates a question (answer) given an answer (question). Significant improvement in model performance is observed empirically on the SQuAD corpus, confirming our hypothesis that the model benefits from jointly learning to perform both tasks. We believe the joint model's novelty offers a new perspective on machine comprehension beyond architectural engineering, and serves as a first step towards autonomous information seeking.


title: A novel method for the absolute energy calibration of large-scale cosmic-ray detectors using radio emission of extensive air showers

abstract: Ultra-high energy cosmic rays impinging onto the atmosphere induce huge cascades of secondary particles. The measurement of the energy radiated by these air showers in form of radio waves enables an accurate measurement of the cosmic-ray energy. Compared to the well-established fluorescence technique, the radio measurements are less dependent on atmospheric conditions and thus potentially reduce the systematic uncertainty in the cosmic-ray energy measurement significantly. Two attractive aspects are that the atmosphere is transparent to MHz radio waves and the radio emission can be calculated from first-principles using classical electrodynamics. This method will be discussed for the Engineering Radio Array (AERA) of the Pierre Auger Cosmic-Ray Observatory. AERA detects radio emission from extensive air showers with energies beyond $10^{17}~$eV in the 30 - 80 MHz frequency band and consists of more than 150 autonomous radio stations covering an area of about 17$~$km$^2$. It is located at the same site as the Auger low-energy detector extensions enabling combinations with various other measurement techniques.


title: End-Fire Silicon Optical Phased Array with Half-Wavelength Spacing

abstract: We demonstrate a one-dimensional optical phased array on an integrated silicon platform for operation at 1.55 microns. Light is emitted end-fire from the chip edge where the waveguides are terminated. The innovative design and high confinement afforded by the silicon waveguides enables lambda/2 spacing (775-nm pitch) at the output. Steering is achieved by inducing a phase shift between the waveguides via integrated thermo-optic heaters. The device forms a beam with a FWHM angular width of 17 degrees, and we demonstrate beam steering over a 64 degrees range.


title: Implementation, demonstration and validation of a user-defined wall-function for direct precipitation fouling in ANSYS Fluent

abstract: In a previous paper (Johnsen et al., 2015) and presentation (Johnsen et al., 2016), we developed and demonstrated a generic modelling framework for the modelling of direct precipitation fouling from multi-component fluid mixtures that become super-saturated at the wall. The modelling concept involves the 1-dimensional transport of the fluid species through the turbulent boundary layer close to the wall. The governing equations include the Reynolds-averaged (RANS) advection-diffusion equations for each fluid species, and the axial momentum and energy equations for the fluid mixture. The driving force for the diffusive transport is the local gradient in the species' chemical potential. Adsorption mechanisms are not modelled per se, but the time-scale of adsorption is reflected in the choice of Dirichlet boundary conditions for the depositing species, at the fluid-solid interface. In this paper, the modelling framework is implemented as a user-defined function (UDF) for the CFD software ANSYS Fluent, to act as a wall boundary condition for mass-transfer to the wall. The subgrid, 1-dimensional formulation of the model reduces the computational cost associated with resolving the fine length-scales at which the boundary-layer mass transfer is determined, and allows for efficient modelling of industry-scale heat exchangers suffering from fouling. The current paper describes the modelling framework, and demonstrates and validates its applicability in a simplified 2D heat exchanger geometry (experimental and detailed CFD modelling data by P\"a\"akk\"onen et al. (2012, 2016)). By tuning the diffusivity, only, good agreement with the experimental data and the detailed CFD model was obtained, in terms of area-averaged deposition rates.


title: Combinatorial identities generated by difference analogs of hyperbolic and trigonometric functions of order $n$

abstract: We naturally obtain some combinatorial identities finding the difference analogs of hyperbolic and trigonometric functions of order $n.$ In particular, we obtain the identities connected with the proved in the paper the addition formulas for these analogs.


title: Metriplectic formalism: friction and much more

abstract: The metriplectic formalism couples Poisson brackets of the Hamiltonian description with metric brackets for describing systems with both Hamiltonian and dissipative components. The construction builds in asymptotic convergence to a preselected equilibrium state. Phenomena such as friction, electric resistivity, thermal conductivity and collisions in kinetic theories are well represented in this framework. In this paper we present an application of the metriplectic formalism of interest for the theory of control: a suitable torque is applied to a free rigid body, which is expressed through a metriplectic extension of its "natural" Poisson algebra. On practical grounds, the effect is to drive the body to align its angular velocity to rotation about a stable principal axis of inertia, while conserving its kinetic energy in the process. On theoretical grounds, this example shows how the non-Hamiltonian part of a metriplectic system may include convergence to a limit cycle, the first example of a non-zero dimensional attractor in this formalism. The method suggests a way to extend metriplectic dynamics to systems with general attractors, e.g. chaotic ones, with the hope of representing bio-physical, geophysical and ecological models.


title: Stability and decay estimates for the marine riser equation in the presence of time dependent data

abstract: In this article we investigate the dynamics of the initial-boundary value problem for the nonlinear marine riser equation in the presence of time dependent boundary conditions at the top end and a time dependent coefficient of the nonlinear drag force. We introduce sufficient conditions on these functions to maintain the structural stability of the system. We deduce their maximum rates of growth to guarantee that the zero solution is globally asymptotically stable.


title: Erasing and Correction of Liquid Metal Printed Electronics Made of Gallium Alloy Ink from the Substrate

abstract: Gallium-based liquid metals have recently been found important in a variety of newly emerging applications such as room temperature metal 3D printing, direct writing of electronics and biomedicine etc. In all these practices, one frequently encounters the situations that a printed circuit or track needs to be corrected or the unwanted parts of the device should be removed as desired. However, few appropriate strategies are currently available to tackle such important issues at this stage. Here we have identified several low cost ways toward this goal by comparatively investigating three typical strategies spanning from mechanical, chemical, to electrochemical principles, for removing the gallium-based liquid metal circuits or thin films. Regarding the mechanical approach, we constructed an eraser for removing the liquid metal thin films. It was shown that ethanol (CH3CH2OH) could serve as a good candidacy material for the mechanical eraser. In the chemical category, we adopted alkalis and acids to remove the finely printed liquid metal circuits and sodium hydroxide (NaOH) solution was particularly revealed to be rather efficient in making a chemical eraser. In the electrochemical strategy, we applied a 15 V voltage to a liquid metal thin film (covered with water) and successfully removed the target metal part. These methods were comparatively evaluated with each of the merits and shortcomings preliminarily clarified in the end. The present work is expected to be important for the increasing applications of the liquid metal enabled additive manufactures.


title: Intelligence of agents produces a structural phase transition in collective behaviour

abstract: Living organisms process information to interact and adapt to their changing environment with the goal of finding food, mates or averting hazards. The structure of their niche has profound repercussions by both selecting their internal architecture and also inducing adaptive responses to environmental cues and stimuli. Adaptive, collective behaviour underpinned by specialized optimization strategies is ubiquitously found in the natural world. This exceptional success originates from the processes of fitness and selection. Here we prove that a universal physical mechanism of a nonequilibrium transition underlies the collective organization of information-processing organisms. As cognitive agents build and update an internal, cognitive representation of the causal structure of their environment, complex patterns emerge in the system, where the onset of pattern formation relates to the spatial overlap of cognitive maps. Studying the exchange of information among the agents reveals a continuous, order-disorder transition. As a result of the spontaneous breaking of translational symmetry, a Goldstone mode emerges, which points at a collective mechanism of information transfer among cognitive organisms. Taken together, the characteristics of this phase transition consolidate different results in cognitive and biological sciences in a universal manner. These finding are generally applicable to the design of artificial intelligent swarm systems that do not rely on centralized control schemes.


title: Information entropy and thermal entropy: apples and oranges

abstract: The frequent misunderstanding of information entropy is pointed out. It is shown that, contrary to fortuitous situations and common beliefs, there is no general interrelation between the information entropy and the thermodynamical entropy. We point out that the change of information entropy during measurement is determined by the resolution of the measurement instrument and these changes do not have a general and clear-cut interrelation with the change of thermal entropy. Moreover, these changes can be separated in space and time. We show classical physical considerations and examples to support this conclusion. Finally, we show that information entropy can violate the Third Law of Thermodynamics which is another indication of major differences from thermal entropy.


title: Extremely late photometry of SN~2011fe

abstract: Type Ia supernovae are widely accepted to be the outcomes of thermonuclear explosions in white dwarf stars. However, many details of these explosions remain uncertain (e.g. the mass, ignition mechanism, and flame speed). Theory predicts that at very late times (beyond 1000 d) it might be possible to distinguish between explosion models. Few very nearby supernovae can be observed that long after the explosion. The Type Ia supernova SN 2011fe located in M101 and along a line of sight with negligible extinction, provides us with the once-in-a-lifetime chance to obtain measurements that may distinguish between theoretical models. In this work, we present the analysis of photometric data of SN 2011fe taken between 900 and 1600 days after explosion with Gemini and HST. At these extremely late epochs theory suggests that the light curve shape might be used to measure isotopic abundances which is a useful model discriminant. However, we show in this work that there are several currently not well constrained physical processes introducing large systematic uncertainties to the isotopic abundance measurement. We conclude that without further detailed knowledge of the physical processes at this late stage one cannot reliably exclude any models on the basis of this dataset.


title: Reconciling the dawn-dusk asymmetry in Mercury's exosphere with the micrometeoroid impact directionality

abstract: Combining dynamical models of dust from Jupiter Family Comets and Halley-type Comets, we demonstrate that the seasonal variation of the dust/meteoroid environment at Mercury is responsible for producing the dawn-dusk asymmetry in Mercury's exosphere observed by the MESSENGER spacecraft. Our latest models, calibrated recently from ground-based and space-borne measurements, provide unprecedented statistics that enable us to study the longitudinal and latitudinal distribution of meteoroids impacting Mercury's surface. We predict that the micrometeoroid impact vaporization source is expected to undergo significant motion on Mercury's surface towards the nightside during Mercury's approach to aphelion and towards the dayside when the planet is approaching the Sun.


title: Enhancing quantum order with fermions by increasing species degeneracy

abstract: One of the challenges for fermionic cold atom experiments in optical lattices is to cool the systems to low enough temperature that they can form quantum degenerate ordered phases. In particular, there has been significant work in trying to find the antiferromagnetic phase transition of the Hubbard model in three dimensions, without success. Here, we attack this problem from a different angle by enhancing the ordering temperature via an increase in the degeneracy of the atomic species trapped in the optical lattice. In addition to developing the general theory, we also discuss some potential systems where one might be able to achieve these results experimentally.


title: Phenomenology of flavorful composite vector bosons in light of $B$ anomalies

abstract: We analyze the flavor structure of composite vector bosons arising in a model of vectorlike technicolor, often called hypercolor (HC), with eight flavors that form a one-family content of HC fermions. Dynamics of the composite vector bosons, referred to as HC rho in this paper, are formulated together with HC pions by the hidden local symmetry (HLS), in a way analogous to QCD vector mesons. Then coupling properties to the standard model (SM) fermions, which respect the HLS gauge symmetry, are described in a way that couplings of the HC rhos to the left-handed SM quarks and leptons are given by a well-defined setup as taking the flavor mixing structures into account. Under the present scenario, we discuss significant bounds on the model from electroweak precision tests, flavor physics, and collider physics. We also try to address B anomalies in processes such as B -> K(*) mu+ mu- and B -> D(*) tau nu, recently reported by LHCb, Belle, (ATLAS, and CMS in part.) Then we find that the present model can account for the anomaly in B -> K(*) mu+ mu- consistently with the other constraints while it predicts no significant deviations in B -> D(*) tau nu from the SM, which can be examined in the future Belle II experiment. The former is archived with the form C9 = -C10 of the Wilson coefficients for effective operators of b -> s mu+ mu-, which has been favored by the recent experimental data. We also investigate current and future experimental limits at the Large Hadron Collider (LHC) and see that possible collider signals come from dijet and ditau, or dimuon resonant searches for the present scenario with TeV mass range. To conclude, the present b -> s mu+ mu- anomaly is likely to imply discovery of new vector bosons in the ditau or dimuon channel in the context of the HC rho model. Our model can be considered as a UV completion of conventional U(1)' model.


title: Using binary statistics in Taurus-Auriga to distinguish between brown dwarf formation processes

abstract: Whether BDs form as stars through gravitational collapse ("star-like") or BDs and some very low-mass stars constitute a separate population which form alongside stars comparable to the population of planets, e.g. through circumstellar disk ("peripheral") fragmentation, is one of the key questions of the star-formation problem. For young stars in Taurus-Auriga the binary fraction is large with little dependence on primary mass above ~0.2Msun, while for BDs it is <10%. We investigate a case in which BDs in Taurus formed dominantly through peripheral fragmentation. The decline of the binary frequency in the transition region between star-like and peripheral formation is modelled. A dynamical population synthesis model is employed in which stellar binary formation is universal. Peripheral objects form separately in circumstellar disks with a distinctive initial mass function (IMF), own orbital parameter distributions for binaries and a low binary fraction. A small amount of dynamical processing of the stellar component is accounted for as appropriate for the low-density Taurus-Auriga embedded clusters. The binary fraction declines strongly between the mass-limits for star-like and peripheral formation. The location of characteristic features and the steepness depend on these mass-limits. Such a trend might be unique to low density regions hosting dynamically unprocessed binary populations. The existence of a strong decline in the binary fraction -- primary mass diagram will become verifiable in future surveys on BD and VLMS binarity in the Taurus-Auriga star forming region. It is a test of the (non-)continuity of star formation along the mass-scale, the separateness of the stellar and BD populations and the dominant formation channel for BDs and BD binaries in regions of low stellar density hosting dynamically unprocessed populations.


title: Extracting the sigma-term from low-energy pion-nucleon scattering

abstract: We present an extraction of the pion-nucleon ($\pi N$) scattering lengths from low-energy $\pi N$ scattering, by fitting a representation based on Roy-Steiner equations to the low-energy data base. We show that the resulting values confirm the scattering-length determination from pionic atoms, and discuss the stability of the fit results regarding electromagnetic corrections and experimental normalization uncertainties in detail. Our results provide further evidence for a large $\pi N$ $\sigma$-term, $\sigma_{\pi N}=58(5)$ MeV, in agreement with, albeit less precise than, the determination from pionic atoms.


title: Renormalization in a Lorentz-violating model and higher-order operators

abstract: The renormalization in a Lorentz-breaking scalar-spinor higher-derivative model involving $\phi^4$ self-interaction and the Yukawa-like coupling is studied. We explicitly de- monstrate that the convergence is improved in comparison with the usual scalar-spinor model, so, the theory is super-renormalizable, with no divergences beyond four loops. We compute the one-loop corrections to the propagators for the scalar and fermionic fields and show that in the presence of higher-order Lorentz invariance violation, the poles that dominate the physical theory, are driven away from the standard on-shell pole mass due to radiatively induced lower dimensional operators. The new operators change the standard gamma-matrix structure of the two-point functions, introduce large Lorentz-breaking corrections and lead to modifications in the renormalization conditions of the theory. We found the physical pole mass in each sector of our model.


title: On the effect of Lyman alpha trapping during the initial collapse of massive black hole seeds

abstract: One viable seeding mechanism for supermassive black holes is the direct gaseous collapse route in pre-galactic dark matter halos, producing objects on the order of $10^4 - 10^6$ solar masses. These events occur when the gas is prevented from cooling below $10^4$ K that requires a metal-free and relatively H$_2$-free medium. The initial collapse cools through atomic hydrogen transitions, but the gas becomes optically thick to the cooling radiation at high densities. We explore the effects ofLyman-$\alpha$ trapping in such a collapsing system with a suite of Monte Carlo radiation transport calculations in uniform density and isotropic cases that are based from a cosmological simulation. Our method includes both non-coherent scattering and two-photon line cooling. We find that Lyman-$\alpha$ radiation is marginally trapped in the parsec-scale gravitationally unstable central cloud, allowing the temperature to increase to 50,000 K at a number density of $3 \times 10^4$ cm$^{-3}$ and increasing the Jeans mass by a factor of five. The effective equation of state changes from isothermal at low densities to have an adiabatic index of 4/3 around the temperature maximum and then slowly retreats back to isothermal at higher densities. Our results suggest that Lyman-$\alpha$ trapping delays the initial collapse by raising the Jeans mass. Afterward the high density core cools back to $10^4$ K that is surrounded by a warm envelope whose inward pressure may alter the fragmentation scales at high densities.


title: Spin-polarized Andreev tunneling through the Rashba chain

abstract: We demonstrate that the selective equal spin Andreev reflection (SESAR) spectroscopy can be used in STM experiments to distinguish the zero-energy Majorana quasiparticles from the ordinary fermionic states of the Rashba chain. Such technique, designed for probing the p-wave superconductivity, could be applied to the intersite pairing of equal-spin electrons in the chain of magnetic Fe atoms deposited on the superconducting Pb substrate. Our calculations of the effective pairing amplitude for individual spin components imply the magnetically polarized Andreev conductance, which can be used to `filter' the Majorana quasiparticles from the ordinary in-gap states, although the pure spin current (i.e., perfect polarization) is impossible.


title: The SLUGGS Survey: Trails of SLUGGS galaxies in a modified spin-ellipticity diagram

abstract: We present radial tracks for four early-type galaxies with embedded intermediate-scale discs in a modified spin-ellipticity diagram. Here, each galaxy's spin and ellipticity profiles are shown as a radial track, as opposed to a single, flux-weighted aperture-dependent value as is common in the literature. The use of a single ellipticity and spin parameter is inadequate to capture the basic nature of these galaxies, which transition from fast to slow rotation as one moves to larger radii where the disc ceases to dominate. After peaking, the four galaxy's radial tracks feature a downturn in both ellipticity and spin with increasing radius, differentiating them from elliptical galaxies, and from lenticular galaxies whose discs dominate at large radii. These galaxies are examples of so-called discy elliptical galaxies, which are a morphological hybrid between elliptical (E) and lenticular (S0) galaxies, and have been designated ES galaxies. The use of spin-ellipticity tracks provides extra structural information about individual galaxies over a single aperture measure. Such tracks provide a key diagnostic for classifying early-type galaxies, particularly in the era of 2D kinematic (and photometric) data beyond one effective radius.


title: Intertwined order and holography: the case of the parity breaking pair density wave

abstract: We present a minimal bottom-up extension of the Chern-Simons bulk action for holographic translational symmetry breaking that naturally gives rise to pair density waves. We construct stationary inhomogeneous black hole solutions in which both the U(1) symmetry and spatially translational symmetry are spontaneously broken at finite temperature and charge density. This novel solution provides a dual description of a superconducting phase intertwined with charge, current and parity orders.


title: Constraining the CO intensity mapping power spectrum at intermediate redshifts

abstract: We compile available constraints on the carbon monoxide (CO) 1-0 luminosity functions and abundances at redshifts 0-3. This is used to develop a data driven halo model for the evolution of the CO galaxy abundances and clustering across intermediate redshifts. It is found that the recent constraints from the CO Power Spectrum Survey ($z \sim 3$; Keating et al. 2016), when combined with existing observations of local galaxies ($z \sim 0$; Keres et al. 2003), lead to predictions which are consistent with the results of smaller surveys at intermediate redshifts ($z \sim 1-2$). We provide convenient fitting forms for the evolution of the CO luminosity - halo mass relation, and estimates of the mean and uncertainties in the CO power spectrum in the context of future intensity mapping experiments.


title: Full-sky Gravitational Lensing Simulation for Large-area Galaxy Surveys and Cosmic Microwave Background Experiments

abstract: We present 108 full-sky gravitational lensing simulation data sets generated by performing multiple-lens plane ray-tracing through high-resolution cosmological $N$-body simulations. The data sets include full-sky convergence and shear maps from redshifts $z=0.05$ to $5.3$ at intervals of $150 \, h^{-1}{\rm Mpc}$ comoving radial distance (corresponding to a redshift interval of $\Delta z \simeq 0.05$ at the nearby universe), enabling the construction of a mock shear catalog for an arbitrary source distribution up to $z=5.3$. The dark matter halos are identified from the same $N$-body simulations with enough mass resolution to resolve the host halos of the Sloan Digital Sky Survey (SDSS) CMASS and Luminous Red Galaxies (LRGs). Angular positions and redshifts of the halos are provided by a ray-tracing calculation, enabling the creation of a mock halo catalog to be used for galaxy-galaxy and cluster-galaxy lensing. The simulation also yields maps of gravitational lensing deflections for a source redshift at the last scattering surface, and we provide 108 realizations of lensed cosmic microwave background (CMB) maps in which the post-Born corrections caused by multiple light scattering are included. We present basic statistics of the simulation data, including the angular power spectra of cosmic shear, CMB temperature and polarization anisotropies, galaxy-galaxy lensing signals for halos, and their covariances. The angular power spectra of the cosmic shear and CMB anisotropies agree with theoretical predictions within $5\%$ up to $\ell = 3000$ (or at an angular scale $\theta > 0.5$ arcsin). The simulation data sets are generated primarily for the ongoing Subaru Hyper Suprime-Cam survey but are freely available for download at http://cosmo.phys.hirosaki-u.ac.jp/takahasi/allsky_raytracing.


title: Analysis of vector boson production within TMD factorization

abstract: We present a comprehensive analysis and extraction of the unpolarized transverse momentum dependent (TMD) parton distribution functions, which are fundamental constituents of the TMD factorization theorem. We provide a general review of the theory of TMD distributions, and present a new scheme of scale fixation. This scheme, called the $\zeta$-prescription, allows to minimize the impact of perturbative logarithms in the large range of scales and does not generate undesired power corrections. Within $\zeta$-prescription we consistently include the perturbatively calculable parts up to next-to-next-to-leading order (NNLO), and perform the global fit of the Drell-Yan and Z-boson production, which include the data of E288, Tevatron and LHC experiments. The nonperturbative part of the TMDs are explored checking a large variety of models. We support the obtained results by a study of theoretical uncertainties, perturbative convergence, and dedicated study of the range of applicability of the TMD factorization theorem. The considered nonperturbative models present significant differences in the fitting behavior, which allow us to clearly disfavor most of them. The numerical evaluations are provided by the "arTeMiDe" code, which is introduced in this work and that can be used for current/future TMD phenomenology.


title: Metallicity Gradient of the Thick Disc Progenitor at High Redshift

abstract: We have developed a novel Markov Chain Mote Carlo (MCMC) chemical "painting" technique to explore possible radial and vertical metallicity gradients for the thick disc progenitor. In our analysis we match an N-body simulation to the data from the Apache Point Observatory Galactic Evolution Experiment (APOGEE) survey. We assume that the thick disc has a constant scale-height and has completed its formation at an early epoch, after which time radial mixing of its stars has taken place. Under these assumptions, we find that the initial radial metallicity gradient of the thick disc progenitor should not be negative, but either flat or even positive, to explain the current negative vertical metallicity gradient of the thick disc. Our study suggests that the thick disc was built-up in an inside-out and upside-down fashion, and older, smaller and thicker populations are more metal poor. In this case, star forming discs at different epochs of the thick disc formation are allowed to have different radial metallicity gradients, including a negative one, which helps to explain a variety of slopes observed in high redshift disc galaxies. This scenario helps to explain the positive slope of the metallicity-rotation velocity relation observed for the Galactic thick disc. On the other hand, radial mixing flattens the slope of an existing gradient.


title: Deconfinement transitions in a generalised XY model

abstract: We find the complete phase diagram of a generalised XY model that includes half-vortices. The model possesses superfluid, pair-superfluid and disordered phases, separated by Kosterlitz-Thouless (KT) transitions for both the half-vortices and ordinary vortices, as well as an Ising-type transition. There also occurs an unusual deconfining phase transition, where the disordered to superfluid transition is of Ising rather than KT type. We show by analytical arguments and extensive numerical simulations that there is a point in the phase diagram where the KT transition line meets the deconfining Ising phase transition. We find that the latter extends into the disordered phase not as a phase transition, but rather solely as a deconfinement transition. It is best understood in the dual height model, where on one side of the transition height steps are bound into pairs while on the other they are unbound. We also extend the phase diagram of the dual model, finding both O(2) loop model and antiferromagnetic Ising transitions.


title: The Metallicity of the Intracluster Medium Over Cosmic Time: Further Evidence for Early Enrichment

abstract: We use Chandra X-ray data to measure the metallicity of the intracluster medium (ICM) in 245 massive galaxy clusters selected from X-ray and Sunyaev-Zel'dovich (SZ) effect surveys, spanning redshifts $0<z<1.2$. Metallicities were measured in three different radial ranges, spanning cluster cores through their outskirts. We explore trends in these measurements as a function of cluster redshift, temperature, and surface brightness "peakiness" (a proxy for gas cooling efficiency in cluster centers). The data at large radii (0.5--1 $r_{500}$) are consistent with a constant metallicity, while at intermediate radii (0.1-0.5 $r_{500}$) we see a late-time increase in enrichment, consistent with the expected production and mixing of metals in cluster cores. In cluster centers, there are strong trends of metallicity with temperature and peakiness, reflecting enhanced metal production in the lowest-entropy gas. Within the cool-core/sharply peaked cluster population, there is a large intrinsic scatter in central metallicity and no overall evolution, indicating significant astrophysical variations in the efficiency of enrichment. The central metallicity in clusters with flat surface brightness profiles is lower, with a smaller intrinsic scatter, but increases towards lower redshifts. Our results are consistent with other recent measurements of ICM metallicity as a function of redshift. They reinforce the picture implied by observations of uniform metal distributions in the outskirts of nearby clusters, in which most of the enrichment of the ICM takes place before cluster formation, with significant later enrichment taking place only in cluster centers, as the stellar populations of the central galaxies evolve.


title: Heat content and horizontal mean curvature on the Heisenberg group

abstract: We identify the short time asymptotics of the sub-Riemannian heat content for a smoothly bounded domain in the first Heisenberg group. Our asymptotic formula generalizes prior work by van den Berg-Le Gall and van den Berg-Gilkey to the sub-Riemannian context, and identifies the first few coefficients in the sub-Riemannian heat content in terms of the horizontal perimeter and the total horizontal mean curvature of the boundary. The proof is probabilistic, and relies on a characterization of the heat content in terms of Brownian motion.


title: An optimal $(\epsilon,\delta)$-approximation scheme for the mean of random variables with bounded relative variance

abstract: Randomized approximation algorithms for many #P-complete problems (such as the partition function of a Gibbs distribution, the volume of a convex body, the permanent of a $\{0,1\}$-matrix, and many others) reduce to creating random variables $X_1,X_2,\ldots$ with finite mean $\mu$ and standard deviation$\sigma$ such that $\mu$ is the solution for the problem input, and the relative standard deviation $|\sigma/\mu| \leq c$ for known $c$. Under these circumstances, it is known that the number of samples from the $\{X_i\}$ needed to form an $(\epsilon,\delta)$-approximation $\hat \mu$ that satisfies $\mathbb{P}(|\hat \mu - \mu| > \epsilon \mu) \leq \delta$ is at least $(2-o(1))\epsilon^{-2} c^2\ln(1/\delta)$. We present here an easy to implement $(\epsilon,\delta)$-approximation $\hat \mu$ that uses $(2+o(1))c^2\epsilon^{-2}\ln(1/\delta)$ samples. This achieves the same optimal running time as other estimators, but without the need for extra conditions such as bounds on third or fourth moments.


title: The very faint hard state of the persistent neutron star X-ray binary SLX 1737-282 near the Galactic centre

abstract: We report on a detailed study of the spectral and temporal properties of the neutron star low mass X-ray binary SLX 1737-282, which is located only ~1degr away from Sgr A. The system is expected to have a short orbital period, even within the ultra-compact regime, given its persistent nature at low X-ray luminosities and the long duration thermonuclear burst that it has displayed. We have analysed a Suzaku (18 ks) observation and an XMM-Newton (39 ks) observation taken 7 years apart. We infer (0.5-10 keV) X-ray luminosities in the range 3-6 x10^35erg s-1, in agreement with previous findings. The spectra are well described by a relatively cool (kTbb = 0.5 keV) black body component plus a Comptonized emission component with {\Gamma} ~1.5-1.7. These values are consistent with the source being in a faint hard state, as confirmed by the ~ 20 per cent fractional root mean square amplitude of the fast variability (0.1 - 7 Hz) inferred from the XMM-Newton data. The electron temperature of the corona is >7 keV for the Suzaku observation, but it is measured to be as low as ~2 keV in the XMM-Newton data at higher flux. The latter is significantly lower than expected for systems in the hard state. We searched for X-ray pulsations and imposed an upper limit to their semi-amplitude of 2 per cent (0.001 - 7 Hz). Finally, we investigated the origin of the low frequency variability emission present in the XMM-Newton data and ruled out an absorption dip origin. This constraint the orbital inclination of the system to 65 degr unless the orbital period is longer than 11 hr (i.e. the length of the XMM-Newton observation).


title: High- and low-energy constraints on top-Higgs couplings

abstract: We study the five chirality-flipping interactions that appear in the top-Higgs sector at leading order in the standard model effective field theory. We consider constraints from collider observables, flavor physics, and electric-dipole-moment experiments. This analysis results in very competitive constraints from indirect observables when one considers a single coupling at a time. In addition, we discuss how these limits are affected in scenarios in which multiple top-Higgs interactions are generated at the scale of new physics.


title: The Bright and the Slow -- GRBs 100724B \& 160509A with high-energy cutoffs at $\lesssim100\;$MeV

abstract: We analyze the prompt emission of GRB 100724B and GRB 160509A, two of the brightest Gamma-Ray Bursts (GRBs) observed by Fermi at $\lesssim{\rm MeV}$ energies but surprisingly faint at $\gtrsim100\;$MeV energies. Time-resolved spectroscopy reveals a sharp high-energy cutoff at energies $E_c\sim20-60\;$MeV for GRB~100724B and $E_c\sim80-150\;$MeV for GRB~160509A. We first characterize phenomenologically the cutoff and its time evolution. We then fit the data to two models where the high-energy cutoff arises from intrinsic opacity to pair production within the source ($\tau_{\gamma\gamma}$): (i) a Band spectrum with $\tau_{\gamma\gamma}$ from the internal-shocks motivated model of Granot et al. (2008), and (ii) the photospheric model of Gill & Thompson (2014). Alternative explanations for the cutoff, such as an intrinsic cutoff in the emitting electron energy distribution, appear to be less natural. Both models provide a good fit to the data with very reasonable physical parameters, providing a direct estimate of bulk Lorentz factors in the range $\Gamma\sim 100-400$, on the lower end of what is generally observed in Fermi GRBs. Surprisingly, their lower cutoff energies $E_c$ compared to other Fermi-LAT GRBs arise not predominantly from the lower Lorentz factors, but also at a comparable level from differences in variability time, luminosity, and high-energy photon index. Finally, particularly low $E_c$ values may prevent detection by Fermi-LAT, thus introducing a bias in the Fermi-LAT GRB sample against GRBs with low Lorentz factors or variability times.


title: Information-entropic analysis of Korteweg--de Vries solitons in the quark-gluon plasma

abstract: Solitary waves propagation of baryonic density perturbations, ruled by the Korteweg--de Vries equation in a mean-field quark-gluon plasma model, are investigated from the point of view of the theory of information. A recently proposed continuous logarithmic measure of information, called configurational entropy, is used to derive the soliton width, defining the pulse, for which the informational content of the soliton spatial profile is more compressed, in the Shannon's sense.


title: Autoignition of Methyl Valerate at Low to Intermediate Temperatures and Elevated Pressures in a Rapid Compression Machine

abstract: Methyl valerate (C$_6$H$_{12}$O$_2$, methyl pentanoate) is a methyl ester and a relevant surrogate component for biodiesel. In this work, we present ignition delays of methyl valerate measured using a rapid compression machine at a range of engine-relevant temperature, pressure, and equivalence ratio conditions. The conditions we have studied include equivalence ratios from 0.25 to 2.0, temperatures between 680 K and 1050 K, and pressures of 15 bar and 30 bar. The ignition delay data demonstrate a negative temperature coefficient region in the temperature range of 720 K-800 K for both $\phi$=2.0, 15 bar and $\phi$=1.0, 30 bar, with two-stage ignition apparent over the narrower temperature ranges of 720 K-760 K for the lower pressure and 740 K-800 K at the higher pressure. In addition, the experimental ignition delay data are compared with simulations using an existing chemical kinetic model from the literature. The simulations with the literature model under-predict the data by factors between 2 and 10 over the entire range of the experimental data. To help determine the possible reasons for the discrepancy between simulations and experiments, a new chemical kinetic model is developed using the Reaction Mechanism Generator (RMG) software. The agreement between the experimental data and the RMG model is improved but still not satisfactory. Directions for future improvement of the methyl valerate model are discussed.


title: What if gravity becomes really repulsive in the future?

abstract: The current acceleration of the Universe is one of the most puzzling issues in theoretical physics nowadays. We are far from giving an answer on this letter to its true nature. Yet, with the observations we have at hand, we analyse the different patterns that the gravitational potential can show in the future. Surprisingly, gravity not only can get weaker in the near future, it can even become repulsive; or equivalently, the gravitational potential may become negative. We show this remark by using one of the simplest phenomenological model we can imagine for dark energy. We have as well reviewed the statefinder approach of these models. For completeness, we have also showed the behaviour of the density contrast of dark matter and dark energy for these simple (yet illustrative models). Our results are displayed at present and how they evolve in the future.


title: High-Pressure Autoignition of Binary Blends of Methanol and Dimethyl Ether

abstract: Reactivity Controlled Compression Ignition (RCCI) is a new advanced engine concept that uses a dual fuel mode of operation to achieve significant improvements in fuel economy and emissions output. The fuels that are typically used in this mode include a low- and a high-reactivity fuel in varying proportions to control ignition timing. As such, understanding the interaction effects during autoignition of binary fuel blends is critical to optimizing these RCCI engines. In this work, we measure the autoignition delays of binary blends of dimethyl ether (C$_2$H$_6$O, DME) and methanol (CH$_4$O, MeOH) in a rapid compression machine. In these experiments, dimethyl ether and methanol function as the high- and low-reactivity fuels, respectively. We considered five fuel blends at varying blending ratios (by mole), including 100% DME-0% MeOH, 75% DME-25% MeOH, 0% DME-0% MeOH, 25% DME-75% MeOH, and 0% DME-100% MeOH. Experiments are conducted at an engine-relevant pressure of 30 bar, for the stoichiometric equivalence ratio. In addition, the experimental results are compared with simulations using a chemical kinetic model for DME/MeOH combustion generated by merging independent, well-validated models for DME and MeOH.


title: Curvature-Squared Multiplets, Evanescent Effects and the U(1) Anomaly in N = 4 Supergravity

abstract: We evaluate one-loop amplitudes of N = 4 supergravity in D dimensions using the double-copy procedure that expresses gravity integrands in terms of corresponding ones in Yang--Mills theory.We organize the calculation in terms of a set of gauge-invariant tensors, allowing us to identify evanescent contributions. Among the latter, we find the matrix elements of supersymmetric completions of curvature-squared operators. In addition, we find that such evanescent terms and the U(1)-anomalous contributions to one-loop N = 4 amplitudes are tightly intertwined. The appearance of evanescent operators in N = 4 supergravity and their relation to anomalies raises the question of their effect on the known four-loop divergence in this theory. We provide bases of gauge-invariant tensors and corresponding projectors useful for Yang--Mills theories as a by-product of our analysis.


title: Visual attention models for scene text recognition

abstract: In this paper we propose an approach to lexicon-free recognition of text in scene images. Our approach relies on a LSTM-based soft visual attention model learned from convolutional features. A set of feature vectors are derived from an intermediate convolutional layer corresponding to different areas of the image. This permits encoding of spatial information into the image representation. In this way, the framework is able to learn how to selectively focus on different parts of the image. At every time step the recognizer emits one character using a weighted combination of the convolutional feature vectors according to the learned attention model. Training can be done end-to-end using only word level annotations. In addition, we show that modifying the beam search algorithm by integrating an explicit language model leads to significantly better recognition results. We validate the performance of our approach on standard SVT and ICDAR'03 scene text datasets, showing state-of-the-art performance in unconstrained text recognition.


title: Random Flag Complexes and Asymptotic Syzygies

abstract: We use the probabilistic method to construct examples of conjectured phenomenon about asymptotic syzygies. In particular, we use the Stanley-Reisner ideals of random flag complexes to construct new examples of Ein and Lazarsfeld's nonvanishing for asymptotic syzygies and of Ein, Erman, and Lazarsfeld's conjectural on the asymptotic normal distribution of Betti numbers.


title: The Dynamics of Truncated Black Hole Accretion Disks I: Viscous, Hydrodynamic Case

abstract: Truncated accretion disks are commonly invoked to explain the spectro-temporal variability from accreting black holes in both small systems, i.e. state transitions in galactic black hole binaries (GBHBs), and large systems, i.e. low-luminosity active galactic nuclei (LLAGNs). In the canonical truncated disk model of moderately low accretion rate systems, gas in the inner region of the accretion disk occupies a hot, radiatively inefficient phase, which leads to a geometrically thick disk, while the gas in the outer region occupies a cooler, radiatively efficient phase that resides in the standard geometrically thin disk. Observationally, there is strong empirical evidence to support this phenomenological model, but a detailed understanding of the dynamics of truncated disks is lacking. We present a well-resolved viscous, hydrodynamic simulation that uses an ad hoc cooling prescription to drive a thermal instability and, hence, produce the first sustained truncated accretion disk. With this simulation, we perform a study of the dynamics, angular momentum transport, and energetics of a truncated disk. We find that time variability introduced by the quasi-periodic transition of gas from efficient cooling to inefficient cooling impacts the evolution of the simulated disk. A consequence of the thermal instability is that an outflow is launched from the hot/cold gas interface which drives large, sub-Keplerian convective cells in the disk atmosphere. The convective cells introduce a viscous $\theta-\phi$ stress that is less than the generic $r-\phi$ viscous stress component, but greatly influences the evolution of the disk. In the truncated disk, we find that the bulk of the accreted gas is in the hot phase.


title: Searching for the heavy custodial fiveplet Higgs in the Georgi-Machacek model at the International Linear Collider

abstract: The Georgi-Machacek (GM) model is one of many Beyond Standard Model scenarios with an extended scalar sector which can group under the custodial $\rm SU(2)_C$ symmetry into a fiveplet, a triplet, and two singlets. The heavy charged custodial fiveplet Higgs $\rm H_5^{\pm}$ are typical particles in GM model which couple to the electroweak gauge bosons, therefore provide a good testing ground for the detection of the $\rm H^\pm W^\mp Z$ vertex. The neutral custodial fiveplet Higgs $\rm H_5^{0}$ in the GM model has the same mass with $\rm H_5^{\pm}$ and couples to the electroweak gauge bosons $\rm W^+ W^-$ and ZZ both. We study the discovery prospects of the exotic scalar bosons $\rm H_5^{\pm}$ and $\rm H_5^{0}$ at the International Linear Collider(ILC) via the vector boson associated production processes, and discuss two different decay modes for both charged and neutral scalars. The discovery potential is discussed. Testing the mass degeneracy of charged and neutral scalar bosons in the GM model is also considered.


title: Ultracold atoms in multiple-radiofrequency dressed adiabatic potentials

abstract: We present the first experimental demonstration of a multiple-radiofrequency dressed potential for the configurable magnetic confinement of ultracold atoms. We load cold $^{87}$Rb atoms into a double well potential with an adjustable barrier height, formed by three radiofrequencies applied to atoms in a static quadrupole magnetic field. Our multiple-radiofrequency approach gives precise control over the double well characteristics, including the depth of individual wells and the height of the barrier, and enables reliable transfer of atoms between the available trapping geometries. We have characterised the multiple-radiofrequency dressed system using radiofrequency spectroscopy, finding good agreement with the eigenvalues numerically calculated using Floquet theory. This method creates trapping potentials that can be reconfigured by changing the amplitudes, polarizations and frequencies of the applied dressing fields, and easily extended with additional dressing frequencies.


title: Asymmetry Observables for Measuring Spin Correlations in Top-Quark Pair Production

abstract: Spin correlations between top-quark and anti-top-quark in pair production are sensitive to physics beyond the standard model. A complete set of observables, that can be used to measure all independent parameters of the production spin density matrix, is necessary to search for anomalous interactions in the production dynamics. In this letter, we provide a set of asymmetry observables for measuring all the density matrix elements. These observables are defined regardless of the production mechanism, and can be specified in any desired reference frame, particularly the laboratory frame in which the first evidence of spin correlation was observed experimentally. Our method can be used to find an optimal Lorentz frame in which the sensitivity to the density matrices is high for given production process and experimental environment, or to test possible correlations among the matrix elements.


title: Monolithic Photoelectrochemical Device for 19% Direct Water Splitting

abstract: Recent rapid progress in efficiencies for solar water splitting by photoelectrochemical devices has enhanced its prospects to enable storable renewable energy. Efficient solar fuel generators all use tandem photoelectrode structures, and advanced integrated devices incorporate corrosion protection layers as well as heterogeneous catalysts. Realization of near thermodynamic limiting performance requires tailoring the energy band structure of the photoelectrode and also the optical and electronic properties of the surface layers exposed to the electrolyte. Here, we report a monolithic device architecture that exhibits reduced surface reflectivity in conjunction with metallic Rh nanoparticle catalyst layers that minimize parasitic light absorption. Additionally, the anatase TiO2 protection layer on the photocathode creates a favorable internal band alignment for hydrogen evolution. An initial solar-to-hydrogen efficiency of 19.3 % is obtained in acidic electrolyte and an efficiency of 18.5 % is achieved at neutral pH condition (under simulated sunlight).


title: Characterization of optimal carbon nanotubes under stretching and validation of the Cauchy-Born rule

abstract: Carbon nanotubes are modeled as point configurations and investigated by minimizing configurational energies including two-and three-body interactions. Optimal configurations are identified with local minima and their fine geometry is fully characterized in terms of lower-dimensional problems. Under moderate tension, we prove the existence of periodic local minimizers, which indeed validates the so-called Cauchy-Born rule in this setting.


title: $\Delta N_{\text{eff}}$ and entropy production from early-decaying gravitinos

abstract: Gravitinos are a fundamental prediction of supergravity, their mass ($m_{G}$) is informative of the value of the SUSY breaking scale, and, if produced during reheating, their number density is a function of the reheating temperature ($T_{\text{rh}}$). As a result, constraining their parameter space provides in turn significant constraints on particles physics and cosmology. We have previously shown that for gravitinos decaying into photons or charged particles during the ($\mu$ and $y$) distortion eras, upcoming CMB spectral distortions bounds are highly effective in constraining the $T_{\text{rh}}-m_{G}$ space. For heavier gravitinos (with lifetimes shorter than a few $\times10^6$ sec), distortions are quickly thermalized and energy injections cause a temperature rise for the CMB bath. If the decay occurs after neutrino decoupling, its overall effect is a suppression of the effective number of relativistic degrees of freedom ($N_{\text{eff}}$). In this paper, we utilize the observational bounds on $N_{\text{eff}}$ to constrain gravitino decays, and hence provide new constaints on gravitinos and reheating. For gravitino masses less than $\approx 10^5$ GeV, current observations give an upper limit on the reheating scale in the range of $\approx 5 \times 10^{10}- 5 \times 10^{11}$GeV. For masses greater than $\approx 4 \times 10^3$ GeV they are more stringent than previous bounds from BBN constraints, coming from photodissociation of deuterium, by almost 2 orders of magnitude.


title: 3D spherical-cap fitting procedure for (truncated) sessile nano- and micro-droplets & -bubbles

abstract: In the study of nanobubbles, nanodroplets or nanolenses immobilised on a substrate, a cross-section of a spherical-cap is widely applied to extract geometrical information from atomic force microscopy (AFM) topographic images. In this paper, we have developed a comprehensive 3D spherical cap fitting procedure (3D-SCFP) to extract morphologic characteristics of complete or truncated spherical caps from AFM images. Our procedure integrates several advanced digital image analysis techniques to construct a 3D spherical cap model, from which the geometrical parameters of the nanostructures are extracted automatically by a simple algorithm. The procedure takes into account all valid data points in the construction of the 3D spherical cap model to achieve high fidelity in morphology analysis. We compare our 3D fitting procedure with the commonly used 2D cross-sectional profile fitting method to determine the contact angle of a complete spherical cap and a truncated spherical cap. The results from 3D-SCFP are consistent and accurate, while 2D fitting is unavoidably arbitrary in selection of the cross-section and has a much lower number of data points on which the fitting can be based, which in addition is biased to the top of the spherical cap. We expect that the developed 3D spherical-cap fitting procedure will find many applications in imaging analysis.


title: FCS and RICS Spectra of Probes in Complex Fluids

abstract: The Fluorescence Correlation Spectroscopy (FCS) spectrum G(t) and Raster Image Correlation Spectroscopy (RICS) spectrum R(t) of dilute diffusing particles are determined by the displacement distribution function P(x,t) of the particles and by the experimental parameters of the associated optical trains. This letter obtains the general relationships between P(x,t) and these spectra. For dilute diffusing molecules in simple liquids, P(x,t) is a Gaussian in the displacement x; the corresponding G(t) is a Lorentzian in (<(x(t))^2>)^(1/2). In complex fluids such as polymer solutions, colloid and protein solutions, and the interior of living cells, P(x,t) may have a non-Gaussian dependence on x, for example an exponential in |x|. We compare theoretical forms for FCS and for RICS spectra of two systems in which P(x,t) is a Gaussian or an exponential in x, but in which the mean-square displacements are precisely equal at all times. If the G(t) and R(t) arising from an exponential P(x,t) are interpreted by using the forms for G(t) and R(t) that are appropriate for a Gaussian P(x,t), the inferred diffusion coefficient may be substantially in error.


title: Stochastic Gradient Monomial Gamma Sampler

abstract: Recent advances in stochastic gradient techniques have made it possible to estimate posterior distributions from large datasets via Markov Chain Monte Carlo (MCMC). However, when the target posterior is multimodal, mixing performance is often poor. This results in inadequate exploration of the posterior distribution. A framework is proposed to improve the sampling efficiency of stochastic gradient MCMC, based on Hamiltonian Monte Carlo. A generalized kinetic function is leveraged, delivering superior stationary mixing, especially for multimodal distributions. Techniques are also discussed to overcome the practical issues introduced by this generalization. It is shown that the proposed approach is better at exploring complex multimodal posterior distributions, as demonstrated on multiple applications and in comparison with other stochastic gradient MCMC methods.


title: Rogers-Shephard and local Loomis-Whitney type inequalities

abstract: We provide functional analogues of the classical geometric inequality of Rogers and Shephard on products of volumes of sections and projections. As a consequence we recover (and obtain some new) functional versions of Rogers-Shephard type inequalities as well as some generalizations of the geometric Rogers-Shephard inequality in the case where the subspaces intersect. These generalizations can be regarded as sharp local reverse Loomis-Whitney inequalities. We also obtain a sharp local Loomis-Whitney inequality.


title: Readings and Misreadings of J. Willard Gibbs Elementary Principles in Statistical Mechanics

abstract: J. Willard Gibbs' Elementary Principles in Statistical Mechanics was the definitive work of one of America's greatest physicists. Gibbs' book on statistical mechanics establishes the basic principles and fundamental results that have flowered into the modern field of statistical mechanics. However, at a number of points, Gibbs' teachings on statistical mechanics diverge from positions on the canonical ensemble found in more recent works, at points where seemingly there should be agreement. The objective of this paper is to note some of these points, so that Gibbs' actual positions are not misrepresented to future generations of students.


title: Direction Dependent Corrections in Polarimetric Radio Imaging I : Characterizing the effects of the primary beam on full Stokes imaging

abstract: Next generation radio telescope arrays are being designed and commissioned to accurately measure polarized intensity and rotation measures across the entire sky through deep, wide-field radio interferometric surveys. Radio interferometer dish antenna arrays are affected by direction-dependent (DD) gains due to both instrumental and atmospheric effects. In this paper we demonstrate the effect of DD errors for parabolic dish antenna array on the measured polarized intensities of radio sources in interferometric images. We characterize the extent of polarimetric image degradation due to the DD gains through wide-band VLA simulations of representative point source simulations of the radio sky at L-Band(1-2GHz). We show that at the 0.5 gain level of the primary beam (PB) there is significant flux leakage from Stokes $I$ to $Q$, $U$ amounting to 10\% of the total intensity. We further demonstrate that while the instrumental response averages down for observations over large parallactic angle intervals, full-polarization DD correction is required to remove the effects of DD leakage. We also explore the effect of the DD beam on the Rotation Measure(RM) signals and show that while the instrumental effect is primarily centered around 0 rad-m$^{-2}$, the effect is significant over a broad range of RM requiring full polarization DD correction to accurately reconstruct RM synthesis signal.


title: UCB Exploration via Q-Ensembles

abstract: We show how an ensemble of $Q^*$-functions can be leveraged for more effective exploration in deep reinforcement learning. We build on well established algorithms from the bandit setting, and adapt them to the $Q$-learning setting. We propose an exploration strategy based on upper-confidence bounds (UCB). Our experiments show significant gains on the Atari benchmark.


title: Unipolar Induction Star-Disk Interaction

abstract: Some critical comments on the prevailing model of star-disk interaction are made, in particular, on the rotating nature of the magnetic field lines and on the application of the magnetohydrodynamic frozen-field theorem to the disk plasma. As an alternative, a unipolar induction model is proposed, where the magnetic field is stationary in space and the stellar unipolar electric field $\vec E_{*}$ on the surface is uploaded to the magnetosphere. Through the Poynting vector, the star and the magnetosphere form a coupled system where the total angular momentum, consisting the mechanical one of the star, the electromagnetic one of the magnetosphere, and the mechanical one of the plasma in the magnetosphere, is conserved. The stellar interaction with the accretion disk is through the projection of the unipolar electric field $\vec E_{*}$ onto the disk via the equipotential field lines, generating a disk current and consequently toroidal fields with opposite signs on both sides of the disk, with return current loops via the stellar surface. As a result, magnetic flux is added to the magnetospheric field in the northern and southern hemispheres with the disk current sheet as the boundary condition. This makes the star-disk system an astrophysical site where intense magnetic fields are generated through the rotational energy of the star. Angular momentum extraction from this star-disk system happens as the magnetic flux in the magnetosphere increases to the point that exceeds the current carrying capacity of the disk, leading to a mega scale magnetic eruption sending Poynting fluxes to space either isotropically or beamed.


title: The chemical composition of the stellar cluster Gaia1: no surprise behind Sirius

abstract: We observed 6 He-clump stars of the intermediate-age stellar cluster Gaia1 with the MIKE/MAGELLAN spectrograph. A possible extra-galactic origin of this cluster, recently discovered thanks to the first data release of the ESA Gaia mission, has been suggested, based on its orbital parameters. Abundances for Fe, alpha, proton- and neutron-capture elements have been obtained. We find no evidence of intrinsic abundance spreads. The iron abundance is solar ([FeI/H]=+0.00 +-0.01; sigma = 0.03 dex). All the other abundance ratios are, by and large, solar-scaled, similar to the Galactic thin disk and open clusters stars of similar metallicity. The chemical composition of Gaia1 does not support an extra-galactic origin for this stellar cluster, that can be considered as a standard Galactic open cluster.


title: Sensitivity of the Cherenkov Telescope Array to the detection of a dark matter signal in comparison to direct detection and collider experiments

abstract: Imaging atmospheric Cherenkov telescopes (IACTs) that are sensitive to potential $\gamma$-ray signals from dark matter (DM) annihilation above $\sim50$ GeV will soon be superseded by the Cherenkov Telescope Array (CTA). CTA will have a point source sensitivity an order of magnitude better than currently operating IACTs and will cover a broad energy range between 20 GeV and 300 TeV. Using effective field theory and simplified models to calculate $\gamma$-ray spectra resulting from DM annihilation, we compare the prospects to constrain such models with CTA observations of the Galactic center with current and near-future measurements at the Large Hadron Collider (LHC) and direct detection experiments. For DM annihilations via vector or pseudoscalar couplings, CTA observations will be able to probe DM models out of reach of the LHC, and, if DM is coupled to standard fermions by a pseudoscalar particle, beyond the limits of current direct detection experiments.


title: Discovery of three strongly lensed quasars in the Sloan Digital Sky Survey

abstract: We present the discovery of 3 quasar lenses in the Sloan Digital Sky Survey (SDSS), selected using two novel photometry-based selection techniques. The J0941+0518 system, with two point sources separated by 5.46" on either side of a galaxy, has source and lens redshifts $z_s = 1.54$ and $z_l = 0.343$. The AO-assisted images of J2211+1929 show two point sources separated by 1.04", corresponding to the same quasar at $z_s = 1.07,$ besides the lens galaxy and Einstein ring. Images of J2257+2349 show two point sources separated by 1.67" on either side of an E/S0 galaxy. The extracted spectra show two images of the same quasar at redshift $z_s = 2.10$. In total, the two selection techniques identified 309 lens candidates, including 47 known lenses, and 6 previously ruled out candidates. 55 of the remaining candidates were observed using NIRC2 and ESI at Keck Observatory, EFOSC2 at the ESO-NTT (La Silla), and SAM and the Goodman spectrograph at SOAR. Of the candidates observed, 3 were confirmed as lenses, 36 were ruled out, and 16 remain inconclusive. Taking into account that we recovered known lenses, this gives us a success rate of at least 50/309 (16%). This initial campaign demonstrates the power of purely photometric selection techniques in finding lensed quasars. Developing and refining these techniques is essential for efficient identification of these rare lenses in ongoing and future photometric surveys.


title: Density Deconvolution for Generalized Skew-Symmetric Distributions

abstract: This paper develops a density deconvolution estimator that assumes the density of interest is a member of the generalized skew-symmetric (GSS) family of distributions. Estimation occurs in two parts: a skewing function, as well as location and scale parameters must be estimated. A kernel method is proposed for estimating the skewing function. The mean integrated square error (MISE) of the resulting GSS deconvolution estimator is derived. Based on derivation of the MISE, two bandwidth estimation methods for estimating the skewing function are also proposed. A generalized method of moments (GMM) approach is developed for estimation of the location and scale parameters. The question of multiple solutions in applying the GMM is also considered, and two solution selection criteria are proposed. The GSS deconvolution estimator is further investigated in simulation studies and is compared to the nonparametric deconvolution estimator. For most simulation settings considered, the GSS estimator has performance superior to the nonparametric estimator.


title: Time-dependent shortest paths in bounded treewidth graphs

abstract: We present a proof that the number of breakpoints in the arrival function between two terminals in graphs of treewidth $w$ is $n^{O(\log^2 w)}$ when the edge arrival functions are piecewise linear. This is an improvement on the bound of $n^{\Theta(\log n)}$ by Foschini, Hershberger, and Suri for graphs without any bound on treewidth. We provide an algorithm for calculating this arrival function using star-mesh transformations, a generalization of the wye-delta-wye transformations.


title: Facial Emotion Detection Using Convolutional Neural Networks and Representational Autoencoder Units

abstract: Emotion being a subjective thing, leveraging knowledge and science behind labeled data and extracting the components that constitute it, has been a challenging problem in the industry for many years. With the evolution of deep learning in computer vision, emotion recognition has become a widely-tackled research problem. In this work, we propose two independent methods for this very task. The first method uses autoencoders to construct a unique representation of each emotion, while the second method is an 8-layer convolutional neural network (CNN). These methods were trained on the posed-emotion dataset (JAFFE), and to test their robustness, both the models were also tested on 100 random images from the Labeled Faces in the Wild (LFW) dataset, which consists of images that are candid than posed. The results show that with more fine-tuning and depth, our CNN model can outperform the state-of-the-art methods for emotion recognition. We also propose some exciting ideas for expanding the concept of representational autoencoders to improve their performance.


title: Commutants of Toeplitz Operators with Harmonic Symbols

abstract: We make a progress towards describing the commutants of Toeplitz operators with harmonic symbols on the Bergman space over the unit disk. Our work greatly generalizes several partial results in the field.


title: Proper-time Quantum Mechanics for Multi-Quark System and Composite-Hadron Spectroscopy

abstract: One of the most important problem in hadron physics is to establish the Lorentz-invariant classification scheme of composite hadrons, extending the framework of non-relativistic quark model. We present an attempt, by developing proper-time $\tau$ quantum mechanics on a multi-quark system in particle frame (with constant boost velocity $\boldsymbol{v}$). We start from the variational method on a classical mechanics action where a constituent quark has Pauli-type $SU(2)_{\sigma}$ spin. Then the $SU(2)_{\mathfrak{m}}$ symmetry, concerning the sign-reversal on quark mass, has arisen with the basic vectors, the normal Dirac spinor with $J^{P}=(1/2)^{+}$ and the chiral one with $J^{P}=(1/2)^{-}$, appearing as a "shadow" of the former. Herewith, the mass reversal between these basic vectors become equivalent to the chirality, which is a symmetry of the standard gauge theory. We describe the role of chirality in hadron spectroscopy and regard it as attribute {$\chi$} of "elementary" hadrons in addition to {$J, P, C$}. A novel feature of our hadron spectroscopy is, in the example of $q\bar{q}$ meson system, that the "Regge trajectories", are given by mass-squared vs. the number of quantum $N$ ; where $M^2 =M_{0}^2 +2N\Omega$ ($N=2n$, $n$ the radial quantum number, $\Omega$ the oscillator quantum), and the intrinsic spin of hadrons $\boldsymbol{J}$ comes only from quark spin $\boldsymbol{S}$, $\boldsymbol{J}=\boldsymbol{S}$. Some phenomenological facts crucial to its validity are pointed out on the light-through-heavy quarkonium system.


title: Feynman-Kac Equations for Reaction and Diffusion Processes

abstract: This paper provides a theoretical framework of deriving the forward and backward Feynman-Kac equations for the distribution of functionals of the path of a particle undergoing both diffusion and chemical reaction. Very general forms of the equations are obtained. Once given the diffusion type and reaction rate, a specific forward or backward Feynman-Kac equation can be obtained. The listed in the paper include the ones for normal/anomalous diffusions and reactions with linear/nonlinear rates. Using the derived equations, we also study the occupation time in half-space, the first passage time to a fixed boundary, and the occupation time in half-space with absorbing or reflecting boundary conditions.


title: Beyond Volume: The Impact of Complex Healthcare Data on the Machine Learning Pipeline

abstract: From medical charts to national census, healthcare has traditionally operated under a paper-based paradigm. However, the past decade has marked a long and arduous transformation bringing healthcare into the digital age. Ranging from electronic health records, to digitized imaging and laboratory reports, to public health datasets, today, healthcare now generates an incredible amount of digital information. Such a wealth of data presents an exciting opportunity for integrated machine learning solutions to address problems across multiple facets of healthcare practice and administration. Unfortunately, the ability to derive accurate and informative insights requires more than the ability to execute machine learning models. Rather, a deeper understanding of the data on which the models are run is imperative for their success. While a significant effort has been undertaken to develop models able to process the volume of data obtained during the analysis of millions of digitalized patient records, it is important to remember that volume represents only one aspect of the data. In fact, drawing on data from an increasingly diverse set of sources, healthcare data presents an incredibly complex set of attributes that must be accounted for throughout the machine learning pipeline. This chapter focuses on highlighting such challenges, and is broken down into three distinct components, each representing a phase of the pipeline. We begin with attributes of the data accounted for during preprocessing, then move to considerations during model building, and end with challenges to the interpretation of model output. For each component, we present a discussion around data as it relates to the healthcare domain and offer insight into the challenges each may impose on the efficiency of machine learning techniques.


title: Emergent Network Modularity

abstract: We introduce a network growth model based on complete redirection: a new node randomly selects an existing target node, but attaches to a random neighbor of this target. For undirected networks, this simple growth rule generates unusual, highly modular networks. Individual network realizations typically contain multiple macrohubs---nodes whose degree scales linearly with the number of nodes $N$. The size of the network "nucleus"---the set of nodes of degree greater than one---grows sublinearly with $N$ and thus constitutes a vanishingly small fraction of the network. The network therefore consists almost entirely of leaves (nodes of degree one) as $N\to\infty$.


title: Limitations in the $2D$ description of the electromagnetic waves propagation in thin dielectric and magnetic layers

abstract: The propagation of electromagnetic waves trapped within dielectric and magnetic layers is considered. The description within the three-dimensional theory is compared with the simplified analysis in two dimensions. Two distinct media configurations with different topology are dealt with: a plane slab and a hollow cylinder. Choosing the appropriate values for the geometrical parameters (layer thickness, radius of the cylinder) and for the electromagnetic properties of the media one can trap exactly one mode corresponding to that obtained within the two-dimensional electromagnetism. However, the symmetry between electric and magnetic fields suggests, that the two versions of the simplified electromagnetism ought to be taken into account. Its usual form is incomplete to describe all modes. It is also found that there is a domain of optimal values of parameters for which the $2D$ model works relatively correctly. In the case of a cylindrical surface we observe, however, several differences which are attributed to the curvature of the layer, and which exclude the propagation of evanescent modes. The two-dimensional electrodynamics, whichever form used, turns out still too poor to describe the so called `hybrid modes' excited in a real layer. The obtained results can be important for proper description of the propagating waves within thin layers for which $3D$ approach is not available due to mathematical complexity, and reducing the layer to a lower-dimensional structure seems the only possible option.


title: Models of the Hydrodynamic Histories of post-AGB Stars. I. Multiflow Shaping of OH231.8+04.2

abstract: We present a detailed hydrodynamic model that matches the present structure of the well-observed preplanetary nebula OH231.8+04.2. The purpose of the model is to present a physically justified and coherent picture of its evolutionary history from about 100 years of the start of the formation of its complex outer structures to the present. We have adopted a set of initial conditions that are heavily constrained by high-quality observations of its present structure and kinematics. The shaping of the nebula occurs while the densities of the flows are "light": less than the surrounding AGB-wind environment. The simulations show that pairs of essentially coeval clumps and sprays of the same extent and density but different outflow speeds sculpted both the pair of thin axial flow "or spine" and the bulbs. The total ejected mass and momentum in the best fit model are surprisingly large -- 3 M_sun and 2.2 x 10^41 gm cm s-1, respectively -- however, these values are reduced by up to a factor of ten in other models that fit the data almost as well. Our ultimate goal is to combine the present model results of masses, momenta, flow speeds, and flow geometries for OH231 with those of other models to be published in the future in order to find common attributes of their ejection histories.


title: Wide Bandwidth, Frequency Modulated Free Electron Laser

abstract: It is shown via theory and simulation that the resonant frequency of a Free Electron Laser may be modulated to obtain an FEL interaction with a frequency bandwidth which is at least an order of magnitude greater than normal FEL operation. The system is described in the linear regime by a summation over exponential gain modes, allowing the amplification of multiple light frequencies simultaneously. Simulation in 3D demonstrates the process for parameters of the UK's CLARA FEL test facility currently under construction. This new mode of FEL operation has close analogies to Frequency Modulation in a conventional cavity laser. This new, wide bandwidth mode of FEL operation scales well for X-ray generation and offers users a new form of high-power FEL output.


title: Degeneration of Kahler-Einstein manifolds of negative scalar curvature

abstract: Let $\pi: \mathcal{X}^* \rightarrow B^*$ be an algebraic family of compact K\"ahler manifolds of complex dimension $n$ with negative first Chern class over a punctured disc $B^*\in \mathbb{C}$. Let $g_t$ be the unique K\"ahler-Einstein metric on $\mathcal{X}_t= \pi^{-1}(t)$. We show that as $t\rightarrow 0$, $(\mathcal{X}_t, g_t)$ converges in pointed Gromov-Hausdorff topology to a unique finite disjoint union of complete metric length spaces $\coprod_{\alpha=1}^\mathcal{A} (Y_\alpha, d_\alpha)$ without loss of volume. Each $(Y_\alpha, d_\alpha)$ is a smooth open K\"ahler-Einstein manifold of complex dimension n outside its closed singular set of Hausdorff dimension no greater than $2n-4$. Furthermore, $\coprod_{\alpha=1}^\mathcal{A} Y_\alpha$ is a quasi-projective variety isomorphic to $\mathcal{X}_0 \setminus LCS(\mathcal{X}_0)$, where $\mathcal{X}_0$ is a projective semi-log canonical model and $LCS(\mathcal{X}_0)$ is the non-log terminal locus of $\mathcal{X}_0$. This is the first step of our approach toward compactification of the analytic geometric moduli space of K\"ahler-Einstein manifolds of negative scalar curvature.


title: Additive decomposition of iterative quantum search operations in the Grover-type algorithm

abstract: In the Grover-type quantum search process a search operator is iteratively applied, say, k times, on the initial database state. We present an additive decomposition scheme such that the iteration process is expressed, in the computational space, as a linear combination of k operators, each of which consists of a single Grover-search followed by an overall phase-rotation. The value of k and the rotation phase are the same as those determined in the framework of the search with certainty. We further show that the final state can be expressed in terms of a single oracle operator of the Grover-search and phase-rotation factors. We discuss how the additive form can be utilized so that it effectively reduces the computational load of the iterative search, and we propose an effective shortcut gate that realizes the same outcome as the iterative search.


title: The Convergence of Markov chain Monte Carlo Methods: From the Metropolis method to Hamiltonian Monte Carlo

abstract: From its inception in the 1950s to the modern frontiers of applied statistics, Markov chain Monte Carlo has been one of the most ubiquitous and successful methods in statistical computing. In that time its development has been fueled by increasingly difficult problems and novel techniques from physics. In this article I will review the history of Markov chain Monte Carlo from its inception with the Metropolis method to today's state-of-the-art in Hamiltonian Monte Carlo. Along the way I will focus on the evolving interplay between the statistical and physical perspectives of the method.


title: Will Break for Productivity: Generalized Symptoms of Cognitive Depletion

abstract: In this work, we address the symptoms of cognitive depletion as they relate to generalized knowledge workers. We unify previous findings within a single analytical model of cognitive depletion. Our purpose is to develop a model that will help us predict when a person has reached a sufficient state of cognitive depletion such that taking a break or some other restorative action will benefit both his or her own wellbeing and the quality of his or her performance. We provide a definition of each symptom in our model as well as the effect it would have on a knowledge worker's ability to work productively. We discuss methods to detect each symptom that do not require self assessment. Understanding symptoms of cognitive depletion provides the ability to support human knowledge workers by reducing the stress involved with cognitive and work overload while maintaining or improving the quality of their performance.


title: Forecasted masses for seven thousand KOIs

abstract: Recent transit surveys have discovered thousands of planetary candidates with directly measured radii, but only a small fraction have measured masses. Planetary mass is crucial in assessing the feasibility of numerous observational signatures, such as radial velocities (RVs), atmospheres, moons and rings. In the absence of a direct measurement, a data-driven, probabilistic forecast enables observational planning and so here we compute posterior distributions for the forecasted mass of approximately seven thousand Kepler Objects of Interest (KOIs). Our forecasts reveal that the predicted RV amplitudes of Neptunian planets are relatively consistent, as a result of transit survey detection bias, hovering around the few m/s level. We find that mass forecasts are unlikely to improve through more precise planetary radii, with the error budget presently dominated by the intrinsic model uncertainty. Our forecasts identify a couple of dozen KOIs near the Terran-Neptunian divide with particularly large RV semi-amplitudes which could be promising targets to follow-up, particularly in the near-IR. With several more transit surveys planned in the near-future, the need to quickly forecast observational signatures is likely to grow and the work here provides a template example of such calculations.


title: Cognitive Depletion in the Wild: a Case Study of NMR Spectroscopy Analysis

abstract: NMR spectroscopy analysis is a detail-oriented analytic feat that typically requires specific domain expertise and hours of concentration. This work presents an ethnographic-style study of this analysis process in the context of evaluating the symptoms of cognitive depletion. The repeated, non-trivial decisions required by and the time-consuming nature of NMR spectroscopy analysis make it an ideal, real-world scenario to study the symptoms of cognitive depletion, its effect on workflow and performance, and potential strategies for mitigating its deleterious effects.


title: Renormalization theory for the Fulde-Ferrell-Larkin-Ovchinnikov states at $T>0$

abstract: Within the renormalization group framework we study the stability of superfluid density wave states, known as Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) phases, with respect to thermal order-parameter fluctuations in two and three-dimensional ($d\in \{2,3\}$) systems. We analyze the renormalization-group flow of the relevant ordering wave-vector $\vec{Q_0}$. The calculation indicates an instability of the FFLO-type states towards either a uniform superfluid or the normal state in $d\in\{2,3\}$ and $T>0$. In $d=2$ this is signaled by $\vec{Q_0}$ being renormalized towards zero, corresponding to the flow being attracted either to the usual Kosterlitz-Thouless fixed-point or to the normal phase. We supplement a solution of the RG flow equations by a simple scaling argument, supporting the generality of the result. The tendency to reduce the magnitude of $\vec{Q_0}$ by thermal fluctuations persists in $d=3$, where the very presence of long-range order is immune to thermal fluctuations, but the effect of attracting $\vec{Q_0}$ towards zero by the flow remains observed at $T>0$.

categories: cond-mat.quant-gas cond-mat.str-el cond-mat.supr-con hep-th


title: The spectrum of a vertex model and related spin one chain sitting in a genus five curve

abstract: We derive the transfer matrix eigenvalues of a three-state vertex model whose weights are based on a $\mathrm{R}$-matrix not of difference form with spectral parameters lying on a genus five curve. We have shown that the basic building blocks for both the transfer matrix eigenvalues and Bethe equations can be expressed in terms of meromorphic functions on an elliptic curve. We discuss the properties of an underlying spin one chain originated from a particular choice of the $\mathrm{R}$-matrix second spectral parameter. We present numerical and analytical evidences that the respective low-energy excitations can be gapped or massless depending on the strength of the interaction coupling. In the massive phase we provide analytical and numerical evidences in favor of an exact expression for the lowest energy gap. We point out that the critical point separating these two distinct physical regimes coincides with the one in which the weights geometry degenerate into union of genus one curves.


title: Quantum control of topological defects in magnetic systems

abstract: Energy-efficient classical information processing and storage based on topological defects in magnetic systems have been studied over past decade. In this work, we introduce a class of macroscopic quantum devices in which a quantum state is stored in a topological defect of a magnetic insulator. We propose non-invasive methods to coherently control and readout the quantum state using ac magnetic fields and magnetic force microscopy, respectively. This macroscopic quantum spintronic device realizes the magnetic analog of the three-level rf-SQUID qubit and is built fully out of electrical insulators with no mobile electrons, thus eliminating decoherence due to the coupling of the quantum variable to an electronic continuum and energy dissipation due to Joule heating. For a domain wall sizes of $10-100$~nm and reasonable material parameters, we estimate qubit operating temperatures in the range of $0.1-1$~K, a decoherence time of about $0.01-1$~$\mu$s, and the number of Rabi flops within the coherence time scale in the range of $10^{2}-10^{4}$.


title: Geometric estimates for complex Monge-Ampere equations

abstract: We prove uniform gradient and diameter estimates for a family of geometric complex Monge-Ampere equations. Such estimates can be applied to study geometric regularity of singular solutions of complex Monge-Ampere equations. We also prove a uniform diameter estimate for collapsing families of twisted Kahler-Einstein metrics on Kahler manifolds of nonnegative Kodaira dimensions.


title: Bose-Einstein Condensation in a Plasmonic Lattice

abstract: Bose-Einstein condensation is a remarkable manifestation of quantum statistics and macroscopic quantum coherence. Superconductivity and superfluidity have their origin in Bose-Einstein condensation. Ultracold quantum gases have provided condensates close to the original ideas of Bose and Einstein, while condensation of polaritons and magnons have introduced novel concepts of non-equilibrium condensation. Here, we demonstrate a Bose-Einstein condensate (BEC) of surface plasmon polaritons in lattice modes of a metal nanoparticle array. Interaction of the nanoscale-confined surface plasmons with a room-temperature bath of dye molecules enables thermalization and condensation in picoseconds. The ultrafast thermalization and condensation dynamics are revealed by an experiment that exploits thermalization under propagation and the open cavity character of the system. A crossover from BEC to usual lasing is realized by tailoring the band structure. This new condensate of surface plasmon lattice excitations has promise for future technologies due to its ultrafast, room-temperature and on-chip nature.


title: Quantifying fermionic decoherence in many-body systems

abstract: Practical measures of electronic decoherence, called distilled purities, that are applicable to many-body systems are introduced. While usual measures of electronic decoherence such as the purity employ the full $N$-particle density matrix which is generally unavailable, the distilled purities are based on the $r$-body reduced density matrices ($r$-RDMs) which are more accessible quantities. The $r$-body distilled purities are derivative quantities of the previously introduced $r$-body reduced purities [I. Franco and H. Appel, J. Chem. Phys. 139, 094109 (2013)] that measure the non-idempotency of the $r$-RDMs. Specifically, the distilled purities exploit the structure of the reduced purities to extract coherences between Slater determinants with integer occupations defined by a given single-particle basis that compose an electronic state. In this way, the distilled purities offer a practical platform to quantify coherences in a given basis that can be used to analyze the quantum dynamics of many-electron systems. Exact expressions for the one-body and two-body distilled purities are presented and the utility of the approach is exemplified via an analysis of the dynamics of oligo-acetylene as described by the Su-Schrieffer-Heeger Hamiltonian. Last, the advantages and limitations of the purity, reduced purity and distilled purity as measures of electronic coherence are discussed.


title: On the $L^p$ boundedness of wave operators for two-dimensional Schr\"odinger operators with threshold obstructions

abstract: Let $H=-\Delta+V$ be a Schr\"odinger operator on $L^2(\mathbb R^2)$ with real-valued potential $V$, and let $H_0=-\Delta$. If $V$ has sufficient pointwise decay, the wave operators $W_{\pm}=s-\lim_{t\to \pm\infty} e^{itH}e^{-itH_0}$ are known to be bounded on $L^p(\mathbb R^2)$ for all $1< p< \infty$ if zero is not an eigenvalue or resonance. We show that if there is an s-wave resonance or an eigenvalue only at zero, then the wave operators are bounded on $L^p(\mathbb R^2)$ for $1 < p<\infty$. This result stands in contrast to results in higher dimensions, where the presence of zero energy obstructions is known to shrink the range of valid exponents $p$.


title: Progressive Boosting for Class Imbalance

abstract: Pattern recognition applications often suffer from skewed data distributions between classes, which may vary during operations w.r.t. the design data. Two-class classification systems designed using skewed data tend to recognize the majority class better than the minority class of interest. Several data-level techniques have been proposed to alleviate this issue by up-sampling minority samples or under-sampling majority samples. However, some informative samples may be neglected by random under-sampling and adding synthetic positive samples through up-sampling adds to training complexity. In this paper, a new ensemble learning algorithm called Progressive Boosting (PBoost) is proposed that progressively inserts uncorrelated groups of samples into a Boosting procedure to avoid loss of information while generating a diverse pool of classifiers. Base classifiers in this ensemble are generated from one iteration to the next, using subsets from a validation set that grows gradually in size and imbalance. Consequently, PBoost is more robust to unknown and variable levels of skew in operational data, and has lower computation complexity than Boosting ensembles in literature. In PBoost, a new loss factor is proposed to avoid bias of performance towards the negative class. Using this loss factor, the weight update of samples and classifier contribution in final predictions are set based on the ability to recognize both classes. Using the proposed loss factor instead of standard accuracy can avoid biasing performance in any Boosting ensemble. The proposed approach was validated and compared using synthetic data, videos from the FIA dataset that emulates face re-identification applications, and KEEL collection of datasets. Results show that PBoost can outperform state of the art techniques in terms of both accuracy and complexity over different levels of imbalance and overlap between classes.


title: Coupling emission from single localized defects in 2D semiconductor to surface plasmon polaritons

abstract: Coupling of an atom-like emitter to surface plasmons provides a path toward significant optical nonlinearity, which is essential in quantum information processing and quantum networks. A large coupling strength requires nanometer-scale positioning accuracy of the emitter near the surface of the plasmonic structure, which is challenging. We demonstrate the coupling of single localized defects in a tungsten diselenide (WSe2) monolayer self-aligned to the surface plasmon mode of a silver nanowire. The silver nanowire induces a strain gradient on the monolayer at the overlapping area, leading to the formation of localized defect emission sites that are intrinsically close to the surface plasmon. We measure a coupling efficiency with a lower bound of 39% from the emitter into the plasmonic mode of the silver nanowire. This technique offers a way to achieve efficient coupling between plasmonic structures and localized defects of 2D semiconductors.